Spatially explicit models of divergence and genome hitchhiking

The role of recombination landscape in species hybridisation and speciation

It is now well recognised that closely related species can hybridize and exchange genetic material, which may promote or oppose adaptation and speciation. In some cases, interspecific hybridisation is very common, making it surprising that species identity is preserved despite active gene exchange. The genomes of most eukaryotic species are highly heterogeneous with regard to gene density, abundance of repetitive DNA, chromatin compactisation etc, which can make certain genomic regions more prone or more resistant to introgression of genetic material from other species. Heterogeneity in local recombination rate underpins many of the observed patterns across the genome (e.g. actively recombining regions are typically gene rich and depleted for repetitive DNA) and it can strongly affect the permeability of genomic regions to interspecific introgression. The larger the region lacking recombination, the higher the chance for the presence of species incompatibility gene(s) in that region, making the entire non- or rarely recombining block impermeable to interspecific introgression. Large plant genomes tend to have highly heterogeneous recombination landscape, with recombination frequently occurring at the ends of the chromosomes and central regions lacking recombination. In this paper we review the relationship between recombination and introgression in plants and argue that large rarely recombining regions likely play a major role in preserving species identity in actively hybridising plant species.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

The mode of speciation during a recent radiation in open-ocean phytoplankton

Despite the enormous ecological importance of marine phytoplankton, surprisingly little is known about how new phytoplankton species originate and evolve in the open ocean, in the absence of apparent geographic barriers that typically act as isolation mechanisms in speciation. To investigate the mechanism of open-ocean speciation, we combined fossil and climatic records from the late Quaternary with genome-wide evolutionary genetic analyses of speciation in the ubiquitous and abundant pelagic coccolithophore genus Gephyrocapsa (including G. huxleyi, formerly known as Emiliania huxleyi). Based on the analysis of 43 sequenced genomes, we report that the best-fitting scenario for all speciation events analyzed included an extended period of complete isolation followed by recent (Holocene) secondary contact, supporting the role of geographic or oceanographic barriers in population divergence and speciation. Consistent with this, fossil data reveal considerable diachroneity of species first occurrence. The timing of all speciation events coincided with glacial phases of glacial-interglacial cycles, suggesting that stronger isolation between the ocean basins and increased segregation of ecological niches during glaciations are important drivers of speciation in marine phytoplankton. The similarity across multiple speciation events implies the generality of this inferred speciation scenario for marine phytoplankton.

Multidimensional divergent selection, local adaptation, and speciation

Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.

Inversions and genomic differentiation after secondary contact: When drift contributes to maintenance, not loss, of differentiation

Due to their effects on reducing recombination, chromosomal inversions may play an important role in speciation by establishing and/or maintaining linked blocks of genes causing reproductive isolation (RI) between populations. This view fits empirical data indicating that inversions typically harbor loci involved in RI. However, previous computer simulations of infinite populations with two to four loci involved in RI implied that, even with gene flux as low as 10^-8 per gamete, per generation between alternative arrangements, inversions may not have large, qualitative advantages over collinear regions in maintaining population differentiation after secondary contact. Here, we report that finite population sizes can help counteract the homogenizing consequences of gene flux, especially when several fitness-related loci reside within the inversion. In these cases, the persistence time of differentiation after secondary contact can be similar to when gene flux is absent and notably longer than the persistence time without inversions. Thus, despite gene flux, population differentiation may be maintained for up to 100,000 generations, during which time new incompatibilities and/or local adaptations might accumulate and facilitate progress toward speciation. How often these conditions are met in nature remains to be determined.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

Demographic history of speciation in a Senecio altitudinal hybrid zone on Mt. Etna

Hybrid zones typically form as a result of species coming into secondary contact, but can also be established in situ as an ecotonal hybrid zone, a situation which has been reported far less frequently. An altitudinal hybrid zone on Mount Etna between two ragwort species (the low elevation Senecio chrysanthemifolius and high elevation S. aethnensis) could potentially represent either of these possibilities. However, a scenario of secondary contact vs. speciation with gene flow has not been explicitly tested. Here, we test these alternatives and demonstrate that the data do not support secondary contact. Furthermore, we report that the previous analyses of speciation history of these species were based on admixed populations, which has led to inflated estimates of ongoing, interspecific gene flow. Our new analyses, based on 'pure' S. aethnensis and S. chrysanthemifolius populations, reveal gene exchange of less than one effective migrant per generation, a level low enough to allow the species to accumulate neutral, genomewide differences. Overall, our results are consistent with a scenario of speciation with gene flow and a divergence time which coincides with the rise of Mt. Etna to altitudes above 2000 m (~150 KY). Further work to quantify the role of adaptation to contrasting environments of high and low altitudes will be needed to support the scenario of recent ecological speciation in this system.

Genomics at the evolving species boundary

Molecular studies on hybridization date back to Dobzhansky who compared chromosomal banding patterns to determine if interspecific gene flow occurred in nature [1]. Now, the advent of high-throughput sequencing provides increasingly fine insights into genomic differentiation between incipient taxa that are changing our view of adaptation and speciation and the links between the two. Empirical data from hybridizing taxa demonstrate highly heterogeneous patterns of genomic differentiation. Although underlining reasons for this heterogeneity are complex, studies of hybridizing taxa offers some of the best insights into the regions of the genome under divergent selection and the role these regions play in species boundaries. The challenge moving forward is to develop a better theoretical framework that fully leverages these powerful natural experiments.

Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae)

BACKGROUND

The evolution of species boundaries and the relative impact of selection and gene flow on genomic divergence are best studied in populations and species pairs exhibiting various levels of divergence along the speciation continuum. We studied species boundaries in Iberian barbels, Barbus and Luciobarbus, a system of populations and species spanning a wide degree of genetic relatedness, as well as geographic distribution and range overlap. We jointly analyze multiple types of molecular markers and morphological traits to gain a comprehensive perspective on the nature of species boundaries in these cyprinid fishes.

RESULTS

Intraspecific molecular and morphological differentiation is visible among many populations. Genomes of all sympatric species studied are porous to gene flow, even if they are not sister species. Compared to their allopatric counterparts, sympatric representatives of different species share alleles and show an increase in all measures of nucleotide polymorphism (S, Hd, K, π and θ). High molecular diversity is particularly striking in L. steindachneri from the Tejo and Guadiana rivers, which co-varies with other sympatric species. Interestingly, different nuclear markers introgress across species boundaries at various levels, with distinct impacts on population trees. As such, some loci exhibit limited introgression and population trees resemble the presumed species tree, while alleles at other loci introgress more freely and population trees reflect geographic affinities and interspecific gene flow. Additionally, extent of introgression decreases with increasing genetic divergence in hybridizing species pairs.

CONCLUSIONS

We show that reproductive isolation in Iberian Barbus and Luciobarbus is not complete and species boundaries are semi-permeable to (some) gene flow, as different species (including non-sister) are exchanging genes in areas of sympatry. Our results support a speciation-with-gene-flow scenario with heterogeneous barriers to gene flow across the genome, strengthening with genetic divergence. This is consistent with observations coming from other systems and supports the notion that speciation is not instantaneous but a gradual process, during which different species are still able to exchange some genes, while selection prevents gene flow at other loci. We also provide evidence for a hybrid origin of a barbel ecotype, L. steindachneri, suggesting that ecology plays a key role in species coexistence and hybridization in Iberian barbels. This ecotype with intermediate, yet variable, molecular, morphological, trophic and ecological characteristics is the local product of introgressive hybridization of L. comizo with up to three different species (with L. bocagei in the Tejo, with L. microcephalus and L. sclateri in the Guadiana). In spite of the homogenizing effects of ongoing gene flow, species can still be discriminated using a combination of morphological and molecular markers. Iberian barbels are thus an ideal system for the study of species boundaries, since they span a wide range of genetic divergences, with diverse ecologies and degrees of sympatry.

Assessing when chromosomal rearrangements affect the dynamics of speciation: implications from computer simulations

Many hypotheses have been put forth to explain the origin and spread of inversions, and their significance for speciation. Several recent genic models have proposed that inversions promote speciation with gene flow due to the adaptive significance of the genes contained within them and because of the effects inversions have on suppressing recombination. However, the consequences of inversions for the dynamics of genome wide divergence across the speciation continuum remain unclear, an issue we examine here. We review a framework for the genomics of speciation involving the congealing of the genome into alternate adaptive states representing species (“genome wide congealing”). We then place inversions in this context as examples of how genetic hitchhiking can potentially hasten genome wide congealing. Specifically, we use simulation models to (i) examine the conditions under which inversions may speed genome congealing and (ii) quantify predicted magnitudes of these effects. Effects of inversions on promoting speciation were most common and pronounced when inversions were initially fixed between populations before secondary contact and adaptation involved many genes with small fitness effects. Further work is required on the role of underdominance and epistasis between a few loci of major effect within inversions. The results highlight five important aspects of the roles of inversions in speciation: (i) the geographic context of the origins and spread of inversions, (ii) the conditions under which inversions can facilitate divergence, (iii) the magnitude of that facilitation, (iv) the extent to which the buildup of divergence is likely to be biased within vs. outside of inversions, and (v) the dynamics of the appearance and disappearance of exceptional divergence within inversions. We conclude by discussing the empirical challenges in showing that inversions play a central role in facilitating speciation with gene flow.

New susceptibility locus for obesity and dyslipidaemia on chromosome 3q22.3

BACKGROUND

The muscle Ras (MRAS) gene resides on chromosome 3q22.3 and encodes a member of the membrane-associated Ras small GTPase proteins, which function as signal transducers in multiple processes including cell growth and differentiation. Its role in cardiovascular disease is not fully understood yet. In a preliminary study in heterozygous familial hypercholesterolaemia, we identified a locus linking the early onset of coronary artery disease (CAD) to chromosome 3q.22 and elected to sequence the MRAS gene using the MegaBACE DNA analysis system. In the present study, we investigated the association of seven single-nucleotide polymorphisms (SNPs) at this locus with CAD and its dyslipidaemia-related risk traits in 4,650 Saudi angiographed individuals using TaqMan assays by the Applied Biosystems real-time Prism 7900HT Sequence Detection System.

RESULTS

Among the studied SNPs, rs6782181 (p = 0.017) and rs9818870T (p = 0.009) were associated with CAD following adjustment for sex, age and other confounding risk factors. The rs6782181_GG also conferred risk for obesity (1,764 cases vs. 2,586 controls) [1.16(1.03-1.30); p = 0.017], hypercholesterolaemia (1,686 vs. 2,744) [1.23(1.02-1.47); p = 0.019], hypertriglyceridaemia (1,155 vs. 3,496) [1.29(1.01-1.45); p = 0.043] and low high-density lipoprotein-cholesterol (lHDL-chol) levels (1,935 vs. 2,401) [1.15(1.02-1.30); p = 0.023] after adjustment. Additionally, rs253662_(CT+TT) [1.16(1.01-1.32); p = 0.030] was associated with lHDL-chol levels. Interestingly, rs253662 (p = 0.014) and rs6782181 (p = 0.019) were protective against acquiring high low-density lipoprotein-cholesterol (hLDL-chol) levels (p = 0.014), while rs1720819 showed similar effects against CAD (p < 0.0001). More importantly, a 7-mer haplotype, ACCTGAC (χ2 = 7.66; p = 0.0056), constructed from the studied SNPs, its 6-mer derivative CCTGAC (χ2 = 6.90; p = 0.0086) and several other shorter derivatives conferred risk for obesity. hLDL-chol was weakly linked to CTAA (χ2 = 3.79; p = 0.052) and CCT (χ2 = 4.32; p = 0.038), while several other haplotypes were protective against both obesity and hLDL-chol level.

CONCLUSION

Our results demonstrate that the genomic locus for the MRAS gene confers risk for CAD, obesity and dyslipidaemia and point to the possible involvement of other genes or regulatory elements at this locus, rather than changes in the M-Ras protein function, in these events.

Genetic hitchhiking and the dynamic buildup of genomic divergence during speciation with gene flow

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. "Divergence hitchhiking" (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. "Genome hitchhiking" (GH) theory emphasizes genome-wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome-wide "congealing" or "coupling" of differentiation (F(ST)) across loci as speciation progressed.