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.

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Mitochondrial genomic analyses provide new insights into the "missing" atp8 and adaptive evolution of Mytilidae

Background

Mytilidae, also known as marine mussels, are widely distributed in the oceans worldwide. Members of Mytilidae show a tremendous range of ecological adaptions, from the species distributed in freshwater to those that inhabit in deep-sea. Mitochondria play an important role in energy metabolism, which might contribute to the adaptation of Mytilidae to different environments. In addition, some bivalve species are thought to lack the mitochondrial protein-coding gene ATP synthase F0 subunit 8. Increasing studies indicated that the absence of atp8 may be caused by annotation difficulties for atp8 gene is characterized by highly divergent, variable length.

Results

In this study, the complete mitochondrial genomes of three marine mussels (Xenostrobus securis, Bathymodiolus puteoserpentis, Gigantidas vrijenhoeki) were newly assembled, with the lengths of 14,972 bp, 20,482, and 17,786 bp, respectively. We annotated atp8 in the sequences that we assembled and the sequences lacking atp8. The newly annotated atp8 sequences all have one predicted transmembrane domain, a similar hydropathy profile, as well as the C-terminal region with positively charged amino acids. Furthermore, we reconstructed the phylogenetic trees and performed positive selection analysis. The results showed that the deep-sea bathymodiolines experienced more relaxed evolutionary constraints. And signatures of positive selection were detected in nad4 of Limnoperna fortunei, which may contribute to the survival and/or thriving of this species in freshwater.

Conclusions

Our analysis supported that atp8 may not be missing in the Mytilidae. And our results provided evidence that the mitochondrial genes may contribute to the adaptation of Mytilidae to different environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08940-8.

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

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

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.

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.

Transoceanic dispersal of the mussel Mytilus galloprovincialis on Japanese tsunami marine debris: An approach for evaluating rafting of a coastal species at sea

Biofouled debris from the 2011 Great East Japan earthquake and tsunami has landed in the Northeast Pacific and along the Hawaiian Islands since 2012. As of 2017, >630 biofouled debris items with >320 living species of algae, invertebrates, and fish have been examined. The invasive mussel Mytilus galloprovincialis was present on >50% of those items. Size, reproduction, and growth of this filter-feeding species were examined to better understand long-distance rafting of a coastal species. The majority of mussels (79%) had developing or mature gametes, and growth rates averaged 0.075±0.018 SE mm/day. Structural and elemental (barium/calcium) analysis of mussel shells generated estimates of growth in coastal waters (mean=1.3 to 25mm total length), which provides an indication of residence times in waters along North America and the Hawaiian Islands prior to landing. Detailed studies of individual species contribute to our understanding of debris as a transport vector and aid efforts to evaluate potential risks associated with marine debris.

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.

Complete female mitochondrial genome of Mytilus chilensis

The controversy surrounding the origin of antitropical distribution of Mytilus mussels and the taxonomic status of southern hemisphere populations remain unsolved, despite the efforts. One of the limiting factors remains the lack of the complete sequences of the representative mitochondrial genomes which would allow their proper comparison with the relatively well-represented northern hemisphere congeneric mussels. To fill this gap we sequenced the representative maternal (F) genome of a native Chilean mussel. The genome is 16,748bp long and structurally identical to the northern hemisphere M. edulis and M. galloprovincialis F genomes. However, the genetic distance from them (≈5%) is twice as high as the maximum distance between them (<2.5%). Thus, the notion that the name M. chilensis should be used for native Chilean Mytilus mussels, with the same rank as M. galloprovincialis and M. edulis is supported.

Sequence motifs associated with paternal transmission of mitochondrial DNA in the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae)

In the majority of metazoans paternal mitochondria represent evolutionary dead-ends. In many bivalves, however, this paradigm does not hold true; both maternal and paternal mitochondria are inherited. Herein, we characterize maternal and paternal mitochondrial control regions of the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae). The maternal control region is 808bp long, while the paternal control region is longer at 2.3kb. We hypothesize that the size difference is due to a combination of repeated duplications within the control region of the paternal mtDNA genome, as well as an evolutionarily ancient recombination event between two sex-associated mtDNA genomes that led to the insertion of a second control region sequence in the genome that is now transmitted via males. In a comparison to other mytilid male control regions, we identified two evolutionarily Conserved Motifs, CMA and CMB, associated with paternal transmission of mitochondrial DNA. CMA is characterized by a conserved purine/pyrimidine pattern, while CMB exhibits a specific 13bp nucleotide string within a stem and loop structure. The identification of motifs CMA and CMB in M. modiolus extends our understanding of Sperm Transmission Elements (STEs) that have recently been identified as being associated with the paternal transmission of mitochondria in marine bivalves.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

The rRNA and tRNA transcripts of maternally and paternally inherited mitochondrial DNAs of Mytilus galloprovincialis suggest presence of a "degradosome" in mussel mitochondria and necessitate the re-annotation of the l-rRNA/CR boundary

Species of the genus Mytilus carry two mitochondrial genomes in obligatory coexistence; one transmitted though the eggs (the F type) and one through the sperm (the M type). We have studied the 3' and 5' ends of rRNA and tRNA transcripts using RT-PCR and RNA circularization techniques in both the F and M genomes of Mytilus galloprovincialis. We have found polyadenylated and non-adenylated transcripts for both ribosomal and transfer RNAs. In all these genes the 5' ends of the transcripts coincided with the first nucleotide of the annotated genes, but the 3' ends were heterogeneous. The l-rRNA 3' end is 47 or 48 nucleotides upstream from the one assigned by a previous annotation, which makes the adjacent first domain (variable domain one, VD1) of the main control region (CR) correspondingly longer. We have observed s-rRNA and l-rRNA transcripts with truncated 3' end and polyadenylated tRNA transcripts carrying the CCA trinucleotide. We have also detected polyadenylated RNA remnants carrying the sequences of the control region, which strongly suggests RNA degradation activity and thus presence of degradosomes in Mytilus mitochondria.

Evidence for a fourteenth mtDNA-encoded protein in the female-transmitted mtDNA of marine Mussels (Bivalvia: Mytilidae)

Background

A novel feature for animal mitochondrial genomes has been recently established: i.e., the presence of additional, lineage-specific, mtDNA-encoded proteins with functional significance. This feature has been observed in freshwater mussels with doubly uniparental inheritance of mtDNA (DUI). The latter unique system of mtDNA transmission, which also exists in some marine mussels and marine clams, is characterized by one mt genome inherited from the female parent (F mtDNA) and one mt genome inherited from the male parent (M mtDNA). In freshwater mussels, the novel mtDNA-encoded proteins have been shown to be mt genome-specific (i.e., one novel protein for F genomes and one novel protein for M genomes). It has been hypothesized that these novel, F- and M-specific, mtDNA-encoded proteins (and/or other F- and/or M-specific mtDNA sequences) could be responsible for the different modes of mtDNA transmission in bivalves but this remains to be demonstrated.

Methodology/Principal Findings

We investigated all complete (or nearly complete) female- and male-transmitted marine mussel mtDNAs previously sequenced for the presence of ORFs that could have functional importance in these bivalves. Our results confirm the presence of a novel F genome-specific mt ORF, of significant length (>100aa) and located in the control region, that most likely has functional significance in marine mussels. The identification of this ORF in five Mytilus species suggests that it has been maintained in the mytilid lineage (subfamily Mytilinae) for ∼13 million years. Furthermore, this ORF likely has a homologue in the F mt genome of Musculista senhousia, a DUI-containing mytilid species in the subfamily Crenellinae. We present evidence supporting the functionality of this F-specific ORF at the transcriptional, amino acid and nucleotide levels.

Conclusions/Significance

Our results offer support for the hypothesis that “novel F genome-specific mitochondrial genes” are involved in key biological functions in bivalve species with DUI.

Homologous recombination between highly diverged mitochondrial sequences: examples from maternally and paternally transmitted genomes

Homologous recombination is restricted to sequences of low divergence. This is attributed to the mismatch repairing system (MMR), which does not allow recombination between sequences that are highly divergent. This acts as a safeguard against recombination between nonhomologous sequences that could result in genome imbalance. Here, we report recombination between maternal and paternal mitochondrial genomes of the sea mussel, whose sequences differ by >20%. We propose that the strict maternal inheritance of the animal mitochondrial DNA and the ensuing homoplasmy has relieved the MMR system of the animal mitochondrion from the pressure to tolerate recombination only among sequences with a high degree of similarity.

The atypical presence of the paternal mitochondrial DNA in somatic tissues of male and female individuals of the blue mussel species Mytilus galloprovincialis

BACKGROUND

In animals mtDNA inheritance is maternal except in certain molluscan bivalve species which have a paternally inherited mitochondrial genome (genome M) along with the standard maternal one (genome F). Normally, the paternal genome occurs in the male gonad, but it can be often found, as a minority, in somatic tissues of males and females. This may happen in two ways. One is through "sperm mtDNA leakage" into somatic tissues, a deviation from the normal situation in which the sperm mtDNA vanishes in females or ends up exclusively in the germ line of males. The other is through "egg heteroplasmy", when the egg contains, in small quantities, the paternal genome in addition to maternal genome.

FINDINGS

To test the two hypotheses, we compared the sequences of one of the most variable domains of the M molecule in a somatic tissue (foot) and in the sperm of ten male and in the foot of ten female individuals of M. galloprovincialis. Presence of the M genome was rarer in the foot of females than males. The M genome in the sperm and in the foot of males was identical.

CONCLUSIONS

Given that the surveyed region differs from individual to individual, the identity of the M genome in the foot and the sperm of males supports strongly the hypothesis that, at least for the tissue examined, the presence of the M genome is due to sperm mtDNA leakage.

Scottish Mytilus trossulus mussels retain ancestral mitochondrial DNA: complete sequences of male and female mtDNA genomes

Mytilus trossulus mussels occur in North America and in the Baltic Sea. Recently genetic markers for the three Mytilus subspecies M. edulis, M. galloprovincialis, and M. trossulus, have been detected at Loch Etive in Scotland suggesting mixed ancestry for this population. Of particular interest is the evidence that M. trossulus occurs at Loch Etive because it had not previously been reported in the British Isles. In the present study, analysis of subspecies-specific diagnostic nuclear DNA markers confirms the presence of a high frequency of mussels with M. trossulus ancestry at Loch Etive. The genetic structure suggests hybridisation at an intermediate stage compared with North American populations, where there is little hybridisation, and Baltic populations where there is extensive introgression. This points strongly against a Baltic origin for Loch Etive M. trossulus. The F and M mitochondrial DNA (mtDNA) genomes of Baltic M. trossulus are similar in sequence to the corresponding genomes in M. edulis and believed to be derived by introgression from that subspecies. Both F and M mtDNA genomes are observed at Loch Etive consistent with the presence of doubly uniparental inheritance. Here we provide the complete sequences of the three M. trossulus mtDNA genomes (one F and two M) from Loch Etive. These genomes are extremely similar to the corresponding genomes from ancestral M. trossulus in America but divergent from the genomes for Baltic M. trossulus. This is the first report of ancestral M. trossulus mtDNA genomes in Europe. The F and M genomes are diverged by 26% in nucleotide sequence, similar to other Mytilus F and M genomes. The gene arrangement in the sequenced genomes is also similar to that in other sequenced Mytilus mtDNA genomes. However the two sequenced M genomes differ by 960bp which is caused by a duplication in the main noncoding region (CR). This duplication has not so far been observed in North American populations of M. trossulus. The coding regions of the Loch Etive genomes have no features suggesting that they are other than functional genomes and have K(a)/K(s) values in coding regions less than one indicative of purifying selection. Estimates of divergence times were made for both genomes and are consistent with invasion of Loch Etive by M. trossulus towards the end of the last glacial period.