Differential patterns of male and female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Trans-Atlantic Distribution and Introgression as Inferred from Single Nucleotide Polymorphism: Mussels Mytilus and Environmental Factors

Large-scale climate changes influence the geographic distribution of biodiversity. Many taxa have been reported to extend or reduce their geographic range, move poleward or displace other species. However, for closely related species that can hybridize in the natural environment, displacement is not the only effect of changes of environmental variables. Another option is subtler, hidden expansion, which can be found using genetic methods only. The marine blue mussels Mytilus are known to change their geographic distribution despite being sessile animals. In addition to natural dissemination at larval phase—enhanced by intentional or accidental introductions and rafting—they can spread through hybridization and introgression with local congeners, which can create mixed populations sustaining in environmental conditions that are marginal for pure taxa. The Mytilus species have a wide distribution in coastal regions of the Northern and Southern Hemisphere. In this study, we investigated the inter-regional genetic differentiation of the Mytilus species complex at 53 locations in the North Atlantic and adjacent Arctic waters and linked this genetic variability to key local environmental drivers. Of seventy-nine candidate single nucleotide polymorphisms (SNPs), all samples were successfully genotyped with a subset of 54 SNPs. There was a clear interregional separation of Mytilus species. However, all three Mytilus species hybridized in the contact area and created hybrid zones with mixed populations. Boosted regression trees (BRT) models showed that inter-regional variability was important in many allele models but did not prevail over variability in local environmental factors. Local environmental variables described over 40% of variability in about 30% of the allele frequencies of Mytilus spp. For the 30% of alleles, variability in their frequencies was only weakly coupled with local environmental conditions. For most studied alleles the linkages between environmental drivers and the genetic variability of Mytilus spp. were random in respect to “coding” and “non-coding” regions. An analysis of the subset of data involving functional genes only showed that two SNPs at Hsp70 and ATPase genes correlated with environmental variables. Total predictive ability of the highest performing models (r² between 0.550 and 0.801) were for alleles that discriminated most effectively M. trossulus from M. edulis and M. galloprovincialis, whereas the best performing allele model (BM101A) did the best at discriminating M. galloprovincialis from M. edulis and M. trossulus. Among the local environmental variables, salinity, water temperature, ice cover and chlorophyll a concentration were by far the greatest predictors, but their predictive performance varied among different allele models. In most cases changes in the allele frequencies along these environmental gradients were abrupt and occurred at a very narrow range of environmental variables. In general, regions of change in allele frequencies for M. trossulus occurred at 8–11 psu, 0–10 °C, 60%–70% of ice cover and 0–2 mg m⁻³ of chlorophyll a, M. edulis at 8–11 and 30–35 psu, 10–14 °C and 60%–70% of ice cover and for M. galloprovincialis at 30–35 psu, 14–20 °C.

Twin introductions by independent invader mussel lineages are both associated with recent admixture with a native congener in Australia

Introduced species can impose profound impacts on the evolution of receiving communities with which they interact. If native and introduced taxa remain reproductively semi-isolated, human-mediated secondary contact may promote genetic exchange across newly created hybrid zones, potentially impacting native genetic diversity and invasive species spread. Here, we investigate the contributions of recent divergence histories and ongoing (post-introduction) gene flow between the invasive marine mussel, Mytilus galloprovincialis, and a morphologically indistinguishable and taxonomically contentious native Australian taxon, Mytilus planulatus. Using transcriptome-wide markers, we demonstrate that two contemporary M. galloprovincialis introductions into south-eastern Australia originate from genetically divergent lineages from its native range in the Mediterranean Sea and Atlantic Europe, where both introductions have led to repeated instances of admixture between introduced and endemic populations. Through increased genome-wide resolution of species relationships, combined with demographic modelling, we validate that mussels sampled in Tasmania are representative of the endemic Australian taxon (M. planulatus), but share strong genetic affinities to M. galloprovincialis. Demographic inferences indicate late-Pleistocene divergence times and historical gene flow between the Tasmanian endemic lineage and northern M. galloprovincialis, suggesting that native and introduced taxa have experienced a period of historical isolation of at least 100,000 years. Our results demonstrate that many genomic loci and sufficient sampling of closely related lineages in both sympatric (e.g. Australian populations) and allopatric (e.g. northern hemisphere Mytilus taxa) ranges are necessary to accurately (a) interpret patterns of intraspecific differentiation and to (b) distinguish contemporary invasive introgression from signatures left by recent divergence histories in high dispersal marine species. More broadly, our study fills a significant gap in systematic knowledge of native Australian biodiversity and sheds light on the intrinsic challenges for invasive species research when native and introduced species boundaries are not well defined.

Local cold adaption increases the thermal window of temperate mussels in the Arctic

Species expand towards higher latitudes in response to climate warming, but the pace of this expansion is related to the physiological capacity to resist cold stress. However, few studies exist that have quantified the level of inter-population local adaptation in marine species freeze tolerance, especially in the Arctic. We investigated the importance of cold adaptation and thermal window width towards high latitudes from the temperate to the Arctic region. We measured upper and lower lethal air temperatures (i.e. LT and LT50) in temperate and Arctic populations of blue mussels (Mytilus edulis), and analysed weather data and membrane fatty acid compositions, following emersion simulations. Both populations had similar upper LT (~38 °C), but Arctic mussels survived 4°C colder air temperatures than temperate mussels (-13 vs. -9°C, respectively), corresponding to an 8% increase in their thermal window. There were strong latitudinal relationships between thermal window width and local air temperatures, indicating Arctic mussels are highly adapted to the Arctic environment where the seasonal temperature span exceeds 60°C. Local adaptation and local habitat heterogeneity thus allow leading-edge M. edulis to inhabit high Arctic intertidal zones. This intraspecific pattern provides insight into the importance of accounting for cold adaptation in climate change, conservation and biogeographic studies.

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.

Genetic diversity and connectivity within Mytilus spp. in the subarctic and Arctic

Climate changes in the Arctic are predicted to alter distributions of marine species. However, such changes are difficult to quantify because information on present species distribution and the genetic variation within species is lacking or poorly examined. Blue mussels, Mytilus spp., are ecosystem engineers in the coastal zone globally. To improve knowledge of distribution and genetic structure of the Mytilus edulis complex in the Arctic, we analyzed 81 SNPs in 534 Mytilus spp. individuals sampled at 13 sites to provide baseline data for distribution and genetic variation of Mytilus mussels in the European Arctic. Mytilus edulis was the most abundant species found with a clear genetic split between populations in Greenland and the Eastern Atlantic. Surprisingly, analyses revealed the presence of Mytilus trossulus in high Arctic NW Greenland (77°N) and Mytilus galloprovincialis or their hybrids in SW Greenland, Svalbard, and the Pechora Sea. Furthermore, a high degree of hybridization and introgression between species was observed. Our study highlights the importance of distinguishing between congener species, which can display local adaptation and suggests that information on dispersal routes and barriers is essential for accurate predictions of regional susceptibility to range expansions or invasions of boreal species in the Arctic.

O mexilhão Perna perna no Brasil: nativo ou exótico?

A mitilicultura no Brasil se baseia no mexilhão Perna perna (Linnaeus, 1758), que é encontrado em todo o litoral brasileiro, sendo especialmente abundante do Espírito Santo a Santa Catarina. Nos últimos anos, foi lançada a hipótese de que P. perna seja uma espécie exótica no litoral brasileiro. A hipótese se baseia na análise da malacofauna presente em sítios arqueológicos. Todas as contestações levantadas não remetem à uma conclusão, pois precisam de estudos específicos, com metodologias claras, aliando arqueologia, ecologia e biologia molecular. O objetivo deste trabalho foi estudar a condição de nativo ou exótico do mexilhão Perna perna no Brasil, a partir de levantamento dos resultados em sítios arqueológicos, de técnicas moleculares e de datação de conchas com C14. A datação indicou que as amostras de Perna perna do sítio arqueológico Rio do Meio/Jurerê, Florianópolis/SC, têm idade de 720±30 e 780±30 anos. O cálculo do tempo de divergência indicou que a separação das populações brasileiras e africanas ocorreu por volta de 200 mil anos. Os resultados apontam a presença da espécie no território brasileiro muito antes do descobrimento do Brasil pelos portugueses no ano de 1500, indicando que P. perna é de fato uma espécie nativa.

Complete male mitochondrial genomes of European Mytilus edulis mussels

Doubly uniparental inheritance (DUI) results in the existence of two gender-specific, divergent mtDNA lineages within a single species. Under DUI, the female genome (F) is transmitted from mothers to the whole offspring, and the male genome (M) is transmitted exclusively from fathers to sons. This system was first described in a marine mussels Mytilus edulis inhabiting European coastal waters, over a decade ago. Despite that, the complete sequence of the M genome from the European M. edulis mussels remained unknown. Here we announce it for the first time. The announcement is based on the two haplotypes isolated from heteroplasmic males of European M. edulis sampled at two moderately distant locations: southern North Sea and western Baltic. The two M genomes are quite similar both in length (16,631 and 16,632 bp) and in sequence (98.3%). Furthermore, both newly sequenced genomes are closely related to the genomes described from Baltic M. trossulus.

Simple identification of mitochondrial lineages in contact zones based on lineage-selective primers

A variety of research projects focus on genetic variation among and within maternal lineages as encompassed by mitochondrial DNA (mtDNA). While mtDNA often differs substantially between species, large differences may also be found within species. The evaluation of such divergent lineages, for example in intraspecific contact zones (hybrid zones), commonly involves sequencing numerous individuals. Large-scale sequencing is both expensive and labour-intensive. Based on sequences from 15 individuals, we devised a simple and quick polymerase chain reaction assay for identification of divergent mtDNA lineages in a secondary contact zone of the side-blotched lizard (Uta stansburiana). The application uses lineage-selective primers to amplify a lineage-diagnostic product, and is based on each group of mtDNA haplotypes being a monophyletic assemblage of haplotypes sharing the same maternal ancestry, deeply divergent from the other group. The assay was tested on a larger sample (n = 147) of specimens from the contact zone, confirming its usefulness in quick and reliable identification of mtDNA lineages. This approach can be modified for other species, provided diagnostic lineage variation is available, and may also be performed in simple laboratory settings while conducting fieldwork.

Molecular population genetics of male and female mitochondrial genomes in subarctic Mytilus trossulus

The doubly uniparental inheritance system allows for the use of two independent mitochondrial genomes for population history investigations. Under this system, two lineages of mitochondrial DNA (mtDNA) exist and males are typically heteroplasmic, having the additional, usually divergent, mitochondrial genome inherited from their male parent. This additional mtDNA typically evolves faster, potentially allowing for insight into more recent events in population history. Few studies did explore this possibility in marine mussels Mytilus showing its usefulness. Recent observations of the Mytilus trossulus mussels who have retained their native mtDNA in European waters posed the question of their origin. Are they part of a population present, but previously undetected, or is this a potentially human mediated, ongoing spread of an invasive species? To tackle this question, we amplified with species-specific primers and sequenced an approximately 1,200-bp-long fragment spanning COIII and ND2 genes from both mitochondrial genomes of mussels sampled at five locations worldwide, covering the whole M. trossulus range. The overall pattern of polymorphisms is compatible with the entirely postglacial history of the whole species, indicating a very deep bottleneck at last glacial maximum, with possible retention of the whole species in a single refugium, and the effective population size of no more than a few thousands. Both analyses of molecular variance and isolation with migration (IM) models point at the West Atlantic as the source of the European M. trossulus mussels, at least the ones who retained their native mtDNA. The hypothesis that this is an ongoing, human-mediated process was considered. To this end, comparison with the well-known case: the introduction of congeneric mussel, Mytilus galloprovincialis, from Mediterranean Sea to Asia was used. This introduction occurred within the last 100 years. The results inferred by the IM model suggest that the timing and structure of transatlantic migration of M. trossulus differs significantly from the M. galloprovincialis case: it is more than 1,000 years old and involves a much larger fraction of the ancestral population. Therefore, most likely, this invasion is not a human-mediated process.

Mitochondrial DNA paradox: sex-specific genetic structure in a marine mussel--despite maternal inheritance and passive dispersal

Background

When genetic structure is identified using mitochondrial DNA (mtDNA), but no structure is identified using biparentally-inherited nuclear DNA, the discordance is often attributed to differences in dispersal potential between the sexes.

Results

We sampled the intertidal rocky shore mussel Perna perna in a South African bay and along the nearby open coast, and sequenced maternally-inherited mtDNA (there is no evidence for paternally-inherited mtDNA in this species) and a biparentally-inherited marker. By treating males and females as different populations, we identified significant genetic structure on the basis of mtDNA data in the females only.

Conclusions

This is the first study to report sex-specific differences in genetic structure based on matrilineally-inherited mtDNA in a passively dispersing species that lacks social structure or sexual dimorphism. The observed pattern most likely stems from females being more vulnerable to selection in habitats from which they did not originate, which also manifests itself in a male-biased sex ratio. Our results have three important implications for the interpretation of population genetic data. First, even when mtDNA is inherited exclusively in the female line, it also contains information about males. For that reason, using it to identify sex-specific differences in genetic structure by contrasting it with biparentally-inherited markers is problematic. Second, the fact that sex-specific differences were found in a passively dispersing species in which sex-biased dispersal is unlikely highlights the fact that significant genetic structure is not necessarily a function of low dispersal potential or physical barriers. Third, even though mtDNA is typically used to study historical demographic processes, it also contains information about contemporary processes. Higher survival rates of males in non-native habitats can erase the genetic structure present in their mothers within a single generation.

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.

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.

Reticulate evolution and marine organisms: the final frontier?

The role that reticulate evolution (i.e., via lateral transfer, viral recombination and/or introgressive hybridization) has played in the origin and adaptation of individual taxa and even entire clades continues to be tested for all domains of life. Though falsified for some groups, the hypothesis of divergence in the face of gene flow is becoming accepted as a major facilitator of evolutionary change for many microorganisms, plants and animals. Yet, the effect of reticulate evolutionary change in certain assemblages has been doubted, either due to an actual dearth of genetic exchange among the lineages belonging to these clades or because of a lack of appropriate data to test alternative hypotheses. Marine organisms represent such an assemblage. In the past half-century, some evolutionary biologists interested in the origin and trajectory of marine organisms, particularly animals, have posited that horizontal transfer, introgression and hybrid speciation have been rare. In this review, we provide examples of such genetic exchange that have come to light largely as a result of analyses of molecular markers. Comparisons among these markers and between these loci and morphological characters have provided numerous examples of marine microorganisms, plants and animals that possess the signature of mosaic genomes.

The flow of antimicrobial peptide genes through a genetic barrier between Mytilus edulis and M. galloprovincialis

We studied the population genetics of two antimicrobial peptide (AMP) loci, called Mytilin B and Mytilus galloprovincialis defensin 2 (MGD2), in the secondary contact mosaic hybrid zone between Mytilus edulis and M. galloprovincialis. The isolation period between the two species was estimated to be approximately 1 million years (range, 0.5 million to 2 million years) long. During this period, coevolution between microbes and the immune system has likely occurred. The secondary contact, which would date back to approximately 25,000 (0-200,000) years, recently allowed these coadaptations to be rearranged through hybridization. Distinctive polymorphisms were uncovered in coding sequences of the two AMP loci such as insertion/deletion of codons or bisubstituted codons. Very low levels of differentiation were observed between populations of the two species at both loci, while other nuclear loci often showed marked structure among the same samples. The absence of population differentiation proved to be the consequence of secondary introgression of highly divergent alleles. While only a few recombinants were observed at the Mytilin B locus, the MGD2 locus showed a high intragenic recombination rate, which increased in the exon coding for the mature peptide. In addition, standard neutrality tests revealed significant deviations from the mutation-drift equilibrium at both loci. These results suggest that either balancing or directional selection is likely to play a role in the evolution of the two AMPs and introgression would be adaptive. However, evidence accumulated at the Mytilin B locus allows neither for identification of the direction of selection nor for any conclusions on whether selection acted directly on the antimicrobial peptide itself. At the MGD2 locus, a spatial variation of polymorphism patterns along the sequence suggests that selection was direct, although the precise nature of the selection (directional vs. balancing) remains unclear. This study concurs with previous reports of an effect of slight selection on AMP genes evolution in other invertebrates, although selection does not necessarily act on the mature peptides.

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.

Ice-age survival of Atlantic cod: agreement between palaeoecology models and genetics

Scant scientific attention has been given to the abundance and distribution of marine biota in the face of the lower sea level, and steeper latitudinal gradient in climate, during the ice-age conditions that have dominated the past million years. Here we examine the glacial persistence of Atlantic cod (Gadus morhua) populations using two ecological-niche-models (ENM) and the first broad synthesis of multi-locus gene sequence data for this species. One ENM uses a maximum entropy approach (<span style='font-size:xx-small;'>Maxent</span> ); the other is a new ENM for Atlantic cod, using ecophysiological parameters based on observed reproductive events rather than adult distribution. Both the ENMs were tested for present-day conditions, then used to hindcast ranges at the last glacial maximum (LGM) ca 21kyr ago, employing climate model data. Although the LGM range of Atlantic cod was much smaller, and fragmented, both the ENMs agreed that populations should have been able to persist in suitable habitat on both sides of the Atlantic. The genetic results showed a degree of trans-Atlantic divergence consistent with genealogically continuous populations on both sides of the North Atlantic since long before the LGM, confirming the ENM results. In contrast, both the ENMs and the genetic data suggest that the Greenland G. morhua population post-dates the LGM.

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.

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.

The unusual system of doubly uniparental inheritance of mtDNA: isn't one enough?

Mitochondria possess their own genetic material (mitochondrial DNA or mtDNA), whose gene products are involved in mitochondrial respiration and oxidative phosphorylation, transcription, and translation. In animals, mitochondrial DNA is typically transmitted to offspring by the mother alone. The discovery of 'doubly uniparental inheritance' (DUI) of mtDNA in some bivalves has challenged the paradigm of strict maternal inheritance (SMI). In this review, we survey recent advances in our understanding of DUI, which is a peculiar system of cytoplasmic DNA inheritance that involves distinct maternal and paternal routes of mtDNA transmission, a novel extension of a mitochondrial gene (cox2), recombination, and periodic 'role-reversals' of the normally male and female-transmitted mitochondrial genomes. DUI provides a unique opportunity for studying nuclear-cytoplasmic genome interactions and the evolutionary significance of different modes of mitochondrial inheritance.

Phylogeography and environmental diversification of a highly adaptable marine amphipod, Gammarus duebeni

Genetic diversity and phylogeographic population structure in the gammarid amphipod, Gammarus duebeni, were investigated across its broad latitudinal distribution in the NE and NW Atlantic by analysis of mitochondrial DNA sequence. Gammarus duebeni has exceptional tolerance of salinity change and inhabits environments ranging from marine to freshwater. The longstanding debate on whether there are distinct marine and freshwater subspecies was assessed by sampling populations from sites characterized by different salinities. Our sequence data demonstrates that there are two major lineages, with little internal geographic structuring. Evidence is provided to suggest a pre-glacial divergence of these two clades, involving segregation between a region historically associated with the freshwater form and the majority of the marine localities on both sides of the Atlantic. A modern contact zone between the marine and freshwater forms is proposed in western Britain.

Two events are responsible for an insertion in a paternally inherited mitochondrial genome of the mussel Mytilus galloprovincialis

Frequent nonhomologous recombination has been previously postulated to explain the 1045-bp insertion in one mitochondrial sperm-transmitted haplotype of Mytilus galloprovincialis. Such recombination would lead to the disruption of gene order and so the existence of a specific mechanism for maintaining the same gene order in both mitochondrial genomes of Mytilus has been proposed. Here the simpler explanation of the observed structure, involving a tandem duplication and a deletion, is presented. Their occasional occurrence in Mytilus mtDNA proves the similarity, not the difference, between animals with and without DUI.

Nearshore fish (Pholis gunnellus) persists across the North Atlantic through multiple glacial episodes

The intertidal biota of the North Atlantic is characterized by two disjunct communities (North American and European) exposed to different climatic regimes during the Pleistocene and in the Holocene. We collect multilocus DNA sequence data from the nearshore fish Pholis gunnellus to help uncover processes determining biogeographical persistence during periodic coastal glaciations. Coalescent-based estimates from the multilocus DNA sequence data suggest that P. gunnellus persisted on both sides of the North Atlantic throughout the last two glacial maxima (> 202,000 years) with little trans-Atlantic gene flow since divergence, very little structure among populations within Europe (Phi(ST) < 0.05) and some structure within the North American coastline (Phi(ST) = 0.0-0.21). Although the ecological flexibility and high local migration of P. gunnellus could have enhanced this species' survival across the Atlantic, logistic regression did not find a significant determinant of trans-Atlantic persistence when considering 12 other North Atlantic phylogeographical studies from the literature.

Divergent Induced Responses to an Invasive Predator in Marine Mussel Populations

Invasive species may precipitate evolutionary change in invaded communities. In southern New England (USA) the invasive Asian shore crab, Hemigrapsus sanguineus, preys on mussels (Mytlius edulis), but the crab has not yet invaded northern New England. We show that southern New England mussels express inducible shell thickening when exposed to waterborne cues from Hemigrapsus, whereas naïve northern mussel populations do not respond. Yet, both populations thicken their shells in response to a long-established crab, Carcinus maenas. Our findings are consistent with the rapid evolution of an inducible morphological response to Hemigrapsus within 15 years of its introduction.

Geographic Variation and Positive Selection on M7 Lysin, an Acrosomal Sperm Protein in Mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely, reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, Mytilus galloprovincialis, and Mytilus trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement (RCD), which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of RCD; neither are there unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection, and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.

From Europe to America: pliocene to recent trans-atlantic expansion of cold-water north atlantic molluscs

Data on the geographical distribution, phylogeny and fossil record of cool-temperate North Atlantic shell-bearing molluscs that live in waters shallower than 100 m depth belong to two biogeographic provinces, one in eastern North America north of Cape Cod, the other in northern Europe. Amphi-Atlantic species, which are found in both provinces, comprise 30.8% of the 402 species in the northeastern Atlantic and 47.3% of the 262 species in the northwestern Atlantic. Some 54.8% of these amphi-Atlantic species have phylogenetic origins in the North Pacific. Comparisons among fossil Atlantic faunas show that amphi-Atlantic distributions became established in the Middle Pliocene (about 3.5 million years ago), and that all represent westward expansions of European taxa to North America. No American taxa spread eastward to Europe without human assistance. These results are in accord with previous phylogeographic studies among populations within several amphi-Atlantic species. Explanations for the unidirectional expansion of species across the Atlantic remain uncertain, but may include smaller size and greater prior extinction of the North American as compared to the European fauna and biased transport mechanisms. Destruction of the European source fauna may jeopardize faunas on both sides of the Atlantic.

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.