Nonneutral evolution and differential mutation rate of gender-associated mitochondrial DNA lineages in the marine mussel Mytilus.

Mussels have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally, and the F type is transmitted maternally. To test hypotheses of the molecular evolution of both mtDNA genomes, 50 nucleotide sequences were obtained for 396 bp of the COIII gene of European populations of Mytilus edulis and the Atlantic and Mediterranean forms of M. galloprovincialis. Analysis based on the proportion of synonymous and nonsynonymous substitutions indicate that mtDNA is evolving in a non-neutral and complex fashion. Previous studies on American mussels demonstrated that the F genome experiences a higher purifying selection and that the M genome evolves faster. Here we show that these patterns also hold in European populations. However, in contrast to American populations, where an excess of replacement substitution between F and M lineages has been reported, a significant excess of replacement polymorphism within mtDNA lineages is observed in European populations of M. galloprovincialis. European populations also show an excess of replacement polymorphism within the F but not within the M genome with respect to American M. trossulus, as well as a consistent pattern of excess of rare variants in both F and M genomes. These results are consistent with a nearly neutral model of molecular evolution and a recent relaxation of selective constraints on European mtDNA. Levels of diversity are significantly higher for the M than F genome, and the M genome also accumulates synonymous and nonsynonymous substitutions at a higher rate, in contrast with earlier reports where no difference for the synonymous rate was observed. It is suggested that a subtle balance between relaxed selection and a higher mutation rate explains the faster evolutionary rate of the M lineage.     

Mitochondrial DNA variation in a species with two mitochondrial genomes: the case of Mytilus galloprovincialis from the Atlantic, the Mediterranean and the Black Sea

We have examined mitochondrial DNA (mtDNA) variation in samples of the mussel Mytilus galloprovincialis from the Black Sea, the Mediterranean and the Spanish Atlantic coast by scoring for presence or absence of cleavage at 20 restriction sites of a fragment of the COIII gene and at four restriction sites of the 16S RNA gene. This species contains two types of mtDNA genomes, one that is transmitted maternally (the F type) and one that is transmitted paternally (the M type). The M genome evolves at a higher rate than the F genome. Normally, females are homoplasmic for an F type and males are heteroplasmic for an F and an M type. Occasionally molecules from the F lineage invade the paternal transmission route, resulting in males that carry two F-type mtDNA genomes. These features of the mussel mtDNA system give rise to a new set of questions when using mtDNA variation in population studies and phylogeny. We show here that the two mtDNA types provide different information with regard to amounts of variation and genetic distances among populations. The F genome exhibits higher degrees of diversity within populations, while the M genome produces higher degrees of differentiation among populations. There is a strong differentiation between the Atlantic and the Black Sea. The Mediterranean samples have intermediate haplotype frequencies, yet are much closer to the Black Sea than to the Atlantic. We conclude that in this species gene flow among the three Seas is restricted and not enough to erase the combined effect of mutation and random drift. In one sample, that from the Black Sea, the majority of males did not contain an M mtDNA type. This suggests that a molecule of the maternal lineage has recently invaded the paternal route and has increased its frequency in the population to the point that the present pool of paternally transmitted mtDNA molecules is highly heterogeneous and cannot be used to read the population's history. This liability of the paternal route means that in species with doubly uniparental inheritance, the maternal lineage provides more reliable information for population and phylogenetic studies.     

Interspecific differentiation and intraspecific substructure in two closely related clupeids with extensive hybridization, Alosa alosa and Alosa fallax

In this study, the evolutionary relationships within and among populations of European shads, Alosa alosa and Alosa fallax , was investigated. Screening of allelic variation across eight allozyme loci and sequencing 448 bp of the mtDNA cytochrome b gene in 14 rivers throughout the range of the species supported that the two taxa were independent lineages (1·3% net nucleotide divergence) despite extensive hybridization. Genetic diversity and structure was considerably higher in A. fallax than A. alosa and the former species revealed evidence of distinct lineages in the Mediterranean and Atlantic basins. A Bayesian clustering approach combined with gill raker counts verified that individuals of the two species could be assigned to their parent group with relatively high confidence. Evaluation of hybridization in the Lima and Mondego Rivers in Portugal provided evidence that introgression is extensive but is not currently obscuring (through hybrid swarming) the diagnosability of the two species.     

Pleistocene separation of mitochondrial lineages of Mytilus spp. mussels from Northern and Southern Hemispheres and strong genetic differentiation among southern populations

Smooth-shelled mussels, Mytilus spp., have an antitropical distribution. In the Northern Hemisphere, the M. edulis complex of species is composed of three genetically well delineated taxa: M. edulis, M. galloprovincialis and M. trossulus. In the Southern Hemisphere, morphological characters, allozymes and intron length polymorphisms suggest that Mytilus spp. populations from South America and Kerguelen Islands are related to M. edulis and those from Australasia to M. galloprovincialis. On the other hand, a phylogeny of the 16S rDNA mitochondrial locus demonstrates a clear distinctiveness of southern mussels and suggests that they are related to Mediterranean M. galloprovincialis. Here, we analysed the faster-evolving cytochrome oxidase subunit I locus. The divergence between haplotypes of populations from the two hemispheres was confirmed and was found to predate the divergence between haplotypes of northern M. edulis and M. galloprovincialis. In addition, strong genetic structure was detected among the southern samples, revealing three genetic entities that correspond to (1) South America and Kerguelen Island, (2) Tasmania, (3) New Zealand. Using the trans-Arctic interchange as a molecular clock calibration, we estimated the time since divergence of populations from the two hemispheres to be between 0.5 million years (MY) and 1.3 MY (average 0.84 MY). The contrasting patterns observed for the nuclear and the organelle genomes suggested two alternative, complex scenarios: two trans-equatorial migrations and the existence of differential barriers to mitochondrial and nuclear gene flow, or a single trans-equatorial migration and a view of the composition of the nuclear genome biased by taxonomic preconception.     

The Effects of Natural Hybridization on the Regulation of Doubly Uniparental Mtdna Inheritance in Blue Mussels (Mytilus Spp.)

Blue mussels in the Mytilus edulis species complex have a doubly uniparental mode of mtDNA inheritance with separate maternal and paternal mtDNA lineages. Female mussels inherit their mtDNA solely from their mother, while males inherit mtDNA from both parents. In the male gonad the paternal mtDNA is preferentially replicated so that only paternal mtDNA is transmitted from fathers to sons. Hybridization is common among differentiated blue mussel taxa; whenever it involves M. trossulus, doubly uniparental mtDNA inheritance is disrupted. We have found high frequencies of males without and females with paternal mtDNA among hybrid mussels produced by interspecific matings between M. galloprovincialis and M. trossulus. In contrast, hybridization between M. galloprovincialis and M. edulis does not affect doubly uniparental inheritance, indicating a difference in the divergence of the mechanisms regulating mtDNA inheritance among the three blue mussel taxa. Our data indicate a high frequency of disrupted mtDNA transmission in F(1) hybrids and suggest that two separate mechanisms, one regulating the transmission of paternal mtDNA to males and another inhibiting the establishment of paternal mtDNA in females, act to regulate doubly uniparental inheritance. We propose a model for the regulation of doubly uniparental inheritance that is consistent with these observations.     

Adaptive gamete-recognition divergence in a hybridizing Mytilus population

Gamete-recognition proteins often evolve rapidly, but it is not known if their divergence occurs within species and corresponds with the evolution of reproductive isolation, or if divergence typically accumulates between already isolated lineages. We examined the evolution of a candidate gamete-recognition protein in several sympatric and allopatric populations of Mytilus blue mussels, species that hybridize in nature. Within a single species, Mytilus galloprovincialis, we found adaptive divergence of Lysin-M7, a sperm acrosomal protein that dissolves the egg vitelline envelope during fertilization. Mytilus galloprovincialis Lysin-M7 alleles group into two distinct clades (termed G and G(D)), and individual alleles in these clades are separated from each other by at least three and up to eleven amino-acid substitutions. Maximum-likelihood estimates of selective pressure (dN/dS =omega) implicate selection in the divergence between M. galloprovincialis Lysin-M7 clades, and within the G(D) clade. Exact tests of population differentiation indicate that the relative frequency of G and G(D) Lysin-M7 alleles differs significantly among M. galloprovincialis populations. Compared with allopatric Mediterranean samples, Lysin-M7 alleles in the G(D) clade are found at elevated frequency in samples from the East Atlantic and California, areas of secondary contact and hybridization between Mytilus species, and Australia, an area of unknown species composition. Adaptive divergence between the alleles most common in allopatry and those found at elevated frequency in samples from sympatry suggests that selection pressures acting in hybridizing populations, likely following Pleistocene secondary contact with M. edulis in the East Atlantic, drove the divergence of Lysin-M7 in M. galloprovincialis.     

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.     

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.     

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.     

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long‐standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well‐studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best‐supported scenario was a long period of allopatric isolation during the first three‐quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.     

Evidence that the large noncoding sequence is the main control region of maternally and paternally transmitted mitochondrial genomes of the marine mussel (Mytilus spp.)

Both the maternal (F-type) and paternal (M-type) mitochondrial genomes of the Mytilus species complex M. edulis/galloprovincialis contain a noncoding sequence between the l-rRNA and the tRNA(Tyr) genes, here called the large unassigned region (LUR). The LUR, which is shorter in M genomes, is capable of forming secondary structures and contains motifs of significant sequence similarity with elements known to have specific functions in the sea urchin and the mammalian control region. Such features are not present in other noncoding regions of the F or M Mytilus mtDNA. The LUR can be divided on the basis of indels and nucleotide variation in three domains, which is reminiscent of the tripartite structure of the mammalian control region. These features suggest that the LUR is the main control region of the Mytilus mitochondrial genome. The middle domain has diverged by only 1.5% between F and M genomes, while the average divergence over the whole molecule is approximately 20%. In contrast, the first domain is among the most divergent parts of the genome. This suggests that different parts of the LUR are under different selection constraints that are also different from those acting on the coding parts of the molecule.     

Multilocus allozyme differentiation of populations of the mussel Mytilus galloprovincialis Lmk. from Moroccan coasts

Along the Moroccan coasts, the systematic status of Mytilus populations have been for a long time uncertain and confused, due to the use of unreliable morphometric criteria. In the present study, allozyme markers reveal the exclusive existence of M. galloprovincialis on Mediterranean and Atlantic coasts. Nei's genetic distances are low and reflect a high gene flow between Atlantic and Mediterranean populations. However, a significant multilocus discontinuity revealed by F-statistics separate southern Atlantic populations from Mediterranean and north Atlantic ones and could be explained by a gene flow breaking because of a larval dispersal decrease, due to a sea surface current direction change from Cap Ghir towards the Canaries archipelago, and probably by differential selection effects in these two geographic areas.     

Evidence for recombination of mtDNA in the marine mussel Mytilus trossulus from the Baltic

A number of studies have claimed that recombination occurs in animal mtDNA, although this evidence is controversial. Ladoukakis and Zouros (2001) provided strong evidence for mtDNA recombination in the COIII gene in gonadal tissue in the marine mussel Mytilus galloprovincialis from the Black Sea. The recombinant molecules they reported had not however become established in the population from which experimental animals were sampled. In the present study, we provide further evidence of the generality of mtDNA recombination in Mytilus by reporting recombinant mtDNA molecules in a related mussel species, Mytilus trossulus, from the Baltic. The mtDNA region studied begins in the 16S rRNA gene and terminates in the cytochrome b gene and includes a major noncoding region that may be analogous to the D-loop region observed in other animals. Many bivalve species, including some Mytilus species, are unusual in that they have two mtDNA genomes, one of which is inherited maternally (F genome) the other inherited paternally (M genome). Two recombinant variants reported in the present study have population frequencies of 5% and 36% and appear to be mosaic for F-like and M-like sequences. However, both variants have the noncoding region from the M genome, and both are transmitted to sperm like the M genome. We speculate that acquisition of the noncoding region by the recombinant molecules has conferred a paternal role on mtDNA genomes that otherwise resemble the F genome in sequence.     

Comparative Genomics of Marine Mussels (Mytilus spp.) Gender Associated mtDNA: Rapidly Evolving atp8

The unusual mode of mitochondrial DNA inheritance, with two separate: maternal (F) and paternal (M) lineages, gives unique opportunities to study the evolution of the mitochondrial genome. This system was first discovered in the marine mussels Mytilus. The three related species: Mytilus edulis, Mytilus galloprovincialis and Mytilus trossulus form a complex in which the divergence of M and F lineages pre-dates the speciation. The complete mitochondrial genomes of both lineages were known for all species except Pacific M. trossulus. Here we report, for the first time, the complete sequences of both mitochondrial genomes of Pacific M. trossulus, filling the gap. While the reported M and F genomes are highly diverged (26%), they have similar organisation. The only difference is the translocation of one tRNA gene into the long, mosaic control region of the F genome. Consistent presence of an ORF which most likely represents the atp8 gene was confirmed in both genomes. The predicted protein has characteristics expected of the functional atp8 even though the M and F versions are markedly different in length. Comparative analysis involving all three species led to the conclusion that the cause of a faster evolution of atp8 and Mytilus mtDNA in general is most likely the Compensation-Draft Feedback process coupled with relatively relaxed selection in the M lineage. Thus, we postulate that the adaptive changes may have played a role in the emergence of highly diverged, barely recognizable atp8 in Mytilus mussels.     

Electrophoretic investigation of systematic relationships in the marine mussels Modiolus modiolus L., Mytilus edulis L., and Mytilus galloprovincialis Lmk. (Mytilidae; Mollusca)

Genetic variation was assayed electrophoretically at 13–16 loci in Modiolus modiolus, Mytilus edulis, and Mytilus galloprovincialis. High genetic distance ( D ) values were observed between Modiolus modiolus and Mytilus edulis (1.516 ± 0.523) and between Modiolus modiolus and Mytilus galloprovincialis (1.564 ± 0.539), whereas the distance between Mytilus edulis and M. galloprovincialis (0.167 ± 0.118) was rather low. The systematic status ot Mytilus edulis and M. galloprovincialis is discussed in relation to these lindings and the genetic distance values are used to estimate divergence times which in turn are compared with paleontological estimates. The observations of high average heterozygosity in Modiolus modiolus, and high correlations of locus heterozygosity between taxa are discussed briefly.     

Quantitation of the male and female types of mitochondrial DNA in a blue mussel, Mytilus galloprovincialis, using real-time polymerase chain reaction assay

The system termed doubly uniparental inheritance (DUI) of mitochondrial transmission to progeny has been reported in Mytilus. Under DUI, it has been thought that males have both paternally (M type) and maternally (F type) transmitted mitochondrial DNA (mtDNA), and females have only F type. However, the presence of M type in females has been reported. To clarify the ratio of M type to F type mtDNA in female and male tissues to further our understanding of mitochondrial transmission, we developed a procedure to measure the copy numbers of the two types of mtDNA in Mytilus galloprovincialis using a real-time polymerase chain reaction assay. The following results were obtained by this method. In females, the copy numbers of M type mtDNA detected in adductor muscle, gonad and eggs were approximately 10 000-fold lower than those of F type. In males, F type dominated in adductor muscle, as in the female tissue. However, copy numbers of M type mtDNA were approximately 1000-fold higher than those of F type in gonad and 100 000-fold higher than those of F type in sperm. We examined the quantity relationship between the two types of mtDNA and the transmission mechanism of mtDNA in M. galloprovincialis.     

An RFLP assay to determine if Mytilus galloprovincialis Lmk. (Mytilidae; Bivalvia) is of Northern or Southern hemisphere origin

Mytilus galloprovincialis is one of three smooth shelled blue mussel species belonging to the Mytilus edulis species complex. Naturally occurring and introduced populations of M. galloprovincialis are widely distributed throughout many regions of the globe. Mytilus galloprovincialis includes morphologically indistinguishable Northern and Southern hemisphere mtDNA lineages that have been separated for ～1 my. To distinguish recently introduced Northern M. galloprovincialis from resident Southern M. galloprovincialis in New Zealand, we developed a 16s rRNA RFLP assay. We compared RFLP assignments of 178 mussels with those generated from a 16s rRNA sequence-estimated phylogeny. All mussels were correctly assigned by the RFLP to their sequence-based phylogenetic placement. This assay allows the rapid identification of Northern and Southern hemisphere M. galloprovincialis and will provide an important tool for monitoring human mediated introductions of otherwise cryptic lineages.     

Macoma balthica in the White and Barents Seas: properties of a widespread marine hybrid swarm (Mollusca: Bivalvia)

A main molecular subdivision in the circumpolar Macoma balthica complex has been described between Atlantic and Pacific taxa. In NE Europe, the clams of the White and Barents Seas, however, show deviant genetic structures. Using allozyme and mitochondrial DNA data, we explore the hypothesis that these deviations result from hybridization between an Atlantic (M. b. rubra) and an invading Pacific (M. b. balthica) lineage. A practically pure Atlantic Macoma extends from France north to the Varanger Peninsula (NE Norway), whereas populations farther east have genetic compositions intermediate between true Atlantic and true Pacific. Admixture estimates range from 32 to 90% Pacific contribution, with a notable deviation in a nearly pure Atlantic outpost in the Mezen Bay (NE White Sea). The pattern of variation is not one of a simple collinear mixing however. Different characters exhibit different degrees of introgression, and the relative introgression varies regionally. Yet, there are practically no interlocus genotypic disequilibria between the diverged loci, which brings out the White Sea-Barents Sea M. balthica as the best-documented marine animal hybrid swarms so far, arisen through amalgamation of genomes previously isolated since pre-Pleistocene times. On top of the main admixture pattern, strong geographical structuring is also seen in characters unrelated to the principal systematic distinction. The persistence of the regional patterns indicates restricted gene flow at the present time, despite the high dispersal potential of the species. The causes of this structuring could be in a complex history of colonization events and features of local hydrography enhancing isolation and divergence of populations.     

Genomic reticulation indicates mixed ancestry in Southern-Hemisphere Mytilus spp. mussels

Previous surveys of allozyme variation in smooth-shell Mytilus spp. mussels have reported the presence in the Southern Hemisphere of both Mytilus edulis and Mytilus galloprovincialis mussels. In the present study, nuclear DNA markers mac-1 and Glu-5 '/ Glu-3 ', both diagnostic for Northern-Hemisphere M. edulis and M. galloprovincialis , were used to further characterize the nuclear genomes of M. edulis from Kerguelen and M. galloprovincialis from Tasmania. Genomic reticulation was observed, with typical M. edulis allelomorphs fixed in both populations at locus mac-1 whereas, at locus Glu-5 '/ Glu-3 ', allelomorphs characteristic of M. galloprovincialis were present in Kerguelen and nearly fixed in Tasmania. Kerguelen mussels had a genome of mixed M. edulis and M. galloprovincialis ancestry without evidence of barriers to merging as shown by Hardy–Weinberg and linkage equilibrium. Tasmanian mussels possessed a predominantly M. galloprovincialis genomic background introgressed by M. edulis allelomorphs at locus mac-1 . Genetic drift superimposed on ancient hybridization and introgression may explain the genomic reticulation observed in both Kerguelen and Tasmanian mussels. There was no evidence of a recent introduction of Northern-Hemisphere M. galloprovincialis or M. edulis to Kerguelen or Tasmania.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 747–754.     

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.