Mitochondrial phylogenomics of the Bivalvia (Mollusca): searching for the origin and mitogenomic correlates of doubly uniparental inheritance of mtDNA

The effect of mitochondrial recombination on fertilization success in blue mussels

The presence of doubly uniparental inheritance (DUI) in bivalves represents a unique mode of mitochondrial transmission, whereby paternal (male-transmitted M-type) and maternal (female-transmitted F-type) haplotypes are transmitted to offspring separately. Male embryos retain both haplotypes, but the M-type is selectively removed from females. Due to the presence of heteroplasmy in males, mtDNA can recombine resulting in a 'masculinized' haplotype referred to as Mf-type. While mtDNA recombination is usually rare, it has been recorded in multiple mussel species across the Northern Hemisphere. Given that mitochondria are the powerhouse of the cell, different mtDNA haplotypes may have different selective advantages under diverse environmental conditions. This may be particularly important for sperm fitness and fertilization success. In this study we aimed to i) determine the presence, prevalence of the Mf-type in Australian blue mussels (Mytilus sp.) and ii) investigate the effect of Mf-mtDNA on sperm performance (a fitness correlate). We found a high prevalence of recombined mtDNA (≈35 %) located within the control region of the mitochondrial genome, which occurred only in specimens that contained Southern Hemisphere mtDNA. The presence of two female mitotypes were identified in the studied mussels, one likely originating from the Northern Hemisphere, and the other either representing the endemic M. planulatus species or introduced genotypes from the Southern Hemisphere. Despite having recombination events present in a third of the studied population, analysis of sperm performance indicated no difference in fertilization success related to mitotype.

Oxidative phosphorylation rather than glycolysis is the primary energy source for sperm motility in the mussels Mytilus edulis

Broadcast-spawning marine mussels rely on high sperm motility for successful fertilization in the dynamic seawater environment. Mitochondria are typically considered the primary source of ATP generation via oxidative phosphorylation (OXPHOS); however, the ATP generation pathways of mussel sperm have not been fully characterized. To better understand the importance of both OXPHOS and glycolysis for mussel sperm function, we conducted experiments inhibiting these pathways in sperm from Mytilus edulis. Our results indicate that oligomycin, an inhibitor of the mitochondrial ATP synthase, immediately decreased sperm motility rate, velocity, and ATP content, while 2-deoxy-d-glucose, a glycolysis inhibitor, had no effect. The OXPHOS inhibitor rotenone also partially reduced sperm motility rate and velocity. Interestingly, no evidence was found for the inhibitors' effects on the content of energy-rich compounds (lipids, carbohydrates, and proteins) in the mussels' sperm, indicating only modest energy demand to fuel sperm motility. Based on these findings, we conclude that OXPHOS is the primary energy source for sperm motility in marine mussels. Our study sheds light on the intricacies of mussel sperm physiology and highlights the importance of understanding the energy requirements for successful fertilization in broadcast-spawning marine invertebrates.

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.

Phylogenomic analyses shed light on the relationships of chiton superfamilies and shell-eye evolution

Mollusca is the second-largest animal phylum with over 100,000 extant species representing eight classes. Across 1000 extant species in the class Polyplacophora, chitons have a relatively constrained morphology but with some notable deviations. Several genera possess "shell eyes", i.e., true eyes with a lens and retina that are embedded within the dorsal shells. The phylogeny of the major chiton clades is mostly well established, in a set of superfamily-level and higher level taxa supported by various approaches, including morphological studies, multiple gene markers, mitogenome-phylogeny, and phylotranscriptomic approaches. However, one critical lineage has remained unclear, namely Schizochiton which was controversially suggested as being the potential independent origin of chiton shell eyes. Here, with the draft genome sequencing of Schizochiton incisus (superfamily Schizochitonoidea) plus assemblies of transcriptome data from other polyplacophorans, we present phylogenetic reconstructions using both mitochondrial genomes and phylogenomic approaches with multiple methods. We found that phylogenetic trees from mitogenomic data are inconsistent, reflecting larger scale confounding factors in molluscan mitogenomes. However, a consistent and robust topology was generated with protein-coding genes using different models and methods. Our results support Schizochitonoidea as the sister group to other Chitonoidea in Chitonina, in agreement with the established classification. Combined with evidence from fossils, our phylogenetic results suggest that the earliest origin of shell eyes is in Schizochitonoidea, and that these structures were also gained secondarily in other genera in Chitonoidea. Our results have generated a holistic review of the internal relationship within Polyplacophora, and a better understanding of the evolution of Polyplacophora.

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.

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

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

Mitochondrial Heteroplasmy and PCR Amplification Bias Lead to Wrong Species Delimitation with High Confidence in the South American and Antarctic Marine Bivalve Aequiyoldia eightsii Species Complex

The species delimitation of the marine bivalve species complex Aequiyoldia eightsii in South America and Antarctica is complicated by mitochondrial heteroplasmy and amplification bias in molecular barcoding. In this study, we compare different data sources (mitochondrial cytochrome c oxidase subunit I (COI) sequences; nuclear and mitochondrial SNPs). Whilst all the data suggest that populations on either side of the Drake Passage belong to different species, the picture is less clear within Antarctic populations, which harbor three distinct mitochondrial lineages (p-dist ≈ 6%) that coexist in populations and in a subset of individuals with heteroplasmy. Standard barcoding procedures lead to amplification bias favoring either haplotype unpredictably and thus overestimate the species richness with high confidence. However, nuclear SNPs show no differentiation akin to the trans-Drake comparison, suggesting that the Antarctic populations represent a single species. Their distinct haplotypes likely evolved during periods of temporary allopatry, whereas recombination eroded similar differentiation patterns in the nuclear genome after secondary contact. Our study highlights the importance of using multiple data sources and careful quality control measures to avoid bias and increase the accuracy of molecular species delimitation. We recommend an active search for mitochondrial heteroplasmy and haplotype-specific primers for amplification in DNA-barcoding studies.

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

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

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.

Are we there yet? Benchmarking low-coverage nanopore long-read sequencing for the assembling of mitochondrial genomes using the vulnerable silky shark Carcharhinus falciformis

Background

Whole mitochondrial genomes are quickly becoming markers of choice for the exploration of within-species genealogical and among-species phylogenetic relationships. Most often, ‘primer walking’ or ‘long PCR’ strategies plus Sanger sequencing or low-pass whole genome sequencing using Illumina short reads are used for the assembling of mitochondrial chromosomes. In this study, we first confirmed that mitochondrial genomes can be sequenced from long reads using nanopore sequencing data exclusively. Next, we examined the accuracy of the long-reads assembled mitochondrial chromosomes when comparing them to a ‘gold’ standard reference mitochondrial chromosome assembled using Illumina short-reads sequencing.

Results

Using a specialized bioinformatics tool, we first produced a short-reads mitochondrial genome assembly for the silky shark C. falciformis with an average base coverage of 9.8x. The complete mitochondrial genome of C. falciformis was 16,705 bp in length and 934 bp shorter than a previously assembled genome (17,639 bp in length) that used bioinformatics tools not specialized for the assembly of mitochondrial chromosomes. Next, low-pass whole genome sequencing using a MinION ONT pocket-sized platform plus customized de-novo and reference-based workflows assembled and circularized a highly accurate mitochondrial genome in the silky shark Carcharhinus falciformis. Indels at the flanks of homopolymer regions explained most of the dissimilarities observed between the ‘gold’ standard reference mitochondrial genome (assembled using Illumina short reads) and each of the long-reads mitochondrial genome assemblies. Although not completely accurate, mitophylogenomics and barcoding analyses (using entire mitogenomes and the D-Loop/Control Region, respectively) suggest that long-reads assembled mitochondrial genomes are reliable for identifying a sequenced individual, such as C. falciformis, and separating the same individual from others belonging to closely related congeneric species.

Conclusions

This study confirms that mitochondrial genomes can be sequenced from long-reads nanopore sequencing data exclusively. With further development, nanopore technology can be used to quickly test in situ mislabeling in the shark fin fishing industry and thus, improve surveillance protocols, law enforcement, and the regulation of this fishery. This study will also assist with the transferring of high-throughput sequencing technology to middle- and low-income countries so that international scientists can explore population genomics in sharks using inclusive research strategies. Lastly, we recommend assembling mitochondrial genomes using specialized assemblers instead of other assemblers developed for bacterial and/or nuclear genomes.

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.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

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.

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.

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.

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

Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates

Bivalves represent valuable taxonomic group for aging studies given their wide variation in longevity (from 1–2 to >500 years). It is well known that aging is associated to the maintenance of Reactive Oxygen Species homeostasis and that mitochondria phenotype and genotype dysfunctions accumulation is a hallmark of these processes. Previous studies have shown that mitochondrial DNA mutation rates are linked to lifespan in vertebrate species, but no study has explored this in invertebrates. To this end, we performed a Bayesian Phylogenetic Covariance model of evolution analysis using 12 mitochondrial protein-coding genes of 76 bivalve species. Three life history traits (maximum longevity, generation time and mean temperature tolerance) were tested against 1) synonymous substitution rates (dS), 2) conservative amino acid replacement rates (Kc) and 3) ratios of radical over conservative amino acid replacement rates (Kr/Kc). Our results confirm the already known correlation between longevity and generation time and show, for the first time in an invertebrate class, a significant negative correlation between dS and longevity. This correlation was not as strong when generation time and mean temperature tolerance variations were also considered in our model (marginal correlation), suggesting a confounding effect of these traits on the relationship between longevity and mtDNA substitution rate. By confirming the negative correlation between dS and longevity previously documented in birds and mammals, our results provide support for a general pattern in substitution rates.

Yes, we can use it: a formal test on the accuracy of low-pass nanopore long-read sequencing for mitophylogenomics and barcoding research using the Caribbean spiny lobster Panulirus argus

Background

Whole mitogenomes or short fragments (i.e., 300–700 bp of the cox1 gene) are the markers of choice for revealing within- and among-species genealogies. Protocols for sequencing and assembling mitogenomes include ‘primer walking’ or ‘long PCR’ followed by Sanger sequencing or Illumina short-read low-coverage whole genome (LC-WGS) sequencing with or without prior enrichment of mitochondrial DNA. The aforementioned strategies assemble complete and accurate mitochondrial genomes but are time consuming and/or expensive. In this study, I first tested whether mitogenomes can be sequenced from long-read nanopore sequencing data exclusively. Second, I explored the accuracy of the long-read assembled genomes by comparing them to a ‘gold’ standard reference mitogenome retrieved from the same individual using Illumina sequencing. Third and lastly, I tested if the long-read assemblies are useful for mitophylogenomics and barcoding research. To accomplish these goals, I used the Caribbean spiny lobster Panulirus argus, an ecologically relevant species in shallow water coral reefs and target of the most lucrative fishery in the greater Caribbean region.

Results

LC-WGS using a MinION ONT device and various de-novo and reference-based assembly pipelines retrieved a complete and highly accurate mitogenome for the Caribbean spiny lobster Panulirus argus. Discordance between each of the long-read assemblies and the reference mitogenome was mostly due to indels at the flanks of homopolymer regions. Although not ‘perfect’, phylogenetic analyses using entire mitogenomes or a fragment of the cox1 gene demonstrated that mitogenomes assembled using long reads reliably identify the sequenced specimen as belonging to P. argus and distinguish it from other related species in the same genus, family, and superorder.

Conclusions

This study serves as a proof-of-concept for the future implementation of in-situ surveillance protocols using the MinION to detect mislabeling in P. argus across its supply chain. Mislabeling detection will improve fishery management in this overexploited lobster. This study will additionally aid in decreasing costs for exploring meta-population connectivity in the Caribbean spiny lobster and will aid with the transfer of genomics technology to low-income countries.

Characterizing transcriptome in female scallop Chlamys farreri provides new insights into the molecular mechanisms of reproductive regulation during ovarian development and spawn

Bivalve mollusks are descendants of an early-Cambrian lineage and have successfully evolved unique strategies for reproduction. Nonetheless, the molecular mechanisms underlying reproductive regulation in mollusks remain to be elucidated. In this study, transcriptomes of ovary at four reproductive stages in female Chlamys farreri were characterized by RNA-Seq. Regarding signaling pathways, ECM-receptor interaction pathway, mTOR signaling pathway, Fanconi anemia pathway, FoxO signaling pathway, Wnt signaling pathway and Hedgehog signaling pathway were enriched during ovarian development processes. In addition, pathways related to energy metabolism such as Nitrogen metabolism and Arachidonic acid metabolism were enriched at spawn stage. Interestingly, Neuroactive ligand-receptor interaction was significantly enriched involved in ovarian development and spawn, and indicated the potential functions of nervous system on reproductive regulation in C. farreri. What's more, this study identified and characterized fourteen genes involved in "sex hormones synthesis and regulation", "ovarian development and spawn" and "maternal immunity" during the four reproductive stages in C. farreri. We determined that CYP17 uniquely affected gamete release by influencing the physiological balance among the steroid hormones and showed that receptors of the 5-HT and GABA neurotransmitters were tightly associated with ovarian maturation. Furthermore, to the best of our knowledge, this is the first study to report the maternal effect gene Zar1 in bivalve mollusks, likewise the maternal immunity genes displayed coordinated and cooperative expression during reproductive periods, which strengthened the environmental adaptation mechanisms of bivalves. Taken together, this study provides the first dynamic transcriptomic analysis of C. farreri at four key reproductive stages, which will assist in revealing the molecular mechanisms underlying bivalves on reproductive regulation in ovarian development and spawn.

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.

Morphology of Mitochondria in Syncytial Annelid Female Germ-Line Cyst Visualized by Serial Block-Face SEM

Mitochondria change their morphology and distribution depending on the metabolism and functional state of a cell. Here, we analyzed the mitochondria and selected structures in female germ-line cysts in a representative of clitellate annelids – the white worm Enchytraeus albidus in which each germ cell has one cytoplasmic bridge that connects it to a common cytoplasmic mass. Using serial block-face scanning electron microscopy (SBEM), we prepared three-dimensional ultrastructural reconstructions of the entire selected compartments of a cyst at the advanced stage of oogenesis, i.e. the nurse cell, cytophore, and cytoplasmic bridges of all 16 cells (15 nurse cells and oocyte). We revealed extensive mitochondrial networks in the nurse cells, cytophore and mitochondria that pass through the cytoplasmic bridges, which indicates that a mitochondrial network can extend throughout the entire cyst. The dynamic hyperfusion state was suggested for such mitochondrial aggregations. We measured the mitochondria distribution and revealed their polarized distribution in the nurse cells and more abundant accumulation within the cytophore compared to the nurse cell. A close association of mitochondrial networks with dispersed nuage material, which seems to be the structural equivalent of a Balbiani body, not described in clitellate annelids so far, was also revealed.

A battle for transmission: the cooperative and selfish animal mitochondrial genomes

The mitochondrial genome is an evolutionarily persistent and cooperative component of metazoan cells that contributes to energy production and many other cellular processes. Despite sharing the same host as the nuclear genome, the multi-copy mitochondrial DNA (mtDNA) follows very different rules of replication and transmission, which translate into differences in the patterns of selection. On one hand, mtDNA is dependent on the host for its transmission, so selections would favour genomes that boost organismal fitness. On the other hand, genetic heterogeneity within an individual allows different mitochondrial genomes to compete for transmission. This intra-organismal competition could select for the best replicator, which does not necessarily give the fittest organisms, resulting in mito-nuclear conflict. In this review, we discuss the recent advances in our understanding of the mechanisms and opposing forces governing mtDNA transmission and selection in bilaterians, and what the implications of these are for mtDNA evolution and mitochondrial replacement therapy.

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.

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.

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.

Genetic diversity analysis of Peking gecko (Gekko swinhonis) in mid-Eastern China based on mitochondrial COI and Cyt b gene sequences

To understand the genetic diversity of Peking gecko (Gekko swinhonis) populations in its endemic region, 60 individuals were sampled from Lushan, Qi, and Linying counties in Henan Province, China. Through PCR amplification and Sanger sequencing, 120 sequences with lengths of 652 bp (COI) and 739 bp (Cyt b) were obtained, and nine haplotypes were detected for each gene. Overall, results indicated that Peking gecko populations in China have high genetic diversity and significant genetic differentiation. This study provides necessary scientific basis for the protection of Peking gecko germplasm resources.

The impacts of warming on the toxicity of carbon nanotubes in mussels

With the increased production and research on nanoparticles, the presence of carbon nanotubes (CNTs) in aquatic systems is very likely to increase. Although it has been shown that CNTs may cause toxicity in marine organisms, to our knowledge, the possible impacts under global temperature increase is still unknown. For this reason, biochemical and physiological impacts induced in Mytilus galloprovincialis due to the presence of functionalized multi-walled CNTs (f-MWCNTs) and increased temperature were investigated in the present study. The mussels exposed to increased temperature alone presented higher metabolic capacity and expenditure of glycogen as an energy resource to fuel up defense mechanisms and thus preventing oxidative damage. Contrarily, organisms exposed to f-MWCNTs alone seemed not stressed enough to demonstrate differences in the metabolism capacity. Furthermore, f-MWCNTs seemed not able to significantly activate their antioxidant and biotransformation enzymes, which in turn may led to oxidative damage in the cells especially when organisms were exposed to a warmer temperature. In fact, at higher temperature, the antioxidant response of organisms exposed to f-MWCNTs was not effective and oxidative damage levels were observed. Nevertheless, no additive or synergetic effects were observed when mussels were exposed to both stressors simultaneously.

Transcriptome analysis of air-breathing land slug, Incilaria fruhstorferi reveals functional insights into growth, immunity, and reproduction

Background

Incilaria (= Meghimatium) fruhstorferi is an air-breathing land slug found in restricted habitats of Japan, Taiwan and selected provinces of South Korea (Jeju, Chuncheon, Busan, and Deokjeokdo). The species is on a decline due to depletion of forest cover, predation by natural enemies, and collection. To facilitate the conservation of the species, it is important to decide on a number of traits related to growth, immunity and reproduction addressing fitness advantage of the species.

Results

The visceral mass transcriptome of I. fruhstorferi was enabled using the Illumina HiSeq 4000 sequencing platform. According to BUSCO (Benchmarking Universal Single-Copy Orthologs) method, the transcriptome was considered complete with 91.8% of ortholog genes present (Single: 70.7%; Duplicated: 21.1%). A total of 96.79% of the raw read sequences were processed as clean reads. TransDecoder identified 197,271 contigs that contained candidate-coding regions. Of a total of 50,230 unigenes, 34,470 (68.62% of the total unigenes) annotated to homologous proteins in the Protostome database (PANM-DB). The GO term and KEGG pathway analysis indicated genes involved in metabolism, phosphatidylinositol signalling system, aminobenzoate degradation, and T-cell receptor signalling pathway. Many genes associated with molluscan innate immunity were categorized under pathogen recognition receptor, TLR signalling pathway, MyD88 dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related. The reproduction-associated unigenes showed homology to protein fem-1, spermatogenesis-associated protein, sperm associated antigen, and testis expressed sequences, among others. In addition, we identified key growth-related genes categorized under somatotrophic axis, muscle growth, chitinases and collagens. A total of 4822 Simple Sequence Repeats (SSRs) were also identified from the unigene sequences of I. fruhstorferi.

Conclusions

This is the first available genomic information for non-model land slug, I. fruhstorferi focusing on genes related to growth, immunity, and reproduction, with additional focus on microsatellites and repeating elements. The transcriptome provides access to greater number of traits of unknown relevance in the species that could be exploited for in-depth analyses of evolutionary plasticity and making informed choices during conservation planning. This would be appropriate for understanding the dynamics of the species on a priority basis considering the ecological, health, and social benefits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5526-3) contains supplementary material, which is available to authorized users.

Towards a global phylogeny of freshwater mussels (Bivalvia: Unionida): Species delimitation of Chinese taxa, mitochondrial phylogenomics, and diversification patterns

The Yangtze River Basin in China is one of the global hotspots of freshwater mussel (order Unionida) diversity with 68 nominal species. Few studies have tested the validity of these nominal species. Some taxa from the Yangtze unionid fauna have not been adequately examined using molecular data and well-positioned phylogenetically with respect to the global Unionida. We evaluated species boundaries of Chinese freshwater mussels, and disentangled their phylogenetic relationships within the context of the global freshwater mussels based on the multi-locus data and complete mitochondrial genomes. Moreover, we produced the time-calibrated phylogeny of Unionida and explored patterns of diversification. COI barcode data suggested the existence of 41 phylogenetic distinct species from our sampled 40 nominal taxa inhabiting the middle and lower reaches of the Yangtze River. Maximum likelihood and Bayesian inference analyses on three loci (COI, 16S, and 28S) and complete mitochondrial genomes showed that the subfamily Unioninae sensu stricto was paraphyletic, and the subfamily Anodontinae should be subsumed under Unioninae. In addition, we described two new tribes (Aculamprotulini tribe nov. and Lepidodesmini tribe nov.) in the subfamily Unioninae and one new genus (Parvasolenaiagen. nov.) in the subfamily Gonideinae. Molecular dating analysis suggested freshwater mussels diversified at 346.1 Mya (HPD = 286.6-409.9). The global diversification rate for Unionida was estimated to be 0.025 species/Myr. Our study found only a single well-supported rate shift in Unionida diversification, occurring at the base of the subfamily Ambleminae. The evolution of active host-attraction may have triggered the burst of speciation in Ambleminae, and the environment and geography of the Mississippi River Basin likely sustained this radiation.

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.

Evaluating the utility of the female-specific mitochondrial f-orf gene for population genetic, phylogeographic and systematic studies in freshwater mussels (Bivalvia: Unionida)

Freshwater mussels (order: Unionida) represent one of the most critically imperilled groups of animals; consequently, there exists a need to establish a variety of molecular markers for population genetics and systematic studies in this group. Recently, two novel mitochondrial protein-coding genes were described in unionoids with doubly uniparental inheritance of mtDNA. These genes are the f-orf in female-transmitted mtDNA and the m-orf in male-transmitted mtDNA. In this study, whole F-type mitochondrial genome sequences of two morphologically similar Lampsilis spp. were compared to identify the most divergent protein-coding regions, including the f-orf gene, and evaluate its utility for population genetic and phylogeographic studies in the subfamily Ambleminae. We also tested whether the f-orf gene is phylogenetically informative at the species level. Our preliminary results indicated that the f-orf gene could represent a viable molecular marker for population- and species-level studies in freshwater mussels.

Distribution and frequency of mitochondrial DNA polymorphisms in blue mussel (Mytilus edulis) populations of southwestern Nova Scotia (Canada)

The Atlantic blue mussel (Mytilus edulis Linnaeus, 1758) exhibits doubly uniparental inheritance of mitochondrial (mt) DNA. Females are usually homoplasmic for a female-transmitted mt genome (the F type) and males are heteroplasmic for an F type and a male-transmitted mt genome (the M type). F types can undergo “role-reversal” events, resulting in new male-transmitted mtDNA genomes known as recently masculinized (RM) types that co-occur in populations with evolutionarily older standard-male (SM) types. Phylogenetic analyses have shown that RM types periodically replace SM types. It has also been shown that sperm with RM mtDNA have greater swimming velocity and more efficient components of the electron transport chain compared to sperm with SM mtDNA, thus leading to the hypothesis that RM sperm may have a selective advantage over SM sperm. The present study examines the distribution of RM and SM mitotypes in male M. edulis (n = 225) from 13 localities in southwestern Nova Scotia (Canada). The SM type was more common in all populations, with the proportion of RM types ranging from 0% to 24.1%. The highest proportion of RM types was observed in an aquaculture operation. Analyses of additional populations are required to evaluate the selective pressures affecting the geographic distribution of RM and SM mitotypes in M. edulis.

Mitogenomics reveals phylogenetic relationships of caudofoveate aplacophoran molluscs

The worm-shaped, shell-less aplacophoran molluscs Caudofoveata and Solenogastres have recently received attention as part of the clade Aculifera, but relationships within these two lineages are still largely unknown. Here, we use complete mitochondrial genomes to shed light on higher-level relationships within Caudofoveata. Mitochondrial genomes have been sequenced for many diverse molluscs, but only two mitochondrial genomes from aplacophoran molluscs (the caudofoveates Scutopus ventrolineatus and Chaetoderma nitidulum) are available to date. We sequenced and assembled complete or near complete mitochondrial genomes of five additional species of Caudofoveata (Falcidens acutargatus, Falcidens halanychi, Scutopus robustus, Psilodens balduri and Spathoderma clenchi) and one species of Solenogastres (Neomenia carinata) for comparison to available mitochondrial genomes of aculiferans. Comparison of mitochondrial gene order among different lineages revealed a highly conserved order of protein coding genes corresponding to the hypothesized ancestral gene order for Mollusca. Unique arrangements of tRNAs were found among lineages of aculiferan molluscs as well as among caudofoveate taxa. Phylogenetic analyses of amino acid sequences for the 13 protein-coding genes recovered a monophyletic Aplacophora. Within Caudofoveata, Chaetodermatidae, but not Limifossoridae, was recovered monophyletic. Taken together, our results suggest that mitochondrial genomes can serve as useful molecular markers for aculiferan phylogenetics, especially for more recent phylogenetic events.

Ancient Occasional Host Switching of Maternally Transmitted Bacterial Symbionts of Chemosynthetic Vesicomyid Clams

Vesicomyid clams in deep-sea chemosynthetic ecosystems harbor sulfur-oxidizing bacteria in their gill epithelial cells. These symbionts, which are vertically transmitted, are species-specific and thought to have cospeciated with their hosts. However, recent studies indicate incongruent phylogenies between some vesicomyid clams and their symbionts, suggesting that symbionts are horizontally transmitted. To more precisely understand the evolution of vesicomyid clams and their symbionts, we compared the evolution of vesicomyid clams and their symbionts through phylogenetic analyses using multi-gene data sets. Many clades in the phylogenetic trees of 13 host species (Abyssogena mariana, Ab. phaseoliformis, Akebiconcha kawamurai, Calyptogena fausta, C. laubieri, C. magnifica, C. nautilei, C. pacifica, Isorropodon fossajaponicum, Phreagena kilmeri, Ph. okutanii, Ph. soyoae, and Pliocardia stearnsii) and their symbionts were well resolved. Six of the 13 host-symbiont pairs (C. fausta, C. magnifica, C. pacifica, Ph. kilmeri, Ph. okutanii, and Ph. soyoae, and their respective symbionts) showed topological congruence. However, the remaining seven pairs (Ak. kawamurai, Ab mariana, Ab. phaseoliformis, C. laubieri, C. nautilei, I. fossajaponicum, and Pl. stearnsii and their corresponding symbionts) showed incongruent topologies, which were supported by the approximately unbiased and Bayes factor tests. Coevolution analyses indicated that six pairs cospeciated, whereas host switching events occurred in the remaining seven pairs. Markedly, multiple host switching events may have occurred in the lineages from the common ancestral symbiont of C. pacifica and C. fausta. Our phylogenetic and coevolution analyses provide additional evidence for host switching during the evolution of vesicomyids.

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.

Genetic structure and diversity of Nodularia douglasiae (Bivalvia: Unionida) from the middle and lower Yangtze River drainage

The Yangtze River drainage in China is among the most species rich rivers for freshwater mussels (order Unionida) on Earth with at least 68 species known. The freshwater mussels of the Yangtze River face a variety of threats with indications that species are declining in abundance and area of occupancy. This study represents the first analyses of the genetic structure and diversity for the common and widespread freshwater mussel Nodularia douglasiae based on microsatellite DNA genotypes and mitochondrial DNA sequences. Phylogenetic analysis a fragment of the COI mitochondrial gene indicated that N. douglasiae collected from across the middle and lower Yangtze River drainage are monophyletic with N. douglasiae from Japan, Russia, and South Korea. The results of the analysis of both the mtDNA and microsatellite datasets indicated that the seven collection locations of N. douglasiae in the middle and lower Yangtze River drainage showed high genetic diversity, significant genetic differentiation and genetic structure, and stable population dynamics over time. Moreover, we found that the connections among tributaries rivers and lakes in the Yangtze River drainage were important in maintaining gene flow among locations that N. douglasiae inhabits. An understanding of the genetic structure and diversity of a widespread species like N. douglasiae could be used as a surrogate to better understand the populations of other freshwater mussel species that are more rare in the Yangtze River drainage. At the same time, these results could provide a basis for the protection of genetic diversity and management of unionid mussels diversity and other aquatic organisms in the system.

Biochemical impacts of Hg in Mytilus galloprovincialis under present and predicted warming scenarios

The interest in the consequences of climate change on the physiological and biochemical functioning of marine organisms is increasing, but the indirect and interactive effects resulting from warming on bioconcentration and responsiveness to pollutants are still poorly explored, particularly in terms of cellular responses. The present study investigated the impacts of Hg in Mytilus galloprovincialis under control (17°C) and warming (21°C) conditions, assessing mussels Hg bioconcentration capacity, metabolic and oxidative status after 14 and 28days of exposure. Results obtained showed greater impacts in mussels exposed for 28days in comparison to 14days of exposure. Furthermore, our findings revealed that the increase in temperature from 17 to 21°C reduced the bioconcentration of Hg by M. galloprovincialis, which may explain higher mortality rates at 17°C in comparison to 21°C. Lower Hg concentration at 21°C in mussels tissue may result from valves closure for longer periods, identified by reduced energy reserves consumption at higher temperature, which in turn might also contributed to higher oxidative stress in organisms exposed to this condition. The highest LPO levels observed in mussels exposed to higher temperatures alone indicate that warming conditions will greatly affect M. galloprovincialis. Furthermore, the present study showed that the impacts induced by the combination of Hg and warming were similar to the ones caused by increased temperature acting alone, mainly due to increased antioxidant defenses in organisms under combined effects of Hg and warming, suggesting that warming was the factor that mostly contributed to oxidative stress in mussels. Although higher mortality was observed in individuals exposed to 17°C and Hg compared to organisms exposed to Hg at 21°C, the oxidative stress induced at higher temperature may generate negative consequences on mussels reproductive and feeding capacity, growth and, consequently, on population maintenance and dynamics.

Hermaphroditic freshwater mussel Anodonta cygnea does not have supranumerary open reading frames in the mitogenome

The complete mitogenome of Anodonta cygnea is 15,613 bp long. This compact, circular molecule contains the set of 37 genes, typical for invertebrate mitogenomes, in the same order and orientation as in maternally inherited genomes of other bivalves from the same subfamily. There are only two unassigned regions longer than 200 bp (266 bp and 274 bp) and no indication of any supranumerary open reading frames.

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

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

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.