Patterns of polymorphism and gene flow of gender-associated mitochondrial DNA lineages in European mussel populations

Mussels of the genus Mytilus have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally and the F type is transmitted maternally. RFLP analysis is used to assess phylogenetic relationships and nucleotide diversity and divergence for both mtDNA genomes in European populations of M. edulis and Atlantic and Mediterranean forms of M. galloprovincialis. Ten restriction endonucleases were used to assay variation in regions of the ND2 and COIII genes for a total of 77 individuals. F and M genomes show a concordant phylogenetic split into two major divergent clades, one specific to Mediterranean M. galloprovincialis and the other containing haplotypes from the three taxa. For both genomes, the geographical distribution of mtDNA variation suggests: (i) extensive levels of mtDNA introgression; (ii) asymmetric mtDNA gene flow from Atlantic to Mediterranean populations; and (iii) recurrent historical hybridization events. Significantly higher mtDNA diversity and divergence are observed for the M than F genome in all three Mytilus taxa, although the evolutionary forces responsible for these differences cannot be resolved. The extensive mtDNA gene flow between European Mytilus taxa conflicts with the restricted mtDNA introgression observed in American mussels , implying geographical variation in the nature of nuclear/mtDNA interactions regulating biparental inheritance.     

Comparative analysis of gender-associated complete mitochondrial genomes in marine mussels (Mytilus spp.)

Marine mussels of the genus Mytilus have an unusual mode of mitochondrial DNA (mtDNA) transmission termed doubly uniparental inheritance (DUI). Female mussels are homoplasmic for the F mitotype, which is inherited maternally, while males are usually heteroplasmic, carrying a mixture of the maternal F mitotype and the paternally inherited M genome. Two classes of M genomes have been observed: "standard" M genomes and "recently masculinized" M genomes. The latter are more similar to F genomes at the sequence level but are transmitted paternally like standard M genomes. In this study we report the complete sequences of two standard male M. edulis and one recently masculinized male M. trossulus mitochondrial genome. A comparative analysis, including the previously sequenced M. edulis F and M. galloprovincialis F and M mtDNAs, reveals that these genomes are identical in gene order, but highly divergent in nucleotide and amino acid sequence. The large amount (>20%) of nucleotide substitutions that fall in coding regions implies that there are several amino acid replacements between the F and M genomes, which likely have an impact on the structural and functional properties of the mitochondrial proteome. Correlation of the divergence rate of different protein-coding genes indicates that mtDNA-encoded proteins of the M genome are still under selective constraints, although less highly than genes of the F genome. The mosaic F/M control region of the masculinized F genome provides evidence for lineage-specific sequences that may be responsible for the different mode of transmission genetics. This analysis shows the value of comparative genomics to better understand the mechanisms of maintenance and segregation of mtDNA sequence variants in mytilid mussels.     

Glacial history of the European marine mussels Mytilus,inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.     

Nonneutral evolution of tandem repeats in the mitochondrial DNA control region of lagomorphs

The mitochondrial DNA of the European rabbit (Oryctolagus cuniculus) contains a tandem array of 153-bp repeats in the vicinity of the replication origin of the H-stand. Variation among molecules in the number of these repeats results in inter- and intraindividual length polymorphism (heteroplasmy). Generally, in an individual, one predominant molecular type is observed, the others representing a low percentage of the mtDNA content. At the tissue level, we observe a particular distribution of this polymorphism in the gonads compared with liver, kidneys, or brain, implying a relationship between the differentiation status of the cells and the types of new mtDNA molecules which appear and accumulate during lifetime. Similar tandem repeats were also found in the mtDNA noncoding region of European hares (Lepus europaeus), a cottontail (Sylvilagus floridanus), and a pika (Ochotona rufescens). The lengths and the sequences of these units evolve rapidly and in a concerted way, but the number of repeats is maintained in a narrow range, and an internal 20-bp segment is highly conserved. Constraints restrict the evolution of the primary sequence of these repeated units, the number of which is probably controlled by a stabilizing selection.      

Phylogenetic evidence for role-reversals of gender-associated mitochondrial DNA in Mytilus (Bivalvia: Mytilidae)

Distinct gender-associated mitochondrial DNA (mtDNA) lineages (i.e., lineages which are transmitted either through males or through females) have been demonstrated in two families of bivalves, the Mytilidae (marine mussels) and the Unionidae (freshwater mussels), which have been separated for more than 400 Myr. The mode of transmission of these M (for male-transmitted) and F (for female-transmitted) molecules has been referred to as doubly uniparental inheritance (DUI), in contrast to standard maternal inheritance (SMI), which is the norm in animals. A previous study suggested that at least three origins of DUI are required to explain the phylogenetic pattern of M and F lineages in freshwater and marine mussels. Here we present phylogenetic evidence based on partial sequences of the cytochrome c oxidase subunit I gene and the 16S RNA gene that indicates the DUI is a dynamic phenomenon. Specifically, we demonstrate that F lineages in three species of Mytilus mussels, M. edulis, M. trossulus, and M. californianus, have spawned separate lineages which are now associated only with males. This process is referred to as "masculinization" of F mtDNA. By extension, we propose that DUI may be a primitive bivalve character and that periodic masculinization events combined with extinction of previously existing M types effectively reset the time of divergence between conspecific gender-associated mtDNA lineages.      

Direct evidence for homologous recombination in mussel (Mytilus galloprovincialis) mitochondrial DNA.

The assumption that animal mitochondrial DNA (mtDNA) does not undergo homologous recombination is based on indirect evidence, yet it has had an important influence on our understanding of mtDNA repair and mutation accumulation (and thus mitochondrial disease and aging) and on biohistorical inferences made from population data. Recently, several studies have suggested recombination in primate mtDNA on the basis of patterns of frequency distribution and linkage associations of mtDNA mutations in human populations, but others have failed to produce similar evidence. Here, we provide direct evidence for homologous mtDNA recombination in mussels, where heteroplasmy is the rule in males. Our results indicate a high rate of mtDNA recombination. Coupled with the observation that mammalian mitochondria contain the enzymes needed for the catalysis of homologous recombination, these findings suggest that animal mtDNA molecules may recombine regularly and that the extent to which this generates new haplotypes may depend only on the frequency of biparental inheritance of the mitochondrial genome. This generalization must, however, await evidence from animal species with typical maternal mtDNA inheritance.     

Origin and evolution of Mytilus mussel satellite DNAs.

A phylogenetic reconstruction based on the amplification of 3 satellite DNAs (stDNAs) was carried out in 1 crustacean species and 15 bivalve species of the subclass Pteriomorphia (10, subfamily Mytilinae; 1, subfamily Litophaginae; 1, subfamily Modiolinae, all belonging to family Mytilidae; 1, family Arcidae; and 2, family Pectinidae). The sequences obtained showed motifs with high similarity to those of A and B boxes of tRNA promoter regions. Dot-blot hybridizations revealed that the 3 stDNAs are present mainly in high copy numbers for each species of the genus Mytilus, whereas for the other species they appear in low copy numbers. Maximum-parsimony trees evidenced a tendency to group Mytilus clones together, and species containing these sequences as a single copy were distributed among the different mytilids. Finally, the possible origin and evolution of these stDNAs is discussed.     

Algal-induced spawning in the marine mussel Mytilus californianus 4

The characteristics of algal-induced spawning of the marine mussel Mytilus californianus were studied. Exposure of mature individuals to culture suspensions of the unicellular alga Pseudoisochrysis paradoxa elicited copious and synchronous release of gametes. Alkaline conditions were necessary to make the animals responsive to the spawning stimulus provided by the algae. The filtered, cell-free fraction of the algal suspensions also stimulated spawning, suggesting that an active principle is secreted into the culture media by the algae.

A key requirement for the development of an efficient and economically viable molluscan mariculture program is the capability to induce copious, synchronous spawning of gravid individuals. This procedure should be easily managed and have no deleterious effects on the gametes, fertilization, or the ontological development of the animals. Miyazaki [1] induced spawning in male oysters by exposing them to an extract of green algae. Recently it was reported that the marine mussels, M. californianus, could be induced to spawn by exposing them to cultures of the marine algae, P. paradoxa [2], This paper summarizes our efforts to characterize this spawning response of mussels to high concentrations of cultured marine algae.     

Evolutionary relationships among the male and female mitochondrial DNA lineages in the Mytilus edulis species complex

A novel form of mitochondrial DNA (mtDNA) inheritance has previously been documented for the blue mussel (Mytilus edulis). Female mussels inherit their mtDNA solely from their mother while males inherit mtDNA from both their mother and their father. In males, the paternal mtDNA is preferentially amplified so that the male gonad is highly enriched for the paternal mtDNA that is then transmitted from fathers to sons. We demonstrate that this mode of mtDNA inheritance also operates in the closely related species M. galloprovincialis and M. trossulus. The evolutionary relationship between the male and female mtDNA lineages is estimated by phylogenetic analysis of 455 nucleotides from the large subunit ribosomal RNA gene. We have found that the male and female lineages are highly divergent; the divergence of these lineages began prior to the speciation of the three species of blue mussels. Further, the separation between the male and female lineages is estimated to have occurred between 5.3 and 5.7 MYA.      

Sperm mitochondrial DNA transmission to both male and female offspring in the blue mussel Mytilus galloprovincialis

In Mytilus mussels, paternal mitochondrial DNA (M type) from sperm is known to be transmitted to offspring. This phenomenon is called doubly uniparental inheritance (DUI). Under DUI, it has been reported that female mussels generally have only maternal mtDNA (F type). In this study, we examined the mode of mtDNA transmission in Mytilus galloprovincialis using M and F type-specific primer sets. The ratio of M and F types were measured in each sample by SNaPshot. The M type was detected in the adductor muscle and female gonad of all females. In unfertilized eggs spawned by 84.6% of females (22/26), M type was also detected. The F type was more abundant than the M type in all females. Although the ratio of M type in females was very low, all females contained the M type. From these results, we propose a new possibility about DUI inheritance. The presence of M type in unfertilized eggs indicates that the M type of eggs may also contribute to M type inheritance.     

Interspecies transfer of female mitochondrial DNA is coupled with role-reversals and departure from neutrality in the mussel Mytilus trossulus.

Mussels of the genus Mytilus have distinct and highly diverged male and female mitochondrial DNA (mtDNA) genomes with separate routes of inheritance. Previous studies of European populations of Mytilus trossulus demonstrated that 33% of males are heteroplasmic for a second mtDNA genome of increased length and that hybridization with Mytilus edulis does not block mtDNA introgression, in contrast to reports for American populations. Here, we demonstrate that the female mtDNA type of M. edulis has replaced the resident female mtDNA type of European M. trossulus. This is supported by COIII sequence data indicating that the female mtDNA of European M. trossulus is very similar to that of M. edulis and that in phylogenetic trees, the mtDNAs of these two species cluster together but separately from American M. trossulus sequences, the latter not being disturbed by introgressive hybridization. We also provide evidence that the mtDNA genome of increased length found in heteroplasmic males of European M. trossulus derives from a recent partition of an introgressed M. edulis female type into the male route of transmission. Neutrality tests reveal that European populations of M. trossulus display an excess of replacement polymorphism within the female mtDNA type with respect to conspecific American populations, as well as a significant excess of rare variants, of a similar magnitude to those previously reported for the invading European M. edulis mtDNA. Results are consistent with a nearly neutral model of molecular evolution and suggest that selection acting on European M. trossulus mtDNA is largely independent of the nuclear genetic background.     

Transition and transversion rate in the evolution of animal mitochondrial DNA

We present a further application of the stochastic model previously described (Lanave et al., 1984, 1985) for measuring the nucleotide substitution rate in the mammalian evolution of the mitochondrial DNA (mtDNA). The applicability of this method depends on the validity of "stationarity conditions" (equal nucleotide frequencies at first, second and third silent codon positions in homologous protein coding genes). In the comparison of homologous sequences satisfying the stationarity condition at the silent sites, only the four codon families (quartets) for which both transitions and transversions are silent at the third position are considered here. This has allowed us to estimate the transition and transversion rates for any pair of species. We have analyzed the third silent codon position of the triplet rat-mouse-cow, of a series of slightly divergent primates and of two Drosophila species. In terms of two external dating input we have then determined the phylogenetic trees for rat, mouse, and cow as well as for a number of primates including man. The phylogenetic tree that we have derived for the triplet rat, mouse and cow agrees with that we had previously determined by analyzing the first, second and third silent codon positions (in both duets and quartets) of mt genes (Lanave et al., 1985). For primates our method leads to the following branching order from the oldest to the most recent: Gibbon, Orangutan, Gorilla, Chimpanzee and Man. In absolute time, fixing the distance Chimpanzee-Man as 5 million years (Myr) we estimate the dating of the divergence nodes as: Gorilla 7 Myr; Orangutan 16 Myr; Gibbon 20 Myr. In all cases analyzed, the transition rate has been found to be substantially higher than the transversion rate. Moreover we have found that the transition/transversion ratio is different in the various lineages. We suggest that this fact is probably related to the nucleotide frequencies at the third silent codon position.     

Molecular evolution and recombination in gender-associated mitochondrial DNAs of the Manila clam Tapes philippinarum

Doubly uniparental inheritance (DUI) provides an intriguing system for addressing aspects of molecular evolution and intermolecular recombination of mitochondrial DNA. For this reason, a large sequence analysis has been performed on Tapes philippinarum (Bivalvia, Veneridae), which has mitochondrial DNA heteroplasmy that is consistent with a DUI. The sequences of a 9.2-kb region (containing 29 genes) from 9 individuals and the sequences of a single gene from another 44 individuals are analyzed. Comparisons suggest that the two sex-related mitochondrial genomes do not experience a neutral pattern of divergence and that selection may act with varying strength on different genes. This pattern of evolution may be related to the long, separate history of M and F genomes within their tissue-specific "arenas." Moreover, our data suggest that recombinants, although occurring in soma, may seldom be transmitted to progeny in T. philippinarum.     

Complete sequences of mitochondrial genomes from the Baltic mussel Mytilus trossulus

Marine mussels Mytilus possess two mitochondrial (mt) genomes, which undergo doubly uniparental inheritance (DUI). Female (F) and male (M) genomes are usually highly diverged at the sequence level. Both genomes contain the same set of metazoan genes (for 12 proteins, 2 rRNAs and 23 tRNAs), both lack the atp8 gene and have two tRNAs for methionine. However, recently recombination between those variants has been reported. Both original F and M mt genomes of M. trossulus were replaced by M. edulis mtDNA in the Baltic populations. Highly diverged M genome occurs rarely in the Baltic mussels. Full sequences of the M genome identified in males (sperm) and F genome in females (eggs) were obtained. Both genomes were diverged by 24% in nucleotide sequence, but had similar nucleotide composition and codon usage bias. Constant domain (CD) of the control region (CR), the tRNA and rRNA genes were the most conserved. The most diverged was the variable domain 1 (VD1) of the control region. The F genome was longer than M by 147 bp. and the main difference was localised in the VD1 region. No recombination was observed in whole mtDNA of both studied variants. Nuclear mitochondrial pseudogenes (numts) have not been found by hybridisation with probes complementary to several fragments of the Baltic M. trossulus mtDNA.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

The marine mites Hyadesia sp. and Copidognathus sp. Associated with the mussel Mytilus galloprovincialis

Two species of marine mites belonging to the families Hyadesiidae and Halacaridae, Hyadesia sp. and Copidognathus sp., respectively, were found associated with the mussel Mytilus galloprovincialis from Baja California in NW México. The first species was found inside the mussel gut with an intensity ranging from one to six mites per mussel and their prevalence was from 20.0 to 46.7%; this species was also found living free in the sediment at a density of 0.7 mite/100 ml. The second species was found on the mantle and gills of the host with an intensity ranging from one to three mites per host and their prevalence was from 3.3 to 6.7%; this species was abundant (4.5 mites/100 ml) and living free in the sediment around mussel clumps. Hyadesia sp. was found alive and attached in the gut of the mussel. A histological analysis revealed this species in the lumen of intestine surrounded by mucus and attached to the epithelial cells of the intestine, where some disorder of epithelial cells was associated. Moreover, this mite may be encapsulated by hemocytes inside the digestive diverticulum, the reproductive follicle, or the connective tissue surrounding the diverticulum. No damages to branches or gills resulting from the presence of Copidognathus sp. were observed. The results suggest that these mites are occasional invaders of mussels; however, as a result of this infestation, Hyadesia sp. may produce damage in the host's tissues. This is the first record of marine mites inside the gut, reproductive follicles, branches, and mantle of a marine bivalve.     

Molecular population genetics of the male and female mitochondrial DNA molecules of the California sea mussel, Mytilus californianus

The presence of two gender-associated mitochondrial genomes in marine mussels provides a unique opportunity to investigate the dynamics of mtDNA evolution without complications inherent in interspecific comparisons. Here, we assess the relative importance of selection, mutation, and differential constraint in shaping the patterns of polymorphism within and divergence between the male (M) and female (F) mitochondrial genomes of the California sea mussel, Mytilus californianus. Partial sequences were obtained from homologous regions of four genes (nad2, cox1, atp6, and nad5) totaling 2307 bp in length. The M and F mtDNA molecules of M. californianus exhibited extensive levels of nucleotide polymorphism and were more highly diverged than observed in other mytilids (overall Tamura-Nei distances >40%). Consistent with previous studies, the M molecule had significantly higher levels of silent and replacement polymorphism relative to F. Both genomes possessed large numbers of singleton and low-frequency mutations that gave rise to significantly negative Tajima's D values. Mutation-rate scalars estimated for silent and replacement mutations were elevated in the M genome but were not sufficient to account for its higher level of polymorphism. McDonald-Kreitman tests were highly significant at all loci due to excess numbers of fixed replacement mutations between molecules. Strong purifying selection was evident in both genomes in keeping the majority of replacement mutations at low population frequencies but appeared to be slightly relaxed in M. Our results suggest that a reduction in selective constraint acting on the M genome remains the best explanation for its greater levels of polymorphism and faster rate of evolution.