Doubly uniparental inheritance (DUI) of mitochondrial DNA in Donax trunculus (Bivalvia: Donacidae) and the problem of its sporadic detection in Bivalvia

The absence of canonical respiratory complex I subunits in male-type mitogenomes of three Donax species

Bivalves are an extraordinary class of animals in which species with a doubly uniparental inheritance (DUI) of mitochondrial DNA have been described. DUI is characterized as a mitochondrial homoplasmy of females and heteroplasmy of male individuals where F-type mitogenomes are passed to the progeny with mother egg cells and divergent M-type mitogenomes are inherited with fathers sperm cells. However, in most cases only male individuals retain divergent mitogenome inherited with spermatozoa. Additionally, in many of bivalves, unique mitochondrial features, like additional genes, gene duplication, gene extensions, mitochondrial introns, and recombination, were observed. In this study, we sequenced and assembled male-type mitogenomes of three Donax species. Comparative analysis of mitochondrial sequences revealed a lack of all seven NADH dehydrogenase subunits as well as the presence of three long additional open reading frames lacking identifiable homology to any of the existing genes.

Insights on Pinna nobilis population genetic structure in the Aegean and Ionian Sea

The fan mussel Pinna nobilis Linnaeus, 1758 is an endemic species of the Mediterranean Sea, protected by international agreements. It is one of the largest bivalves in the world, playing an important role in the benthic communities; yet it has been recently characterized as Critically Endangered by the IUCN, due to mass mortality events. In this context, the assessment of the genetic variation of the remaining P. nobilis populations and the evaluation of connectivity among them are crucial elements for the conservation of the species. For this purpose, samples were collected from six regions of the Eastern Mediterranean Sea; the Islands of Karpathos, Lesvos and Crete; the Chalkidiki and Attica Peninsulas; and the Amvrakikos Gulf. Sampling was performed either by collecting tissue from the individuals or by using a non-invasive method, i.e., by scraping the inside of their shells aiming to collect their mucus and thus avoid stress induction to them. Conventional molecular techniques with the use of the COI and 16S rRNA mitochondrial markers were selected for the depiction of the intra-population genetic variability. The analyses included 104 samples from the present study and publicly available sequences of individuals across the whole Mediterranean Sea. The results of this work (a) suggest the use of eDNA as an efficient sampling method for protected bivalves and (b) shed light to the genetic structure of P. nobilis population in the Eastern Mediterranean; this latter knowledge might prove to be fundamental for the species conservation and hence the ecosystem resilience. The haplotype analyses reinforced the evidence that there is a certain degree of connectivity among the distinct regions of the Mediterranean; yet there is evidence of population distinction within the basin, namely between the Western and the Eastern basins. The combination of both genetic markers in the same analysis along with the inclusion of a large number of individuals produced more robust results, revealing a group of haplotypes being present only in the Eastern Mediterranean and providing insights for the species' most suitable conservation management.

A tale of two paths: The evolution of mitochondrial recombination in bivalves with doubly uniparental inheritance

In most animals, mitochondrial DNA is strictly maternally inherited and non-recombining. One exception to this pattern is called doubly uniparental inheritance (DUI), a phenomenon involving the independent transmission of female and male mitochondrial genomes. DUI is known only from the molluskan class Bivalvia. The phylogenetic distribution of male-transmitted mitochondrial DNA (M mtDNA) in bivalves is consistent with several evolutionary scenarios, including multiple independent gains, losses, and varying degrees of recombination with female-transmitted mitochondrial DNA (F mtDNA). In this study, we use phylogenetic methods to test M mtDNA origination hypotheses and infer the prevalence of mitochondrial recombination in bivalves with DUI. Phylogenetic modeling using site concordance factors supported a single origin of M mtDNA in bivalves coupled with recombination acting over long evolutionary timescales. Ongoing mitochondrial recombination is present in Mytilida and Venerida, which results in a pattern of concerted evolution of F mtDNA and M mtDNA. Mitochondrial recombination could be favored to offset the deleterious effects of asexual inheritance and maintain mitonuclear compatibility across tissues. Cardiida and Unionida have gone without recent recombination, possibly due to an extension of the COX2 gene in male mitochondrial DNA. The loss of recombination could be connected to the role of M mtDNA in sex determination or sexual development. Our results support that recombination events may occur throughout the mitochondrial genomes of DUI species. Future investigations may reveal more complex patterns of inheritance of recombinants, which could explain the retention of signal for a single origination of M mtDNA in protein-coding genes.

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.

No evidence of DUI in the Mediterranean alien species Brachidontes pharaonis (P. Fisher, 1870) despite mitochondrial heteroplasmy

Two genetically different mitochondrial haplogroups of Brachidontes pharaonis (p-distance 6.8%) have been identified in the Mediterranean Sea. This hinted at a possible presence of doubly uniparental inheritance in this species. To ascertain this possibility, we sequenced two complete mitogenomes of Brachidontes pharaonis mussels and performed a qPCR analysis to measure the relative mitogenome copy numbers of both mtDNAs. Despite the presence of two very similar regions composed entirely of repetitive sequences in the two haplogroups, no recombination between mitogenomes was detected. In heteroplasmic individuals, both mitogenomes were present in the generative tissues of both sexes, which argues against the presence of doubly uniparental inheritance in this species.

The complete mitochondrial genome of Mactra chinensis (Bivalvia: Macridae)

The structure and composition of the mitogenome of a bivalve mollusk, denoted as Mactra sp. (MT780813), has been obtained. The genome has a variable organization: it includes 12 protein-coding genes, 28 tRNAs, and two rRNA genes. Its content is sufficiently different from that of nearest specimen, accessed from GenBank, supposedly belonging to the other gender. All genes are encoded on the "+"-strand. All protein-coding genes are initiated with ATG codon. Analysis confirms the close topological position of the GenBank Mactra chinensis (KJ754823) and our M. sp. specimen on gene tree. Above data suggesting female- vs. male-type mitogenomes or cryptic species presence.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida)

BACKGROUND

Supernumerary ORFan genes (i.e., open reading frames without obvious homology to other genes) are present in the mitochondrial genomes of gonochoric freshwater mussels (Bivalvia: Unionida) showing doubly uniparental inheritance (DUI) of mitochondria. DUI is a system in which distinct female-transmitted and male-transmitted mitotypes coexist in a single species. In families Unionidae and Margaritiferidae, the transition from dioecy to hermaphroditism and the loss of DUI appear to be linked, and this event seems to affect the integrity of the ORFan genes. These observations led to the hypothesis that the ORFans have a role in DUI and/or sex determination. Complete mitochondrial genome sequences are however scarce for most families of freshwater mussels, therefore hindering a clear localization of DUI in the various lineages and a comprehensive understanding of the influence of the ORFans on DUI and sexual systems. Therefore, we sequenced and characterized eleven new mitogenomes from poorly sampled freshwater mussel families to gather information on the evolution and variability of the ORFan genes and their protein products.

RESULTS

We obtained ten complete plus one almost complete mitogenome sequence from ten representative species (gonochoric and hermaphroditic) of families Margaritiferidae, Hyriidae, Mulleriidae, and Iridinidae. ORFan genes are present only in DUI species from Margaritiferidae and Hyriidae, while non-DUI species from Hyriidae, Iridinidae, and Mulleriidae lack them completely, independently of their sexual system. Comparisons among the proteins translated from the newly characterized ORFans and already known ones provide evidence of conserved structures, as well as family-specific features.

CONCLUSIONS

The ORFan proteins show a comparable organization of secondary structures among different families of freshwater mussels, which supports a conserved physiological role, but also have distinctive family-specific features. Given this latter observation and the fact that the ORFans can be either highly mutated or completely absent in species that secondarily lost DUI depending on their respective family, we hypothesize that some aspects of the connection among ORFans, sexual systems, and DUI may differ in the various lineages of unionids.

Sequence characterization of the 5S ribosomal DNA and the internal transcribed spacer (ITS) region in four European Donax species (Bivalvia: Donacidae)

Background

The whole repeat unit of 5S rDNA and the internal transcribed spacer (ITS) of four European Donax species were analysed. After amplifying, cloning and sequencing several 5S and ITS units, their basic features and their variation were described. The phylogenetic usefulness of 5S and ITS sequences in the inference of evolutionary relationships among these wedge clams was also investigated.

Results

The length of the 5S repeat presented little variation among species, except D. trunculus that differed from the rest of the Donax species in 170–210 bp. The deduced coding region covered 120 bp, and showed recognizable internal control regions (ICRs) involved in the transcription. The length of non-transcribed spacer region (NTS) ranged from 157 bp to 165 bp in Donax trunculus and from 335 bp to 367 bp in the other three species. The conservation degree of transcriptional regulatory regions was analysed revealing a conserved TATA-like box in the upstream region. Regarding ITS sequences, the four Donax species showed slight size differences among clones due to the variation occurring in the ITS1 and ITS2, except Donax variegatus did not display size differences in the ITS2. The total length of the ITS sequence ranged between 814 and 1014 bp. Resulting phylogenetic trees display that the two ribosomal DNA regions provide well-resolved phylogenies where the four European Donax species form a single clade receiving high support in nodes. The topology obtained with 5S sequences was in agreement with Donax evolutionary relationships inferred from several sequences of different nature in previous studies.

Conclusions

This is not only a basic research work, where new data and new knowledge is provided about Donax species, but also have allowed the authentication of these wedge clams and offers future applications to provide other genetic resources.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0684-x) contains supplementary material, which is available to authorized users.

Lose it or keep it: (how bivalves can provide) insights into mitochondrial inheritance mechanisms

The strictly maternal inheritance (SMI) is a pattern of mitochondrial inheritance observed across the whole animal kingdom. However, some interesting exceptions are known for the class Bivalvia, in which several species show an unusual pattern called doubly uniparental inheritance (DUI) whose outcome is a heteroplasmic pool of mtDNA in males. Even if DUI has been studied for long, its molecular basis has not been established yet. The aim of this work is to select classes of proteins known to be involved in the maintenance of SMI and to compare their features in two clam species differing for their mitochondrial inheritance mechanism, that is, the SMI species Ruditapes decussatus and the DUI species Ruditapes philippinarum. Data have been obtained from the transcriptomes of male and female ripe gonads of both species. Our analysis focused on nucleases and polymerases, ubiquitination and ubiquitin-like modifier pathways, and proteins involved in autophagy and mitophagy. For each protein group of interest, transcription bias (male or female), annotation, and mitochondrial targeting (when appropriate) were assessed. We did not find evidence supporting a role of nucleases/polymerases or autophagic machinery in the enforcement of SMI in R. decussatus. On the other hand, ubiquitinating enzymes with the expected features have been retrieved, providing us with two alternative testable models for mitochondrial inheritance mechanisms at the molecular level.

Gender-Associated Mitochondrial DNA Heteroplasmy in Somatic Tissues of the Endangered Freshwater Mussel Unio crassus (Bivalvia: Unionidae): Implications for Sex Identification and Phylogeographical Studies

Some bivalve species possess two independent mitochondrial DNA lineages: maternally (F-type) and paternally (M-type) inherited. This phenomenon is called doubly uniparental inheritance. It is generally agreed that F-type mtDNA is typically present in female somatic and gonadal tissues as well as in male somatic tissues, whereas the M-type mtDNA occurs only in male germ line and gonadal tissue. In the present study, the mtDNA heteroplasmy (for both F and M genomes) in male somatic tissues of Unio crassus (Philipsson, 1788), species threatened with extinction, has been confirmed. Taking advantage from the presence of Mcox1 marker only in male somatic tissues, we developed a new method of sex identification in this endangered species, using nondestructive tissue sampling. Probability of correct sex identification was estimated at 97.5%. The present study is the first report on gender-associated mitochondrial DNA heteroplasmy in male somatic tissues of thick-shelled river mussel and first approach to U. crassus sex identification at molecular level. Our study also confirmed the utility of paternally inherited Mcox1 gene fragment as a complementary molecular tool for resolving phylogeographical relationships among populations of thick-shelled river mussel.

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.

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.

Implications for management and conservation of the population genetic structure of the wedge clam Donax trunculus across two biogeographic boundaries

In a resource management perspective, the understanding of the relative influence of the physical factors on species connectivity remains a major challenge and is also of great ecological and conservation biology interest. Despite the overfishing threat on the wedge clam Donax trunculus in Europe, relatively little information is known about its population genetic structure and connectivity and their consequences on conservation policies. We employed 16 microsatellite loci to characterise the genetic diversity and population structure of D. trunculus. A total of 514 samples from seven different localities along the Atlantic-Mediterranean transition, from the Atlantic (Gulf of Cádiz) to the north-western Mediterranean were genotyped. The analysis of the population genetic structure displayed a clear distinction along the Atlantic-Mediterranean transition with different clusters in the Atlantic Ocean, the Alboran Sea and the northwestern Mediterranean. Consequently, we recommend that these three areas should be considered as different management units. We showed that all populations seem to be at high long-term risk of extinction with the exception of the protected Doñana National Park population which still seems to have evolutionary potential. Therefore, our results emphasized the necessity of protection of this economic resource and the validity of molecular tools to evaluate the population dynamics.

Pursuing the quest for better understanding the taxonomic distribution of the system of doubly uniparental inheritance of mtDNA

There is only one exception to strict maternal inheritance of mitochondrial DNA (mtDNA) in the animal kingdom: a system named doubly uniparental inheritance (DUI), which is found in several bivalve species. Why and how such a radically different system of mitochondrial transmission evolved in bivalve remains obscure. Obtaining a more complete taxonomic distribution of DUI in the Bivalvia may help to better understand its origin and function. In this study we provide evidence for the presence of sex-linked heteroplasmy (thus the possible presence of DUI) in two bivalve species, i.e., the nuculanoid Yoldia hyperborea(Gould, 1841)and the veneroid Scrobicularia plana(Da Costa,1778), increasing the number of families in which DUI has been found by two. An update on the taxonomic distribution of DUI in the Bivalvia is also presented.

From writing to reading the encyclopedia of life

Prologue 'As the study of natural science advances, the language of scientific description may be greatly simplified and abridged. This has already been done by Linneaus and may be carried still further by other invention. The descriptions of natural orders and genera may be reduced to short definitions, and employment of signs, somewhat in the manner of algebra, instead of long descriptions. It is more easy to conceive this, than it is to conceive with what facility, and in how short a time, a knowledge of all the objects of natural history may ultimately be acquired; and that which is now considered learning and science, and confined to a few specially devoted to it, may at length be universally possessed in every civilized country and in every rank of life'. J. C. Louden 1829. Magazine of natural history, vol. 1: This article is part of the themed issue 'From DNA barcodes to biomes'.

The Complete Female- and Male-Transmitted Mitochondrial Genome of Meretrix lamarckii

Bivalve mitochondrial genomes show many uncommon features, like additional genes, high rates of gene rearrangement, high A-T content. Moreover, Doubly Uniparental Inheritance (DUI) is a distinctive inheritance mechanism allowing some bivalves to maintain and transmit two separate sex-linked mitochondrial genomes. Many bivalve mitochondrial features, such as gene extensions or additional ORFs, have been proposed to be related to DUI but, up to now, this topic is far from being understood. Several species are known to show this unusual organelle inheritance but, being widespread only among Unionidae and Mytilidae, DUI distribution is unclear. We sequenced and characterized the complete female- (F) and male-transmitted (M) mitochondrial genomes of Meretrix lamarckii, which, in fact, is the second species of the family Veneridae where DUI has been demonstrated so far. The two mitochondrial genomes are comparable in length and show roughly the same gene content and order, except for three additional tRNAs found in the M one. The two sex-linked genomes show an average nucleotide divergence of 16%. A 100-aminoacid insertion in M. lamarckii M-cox2 gene was found; moreover, additional ORFs have been found in both F and M Long Unassigned Regions of M. lamarckii. Even if no direct involvement in DUI process has been demonstrated so far, the finding of cox2 insertions and supernumerary ORFs in M. lamarckii both strengthens this hypothesis and widens the taxonomical distribution of such unusual features. Finally, the analysis of inter-sex genetic variability shows that DUI species form two separate clusters, namely Unionidae and Mytilidae+Veneridae; this dichotomy is probably due to different DUI regimes acting on separate taxa.

Early replication dynamics of sex-linked mitochondrial DNAs in the doubly uniparental inheritance species Ruditapes philippinarum (Bivalvia Veneridae)

Mitochondrial homoplasmy, which is maintained by strictly maternal inheritance and a series of bottlenecks, is thought to be an adaptive condition for metazoans. Doubly uniparental inheritance (DUI) is a unique mode of mitochondrial transmission found in bivalve species, in which two distinct mitochondrial genome (mtDNA) lines are present, one inherited through eggs (F) and one through sperm (M). During development, the two lines segregate in a sex- and tissue-specific manner: females lose M during embryogenesis, whereas males actively segregate it in the germ line. These two pivotal events are still poorly characterized. Here we investigated mtDNA replication dynamics during embryogenesis and pre-adulthood of the venerid Ruditapes philippinarum using real-time quantitative PCR. We found that both mtDNAs do not detectably replicate during early embryogenesis, and that the M line might be lost from females around 24 h of age. A rise in mtDNA copy number was observed before the first reproductive season in both sexes, with the M mitochondrial genome replicating more than the F in males, and we associate these boosts to the early phase of gonad production. As evidence indicates that DUI relies on the same molecular machine of mitochondrial maternal inheritance that is common in most animals, our data are relevant not only to DUI but also to shed light on how differential segregations of mtDNA variants, in the same nuclear background, may be controlled during development.

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.

A Distinct Mitochondrial Genome with DUI-Like Inheritance in the Ocean Quahog Arctica islandica

Mitochondrial DNA (mtDNA) is strictly maternally inherited in metazoans. The major exception to this rule has been found in many bivalve species which allow the presence of different sex-linked mtDNA molecules. This mechanism, named doubly uniparental inheritance (DUI), is characterized by the presence of two mtDNAs: The female mtDNA is found in somatic tissue and female gonads, whereas the male mtDNA is usually found in male gonads and sperm. In this study we highlight the existence of two divergent mitochondrial haplotypes with a low genetic difference around 6–8% in Arctica islandica, a long-lived clam belonging to the Arcticidae, a sister group to the Veneridae in which DUI has been found. Phylogenetic analysis on cytochrome b and 16S sequences from somatic and gonadic tissues of clams belonging to different populations reveals the presence of the “divergent” type in male gonads only and the “normal” type in somatic tissues and female gonads. This peculiar segregation of divergent mtDNA types speaks for the occurrence of the DUI mechanism in A. islandica. This example also highlights the difficulties to assess the presence of such particular mitochondrial inheritance system and underlines the possible misinterpretations in phylogeographic and phylogenetic studies of bivalve species linked to the presence of two poorly differentiated mitochondrial genomes.

Extensive mitochondrial heteroplasmy in hybrid water frog (Pelophylax spp.) populations from Southeast Europe

Water frogs of the genus Pelophylax (previous Rana) species have been much studied in Europe for their outstanding reproductive mechanism in which sympatric hybridization between genetically distinct parental species produces diverse genetic forms of viable hybrid animals. The most common hybrid is P. esculentus that carries the genomes of both parental species, P. ridibundus and P. lessonae, but usually transfers the whole genome of only one parent to its offsprings (hybridogenesis). The evolutionary cost of transfer of the intact genome and hence the hemiclonal reproduction is the depletion of heterozygosity in the hybrid populations. Pelophylax esculentus presents an excellent example of the long‐term sustained hybridization and hemiclonal reproduction in which the effects of the low genetic diversity are balanced through the novel mutations and periodic recombinations. In this study, we analyzed the mitochondrial (mt) and microsatellites DNA variations in hybrid Pelophylax populations from southern parts of the Pannonian Basin and a north–south transect of the Balkan Peninsula, which are home for a variety of Pelophylax genetic lineages. The mtDNA haplotypes found in this study corresponded to P. ridibundus and P. epeiroticus of the Balkan – Anatolian lineage (ridibundus–bedriagae) and to P. lessonae and a divergent lessonae haplotype of the lessonae lineage. The mtDNA genomes showed considerable intraspecific variation and geographic differentiation. The Balkan wide distributed P. ridibundus was found in all studied populations and its nuclear genome, along with either the lessonae or the endemic epeiroticus genome, in all hybrids. An unexpected finding was that the hybrid populations were invariably heteroplasmic, that is, they contained the mtDNA of both parental species. We discussed the possibility that such extensive heteroplasmy is a result of hybridization and it comes from regular leakage of the paternal mtDNA from a sperm of one species that fertilizes eggs of another. In this case, the mechanisms that protect the egg from heterospecific fertilization and further from the presence of sperm mtDNA could become compromised due to their differences and divergence at both, mitochondrial and nuclear DNA. The heteroplasmy once retained in the fertilized egg could be transmitted by hybrid backcrossing to the progeny and maintained in a population over generations. The role of interspecies and heteroplasmic hybrid animals due to their genomic diversity and better fitness compare to the parental species might be of the special importance in adaptations to miscellaneous and isolated environments at the Balkan Peninsula.

Complete mitochondrial genome of the soft-shell clam Mya arenaria

We have sequenced and characterized the complete mitochondrial genome of the soft-shell clam, Mya arenaria, an important organism for environmental toxicology and aquaculture. Mya arenaria is located in the taxonomic order Myoida, which lacks any member with a completely annotated mitogenome. The M. arenaria mitochondrial genome is 17 947 bp in length. Like most marine bivalves, the circular mitogenome codes entirely on the heavy strand, with no introns. As with other bivalves, the gene order of the mitochondrion is highly rearranged. The mitogenome contains 12 protein-coding genes but ATP8 is missing, consistent with about half of all bivalve genera. Twenty-three tRNAs were identified. Phylogenetic analysis shows that M. arenaria is related most closely with the bivalves Sinonovacula constricta, and Moerella iridescens, of the infraclass Euheterodonta (unassigned). This, along with the close grouping of the phylogenetic trees, confirms a close tie between Myoida and Euheterodonta (unassigned).

Presence of two mitochondrial genomes in the mytilid Perumytilus purpuratus: Phylogenetic evidence for doubly uniparental inheritance

This study presents evidence, using sequences of ribosomal 16S and COI mtDNA, for the presence of two mitochondrial genomes in Perumytilus purpuratus. This may be considered evidence of doubly uniparental mtDNA inheritance. The presence of the two types of mitochondrial genomes differentiates females from males. The F genome was found in the somatic and gonadal tissues of females and in the somatic tissues of males; the M genome was found in the gonads and mantle of males only. For the mitochondrial 16S region, ten haplotypes were found for the F genome (nucleotide diversity 0.004), and 7 haplotypes for the M genome (nucleotide diversity 0.001), with a distance Dxy of 0.125 and divergence Kxy of 60.33%. For the COI gene 17 haplotypes were found for the F genome (nucleotide diversity 0.009), and 10 haplotypes for the M genome (nucleotide diversity 0.010), with a genetic distance Dxy of 0.184 and divergence Kxy of 99.97%. Our results report the presence of two well-differentiated, sex-specific types of mitochondrial genome (one present in the male gonad, the other in the female gonad), implying the presence of DUI in P. purpuratus. These results indicate that care must be taken in phylogenetic comparisons using mtDNA sequences of P. purpuratus without considering the sex of the individuals.

Searching for doubly uniparental inheritance of mtDNA in the apple snail Pomacea diffusa

Doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) is an exceptional mode of mtDNA transmission, restricted so far to the class of bivalves. We searched for DUI outside bivalves using the apple snail Pomacea diffusa. It was an appropriate candidate to search for DUI for three reasons; it belongs to gastropods, which is the closest sister group to bivalves, it is gonochoristic and it has a strong sex bias in the progeny of different female individuals. These phenomena (gonochorism and sex-biased progeny) are also found in species with DUI. We searched for heteroplasmy in males and for high sequence divergence among mtDNA sequences obtained from male and female gonads. All sequences examined were identical. These data suggest that the mtDNA in P. diffusa is maternally transmitted and DUI remains an exclusive characteristic of bivalves.

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.

Complete male mitochondrial genome of Anodonta anatina (Mollusca: Unionidae)

Anodonta anatina is a freshwater mussel of the family Unionidae. These mussels have a unique mitochondria inheritance system named doubly uniparental inheritance (DUI). Under DUI males have two, potentially very divergent mitochondrial genomes: F-type inherited from mother and M-type inherited from father. F-type is present in soma whereas M-type is present in gonadal tissues and sperm. Here we report two M-type sequences of complete mitochondrial genomes from Anodonta anatina. They are 16,906 bp long and their sequences are similar (0.1% divergence). The genome organization is identical to the other Unionidean M-type genomes published to date. There are 38 genes, including the recently described M-type specific M ORF. The presence of tRNA-like repeat in one of the noncoding regions, suggests that the control region is located in this area. Nucleotide composition is quite extreme, with AT content (66.2%) higher than in any other of the six published Unionidean M genomes.

Patterns of DNA barcode variation in Canadian marine molluscs

BACKGROUND

Molluscs are the most diverse marine phylum and this high diversity has resulted in considerable taxonomic problems. Because the number of species in Canadian oceans remains uncertain, there is a need to incorporate molecular methods into species identifications. A 648 base pair segment of the cytochrome c oxidase subunit I gene has proven useful for the identification and discovery of species in many animal lineages. While the utility of DNA barcoding in molluscs has been demonstrated in other studies, this is the first effort to construct a DNA barcode registry for marine molluscs across such a large geographic area.

METHODOLOGY/PRINCIPAL FINDINGS

This study examines patterns of DNA barcode variation in 227 species of Canadian marine molluscs. Intraspecific sequence divergences ranged from 0-26.4% and a barcode gap existed for most taxa. Eleven cases of relatively deep (>2%) intraspecific divergence were detected, suggesting the possible presence of overlooked species. Structural variation was detected in COI with indels found in 37 species, mostly bivalves. Some indels were present in divergent lineages, primarily in the region of the first external loop, suggesting certain areas are hotspots for change. Lastly, mean GC content varied substantially among orders (24.5%-46.5%), and showed a significant positive correlation with nearest neighbour distances.

CONCLUSIONS/SIGNIFICANCE

DNA barcoding is an effective tool for the identification of Canadian marine molluscs and for revealing possible cases of overlooked species. Some species with deep intraspecific divergence showed a biogeographic partition between lineages on the Atlantic, Arctic and Pacific coasts, suggesting the role of Pleistocene glaciations in the subdivision of their populations. Indels were prevalent in the barcode region of the COI gene in bivalves and gastropods. This study highlights the efficacy of DNA barcoding for providing insights into sequence variation across a broad taxonomic group on a large geographic scale.

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.

Mitochondrial DNA transcription levels during spermatogenesis and early development in doubly uniparental inheritance of the mitochondrial DNA system of the blue mussel Mytilus galloprovincialis

In some species of bivalve, there are two highly diverged mitochondrial genomes, one found in all individuals (F type) and the other normally in males only (M type). In Mytilus, a maternally-dependent sex ratio of the progeny has been reported. Some females almost exclusively produce daughters, while others produce a high proportion of sons. We previously reported that in M. galloprovincialis, M type mtDNA copy number may be maintained during spermatogenesis and the development of larvae of male-biased mothers to sustain the doubly uniparental inheritance system. In this study, we investigated transcription levels of M type mtDNA before and after fertilization to understand its function in the germ line. First, we quantified transcription levels of M type mtDNA in testicular cells dissected using laser-capture micro-dissection. The transcription levels of M type mtDNA were not significantly different between spermatogonia and spermatocytes versus spermatids and spermatozoa. Next, we examined differences in transcription levels of M type mtDNA between larvae from male-biased and female-biased mothers. The transcription levels of M type mtDNA significantly increased 24 and 48 h after fertilization in male-biased crosses. By contrast, transcription levels significantly decreased in female-biased crosses. These results suggest M type mtDNA may play a role in early germ line formation.

Proteomic analysis of eggs from Mytilus edulis females differing in mitochondrial DNA transmission mode

Many bivalves have an unusual mechanism of mitochondrial DNA (mtDNA) inheritance called doubly uniparental inheritance (DUI) in which distinctly different genomes are inherited through the female (F genome) and male (M genome) lineages. In fertilized eggs that will develop into male embryos, the sperm mitochondria remain in an aggregation, which is believed to be delivered to the primordial germ cells and passed to the next generation through the sperm. In fertilized eggs that will develop into female embryos, the sperm mitochondria are dispersed throughout the developing embryo and make little if any contribution to the next generation. The frequency of embryos with the aggregated or dispersed mitochondrial type varies among females. Previous models of DUI have predicted that maternal nuclear factors cause molecular differences among unfertilized eggs from females producing embryos with predominantly dispersed or aggregated mitochondria. We test this hypothesis using females of each of the two types from a natural population. We have found small, yet detectable, differences of the predicted type at the proteome level. We also provide evidence that eggs of females giving the dispersed pattern have consistently lower expression for different proteasome subunits than eggs of females giving the aggregated pattern. These results, combined with those of an earlier study in which we used hatchery lines of Mytilus, and with a transcriptomic study in a clam that has the DUI system of mtDNA transmission, reinforce the hypothesis that the ubiquitin-proteasome system plays a key role in the mechanism of DUI and sex determination in bivalves. We also report that eggs of females giving the dispersed pattern have higher expression for arginine kinase and enolase, enzymes involved in energy production, whereas ferritin, which is involved in iron homeostasis, has lower expression. We discuss these results in the context of genetic models for DUI and suggest experimental methods for further understanding the role of these proteins in DUI.

Diversity of mitochondrial genome organization

In this review, we discuss types of mitochondrial genome structural organization (architecture), which includes the following characteristic features: size and the shape of DNA molecule, number of encoded genes, presence of cryptogenes, and editing of primary transcripts.

Structure, transcription, and variability of metazoan mitochondrial genome: perspectives from an unusual mitochondrial inheritance system

Despite its functional conservation, the mitochondrial genome (mtDNA) presents strikingly different features among eukaryotes, such as size, rearrangements, and amount of intergenic regions. Nonadaptive processes such as random genetic drift and mutation rate play a fundamental role in shaping mtDNA: the mitochondrial bottleneck and the number of germ line replications are critical factors, and different patterns of germ line differentiation could be responsible for the mtDNA diversity observed in eukaryotes. Among metazoan, bivalve mollusc mtDNAs show unusual features, like hypervariable gene arrangements, high mutation rates, large amount of intergenic regions, and, in some species, an unique inheritance system, the doubly uniparental inheritance (DUI). The DUI system offers the possibility to study the evolutionary dynamics of mtDNAs that, despite being in the same organism, experience different genetic drift and selective pressures. We used the DUI species Ruditapes philippinarum to study intergenic mtDNA functions, mitochondrial transcription, and polymorphism in gonads. We observed: 1) the presence of conserved functional elements and novel open reading frames (ORFs) that could explain the evolutionary persistence of intergenic regions and may be involved in DUI-specific features; 2) that mtDNA transcription is lineage-specific and independent from the nuclear background; and 3) that male-transmitted and female-transmitted mtDNAs have a similar amount of polymorphism but of different kinds, due to different population size and selection efficiency. Our results are consistent with the hypotheses that mtDNA evolution is strongly dependent on the dynamics of germ line formation, and that the establishment of a male-transmitted mtDNA lineage can increase male fitness through selection on sperm function.

Mitochondrial DNA paradox: sex-specific genetic structure in a marine mussel--despite maternal inheritance and passive dispersal

Background

When genetic structure is identified using mitochondrial DNA (mtDNA), but no structure is identified using biparentally-inherited nuclear DNA, the discordance is often attributed to differences in dispersal potential between the sexes.

Results

We sampled the intertidal rocky shore mussel Perna perna in a South African bay and along the nearby open coast, and sequenced maternally-inherited mtDNA (there is no evidence for paternally-inherited mtDNA in this species) and a biparentally-inherited marker. By treating males and females as different populations, we identified significant genetic structure on the basis of mtDNA data in the females only.

Conclusions

This is the first study to report sex-specific differences in genetic structure based on matrilineally-inherited mtDNA in a passively dispersing species that lacks social structure or sexual dimorphism. The observed pattern most likely stems from females being more vulnerable to selection in habitats from which they did not originate, which also manifests itself in a male-biased sex ratio. Our results have three important implications for the interpretation of population genetic data. First, even when mtDNA is inherited exclusively in the female line, it also contains information about males. For that reason, using it to identify sex-specific differences in genetic structure by contrasting it with biparentally-inherited markers is problematic. Second, the fact that sex-specific differences were found in a passively dispersing species in which sex-biased dispersal is unlikely highlights the fact that significant genetic structure is not necessarily a function of low dispersal potential or physical barriers. Third, even though mtDNA is typically used to study historical demographic processes, it also contains information about contemporary processes. Higher survival rates of males in non-native habitats can erase the genetic structure present in their mothers within a single generation.