Evidence for a slightly deleterious effect of intron polymorphisms at the EF1alpha gene in the deep-sea hydrothermal vent bivalve Bathymodiolus

Divergent ecological histories of two sister Antarctic krill species led to contrasted patterns of genetic diversity in their heat-shock protein (hsp70) arsenal

The Arctic and the Antarctic Peninsula are currently experiencing some of the most rapid rates of ocean warming on the planet. This raises the question of how the initial adaptation to extreme cold temperatures was put in place and whether or not directional selection has led to the loss of genetic variation at key adaptive systems, and thus polar species’ (re)adaptability to higher temperatures. In the Southern Ocean, krill represents the most abundant fauna and is a critical member at the base of the Antarctic food web. To better understand the role of selection in shaping current patterns of polymorphisms, we examined genetic diversity of the cox‐1 and hsp70 genes by comparing two closely related species of Euphausiid that differ in ecology. Results on mtcox‐1 agreed with previous studies, indicating high and similar effective population sizes. However, a coalescent‐based approach on hsp70 genes highlighted the role of positive selection and past demographic changes in their recent evolution. Firstly, some form of balancing selection was acting on the inducible isoform C, which reflected the maintenance of an ancestral adaptive polymorphism in both species. Secondly, E. crystallorophias seems to have lost most of its hsp70 diversity because of a population crash and/or directional selection to cold. Nonsynonymous diversities were always greater in E. superba, suggesting that it might have evolved under more heterogeneous conditions. This can be linked to species’ ecology with E. superba living in more variable pelagic conditions, while E. crystallorophias is strictly associated with continental shelves and sea ice.

Species distribution and population connectivity of deep-sea mussels at hydrocarbon seeps in the Gulf of Mexico

Hydrocarbon seepage is widespread and patchy in the Gulf of Mexico, and six species of symbiont containing bathymodiolin mussels are found on active seeps over wide and overlapping depth and geographic ranges. We use mitochondrial genes to discriminate among the previously known and a newly discovered species and to assess the connectivity among populations of the same species in the northern Gulf of Mexico (GoM). Our results generally validate the morphologically based distribution of the three previously known GoM species of Bathymodiolus, although we found that approximately 10% of the morphologically based identifications were incorrect and this resulted in some inaccuracies with respect to their previously assigned depth and geographical distribution patterns. These data allowed us to confirm that sympatry of two species of Bathymodiolus within a single patch of mussels is common. A new species of bathymodiolin, Bathymodiolus sp. nov., closely related to B. heckerae was also discovered. The two species live at the same depths but have not been found in sympatry and both have small effective population sizes. We found evidence for genetic structure within populations of the three species of Bathymodiolinae for which we had samples from multiple sites and suggest limited connectivity for populations at some sites. Despite relatively small sample sizes, genetic diversity indices suggest the largest population sizes for B. childressi and Tamu fisheri and the smallest for B. heckerae and B. sp. nov. among the GoM bathymodiolins. Moreover, we detected an excess of rare variants indicating recent demographic changes and population expansions for the four species of bathymodiolins from the Gulf of Mexico.

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.

A new barrier to dispersal trapped old genetic clines that escaped the Easter Microplate tension zone of the Pacific vent mussels

Comparative phylogeography of deep-sea hydrothermal vent species has uncovered several genetic breaks between populations inhabiting northern and southern latitudes of the East Pacific Rise. However, the geographic width and position of genetic clines are variable among species. In this report, we further characterize the position and strength of barriers to gene flow between populations of the deep-sea vent mussel Bathymodiolus thermophilus. Eight allozyme loci and DNA sequences of four nuclear genes were added to previously published sequences of the cytochrome c oxidase subunit I gene. Our data confirm the presence of two barriers to gene flow, one located at the Easter Microplate (between 21°33′S and 31°S) recently described as a hybrid zone, and the second positioned between 7°25′S and 14°S with each affecting different loci. Coalescence analysis indicates a single vicariant event at the origin of divergence between clades for all nuclear loci, although the clines are now spatially discordant. We thus hypothesize that the Easter Microplate barrier has recently been relaxed after a long period of isolation and that some genetic clines have escaped the barrier and moved northward where they have subsequently been trapped by a reinforcing barrier to gene flow between 7°25′S and 14°S.

Elongation factor-1α, a putative single-copy nuclear gene, has divergent sets of paralogs in an arachnid

Identification of paralogy in candidate nuclear loci is an important prerequisite in phylogenetics and statistical phylogeography, but one that is often overlooked. One marker commonly assumed to be a single-copy gene and claimed to harbor great utility for inferring recent divergences is elongation factor-1α (EF-1α). To test this hypothesis, we systematically cloned EF-1α in three disjunct populations of the harvestman Metasiro americanus. Here we show that EF-1α has a large number of paralogs in this species. The paralogs do not evolve in a concerted manner, and the paralogs diverged prior to the population divergence. Moreover, the paralogs of M. americanus are not comparable to the highly divergent EF-1α paralogs found in bees and spiders, which are easily recognized and separated through the use of specific primers. We demonstrate statistically that our detection of paralogs cannot be attributed to amplification error. The presence of EF-1α paralogs in M. americanus prevents its use in statistical phylogeography, and the presence of out-paralogs argues against its use in phylogenetic inference among recently diverged clades. These data contradict the common assumption that EF-1α is for most or all taxa a single-copy gene, or that it has a small number of paralogs that are homogenized through gene conversion, unequal crossing over, or other processes.

Phylogeography of the red coral (Corallium rubrum): inferences on the evolutionary history of a temperate gorgonian

The red coral Corallium rubrum (Cnidaria, Octocorallia) is an exploited, long-lived sessile species from the Mediterranean Sea and the adjacent coastline in the Atlantic Ocean. Surveys of genetic variation using microsatellites have shown that populations of C. rubrum are characterized by strong differentiation at the local scale but a study of the phylogeography of this species was still lacking. Here, we used seven polymorphic microsatellite loci, together with sequence data from an intron of the elongation factor 1 (EF1) gene, to investigate the genetic structure of C. rubrum across its geographical range in the western Mediterranean Sea and in the Adriatic Sea. The EF1 sequences were also used to analyse the consequences of demographic fluctuations linked with past environmental change. Clustering analysis with microsatellite loci highlighted three to seven genetic groups with the distinction of North African and Adriatic populations; this distinction appeared significant with AMOVA and differentiation tests. Microsatellite and EF1 data extended the isolation by distance pattern previously observed for this species at the western Mediterranean scale. EF1 sequences confirmed the genetic differentiation observed between most samples with microsatellites. A statistical parsimony network of EF1 haplotypes provided no evidence of high sequence divergence among regions, suggesting no long-term isolation. Selective neutrality tests on microsatellites and EF1 were not significant but should be interpreted with caution in the case of EF1 because of the low sample sizes for this locus. Our results suggest that recent Quaternary environmental fluctuations had a limited impact on the genetic structure of C. rubrum.

Speciation in the deep sea: multi-locus analysis of divergence and gene flow between two hybridizing species of hydrothermal vent mussels

BACKGROUND

Reconstructing the history of divergence and gene flow between closely-related organisms has long been a difficult task of evolutionary genetics. Recently, new approaches based on the coalescence theory have been developed to test the existence of gene flow during the process of divergence. The deep sea is a motivating place to apply these new approaches. Differentiation by adaptation can be driven by the heterogeneity of the hydrothermal environment while populations should not have been strongly perturbed by climatic oscillations, the main cause of geographic isolation at the surface.

METHODOLOGY/PRINCIPAL FINDING

Samples of DNA sequences were obtained for seven nuclear loci and a mitochondrial locus in order to conduct a multi-locus analysis of divergence and gene flow between two closely related and hybridizing species of hydrothermal vent mussels, Bathymodiolus azoricus and B. puteoserpentis. The analysis revealed that (i) the two species have started to diverge approximately 0.760 million years ago, (ii) the B. azoricus population size was 2 to 5 time greater than the B. puteoserpentis and the ancestral population and (iii) gene flow between the two species occurred over the complete species range and was mainly asymmetric, at least for the chromosomal regions studied.

CONCLUSIONS/SIGNIFICANCE

A long history of gene flow has been detected between the two Bathymodiolus species. However, it proved very difficult to conclusively distinguish secondary introgression from ongoing parapatric differentiation. As powerful as coalescence approaches could be, we are left by the fact that natural populations often deviates from standard assumptions of the underlying model. A more direct observation of the history of recombination at one of the seven loci studied suggests an initial period of allopatric differentiation during which recombination was blocked between lineages. Even in the deep sea, geographic isolation may well be a crucial promoter of speciation.

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.