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

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

SmithHunter: a workflow for the identification of candidate smithRNAs and their targets

BACKGROUND

SmithRNAs (Small MITochondrial Highly-transcribed RNAs) are a novel class of small RNA molecules that are encoded in the mitochondrial genome and regulate the expression of nuclear transcripts. Initial evidence for their existence came from the Manila clam Ruditapes philippinarum, where they have been described and whose activity has been biologically validated through RNA injection experiments. Current evidence on the existence of these RNAs in other species is based only on small RNA sequencing. As a preliminary step to characterize smithRNAs across different metazoan lineages, a dedicated, unified, analytical workflow is needed.

RESULTS

We propose a novel workflow specifically designed for smithRNAs. Sequence data (from small RNA sequencing) uniquely mapping to the mitochondrial genome are clustered into putative smithRNAs and prefiltered based on their abundance, presence in replicate libraries and 5' and 3' transcription boundary conservation. The surviving sequences are subsequently compared to the untranslated regions of nuclear transcripts based on seed pairing, overall match and thermodynamic stability to identify possible targets. Ample collateral information and graphics are produced to help characterize these molecules in the species of choice and guide the operator through the analysis. The workflow was tested on the original Manila clam data. Under basic settings, the results of the original study are largely replicated. The effect of additional parameter customization (clustering threshold, stringency, minimum number of replicates, seed matching) was further evaluated.

CONCLUSIONS

The study of smithRNAs is still in its infancy and no dedicated analytical workflow is currently available. At its core, the SmithHunter workflow builds over the bioinformatic procedure originally applied to identify candidate smithRNAs in the Manila clam. In fact, this is currently the only evidence for smithRNAs that has been biologically validated and, therefore, the elective starting point for characterizing smithRNAs in other species. The original analysis was readapted using current software implementations and some minor issues were solved. Moreover, the workflow was improved by allowing the customization of different analytical parameters, mostly focusing on stringency and the possibility of accounting for a minimal level of genetic differentiation among samples.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Signatures of Extreme Longevity: A Perspective from Bivalve Molecular Evolution

Among Metazoa, bivalves have the highest lifespan disparity, ranging from 1 to 500+ years, making them an exceptional testing ground to understand mechanisms underlying aging and the evolution of extended longevity. Nevertheless, comparative molecular evolution has been an overlooked approach in this instance. Here, we leveraged transcriptomic resources spanning 30 bivalve species to unravel the signatures of convergent molecular evolution in four long-lived species: Margaritifera margaritifera, Elliptio complanata, Lampsilis siliquoidea, and Arctica islandica (the latter represents the longest-lived noncolonial metazoan known so far). We applied a comprehensive approach-which included inference of convergent dN/dS, convergent positive selection, and convergent amino acid substitution-with a strong focus on the reduction of false positives. Genes with convergent evolution in long-lived bivalves show more physical and functional interactions to each other than expected, suggesting that they are biologically connected; this interaction network is enriched in genes for which a role in longevity has been experimentally supported in other species. This suggests that genes in the network are involved in extended longevity in bivalves and, consequently, that the mechanisms underlying extended longevity are-at least partially-shared across Metazoa. Although we believe that an integration of different genes and pathways is required for the extended longevity phenotype, we highlight the potential central roles of genes involved in cell proliferation control, translational machinery, and response to hypoxia, in lifespan extension.

Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes

Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00178-x.

Mitochondrial genomic analyses provide new insights into the "missing" atp8 and adaptive evolution of Mytilidae

Background

Mytilidae, also known as marine mussels, are widely distributed in the oceans worldwide. Members of Mytilidae show a tremendous range of ecological adaptions, from the species distributed in freshwater to those that inhabit in deep-sea. Mitochondria play an important role in energy metabolism, which might contribute to the adaptation of Mytilidae to different environments. In addition, some bivalve species are thought to lack the mitochondrial protein-coding gene ATP synthase F0 subunit 8. Increasing studies indicated that the absence of atp8 may be caused by annotation difficulties for atp8 gene is characterized by highly divergent, variable length.

Results

In this study, the complete mitochondrial genomes of three marine mussels (Xenostrobus securis, Bathymodiolus puteoserpentis, Gigantidas vrijenhoeki) were newly assembled, with the lengths of 14,972 bp, 20,482, and 17,786 bp, respectively. We annotated atp8 in the sequences that we assembled and the sequences lacking atp8. The newly annotated atp8 sequences all have one predicted transmembrane domain, a similar hydropathy profile, as well as the C-terminal region with positively charged amino acids. Furthermore, we reconstructed the phylogenetic trees and performed positive selection analysis. The results showed that the deep-sea bathymodiolines experienced more relaxed evolutionary constraints. And signatures of positive selection were detected in nad4 of Limnoperna fortunei, which may contribute to the survival and/or thriving of this species in freshwater.

Conclusions

Our analysis supported that atp8 may not be missing in the Mytilidae. And our results provided evidence that the mitochondrial genes may contribute to the adaptation of Mytilidae to different environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08940-8.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

Exploring mitogenome evolution in Branchiopoda (Crustacea) lineages reveals gene order rearrangements in Cladocera

The class Branchiopoda, whose origin dates back to Cambrian, includes ~ 1200 species which mainly occupy freshwater habitats. The phylogeny and systematics of the class have been debated for long time, until recent phylogenomic analyses allowed to better clarify the relationships among major clades. Based on these data, the clade Anostraca (fairy and brine shrimps) is sister to all other branchiopods, and the Notostraca (tadpole shrimps) results as sister group to Diplostraca, which includes Laevicaudata + Spinicaudata (clam shrimps) and Cladoceromorpha (water fleas + Cyclestherida). In the present analysis, thanks to an increased taxon sampling, a complex picture emerges. Most of the analyzed mitogenomes show the Pancrustacea gene order while in several other taxa they are found rearranged. These rearrangements, though, occur unevenly among taxa, most of them being found in Cladocera, and their taxonomic distribution does not agree with the phylogeny. Our data also seems to suggest the possibility of potentially homoplastic, alternative gene order within Daphniidae.

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.

Mitogenome sequence comparison in the endangered congeneric Pinna nobilis and Pinna rudis bivalves

BACKGROUND

The pen shells Pinna nobilis and Pinna rudis are large wedge-shaped bivalve molluscs. Both species are threatened by different anthropogenic pressures. In the last few years, P. nobilis populations have significantly reduced due to massive mortality events. The complete mitochondrial DNA sequences of these congeneric species have been determined and compared for the first time.

RESULTS

The mitogenome sequences of P. nobilis and P. rudis were 18,919 bp and 18,264 bp in length, respectively. Each mitogenome is composed of 12 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA (tRNAs) genes and non-coding regions. A putative Adenosine Triphosphate synthase subunit 8 gene could only be proposed for P. nobilis. Both newly sequenced mitogenomes present a conserved gene order between them, comparable to the closely related Atrina pectinata, but global arrangement greatly differs from other available bivalve mitochondrial sequences. Multiple copies of tRNA-Cys were identified, located in different positions probably due to mechanisms of mitochondrial genome rearrangements, and detected 2 and 3 times in P. rudis and in P. nobilis, respectively.

CONCLUSION

A close relationship was shown between Pinna species and Atrina pectinata and a consistent clustering showing a monophyletic origin of Pinnidae family sequences was evidenced. The mitochondrial genomes will provide a valuable genetic resource for further studies on population genetics and species identification.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

The first two complete mitogenomes of the order Apodida from deep-sea chemoautotrophic environments: New insights into the gene rearrangement, origin and evolution of the deep-sea sea cucumbers

The deep-sea ecosystem is considered as the largest and most remote biome of the world. It is meaningful and important to elucidate the life origins by exploring the origin and adaptive genetic mechanisms of the large deep-sea organisms. Sea cucumbers (Holothuroidea) are abundant and economically important group of echinoderms, living from the shallow-waters to deep-sea. In this study, we present the mitochondrial genomes of the sea cucumber Chiridota heheva and Chiridota sp. collected from the deep-sea cold seep and hydrothermal vent, respectively. This is the first reported mitochondrial genomes from the order Apodida. The mitochondrial genomes of C. heheva (17,200 bp) and Chiridota sp. (17,199 bp) display novel gene arrangements with the first protein-coding gene rearrangements in the class Holothuroidea. Bases composition analysis showed that the A + T content of deep-sea holothurians were significantly higher than that of the shallow-water groups. We compared the arrangement of genes from the 24 available holothurian mitogenomes and found that the transposition, reverse transposition and tandem-duplication-random-losses (TDRL) may be involved in the evolution of mitochondrial gene arrangements in Holothuroidea. Phylogenetic analysis revealed that the Apodida clustered with Elasipodida, forming two basal deep-sea holothurian clades. The divergence between the deep-sea and shallow-water holothurians was located at 386.93 Mya, during the Late Devonian. Mitochondrial protein-coding genes of deep-sea holothurians underwent relaxed purifying selection. There are 57 positive selected amino acids sites for some mitochondrial genes of the three deep-sea clades, implying they may involve in the adaption of deep-sea sea cucumbers.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Comparative mitogenomic analysis of the superfamily Tellinoidea (Mollusca: Bivalvia): Insights into the evolution of the gene rearrangements

The superfamily Tellinoidea is widespread and contains approximately 180 living species, which is one of the most diverse and representative groups among the bivalves. In order to extend our knowledge on evolution of tellinoidean species, we newly determined five tellinoidean mitochondrial genomes (mitogenomes). The newly determined mitogenome vary in size from 16,333 to 16,986 bp. The results show that the genome size and genome organization are conserved in tellinoideans. However, gene arrangement and the location of the major non-coding region (NCR) show diversity. The atp8 gene presents in all the five new mitogenomes. Two trnK and trnP genes were detected in Gari togata mitogenome. Phylogenetic analysis supports the monophyly of Tellinoidea, however, it's family Psammobiidae is polyphyletic. CREx analysis suggests that the gene order of Nuttallia olivacea is assumed as the most primitive condition of Tellinoidea. We map the gene order onto the phylogeny and infers the possible gene rearrangement scenarios among tellinoidean mitogenomes. The mitochondrial gene rearrangement is a useful information that help reassessing the phylogeny of Tellinoidea. Phylogenetic relationship and gene arrangement analyses suggest that a careful review for the current taxonomy of the family Psammobiidae is required.

Horizontal transmission and recombination maintain forever young bacterial symbiont genomes

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation.

Clues of in vivo nuclear gene regulation by mitochondrial short non-coding RNAs

Gene expression involves multiple processes, from transcription to translation to the mature, functional peptide, and it is regulated at multiple levels. Small RNA molecules are known to bind RNA messengers affecting their fate in the cytoplasm (a process generically termed ‘RNA interference’). Such small regulatory RNAs are well-known to be originated from the nuclear genome, while the role of mitochondrial genome in RNA interference was largely overlooked. However, evidence is growing that mitochondrial DNA does provide the cell a source of interfering RNAs. Small mitochondrial highly transcribed RNAs (smithRNAs) have been proposed to be transcribed from the mitochondrion and predicted to regulate nuclear genes. Here, for the first time, we show in vivo clues of the activity of two smithRNAs in the Manila clam, Ruditapes philippinarum. Moreover, we show that smithRNAs are present and can be annotated in representatives of the three main bilaterian lineages; in some cases, they were already described and assigned to a small RNA category (e.g., piRNAs) given their biogenesis, while in other cases their biogenesis remains unclear. If mitochondria may affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for them to interact with the nucleus and makes metazoan mitochondrial DNA a much more complex genome than previously thought.

Resolving the relationships of clams and cockles: dense transcriptome sampling drastically improves the bivalve tree of life

Bivalvia has been the subject of extensive recent phylogenetic work to attempt resolving either the backbone of the bivalve tree using transcriptomic data, or the tips using morpho-anatomical data and up to five genetic markers. Yet the first approach lacked decisive taxon sampling and the second failed to resolve many interfamilial relationships, especially within the diverse clade Imparidentia. Here we combine dense taxon sampling with 108 deep-sequenced Illumina-based transcriptomes to provide resolution in nodes that required additional study. We designed specific data matrices to address the poorly resolved relationships within Imparidentia. Our results support the overall backbone of the bivalve tree, the monophyly of Bivalvia and all its main nodes, although the monophyly of Protobranchia remains less clear. Likewise, the inter-relationships of the six main bivalve clades were fully supported. Within Imparidentia, resolution increases when analysing Imparidentia-specific matrices. Lucinidae, Thyasiridae and Gastrochaenida represent three early branches. Gastrochaenida is sister group to all remaining imparidentians, which divide into six orders. Neoheterodontei is always fully supported, and consists of Sphaeriida, Myida and Venerida, with the latter now also containing Mactroidea, Ungulinoidea and Chamidae, a family particularly difficult to place in earlier work. Overall, our study, by using densely sampled transcriptomes, provides the best-resolved bivalve phylogeny to date.

Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria

In animals, strictly maternal inheritance (SMI) of mitochondria is the rule, but one exception (doubly uniparental inheritance or DUI), marked by the transmission of sex-specific mitogenomes, has been reported in bivalves. Associated with DUI is a frequent modification of the mitochondrial cox2 gene, as well as additional sex-specific mitochondrial genes not involved in oxidative phosphorylation. With the exception of freshwater mussels (for 3 families of the order Unionida), these DUI-associated features have only been shown in few species [within Mytilidae (order Mytilida) and Veneridae (order Venerida)] because of the few complete sex-specific mitogenomes published for these orders. Here, we present the complete sex-specific mtDNAs of two recently-discovered DUI species in two families of the order Venerida, Scrobicularia plana (Semelidae) and Limecola balthica (Tellinidae). These species display the largest differences in genome size between sex-specific mitotypes in DUI species (>10 kb), as well as the highest mtDNA divergences (sometimes reaching >50%). An important in-frame insertion (>3.5 kb) in the male cox2 gene is partly responsible for the differences in genome size. The S. plana cox2 gene is the largest reported so far in the Kingdom Animalia. The mitogenomes may be carrying sex-specific genes, indicating that general mitochondrial features are shared among DUI species.

Linking paternally inherited mtDNA variants and sperm performance

Providing robust links between mitochondrial genotype and phenotype is of major importance given that mitochondrial DNA (mtDNA) variants can affect reproductive success. Because of the strict maternal inheritance (SMI) of mitochondria in animals, haplotypes that negatively affect male fertility can become fixed in populations. This phenomenon is known as 'mother's curse'. Doubly uniparental inheritance (DUI) of mitochondria is a stable exception in bivalves, which entails two mtDNA lineages that evolve independently and are transmitted separately through oocytes and sperm. This makes the DUI mitochondrial lineages subject to different sex-specific selective sieves during mtDNA evolution, thus DUI is a unique model to evaluate how direct selection on sperm mitochondria could contribute to male reproductive fitness. In this study, we tested the impact of mtDNA variants on sperm performance and bioenergetics in DUI and SMI species. Analyses also involved measures of sperm performance following inhibition of main energy pathways and sperm response to oocyte presence. Compared to SMI, DUI sperm exhibited (i) low speed and linearity, (ii) a strict OXPHOS-dependent strategy of energy production, and (iii) a partial metabolic shift towards fermentation following egg detection. Discussion embraces the adaptive value of mtDNA variation and suggests a link between male-energetic adaptation, fertilization success and paternal mitochondria preservation. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Sperm competitive advantage of a rare mitochondrial haplogroup linked to differential expression of mitochondrial oxidative phosphorylation genes

Maternal inheritance of mitochondria creates a sex-specific selective sieve through which mitochondrial mutations harmful to males but not females accumulate and contribute to sexual differences in longevity and disease susceptibility. Because eggs and sperm are under disruptive selection, sperm are predicted to be particularly vulnerable to the genetic load generated by maternal inheritance, yet evidence for mitochondrial involvement in male fertility is limited and controversial. Here, we exploit the coexistence of two divergent mitochondrial haplogroups (A and B2) in a Neotropical arachnid to investigate the role of mitochondria in sperm competition. DNA profiling demonstrated that B2-carrying males sired more than three times as many offspring in sperm competition experiments than A males, and this B2 competitive advantage cannot be explained by female mitochondrial haplogroup or male nuclear genetic background. RNA-Seq of testicular tissues implicates differential expression of mitochondrial oxidative phosphorylation (OXPHOS) genes in the B2 competitive advantage, including a 22-fold upregulation of atp8 in B2 males. Previous comparative genomic analyses have revealed functionally significant amino acid substitutions in differentially expressed genes, indicating that the mitochondrial haplogroups differ not only in expression but also in DNA sequence and protein functioning. However, mitochondrial haplogroup had no effect on sperm number or sperm viability, and, when females were mated to a single male, neither male haplogroup, female haplogroup nor the interaction between male/female haplogroup significantly affected female reproductive success. Our findings therefore suggest that mitochondrial effects on male reproduction may often go undetected in noncompetitive contexts and may prove more important in nature than is currently appreciated.

Arm-less mitochondrial tRNAs conserved for over 30 millions of years in spiders

BACKGROUND

In recent years, Next Generation Sequencing (NGS) has accelerated the generation of full mitogenomes, providing abundant material for studying different aspects of molecular evolution. Some mitogenomes have been observed to harbor atypical sequences with bizarre secondary structures, which origins and significance could only be fully understood in an evolutionary framework.

RESULTS

Here we report and analyze the mitochondrial sequences and gene arrangements of six closely related spiders in the sister genera Parachtes and Harpactocrates, which belong to the nocturnal, ground dwelling family Dysderidae. Species of both genera have compacted mitogenomes with many overlapping genes and strikingly reduced tRNAs that are among the shortest described within metazoans. Thanks to the conservation of the gene order and the nucleotide identity across close relatives, we were able to predict the secondary structures even on arm-less tRNAs, which would be otherwise unattainable for a single species. They exhibit aberrant secondary structures with the lack of either DHU or TΨC arms and many miss-pairings in the acceptor arm but this degeneracy trend goes even further since at least four tRNAs are arm-less in the six spider species studied.

CONCLUSIONS

The conservation of at least four arm-less tRNA genes in two sister spider genera for about 30 myr suggest that these genes are still encoding fully functional tRNAs though they may be post-transcriptionally edited to be fully functional as previously described in other species. We suggest that the presence of overlapping and truncated tRNA genes may be related and explains why spider mitogenomes are smaller than those of other invertebrates.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Comparative analysis of the mitochondrial genomes of Orthonectida: insights into the evolution of an invertebrate parasite species

Among invertebrates, only a few groups still have uncertain phylogenetic position, Orthonectida, a small group of rare multi-cellular parasites of marine invertebrates, being one of them. Recent molecular and morphological findings suggest that orthonectids belong to Lophotrochozoa and are close to Annelida. Nevertheless, phylogenetic relationships between orthonectids and annelids are unclear, and the phylogeny within the group itself has never been studied. Sequencing of mitochondrial genomes is used here to clarify this issue. Complete mt genomes of the orthonectids Intoshia variabili and Rhopalura litoralis were characterized and compared with Intoshia linei mt genome. Our results show that Orthonectida mt genomes have undergone reduction and gene loss, and that they have complicated organization revealed in strand asymmetry in nucleotide composition, in some features of intergenic non-coding regions, tRNA duplication and folding. Moreover, all species of Orthonectida have a unique gene order with complicated rearrangement landscape. Significant differences in mitochondrial genomes in the three orthonectid species could be explained by the fact that their host species belong to different taxa (flat worms, nemertines and gastropods). Among the analyzed mt genomes of Orthonectida, I. linei possesses the closest gene order to the ancestral genome. All Orthonectida species are monophyletic, and in the phylogenetic tree are close to Pleistoannelida, and specifically, to Clitellata.

Relation between mitochondrial DNA hyperdiversity, mutation rate and mitochondrial genome evolution in Melarhaphe neritoides (Gastropoda: Littorinidae) and other Caenogastropoda

Mitochondrial DNA hyperdiversity is primarily caused by high mutation rates (µ) and has potential implications for mitogenome architecture and evolution. In the hyperdiverse mtDNA of Melarhaphe neritoides (Gastropoda: Littorinidae), high mutational pressure generates unusually large amounts of synonymous variation, which is expected to (1) promote changes in synonymous codon usage, (2) reflect selection at synonymous sites, (3) increase mtDNA recombination and gene rearrangement, and (4) be correlated with high mtDNA substitution rates. The mitogenome of M. neritoides was sequenced, compared to closely related littorinids and put in the phylogenetic context of Caenogastropoda, to assess the influence of mtDNA hyperdiversity and high µ on gene content and gene order. Most mitogenome features are in line with the trend in Mollusca, except for the atypical secondary structure of the methionine transfer RNA lacking the TΨC-loop. Therefore, mtDNA hyperdiversity and high µ in M. neritoides do not seem to affect its mitogenome architecture. Synonymous sites are under positive selection, which adds to the growing evidence of non-neutral evolution at synonymous sites. Under such non-neutrality, substitution rate involves neutral and non-neutral substitutions, and high µ is not necessarily associated with high substitution rate, thus explaining that, unlike high µ, a high substitution rate is associated with gene order rearrangement.

Structure-Related Differences between Cytochrome Oxidase I Proteins in a Stable Heteroplasmic Mitochondrial System

Many bivalve species have two types of mitochondrial DNA passed independently through the female line (F genome) and male line (M genome). Here we study the cytochrome oxidase I protein in such bivalve species and provide evidence for differences between the F and M proteins in amino acid property values, particularly relating to hydrophobicity and helicity. The magnitude of these differences varies between different regions of the protein and the change from the ancestor is most marked in the M protein. The observed changes occur in parallel and in the same direction in the different species studied. Two possible causes are considered, first relaxation of purifying selection with drift and second positive selection. These may operate in different ways in different regions of the protein. Many different amino acid substitutions contribute in a small way to the observed variation, but substitutions involving alanine and serine have a quantitatively large effect. Some of these substitutions are potential targets for phosphorylation and some are close to residues of functional importance in the catalytic mechanism. We propose that the observed changes in the F and M proteins might contribute to functional differences between them relating to ATP production and mitochondrial membrane potential with implications for sperm function.

Comparative analyses of the complete mitochondrial genomes of Dosinia clams and their phylogenetic position within Veneridae

Mitochondrial genomes have proved to be a powerful tool in resolving phylogenetic relationship. In order to understand the mitogenome characteristics and phylogenetic position of the genus Dosinia, we sequenced the complete mitochondrial genomes of Dosinia altior and Dosinia troscheli (Bivalvia: Veneridae), compared them with that of Dosinia japonica and established a phylogenetic tree for Veneridae. The mitogenomes of D. altior (17,536 bp) and D. troscheli (17,229 bp) are the two smallest in Veneridae, which include 13 protein-coding genes, 2 ribosomal RNA genes, 22 tRNA genes, and non-coding regions. The mitogenomes of the Dosinia species are similar in size, gene content, AT content, AT- and GC- skews, and gene arrangement. The phylogenetic relationships of family Veneridae were established based on 12 concatenated protein-coding genes using maximum likelihood and Bayesian analyses, which supported that Dosininae and Meretricinae have a closer relationship, with Tapetinae being the sister taxon. The information obtained in this study will contribute to further understanding of the molecular features of bivalve mitogenomes and the evolutionary history of the genus Dosinia.

Reclassification of Lamprotula rochechouartii as Margaritifera rochechouartiicomb. nov. (Bivalvia: Margaritiferidae) revealed by time-calibrated multi-locus phylogenetic analyses and mitochondrial phylogenomics of Unionoida

The family Margaritiferidae encompasses 12 valid species, which are distributed widely but disjunctively in the Northern Hemisphere. A lack of a well resolved and temporally calibrated phylogenetic framework of Margaritiferidae has made it difficult to discuss the evolutionary pattern and process. Phylogenetic relationships between five major clades, which were revealed in earlier studies, remain elusive and unresolved. Lamprotula rochechouartii has long been classified within the family Unionidae based on shell morphology, but our preliminary molecular study on this species made us hypothesize that it has an affinity with margaritiferids. Hence, five loci (COI, 16S, 18S, 28S and histone H3) were used to investigate the phylogenetic position of L. rochechouartii and intra-familial relationships within Margaritiferidae using various partitioning strategies. Moreover, two mitochondrial genomes were newly obtained to further resolve and validate the five-clade relationships within Margaritiferidae in a broad view of Unionoida evolution. Both five-gene and mitogenome datasets strongly advocated treating Lamprotula rochechouartii as Margaritifera rochechouartiicomb. nov. Maximum likelihood and Bayesian inference analyses using partitioned five-gene dataset resulted in various topologies, but five well-supported clades were obtained. The most probable cladistic relationships generated by five-gene dataset analyses were identical to subsequent whole mitogenome analyses except the position of M. monodonta. M. rochechouartii and M. laosensis had a well-supported sister relationship and formed a basal clade splitting from the rest of the family. Based on six reliable fossils, crown age of the extant Margaritiferidae was estimated during the Late Cretaceous at 88.3 Ma (95% HPD = 66.2-117.4). But we hypothesized a much earlier origin of this family due to the Permian stem age (mean = 257 Ma, 95% HPD = 230.0-296.0) and a high extinction rate in the whole order. Biogeographic scenarios supported a Laurasian origin of extant Margaritiferidae during the Late Cretaceous, and suggested that Asian margaritiferids may have had two origins, having either Asia (M. rochechouartii, M. laosensis) or North America (M. dahurica, M. laevis, and M. middendorffi) as ancestral. The newly added Margaritiferidae species M. rochechouartii expands our recognized distribution range of modern margaritiferids. Our results indicate that whole mitogenome sequences can be used to reconstruct robust phylogenetic relationships for freshwater mussels, especially with the help of adding M-type mitogenomes.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

The complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus (Bivalvia, Veneridae)

Despite the large number of animal complete mitochondrial genomes currently available in public databases, knowledge about mitochondrial genomics in invertebrates is uneven. This paper reports, for the first time, the complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus, also known as the European clam. Ruditapes decussatus is morphologically and ecologically similar to the Manila clam Ruditapes philippinarum, which has been recently introduced for aquaculture in the very same habitats of Ruditapes decussatus, and that is replacing the native species. Currently the production of the European clam is almost insignificant, nonetheless it is considered a high value product, and therefore it is an economically important species, especially in Portugal, Spain and Italy. In this work we: (i) assembled Ruditapes decussatus mitochondrial genome from RNA-Seq data, and validated it by Sanger sequencing; (ii) analyzed and characterized the Ruditapes decussatus mitochondrial genome, comparing its features with those of other venerid bivalves; (iii) assessed mitochondrial sequence polymorphism (SP) and copy number variation (CNV) of tandem repeats across 26 samples. Despite using high-throughput approaches we did not find evidence for the presence of two sex-linked mitochondrial genomes, typical of the doubly uniparental inheritance of mitochondria, a phenomenon known in ∼100 bivalve species. According to our analyses, Ruditapes decussatus is more genetically similar to species of the Genus Paphia than to the congeneric Ruditapes philippinarum, a finding that bolsters the already-proposed need of a taxonomic revision. We also found a quite low genetic variability across the examined samples, with few SPs and little variability of the sequences flanking the control region (Largest Unassigned Regions (LURs). Strikingly, although we found low nucleotide variability along the entire mitochondrial genome, we observed high levels of length polymorphism in the LUR due to CNV of tandem repeats, and even a LUR length heteroplasmy in two samples. It is not clear if the lack of genetic variability in the mitochondrial genome of Ruditapes decussatus is a cause or an effect of the ongoing replacement of Ruditapes decussatus with the invasive Ruditapes philippinarum, and more analyses, especially on nuclear sequences, are required to assess this point.

SmithRNAs: Could Mitochondria "Bend" Nuclear Regulation?

Typically, animal mitochondria have very compact genomes, with few short intergenic regions, and no introns. Hence, it may seem that there is little space for unknown functions in mitochondrial DNA (mtDNA). However, mtDNA can also operate through RNA interference, as small non coding RNAs (sncRNAs) produced by mtDNA have already been proposed for humans. We sequenced sncRNA libraries from isolated mitochondria of Ruditapes philippinarum (Mollusca Bivalvia) gonads, a species with doubly uniparental inheritance of mitochondria, and identified several putative sncRNAs of mitochondrial origin. Some sncRNAs are transcribed by intergenic regions that form stable stem-hairpin structures, which makes them good miRNA-like candidates. We decided to name them small mitochondrial highly-transcribed RNAs (smithRNAs). Many concurrent data support that we have recovered sncRNAs of mitochondrial origin that might be involved in gonad formation and able to affect nuclear gene expression. This possibility has been never suggested before. If mtDNA can affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for it to interact with the nucleus, and makes metazoan mtDNA a much more complex genome than previously thought.

Evolution of Rosaceae Fruit Types Based on Nuclear Phylogeny in the Context of Geological Times and Genome Duplication

Fruits are the defining feature of angiosperms, likely have contributed to angiosperm successes by protecting and dispersing seeds, and provide foods to humans and other animals, with many morphological types and important ecological and agricultural implications. Rosaceae is a family with ∼3000 species and an extraordinary spectrum of distinct fruits, including fleshy peach, apple, and strawberry prized by their consumers, as well as dry achenetum and follicetum with features facilitating seed dispersal, excellent for studying fruit evolution. To address Rosaceae fruit evolution and other questions, we generated 125 new transcriptomic and genomic datasets and identified hundreds of nuclear genes to reconstruct a well-resolved Rosaceae phylogeny with highly supported monophyly of all subfamilies and tribes. Molecular clock analysis revealed an estimated age of ∼101.6 Ma for crown Rosaceae and divergence times of tribes and genera, providing a geological and climate context for fruit evolution. Phylogenomic analysis yielded strong evidence for numerous whole genome duplications (WGDs), supporting the hypothesis that the apple tribe had a WGD and revealing another one shared by fleshy fruit-bearing members of this tribe, with moderate support for WGDs in the peach tribe and other groups. Ancestral character reconstruction for fruit types supports independent origins of fleshy fruits from dry-fruit ancestors, including the evolution of drupes (e.g., peach) and pomes (e.g., apple) from follicetum, and drupetum (raspberry and blackberry) from achenetum. We propose that WGDs and environmental factors, including animals, contributed to the evolution of the many fruits in Rosaceae, which provide a foundation for understanding fruit evolution.

Evolution of sex-dependent mtDNA transmission in freshwater mussels (Bivalvia: Unionida)

Doubly uniparental inheritance (DUI) describes a mode of mtDNA transmission widespread in gonochoric freshwater mussels (Bivalvia: Palaeoheterodonta: Unionida). In this system, both female- and male-transmitted mtDNAs, named F and M respectively, coexist in the same species. In unionids, DUI is strictly correlated to gonochorism and to the presence of the atypical open reading frames (ORFans) F-orf and M-orf, respectively inside F and M mtDNAs, which are hypothesized to participate in sex determination. However, DUI is not found in all three Unionida superfamilies (confirmed in Hyrioidea and Unionoidea but not in Etherioidea), raising the question of its origin in these bivalves. To reconstruct the co-evolution of DUI and of ORFans, we sequenced the mtDNAs of four unionids (two gonochoric with DUI, one gonochoric and one hermaphroditic without DUI) and of the related gonochoric species Neotrigonia margaritacea (Palaeoheterodonta: Trigoniida). Our analyses suggest that rearranged mtDNAs appeared early during unionid radiation, and that a duplicated and diverged atp8 gene evolved into the M-orf associated with the paternal transmission route in Hyrioidea and Unionoidea, but not in Etherioidea. We propose that novel mtDNA-encoded genes can deeply influence bivalve sex determining systems and the evolution of the mitogenomes in which they occur.

Burrowers from the Past: Mitochondrial Signatures of Ordovician Bivalve Infaunalization

Bivalves and gastropods are the two largest classes of extant molluscs. Despite sharing a huge number of features, they do not share a key ecological one: gastropods are essentially epibenthic, although most bivalves are infaunal. However, this is not the ancestral bivalve condition; Cambrian forms were surface crawlers and only during the Ordovician a fundamental infaunalization process took place, leading to bivalves as we currently know them. This major ecological shift is linked to the exposure to a different redox environoments (hypoxic or anoxic) and with the Lower Devonian oxygenation event. We investigated selective signatures on bivalve and gastropod mitochondrial genomes with respect to a time calibrated mitochondrial phylogeny by means of dN/dS ratios. We were able to detect 1) a major signal of directional selection between the Ordovician and the Lower Devonian for bivalve mitochondrial Complex I, and 2) an overall higher directional selective pressure on bivalve Complex V with respect to gastropods. These and other minor dN/dS patterns and timings are discussed, showing that the Ordovician infaunalization event left heavy traces in bivalve mitochondrial genomes.