Characterization of the complete mitochondrial genome of Hydrotaea spinigera (Diptera: Muscidae) with phylogenetic implications

Hydrotaea spinigera Stein is a muscid species with high forensic importance. To further supplement the genome-level features of related species, the complete mitochondrial genome of H. spinigera is amplified, sequenced, annotated, analyzed, and compared with 13 other species of the family Muscidae. This mitogenome is 15,517?bp in length, with a standard set of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a typical non-coding control region. The gene organization, base composition, and secondary structures of tRNA genes of this mitogenome are similar to other dipteran species. Phylogenetic analyses based on all the mitochondrial PCGs and rRNA genes for 14 muscid and one outgroup taxa are performed using Bayesian inference (BI) and maximum likelihood (ML) analysis. The inferred trees indicate that the Muscidae, the subfamily Muscinae, and the tribes Azeliini and Reinwardtiini are monophyletic, whereas the monophyly of the subfamily Azeliinae and the tribe Muscini is not supported.     

The complete chloroplast genome sequence of the endangered species Syringa pinnatifolia (Oleaceae)

Syringa pinnatifolia is an endangered endemic species in China with important ornamental and medicinal value, and it needs urgent protection. Here, we report the complete chloroplast (cp) genome structure of S. pinnatifolia and its evolution is inferred through comparative studies with related species. The S. pinnatifolia cp genome was 155 326 bp and contained a large single copy region (LSC) of 86 167 bp and a small single copy region (SSC) of 17 775 bp, as well as a pair of inverted repeat regions (IRs) of 25 692 bp. A total of 113 unique genes were annotated, including 79 protein‐coding genes, 30 tRNA genes and four rRNA genes. The GC content of the S. pinnatifolia cp genome was 37.9%, and the corresponding values in the LSC, SSC and IR regions were 36.0, 32.1, 43.2% respectively. Repetitive sequences analysis revealed that the S. pinnatifolia cp genome contained 38 repeats. Microsatellite marker detection analysis identified 253 simple sequence repeats (SSRs), which provides opportunities for future studies of the population genetics and phylogenetic relationships of Syringa. Phylogenetic analysis of 29 selected cp genomes revealed that S. pinnatifolia is closely related to Syringa vulgaris and all 27 Lamiales species formed a clade separate from the two outgroup species. This newly characterized S. pinnatifolia chloroplast genome will provide a useful genomic resource of phylogenetic inference and the development of more genetic markers for species discrimination and population studies in the genus Syringa.     

Complete mitochondrial genome of Dacus conopsoides (Insecta: Tephritidae) with tRNA gene duplication and molecular phylogeny of Dacini tribe

To date there is only a single report on the complete mitochondrial genome of the Dacus fruit flies. We report here the whole mitogenome of Dacus conopsoides with first report of tRNA gene duplication in tephritid fruit flies determined using next-generation sequencing and discuss the molecular phylogeny of Dacini tribe. It had a total length of 15,852 bp, comprising 13 protein coding genes, 2 rRNA genes, 23 tRNA genes, and a non-coding region (A + T-rich control region). The 65-bp trnF gene was duplicated, and the 68-bp trnE gene was partially duplicated resulting in a 31-bp pseudogene. The cloverleaf structure for trnN, trnH, and trnF lacked the TΨC-loop, while trnS lacked the D-stem. The start codons for the protein coding genes included 6 ATG, 3 ATC, 2 ATA, and 1 each of ATT and TCG. Seven PCGs had TAA stop codon, two had TAG and four had incomplete T stop codon. Molecular phylogeny based on 15 mt-genes (13 PCGs +2 rRNA genes) and 30 taxa of Tephritidae indicated D. conopsoides forming a monophyletic sister group with D. longicornis supported by high bootstrap value. The lineage containing also the monophyletic genus Zeugodacus. The Dacini and Ceratitidini tribes of the subfamily Dacinae were monophyletic but the subfamilies Dacinae and Trypetinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

The complete mitochondrial genome of Rhynchocypris oxycephalus (Teleostei: Cyprinidae) and its phylogenetic implications

Rhynchocypris oxycephalus (Teleostei: Cyprinidae) is a typical small cold water fish, which is distributed widely and mainly inhabits in East Asia. Here, we sequenced and determined the complete mitochondrial genome of R. oxycephalus and studied its phylogenetic implication. R. oxycephalus mitogenome is 16,609 bp in length (GenBank accession no.: MH885043), and it contains 13 protein‐coding genes (PCGs), two rRNA genes, 22 tRNA genes, and two noncoding regions (the control region and the putative origin of light‐strand replication). 12 PCGs started with ATG, while COI used GTG as the start codon. The secondary structure of tRNA‐Ser (AGN) lacks the dihydrouracil (DHU) arm. The control region is 943bp in length, with a termination‐associated sequence, six conserved sequence blocks (CSB‐1, CSB‐2, CSB‐3, CSB‐D, CSB‐E, CSB‐F), and a repetitive sequence. Phylogenetic analysis was performed with maximum likelihood and Bayesian methods based on the concatenated nucleotide sequence of 13 PCGs and the complete sequence without control region, and the result revealed that the relationship between R. oxycephalus and R. percnurus is closest, while the relationship with R. kumgangensis is farthest. The genus Rhynchocypris is revealed as a polyphyletic group, and R. kumgangensis had distant relationship with other Rhynchocypris species. In addition, COI and ND2 genes are considered as the fittest DNA barcoding gene in genus Rhynchocypris. This work provides additional molecular information for studying R. oxycephalus conservation genetics and evolutionary relationships.     

Comparative Analysis of the Mitochondrial Genomes of Callitettixini Spittlebugs (Hemiptera: Cercopidae) Confirms the Overall High Evolutionary Speed of the AT-Rich Region but Reveals the Presence of Short Conservative Elements at the Tribal Level

The present study compares the mitochondrial genomes of five species of the spittlebug tribe Callitettixini (Hemiptera: Cercopoidea: Cercopidae) from eastern Asia. All genomes of the five species sequenced are circular double-stranded DNA molecules and range from 15,222 to 15,637 bp in length. They contain 22 tRNA genes, 13 protein coding genes (PCGs) and 2 rRNA genes and share the putative ancestral gene arrangement of insects. The PCGs show an extreme bias of nucleotide and amino acid composition. Significant differences of the substitution rates among the different genes as well as the different codon position of each PCG are revealed by the comparative evolutionary analyses. The substitution speeds of the first and second codon position of different PCGs are negatively correlated with their GC content. Among the five species, the AT-rich region features great differences in length and pattern and generally shows a 2–5 times higher substitution rate than the fastest PCG in the mitochondrial genome, atp8. Despite the significant variability in length, short conservative segments were identified in the AT-rich region within Callitettixini, although absent from the other groups of the spittlebug superfamily Cercopoidea.     

Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species

Magnolia grandiflora is an important medicinal,ornamental and horticultural plant species.The chloroplast(cp) genome of M.grandiflora was sequenced using a 454 sequencing platform and the genome structure was compared with other related species.The complete cp genome of M.grandiflora was 159623 bp in length and contained a pair of inverted repeats(IR) of 26563 bp separated by large and small single copy(LSC,SSC) regions of 87757 and 18740 bp,respectively.A total of 129 genes were successfully annotated,18 of which included introns.The identity,number and GC content of M.grandiflora cp genes were similar to those of other Magnoliaceae species genomes.Analysis revealed 218 simple sequence repeat(SSR) loci,most composed of A or T,contributing to a bias in base composition.The types and abundances of repeat units in Magnoliaceae species were relatively conserved and these loci will be useful for developing M.grandiflora cp genome vectors.In addition,results indicated that the cp genome size in Magnoliaceae species and the position of the IR border were closely related to the length of the ycf1 gene.Phylogenetic analyses based on 66 shared genes from 30 species using maximum parsimony(MP) and maximum likelihood(ML) methods provided strong support for the phylogenetic position of Magnolia.The availability of the complete cp genome sequence of M.grandiflora provides valuable information for breeding of desirable varieties,cp genetic engineering,developing useful molecular markers and phylogenetic analyses in Magnoliaceae.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

Complete sequence and organization of the cucumber (Cucumis sativus L. cv. Baekmibaekdadagi) chloroplast genome

The nucleotide sequence of the cucumber (Cucumis sativus L. cv. Baekmibaekdadagi) chloroplast genome was completed (DQ119058). The circular double-stranded DNA, consisting of 155,527 bp, contained a pair of inverted repeat regions (IRa and IRb) of 25,187 bp each, which were separated by small and large single copy regions of 86,879 and 18,274 bp, respectively. The presence and relative positions of 113 genes (76 peptide-encoding genes, 30 tRNA genes, four rRNA genes, and three conserved open reading frames) were identified. The major portion (55.76%) of the C. sativus chloroplast genome consisted of gene-coding regions (49.13% protein coding and 6.63% RNA regions; 27.81% LSC, 9.46% SSC and 18.49% IR regions), while intergenic spacers (including 20 introns) made up 44.24%. The overall G-C content of C. sativus chloroplast genome was 36.95%. Sixteen genes contained one intron, while two genes had two introns. The expansion/contraction manner of IR at IRb/LSC and IR/SSC border in Cucumis was similar to that of Lotus and Arabidopsis, and the manner at IRa/LSC was similar to Lotus and Nicotiana. In total, 56 simple sequence repeats (more than 10 bases) were identified in the C. sativus chloroplast genome.     

The complete structure of the cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.) chloroplast genome: Its composition and comparative analysis

The complete nucleotide sequence of the cucumber (C. sativus L. var. Borszczagowski) chloroplast genome has been determined. The genome is composed of 155,293 bp containing a pair of inverted repeats of 25,191 bp, which are separated by two single-copy regions, a small 18,222-bp one and a large 86,688-bp one. The chloroplast genome of cucumber contains 130 known genes, including 89 protein-coding genes, 8 ribosomal RNA genes (4 rRNA species), and 37 tRNA genes (30 tRNA species), with 18 of them located in the inverted repeat region. Of these genes, 16 contain one intron, and two genes and one ycf contain 2 introns. Twenty-one small inversions that form stem-loop structures, ranging from 18 to 49 bp, have been identified. Eight of them show similarity to those of other species, while eight seem to be cucumber specific. Detailed comparisons of ycf2 and ycf15, and the overall structure to other chloroplast genomes were performed.     

三种厚朴叶绿体基因组的比较研究

为了深入发掘日本厚朴、厚朴、凹叶厚朴叶绿体基因组差异,筛选厚朴优良性状候选基因,开展三种厚朴的分子遗传研究,该文利用Illumina HiSeq高通量测序平台首次对日本厚朴叶绿体进行测序、组装,并与已有的厚朴、凹叶厚朴叶绿体基因组共同注释,获得三个物种叶绿体基因图谱,筛选出三个基因组中的差异基因,又与同科中11个亲缘物种进行叶绿体基因组比对,构建NJ遗传树。结果表明:(1)日本厚朴叶绿体基因组的Clean Reads为19 791 019,Q30为91.33%,组装后基因组全长160 051 bp,GC含量为39.2%,含tRNA 37个,rRNA 8个。(2)比对分析发现三种厚朴具有相似的IR、LSC和SSC结构,以及GC含量和tRNA数量,但编码基因种类和数量、内含子和外显子的数量和结构等存在差异。(3)日本厚朴的功能基因数目较厚朴、凹叶厚朴分别多6个和4个,主要分布于LSC区和IR区,涉及核糖体大亚基、核糖体小亚基和未知功能基因类群。(4)系统发育分析结果进一步显示日本厚朴与凹叶厚朴亲缘关系较近,其次是厚朴。该研究表明日本厚朴具有更丰富的叶绿体基因组结构、组成和变异特征,是其适应高纬度地区弱光、低温环境的分子机制,这为厚朴类优良品种的分子选育提供有力的指导。     

De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an edible root vegetable crop that is cultivated worldwide and whose genome has been sequenced. Here we report the complete nucleotide sequence of the radish cultivar WK10039 chloroplast (cp) genome, along with a de novo assembly strategy using whole genome shotgun sequence reads obtained by next generation sequencing. The radish cp genome is 153,368 bp in length and has a typical quadripartite structure, composed of a pair of inverted repeat regions (26,217 bp each), a large single copy region (83,170 bp), and a small single copy region (17,764 bp). The radish cp genome contains 87 predicted protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Sequence analysis revealed the presence of 91 simple sequence repeats (SSRs) in the radish cp genome.     

The complete mitochondrial genome of Eriocheir japonica sinensis (Decapoda: Varunidae) and its phylogenetic analysis

The complete mitochondrial genome (mitogenome) can provide novel insights into understanding the mechanisms underlying mitogenome evolution. In this study, the complete mitogenome of Eriocheir japonica sinensis (Decapoda: Varunidae) was determined to be 16,378 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a D-loop region. The AT skew of the E. j. sinensis mitogenome was slightly negative (−0.016), indicating a higher number of T compared with A nucleotides. The nucleotide composition of the mitogenome was also biased toward A + T nucleotides (71.6%). All PCGs were initiated by ATN codons. Eight of the 13 PCGs harbored the incomplete termination codon by T, or TA. All other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The D-loop region of the E. j. sinensis mitogenome was 918 bp in length. Based on 13 PCGs, phylogenetic analysis confirmed the placement of E. j. sinensis within the Varunidae.     

The complete mitochondrial genome of <Emphasis Type="Italic">Damora sagana</Emphasis> and phylogenetic analyses of the family Nymphalidae

The monotypic genus Damora (Nymphalidae, Heliconiinae) contains a single species, Damora sagana, which is widely distributed across southern China. Herein, its complete mitogenome was sequenced to further understand lepidopteran mitogenome characteristics, reconstruct the nymphalid family phylogeny, and infer the subdivision of Heliconiinae species. The circular mitogenome was 15,151 bp long, abundant in A and T, and comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and one control region with a gene arrangement typical of lepidopteran mitogenomes. ATN codons initiated all PCGs, except cytochrome c oxidase subunit 1 (COX1), which was initiated by a CGA sequence as has been observed in other lepidopterans. Three PCGs (COX1, COX2 and ND4) employed a single T termination signal, whereas others had the typical complete termination codon (TAA). All tRNA genes were folded into the typical cloverleaf structure except for tRNA-Ser (AGN). The A+T-rich region included the conserved motif ‘ATAGA’ followed by a 17 bp poly-T stretch, which was also observed in tribe Argynnini mitogenomes. A phylogenetic tree was constructed via multiple methods using the 13 PCGs data of D. sagana and other available mitogenomes of nymphalid species. All three phylogenetic trees yielded the same topology. These results were consistent with those from previous studies of most major nymphalid groups, except those regarding tribe subdivision in certain subfamilies such as Argynnini + (Acraeini + Heliconiini) for Heliconiine. Furthermore, our analyses identified that the genus Cethosia was grouped with the genus Acraea composing the tribe Acraeini with strong support.     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).     

The complete chloroplast genome sequence of the Chinese endemic species Sorbus setschwanensis (Rosaceae) and its phylogenetic analysis

Sorbus setschwanensis Koehne is a pinnate-leaved Sorbus s.str. species endemic to China with narrow distribution and intriguing phylogeny that needs wider attention. In this paper, the complete chloroplast (cp) genome of S. setschwanensis is reported, and its phylogenetic position is analyzed. The complete cp genome of S. setschwanensis is 160 064 bp in size with 36.50% GC content. It has a typical quadripartite structure including a pair of inverted repeat regions (IRs) of 26 378 bp that separates a large single copy (LSC) region of 86 013 bp and a small single copy (SSC) region of 19 295 bp. The cp genome encodes 108 genes, comprising 76 protein-coding genes, 28 tRNA genes and 4 rRNA genes. Additionally, 52 simple sequence repeats (SSRs) and 43 dispersed repeats were identified. Comparison of the whole cp genome with those of other Sorbus species showed an overall high degree of sequence similarity, but there are six highly variable regions (trnR-atpA, petN-psbM, ndhC-trnV, trnE-trnT, trnT-trnL and rpl32-trnL) located in intergenic spacers that may be useful as molecular markers in future population genetic and phylogenetic studies in the genus. Phylogenetic analyses based on 108 coding genes from 25 species in Rosaceae revealed that S. setschwanensis is nested within Sorbus sect. Sorbus together with other pinnately leaved species, but does not form a sister lineage to S. rufopilosa belonging to the same series Multijugae. Thus, the systematic position of S. setschwanensis and relationships with other species in the genus needs to be further studied.     

Characterization and phylogenetic implication of complete mitochondrial genome of the Cataglyphis aenescens (Hymenoptera: Formicidae): Genomic comparisons in Formicinae

Ants are the most widely distributed social insects in terrestrial ecosystems with the largest number of species and quantities; their ecological function is very significant and they have certain medicinal value. Because mitochondria have their own DNA (mtDNA) and genetic system, so the complete mitochondrial genome has been used more and more in molecular phylogeny, population genetics, molecular diagnosis, and evolution of insects. In this study, the mitochondrial genome of Cataglyphis aenescens (Nylander, 1849) (Hymenoptera: Formicidae) was sequenced and annotated. The results showed that the total length of the mitochondrial genome of C. aenescens (Nylander, 1849) was 17,197 bp; it contained 37 genes, including 13 protein coding genes (PCGs), 22 tRNA genes (tRNAs), 2 rRNA genes (rRNAs) and a non-coding region. All PCGs of C. aenescens were initiated with ATN codons and terminated with the TAA codon. The overall nucleotide composition of C. aenescens was AT-biased (81.74%). Cataglyphis and Formica clustered together. C. aenescens is closely related to Formica sinae, Formica. fusca and Formica selysi. This work enhances the genetic data of Formicidae and contributes to our understanding of their phylogenic relationship, evolution, and utilization.     

Complete mitochondrial genome of the wild silkmoth,Saturnia boisduvalii (Lepidoptera: Saturniidae)

We sequenced mitochondrial genome (mitogenome) of the wild silkmoth, Saturnia boisduvalii (Lepidoptera: Saturniidae), which occurs in mainland Korea, and compared it with other species in Bombycoidea to characterize the genomic evolution of the superfamily. We found that the composition and arrangement of genes in the 15,257‐bp S. boisduvalii genome are typical of the majority of Lepidoptera, and the genome is biased toward A/T nucleotides, as previously reported. Comparison of individual gene divergence among bombycoid species showed that ND6 was most variable (p‐distance = 0.21), whereas COI and COII were most conserved, indicating that of all the protein‐coding genes (PCGs) ND6 appear to have evolved most rapidly. Thus, other PCGs beside COI are potential alternative markers, where scrutinized discrimination among species is required.     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

青藏高原特有种-藏扇穗茅叶绿体基因组测序及序列分析

藏扇穗茅（Littledalea tibetica）是禾本科（Poaceae）雀麦族（Bromeae）中一个具有重要生态价值的多年生高山特有种,主要分布于青藏高原及其毗邻地区。本文采用基于第二代高通量测序平台的Illumina MiSeq技术,对青藏高原特有种—藏扇穗茅进行了叶绿体基因组测序,首次建立了雀麦族物种的标准测序流程;同时,以其近缘物种—黑麦草（Lolium perenne）的叶绿体基因组序列作为参考,组装获得它的叶绿体基因组序列。结果表明,藏扇穗茅叶绿体基因组序列全长136 852 bp,GC含量为38.5%,呈典型的四段式结构,其中大（LSC）、小（SSC）单拷贝区大小分别为80 970和12 876 bp,反向互补重复区（IR）大小为21 503 bp,共注释得到141个基因,包含95个蛋白编码基因、38个tRNA基因和8个rRNA基因,主要分布于大单拷贝区和小单拷贝区。同时,基于藏扇穗茅和其它30种禾本科植物叶绿体基因全序列构建的系统发育树显示,藏扇穗茅与早熟禾亚科中小麦族植物亲缘关系较近。     

The first complete mitochondrial genome of the Indian Tent Turtle,Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

The characterization of a complete mitogenome is widely used in genomics studies for systematics and evolutionary research. However, the sequences and structural motifs contained within the mitogenome of Testudines taxa have rarely been examined. The present study decodes the first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (16,657 bp) by using next‐generation sequencing. This denovo assembly encodes 37 genes: 13 protein‐coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA, and one control region (CR). Most of the genes were encoded on majority strand, except for one PCG (NADH dehydrogenase subunit 6) and eight tRNAs. Most of the PCGs were started with an ATG initiation codon, except for Cytochrome oxidase subunit 1 with “GTG” and NADH dehydrogenase subunit 5 with “ATA.” The termination codons, “TAA” and “AGA” were observed in two subunits of NADH dehydrogenase gene. The relative synonymous codon usage analysis revealed the maximum abundance of alanine, isoleucine, leucine, and threonine. The nonsynonymous/synonymous ratios were <1 in all PCGs, which indicates strong negative selection among all Geoemydid species. The study also found the typical cloverleaf secondary structure in most of the tRNA genes, except for serine with the lack of the conventional DHU arm. The comparative study of Geoemydid mitogenomes revealed the occurrence of tandem repeats was frequent in the 3′ end of CR. Further, two copies of a unique tandem repeat “TTCTCTTT” were identified in P. tentoria. The Bayesian and maximum‐likelihood phylogenetic trees using concatenation of 13 PCGs revealed the close relationships of P. tentoria with Batagur trivittata in the studied dataset. All the Geoemydid species showed distinct clustering with high bootstrap support congruent with previous evolutionary hypotheses. We suggest that the generations of more mitogenomes of Geoemydid species are required, to improve our understanding of their in‐depth phylogenetic and evolutionary relationships.