CN117133351B - 一种优化的线粒体基因重排量化方法 - Google Patents
一种优化的线粒体基因重排量化方法 Download PDFInfo
- Publication number
- CN117133351B CN117133351B CN202311378396.9A CN202311378396A CN117133351B CN 117133351 B CN117133351 B CN 117133351B CN 202311378396 A CN202311378396 A CN 202311378396A CN 117133351 B CN117133351 B CN 117133351B
- Authority
- CN
- China
- Prior art keywords
- gene
- gene sequence
- array
- sequence
- rearrangement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000008707 rearrangement Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 28
- 108020005196 Mitochondrial DNA Proteins 0.000 title claims abstract description 22
- 238000011002 quantification Methods 0.000 title claims description 7
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 272
- 230000002438 mitochondrial effect Effects 0.000 claims abstract description 26
- 238000013139 quantization Methods 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract description 4
- 238000003491 array Methods 0.000 claims abstract description 3
- 101150076211 TH gene Proteins 0.000 claims description 18
- 101150077543 st gene Proteins 0.000 claims description 18
- 241001465754 Metazoa Species 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 5
- 108091008053 gene clusters Proteins 0.000 description 7
- 241000251468 Actinopterygii Species 0.000 description 5
- 241000271566 Aves Species 0.000 description 5
- 208000037516 chromosome inversion disease Diseases 0.000 description 5
- 241000124008 Mammalia Species 0.000 description 4
- 230000010076 replication Effects 0.000 description 4
- 241000270322 Lepidosauria Species 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 108091036078 conserved sequence Proteins 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000006798 recombination Effects 0.000 description 3
- 238000005215 recombination Methods 0.000 description 3
- 241000894007 species Species 0.000 description 3
- 241000938605 Crocodylia Species 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 101100047216 Mus musculus Trnp1 gene Proteins 0.000 description 1
- 101150117428 NAD6 gene Proteins 0.000 description 1
- 108091008109 Pseudogenes Proteins 0.000 description 1
- 102000057361 Pseudogenes Human genes 0.000 description 1
- 241000251539 Vertebrata <Metazoa> Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
Landscapes
- Bioinformatics & Cheminformatics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Chemical & Material Sciences (AREA)
- Molecular Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Bioinformatics & Computational Biology (AREA)
- Analytical Chemistry (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于分子生物学和生物信息学技术领域,公开了一种优化的线粒体基因重排量化方法,输入一条基准基因序列和待量化基因序列;待量化基因序列逐列与基准基因序列进行异或运算,将运算结果按差异程度由高到低存储到划分数组中;取划分数组中差异程度最高的N个位置将基准基因序列和每个待量化基因序列分别划分为N+1个数组段;指定滑动窗口的距离,使用滑动窗口技术在数组段内滑动,在滑动窗口内对每一段划分后的数组段分别计算基因得分;根据基因得分更新相对重排频率值数组和重排频率累加值数组并输出。本发明可以量化不同线粒体基因簇内的重排事件,从而显著地展示出不同类群线粒体基因组结构的差异性,具有较高的准确性。
Description
技术领域
本发明属于分子生物学和生物信息学技术领域,更具体地说,涉及一种优化的线粒体基因重排量化方法。
背景技术
动物线粒体基因组因其结构紧凑,被认为是研究基因组的窗口。从8000种后生动物的线粒体基因组的统计分析中发现,线粒体基因顺序受到选择限制,可能影响线粒体DNA调控,基因重排的发生频率和参与重排的基因簇(gene cluster)在不同类群中变化较大,还有个别类群具有特异性的基因顺序。正确地识别基因顺序以及评估相应变量,在揭示基因组的演化规律、探索基因组功能、研究分子演化过程等方面具有重要的科学价值,也为基因顺序变化的模式及其机制提供可靠的基础数据。
目前普遍被接受的线粒体基因重排机制是基于后口动物提出的“串联复制,随机丢失”(tandem duplication random loss,TDRL)模型,随着更多的线粒体基因重排被识别,重组(recombination)再次被认为或是重排的重要机制。
虽然先前的相关研究都取得了一定成果,但是也存在着一些不足:这些研究没有考虑到基因重排事件中的热点特异性,在量化过程中对于高频重排基因簇和保守基因簇没有给予不同权重,在重排复杂度较高的类群中表现欠佳,不利于基因组比较工作的开展。
专利公开号CN107885977A公布了一种用于检测动物类群线粒体基因组重排的方法,该方法是一种新的量化线粒体基因组重排的方法,它视单基因为一个结构单位,但计算其相邻的两个基因的变化。通过全局等价量化分类单元中的多个重排特征。除了对特定分类群的重排频率(RF)予以计分之外,还产出线粒体基因组中单个基因的重排分数(RS)以检测分类单元内的保守基因簇和高频重排基因簇。目前,倒置的发生率常被用作线粒体内重组发生的证据,但qMGR算法未能考虑倒置,这导致其低估了复杂度较高的基因重排。此外,该方法虽然能够检测出高频重排的基因簇,但并未进一步对保守基因簇和高频重排基因簇区分比较。
但qMGR因其算法的固有问题,导致其高估了重排基因两侧基因的重排频率,低估了单侧未发生重排的基因的重排频率;同时无法有效处理基因重复、假基因或多控制区等问题。也因此,该算法无法有效评价强差异的基因重排类型的重排分数。
发明内容
针对上述问题,本发明提供了一种优化的线粒体基因重排量化方法,允许对重排强度各异的基因簇给予不同权重,可以量化不同线粒体基因簇内的移位、倒置和基因洗牌等多种重排事件,从而显著地展示出不同类群线粒体基因组结构的差异性,在揭示基因组的演化规律、探索基因组功能、研究分子演化过程等方面具有重要的科学价值。
本发明通过下述技术方案来实现:一种优化的线粒体基因重排量化方法,包括以下步骤:
步骤S1.输入一条基准基因序列和待量化基因序列;
步骤S2. 待量化基因序列逐列与基准基因序列进行异或运算,统计待量化基因序列每个位置的基因与基准基因序列对应位置的差异程度,并将运算结果按差异程度由高到低存储到划分数组中;
步骤S3.取划分数组中差异程度最高的N个位置表示重排次数最高的N个基因,将基准基因序列和每个待量化基因序列分别划分为N+1个数组段;
步骤S4.指定滑动窗口的距离,使用滑动窗口技术在数组段内滑动,在滑动窗口内对每一段划分后的数组段分别计算得分;在滑动窗口滑动的过程中,在每个滑动窗口内寻找是否有与B[i]相同的基因S[j][i],如果存在,就按照如下规则进行计分,如果不存在则直接计得分为2分;计算得分的规则为:
;
式中,Xji表示第j条待量化基因序列中第i个基因的得分;B[i]表示基准基因序列第i个位置上的基因;B[i-1]表示基准基因序列第i-1个位置上的基因;B[i+1]表示基准基因序列第i+1个位置上的基因;S[j][i]表示第j条待量化基因序列中第i个位置的基因;S[j][i-1]表示第j条待量化基因序列中第i-1个位置的基因;S[j][i+1]表示第j条待量化基因序列中第i+1个位置的基因;
步骤S5.根据计算的Xji,按如下方式更新相对重排频率值数组BScore和重排频率累加值数组SScore并输出:
;
;
其中,BScore[i]表示待量化基因序列中第i个基因的相对重排频率值, SScore[j]表示第j条待量化基因序列的重排频率累加值;m表示待量化基因序列的数量,n为基准基因序列的基因数量。
具体地,步骤S2的运算公式如下:
;
其中,D[i]表示划分数组第i个位置上的基因,xor表示异或运算。
具体地,选择N为5,取划分数组中差异程度最高的a,b,c,d,e 共5个位置,将基准基因序列和每个待量化基因序列分别划分为6个数组段:
基准基因序列=B[1,a)+B[a,b)+B[b,c)+B[c,d)+B[d,e)+B[e,n];
第j条待量化基因序列=S[j][1,a)+S[j][a,b)+S[j][b,c)+ S[j][c,d)+ S[j][d,e)+ S[j][e,n];
其中,B[1,a)表示基准基因序列第1个基因到第a-1个基因的数组段,B[a,b)表示基准基因序列第a个基因到第b-1个基因的数组段,B[b,c)表示基准基因序列第b个基因到第c-1个基因的数组段,B[c,d)表示基准基因序列第c个基因到第d-1个基因的数组段,B[d,e) 表示基准基因序列第d个基因到第e-1个基因的数组段,B[e,n] 表示基准基因序列第e个基因到第n基因的数组段;S[j][1,a)表示第j条待量化基因序列第1个基因到第a-1个基因的数组段,S[j][a,b)表示第j条待量化基因序列第a个基因到第b-1个基因的数组段,S[j][b,c)表示第j条待量化基因序列第b个基因到第c-1个基因的数组段, S[j][c,d)表示第j条待量化基因序列第c个基因到第d-1个基因的数组段, S[j][d,e) 表示第j条待量化基因序列第d个基因到第e-1个基因的数组段, S[j][e,n] 表示第j条待量化基因序列第e个基因到第n基因的数组段。
具体地,利用Python语言对动物线粒体基因组中每一个基因的基因序列进行分析提取,获取每条基因组上的所有的基因的排列信息,以每个类群中典型的动物线粒体基因组排布为基准基因序列。
本发明通过划分数组来基因序列分段,并使用滑动窗口来计算每个基因的得分,可以量化不同线粒体基因簇内的移位、倒置和基因洗牌等多种重排事件,从而显著地展示出不同类群线粒体基因组结构的差异性,具有较高的准确性。通过对一组待量化的基因序列填充间隙来完成比对,以更大程度的在最大简约原则下复现这组数据的同源性。
附图说明
图1为滑动窗口分区示意图,图中字母A、B、C、D、E、F、G、H、I、J均为基因名称。
具体实施方式
下面结合实施例对本发明进一步详细说明。
从线粒体基因组公共数据库中获取不同类型物种的所有线粒体基因组信息,如在Genebank数据库(https://www .ncbi .nlm .nih .gov/nuccore)中下载某一特定动物线粒体基因组的全序列数据,然后对数据进行处理:利用编程语言Python语言分析提取相关数据,得到线粒体基因组中每一个基因的位置信息,进而获取每条基因组上的所有基因的位置排序,为后期线粒体基因重排做准备。
本发明的一种优化的线粒体基因重排量化方法,包括以下步骤:
步骤S1.输入一条基准基因序列和待量化基因序列;
步骤S2. 待量化基因序列逐列与基准基因序列进行异或运算,统计待量化基因序列每个位置的基因与基准基因序列对应位置的差异程度,并将运算结果按差异程度由高到低存储到划分数组中;
对每个位置上的基因,将待量化基因序列与基准基因序列相比较,与基准基因不同的个数记录下来,如基准基因序列中第五个位置上的基因为trnW,对第五列的所有基因进行纵向比较,将不是trnW的基因数量记录在划分数组的第五个位置上,划分数组的第五个位置上记录的基因数量就代表待量化基因序列第五个位置的基因与基准基因序列差异程度,对每一列都进行同样的操作,进而初步统计出每个位置基因的差异程度,并将运算结果按差异程度由高到低存储到划分数组中;
;
其中,D[i]表示划分数组第i个位置上的基因,基准基因序列是一维数组,B[i]表示基准基因序列第i个位置上的基因, 待量化基因序列是二维数组,m表示待量化基因序列的数量,S[j][i]表示第j条待量化基因序列中第i个位置的基因,xor表示异或运算。
步骤S3.取划分数组中差异程度最高的N个位置表示重排次数最高的N个基因,将基准基因序列和每个待量化基因序列分别划分为N+1个数组段。图1中N为2,划分为3个数组段,类似地,本实施例选择N为5,取划分数组中差异程度最高的5个位置(a,b,c,d,e),将基准基因序列和每个待量化基因序列分别划分为6个数组段。
基准基因序列=B[1,a)+B[a,b)+B[b,c)+B[c,d)+B[d,e)+B[e,n];
第j条待量化基因序列=S[j][1,a)+S[j][a,b)+S[j][b,c)+ S[j][c,d)+ S[j][d,e)+ S[j][e,n];
其中,B[1,a)表示基准基因序列第1个基因到第a-1个基因的数组段,B[a,b)表示基准基因序列第a个基因到第b-1个基因的数组段,B[b,c)表示基准基因序列第b个基因到第c-1个基因的数组段,B[c,d)表示基准基因序列第c个基因到第d-1个基因的数组段,B[d,e) 表示基准基因序列第d个基因到第e-1个基因的数组段,B[e,n] 表示基准基因序列第e个基因到第n基因的数组段,n为基准基因序列的基因数量;S[j][1,a)表示第j条待量化基因序列第1个基因到第a-1个基因的数组段,S[j][a,b)表示第j条待量化基因序列第a个基因到第b-1个基因的数组段,S[j][b,c)表示第j条待量化基因序列第b个基因到第c-1个基因的数组段, S[j][c,d)表示第j条待量化基因序列第c个基因到第d-1个基因的数组段, S[j][d,e) 表示第j条待量化基因序列第d个基因到第e-1个基因的数组段, S[j][e,n] 表示第j条待量化基因序列第e个基因到第n基因的数组段。
步骤S4.指定滑动窗口的距离,使用滑动窗口技术在数组段内滑动,在滑动窗口内对每一段划分后的数组段(B[1,a)、B[a,b)、B[b,c)、B[c,d)、B[d,e)、B[e,n]、S[j][1,a)等)分别计算得分,如指定滑动窗口距离为3,对数组段B[1,a)进行操作,则每一段滑动窗口分别为:B[1,3)、B[3,6)……B[a-3,a)。
在滑动窗口滑动的过程中,在每个滑动窗口内寻找是否有与B[i]相同的基因S[j][i],如果存在,就按照如下规则进行计分,如果不存在则直接计得分为2分;计算得分的规则为:
;
式中,Xji表示第j条待量化基因序列中第i个基因的得分;B[i-1]表示基准基因序列第i-1个位置上的基因;B[i+1]表示基准基因序列第i+1个位置上的基因; S[j][i-1]表示第j条待量化基因序列中第i-1个位置的基因;S[j][i+1]表示第j条待量化基因序列中第i+1个位置的基因;
得分的规则实际是:
相对于基准基因序列,如果该基因左右两侧相邻基因不变,则不给该基因加分;
如果左侧或是右侧相邻基因发生了改变,则给该基因加1分;
如果被考察基因出现缺失,重复,转置或是其两侧相邻的基因均发生了改变,则给该基因加2分。
步骤S5.根据计算的Xji,按如下方式更新相对重排频率值数组BScore和重排频率累加值数组SScore并输出:
;
;
其中,BScore[i]表示待量化基因序列中第i个基因的相对重排频率值, SScore[j]表示第j条待量化基因序列的重排频率累加值;m表示待量化基因序列的数量,n为基准基因序列的基因数量。
对类群中每条待量化基因序列中的每个基因进行上述方法扫描,从而就可以获得该类群线粒体基因组中各个基因的相对重排频率值和重排频率累加值。
依据本发明中所提出的方法,本发明进行了两栖类、鱼类、爬行类、鸟类、哺乳类动物线粒体基因重排的研究,实验步骤如下:
(1)重新运算已发表的脊椎动物线粒体基因重排比较研究中的源数据。
(2)利用Python语言就上述数据包进行清洗。
(3)利用本发明提出的方法,分别计算各分类单元的单个基因的重排频率。再对不同基因重排类型的量化分数进行分析研究。
(4)比较qMGR和本发明提出的方法的量化结果,获得的新发现如下:
(1)从NCBI网站下载两栖类、鱼类、爬行类、鸟类、哺乳类动物动物线粒体基因组全序列数据,共计2834个物种信息(2023年6月下载)。
(2)利用Python语言对上述动物线粒体基因组中每一个基因的基因序列进行分析提取,获取每条基因组上的所有的基因的排列信息,以每个类群中典型的动物线粒体基因组排布为基准基因序列。
(3)利用本发明提出的方法,计算出单个基因重排频率。再对各个物种的量化分数进行分析研究。
(4)在对两栖动物线粒体基因组组成中的37个基因和CR区(ControlRegion区)的基因重排情况进行量化计算,获得的新发现如下:
①去除离群值后,五个分类群的重排得分(RS)依次为两栖类 > 鱼类 > 哺乳类 >爬行类 > 鸟类;类群内差异度分化强度依次为哺乳类> 鱼类 > 两栖类 > 爬行类 > 鸟类。
②鱼类线粒体基因组的重排类型多样,保守性最低。重排得分(RS)均值在五个分类群中最高。重排频率(RF)最高的基因是trnQ、trnM和trnE,重排频率最高的蛋白质编码基因为nad6。
③两栖类线粒体基因组完全没有反转,但存在单个基因或基因簇的长距离位移。重排频率(RF)较高的基因有trnM、trnL1和trnP,蛋白质编码基因的重排频率均不高。
④爬行类线粒体基因组的重排类型主要包括反转、长距离位移、基因复制等。重排频率(RF)较高的基因有发生复制的trnL2,和重排高频发生的控制区上游若干基因。
⑤鸟类线粒体基因组仅在控制区上游段落发生基因复制,其余基因均未发生任何重排。但需要注意鸟类基因组存在大量源数据的注释错误,需要清洗数据后再重新量化。
⑥哺乳动物线粒体基因组仅在trnW至trnC这一基因簇中发生位移变化,反转基因仅局限于tRNA。
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该本发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
Claims (7)
1.一种优化的线粒体基因重排量化方法,其特征在于,包括以下步骤:
步骤S1.输入一条基准基因序列和待量化基因序列;
步骤S2. 待量化基因序列逐列与基准基因序列进行异或运算,统计待量化基因序列每个位置的基因与基准基因序列对应位置的差异程度,并将运算结果按差异程度由高到低存储到划分数组中;
步骤S3.取划分数组中差异程度最高的N个位置表示重排次数最高的N个基因,将基准基因序列和每个待量化基因序列分别划分为N+1个数组段;
步骤S4.指定滑动窗口的距离,使用滑动窗口技术在数组段内滑动,在滑动窗口内对每一段划分后的数组段分别计算得分;在滑动窗口滑动的过程中,在每个滑动窗口内寻找是否有与B[i]相同的基因S[j][i],如果存在,就按照如下规则进行计分,如果不存在则直接计得分为2分;计算得分的规则为:
;
式中,Xji表示第j条待量化基因序列中第i个基因的得分;B[i]表示基准基因序列第i个位置上的基因;B[i-1]表示基准基因序列第i-1个位置上的基因;B[i+1]表示基准基因序列第i+1个位置上的基因;S[j][i]表示第j条待量化基因序列中第i个位置的基因;S[j][i-1]表示第j条待量化基因序列中第i-1个位置的基因;S[j][i+1]表示第j条待量化基因序列中第i+1个位置的基因;
步骤S5.根据计算的Xji,按如下方式更新相对重排频率值数组BScore和重排频率累加值数组SScore并输出:
;
;
其中,BScore[i]表示待量化基因序列中第i个基因的相对重排频率值, SScore[j]表示第j条待量化基因序列的重排频率累加值;m表示待量化基因序列的数量,n为基准基因序列的基因数量。
2.根据权利要求1所述的一种优化的线粒体基因重排量化方法,其特征在于,步骤S2的运算公式如下:
;
其中,D[i]表示划分数组第i个位置上的基因,xor表示异或运算。
3.根据权利要求1所述的一种优化的线粒体基因重排量化方法,其特征在于,选择N为5,取划分数组中差异程度最高的a,b,c,d,e 共5个位置,将基准基因序列和每个待量化基因序列分别划分为3个数组段。
4.根据权利要求3所述的一种优化的线粒体基因重排量化方法,其特征在于,对基准基因序列和待量化基因序列的划分方式如下:
基准基因序列=B[1,a)+B[a,b)+B[b,c)+B[c,d)+B[d,e)+B[e,n];
第j条待量化基因序列=S[j][1,a)+S[j][a,b)+S[j][b,c)+ S[j][c,d)+ S[j][d,e)+S[j][e,n];
其中,B[1,a)表示基准基因序列第1个基因到第a-1个基因的数组段,B[a,b)表示基准基因序列第a个基因到第b-1个基因的数组段,B[b,c)表示基准基因序列第b个基因到第c-1个基因的数组段,B[c,d)表示基准基因序列第c个基因到第d-1个基因的数组段,B[d,e) 表示基准基因序列第d个基因到第e-1个基因的数组段,B[e,n] 表示基准基因序列第e个基因到第n基因的数组段;S[j][1,a)表示第j条待量化基因序列第1个基因到第a-1个基因的数组段,S[j][a,b)表示第j条待量化基因序列第a个基因到第b-1个基因的数组段,S[j][b,c)表示第j条待量化基因序列第b个基因到第c-1个基因的数组段, S[j][c,d)表示第j条待量化基因序列第c个基因到第d-1个基因的数组段, S[j][d,e) 表示第j条待量化基因序列第d个基因到第e-1个基因的数组段, S[j][e,n] 表示第j条待量化基因序列第e个基因到第n基因的数组段。
5.根据权利要求1所述的一种优化的线粒体基因重排量化方法,其特征在于,利用Python语言对动物线粒体基因组中每一个基因的基因序列进行分析提取,获取每条基因组上的所有的基因的排列信息,以每个类群中典型的动物线粒体基因组排布为基准基因序列。
6.根据权利要求1所述的一种优化的线粒体基因重排量化方法,其特征在于,步骤S3中选择N为2。
7.根据权利要求1所述的一种优化的线粒体基因重排量化方法,其特征在于,步骤S4中指定滑动窗口距离为3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311378396.9A CN117133351B (zh) | 2023-10-24 | 2023-10-24 | 一种优化的线粒体基因重排量化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311378396.9A CN117133351B (zh) | 2023-10-24 | 2023-10-24 | 一种优化的线粒体基因重排量化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117133351A CN117133351A (zh) | 2023-11-28 |
CN117133351B true CN117133351B (zh) | 2024-01-23 |
Family
ID=88863091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311378396.9A Active CN117133351B (zh) | 2023-10-24 | 2023-10-24 | 一种优化的线粒体基因重排量化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117133351B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2863274A1 (fr) * | 2003-12-05 | 2005-06-10 | Commissariat Energie Atomique | Procede d'evaluation quantitative d'un rearrangement ou d'une recombinaison genetique ciblee d'un individu et ses applications. |
CN104094120A (zh) * | 2011-12-08 | 2014-10-08 | 凡弗3基因组有限公司 | 包含mdm2的双微染色体及其方法 |
CN107885977A (zh) * | 2017-11-30 | 2018-04-06 | 淮南师范学院 | 一种用于检测动物类群线粒体基因组重排的方法 |
CN114395630A (zh) * | 2021-02-04 | 2022-04-26 | 湖南农业大学 | 基于高通量测序寄生虱线粒体基因组组装的方法和应用 |
-
2023
- 2023-10-24 CN CN202311378396.9A patent/CN117133351B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2863274A1 (fr) * | 2003-12-05 | 2005-06-10 | Commissariat Energie Atomique | Procede d'evaluation quantitative d'un rearrangement ou d'une recombinaison genetique ciblee d'un individu et ses applications. |
CN104094120A (zh) * | 2011-12-08 | 2014-10-08 | 凡弗3基因组有限公司 | 包含mdm2的双微染色体及其方法 |
CN107885977A (zh) * | 2017-11-30 | 2018-04-06 | 淮南师范学院 | 一种用于检测动物类群线粒体基因组重排的方法 |
CN114395630A (zh) * | 2021-02-04 | 2022-04-26 | 湖南农业大学 | 基于高通量测序寄生虱线粒体基因组组装的方法和应用 |
Non-Patent Citations (2)
Title |
---|
Dynamic comparisons of high-resolution expression profiles highlighting mitochondria-related genes between in vivo and in vitro fertilized early mouse embryos;Ren, LK等;《HUMAN REPRODUCTION》;全文 * |
原发性干燥综合征并发非霍奇金淋巴瘤患者唇腺组织BCL-6基因重排的检测及意义;杨敏;梁胜翔;李宝贞;石宇红;许佳;莫汉有;;中国当代医药(34);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117133351A (zh) | 2023-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Waring et al. | Automated machine learning: Review of the state-of-the-art and opportunities for healthcare | |
Hu et al. | A coherent pattern mining algorithm based on all contiguous column bicluster | |
US20230207061A1 (en) | Pathogenicity language model | |
Rashed et al. | Accelerating DNA pairwise sequence alignment using FPGA and a customized convolutional neural network | |
Kasap et al. | Design and implementation of an FPGA-based core for gapped BLAST sequence alignment with the two-hit method | |
Yoshida et al. | Tempo and mode in karyotype evolution revealed by a probabilistic model incorporating both chromosome number and morphology | |
CN117133351B (zh) | 一种优化的线粒体基因重排量化方法 | |
Hafemeister et al. | Classifying short gene expression time-courses with Bayesian estimation of piecewise constant functions | |
Church et al. | Phylotranscriptomics reveals discordance in the phylogeny of Hawaiian Drosophila and Scaptomyza (Diptera: Drosophilidae) | |
Mäkinen et al. | Genome-scale algorithm design: bioinformatics in the era of high-throughput sequencing | |
Du et al. | Deep multi-label joint learning for RNA and DNA-binding proteins prediction | |
Miralavy et al. | Using genetic programming to predict and optimize protein function | |
Dong et al. | Joint gene network construction by single-cell rna sequencing data | |
Rorabaugh et al. | High frequency accuracy and loss data of random neural networks trained on image datasets | |
CN115249538B (zh) | 一种基于异构图生成对抗网络的lncRNA-疾病关联预测模型的构建方法 | |
Nguyen et al. | Efficient agglomerative hierarchical clustering for biological sequence analysis | |
Rabea et al. | A fast algorithm for constructing suffix arrays for DNA alphabets | |
Ullrich et al. | oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data | |
Wang et al. | AFSBN: A method of artificial fish swarm optimizing bayesian network for epistasis detection | |
Lee et al. | Protein secondary structure prediction using BLAST and exhaustive RT-RICO, the search for optimal segment length and threshold | |
CN107885977A (zh) | 一种用于检测动物类群线粒体基因组重排的方法 | |
Khitmoh et al. | Stretch Profile: A pruning technique to accelerate DNA sequence search | |
Behera et al. | Evolutionary computation approach to enhance protein multiple sequence alignments | |
Wegmann | Technical and Algorithmic Optimization of PaPaRa | |
Haeusler | Prof. Marcus Vinicius Soledade Poggi de Aragao |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |