Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution.

BACKGROUND: Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics. RESULTS: We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads. CONCLUSIONS: We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

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

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

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 dN /dS 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.

SmithRNAs: Could Mitochondria "Bend" Nuclear Regulation?

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

Germ cell line during the seasonal sexual rest of clams: finding niches of cells for gonad renewal.

Reconstitution and renewal of tissues are key topics in developmental biology. In this brief work, we analyzed the wintry spent phase of the reproductive cycle in the Manila clam Ruditapes philippinarum (Bivalvia, Veneridae) in order to study the gonad rebuilding that in this species occurs at the beginning of the warmer months. We labeled VASA homolog protein-a germ cell marker-and compared the histological observations of the spent phase with those of the previously analyzed gametogenic phase. In R. philippinarum, during the reproductive season, most of the body mass is represented by sack-like structures (acini) full of developing gametes. In that period, VASA-stained cells are present at the basal pole of the gut epithelium, in the connective tissue, and around the acini. We here show that during the spent phase large portions of the intestine lack such cell type, except for some areas showing a few faintly VASA-stained cells. Cells with similar nuclear morphology are present among loosely organized cells of connective tissue, sometimes as single units, sometimes in small groups, rarely partially organized in primordial gonadic structures. These observations match the findings of RNA-targeting studies that during the spent phase identified the source of bivalve germ cells within the connective tissue in the form of quiescent units and add new information on the possible maintenance of VASA-stained, multipotent cells among the batiprismatic cells of the intestine during the whole life span of these bivalves.

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.

A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs.

Recent data from mitochondrial genomics and proteomics research demonstrate the existence of several atypical mitochondrial protein-coding genes (other than the standard set of 13) and the involvement of mtDNA-encoded proteins in functions other than energy production in several animal species including humans. These results are of considerable importance for evolutionary and cellular biology because they indicate that animal mtDNAs have a larger functional repertoire than previously believed. This review summarizes recent studies on animal species with a non-standard mitochondrial functional repertoire and discusses how these genetic novelties represent promising candidates for studying the role of the mitochondrial genome in speciation.

VASA expression suggests shared germ line dynamics in bivalve molluscs.

Germ line segregation can occur during embryogenesis or after embryogenesis completion, with multipotent cells able to give rise to both germ and somatic cells in the developing juvenile or even in adulthood. These undifferentiated cells, in some animals, are self-renewing stem cells. In all these cell lineages, the same set of genes, among which vasa, appears to be expressed. We traced VASA expression during the peculiar gonad rebuilding of bivalves to verify its presence from undifferentiated germ cells to mature gametes in an animal taxon in which the mechanism of germ line establishment is still under investigation. We utilized antibodies produced against VASPH, VASA homolog of Ruditapes philippinarum (Subclass Heterodonta), to compare the known expression pattern of R. philippinarum to two species of the Subclass Pteriomorphia, Anadara kagoshimensis and Crassostrea gigas, and another species of the Subclass Heterodonta, Mya arenaria. The immunohistological data obtained support a conserved mechanism of proliferation of "primordial stem cells" among the simple columnar epithelium of the gut, as well as in the connective tissue, contributing to the seasonal gonad reconstitution. Given the taxonomic separation of the analyzed species, we suggest that the process could be shared in bivalve molluscs. The presence of germ cell precursors in the gut epithelium appears to be a feature in common with model organisms, such as mouse, fruit fly, and human. Thus, the comparative study of germ line establishment can add details on bivalve development, but can also help to clarify the role that VASA plays during germ cell specification.

Evidence for somatic transcription of male-transmitted mitochondrial genome in the DUI species Ruditapes philippinarum (Bivalvia: Veneridae).

In species with doubly uniparental inheritance (DUI), males are heteroplasmic for two sex-linked mitochondrial genomes (M- and F-mtDNA). While a role of M-mtDNA in male gametogenesis and sperm function is evident, there is an ongoing debate on whether it is transcribed or not in male soma. In this work we report a qPCR analysis in the DUI species Ruditapes philippinarum, showing that M-mtDNA is transcribed in somatic tissues. We observed a correlation between DNA copy numbers of the two analyzed genes, cytochrome b and a novel male-specific mitochondrial gene thought to be involved in DUI (orf21), and between their transcription levels. No correlation between a transcript and its DNA copy number was found, supporting the existence of complex regulatory mechanisms of mitochondrial transcription. We found the highest amount of mtDNA and mtRNA in gonads, likely due to the intense cell proliferation and high energy request for gametogenesis, while the observed variation among specimens is probably related to their different stages of gonad development. Finally, orf21 showed a highly variable transcription in advanced stages of gametogenesis. We hypothesize a differential storage of orf21 transcripts in spermatozoa, representing different paternal contributions to progeny, possibly leading to different developmental outcomes. A transcriptional activity does not necessarily imply the translation of M-mtDNA genes, and studies on mitochondrial proteins and their localization are needed to definitively assess the functioning of male-transmitted mitochondria in male soma. All that considered, the male soma of DUI species may represent an intriguing experimental model to study cytoplasmic genetic conflicts.

De Novo assembly of the Manila clam Ruditapes philippinarum transcriptome provides new insights into expression bias, mitochondrial doubly uniparental inheritance and sex determination.

Males and females share the same genome, thus, phenotypic divergence requires differential gene expression and sex-specific regulation. Accordingly, the analysis of expression patterns is pivotal to the understanding of sex determination mechanisms. Many bivalves are stable gonochoric species, but the mechanism of gonad sexualization and the genes involved are still unknown. Moreover, during the period of sexual rest, a gonad is not present and sex cannot be determined. A mechanism associated with germ line differentiation in some bivalves, including the Manila clam Ruditapes philippinarum, is the doubly uniparental inheritance (DUI) of mitochondria, a variation of strict maternal inheritance. Two mitochondrial lineages are present, one transmitted through eggs and the other through sperm, as well as a mother-dependent sex bias of the progeny. We produced a de novo annotation of 17,186 transcripts from R. philippinarum and compared the transcriptomes of males and females and identified 1,575 genes with strong sex-specific expression and 166 sex-specific single nucleotide polymorphisms, obtaining preliminary information about genes that could be involved in sex determination. Then we compared the transcriptomes between a family producing predominantly females and a family producing predominantly males to identify candidate genes involved in regulation of sex-specific aspects of DUI system, finding a relationship between sex bias and differential expression of several ubiquitination genes. In mammalian embryos, sperm mitochondria are degraded by ubiquitination. A modification of this mechanism is hypothesized to be responsible for the retention of sperm mitochondria in male embryos of DUI species. Ubiquitination can additionally regulate gene expression, playing a role in sex determination of several animals. These data enable us to develop a model that incorporates both the DUI literature and our new findings.

Nuclear genes with sex bias in Ruditapes philippinarum (Bivalvia, veneridae): Mitochondrial inheritance and sex determination in DUI species.

Mitochondria are inherited maternally in most metazoans, but in bivalves with Doubly Uniparental Inheritance (DUI) a mitochondrial lineage is transmitted through eggs (F-type), and another through sperm (M-type). In DUI species, a sex-ratio distortion of the progeny was observed: some females produce a female-biased offspring (female-biased family), others a male-biased progeny (male-biased family), and others a 50:50 sex-ratio. A peculiar segregation pattern of M-type mitochondria in DUI organisms appears to be correlated with the sex bias of these families. According to a proposed model for the inheritance of M-type mitochondria in DUI, the transmission of sperm mitochondria is controlled by three nuclear genes, named W, X, and Z. An additional S gene with different dosage effect would be involved in sex determination. In this study, we analyzed structure and localization of three transcripts (psa, birc, and anubl1) with specific sex and family biases in the Manila clam Ruditapes philippinarum. In situ hybridization confirmed the localization of these transcripts in gametogenic cells. In other animals, homologs of these genes are involved in reproduction and ubiquitination. We hypothesized that these genes may have a role in sex determination and could also be responsible for the maintenance/degradation of spermatozoon mitochondria during embryo development of the DUI species R. philippinarum, so that we propose them as candidate factors of the W/X/Z/S system.

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.

Signatures of Extreme Longevity: A Perspective from Bivalve Molecular Evolution.

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

Strict sex-specific mtDNA segregation in the germ line of the DUI species Venerupis philippinarum (Bivalvia: Veneridae).

Doubly Uniparental Inheritance (DUI) is one of the most striking exceptions to the common rule of standard maternal inheritance of metazoan mitochondria. In DUI, two mitochondrial genomes are present, showing different transmission routes, one through eggs (F-type) and the other through sperm (M-type). In this paper, we report results from a multiplex real-time quantitative polymerase chain reaction analysis on the Manila clam Venerupis philippinarum (formerly Tapes philippinarum). We quantified M- and F-types in somatic tissues, gonads, and gametes. Nuclear and external reference sequences were used, and the whole experimental process was designed to avoid any possible cross-contamination. In most male somatic tissues, the M-type is largely predominant: This suggests that the processes separating sex-linked mitochondrial DNAs (mtDNAs) in somatic tissues are less precise than in other DUI species. In the germ line, we evidenced a strict sex-specific mtDNA segregation because both sperm and eggs do carry exclusively M- and F-types, respectively, an observation that is in contrast with a previous analysis on Mytilus galloprovincialis. More precisely, whereas two mtDNAs are present in the whole gonad, only the sex-specific one is detected in gametes. Because of this, we propose that the mtDNA transmission is achieved through a three-checkpoint process in V. philippinarum. The cytological mechanisms of male mitochondria segregation in males and degradation in females during the embryo development (here named Checkpoint #1 and Checkpoint #2) are already well known for DUI species; a Checkpoint #3 would act when primordial germ cells (PGCs) are first formed and would work in both males and females. We believe that Checkpoint #3 is a mere variation of the "mitochondrial bottleneck" in species with standard maternal inheritance, established when their PGCs separate during embryo cleavage.

Sex-linked mitochondrial behavior during early embryo development in Ruditapes philippinarum (Bivalvia Veneridae) a species with the Doubly Uniparental Inheritance (DUI) of mitochondria.

In most metazoans mitochondria are inherited maternally. However, in some bivalve molluscs, two mitochondrial lineages are present: one transmitted through females (F-type), the other through males (M-type). This unique system is called Doubly Uniparental Inheritance (DUI) of mitochondria. In DUI species, M-type mitochondria have to invade the germ line of male embryos during development, otherwise sperm would transmit F-type mtDNA and DUI would fail. The mechanisms by which sperm mitochondria enter the germ line are still unknown. To address this question, we traced the movement of spermatozoon mitochondria (M-type) in embryos of the DUI species Ruditapes philippinarum by fertilizing eggs with sperm stained with the mitochondrial-specific vital dye MitoTracker Green. As in Mytilus DUI species, in R. philippinarum the distribution of sperm mitochondria follows two different patterns: an aggregated one in which these organelles locate near the first cleavage furrow, and a dispersed one in which sperm mitochondria are scattered. The presence of the two mitochondrial patterns in these taxa, together with their absence in species with Strictly Maternal Inheritance (SMI), confirms that their occurrence is related to DUI. Moreover, a Real-Time qPCR analysis showed that neither M-type nor F-type mitochondria undergo replication boosts in the earliest embryo development. This is the first study on sex-linked mtDNA copy number carried out by qPCR analysis on embryos of a DUI species and the first time the segregation patterns of sperm mitochondria are described in a DUI system other than Mytilus.

Repeat Age Decomposition Informs an Ancient Set of Repeats Associated With Coleoid Cephalopod Divergence.

In comparison with other molluscs and bilaterians, the genomes of coleoid cephalopods (squid, cuttlefish, and octopus) sequenced so far show remarkably different genomic organization that presumably marked the early evolution of this taxon. The main driver behind this genomic rearrangement remains unclear. About half of the genome content in coleoids is known to consist of repeat elements; since selfish DNA is one of the powerful drivers of genome evolution, its pervasiveness could be intertwined with the emergence of cephalopod-specific genomic signatures and could have played an important role in the reorganization of the cephalopod genome architecture. However, due to abundant species-specific repeat expansions, it has not been possible so far to identify the ancient shared set of repeats associated with coleoid divergence. By means of an extensive repeat element re-evaluation and annotation combined with network sequence divergence approaches, we are able to identify and characterize the ancient repeat complement shared by at least four coleoid cephalopod species. Surprisingly, instead of the most abundant elements present in extant genomes, lower-copy-number DNA and retroelements were most associated with ancient coleoid radiation. Furthermore, evolutionary analysis of some of the most abundant families shared in Octopus bimaculoides and Euprymna scolopes disclosed within-family patterns of large species-specific expansions while also identifying a smaller shared expansion in the coleoid ancestor. Our study thus reveals the apomorphic nature of retroelement expansion in octopus and a conserved complement composed of several DNA element types and fewer LINE families.

Lack of transcriptional coordination between mitochondrial and nuclear oxidative phosphorylation genes in the presence of two divergent mitochondrial genomes.

In most eukaryotes, oxidative phosphorylation (OXPHOS) is the main energy production process and it involves both mitochondrial and nuclear genomes. The close interaction between the two genomes is critical for the coordinated function of the OXPHOS process. Some bivalves show doubly uniparental inheritance (DUI) of mitochondria, where two highly divergent mitochondrial genomes, one inherited through eggs (F-type) and the other through sperm (M-type), coexist in the same individual. However, it remains a puzzle how nuclear OXPHOS genes coordinate with two divergent mitochondrial genomes in DUI species. In this study, we compared transcription, polymorphism, and synonymous codon usage in the mitochondrial and nuclear OXPHOS genes of the DUI species Ruditapes philippinarum using sex- and tissue-specific transcriptomes. Mitochondrial and nuclear OXPHOS genes showed different transcription profiles. Strong co-transcription signal was observed within mitochondrial (separate for F- and M-type) and within nuclear OXPHOS genes but the signal was weak or absent between mitochondrial and nuclear OXPHOS genes, suggesting that the coordination between mitochondrial and nuclear OXPHOS subunits is not achieved transcriptionally. McDonald-Kreitman and frequency-spectrum based tests indicated that M-type OXPHOS genes deviated significantly from neutrality, and that F-type and M-type OXPHOS genes undergo different selection patterns. Codon usage analysis revealed that mutation bias and translational selection were the major factors affecting the codon usage bias in different OXPHOS genes, nevertheless, translational selection in mitochondrial OXPHOS genes appears to be less efficient than nuclear OXPHOS genes. Therefore, we speculate that the coordination between OXPHOS genes may involve post-transcriptional/translational regulation.

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.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae).

BACKGROUND: Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species. RESULTS: Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer. CONCLUSIONS: The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

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