Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous–Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India–Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.     

Complete mitochondrial genome sequence of Bonasa sewerzowi(Galliformes: Phasianidae) and phylogenetic analysis

Bonasa sewerzowi, the smallest and most southerly distributed grouse species in the world, is a bird endemic to China. The population of B. sewerzowi had shown a declining trend, which made it to be the endangered species in the China Red Data Book and Category I of nationally protected animals. So far, however, most studies about this species were mainly focused on the morphological and ecological aspects. In order to further study the feature of B. sewerzowi, the complete mitochondrial genome(mitogenome) of B. sewerzowi was sequenced by Illumina Hiseq 2000 high-throughput sequencing. Then, we focused on comparative genomics of two Bonasa species to find their characteristics. Finally, phylogenetic position of Bonasa was made based on the mitogenome dataset. Our results revealed that:(1) the mitogenome of B. sewerzowi, consisting of 16 658 bp, displayed typical genome organization and gene order found in other previously determined Galliformes mitogenomes;(2) the structure and composition of mitogenomes were similar between B. sewerzowi and B. bonasia;(3) the monophyly of Bonasa was well supported, which had a closer phylogenetic relationship with Meleagris gallopavo.     

访花昆虫野蚜蝇线粒体基因组结构分析

野蚜蝇Syrphus torvus(Osten-Sacken)成虫具访花习性,是重要的传粉昆虫;幼虫捕食蚜虫,是蚜虫重要天敌之一.本研究通过高通量测序、组装拼接获得了完整的野蚜蝇线粒体基因组全序列(GenBank登录号:MW074962),其序列全长为16 444 bp,包含13个蛋白编码基因(PCGs)、22个tRN...     

Molecular characterization of the complete mitochondrial genome sequence of Indian wild pig (Sus scrofa cristatus)

The Indian wild pig is a sub-species (Sus scrofa cristatus) which is different from the other pig breeds and is protected under Schedule-III of the Indian Wildlife (Protection) Act, 1972. In this study, complete mitogenome of two Indian wild pigs was sequenced and characterized by shotgun sequencing and de novo assembly, which revealed sequence size of 16,738 and 16,251?bp, respectively, (Accession no. MG725630 and MG725631). The mitogenome sequence in this study displayed 98% homology with previously reported mitogenome of pigs from different parts of the world. Mitogenome analysis by MITOS Web server revealed similarity of gene organization with the other vertebrates (13 protein-coding, 22 tRNAs, 2 rRNAs genes, and a control region). The mitogenomic sequences of Indian wild pig maintained a separate clade in the phylogenetic tree constructed by using 62 whole mitogenome sequences across the world. The phylogeny derived from mitogenomic sequences revealed distinct separate European–American and Asiatic pig clades. It was concluded that whole mitogenome sequencing using NGS without designing mitogenome-specific primer for amplification, is possible thereby reducing the cost and labor. This study is the first report of complete sequence of mitogenome of Indian wild pig.     

甜果螨全线粒体基因组测序与分析

【目的】测定和分析甜果螨Carpoglyphus lactis线粒体基因组全序列,并在线粒体基因组水平探讨其在真螨总目(Acariformes)中的系统发育地位,为真螨总目分类及果螨科线粒体基因组研究提供科学依据。【方法】挑取实验室饲养的甜果螨成螨,用传统的酚氯仿抽提法和试剂盒提取法提取甜果螨基因组DNA。然后采用节肢动物或螨类线粒体基因的通用引物PCR扩增出甜果螨线粒体基因cox1,cob,rrnS和nad4-nad5的部分序列;再设计种特异性引物进行Long-PCR扩增和步移法测序,测出甜果螨线粒体基因组全序列。应用SeqMan, SEQUIN 9.0和tRNAscan等生物信息学软件,对甜果螨线粒体基因组的基因结构等进行生物信息学分析。最后基于17种真螨总目螨类的蛋白质编码基因,采用最大似然法构建系统发育树。【结果】甜果螨线粒体全基因组总长为14 060 bp(GenBank登录号:MN073839),为典型的闭合双链DNA分子,共由37个基因组成,包括13个蛋白质编码基因(PCGs)、22个tRNA基因和2个rRNA基因;甜果螨线粒体基因组还包括1个大的非编码区(large n...     

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).     

Characterization of the first mitochondrial genome of a catenulid flatworm: Stenostomum leucops (Platyhelminthes)

The first complete mitochondrial genome of a catenulid, Stenostomum leucops, was characterized. Illumina sequencing and 90 813 reads were utilized in the assembly, producing one contig with an average coverage of 1118×. The length of this genome is 15 742 bp with 12 protein‐coding, two rRNA and 22 tRNA genes. Although the atp8 gene is absent in other Platyhelminthes, a highly divergent putative atp8 gene was found in S. leucops. In contrast to other Platyhelminthes, the mitochondrial genes of S. leucops are encoded on both strands. The gene order in the S. leucops mitogenome is very divergent from those observed in other Platyhelminthes, showing only small blocks of synteny. With AAA as the codon for lysine S. leucops shows the probable plesiomorphic condition, whereas Rhabditophora possess the derived GAA. This evolutionary transition is correlated with changes in the respective anticodons in trnK. It remains unclear whether the absence of the D arm and loop in trnS1 is a convergence in Catenulida and Neodermata.     

Grayling (Thymallinae) phylogeny within salmonids: complete mitochondrial DNA sequences of Thymallus arcticus and Thymallus thymallus

The phylogenetic relationships among the three subfamilies (Salmoninae, Coregoninae and Thymallinae) in the Salmonidae have not been addressed extensively at the molecular level. In this study, the whole mitochondrial genomes of two Thymallinae species, Thymallus arcticus and Thymallus thymallus were sequenced, and the published mitochondrial genome sequences of other salmonids were used for Bayesian and maximum‐likelihood phylogenetic analyses. These results support an ancestral Coregoninae, branching within the Salmonidae, with Thymallinae as the sister group to Salmoninae.     

Complete mitochondrial genome sequence of the Indian pipistrelle Pipistrellus coromandra (Vespertilioninae)

We characterized the complete mitogenome of Pipistrellus coromandra (Indian pipistrelle) for comparative analysis of mitogenomes and for resolving the phylogenetic relationship of four tribes in the subfamily Vespertilioninae. The mitogenome size of P. coromandra was 17,153?bp, with a control region and a typical set of 37 mitochondrial genes. The nucleotide composition of the P. coromandra mitogenome showed an AT bias with a nucleotide composition of 33.5% A, 30.7% T, 13.3% G, and 22.5% C. The mitochondrial protein-coding genes in P. coromandra use the standard start codon (ATN), two stop codons (TAA and AGA), and two incomplete stop codons (TA- and T--). The intertribal relationship of four tribes was highly resolved from the phylogenetic analysis of mitogenome sequences.     

Draft genome of the American Eel (Anguilla rostrata)

Freshwater eels (Anguilla sp.) have large economic, cultural, ecological and aesthetic importance worldwide, but they suffered more than 90% decline in global stocks over the past few decades. Proper genetic resources, such as sequenced, assembled and annotated genomes, are essential to help plan sustainable recoveries by identifying physiological, biochemical and genetic mechanisms that caused the declines or that may lead to recoveries. Here, we present the first sequenced genome of the American eel. This genome contained 305 043 contigs (N50 = 7397) and 79 209 scaffolds (N50 = 86 641) for a total size of 1.41 Gb, which is in the middle of the range of previous estimations for this species. In addition, protein‐coding regions, including introns and flanking regions, are very well represented in the genome, as 95.2% of the 458 core eukaryotic genes and 98.8% of the 248 ultra‐conserved subset were represented in the assembly and a total of 26 564 genes were annotated for future functional genomics studies. We performed a candidate gene analysis to compare three genes among all three freshwater eel species and, congruent with the phylogenetic relationships, Japanese eel (A. japanica) exhibited the most divergence. Overall, the sequenced genome presented in this study is a crucial addition to the presently available genetic tools to help guide future conservation efforts of freshwater eels.     

The complete mitochondrial genome of Pardosa pusiola (Araneae,Lycosidae) and its phylogenetic implications

The mitochondrial genome (mitogenome) is useful for identification and phylogenetic analyses among arthropods, but there are no sufficient mitogenome data for wolf spiders. To enrich the mitogenome database of wolf spiders, the complete mitogenome of Pardosa pusiola was sequenced by high-throughput sequencing. It is 14,284 bp, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region (CR). It represents a high bias toward A and T nucleotides with an A + T content of 76.49%. The mitogenome exhibited a negative AT skew (−0.13) and a positive GC skew (0.32). Most PCGs started with ATN codons and ended with TAA, TAG, or an incomplete T. In addition, most tRNAs had aberrant secondary structures with the absence of DHU arm or TΨC arm. Analysis performed with CREx software demonstrated that large-scale rearrangements of tRNAs were observed in the mitogenome of P. pusiola as compared with the putative ancestral mitogenome. The Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees based on the 13 PCGs of 25 spiders had the same topology, which could be presented as (Araneidae + (Agelenidae + (Dictynidae + Desidae)) + (Salticidae + (Thomisidae + (Oxyopidae + (Pisauridae + Lycosidae))))). This study offers a useful genetic resource for the taxonomy and phylogeny of spiders.     

变灰青霉线粒体基因组特征及系统发育分析

本研究对一株分离自细虫草上的变灰青霉SFY00C3菌株线粒体基因组进行测定,分析其组成特征,并探究其与青霉属真菌的系统发育关系。结果表明,SFY00C3的线粒体基因组是一条长度为28 301 bp的环状DNA分子,共编码42个基因(15个蛋白编码基因、2个rRNA基因和25个tRNA基因),其碱基组成有显著的A-T偏好性,25个tRNA基因均可形成典型三叶草结构,并存在32处G-U错配。通过青霉属物种间共线性分析发现其线粒体基因组发生了基因重排;共有的14个蛋白编码基因的Ka/Ks值均小于1,表明受到了纯化选择压力的影响;系统发育分析表明：SFY00C3在青霉属中是一个独立的分支,应该是6种青霉祖先的姊妹。本研究丰富了变灰青霉的线粒体基因组序列信息,为青霉属的系统发育、资源保护及遗传多样性研究提供基础数据。     

Mitochondrial phylogeography of the Japanese pond turtle,Mauremys japonica (Testudines,Geoemydidae)

The phylogeography of the Japanese pond turtle, Mauremys japonica (Temminck and Schlegel, 1835), which is an endemic species in Japan, was studied by analyzing the variation in two mitochondrial DNA sequences, the cytochrome b gene and the control region. These analyses suggest that M. japonica comprises two major groups. The first one was found in the more eastern region, eastern Honshu Island and Shikoku Island, while the second was found in a western region, Kyushu Island and the Chugoku District (the westernmost part of Honshu Island). The boundary between the two groups is located in the Chugoku District. The nucleotide and haplotype diversities were very low, and these low diversities seem to have been caused by a bottleneck in the last glacial age. These results suggest that this species survived the last glacial period in two refugia, one in the central part of Honshu Island and the other one in Kyushu Island. Subsequently, population expansion took place in the postglacial period, and the groups from the two refugia extended their distribution ranges to the present boundary in the Chugoku District which represents a secondary contact zone.     

The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications

We determined the nucleotide sequence of the complete mitochondrial genome of the nematode species Anisakis simplex. The genome is circular, 13,916 bp in size and conforms to the general characteristics of nematode mitochondrial DNAs. The gene arrangement of A. simplex is the same as that of Ascaris suum and almost identical to those of rhabditid species with a minor exception concerning the relative position of the AT-rich and non-coding regions and radically different from those of spirurid species. Along with comparisons of gene arrangement, phylogenetic analyses (maximum parsimony, neighbour joining and maximum likelihood methods) based on concatenated amino acid sequences of 12 protein-coding genes from 13 nematode species provided strong support for the sister-group relationship between Ascaridida and Rhabditida. The Shimodaira-Hasegawa and Templeton's tests both rejected the alternative hypothesis of a closer relationship between Ascaridida and Spirurida. These results contradicted the traditional view of nematode classification and a recent molecular phylogenetic study of 18S rDNA data that assigned Ascaridida and Spirurida as being a sister-group. Mapping of gene arrangement across the phylogenetic tree lead to the assumption that the conserved gene arrangement found in Ascaridida-Rhabditida members might have been acquired after the most recent common ancestor of ascaridid/rhabditid members branched off from the basal stock of the rhabditid lineage.     

The complete mitochondrial genome of Plesiastrea versipora (Scleractinia,Plesiastreidae) sheds light on its phylogeny and taxonomy of the family Plesiastreidae

The genus Plesiastrea used to be a member of the traditional family Faviidae, falling into the challenging ‘Bigmessidae’ clade, and was re-established until recent molecular phylogenies published. The entire mitogenome of the symbiotic coral Plesiastrea versipora (Lamarck, 1816), the type species of the family Plesiastreidae, was sequenced. The length of the mitochondrial genome is 15,320 bp and it includes thirteen protein-coding genes (PCGs), two rRNAs and two tRNAs. The nucleotide composition of GC is 32%. We perform phylogenetic reconstruction based on maximum likelihood (ML) and Bayesian analysis(BI) using all PCGs. Our result indicates that P. versipora clusters closely with species which belong to Mussidae, Merulinidae and Lobophylliidae. Our phylogenetic analyses provide solid evidence for phylogenetic placement of P. versipora and the evolutionary relationships among different families within the traditional robust clade of Scleractinia. In addition, the mitogenome data provide useful information for further molecular systematic investigations on Plesiastreidae as well as conservation biology research of P. versipora.     

The mitogenome of freshwater loach Homatula laxiclathra (Teleostei: Nemacheilidae) with phylogenetic analysis of Nemacheilidae

The complete mitogenome can provide valuable genetic information to reconstruct relationships between species. In this study, we sequenced a stone loach, Homatula laxiclathra (Teleostei: Nemacheilidae), which is found in the northern region of the Qinling Mountains in China. The size of the H. laxiclathra mitogenome is 16,570 bp, which contains 37 typical mitochondrial genes including 13 protein‐coding genes, 22 transfer RNAs, two ribosomal RNAs, and a control region (D‐loop) with a total AT content of 55.8%. This is similar to other Nemacheilidae sequences published in GenBank. Furthermore, a mito‐phylogenomic analysis of 46 Nemacheilidae species places H. laxiclathra in a robust monophyletic Homatula cluster with other Homatula species. Our results contribute toward a better understanding of a true phylogeny of these species based on large‐scale taxonomic samplings as well as to help grasp the evolution of fish mitogenomes.     

Evolutionary implications of analyses of complete mitochondrial genomes across order Zoantharia (Cnidaria: Hexacorallia)

Cnidarians are early-diverging metazoans, but evolutionary aspects of some taxa are still poorly understood, as in the order Zoantharia (Anthozoa: Hexacorallia). Zoantharians have been divided into two suborders based on the arrangement of the fifth septae as complete (Macrocnemina) or incomplete (Brachycnemina). Previous molecular phylogenetic analyses have indicated the need for re-evaluation as Macrocnemina has been found to be paraphyletic. Despite many phylogenetic studies, the recovery of complete mitochondrial genomes (mt-genomes) for systematic and evolutionary studies of zoantharians has been limited. The present study represents the first to sequence the complete mt-genomes of members of eight of nine zoantharian families. Although all examined mt-genomes had the same gene order arrangement, there were variations among mt-genomes' sizes, nucleotide substitution rates, and introns. Only two species did not have the cox1 intron, which harbors a gene coding a homing endonuclease of the LAGLIDADG type. Our mitogenomic analyses also showed relatively high nucleotide diversity in mt-DNA regions other than the standard regions traditionally considered for DNA barcoding of this group. Phylogenetic analyses using 13 mt-genome protein-coding genes recovered a fully resolved tree with clear separation between macrocnemic representatives. Ancestral state reconstruction analyses revealed three main transitions in arrangement of the marginal musculature through the evolutionary history of the order. An “early” transition from reticulate mesogleal to a cteniform endodermal arrangement was followed by transitions that occurred in the common ancestor of the Brachycnemina and family Hydrozoanthidae. Our results indicate the need for clarification of higher-level phylogeny and taxonomy of Zoantharia.