Novel protein genes in animal mtDNA: a new sex determination system in freshwater mussels (Bivalvia: Unionoida)?

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

Integrative phylogenetic, phylogeographic and morphological characterisation of the Unio crassus species complex reveals cryptic diversity with important conservation implications

The global decline of freshwater mussels and their crucial ecological services highlight the need to understand their phylogeny, phylogeography and patterns of genetic diversity to guide conservation efforts. Such knowledge is urgently needed for Unio crassus, a highly imperilled species originally widespread throughout Europe and southwest Asia. Recent studies have resurrected several species from synonymy based on mitochondrial data, revealing U. crassus to be a complex of cryptic species. To address long-standing taxonomic uncertainties hindering effective conservation, we integrate morphometric, phylogenetic, and phylogeographic analyses to examine species diversity within the U. crassus complex across its entire range. Phylogenetic analyses were performed using cytochrome c oxidase subunit I (815 specimens from 182 populations) and, for selected specimens, whole mitogenome sequences and Anchored Hybrid Enrichment (AHE) data on ∼ 600 nuclear loci. Mito-nuclear discordance was detected, consistent with mitochondrial DNA gene flow between some species during the Pliocene and Pleistocene. Fossil-calibrated phylogenies based on AHE data support a Mediterranean origin for the U. crassus complex in the Early Miocene. The results of our integrative approach support 12 species in the group: the previously recognised Unio bruguierianus, Unio carneus, Unio crassus, Unio damascensis, Unio ionicus, Unio sesirmensis, and Unio tumidiformis, and the reinstatement of five nominal taxa: Unio desectusstat. rev., Unio gontieriistat. rev., Unio mardinensisstat. rev., Unio nanusstat. rev., and Unio vicariusstat. rev. Morphometric analyses of shell contours reveal important morphospace overlaps among these species, highlighting cryptic, but geographically structured, diversity. The distribution, taxonomy, phylogeography, and conservation of each species are succinctly described.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Mitonuclear Sex Determination? Empirical Evidence from Bivalves

Genetic elements encoded in nuclear DNA determine the sex of an individual in many animals. In certain bivalve lineages that possess doubly uniparental inheritance (DUI), mitochondrial DNA (mtDNA) has been hypothesized to contribute to sex determination. In these cases, females transmit a female mtDNA to all offspring, while male mtDNA (M mtDNA) is transmitted only from fathers to sons. Because M mtDNA is inherited in the same way as Y chromosomes, it has been hypothesized that mtDNA may be responsible for sex determination. However, the role of mitochondrial and nuclear genes in sex determination has yet to be validated in DUI bivalves. In this study, we used DNA, RNA, and mitochondrial short noncoding RNA (sncRNA) sequencing to explore the role of mitochondrial and nuclear elements in the sexual development pathway of the freshwater mussel Potamilus streckersoni (Bivalvia: Unionida). We found that the M mtDNA sheds a sncRNA partially within a male-specific mitochondrial gene that targets a pathway hypothesized to be involved in female development and mitophagy. RNA-seq confirmed the gene target was significantly upregulated in females, supporting a direct role of mitochondrial sncRNAs in gene silencing. These findings support the hypothesis that M mtDNA inhibits female development. Genome-wide patterns of genetic differentiation and heterozygosity did not support a nuclear sex-determining region, although we cannot reject that nuclear factors are involved with sex determination. Our results provide further evidence that mitochondrial loci contribute to diverse, nonrespiratory functions and additional insights into an unorthodox sex-determining system.

Validation of the male-specific ORF of the paternally-transmitted mtDNA in Mytilus edulis as a protein-coding gene

There appears to be an additional set of sex-specific mtDNA-encoded proteins in bivalve species with doubly uniparental mitochondrial inheritance that may be involved in the transmission of the female and male mitogenomes. In the marine mussel Mytilus edulis, the translation of the female-specific open reading frame (F-ORF) was demonstrated but the translation of the male-specific ORF (M-ORF) remains to be shown. Here we validate the male-specific ORF of the paternal mitogenome in M. edulis as a protein-coding gene. The M-ORF protein was detected only in male gonads and localized in sperm mitochondria and acrosome, suggesting that it is involved in a key sperm function in Mytilus edulis.

Mitonuclear sex determination? Empirical evidence from bivalves

Genetic elements encoded in nuclear DNA determine the sex of an individual in many animals. In bivalves, however, mitochondrial DNA (mtDNA) has been hypothesized to contribute to sex determination in lineages that possess doubly uniparental inheritance (DUI). In these cases, females transmit a female mtDNA (F mtDNA) to all offspring, while male mtDNA (M mtDNA) is transmitted only from fathers to sons. Because M mtDNA is inherited in the same way as Y chromosomes, it has been hypothesized that mtDNA may be responsible for sex determination. However, the role of mitochondrial and nuclear genes in sex determination has yet to be validated in DUI bivalves. In this study, we used DNA, RNA, and mitochondrial short non-coding RNA (sncRNA) sequencing to explore the role of mitochondrial and nuclear elements in the sexual development pathway of the freshwater mussel Potamilus streckersoni (Bivalvia: Unionida). We found that the M mtDNA shed a sncRNA partially within a male-specific mitochondrial gene that targeted pathways hypothesized to be involved in female development and mitophagy. RNA-seq confirmed the gene target was significantly upregulated in females, supporting a direct role of mitochondrial sncRNAs in gene silencing. These findings support the hypothesis that M mtDNA inhibits female development. Genome-wide patterns of genetic differentiation and heterozygosity did not support a nuclear sex determining region, although we cannot reject that nuclear factors are involved with sex determination. Our results provide further evidence that mitochondrial loci contribute to diverse, non-respiratory functions and provide a first glimpse into an unorthodox sex determining system.

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

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.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BACKGROUND

Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.

RESULTS

We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.

CONCLUSIONS

The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

The male-type mitochondrial genome of the freshwater mussel Potamilus streckersoni Smith, Johnson, Inoue, Doyle, & Randklev, 2019 (Bivalvia: Unionidae)

The global decline of freshwater mussels emphasizes the need to establish genetic resources to better understand their biology, including a unique mitochondrial biology known as doubly uniparental inheritance. In this study, we present the complete male-type (M-type) mitochondrial genome of the freshwater mussel, Potamilus streckersoni Smith, Johnson, Inoue, Doyle, & Randklev, 2019. The M-type mtDNA is approximately 16 kilobases and contains 22 tRNAs, two rRNAs, and 14 protein-coding genes, including a male-specific open reading frame. Read coverage revealed that M-type mtDNA was more abundant than female-type mtDNA in male gonadal tissue, with respect to a non-spawning male individual. Novel mitogenomes were resolved within previously described sex-specific monophyletic clades across the subfamily Ambleminae. The availability of high-quality nuclear and mitochondrial genomic data for P. streckersoni makes it a model for future research into the potential role of mtDNA in sex determination or sexual development in freshwater mussels.

The longest mitochondrial protein in metazoans is encoded by the male-transmitted mitogenome of the bivalve Scrobicularia plana

Cytochrome c oxidase subunit II (COX2) is one of the three mitochondrially encoded proteins of the complex IV of the respiratory chain that catalyses the reduction of oxygen to water. The cox2 gene spans about 690 base pairs in most animal species and produces a protein composed of approximately 230 amino acids. We discovered an extreme departure from this pattern in the male-transmitted mitogenome of the bivalve Scrobicularia plana with doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA), which possesses an important in-frame insertion of approximately 4.8 kb in its cox2 gene. This feature—an enlarged male cox2 gene—is found in many species with DUI; the COX2 protein can be up to 420 amino acids long. Through RT-PCRs, immunoassays and comparative genetics, the evolution and functionality of this insertion in S. plana were characterized. The in-frame insertion is conserved among individuals from different populations and bears the signature of purifying selection seemingly indicating maintenance of functionality. Its transcription and translation were confirmed: this gene produces a polypeptide of 1892 amino acids, making it the largest metazoan COX2 protein known to date. We hypothesize that these extreme modifications in the COX2 protein affect the metabolism of mitochondria containing the male-transmitted mtDNA in Scrobicularia plana.

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.

Did doubly uniparental inheritance (DUI) of mtDNA originate as a cytoplasmic male sterility (CMS) system?

Animal and plant species exhibit an astonishing diversity of sexual systems, including environmental and genetic determinants of sex, with the latter including genetic material in the mitochondrial genome. In several hermaphroditic plants for example, sex is determined by an interaction between mitochondrial cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Specifically, CMS involves aberrant mitochondrial genes that prevent pollen development and specific nuclear genes that restore it, leading to a mixture of female (male-sterile) and hermaphroditic individuals in the population (gynodioecy). Such a mitochondrial-nuclear sex determination system is thought to be rare outside plants. Here, we present one possible case of CMS in animals. We hypothesize that the only exception to the strict maternal mtDNA inheritance in animals, the doubly uniparental inheritance (DUI) system in bivalves, might have originated as a mitochondrial-nuclear sex-determination system. We document and explore similarities that exist between DUI and CMS, and we propose various ways to test our hypothesis.

Mussel Unio douglasiae MG from the Chihe River: mitogenome description and phylogenetic analysis

The genus Unio is one of the widespread freshwater bivalves. To date, its intra-generic phylogeny remains controversial and therefore the mitochondrial genome data is needed. Here, we report the complete mitogenome of Unio douglasiae MG that is distributed in the Chihe River, a branch of Huaihe River, East China. This mitochondrial genome is 15,764 base pair in total length. It consists of 37 genes: 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes (12S and 16S). The base composition was 38.38% for A, 26.48% for T, 23.17% for C, and 11.98% for G, showing an obvious bias of higher A + T content (64.86%) than the G + C content (35.14%). Phylogenetic analysis showed that U. douglasiae MG is clustered with other Unio and Nodularia mollusks in the family Unionidae. These results showed that combine with morphological techniques, the mitogenome can provide useful information to further understanding of the genetics, systematics, and conservation of this endangered species.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Presence of male mitochondria in somatic tissues and their functional importance at the whole animal level in the marine bivalve Arctica islandica

Metazoans normally possess a single lineage of mitochondria inherited from the mother (♀-type mitochondria) while paternal mitochondria are absent or eliminated in fertilized eggs. In doubly uniparental inheritance (DUI), which is specific to the bivalve clade including the ocean quahog, Arctica islandica, ♂-type mitochondria are retained in male gonads and, in a few species, small proportions of ♂-type mitochondria co-exist with ♀-type in somatic tissues. To the best of our knowledge, we report, for the first time in metazoan, the natural occurrence of male and female individuals with exclusively ♂-type mitochondria in somatic tissues of the bivalve A. islandica. Mitochondrial genomes differ by ~5.5% at DNA sequence level. Exclusive presence of ♂-type mitochondria affects mitochondrial complexes partially encoded by mitochondrial genes and leads to a sharp drop in respiratory capacity. Through a combination of whole mitochondrial genome sequencing and molecular assays (gene presence and expression), we demonstrate that 1) 11% of individuals of an Icelandic population appear homoplasmic for ♂-type mitochondria in somatic tissues, 2) ♂-type mitochondrial genes are transcribed and 3) individuals with ♂-type mitochondria in somatic cells lose 30% of their wild-type respiratory capacity. This mitochondrial pattern in A. islandica is a special case of DUI, highlighted in individuals from both sexes with functional consequences at cellular and conceivably whole animal level.

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.

Bioenergetic consequences of sex-specific mitochondrial DNA evolution

Doubly uniparental inheritance (DUI) represents a notable exception to the general rule of strict maternal inheritance (SMI) of mitochondria in metazoans. This system entails the coexistence of two mitochondrial lineages (F- and M-type) transmitted separately through oocytes and sperm, thence providing an unprecedented opportunity for the mitochondrial genome to evolve adaptively for male functions. In this study, we explored the impact of a sex-specific mitochondrial evolution upon gamete bioenergetics of DUI and SMI bivalve species, comparing the activity of key enzymes of glycolysis, fermentation, fatty acid metabolism, tricarboxylic acid cycle, oxidative phosphorylation and antioxidant metabolism. Our findings suggest reorganized bioenergetic pathways in DUI gametes compared to SMI gametes. This generally results in a decreased enzymatic capacity in DUI sperm with respect to DUI oocytes, a limitation especially prominent at the terminus of the electron transport system. This bioenergetic remodelling fits a reproductive strategy that does not require high energy input and could potentially link with the preservation of the paternally transmitted mitochondrial genome in DUI species. Whether this phenotype may derive from positive or relaxed selection acting on DUI sperm is still uncertain.

The complete mitogenome of the inarticulate brachiopod Glottidia pyramidata reveals insights into gene order variation, deviant ATP8 and mtORFans in the Brachiopoda

Brachiopods are a clade of marine organisms with a tremendously diverse and abundant fossil record but with fewer than 500 species extant today. Even if a better understanding of their biology and genetics could help to test hypotheses about their impressive decline, knowledge of genetics and evolutionary genomics in extant brachiopods is very poor. Here, we present the complete mitochondrial genome sequence of the inarticulate Glottidia pyramidata, an eastern North American extant representative of the phylum Brachiopoda. Besides the general characteristics of the sequenced mitogenome, we present its unusual features such as deviant ATP8 protein sequence and supernumerary ORFs, and also unique gene order, considering the available genome sequences of other brachiopod species.

A proposed method for analyzing molecular signatures to detect hermaphroditism in freshwater mussels: a case study using the eastern floater (Pyganodon cataracta)

Many freshwater mussels (order Unionida) have an unusual system of doubly uniparental inheritance (DUI) of mitochondrial (mt) DNA. In species with DUI, males possess a female-transmitted (F-type) mt genome and a male-transmitted (M-type) mt genome. These genomes contain non-canonical open reading frame (orf) genes, referred to as f-orf and m-orf, present in F and M mt genomes, respectively. These genes have been implicated in sexual development in Unionida. When gonochoric species become hermaphroditic, which has happened several times in Unionida, they lose their M-type mt genome and f-orf genes evolve dramatically. Resulting F-ORF proteins are highly divergent in terms of primary nucleotide sequence, inferred amino acids, and hydrophobic properties; these genes (and proteins) are referred to as hermaphroditic orfs or h-orfs (and H-ORFs). We investigated patterns of hydrophobicity divergence for H-ORF proteins in hermaphrodites versus F-ORF proteins in closely related gonochoric species against cytochrome c oxidase subunit 1 (cox1) divergences. This approach was used to assess whether cryptic hermaphrodites can be detected. Although we did not detect evidence for the recent transition of any populations of eastern floaters (Pyganodon cataracta (Say, 1817)) to hermaphroditism, our analyses demonstrate that molecular signatures in mtDNA can be used to detect hermaphroditism in freshwater mussels.

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.

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

Molluscan mitochondrial genomes break the rules

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

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

Semimytilus algosus: first known hermaphroditic mussel with doubly uniparental inheritance of mitochondrial DNA

Doubly uniparental inheritance (DUI) of mitochondrial DNA is a rare phenomenon occurring in some freshwater and marine bivalves and is usually characterized by the mitochondrial heteroplasmy of male individuals. Previous research on freshwater Unionida mussels showed that hermaphroditic species do not have DUI even if their closest gonochoristic counterparts do. No records showing DUI in a hermaphrodite have ever been reported. Here we show for the first time that the hermaphroditic mussel Semimytilus algosus (Mytilida), very likely has DUI, based on the complete sequences of both mitochondrial DNAs and the distribution of mtDNA types between male and female gonads. The two mitogenomes show considerable divergence (34.7%). The presumably paternal M type mitogenome dominated the male gonads of most studied mussels, while remaining at very low or undetectable levels in the female gonads of the same individuals. If indeed DUI can function in the context of simultaneous hermaphroditism, a change of paradigm regarding its involvement in sex determination is needed. It is apparently associated with gonadal differentiation rather than with sex determination in bivalves.

Metabolic remodelling associated with mtDNA: insights into the adaptive value of doubly uniparental inheritance of mitochondria

Mitochondria produce energy through oxidative phosphorylation (OXPHOS), which depends on the expression of both nuclear and mitochondrial DNA (mtDNA). In metazoans, a striking exception from strictly maternal inheritance of mitochondria is doubly uniparental inheritance (DUI). This unique system involves the maintenance of two highly divergent mtDNAs (F- and M-type, 8-40% of nucleotide divergence) associated with gametes, and occasionally coexisting in somatic tissues. To address whether metabolic differences underlie this condition, we characterized the OXPHOS activity of oocytes, spermatozoa, and gills of different species through respirometry. DUI species express different gender-linked mitochondrial phenotypes in gametes and partly in somatic tissues. The M-phenotype is specific to sperm and entails (i) low coupled/uncoupled respiration rates, (ii) a limitation by the phosphorylation system, and (iii) a null excess capacity of the final oxidases, supporting a strong control over the upstream complexes. To our knowledge, this is the first example of a phenotype resulting from direct selection on sperm mitochondria. This metabolic remodelling suggests an adaptive value of mtDNA variations and we propose that bearing sex-linked mitochondria could assure the energetic requirements of different gametes, potentially linking male-energetic adaptation, mitotype preservation and inheritance, as well as resistance to both heteroplasmy and ageing.

Natural Heteroplasmy and Mitochondrial Inheritance in Bivalve Molluscs

Heteroplasmy is the presence of more than one type of mitochondrial genome within an individual, a condition commonly reported as unfavorable and affecting mitonuclear interactions. So far, no study has investigated heteroplasmy at protein level, and whether it occurs within tissues, cells, or even organelles. The only known evolutionarily stable and natural heteroplasmic system in Metazoa is the Doubly Uniparental Inheritance (DUI)-reported so far in ∼100 bivalve species-in which two mitochondrial lineages are present: one transmitted through eggs (F-type) and the other through sperm (M-type). Because of such segregation, mitochondrial oxidative phosphorylation proteins reach a high amino acid sequence divergence (up to 52%) between the two lineages in the same species. Natural heteroplasmy coupled with high sequence divergence between F- and M-type proteins provides a unique opportunity to study their expression and assess the level and extent of heteroplasmy. Here, for the first time, we immunolocalized F- and M-type variants of three mitochondrially-encoded proteins in the DUI species Ruditapes philippinarum, in germline and somatic tissues at different developmental stages. We found heteroplasmy at organelle level in undifferentiated germ cells of both sexes, and in male soma, whereas gametes were homoplasmic: eggs for the F-type and sperm for the M-type. Thus, during gametogenesis, only the sex-specific mitochondrial variant is maintained, likely due to a process of meiotic drive. We examine the implications of our results for DUI proposing a revised model, and we discuss interactions of mitochondria with germ plasm and their role in germline development. Molecular and phylogenetic evidence suggests that DUI evolved from the common Strictly Maternal Inheritance, so the two systems likely share the same underlying molecular mechanism, making DUI a useful system for studying mitochondrial biology.

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.

Mesozoic mitogenome rearrangements and freshwater mussel (Bivalvia: Unionoidea) macroevolution

Using a new fossil-calibrated mitogenome-based approach, we identified macroevolutionary shifts in mitochondrial gene order among the freshwater mussels (Unionoidea). We show that the early Mesozoic divergence of the two Unionoidea clades, Margaritiferidae and Unionidae, was accompanied by a synchronous split in the gene arrangement in the female mitogenome (i.e., gene orders MF1 and UF1). Our results suggest that this macroevolutionary jump was completed within a relatively short time interval (95% HPD 201-226 Ma) that coincided with the Triassic-Jurassic mass extinction. Both gene orders have persisted within these clades for ~200 Ma. The monophyly of the so-called "problematic" Gonideinae taxa was supported by all the inferred phylogenies in this study using, for the first time, the M- and F-type mitogenomes either singly or combined. Within Gonideinae, two additional splits in the gene order (UF1 to UF2, UF2 to UF3) occurred in the Mesozoic and have persisted for ~150 and ~100 Ma, respectively. Finally, the mitogenomic results suggest ancient connections between freshwater basins of East Asia and Europe near the Cretaceous-Paleogene boundary, probably via a continuous paleo-river system or along the Tethys coastal line, which are well supported by at least three independent but almost synchronous divergence events.

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.

Phylogeny of European Anodontini (Bivalvia: Unionidae) with a redescription of Anodonta exulcerata

Freshwater bivalves are highly threatened and globally declining due to multiple anthropogenic impacts, making them important conservation targets. Because conservation policies and actions generally occur at the species level, accurate species identification and delimitation is critical. A recent phylogenetic study of Italian mussel populations revalidated an Anodonta species bringing the number of known European Anodontini from three to four species. The current study contributes to the clarification of the taxonomy and systematics of European Anodontini, using a combination of molecular, morphological and anatomical data, and constructs phylogenies based on complete mitogenomes. A redescription of A. exulcerata and a comparative analysis of morphological and anatomical characters with respect to the other two species of Anodonta present in the area are provided. No reliable diagnostic character has emerged from comparative analysis of the morphometric characters of 109 specimens from 16 sites across the Italian peninsula. In fact, the discriminant analysis resulted in a greater probability of correct assignment to the site of origin than to the species. This confirms the difficulties of an uncritical application of visual characters for the delimitation of species, especially for Anodontinae.

Putative Mitochondrial Sex Determination in the Bivalvia: Insights From a Hybrid Transcriptome Assembly in Freshwater Mussels

Bivalves exhibit an astonishing diversity of sexual systems, with genetic and environmental determinants of sex, and possibly the only example of mitochondrial genes influencing sex determination pathways in animals. In contrast to all other animal species in which strict maternal inheritance (SMI) of mitochondria is the rule, bivalves possess a system known as doubly uniparental inheritance (DUI) of mitochondria in which maternal and paternal mitochondria (and their corresponding female-transmitted or F mtDNA and male-transmitted or M mtDNA genomes) are transmitted within a species. Species with DUI also possess sex-associated mtDNA-encoded proteins (in addition to the typical set of 13), which have been hypothesized to play a role in sex determination. In this study, we analyzed the sex-biased transcriptome in gonads of two closely-related freshwater mussel species with different reproductive and mitochondrial transmission modes: the gonochoric, DUI species, Utterbackia peninsularis, and the hermaphroditic, SMI species, Utterbackia imbecillis. Through comparative analysis with other DUI and non-DUI bivalve transcriptomes already available, we identify common male and female-specific genes, as well as SMI and DUI-related genes, that are probably involved in sex determination and mitochondrial inheritance in this animal group. Our results contribute to the understanding of what could be the first animal sex determination system involving the mitochondrial genome.

Differences Between the Raw Material and the Products of de Novo Gene Birth Can Result from Mutational Biases

Proteins are among the most important constituents of biological systems. Because all protein-coding genes have a noncoding ancestral form, the properties of noncoding sequences and how they shape the birth of novel proteins may influence the structure and function of all proteins. Differences between the properties of young proteins and random expectations from noncoding sequences have previously been interpreted as the result of natural selection. However, interpreting such deviations requires a yet-unattained understanding of the raw material of de novo gene birth and its relation to novel functional proteins. We mathematically show that the average properties and selective filtering of the "junk" polypeptides of which this raw material is composed are not the only factors influencing the properties of novel functional proteins. We find that in some biological scenarios, they also depend on the variance of the properties of junk polypeptides and their correlation with the rate of allelic turnover, which may itself depend on mutational biases. This suggests for instance that any property of polypeptides that accelerates their exploration of the sequence space could be overrepresented in novel functional proteins, even if it has a limited effect on adaptive value. To exemplify the use of our general theoretical results, we build a simple model that predicts the mean length and mean intrinsic disorder of novel functional proteins from the genomic GC content and a single evolutionary parameter. This work provides a theoretical framework that can guide the prediction and interpretation of results when studying the de novo emergence of protein-coding genes.

[The underestimated coding potential of mitochondrial DNA]

Mitochondria are ancient organelles that emerged from the endosymbiosis of free-living proto-bacteria. They still retain a semi-autonomous genetic system with a small genome. Mitochondrial DNA (mtDNA) codes for 13 essential proteins for the production of ATP, the sequences of which are relatively conserved across Metazoans. The discovery of additional mitochondria-derived peptides (MDPs) indicates an underestimated coding potential. Humanin, an anti-apoptotic peptide, is likely independently transcribed from within the 16S rRNA gene, as are recently described SHLPs. MOTS-c, discovered in silico, has been demonstrated to be involved in metabolism and insulin sensitivity. Gau, is a positionally conserved open reading frame (ORF) sequence found in the antisense strand of the COX1 gene and its corresponding peptide is strictly colocalized with mitochondrial markers. In bivalves with doubly uniparental inheritance of mtDNA, male and female mtDNAs each carry a separate additional gene possibly involved in sex determination. Other MDPs likely exist and their investigation will shed light on the underestimated functional repertoire of mitochondria.

Analysis of mitochondrial genomes resolves the phylogenetic position of Chinese freshwater mussels (Bivalvia, Unionidae)

The Yangtze River basin is one of the most species-rich regions for freshwater mussels on Earth, but is gravely threatened by anthropogenic activities. However, conservation planning and management of mussel species has been hindered by a number of taxonomic uncertainties. In order to clarify the taxonomic status and phylogenetic position of these species, mitochondrial genomes of four species (Acuticostachinensis, Schistodesmuslampreyanus, Cuneopsisheudei and Cuneopsiscapitatus) were generated and analyzed along with data from 43 other mitogenomes. The complete F-type mitogenomes of A.chinensis, S.lampreyanus, C.heudei, and C.capitatus are 15652 bp, 15855 bp, 15892 bp, and 15844 bp, respectively, and all four F-type mitogenomes have the same pattern of gene arrangement. ML and BI trees based on the mitogenome dataset are completely congruent, and indicate that the included Unionidae belong to three subfamilies with high bootstrap and posterior probabilities, i.e., Unioninae (Aculamprotula, Cuneopsis, Nodularia, and Schistodesmus), Anodontinae (Cristaria, Arconaia, Acuticosta, Lanceolaria, Anemina, and Sinoanodonta), and Gonideinae (Ptychorhynchus, Solenaia, Lamprotula, and Sinohyriopsis). Results also indicate that A.chinensis has affinities with Arconaialanceolata and Lanceolariagrayii and is a member of the subfamily Anodontinae.

Genome Survey of the Freshwater Mussel Venustaconcha ellipsiformis (Bivalvia: Unionida) Using a Hybrid De Novo Assembly Approach

Freshwater mussels (Bivalvia: Unionida) serve an important role as aquatic ecosystem engineers but are one of the most critically imperilled groups of animals. Here, we used a combination of sequencing strategies to assemble and annotate a draft genome of Venustaconcha ellipsiformis, which will serve as a valuable genomic resource given the ecological value and unique “doubly uniparental inheritance” mode of mitochondrial DNA transmission of freshwater mussels. The genome described here was obtained by combining high-coverage short reads (65× genome coverage of Illumina paired-end and 11× genome coverage of mate-pairs sequences) with low-coverage Pacific Biosciences long reads (0.3× genome coverage). Briefly, the final scaffold assembly accounted for a total size of 1.54 Gb (366,926 scaffolds, N50 = 6.5 kb, with 2.3% of “N” nucleotides), representing 86% of the predicted genome size of 1.80 Gb, while over one third of the genome (37.5%) consisted of repeated elements and >85% of the core eukaryotic genes were recovered. Given the repeated genetic bottlenecks of V. ellipsiformis populations as a result of glaciations events, heterozygosity was also found to be remarkably low (0.6%), in contrast to most other sequenced bivalve species. Finally, we reassembled the full mitochondrial genome and found six polymorphic sites with respect to the previously published reference. This resource opens the way to comparative genomics studies to identify genes related to the unique adaptations of freshwater mussels and their distinctive mitochondrial inheritance mechanism.

Mitogenomics of Perumytilus purpuratus (Bivalvia: Mytilidae) and its implications for doubly uniparental inheritance of mitochondria

Animal mitochondria are usually inherited through the maternal lineage. The exceptional system allowing fathers to transmit their mitochondria to the offspring exists in some bivalves. Its taxonomic spread is poorly understood and new mitogenomic data are needed to fill the gap. Here, we present for the first time the two divergent mitogenomes from Chilean mussel Perumytilus purpuratus. The existence of these sex-specific mitogenomes confirms that this species has the doubly uniparental inheritance (DUI) of mitochondria. The genetic distance between the two mitochondrial lineages in P. purpuratus is not only much bigger than in the Mytilus edulis species complex but also greater than the distance observed in Musculista senhousia, the only other DUI-positive member of the Mytilidae family for which both complete mitochondrial genomes were published to date. One additional, long ORF (open reading frame) is present exclusively in the maternal mitogenome of P. purpuratus. This ORF evolves under purifying selection, and will likely be a target for future DUI research.