Deciphering the Link between Doubly Uniparental Inheritance of mtDNA and Sex Determination in Bivalves: Clues from Comparative Transcriptomics

The transcriptome assembly of the European freshwater mussel Unio elongatulus C. Pfeiffer, 1825

Freshwater mussels of the order Unionida are a global conservation concern. Species of this group are strictly freshwater, sessile, slow-growing animals and, extremely sensitive to environmental changes. Human-mediated changes in freshwater habitats are imposing enormous pressure on the survival of freshwater mussels. Although a few flagship species are protected in Europe, other highly imperilled species receive much less attention. Moreover, knowledge about biology, ecology, and evolution and proper conservation assessments of many European species are still sparse. This knowledge gap is further aggravated by the lack of genomic resources available, which are key tools for conservation. Here we present the transcriptome assembly of Unio elongatulus C. Pfeiffer, 1825, one of the least studied European freshwater mussels. Using the individual sequencing outputs from eight physiologically representative mussel tissues, we provide an annotated panel of tissue-specific Relative Gene Expression profiles. These resources are pivotal to studying the species' biological and ecological features, as well as helping to understand its vulnerability to current and future threats.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

PacBio Hi-Fi genome assembly of the Iberian dolphin freshwater mussel Unio delphinus Spengler, 1793

Mussels of order Unionida are a group of strictly freshwater bivalves with nearly 1,000 described species widely dispersed across world freshwater ecosystems. They are highly threatened showing the highest record of extinction events within faunal taxa. Conservation is particularly concerning in species occurring in the Mediterranean biodiversity hotspot that are exposed to multiple anthropogenic threats, possibly acting in synergy. That is the case of the dolphin freshwater mussel Unio delphinus Spengler, 1793, endemic to the western Iberian Peninsula with recently strong population declines. To date, only four genome assemblies are available for the order Unionida and only one European species. We present the first genome assembly of Unio delphinus. We used the PacBio HiFi to generate a highly contiguous genome assembly. The assembly is 2.5 Gb long, possessing 1254 contigs with a contig N50 length of 10 Mbp. This is the most contiguous freshwater mussel genome assembly to date and is an essential resource for investigating the species' biology and evolutionary history that ultimately will help to support conservation strategies.

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.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

The gill transcriptome of threatened European freshwater mussels

Genomic tools applied to non-model organisms are critical to design successful conservation strategies of particularly threatened groups. Freshwater mussels of the Unionida order are among the most vulnerable taxa and yet almost no genetic resources are available. Here, we present the gill transcriptomes of five European freshwater mussels with high conservation concern: Margaritifera margaritifera, Unio crassus, Unio pictorum, Unio mancus and Unio delphinus. The final assemblies, with N50 values ranging from 1069–1895 bp and total BUSCO scores above 90% (Eukaryote and Metazoan databases), were structurally and functionally annotated, and made available. The transcriptomes here produced represent a valuable resource for future studies on these species’ biology and ultimately guide their conservation.

Measurement(s)	transcriptomics
Technology Type(s)	Illumina sequencing
Sample Characteristic - Organism	Margaritifera margaritifera • Unio crassus • Unio delphinus • Unio mancus • Unio pictorum
Sample Characteristic - Location	Europe

Comparative phylotranscriptomics reveals putative sex differentiating genes across eight diverse bivalve species

Mollusks, especially bivalves, exhibit a great diversity of sex determining mechanisms, including both genetic and environmental sex determination. Some bivalve species can be gonochoristic (separate sexes), while others are hermaphroditic (sequential or simultaneous). Several models have been proposed for specific bivalve species, utilizing information gained from gene expression data, as well as limited RAD-seq data (e.g., from Crassostrea gigas). However, these mechanisms are not as well studied as those in model organisms (e.g., Mus musculus, Drosophila melanogaster, Caenorhabditis elegans) and many genes involved in sex differentiation are not well characterized. We used phylotranscriptomics to better understand which possible sex differentiating genes are in bivalves and how these genes relate to similar genes in diverse phyla. We collected RNAseq data from eight phylogenetically diverse bivalve species: Argopecten irradians, Ensis directus, Geukensia demissa, Macoma tenta, Mercenaria mercenaria, Mya arenaria, Mytilus edulis, and Solemya velum. Using these data, we assembled representative transcriptomes for each species. We then searched for candidate sex differentiating genes using BLAST and confirmed the identity of nine genes using phylogenetics analyses from nine phyla. To increase the confidence of identification, we included ten bivalve genomes in our analyses. From the analysis of doublesex and mab-3 related transcription factor (DMRT) genes, we confirmed the identify of a Mollusk-specific sex determining DMRT gene: DMRT1L. Based on gene expression data from M. edulis and previous research, DMRT1L and FoxL2 are key genes for male and female development, respectively.

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.

Pathophysiological Responses of Pinna nobilis Individuals Enlightens the Etiology of Mass Mortality Situation in the Mediterranean Populations

Due to the rapid decrease of Pinna nobilis populations during the previous decades, this bivalve species, endemic in the Mediterranean Sea, is characterized as 'critically endangered'. In addition to human pressures, various pathogen infections have resulted in extended reduction, even population extinction. While Haplosporidium pinnae is characterized as one of the major causative agents, mass mortalities have also been attributed to Mycobacterium sp. and Vibrio spp. Due to limited knowledge concerning the physiological response of infected P. nobilis specimens against various pathogens, this study's aim was to investigate to pathophysiological response of P. nobilis individuals, originating from mortality events in the Thermaikos Gulf and Lesvos and Limnos islands (Greece), and their correlation to different potential pathogens detected in the diseased animals. In isolated tissues, several cellular stress indicators of the heat shock and immune response, apoptosis and autophagy, were examined. Despite the complexity and limitations in the study of P. nobilis mortality events, the present investigation demonstrates the cumulative negative effect of co-infection additionally with H. pinnae in comparison to the non-presence of haplosporidian parasite. In addition, impacts of global climate change affecting physiological performance and immune responses result in more vulnerable populations in infectious diseases, a phenomenon which may intensify in the future.

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'.

Mechanisms for sperm mitochondrial removal in embryos

BACKGROUND

Different animal species have different characteristics regarding the transmission of mitochondrial DNA. While some species have biparental mitochondrial inheritance, others have developed pathways to remove paternal mtDNA. These pathways guarantee the uniparental mitochondrial inheritance, so far well known in mammals, avoiding heteroplasmy, which may have the potential to cause certain mitochondrial diseases in the offspring.

SCOPE OF REVIEW

This review aims to address the main mechanisms that involve mitochondrial degradation in different animal species, as well as to describe what is present in the literature on the mechanisms involved in mitochondrial inheritance.

MAJOR CONCLUSIONS

Two theories are proposed to explain the uniparental inheritance of mtDNA: (i) active degradation, where mechanisms for paternal mitochondrial DNA elimination involve mitochondrial degradation pathway by autophagy and, in some species, may also involve the endocytic degradation pathway; and (ii) passive dilution, where the paternal mitochondria are diluted in the cells of the embryo according to cell division, until becoming undetectable.

GENERAL SIGNIFICANCE

This work brings a wide review of the already published evidence on mitochondrial inheritance in the animal kingdom and the possible mechanisms to mtDNA transmission already described in literature.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

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.

[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.

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.