Complete mitochondrial genomes of three sea basses Lateolabrax (Perciformes,Lateolabracidae) species: Genome description and phylogenetic considerations

Genus Lateolabrax consists of three species, Japanese sea bass Lateolabrax japonicus, spotted sea bass Lateolabrax maculatus and blackfin sea bass Lateolabrax latus. The complete mitochondrial DNA (mtDNA) of the three sea basses were amplified and sequenced to characterize and discuss their phylogenetic relationships. The length of mitogenomes was 16,593 bp, 16,479 bp and 16,600 bp, respectively, and all of them consisted of 13 protein-coding genes, 2 ribosomal RNA (rRNA), 22 transfer RNA (tRNA) and a control region, which are typical for mtDNA of vertebrate. Most genes were encoded on the H-strand, except for the ND6 and eight tRNA genes encoding on the L-strand. A significant variation among the three species was detected in length of the control region. Phylogenetic relationship among the three species was constructed based on the datasets, including the 12 protein-coding genes (except ND6 gene), 22 tRNA and 2 rRNA sequences. The results supported the sister taxon between L. japonicus and L. maculatus. The genetic resources reported here are useful for further studies in taxonomy and phylogeny of the three sea basses and related species.     

Molecular phylogenetics and biogeography of Lepus in Eastern Asia based on mitochondrial DNA sequences

In spite of several classification attempts among taxa of the genus Lepus, phylogenetic relationships still remain poorly understood. Here, we present molecular genetic evidence that may resolve some of the current incongruities in the phylogeny of the leporids. The complete mitochondrial cytb, 12S genes, and parts of ND4 and control region fragments were sequenced to examine phylogenetic relationships among Chinese hare taxa and other leporids throughout the World using maximum parsimony, maximum likelihood, and Bayesian phylogenetic reconstruction approaches. Using reconstructed phylogenies, we observed that the Chinese hare is not a single monophyletic group as originally thought. Instead, the data infers that the genus Lepus is monophyletic with three unique species groups: North American, Eurasian, and African. Ancestral area analysis indicated that ancestral Lepus arose in North America and then dispersed into Eurasia via the Bering Land Bridge eventually extending to Africa. Brooks Parsimony analysis showed that dispersal events followed by subsequent speciation have occurred in other geographic areas as well and resulted in the rapid radiation and speciation of Lepus. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches estimated the divergence time between the three major groups within Lepus. The genus appears to have arisen approximately 10.76 MYA (+/-0.86 MYA), with most speciation events occurring during the Pliocene epoch (5.65+/-1.15 MYA approximately 1.12 +/- 0.47 MYA).     

Phylogenetic relationships of the North American sturgeons (order Acipenseriformes) based on mitochondrial DNA sequences

The evolutionary relationships of the extant species within the order Acipenseriformes are not well understood. Nucleotide sequences of four mitochondrial genes (12S rRNA, COII, tRNA(Phe), and tRNA(Asp) genes) in North American sturgeon and paddlefish were examined to reconstruct a phylogeny. Analysis of the combined gene sequences suggests a basal placement of the paddlefish with regard to the sturgeons. Nucleotide sequences of all four genes for the three Scaphirhynchus species were identical. The position of Scaphirhynchus based on our data was uncertain. Within the genus Acipenser, the two Acipenser oxyrinchus subspecies were very similar in sequence and found to be basal to the remaining Acipenser species examined. Based on our data, Acipenser transmontanus and Acipenser medirostris were sister taxa, as were Acipenser fulvescens and Acipenser brevirostrum. Comparison of our results with hypotheses of sturgeon relationships proposed by previous authors is presented. The sequence data presented here are phylogenetically useful and provide a solid foundation of genetic information for the North American Acipenseriformes that can be expanded to include Eurasian species to provide a global picture of sturgeon evolution.     

Mutation rates and intraspecific divergence of the mitochondrial genome of Pristionchus pacificus

Evolutionary reconstruction of the natural history of an organism ultimately requires knowledge about the development, population genetics, ecology, and phylogeny of the species. Such investigations would benefit from studies of mutational processes because mutations are the source of natural variation. The nematode Pristionchus pacificus has been developed as a model organism in evolutionary biology by comparing its development with Caenorhabditis elegans. Pristionchus pacificus and related species are associated with scarab beetles, and their ecology and phylogeny are well known. More than 200 P. pacificus isolates from all over the world are available for this cosmopolitan species. We generated mutation accumulation (MA) lines in P. pacificus to study spontaneous mutation rates in the mitochondrial genome and compared mutation rate estimates with natural variation between nine representative isolates of the species. The P. pacificus mitochondrial genome is 15,955 bp in length and is typical for nematodes. Pristionchus pacificus has all known mitochondrial genes and contains an unusual suppressor transfer RNA (tRNA) for the codon UAA. This has most likely influenced the spectrum of observable mutations because 6 of 12 mutations found in the 82 MA lines analyzed are nonsense mutations that can be suppressed by the suppressor tRNA. The overall mutation rate in P. pacificus is 7.6 × 10?? per site per generation and is less than one order of magnitude different from estimates in C. elegans and Drosophila. Using this mutation rate estimate in a comparison of the mitochondrial genome of nine P. pacificus isolates, we calculate the minimum time to the most recent common ancestor at 10?-10? generations. The combination of mutation rate analysis with intraspecific divergence provides a powerful tool for the reconstruction of the natural history of P. pacificus, and we discuss the ecological implication of these findings.     

Mitochondrial tRNA mutations in Chinese children with tic disorders

Aim: To conduct the clinical, genetic, and molecular characterization of 494 Han Chinese subjects with tic disorders (TD).Methods: In the present study, we performed the mutational analysis of 22 mitochondrial tRNA genes in a large cohort of 494 Han Chinese subjects with TD via Sanger sequencing. These variants were then assessed for their pathogenic potential via phylogenetic, functional, and structural analyses.Results: A total of 73 tRNA gene variants (49 known and 24 novel) on 22 tRNA genes were identified. Among these, 18 tRNA variants that were absent or present in <1% of 485 Chinese control patient samples were localized to highly conserved nucleotides, or changed the modified nucleotides, and had the potential structural to alter tRNA structure and function. These variants were thus considered to be TD-associated mutations. In total, 25 subjects carried one of these 18 putative TD-associated tRNA variants with the total prevalence of 4.96%.Limitations: The phenotypic variability and incomplete penetrance of tic disorders in pedigrees carrying these tRNA mutations suggested the involvement of modifier factors, such as nuclear encoded genes associated mitochondrion, mitochondrial haplotypes, epigenetic, and environmental factors.Conclusion: Our data provide the evidence that mitochondrial tRNA mutations are the important causes of tic disorders among Chinese population. These findings also advance current understanding regarding the clinical relevance of tRNA mutations, and will guide future studies aimed at elucidating the pathophysiology of maternal tic disorders.     

The complete mitochondrial genome of Eterusia aedea (Lepidoptera,Zygaenidae) and comparison with other zygaenid moths

Zygaenidae comprises >1036 species, including many folivorous pests in agriculture. In the present study, the complete mitochondrial genome (mitogenome) of a major pest of tea trees, Eterusia aedea was determined. The 15,196-bp circular genome contained the common set of 37 mitochondrial genes (including 13 protein-coding genes, two rRNA genes, and 22 tRNA genes) and exhibited the similar genomic features to reported Zygaenidae mitogenome. Comparative analyses of Zygaenidae mitogenomes showed a typical evolutionary trend of lepidopteran mitogenomes. In addition, we also investigated the gene order of lepidopteran mitogenomes and proposed that the novel gene order trnA-trnR-trnN-trnE-trnS-trnF from Zygaenidae and Gelechiidae and most other gene rearrangements of this tRNA cluster evolved independently. Finally, the mitogenomic phylogeny of Lepidoptera was reconstructed based on multiple mitochondrial datasets. And all the phylogenetic results revealed the sister relationships of Cossoidea and Zygaenoidea with both BI and ML methods, which is the first stable mitogenomic evidence for this clade.     

Linking the mitochondrial genotype to the organismal phenotype

One of the grand challenges of the postgenomics era is to mechanistically link the genotype with the phenotype. Here, we consider the link between the mitochondrial genotype and the organismal phenotype that is provided by bioenergetic studies of the electron transport chain. That linkage is pertinent for the fields of molecular ecology and phylogeography as it tests if, and potentially how, natural selection can influence the evolutionary and demographic past of both populations and species. We introduce the mitochondrial genotype in terms of mitochondrial‐encoded genes, nuclear‐encoded genes that produce structural proteins imported into the mitochondria, and mitochondrial DNA–nuclear interactions. We then review the potential for quaternary structure modelling to predict the functional consequence of specific naturally occurring mutations. We discuss how the energy‐producing reactions of oxidative phosphorylation can be used to provide a mechanistic biochemical link between genotype and phenotype. Experimental manipulations can then be used to test the functional consequences of specific mutations in multiple genetic backgrounds. Finally, we examine how mitochondria can influence the organismal mitochondrial phenotype using the examples of lifespan, fertility and starvation resistance and discuss how mitochondria may be involved in establishing both the upper and lower thermal limits of organisms. We conclude that mitochondrial DNA mutations can be important in determining aspects of organism life history. The question that remains to be resolved is how common are these adaptive mutations?     

Use of RAD sequencing for delimiting species

RAD-tag sequencing is a promising method for conducting genome-wide evolutionary studies. However, to date, only a handful of studies empirically tested its applicability above the species level. In this communication, we use RAD tags to contribute to the delimitation of species within a diverse genus of deep-sea octocorals, Chrysogorgia, for which few classical genetic markers have proved informative. Previous studies have hypothesized that single mitochondrial haplotypes can be used to delimit Chrysogorgia species. On the basis of two lanes of Illumina sequencing, we inferred phylogenetic relationships among 12 putative species that were delimited using mitochondrial data, comparing two RAD analysis pipelines (Stacks and PyRAD). The number of homologous RAD loci decreased dramatically with increasing divergence, as >70% of loci are lost when comparing specimens separated by two mutations on the 700-nt long mitochondrial phylogeny. Species delimitation hypotheses based on the mitochondrial mtMutS gene are largely supported, as six out of nine putative species represented by more than one colony were recovered as discrete, well-supported clades. Significant genetic structure (correlating with geography) was detected within one putative species, suggesting that individuals characterized by the same mtMutS haplotype may belong to distinct species. Conversely, three mtMutS haplotypes formed one well-supported clade within which no population structure was detected, also suggesting that intraspecific variation exists at mtMutS in Chrysogorgia. Despite an impressive decrease in the number of homologous loci across clades, RAD data helped us to fine-tune our interpretations of classical mitochondrial markers used in octocoral species delimitation, and discover previously undetected diversity.     

Mitochondrial tRNA cleavage by tRNA-targeting ribonuclease causes mitochondrial dysfunction observed in mitochondrial disease

Mitochondrial DNA (mtDNA) is a genome possessed by mitochondria. Since reactive oxygen species (ROS) are generated during aerobic respiration in mitochondria, mtDNA is commonly exposed to the risk of DNA damage. Mitochondrial disease is caused by mitochondrial dysfunction, and mutations or deletions on mitochondrial tRNA (mt tRNA) genes are often observed in mtDNA of patients with the disease. Hence, the correlation between mt tRNA activity and mitochondrial dysfunction has been assessed. Then, cybrid cells, which are constructed by the fusion of an enucleated cell harboring altered mtDNA with a ρ⁰ cell, have long been used for the analysis due to difficulty in mtDNA manipulation. Here, we propose a new method that involves mt tRNA cleavage by a bacterial tRNA-specific ribonuclease. The ribonuclease tagged with a mitochondrial-targeting sequence (MTS) was successfully translocated to the mitochondrial matrix. Additionally, mt tRNA cleavage, which resulted in the decrease of cytochrome c oxidase (COX) activity, was observed.     

基于线粒体基因组序列的鞘翅目肉食亚目水生类群系统发育分析

【目的】明确肉食亚目(Adephaga)水生类群线粒体基因组的基本特征,并基于线粒体基因组序列分析肉食亚目水生类群的系统发育关系。【方法】基于Illumina HiSeq X Ten测序技术测定了圆鞘隐盾豉甲Dineutus mellyi和齿缘龙虱Eretes sticticus的线粒体全基因组序列,对其进行了基因注释,并对其tRNA基因二级结构进行了预测分析。加上已公布的鞘翅目(Coleoptera)肉食亚目水生类群17个种的线粒体基因组序列,对该类群共19个种线粒体的蛋白质编码基因(protein-coding genes, PCGs)开展了比较基因组学分析,包括AT含量、密码子偏好性、选择压力等。基于13个PCGs的氨基酸序列和核苷酸序列,利用最大似然法(ML)和贝叶斯法(BI)分别构建鞘翅目肉食亚目水生类群的系统发育关系,并通过FcLM分析进一步评估伪龙虱科(Noteridae)和瀑甲科(Meruidae)的系统发育位置。【结果】圆鞘隐盾豉甲和齿缘龙虱的线粒体基因组全长分别为16 123 bp(GenBank登录号: MN781126)和16 196 bp(GenBank登录号: MN781132),都包含13个PCGs、22个tRNA基因、2个rRNA基因和1个D-loop区(控制区)。19个肉食亚目水生类群线粒体基因组PCGs的碱基组成都呈现A+T偏好性,在密码子使用上也都偏向于使用富含A+T的密码子;在进化过程中13个PCGs的进化模式相同,都受到纯化选择。基于线粒体基因组13个PCGs的氨基酸序列的肉食亚目水生类群的系统发育关系为(豉甲科Gyrinidae+(沼梭甲科Haliplidae+((壁甲科Aspidytidae+(两栖甲科Amphizoidae+龙虱科Dytiscidae))+(水甲科Hygrobiidae+(瀑甲科Meruidae+伪龙虱科Noteridae)))))。【结论】研究结果表明,豉甲科是肉食亚目水生类群的基部类群,接下来是沼梭甲科和龙虱总科;伪龙虱科和瀑甲科互为姐妹群,并一起作为龙虱总科内部的一个分支;两栖甲科与龙虱科具有更近的亲缘关系。     

The complete mitochondrial genome of Caprella scaura (Crustacea, Amphipoda, Caprellidea), with emphasis on the unique gene order pattern and duplicated control region

The nucleotide and amino acid sequences and the gene order of the mitochondrial genome are highly informative for studying phylogeny, population genetics, and phylogeography. This study determined the complete mitochondrial genome of the caprellid species Caprella scaura. The mitochondrial genome of C. scaura has a total length of 15,079 bp, with an AT content of 66.43%. The mitochondrial genome contains typical gene components, including 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. In comparison with the mitochondrial genome of a gammarid, some distinct characteristics were found. For example, the order of the two conserved gene blocks is inverted between Gammaridea and C. scaura. In addition, two copies of almost identical control regions were found in the mitochondrial genome of C. scaura. These unique characteristics will be useful for determining the evolutionary history of the Caprellidea.     

Hearing impairment-associated KARS mutations lead to defects in aminoacylation of both cytoplasmic and mitochondrial tRNALys

Aminoacyl-tRNA synthetases(aaRSs) are ubiquitously expressed, essential enzymes, synthesizing aminoacyl-tRNAs for protein synthesis. Functional defects of aaRSs frequently cause various human disorders. Human KARS encodes both cytosolic and mitochondrial lysyl-tRNA synthetases(LysRSs). Previously, two mutations(c.1129 GA and c.517 TC) were identified that led to hearing impairment; however, the underlying biochemical mechanism is unclear. In the present study, we found that the two mutations have no impact on the incorporation of LysRS into the multiple-synthetase complex in the cytosol, but affect the cytosolic LysRS level, its tertiary structure, and cytosolic tRNA aminoacylation in vitro. As for mitochondrial translation, the two mutations have little effect on the steady-state level, mitochondrial targeting, and tRNA binding affinity of mitochondrial LysRS. However, they exhibit striking differences in charging mitochondrial tRNA~(Lys), with the c.517TC mutant being completely deficient in vitro and in vivo. We constructed two yeast genetic models, which are powerful tools to test the in vivo aminoacylation activity of KARS mutations at both the cytosolic and mitochondrial levels. Overall, our data provided biochemical insights into the potentially molecular pathological mechanism of KARS c.1129GA and c.517TC mutations and provided yeast genetic bases to investigate other KARS mutations in the future.     

New insights into mitogenomic phylogeny and copy number in eight indigenous sheep populations based on the ATP synthase and cytochrome c oxidase genes

The origins and phylogeny of different sheep breeds has been widely studied using polymorphisms within the mitochondrial hypervariable region. However, little is known about the mitochondrial DNA (mtDNA) content and phylogeny based on mtDNA protein-coding genes. In this study, we assessed the phylogeny and copy number of the mtDNA in eight indigenous (population size, n=184) and three introduced (n=66) sheep breeds in China based on five mitochondrial coding genes (COX1, COX2, ATP8, ATP6 and COX3). The mean haplotype and nucleotide diversities were 0.944 and 0.00322, respectively. We identified a correlation between the lineages distribution and the genetic distance, whereby Valley-type Tibetan sheep had a closer genetic relationship with introduced breeds (Dorper, Poll Dorset and Suffolk) than with other indigenous breeds. Similarly, the Median-joining profile of haplotypes revealed the distribution of clusters according to genetic differences. Moreover, copy number analysis based on the five mitochondrial coding genes was affected by the genetic distance combining with genetic phylogeny; we also identified obvious non-synonymous mutations in ATP6 between the different levels of copy number expressions. These results imply that differences in mitogenomic compositions resulting from geographical separation lead to differences in mitochondrial function.     

The ubiquitous mountain hare mitochondria: multiple introgressive hybridization in hares, genus Lepus

Climatic oscillations during the glaciations forced dramatic changes in species distributions, such that some presently temperate regions were alternately occupied by temperate and arctic species. These species could have met and hybridized during climatic transitions. This phenomenon happened for three hare species present in Iberia (Lepus granatensis, Lepus europaeus and Lepus castroviejoi), which display high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species presently extinct in Iberia. Here, we extend our previous geographical survey to determine whether the distribution of this mtDNA lineage extends beyond the northern half of the Iberian Peninsula, where it is found at high frequencies. We also review the taxonomy, distribution and molecular phylogeny of the genus Lepus. The phylogenetic inference reveals the presence of L. timidus-like mtDNA in several other hare species in Asia and North America, suggesting that the mitochondrial introgression observed in Iberia might be generalized. Comparison with the available nuclear gene phylogenies suggests that introgression could have happened repeatedly, possibly during different climatic transitions. We discuss demographic and adaptive scenarios that could account for the repetition in time and space of this spectacular phenomenon and suggest ways to improve our understanding of its determinants and consequences. Such high levels of introgressive hybridization should discourage attempts to revise hare taxonomy based solely on mtDNA.     

Automatic sequencing of mitochondrial tRNA genes in patients with mitochondrial encephalomyopathy

We have investigated nine children with infantile onset of mitochondrial myopathy and two adults with myoclonus epilepsy and ragged-red fibers (MERRF) and chronic progressive external ophthalmoplegia (CPEO), respectively. These patients lacked any of the previously known pathogenic tRNA mutations. Southern blot analysis of muscle mtDNA revealed no deletions. The tRNA genes of muscle mtDNA were sequenced. Restriction enxyme analysis of PCR fragments was performed to verify the presence of the mutations identified by automatic sequencing. Several tRNA mutations were found, but they were all homoplasmic. Furthermore, the mutations were either present in controls or did not change nucleotides conserved between species. This strongly suggests that none of the tRNA mutations identified in the 11 patients with mitochondrial encephalomyopathy was pathogenic. It can thus be concluded that mitochondrial tRNA mutations and mtDNA deletions probably are an infrequent cause of mitochondrial disorders in infants. Patients with MERRF and CPEO may lack both pathogenic point mutations of tRNA genes and deletions of mtDNA.     

Hares on ice: phylogeography and historical demographics of Lepus arcticus, L. othus, and L. timidus (Mammalia: Lagomorpha)

Phylogeographical investigations of arctic organisms provide spatial and temporal frameworks for interpreting the role of climate change on biotic diversity in high-latitude ecosystems. Phylogenetic analyses were conducted on 473 base pairs of the mitochondrial control region in 192 arctic hares (Lepus arcticus, Lepus othus, Lepus timidus) and two individual Lepus townsendii. The three arctic hare species are closely related. All L. othus individuals form one well-supported clade, L. arcticus individuals form two well-supported clades, and L. timidus individuals are scattered throughout the phylogeny. Arctic hare distribution was altered dramatically following post-Pleistocene recession of continental ice sheets. We tested for genetic signatures of population expansion for hare populations now found in deglaciated areas. Historical demographic estimates for 12 arctic hare populations from throughout their range indicate that L. arcticus and L. othus persisted in two separate North American arctic refugia (Beringia and High Canadian Arctic) during glacial advances of the Pleistocene, while the high genetic diversity in L. timidus likely reflects multiple Eurasian refugia.     

Analysis of the mitochondrial DNA from patients with Wolfram (DIDMOAD) syndrome

Wolfram or DIDMOAD (Diabetes Insipidus, Diabetes Mellitus, Optic Atrophy and Deafness) syndrome, which has long been known as an autosomal-recessive disorder, has recently been proposed to be a mitochondrial-mediated disease with either a nuclear or a mitochondrial genetic background. The phenotypic characteristics of the syndrome resemble those found in other mitochondrial (mt)DNA mediated disorders such as Leber's hereditary optic neuropathy (LHON) or maternally inherited diabetes and deafness (MIDD). Therefore, we looked for respective mtDNA alterations in blood samples from 7 patients with DIDMOAD syndrome using SSCP-analysis of PCR-amplified fragments, encompassing all mitochondrial ND and tRNA genes, followed by direct sequencing. Subsequently, we compared mtDNA variants identified in this disease group with those detected in a group of LHON patients (n = 17) and in a group of 69 healthy controls. We found that 4/7 (57%) DIDMOAD patients harbored a specific set of point mutations in tRNA and ND genes including the so-called class II or secondary LHON mutations at nucleotide positions (nps) 4216 and 4917 (haplogroup B). In contrast, LHON-patients were frequently (10/17, 59%) found in association with another cluster of mtDNA variants including the secondary LHON mutations at nps 4216 and 13708 and further mtDNA polymorphisms in ND genes (haplogroup A), overlapping with haplogroup B only by variants at nps 4216 and 11251. The frequencies of both haplogroups were significantly lower in the control group versus the respective disease groups. We propose that haplogroup B represents a susceptibility factor for DIDMOAD which, by interaction with further exogeneous or genetic factors, might increase the risk for disease. (Mol Cell Biochem 174: 209–213, 1997)     

两种毛蚜线粒体基因组比较分析

【目的】线粒体基因组分析已被应用于昆虫系统发育研究。本研究以蚜科Aphididae重要类群毛蚜亚科物种为代表,测定并比较分析了该类蚜虫的线粒体基因组特征,探讨了基于线粒体基因组信息的蚜虫系统发育关系重建。【方法】以毛蚜亚科三角枫多态毛蚜Periphyllus acerihabitans Zhang和针茅小毛蚜Chaetosiphella stipae Hille Ris Lambers,1947为研究对象,利用长短PCR相结合的方法测定线粒体基因组的序列,分析了基因组的基本特征;基于在线t RNAscan-SE Search Server搜索方法预测了t RNA的二级结构;基于12个物种(本研究获得的2个物种和10个Gen Bank上下载的物种数据)的蛋白编码基因(PCGs)序列,利用最大似然法和贝叶斯法重建了蚜科的系统发育关系。【结果】两种毛蚜均获得了约94%的线粒体基因组数据,P.acerihabitans获得了14 908 bp,控制区为1 205 bp;C.stipae获得了13 893 bp,控制区为609 bp。两种毛蚜同时获得33个基因,包含接近完整的13个蛋白编码基因(PCGs)(nad5不完整),18个tRNA,2个rRNA基因;ka/ks值表明,C.stipae的进化速率更快。从基因组组成、基因排列顺序、核苷酸组成分析、密码子使用情况、t RNA二级结构等特征来分析,两种蚜虫线粒体基因组基本特征相似。系统发育重建结果表明毛蚜亚科、蚜亚科的单系性得到了支持,毛蚜亚科位于蚜科的基部位置。【结论】两种毛蚜线粒体基因组的基本特征相似,符合蚜虫线粒体基因组的一般特征,两种线粒体基因组的长度差异主要来自控制区长度的不同;系统发育重建支持毛蚜亚科与蚜亚科的单系性,毛蚜亚科位于蚜科较为基部的位置。研究结果为蚜虫类系统发育重建提供了参考。     

Comparative mitochondrial genomic analyses of three chemosynthetic vesicomyid clams from deep‐sea habitats

Vesicomyid clams of the subfamily Pliocardinae are among the dominant chemosymbiotic bivalves found in sulfide‐rich deep‐sea habitats. Plastic morphologies and present molecular data could not resolve taxonomic uncertainties. The complete mitochondrial (mt) genomes will provide more data for comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group, and help to clarify generic classifications. In this study, we analyze the features and evolutionary dynamics of mt genomes from three Archivesica species (Archivesica sp., Ar. gigas and Ar. pacifica) pertaining to subfamily Pliocardinae. Sequence coverage is nearly complete for the three newly sequenced mt genomes, with only the control region and some tRNA genes missing. Gene content, base composition, and codon usage are highly conserved in these pliocardiin species. Comparative analysis revealed the vesicomyid have a relatively lower ratio of Ka/Ks, and all 13 protein‐coding genes (PGCs) are under strong purifying selection with a ratio of Ka/Ks far lower than one. Minimal changes in gene arrangement among vesicomyid species are due to the translocation trnaG in Isorropodon fossajaponicum. Additional tRNA genes were detected between trnaG and nad2 in Abyssogena mariana (trnaL3), Ab. phaseoliformis (trnaS3), and Phreagena okutanii (trnaM2), and display high similarity to other pliocardiin sequences at the same location. Single base insertion in multiple sites of this location could result in new tRNA genes, suggesting a possible tRNA arising from nongeneic sequence. Phylogenetic analysis based on 12 PCGs (excluding atp8) supports the monophyly of Pliocardiinae. These nearly complete mitogenomes provide relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group of chemosymbiotic bivalves.