Complete sequences of the highly rearranged molluscan mitochondrial genomes of the Scaphopod Graptacme eborea and the bivalve Mytilus edulis

Molluscan systematics: historical perspectives and the way ahead

Mollusca, the second-most diverse animal phylum, is estimated to have over 100,000 living species with great genetic and phenotypic diversity, a rich fossil record, and a considerable evolutionary significance. Early work on molluscan systematics was grounded in morphological and anatomical studies. With the transition from oligo gene Sanger sequencing to cutting-edge genomic sequencing technologies, molecular data has been increasingly utilised, providing abundant information for reconstructing the molluscan phylogenetic tree. However, relationships among and within most major lineages of Mollusca have long been contentious, often due to limited genetic markers, insufficient taxon sampling and phylogenetic conflict. Fortunately, remarkable progress in molluscan systematics has been made in recent years, which has shed light on how major molluscan groups have evolved. In this review of molluscan systematics, we first synthesise the current understanding of the molluscan Tree of Life at higher taxonomic levels. We then discuss how micromolluscs, which have adult individuals with a body size smaller than 5 mm, offer unique insights into Mollusca's vast diversity and deep phylogeny. Despite recent advancements, our knowledge of molluscan systematics and phylogeny still needs refinement. Further advancements in molluscan systematics will arise from integrating comprehensive data sets, including genome-scale data, exceptional fossils, and digital morphological data (including internal structures). Enhanced access to these data sets, combined with increased collaboration among morphologists, palaeontologists, evolutionary developmental biologists, and molecular phylogeneticists, will significantly advance this field.

Mitochondrial epigenetics brings new perspectives on doubly uniparental inheritance in bivalves

Mitochondrial epigenetics, particularly mtDNA methylation, is a flourishing field of research. MtDNA methylation appears to play multiple roles, including regulating mitochondrial transcription, cell metabolism and mitochondrial inheritance. In animals, bivalves with doubly uniparental inheritance (DUI) of mitochondria are the exception to the rule of maternal mitochondrial inheritance since DUI also involve a paternal mtDNA transmitted from the father to sons. The mechanisms underlying DUI are still unknown, but mtDNA methylation could play a role in its regulation. Here, we investigated mtDNA methylation levels and machinery in gonads of the mussel Mytilus edulis using methods based on antibodies, enzymatic cleavage and methylome sequencing. Our results confirm the presence in mitochondria of methylated cytosines and adenines and methyltransferases and unveil a more variable cytosine methylation state among males than females. Also, spermatid mtDNA is always methylated, while only few spermatozoa present methylated mtDNA suggesting a relation between cytosine methylation and development stage of male gametes. We propose that mtDNA methylation could play a role in the different fates of the parental mtDNAs in male and female embryos in M. edulis. Our study provides novel insights into the epigenetic landscape of bivalve mtDNA and highlights the multiple roles of mtDNA methylation in animals.

Complementing aculiferan mitogenomics: comparative characterization of mitochondrial genomes of Solenogastres (Mollusca, Aplacophora)

BACKGROUND

With the advances in high-throughput sequencing and bioinformatic pipelines, mitochondrial genomes have become increasingly popular for phylogenetic analyses across different clades of invertebrates. Despite the vast rise in available mitogenomic datasets of molluscs, one class of aplacophoran molluscs - Solenogastres (or Neomeniomorpha) - is still neglected.

RESULTS

Here, we present six new mitochondrial genomes from five families of Solenogastres (Amphimeniidae, Gymnomeniidae, Proneomeniidae, Pruvotinidae, Simrothiellidae), including the first complete mitogenomes, thereby now representing three of the four traditional orders. Solenogaster mitogenomes are variable in size (ranging from approximately 15,000 bp to over 17,000 bp). The gene order of the 13 protein coding genes and two rRNA genes is conserved in three blocks, but considerable variation occurs in the order of the 22 tRNA genes. Based on phylogenetic analyses and reconstruction of ancestral mitochondrial genomes of Aculifera, the position of (1) trnD gene between atp8 and atp6, (2) trnT and P genes between atp6 and nad5, and (3) trnL1 gene between G and E, resulting in a 'MCYWQGL1E'-block of tRNA genes, are all three considered synapomorphies for Solenogastres. The tRNA gene block 'KARNI' present in Polyplacophora and several conchiferan taxa is dissolved in Solenogastres.

CONCLUSION

Our study shows that mitogenomes are suitable to resolve the phylogenetic relationships among Aculifera and within Solenogastres, thus presenting a cost and time efficient compromise to approach evolutionary history in these clades.

The Mitochondrial Genome of Ylistrum japonicum (Bivalvia, Pectinidae) and Its Phylogenetic Analysis

The Ylistrum japonicum is a commercially valuable scallop known for its long-distance swimming abilities. Despite its economic importance, genetic and genomic research on this species is limited. This study presents the first complete mitochondrial genome of Y. japonicum. The mitochondrial genome is 19,475 bp long and encompasses 13 protein-coding genes, three ribosomal RNA genes, and 23 transfer RNA genes. Two distinct phylogenetic analyses were used to explore the phylogenetic position of the Y. japonicum within the family Pectinidae. Based on one mitochondrial phylogenetic analysis by selecting 15 Pectinidae species and additional outgroup taxa and one single gene phylogenetic analysis by 16S rRNA, two phylogenetic trees were constructed to provide clearer insights into the evolutionary placement of Y. japonicum within the family Pectinidae. Our analysis reveals that Ylistrum is a basal lineage to the Pectininae clade, distinct from its previously assigned tribe, Amusiini. This study offers critical insights into the genetic makeup and evolutionary history of Y. japonicum, enhancing our knowledge of this economically vital species.

The effect of mitochondrial recombination on fertilization success in blue mussels

The presence of doubly uniparental inheritance (DUI) in bivalves represents a unique mode of mitochondrial transmission, whereby paternal (male-transmitted M-type) and maternal (female-transmitted F-type) haplotypes are transmitted to offspring separately. Male embryos retain both haplotypes, but the M-type is selectively removed from females. Due to the presence of heteroplasmy in males, mtDNA can recombine resulting in a 'masculinized' haplotype referred to as Mf-type. While mtDNA recombination is usually rare, it has been recorded in multiple mussel species across the Northern Hemisphere. Given that mitochondria are the powerhouse of the cell, different mtDNA haplotypes may have different selective advantages under diverse environmental conditions. This may be particularly important for sperm fitness and fertilization success. In this study we aimed to i) determine the presence, prevalence of the Mf-type in Australian blue mussels (Mytilus sp.) and ii) investigate the effect of Mf-mtDNA on sperm performance (a fitness correlate). We found a high prevalence of recombined mtDNA (≈35 %) located within the control region of the mitochondrial genome, which occurred only in specimens that contained Southern Hemisphere mtDNA. The presence of two female mitotypes were identified in the studied mussels, one likely originating from the Northern Hemisphere, and the other either representing the endemic M. planulatus species or introduced genotypes from the Southern Hemisphere. Despite having recombination events present in a third of the studied population, analysis of sperm performance indicated no difference in fertilization success related to mitotype.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.)

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

Mitochondrial Genome Evolution in Annelida-A Systematic Study on Conservative and Variable Gene Orders and the Factors Influencing its Evolution

The mitochondrial genomes of Bilateria are relatively conserved in their protein-coding, rRNA, and tRNA gene complement, but the order of these genes can range from very conserved to very variable depending on the taxon. The supposedly conserved gene order of Annelida has been used to support the placement of some taxa within Annelida. Recently, authors have cast doubts on the conserved nature of the annelid gene order. Various factors may influence gene order variability including, among others, increased substitution rates, base composition differences, structure of noncoding regions, parasitism, living in extreme habitats, short generation times, and biomineralization. However, these analyses were neither done systematically nor based on well-established reference trees. Several focused on only a few of these factors and biological factors were usually explored ad-hoc without rigorous testing or correlation analyses. Herein, we investigated the variability and evolution of the annelid gene order and the factors that potentially influenced its evolution, using a comprehensive and systematic approach. The analyses were based on 170 genomes, including 33 previously unrepresented species. Our analyses included 706 different molecular properties, 20 life-history and ecological traits, and a reference tree corresponding to recent improvements concerning the annelid tree. The results showed that the gene order with and without tRNAs is generally conserved. However, individual taxa exhibit higher degrees of variability. None of the analyzed life-history and ecological traits explained the observed variability across mitochondrial gene orders. In contrast, the combination and interaction of the best-predicting factors for substitution rate and base composition explained up to 30% of the observed variability. Accordingly, correlation analyses of different molecular properties of the mitochondrial genomes showed an intricate network of direct and indirect correlations between the different molecular factors. Hence, gene order evolution seems to be driven by molecular evolutionary aspects rather than by life history or ecology. On the other hand, variability of the gene order does not predict if a taxon is difficult to place in molecular phylogenetic reconstructions using sequence data or not. We also discuss the molecular properties of annelid mitochondrial genomes considering canonical views on gene evolution and potential reasons why the canonical views do not always fit to the observed patterns without making some adjustments. [Annelida; compositional biases; ecology; gene order; life history; macroevolution; mitochondrial genomes; substitution rates.].

Evidence of multiple genome duplication events in Mytilus evolution

Background

Molluscs remain one significantly under-represented taxa amongst available genomic resources, despite being the second-largest animal phylum and the recent advances in genomes sequencing technologies and genome assembly techniques. With the present work, we want to contribute to the growing efforts by filling this gap, presenting a new high-quality reference genome for Mytilus edulis and investigating the evolutionary history within the Mytilidae family, in relation to other species in the class Bivalvia.

Results

Here we present, for the first time, the discovery of multiple whole genome duplication events in the Mytilidae family and, more generally, in the class Bivalvia. In addition, the calculation of evolution rates for three species of the Mytilinae subfamily sheds new light onto the taxa evolution and highlights key orthologs of interest for the study of Mytilus species divergences.

Conclusions

The reference genome presented here will enable the correct identification of molecular markers for evolutionary, population genetics, and conservation studies. Mytilidae have the capability to become a model shellfish for climate change adaptation using genome-enabled systems biology and multi-disciplinary studies of interactions between abiotic stressors, pathogen attacks, and aquaculture practises.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08575-9.

Comparative Analysis of Four Complete Mitochondrial Genomes of Epinephelidae (Perciformes)

Groupers are commercial, mainly reef-associated fishes, classified in the family Epinephelidae (Perciformes). This study first sequenced the complete mitogenomes of Cephalopholis leopardus, Cephalopholis spiloparaea, Epinephelus amblycephalus, and Epinephelus hexagonatus. The lengths of the four Epinephelidae mitogenomes ranged from 16,585 base pair (bp) to 16,872 bp with the typical gene order. All tRNA genes had a typical cloverleaf structure, except the tRNA-Ser (AGY) gene which was lacking the entire dihydrouridine arm. The ratio of nonsynonymous substitution (Ka) and synonymous substitution (Ks) indicated that four groupers were suffering a purifying selection. Phylogenetic relationships were reconstructed by Bayesian inference (BI) and maximum likelihood (ML) methods based on all mitogenomic data of 41 groupers and 2 outgroups. The identical topologies result with high support values showed that Cephalopholis and Epinephelus are not monophyletic genera. Anyperodon and Cromileptes clustered to Epinephelus. Aethaloperca rogaa and Cephalopholis argus assembled a clad. Cephalopholis leopardus, C. spiloparaea, and Cephalopholis miniata are also in a clade. Epinephelushexagonatus is close to Epinephelus tauvina and Epinephelus merra, and E. amblycephalus is a sister group with Epinephelus stictus. More mitogenomic data from Epinephelidae species are essential to understand its taxonomic status with the family Serranidae.

Complex mitogenomic rearrangements within the Pectinidae (Mollusca: Bivalvia)

BACKGROUND

Scallops (Bivalvia: Pectinidae) present extraordinary variance in both mitochondrial genome size, structure and content, even when compared to the extreme diversity documented within Mollusca and Bivalvia. In pectinids, mitogenome rearrangements involve protein coding and rRNA genes along with tRNAs, and different genome organization patterns can be observed even at the level of Tribes. Existing pectinid phylogenies fail to resolve some relationships in the family, Chlamydinae being an especially problematic group.

RESULTS

In our study, we sequenced, annotated and characterized the mitochondrial genome of a member of Chlamydinae, Mimachlamys varia-a species of commercial interest and an effective bioindicator-revealing yet another novel gene arrangement in the Pectinidae. The phylogeny based on all mitochondrial protein coding and rRNA genes suggests the paraphyly of the Mimachlamys genus, further commending the taxonomic revision of the classification within the Chlamydinae subfamily. At the scale of the Pectinidae, we found that 15 sequence blocks are involved in mitogenome rearrangements, which behave as separate units.

CONCLUSIONS

Our study reveals incongruities between phylogenies based on mitochondrial protein-coding versus rRNA genes within the Pectinidae, suggesting that locus sampling affects phylogenetic inference at the scale of the family. We also conclude that the available taxon sampling does not allow for understanding of the mechanisms responsible for the high variability of mitogenome architecture observed in the Pectinidae, and that unraveling these processes will require denser taxon sampling.

Morphological and complete mitogenomic characterisation of the acanthocephalan Polymorphus minutus infecting the duck Anas platyrhynchos

Morphological characteristics of the acanthocephalan Polymorphus minutus (Goeze, 1782), which was collected from the duck Anas platyrhynchos Linnaeus in the Czech Republic, are described. The mitochondrial (mt) genome of P. minutus was sequenced, with a total length of 14,149 bp, comprising 36 genes including 12 protein coding genes (PCGs), 22 transfer RNA (tRNA) genes and two ribosomal RNA genes (rrnL and rrnS). This genome is similar to the mt genomes of other syndermatan species. All these genes were encoded on the same DNA strand and in the same orientation. The overall nucleotide composition of the P. minutus mt genome was 38.2% T, 27.3% G, 26.2% A, and 8.3% C. The amino acid sequences of 12 PCGs for mt genomes of 28 platyzoans, including P. minutus, were used for phylogenetic analysis, and the resulting topology recovers P. minutus as sister to Southwellina hispida (Van Cleave, 1925), and the two taxa form a sister clade to Centrorhynchus aluconis (Müller, 1780) and Plagiorhynchus transversus (Rudolphi, 1819), which are all species in the Palaeacanthocephala, thus supporting the monophyly of this class.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Phylogenetic Relationships and Adaptation in Deep-Sea Mussels: Insights from Mitochondrial Genomes

Mitochondrial genomes (mitogenomes) are an excellent source of information for phylogenetic and evolutionary studies, but their application in marine invertebrates is limited. In the present study, we utilized mitogenomes to elucidate the phylogeny and environmental adaptation in deep-sea mussels (Mytilidae: Bathymodiolinae). We sequenced and assembled seven bathymodioline mitogenomes. A phylogenetic analysis integrating the seven newly assembled and six previously reported bathymodioline mitogenomes revealed that these bathymodiolines are divided into three well-supported clades represented by five Gigantidas species, six Bathymodiolus species, and two "Bathymodiolus" species, respectively. A Common interval Rearrangement Explorer (CREx) analysis revealed a gene order rearrangement in bathymodiolines that is distinct from that in other shallow-water mytilids. The CREx analysis also suggested that reversal, transposition, and tandem duplications with subsequent random gene loss (TDRL) may have been responsible for the evolution of mitochondrial gene orders in bathymodiolines. Moreover, a comparison of the mitogenomes of shallow-water and deep-sea mussels revealed that the latter lineage has experienced relaxed purifying selection, but 16 residues of the atp6, nad4, nad2, cob, nad5, and cox2 genes have underwent positive selection. Overall, this study provides new insights into the phylogenetic relationships and mitogenomic adaptations of deep-sea mussels.

Complete mitochondrial genomes of Epeorus carinatus and E. dayongensis (Ephemeroptera: Heptageniidae): Genomic comparison and phylogenetic inference

The current research on Ephemeroptera is mainly based on its morphology, since only small numbers of mitogenomes have been reported. In this study, the mitogenomes of Epeorus carinatus (15,338 bp) and E. dayongensis (15,609 bp) were sequenced, annotated and compared to genome data from congeners. Both mitogenomes had 23 tRNA genes including standard 22 and one extra tRNA^Met. The duplicated tRNA^Met gene had been found in other heptageniid species except Paegniodes cupulatus, suggesting it could be used as a molecular synapomorphy for partial Heptageniidae. The phylogenetic analyses based on Bayesian Inference (BI) and Maximum Likelihood (ML) showed that Heptageniidae was monophyletic and the relationships among known Epeorus species were ((E. carinatus + E. herklotsi) + (E. dayongensis + E. sp. 1)), which implied the focal species E. carinatus and E. dayongensis should be grouped into different subgenera.

Complete mitochondrial genomes of two scaphopod molluscs

Complete mitochondrial genomes were determined for two scaphopod molluscs: the dentaliid Antalis entalis and an unidentified Antarctic gadilid. Both genomes are complete except, in Gadilida sp. indet., a short stretch of nad5 was undetermined and trnR could not be annotated. Organization of the Gadilida sp. genome is nearly identical to that previously reported for the gadilid Siphonodentalium whereas trnK, nad5, trnD, nad4, and nad4l are transposed to the opposite strand in the previously published Graptacme genome relative to that of Antalis. Phylogenetic analysis of the 13 protein-coding and 2 rRNA genes recovered Scaphopoda, Gadilida, and Dentaliida monophyletic with maximal support.

A mitochondrial genome phylogeny of Mytilidae (Bivalvia: Mytilida)

The family Mytilidae is a family of bivalve mussels that are distributed worldwide in diverse marine habitats. Within the family, classification systems and phylogenetic relationships among subfamilies remain not yet fully resolved. In this study, we newly determined 9 mitochondrial genome sequences from 7 subfamilies: Bathymodiolus thermophilus (Bathymodiolinae), Modiolus nipponicus (Modiolinae), Lithophaga curta (the first representative of Lithophaginae), Brachidontes mutabilis (Brachidontinae), Mytilisepta virgata (Brachidontinae), Mytilisepta keenae (Brachidontinae), Crenomytilus grayanus (Mytilinae), Gregariella coralliophaga (Crenellinae), and Septifer bilocularis (the first representative of Septiferinae). Phylogenetic trees using maximum likelihood and Bayesian inference methods for 28 mitochondrial genomes (including 19 previously published sequences) showed two major clades with high support values: Clade 1 ((Bathymodiolinae + Modiolinae) + (Lithophaginae + Limnoperninae)) and Clade 2 (((Mytilinae + Crenellinae) + Septiferinae) + Brachidontinae). The position of the genus Lithophaga (representing Lithophaginae) differed from a previously published molecular phylogeny. Divergence time analysis with a molecular clock indicated that lineage splitting among the major subfamilies of Mytilidae (including the habitat transition from marine to freshwater environments by ancestral Limnoperninae) occurred in the Mesozoic period, coinciding with high diversification rates of marine fauna during that time. This is the first mitochondrial genome-based phylogenetic study of the Mytilidae that covers nearly all subfamily members, excluding the subfamily Dacrydiinae.

A novel gene arrangement among the Stylommatophora by the complete mitochondrial genome of the terrestrial slug Meghimatium bilineatum (Gastropoda, Arionoidea)

Stylommatophora is a main clade of Gastropoda that encompasses approximately 112 gastropod families and may exceed a total of 30,000 species. Twenty-four complete stylommatophoran mitogenomes have been sequenced to date, yet our understanding of mitochondrial evolution in stylommatophorans is still in its infancy. To further expand the set of available mitogenomes, we sequenced the mitogenome of Meghimatium bilineatum (Arionoidea: Philomycidae), a widespread land slug in East Asia. This is the first report on a mitogenome of the superfamily Arionoidea, and indeed on a terrestrial slug. The mitogenome of Meghimatium bilineatum comprises 13,972 bp and exhibits a novel, highly distinctive gene arrangement among the Stylommatophora. Phylogenetic reconstructions based on the sequences of all protein-coding genes consistently recovered Meghimatium bilineatum as sister-group of the Succineidae. A phylogenetic reconstruction based on gene order, however, suggested a highly divergent tree topology, which is less credible when taking into account prior knowledge of stylommatophoran relationships. Our CREx (Common interval Rearrangement Explorer) analysis suggested that three successive events of tandem duplication random loss (TDRL) best explain the evolutionary process of gene order rearrangement in Meghimatium bilineatum from an ancestral stylommatophoran mitogenome. The present example offers new insights into the mechanisms of mitogenome rearrangements in gastropods at large and into the usefulness of mitogenomic gene order as a phylogenetic marker.

Mitochondrial genome of Argopecten irradians reveals higher-level phylogenetic relationships in Anisomyaria

The complete mitochondrial genome of Argopecten irradians strain Zhongkehong was sequenced and annotated: it is 16,212 bp in length and contains twelve protein-coding genes (atp8 is absent, as in most species in Anisomyaria), two ribosomal RNA genes, and 21 transfer RNA genes (trnS is absent and there are two copies of trnF). The heavy strand has an overall A + T content of 57.3%; GC and AT skews are 0.249 and -0.262, respectively, indicating more Gs and more Ts than Cs and As. Phylogenetic analysis based on Bayesian Inference and Maximum Likelihood of the twelve protein-coding genes shows that A. irradians has close relationships with A. purpuratus and A. ventricosus; this indicated that A. irradians belongs to the Pectinidae family. The Pectinidae was sister to (Ostreidae + Mytilidae). This work provides general information on the evolution of cultured scallops.

Complete paternally inherited mitogenomes of two freshwater mussels Unio pictorum and Sinanodonta woodiana (Bivalvia: Unionidae)

Freshwater bivalves from the family Unionidae usually have two very divergent mitogenomes, inherited according to the doubly uniparental model. The early divergence of these two mitogenomic lineages gives a unique opportunity to use two mitogenomic data sets in a single phylogenetic context. However, the number of complete sequences of the maternally inherited mitogenomes of these animals available in GenBank greatly exceeds that of the paternally inherited mitogenomes. This is a problem for phylogenetic reconstruction because it limits the use of both mitogenomic data sets. Moreover, since long branch attraction phenomenon can bias reconstructions if only a few but highly divergent taxa are considered, the shortage of the faster evolving paternally inherited mitogenome sequences is a real problem. Here we provide, for the first time, complete sequences of the M mitogenomes sampled from Polish populations of two species: native Unio pictorum and invasive Sinanodonta woodiana. It increases the available set of mitogenomic pairs to 18 species per family, and allows unambiguous reconstruction of phylogenetic relationships among them. The reconstructions based on M and F mitogenomes which were separated for many millions of years, and subject to differing evolutionary dynamics, are fully congruent.

Complete mitochondrial genome of the first deep-sea spongicolid shrimp Spongiocaris panglao (Decapoda: Stenopodidea): Novel gene arrangement and the phylogenetic position and origin of Stenopodidea

Stenopodidea Claus, 1872 (Crustacea: Decapoda) is one of the major groups of decapods crustaceans. Hitherto, only one complete mitochondrial genome (mitogenome) from the family Stenopodidae is available for the infraorder Stenopodidea. Here, we determined the complete mitogenome of Spongiocaris panglao de Grave and Saito, 2016 using Illumina sequencing, representing the first species from the family Spongicolidae. The 15,909 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. Although the overall genome organization is typical for metazoans, the mitogenome of S. panglao shows some derived characters. A + T content of 77.42% in S. pamglao mitogenome is second-highest among the dacapods described to date. The trnR gene exhibit modified secondary structure with the TψC loop completely missing, which might be a putative autapomorphy of S. pamglao mitogenome. Compared with the shallow-water stenopodidean species S. hispidus, the control region of S. pamglao exhibits three characteristics: larger size, higher A + T content, and more tandem repeat sequences. The gene order exhibited difference from the ancestral mitogenome pattern of the Pancrustacea, with 5 tRNA genes rearrangement. The result from BI was agreed with most morphological characters and molecular evidences, revealing that Stenopodidea and Reptantia had the closest relationship, as the sister group of Caridea. Still, the alternative hypothesis supported from ML topology cannot be completely rejected based on the current data. Estimated times revealed that the two stenopodideans families Stenopodidae and Spongicolidae diverged from each other around 122 Mya. The divergence time of spongicolid shrimp is in good agreement with the origin of their hexactinellid hosts (78-144 Mya).

Actively transcribed and expressed atp8 gene in Mytilus edulis mussels

Background

Animal mitochondrial genomes typically encode 37 genes: 13 proteins, 22 tRNAs and two rRNAs. However, many species represent exceptions to that rule. Bivalvia along with Nematoda and Platyhelminthes are often suspected to fully or partially lack the ATP synthase subunit 8 (atp8) gene. This raises the question as to whether they are really lacking this gene or is this maybe an annotation problem? Among bivalves, Mytilus edulis has been inferred to lack an ATP8 gene since the characterization of its mitochondrial genome in 1992. Even though recent bioinformatic analyses suggested that atp8 is present in Mytilus spp., due to high divergence in predicted amino acid sequences, the existence of a functional atp8 gene in this group remains controversial.

Results

Here we demonstrate that M. edulis mitochondrial open reading frames suggested to be atp8 (in male and female mtDNAs) are actively translated proteins. We also provide evidence that both proteins are an integral part of the ATP synthase complex based on in-gel detection of ATP synthase activity and two-dimensional Blue-Native and SDS polyacrylamide electrophoresis.

Conclusion

Many organisms (e.g., Bivalvia along with Nematoda and Platyhelminthes) are considered to be lacking certain mitochondrial genes often only based on poor similarity between protein coding gene sequences in genetically closed species. In some situations, this may lead to the inference that the ATP8 gene is absent, when it is in fact present, but highly divergent. This shows how important complementary role protein-based approaches, such as those in the present study, can provide to bioinformatic, genomic studies (i.e., ability to confirm the presence of a gene).

Mitochondrial gene order evolution in Mollusca: Inference of the ancestral state from the mtDNA of Chaetopleura apiculata (Polyplacophora, Chaetopleuridae)

The mitochondrial genome architecture of polyplacophorans has been usually regarded as being very ancient in comparison to all mollusks. However, even if some complete chiton mtDNAs have been recently sequenced, thorough studies of their evolution are lacking. To further expand the set of complete chiton mtDNAs and perform such analysis, we sequenced the mitochondrial genome of the Eastern beaded chiton Chaetopleura apiculata (Chaetopleuridae) using next-generation sequencing. With mitochondrial sequences from all available chiton mtDNAs, we also built a phylogeny on which we reconstructed the evolution of gene arrangement in this class. The arrangement of C. apiculata proved to be the most primitive known so far for polyplacophorans. Comparing this gene order to those of other molluscan classes, we found that it most probably is the original gene order of the last common ancestor to all extant Mollusca. The ancient mitochondrial genome organization of C. apiculata is an important information that may help reconstructing the phylogeny of Mollusca and their relationship with other lophotrochozoans.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

Curious bivalves: Systematic utility and unusual properties of anomalodesmatan mitochondrial genomes

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low.

Complete female mitochondrial genome of Mytilus chilensis

The controversy surrounding the origin of antitropical distribution of Mytilus mussels and the taxonomic status of southern hemisphere populations remain unsolved, despite the efforts. One of the limiting factors remains the lack of the complete sequences of the representative mitochondrial genomes which would allow their proper comparison with the relatively well-represented northern hemisphere congeneric mussels. To fill this gap we sequenced the representative maternal (F) genome of a native Chilean mussel. The genome is 16,748bp long and structurally identical to the northern hemisphere M. edulis and M. galloprovincialis F genomes. However, the genetic distance from them (≈5%) is twice as high as the maximum distance between them (<2.5%). Thus, the notion that the name M. chilensis should be used for native Chilean Mytilus mussels, with the same rank as M. galloprovincialis and M. edulis is supported.

Monoplacophoran mitochondrial genomes: convergent gene arrangements and little phylogenetic signal

BACKGROUND

Although recent studies have greatly advanced understanding of deep molluscan phylogeny, placement of some taxa remains uncertain as different datasets support competing class-relationships. Traditionally, morphologists have placed Monoplacophora, a group of morphologically simple, limpet-like molluscs as sister group to all other conchiferans (shelled molluscs other than Polyplacophora), a grouping that is supported by the latest large-scale phylogenomic study that includes Laevipilina. However, molecular datasets dominated by nuclear ribosomal genes support Monoplacophora + Polyplacophora (Serialia). Here, we evaluate the potential of mitochondrial genome data for resolving placement of Monoplacophora.

RESULTS

Two complete (Laevipilina antarctica and Vema ewingi) and one partial (Laevipilina hyalina) mitochondrial genomes were sequenced, assembled, and compared. All three genomes show a highly similar architecture including an unusually high number of non-coding regions. Comparison of monoplacophoran gene order shows a gene arrangement pattern not previously reported; there is an inversion of one large gene cluster. Our reanalyses of recently published polyplacophoran mitogenomes show, however, that this feature is also present in some chiton species. Maximum Likelihood and Bayesian Inference analyses of 13 mitochondrial protein-coding genes failed to robustly place Monoplacophora and hypothesis testing could not reject any of the evaluated placements of Monoplacophora.

CONCLUSIONS

Under both serialian or aculiferan-conchiferan scenarios, the observed gene cluster inversion appears to be a convergent evolution of gene arrangements in molluscs. Our phylogenetic results are inconclusive and sensitive to taxon sampling. Aculifera (Polyplacophora + Aplacophora) and Conchifera were never recovered. However, some analyses recovered Serialia (Monoplacophora + Polyplacophora), Diasoma (Bivalvia + Scaphopoda) or Pleistomollusca (Bivalvia + Gastropoda). Although we could not shed light on deep evolutionary traits of Mollusca we found unique patterns of gene arrangements that are common to monoplacophoran and chitonine polyplacophoran species but not to acanthochitonine Polyplacophora. Complete mitochondrial genome of Laevipilina antarctica.

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.

Mitochondrial Genomes of Kinorhyncha: trnM Duplication and New Gene Orders within Animals

Many features of mitochondrial genomes of animals, such as patterns of gene arrangement, nucleotide content and substitution rate variation are extensively used in evolutionary and phylogenetic studies. Nearly 6,000 mitochondrial genomes of animals have already been sequenced, covering the majority of animal phyla. One of the groups that escaped mitogenome sequencing is phylum Kinorhyncha-an isolated taxon of microscopic worm-like ecdysozoans. The kinorhynchs are thought to be one of the early-branching lineages of Ecdysozoa, and their mitochondrial genomes may be important for resolving evolutionary relations between major animal taxa. Here we present the results of sequencing and analysis of mitochondrial genomes from two members of Kinorhyncha, Echinoderes svetlanae (Cyclorhagida) and Pycnophyes kielensis (Allomalorhagida). Their mitochondrial genomes are circular molecules approximately 15 Kbp in size. The kinorhynch mitochondrial gene sequences are highly divergent, which precludes accurate phylogenetic inference. The mitogenomes of both species encode a typical metazoan complement of 37 genes, which are all positioned on the major strand, but the gene order is distinct and unique among Ecdysozoa or animals as a whole. We predict four types of start codons for protein-coding genes in E. svetlanae and five in P. kielensis with a consensus DTD in single letter code. The mitochondrial genomes of E. svetlanae and P. kielensis encode duplicated methionine tRNA genes that display compensatory nucleotide substitutions. Two distant species of Kinorhyncha demonstrate similar patterns of gene arrangements in their mitogenomes. Both genomes have duplicated methionine tRNA genes; the duplication predates the divergence of two species. The kinorhynchs share a few features pertaining to gene order that align them with Priapulida. Gene order analysis reveals that gene arrangement specific of Priapulida may be ancestral for Scalidophora, Ecdysozoa, and even Protostomia.

Comparative Large-Scale Mitogenomics Evidences Clade-Specific Evolutionary Trends in Mitochondrial DNAs of Bivalvia

Despite the figure of complete bivalve mitochondrial genomes keeps growing, an assessment of the general features of these genomes in a phylogenetic framework is still lacking, despite the fact that bivalve mitochondrial genomes are unusual under different aspects. In this work, we constructed a dataset of one hundred mitochondrial genomes of bivalves to perform the first systematic comparative mitogenomic analysis, developing a phylogenetic background to scaffold the evolutionary history of the class' mitochondrial genomes. Highly conserved domains were identified in all protein coding genes; however, four genes (namely, atp6, nad2, nad4L, and nad6) were found to be very divergent for many respects, notwithstanding the overall purifying selection working on those genomes. Moreover, the atp8 gene was newly annotated in 20 mitochondrial genomes, where it was previously declared as lacking or only signaled. Supernumerary mitochondrial proteins were compared, but it was possible to find homologies only among strictly related species. The rearrangement rate on the molecule is too high to be used as a phylogenetic marker, but here we demonstrate for the first time in mollusks that there is correlation between rearrangement rates and evolutionary rates. We also developed a new index (HERMES) to estimate the amount of mitochondrial evolution. Many genomic features are phylogenetically congruent and this allowed us to highlight three main phases in bivalve history: the origin, the branching of palaeoheterodonts, and the second radiation leading to the present-day biodiversity.

The complete mitochondrial genome of the golden mussel Limnoperna fortunei and comparative mitogenomics of Mytilidae

Here we describe the mitochondrial genome of the golden mussel Limnoperna fortunei, an Asian bivalve which has become one of the most aggressive invasive species in Japan and South America. The mitochondrial genome of L. fortunei does not present conserved gene arrangement when compared to the other Mytilidae species suggesting a high degree of gene recombination in the mitochondria of this clade. In addition, the golden mussel mitogenome encodes two copies of tRNA(Lys) and presents a putative pseudogene for the atp8 gene sequence that encodes a 27 amino acid peptide containing an in-frame stop codon. The presence of this pseudogene raises the question as to whether atp8 is encoded in some bivalve mitochondrial genomes or not. The phylogenetic analysis of all complete mitochondrial genomes available from Mytilidae mussels confirmed the close evolutionary relationships among bivalves from the genus Mytilys and placed L. fortunei coming from a more ancestral branch on the family. The supermatrix phylogeny described used the concatenation of all 12 genes from the mitochondria and disputed the monophyly of the genus Perna, as Perna perna was shown to be more closely related to Brachidontes exustus than to Perna viridis. The comparative analysis of mitogenome synteny also confirmed the polyphyly of the genus Perna. The complete and annotated mitogenome has been published in GenBank under the accession number KP756905.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

The mitochondrial genome of the land snail Camaenacicatricosa (Müller, 1774) (Stylommatophora, Camaenidae): the first complete sequence in the family Camaenidae

The complete mitochondrial (mt) genome of the snail Camaenacicatricosa (Müller, 1774) has been sequenced and annotated in this study. The entire circular genome is 13,843 bp in size and represents the first camaenid mt genome, with content of 31.9%A, 37.9%T, 13.5%C and 16.7%G. Gene content, codon usage and base organization show similarity to a great extent to the sequenced mt genome from Stylommatophora, whereas, gene order is different from them, especially the positions of tRNA(Cys) , tRNA(Phe) , COII, tRNA(Asp) , tRNA(Gly) , tRNA(His) and tRNA(Trp) . All protein coding genes use standard initiation codons ATN except for COII with GTG as start signal. Conventional stop codons TAA and TAG have been assigned to all protein coding genes. All tRNA genes possess the typical clover leaf structure, but the TψC arm of tRNA(Asp) and dihydrouridine arm of tRNA(Ser(AGN)) only form a simple loop. Shorter intergenic spacers have been found in this mt genome. Phylogenetic study based on protein coding genes shows close relationship of Camaenidae and Bradybaenidae. The presented phylogeny is consistent with the monophyly of Stylommatophora.

Scorched mussels (BIVALVIA: MYTILIDAE: BRACHIDONTINAE) from the temperate coasts of South America: phylogenetic relationships, trans-Pacific connections and the footprints of Quaternary glaciations

This study addresses aspects of the phylogeny and phylogeography of scorched mussels (BIVALVIA: MYTILIDAE: BRACHIDONTINAE) from southern South America (Argentina and Chile), as well as their ecophylogenetic implications. Relationships were inferred from sequences of two nuclear (28S and 18S) and one mitochondrial (COI) genes, using Bayesian and maximum likelihood analyses. Our results indicate that the monophyletic BRACHIDONTINAE include three well supported clades: [i] Brachidontes Swainson (=Hormomya Mörch), [ii] Ischadium Jukes-Browne+Geukensia van de Poel, and [iii] Austromytilus Laseron+Mytilisepta Habe (usually considered a member of the SEPTIFERINAE)+Perumytilus Olsson. Species of clade [iii] are distributed along the temperate coasts of the Pacific Ocean. Available evidence supports divergence between Austromytilus (Australia) and Perumytilus (South American) following the breakup of Australian, Antarctic and South American shelves. Four brachidontins occur in southern South America: Brachidontes rodriguezii (d'Orbigny), B. granulatus (Hanley), and two genetically distinct clades of Perumytilus. The latter are confined to the Chile-Peru (North Clade) and Magellanic (South Clade) Biogeographic Provinces, respectively warm- and cold-temperate. The South Clade is the only brachidontin restricted to cold-temperate waters. Biogeographic considerations and the fossil record prompted the hypothesis that the South Clade originated from the North Clade by incipient peripatric differentiation, followed by isolation during the Quaternary glaciations, genetic differentiation in the non-glaciated coasts of eastern Patagonia, back-expansion over southern Chile following post-LGM de-glaciation, and development of a secondary contact zone between the two clades in south-central Chile. Evidence of upper Pleistocene expansion of the South Clade parallels similar results on other organisms that have colonized coastal ecosystems from eastern Patagonia since the LGM, apparently occupying free ecological space. We emphasize that the assembly of communities cannot be explained solely in terms of environmental drivers, as history also matters.

The complete mitochondrial genome of Scutopus ventrolineatus (Mollusca: Chaetodermomorpha) supports the Aculifera hypothesis

Background

With more than 100000 living species, mollusks are the second most diverse metazoan phylum. The current taxonomic classification of mollusks recognizes eight classes (Neomeniomorpha, Chaetodermomorpha, Polyplacophora, Monoplacophora, Cephalopoda, Gastropoda, Bivalvia, and Scaphopoda) that exhibit very distinct body plans. In the past, phylogenetic relationships among mollusk classes have been contentious due to the lack of indisputable morphological synapomorphies. Fortunately, recent phylogenetic analyses based on multi-gene data sets are rendering promising results. In this regard, mitochondrial genomes have been widely used to reconstruct deep phylogenies. For mollusks, complete mitochondrial genomes are mostly available for gastropods, bivalves, and cephalopods, whereas other less-diverse lineages have few or none reported.

Results

The complete DNA sequence (14662 bp) of the mitochondrial genome of the chaetodermomorph Scutopus ventrolineatus Salvini-Plawen, 1968 was determined. Compared with other mollusks, the relative position of protein-coding genes in the mitochondrial genome of S. ventrolineatus is very similar to those reported for Polyplacophora, Cephalopoda and early-diverging lineages of Bivalvia and Gastropoda (Vetigastropoda and Neritimorpha; but not Patellogastropoda). The reconstructed phylogenetic tree based on combined mitochondrial and nuclear sequence data recovered monophyletic Aplacophora, Aculifera, and Conchifera. Within the latter, Cephalopoda was the sister group of Gastropoda and Bivalvia + Scaphopoda.

Conclusions

Phylogenetic analyses of mitochondrial sequences showed strong among-lineage rate heterogeneity that produced long-branch attraction biases. Removal of long branches (namely those of bivalves and patellogastropods) ameliorated but not fully resolved the problem. Best results in terms of statistical support were achieved when mitochondrial and nuclear sequence data were concatenated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0197-9) contains supplementary material, which is available to authorized users.

Comparative mitogenomic analysis reveals cryptic species: A case study in Mactridae (Mollusca: Bivalvia)

The Chinese surf clam Mactra chinensis Philippi, 1846 is a commercially important marine bivalve belonging to the family Mactridae (Mollusca: Bivalvia). In this study, the M. chinensis mitochondrial genomic features are analyzed. The genome has 34 genes on the same strand, lacking atp8 and both trnS (trnS1 and trnS2) as compared with the typical gene content of metazoan mitochondrial genomes. The A+T content of M. chinensis mitochondrial genome is 63.72%, which is slightly lower than that of M. veneriformis (67.59%) and Coelomactra antiquata (64.33% and 64.14% for the samples from Ri Zhao, Shandong Province, and Zhang Zhou, Fujian Province, China, respectively) in the same family. There are 22 NCRs in the M. chinensis mitochondrial genome, accounting for 12.91% of the genome length. The longest NCR (1,075bp in length) is located between trnT and trnQ. A TRS (127bp×8.15) accounts for 96.3% (1,035/1,075) of this NCR. The occurrence of TRS in NCR is shared by the two Mactra mitochondrial genomes, but is not found in the two Coelomactra mitochondrial genomes. A phylogenetic tree constructed based on 12 PCGs of 25 bivalve mitochondrial genomes shows that all seven genera (Mactra, Coelomactra, Paphia, Meretrix, Solen, Mytilus, and Crassostrea) constitute monophyletic groups with very high support values. Pairwise genetic distance analyses indicate that the genetic distance of C. antiquata from the two localities is 0.084, which is greater than values between congeneric species, such as those in Mactra, Mytilus, Meretrix, and Crassostrea. The results show that the C. antiquata from the two localities represent cryptic species.

Complete mitochondrial genome of the giant African snail, Achatina fulica (Mollusca: Achatinidae): a novel location of putative control regions (CR) in the mitogenome within Pulmonate species

The whole sequence (15,057 bp) of the mitochondrial DNA (mtDNA) of the terrestrial snail Achatina fulica (order Stylommatophora) was determined. The mitogenome, as the typical metazoan mtDNA, contains 13 protein-coding genes (PCG), 2 ribosomal RNA genes (rRNA) and 22 transfer RNA genes (tRNA). The tRNA genes include two trnS without standard secondary structure. Interestingly, among the known mitogenomes of Pulmonata species, we firstly characterized an unassigned lengthy sequence (551 bp) between the cox1 and the trnV which may be the CR for the sake of its AT bases usage bias (65.70%) and potential hairpin structure.

Physella acuta: atypical mitochondrial gene order among panpulmonates (Gastropoda)

Mitochondrial (mt) sequences are frequently used for phylogenetic reconstruction and for identification of species of molluscs. This study expands the phylogenetic range of Hygrophila (Panpulmonata) for which such sequence data are available by characterizing the full mt genome of the invasive freshwater snail Physella acuta (Physidae). The mt genome sequences of two P. acuta isolates from Stubblefield Lake, New Mexico, USA, differed in length (14,490 vs 14,314 bp) and showed 11.49% sequence divergence, whereas ITS1 and ITS2 sequences from the nuclear genome differed by 1.75%. The mt gene order of P. acuta (cox1, P, nad6, nad5, nad1, D, F, cox2, Y, W, nad4L, C, Q, atp6, R, E, rrnS, M, T, cox3, I, nad2, K, V, rrnL, L1, A, cytb, G, H, L2, atp8, N, nad2, S1, S2, nad4) differs considerably from the relatively conserved gene order within Panpulmonata. Phylogenetic trees show that the 13 protein-encoding mt gene sequences (equivalent codons) of P. acuta group according to gastropod phylogeny, yet branch lengths and dN/dS ratios for P. acuta indicate elevated amino acid substitutions relative to other gastropods. This study indicates that mt sequences of P. acuta are phylogenetically informative despite a considerable intraspecific divergence and the atypical gene order in its mt genome.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

The rRNA and tRNA transcripts of maternally and paternally inherited mitochondrial DNAs of Mytilus galloprovincialis suggest presence of a "degradosome" in mussel mitochondria and necessitate the re-annotation of the l-rRNA/CR boundary

Species of the genus Mytilus carry two mitochondrial genomes in obligatory coexistence; one transmitted though the eggs (the F type) and one through the sperm (the M type). We have studied the 3' and 5' ends of rRNA and tRNA transcripts using RT-PCR and RNA circularization techniques in both the F and M genomes of Mytilus galloprovincialis. We have found polyadenylated and non-adenylated transcripts for both ribosomal and transfer RNAs. In all these genes the 5' ends of the transcripts coincided with the first nucleotide of the annotated genes, but the 3' ends were heterogeneous. The l-rRNA 3' end is 47 or 48 nucleotides upstream from the one assigned by a previous annotation, which makes the adjacent first domain (variable domain one, VD1) of the main control region (CR) correspondingly longer. We have observed s-rRNA and l-rRNA transcripts with truncated 3' end and polyadenylated tRNA transcripts carrying the CCA trinucleotide. We have also detected polyadenylated RNA remnants carrying the sequences of the control region, which strongly suggests RNA degradation activity and thus presence of degradosomes in Mytilus mitochondria.

The complete mitochondrial genome of a turbinid vetigastropod from MiSeq Illumina sequencing of genomic DNA and steps towards a resolved gastropod phylogeny

A need to increase sampling of mitochondrial genomes for Vetigastropoda has been identified as an important step towards resolving relationships within the Gastropoda. We used shotgun sequencing of genomic DNA, using an Illumina MiSeq, to obtain the first mitochondrial genome for the vetigastropod family Turbinidae, doubling the number of genomes for the species-rich superfamily Trochoidea. This method avoids the necessity of finding suitable primers for long PCRs or primer-walking amplicons, resulting in a timely and cost-effective method for obtaining whole mitochondrial genomes from ethanol-preserved tissue samples. Bayesian analysis of amino acid variation for all available gastropod genomes including the new turbinid mtgenome produced a well resolved tree with high nodal support for most nodes. Major clades within Gastropoda were recovered with strong support, with the exception of Littorinimorpha, which was polyphyletic. We confirm here that mitogenomics is a useful tool for molluscan phylogenetics, especially when using powerful new models of amino acid evolution, but recognise that increased taxon sampling is still required to resolve existing differences between nuclear and mitochondrial gene trees.

Complete female mitochondrial genome of Anodonta anatina (Mollusca: Unionidae): confirmation of a novel protein-coding gene (F ORF)

Freshwater mussels are among animals having two different, gender-specific mitochondrial genomes. We sequenced complete female mitochondrial genomes from five individuals of Anodonta anatina, a bivalve species common in palearctic ecozone. The length of the genome was variable: 15,637-15,653 bp. This variation was almost entirely confined to the non-coding parts, which constituted approximately 5% of the genome. Nucleotide diversity was moderate, at 0.3%. Nucleotide composition was typically biased towards AT (66.0%). All genes normally seen in animal mtDNA were identified, as well as the ORF characteristic for unionid mitochondrial genomes, bringing the total number of genes present to 38. If this additional ORF does encode a protein, it must evolve under a very relaxed selection since all substitutions within this gene were non-synonymous. The gene order and structure of the genome were identical to those of all female mitochondrial genomes described in unionid bivalves except the Gonideini.

A comprehensive analysis of bilaterian mitochondrial genomes and phylogeny

About 2800 mitochondrial genomes of Metazoa are present in NCBI RefSeq today, two thirds belonging to vertebrates. Metazoan phylogeny was recently challenged by large scale EST approaches (phylogenomics), stabilizing classical nodes while simultaneously supporting new sister group hypotheses. The use of mitochondrial data in deep phylogeny analyses was often criticized because of high substitution rates on nucleotides, large differences in amino acid substitution rate between taxa, and biases in nucleotide frequencies. Nevertheless, mitochondrial genome data might still be promising as it allows for a larger taxon sampling, while presenting a smaller amount of sequence information. We present the most comprehensive analysis of bilaterian relationships based on mitochondrial genome data. The analyzed data set comprises more than 650 mitochondrial genomes that have been chosen to represent a profound sample of the phylogenetic as well as sequence diversity. The results are based on high quality amino acid alignments obtained from a complete reannotation of the mitogenomic sequences from NCBI RefSeq database. However, the results failed to give support for many otherwise undisputed high-ranking taxa, like Mollusca, Hexapoda, Arthropoda, and suffer from extreme long branches of Nematoda, Platyhelminthes, and some other taxa. In order to identify the sources of misleading phylogenetic signals, we discuss several problems associated with mitochondrial genome data sets, e.g. the nucleotide and amino acid landscapes and a strong correlation of gene rearrangements with long branches.