麦鱼线粒体基因组序列测定与结构特征分析

通过PCR扩增和一代测序技术获得麦鱼(Rhinogobius sp.)线粒体基因组全序列，并对其序列结构进行了分析。结果表明：1)麦鱼线粒体基因组全长16 511 bp,其碱基G+C的含量为47.7%,A+T的含量为52.3%,呈现AT偏好性；2)共含有蛋白质编码基因13个，rRNA基因2个，tRNA基因22个和D-loop区，除ND6、tRNA-Gln、tRNA-Tyr、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Ser、tRNA-Glu、tRNA-Pro基因在L链上编码之外，其他基因均在H链上进行编码；3)基于线粒体基因组全序列和Cytb基因，分别构建了鰕虎鱼类中部分属的进化树，其表明麦鱼与波氏吻鰕虎鱼、褐吻鰕虎鱼、子陵吻鰕虎鱼聚在一个分支上，均属于吻鰕虎鱼属，并且其与波氏吻鰕虎鱼的亲缘关系最近。     

黄带拟鲹线粒体基因组测序及鲹科鱼类系统发育分析

摘要: 为深入开展黄带拟鲹种质鉴定、分类及遗传进化等研究,通过二代高通量测序与生物信息学分析,揭示了黄带拟鲹(Pseudocaranx dentex)线粒体基因组全序列基因结构及鲹科鱼类系统发育特征。结果显示,黄带拟鰺线粒体基因组为典型的环状DNA结构,序列全长为16 570 bp,碱基组成分别为A(27.2 %)、G(17.18%)、C(30.24%)和T(25.38%),包括13种蛋白编码基因,22个tRNA基因,2个rRNA基因,除ND6、tRNAGln、tRNAAla、tRNAAsn、tRNACys、tRNATyr、tRNASer、tRNAGlu、tRNAPro外,其余基因均在H链上编码,且基因组存在多处重叠区域。除CO Ⅰ和ND5的起始密码子分别为ATC和ATA外,其余11个蛋白编码基因的起始密码子均为ATG,以典型的TAA和TAG为终止密码子,在Cyt b中存在不完全终止密码子T;黄带拟鲹线粒体基因组全序列与蛋白编码基因的A+T含量分别为52.58%和51.55%,非编码控制区(D-loop)A+T富含61.69%,具有明显的AT偏好性。除tRNASer-(GCT)外,其余21个tRNA均含典型三叶草二级结构。为进一步研究黄带拟鲹系统进化特性,通过与18种鲹科鱼类线粒体基因组全序列及16S rRNA基因构建系统进化树显示,黄带拟鲹与高体若鲹同属一个分支,亲缘关系最近。结果可为黄带拟鲹物种鉴别、遗传多样性评价与保护提供技术依据。     

鲂属鱼类线粒体基因组的比较及其系统发育分析

基于GenBank中团头鲂线粒体基因组全序列和三角鲂、厚颌鲂、广东鲂的部分线粒体基因组序列,设计引物扩增出三角鲂、厚颌鲂和广东鲂3种鱼线粒体基因组全序列,同时对4种鲂属鱼类线粒体基因组全序列进行了比较分析。结果表明,4种鲂属鱼类线粒体基因组基因排列顺序完全相同,排列紧密,均包含13个蛋白质编码基因、22个tRNA、2个rRNA、1个非编码控制区(D-loop区)和1个轻链复制起始区(OL区)。除ND6和8个tRNA在L链上编码外,其余的基因均在H链上编码。4种鲂属线粒体基因组13个蛋白质编码基因中,均呈现出较强的A+T偏向性和C碱基偏好。全序列比对结果显示,共有758个变异位点,其中非简约性信息位点有691个,占总变异位点的91.16%,简约性信息位点有67个,仅占总变异位点的8.84%。22个tRNA基因中只有11个存在种间变异,共23个变异位点,主要发生在tRNA三叶草结构的TΨC和DHU臂环上。13个蛋白质编码基因中共检测出626个变异位点,这些变异主要发生在密码子第三位,占总变异位点的82.59%,其中变异位点数最多的是Cyt b基因,达84个,其次是ND 4基因(83个)。因此,Cyt b和ND4基因可作为备选的分子标记,用于鲂属群体间的遗传学研究。基于4种鲂属鱼类线粒体基因组全序列构建的ML树和BI树均显示,三角鲂与厚颌鲂的亲缘关系最近,团头鲂与它们的亲缘关系相对较近,而广东鲂与前述3种鲂属鱼类的亲缘关系均较远。     

鳀科鱼类线粒体全基因组序列结构特征及系统发育信息分析

为了解鳀科鱼类的线粒体全基因组序列结构特征及系统发育信息,以期为进化遗传学研究和分子标记的选取提供参考依据,对已知的10种鳀科鱼类的线粒体全基因组进行分析。结果显示:1)鳀科线粒体基因组全序列长度在16 660 bp到17 069 bp之间,基因组的结构和基因排列顺序与其它硬骨鱼类一致。2)比对后获得一致序列长度为15 704 bp(不含D-loop),其中变异位点5 570个,占所有位点数的35.5%。在编码基因中,序列变异程度和Kimura双参数遗传距离最大的是ND6基因(分别是47.5%和0.276),最小的是tRNA拼接序列(分别为18.7%和0.072)。3)基于Ka-Ks的Z检验和Tajima’s D检验表明蛋白质编码基因主要受到净化选择(即负选择,purifying selection)的作用;其中12个蛋白质编码基因(ND6除外)有强烈的净化选择(Ka/Ks<1),而ND6基因受正选择(positive selection)影响较大(Ka/Ks>1)。4)ND4、ND2和Cytb是进行鳀科鱼类系统发育分析的较为理想的分子标记。     

瘤背石磺线粒体基因组全序列分析

采用LA-PCR技术对瘤背石磺线粒体基因组全序列进行了测定和分析。结果表明,瘤背石磺线粒体基因组序列全长13 957 bp,由22个tRNA、2个rRNA、13个蛋白编码基因和19个长度为2~138 bp的非编码区组成。4个蛋白质编码基因和8个tRNA基因从L链编码,其余基因均从H链编码。蛋白质基因的起始密码子,除ND2为TTG以外,均为典型的起始密码子ATN。COⅢ和Cytb基因使用了不完全终止密码子T,其余基因均使用典型的TAA或TAG。预测了22个tRNA基因的二级结构,发现tRNA^Ser缺少DHU臂,tRNA^Ser和tRNA^Thr的反密码子环上有9个碱基,而不是通常的7个碱基。最长的非编码区含有类似于tRNA的二级结构。基于线粒体基因组编码的13个蛋白质的氨基酸序列,用NJ、MP、ME和UPGMA法构建系统进化树。分析6种软体动物之间的亲缘关系,结果与传统的系统分类基本一致。研究初步确定瘤背石磺与平疣桑椹石磺的亲缘关系比与凯尔特石磺的亲缘关系近。     

拟腹吸鳅属鱼类线粒体基因组特征与系统发育分析

为了解拟腹吸鳅属(Pseudogastromyzon)鱼类的系统发育特征,对已获得的11种该属鱼类进行线粒体基因组的测序与分析,结果表明,其线粒体基因组的长度在16 560～16 574 bp之间,且基因组成和排列模式与大多数脊椎动物线粒体基因组基本一致;各线粒体基因组序列中A+T的含量均超过50%,并存在反G偏倚现象。通过对11种拟腹吸鳅属鱼类线粒体基因组13个蛋白质编码基因与控制区序列的比较发现：在COⅠ基因中有一段6 bp的碱基插入和缺失;D-loop基因的进化速率介于蛋白质编码基因之间,ND4、ND5、Cyt b和D-loop等基因适合作为本属鱼类种间系统发育研究的分子标记。基于COⅠ基因序列的11种拟腹吸鳅属鱼类种间遗传距离以及基于线粒体基因组序列构建的11个种系统发育树研究表明,九龙江拟腹吸鳅(P.fasciatus jiulongjiangensis)与梅花山拟腹吸鳅(P.meihuashanensis)应为拟腹吸鳅(P.fasciatus fasciatus)的同物异名。9个有效种可以分为两大类：颏部吸附器为品字型的3种属于品唇鳅亚属(Labigastromyzon),...     

基于线粒体全基因组结构的鲳属鱼类分子分类研究

为明确鲳属鱼类的线粒体全基因组序列结构、寻找有效的种间鉴别分子标记,并剖析其系统发育信息,从而为鲳属鱼类的分类学、进化遗传学和种质鉴定提供参考依据,测定了鲳属5种鱼类(中国鲳Pampus chinensis、灰鲳P.cinereus、银鲳P.argenteus、珍鲳P.minor、翎鲳P.punctatissimus)的线粒体全基因组序列,并结合NCBI数据库中已有的鲳属鱼类线粒体全基因组序列进行分析。结果显示:1)鲳属线粒体基因组全序列长度在16551 bp到18062 bp之间,其基因组结构与大多数硬骨鱼类一致;比较特别的是,银鲳和镰鲳具有两个tRNA;结构,珍鲳出现ND1基因重排以及两个重复的D-loop结构。2)银鲳与镰鲳的D-loop区出现重复CSB3结构,灰鲳与翎鲳D-loop区出现重复TAS结构。3)ND2基因条形码及其微型条形码均可将鲳属鱼类区分开,可作为有效的鲳属鉴别分子标记。4)东海区的银鲳与镰鲳应为同一物种,灰鲳可能为翎鲳的同物异名。     

滁州鲫(Carassius auratus)线粒体全基因组序列分析及系统进化

为探讨天然三倍体滁州鲫的系统进化地位,采用直接测序法获得滁州鲫线粒体基因组。其序列全长为16581 bp,碱基组成为31.6% A、26.2% T、16.1% G和26.1% C,包括13个蛋白质基因、22个tRNA基因、2个rRNA基因和1个非编码区,各基因的位置及组成与已公布的鲤科鱼类一致。除tRNA-Ser (AGY)外,其他21个tRNA的二级结构均具有典型的三叶草结构;13个蛋白编码基因中,除COⅠ起始密码子为GTG外,其余均以ATG为起始密码子;COⅡ、ND3、ND4和Cytb基因的终止密码子为不完整的T,其他9个基因均具有完整的终止密码子TAA或TAG。序列分析表明,滁州鲫与其他鲫属鱼类(方正银鲫A系和D系、鲫、淇河鲫、萍乡肉红鲫、黑鲫、日本白鲫和日本银鲫)在线粒体基因组上均具有较高的序列同源性(>94%)。以鲤(Cyprinus carpio haematopterus)为外类群,基于线粒体13个蛋白质基因的核苷酸与氨基酸序列构建上述鲫属鱼类的系统进化树,结果显示,滁州鲫与方正银鲫亲缘关系最近,与黑鲫最远。综合以上研究结果,认为滁州鲫应为银鲫亚种的一个地方种群。     

笛鲷属线粒体基因组编码序列的系统进化分析能力评估

运用通用引物长距PCR(Long-PCR)和常规PCR相结合的方法测定了鲈形目笛鲷科的4种笛鲷属鱼类(孟加拉笛鲷、四带笛鲷、千年笛鲷和马拉巴笛鲷)和军曹鱼科的军曹鱼线粒体DNA基因组全序列(GenBank序列号分别为FJ171339,FJ416614,FJ824741,FJ824742 和NC_011219),得出所用全序列测定体系的方法通用性较强,操作简单。线粒体基因组的比对分析表明,测定的mtDNA基因组的绝大部分区段与GenBank中现有的脊椎动物的序列有较高的同源性。以军曹鱼外群结合GenBank中近缘笛鲷鱼类(勒氏笛鲷、蓝点笛鲷和黑带鳞鳍梅鲷)进行的聚类分析中,勒氏笛鲷与黑带鳞鳍梅鲷的聚类关系近于同属物种,与形态学分类存在矛盾。通过单个基因和基因拼接序列比较,综合考虑邻位连接法构建系统进化树的置信度和序列的信息量,对13种蛋白质编码基因在属内种间的系统进化分析能力进行了评估,将基因分成不同的4等:很好的为ATPase6和cox2;好的序列为ND2和cox1;差的为ND6、ND3和ATPase8;包括Cyt b在内的其余6种基因为中等。分析还揭示出序列长度的增加可以提高系统进化树的置信度,且属内物种间比较时序列长度的影响小于高级阶元。     

Cyt b和12S rRNA基因条形码在灯笼鱼科鱼类物种鉴定中的应用

灯笼鱼科鱼类种类繁多, 且同属鱼类形态学相近, 因此利用分子标记对灯笼鱼进行准确的物种鉴定具有重要价值。为探讨线粒体细胞色素 b 基因(Cyt b)和 12S rRNA 基因在灯笼鱼科物种鉴定中的适用性, 对西北太平洋采集的 56 尾灯笼鱼进行扩增, 并进行序列对比与系统发育分析。研究表明, 采集的样本包括 6 种灯笼鱼, 分别为瓦氏角灯鱼(Ceratoscopelus warmingii)、长体标灯鱼(Symbolophorus californiensis)、粗鳞灯笼鱼(Myctophum asperum)、 细泰勒灯鱼(Tarletonbeania crenularis)、日本背灯鱼(Notoscopelus japonicus)以及某背灯鱼属鱼类(Notoscopelus sp.)。 核苷酸多态性分析显示, 基于 Cyt b 基因的种内与种间遗传距离比基于 12S rRNA 基因的更大。比较灯笼鱼科 2 种基因序列的结构特征, 发现 Cyt b 基因的种间平均遗传距离是种内平均遗传距离的 25 倍, 12S rRNA 基因的种间平均遗传距离是种内平均遗传距离的 26 倍, 均符合作为 DNA 条形码的基本要求。系统进化分析显示, 每种灯笼鱼均能形成独立分支, 2 个基因均能对 6 种灯笼鱼类进行鉴别; 但在 Cyt b 基因构建的进化树中, 每种鱼类能更好与数据库中已有的序列进行聚类。综上所述, Cyt b 和 12S rRNA 作为 DNA 条形码可以有效地对灯笼鱼科鱼类物种进行鉴定, 且 Cyt b 基因在系统进化关系的研究上具有更高的适用性。     

青鱼线粒体ND4L的SNP筛选及在两种体色群体中的分布

为了分析青鱼线粒体蛋白编码基因与体色性状的相关性,从而筛选出与体色相关的SNP分子标记。本研究选取了5尾广东佛山灰色青鱼和5尾扬州邗江正常体色黑色青鱼群体的皮肤组织,通过qPCR检测13个线粒体蛋白质编码基因在两种颜色群体皮肤组织中的表达量,并进行显著性分析。对表达差异最明显的蛋白编码基因CDS区序列设计引物,并选取78尾佛山灰色青鱼样本和92尾邗江黑色青鱼样本进行测序,根据测序峰图筛选出SNP位点。结果发现：13个线粒体蛋白质编码基因中COII基因表达量最高,ND5基因表达量最低,ND4L基因在两种体色青鱼中的表达差异最为明显。在对170尾青鱼样本的线粒体ND4L基因上共检测到了2个SNP位点,位于编码区,且都为同义突变。卡方检验表明,252bp(C/T)处的SNP位点在两种体色青鱼群体中的等位基因频率和基因型频率都有极显著的差异(p＜0.01),243bp(A/G)处的SNP位点在两种体色青鱼群体中的等位基因频率和基因型频率差异不显著(p＞0.05),因此线粒体ND4L基因中的SNP位点C252T与青鱼体色显著相关。本研究中成功找到一个与青鱼体色相关的SNP分子标记,可以在青鱼的养殖和育种中提供理论分子层面的指导和帮助。     

黄条鰤线粒体全基因组测序及结构特征分析

通过二代测序、软件拼接获得中国黄海海域黄条鰤(Seriola aureovittata)线粒体基因组全序列。其序列全长为16609 bp,碱基组成分别为A(26.68%)、G(17.84%)、C(30.12%)和T(25.36%);共有13条蛋白编码基因, 22个tRNA基因, 2个rRNA基因,除NAD6、trnQ、trnA、trnN、trnC、trnY、trnS、trnE、trnP外,其余基因均在H链上编码。黄条鰤线粒体基因组全序列与蛋白编码基因的A+T含量分别为52.05%和51.085%,具有明显的AT偏好性。线粒体基因中存在2个散在重复序列,分别位于NAD1基因序列正义链的中游和COX2基因序列反义链的上游。在其22个tRNA基因中,除了tRNAGly外,均具有典型的三叶草二级结构。黄条鰤线粒体基因组的蛋白编码基因起止位点与密码子除COX1、NAD5外,均与日本海域出产的黄条鰤(Seriola lalandi)完全吻合,且COX1、NAD5基因皆短于日本黄条鰤;两者间存在一定的遗传差异。基于18种隶属于13属的鲹科鱼类线粒体基因组全序列构建系统发生树,可知小甘鲹(Seriolina nigrofasciata)、黄条鰤、五条鰤(Seriola quinqueradiata)、高体鰤(Seriola dumerili)、长鳍鰤(Seriola rivoliana)同属一近支,且黄条鰤与五条鰤亲缘关系最近,与小甘鲹进化距离较远。     

仿刺参线粒体全基因组序列结构及比较研究

利用已构建的仿刺参cDNA文库得到的线粒体DNA(mtDNA)相关基因序列设计扩增引物,测定了大连仿刺参线粒体基因组全序列,并对其进行了基因构成和进化分析。仿刺参线粒体基因组序列长16 109bp,其基因构成与其他后口动物基本一致,包括37个基因(2个rRNA基因、22个tRNA基因和13个蛋白质编码基因)和3个主要的非编码区。在其37个基因中,ND6、tRNASer(AGN)、tRNAGln、tRNAAla、tRNAVal、TrnaAsp位于L链上,其余均位于H链上。在13个蛋白质编码基因中,除ND1的起始密码子为GTG外,其余均以ATG作为起始密码子;除Cytb以"T"作为终止密码子外,其他蛋白质基因均具有完全的终止密码子,且在已知的棘皮动物线粒体蛋白质基因中,部分基因的起始和终止密码子表现出一定的纲内特异性。比较分析了大连、青岛、威海仿刺参线粒体基因组,三者的基因组成和排列相同,碱基组成相近,蛋白质编码基因的起始和终止密码子完全一致,但存在核苷酸和氨基酸序列的差异。三者的控制区序列存在多个插入/缺失和SNP位点。根据COI、Cytb和ND4计算了三者之间的遗传距离为0.006～0.018,遗传距离分析和系统进化关系分析都显示青岛仿刺参和威海仿刺参的关系较大连仿刺参更近。     

Complete nucleotide sequence and variation of mitochondrial DNA from 10 individuals of walleye pollock, Theragra chalcogramma

ABSTRACT:   The complete mitochondrial genome sequence of 10 walleye pollocks, Theragra chalcogramma , from the Japan Sea and Bering Sea was determined. The 16 568–16 571 bp genome contains the same 37 mitochondrial structural genes (two ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes) as found in all other vertebrates analyzed, in an organization identical to that of other bony fish. The major non-coding region had several conserved sequence features. Nucleotide variations of ND1, ND5, and control region were high, and these regions appear to be good candidates for high-resolution markers in population studies.     

Decadal change in abundance of surface migratory myctophid fishes in the Kuroshio region from 1957 to 1994

Decadal change in abundance of surface migratory myctophid fishes was examined in the Kuroshio region of the western North Pacific for the 35 yr from 1957 to 1994 in relation to the Kuroshio regime shift, indicated by changes in its transport volume and sea surface temperature in the early 1970s. Each year, samples were collected from the 0‐ to 1‐m depth layer at night, from January to March, at 83 stations on average. Ten species belonging to five genera occurred. Juveniles [standard length (SL) ≤ 40 mm] and adults (SL > 40 mm) of each species were analysed. In juveniles, Myctophum asperum, M. nitidulum, Symbolophorus evermanni and Centrobranchus brevirostris accounted for 91.1% of the total catch (7419 individuals), while M. nitidulum contributed 82.4% of the total adult catch (494 individuals). The abundance of juvenile S. evermanni, C. brevirostris and M. nitidulum showed prominent peaks in 1971, together with adult M. nitidulum. These peaks corresponded to extraordinarily high zooplankton abundance in the Kuroshio region in 1971, suggesting high recruitment success and low adult mortality in surface migratory myctophids. Smaller changes in abundance in other years were species‐specific, and the fluctuation patterns were mainly explained by year‐to‐year fluctuations in the strength of the Kuroshio and the geographical distribution of each species, i.e. an increased Kuroshio flow volume resulted in an increase in tropical species and a decrease in temperate species, but no change in the pan‐Kuroshio species.     

长江流域4个野生大眼鳜群体的遗传多样性分析

为探讨长江流域大眼鳜(Siniperca kneri)野生群体的遗传多样性和遗传结构,为长江大眼鳜野生群体资源的有效管理和合理开发利用提供基础科学数据支撑,本研究基于线粒体细胞色素 b (cytochrome b, Cyt b)基因及控制区(D-loop)全序列,对 4 个群体(赤水河群体-CS、南京群体-NJ、岳阳群体-YY、宜昌群体-YC)共计 79 尾个体进行了分析。结果如下:(1)获得了长 1141 bp 的 Cyt b 基因全序列,共检测到 26 种单倍型。单倍型多样性(haplotype diversity,Hd)指数为 0.685 (CS)~0.924 (YC),核苷酸多样性(nucleotide diversity,π)指数为 0.176%(CS)~0.285%(YY)。群体内与群体间的遗传距离均在 0.002~0.003。(2)获得了长度为 834~840 bp 的 D-loop 全序列,在 79 尾个体中共检测到46 种单倍型。单倍型多样性指数为 0.754 (CS)~0.990 (NJ),核苷酸多样性指数为 0.548%(CS)~1.412%(YC)。群体内遗传距离在 0.005 (CS)~0.014 (YC),群体间遗传距离在 0.008~0.012,提示 4 个野生群体均尚未达到亚种分化水平。分子方差分析(AMOVA)显示群体内的变异是总变异的主要来源,长江流域大眼鳜野生群体无显著遗传结构差异。     

Complete nucleotide sequences of mitochondrial DNA of long-finned squid Loligo edulis

The complete nucleotide sequence of the mitochondrial genome of long-finned squid Loligo edulis forma kensaki (L. edulis f. kensaki, Teuthida: Myopsida) was determined. It contains 13 protein-coding genes, 12SrRNA, 16SrRNA, and the 22 tRNA genes, and three long noncoding regions. We also determined the mitochondrial genome nucleotide sequence of seasonal brood samples of L. edulis f. budo collected from the western part of the Sea of Japan in the autumn to winter seasons, which showed morphologic discrepancies with L. edulis f. kensaki. Comparison of the two complete sets of mitochondrial DNA data showed that the similarity of the two sequences was as high as 99.9?%, with only eight nucleotide substitutions in coding regions. The complete mitochondrial genome of Loligo edulis f. kensaki from Nagasaki Pref. is 17360?bp and L. edulis f. budo from the Shimane Pref. is 17351?bp. The nucleotide sequences and haplotype patterns in ND1 suggest that L. edulis f. kensaki and L. edulis f. budo share a single genetic background.     

Evaluation of a Fast Method Based on the Presence of Two Restriction Sites in the Mitochondrial ND5 (mt ND5) Gene for the Identification of Scomber Species

The purpose of this work was to evaluate the suitability of a method based on the presence of two restriction sites (for Hae III and Hindf I) in the mitochondrial NADH dehydrogenase subunit 5 (mt ND5) gene to identify Scomber species. The evaluation was performed on 144 reference and market samples by sequencing of the entire 505-bp fragment of the mt ND5 gene and of a 464-bp fragment of the Kocher fragment of the cytochrome b gene (mt Cytb). Sequence analysis of any of the two fragments allows the identification of each of the four Scomber species, but S. japonicus and S. colias had the same restriction sites at the ND5 amplicon and would not have been differentiated by this analysis. Similarly, loss of the Hae III site in some S. scombrus individuals would have misidentified them as not being Scomber. All the market products were correctly labeled except one acquired in Spain labeled as originating in the Atlantic and containing S. japonicus.     

Complete mitochondrial DNA sequence of the Japanese flying fish <Emphasis Type="Italic">Cypselurus hiraii</Emphasis>

ABSTRACT:   In this study, to develop a technique that enables authentication of processed seafood, the complete nucleotide sequence of the mitochondrial genome for the Japanese flying fish Cypselurus hiraii was determined. Three segments spanning the entire genome were amplified using polymerase chain reaction, and products were subsequently used as templates for direct sequencing with 60 primers. The genome (16 528 base pairs) was found to contain the same 37 genes (two ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes) as those found in other vertebrate mitochondrial genomes, with the gene order being identical to that typical of vertebrates. A major non-coding region between the tRNA^Pro and tRNA^Phe genes (868 base pairs) appears to be the control (D-loop) region, as it has several conserved blocks characteristic of control regions.