Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

Patterns of nucleotide diversity in two species of Mimulus are affected by mating system and asymmetric introgression

The evolutionary transition from outcrossing to self-fertilization has far-reaching implications for patterns of intraspecific genetic diversity and the potential for speciation. Using DNA sequence variation at two nuclear loci, we examined the divergence history of two closely related species of Mimulus. To investigate the effects of mating system and introgressive hybridization on the outcrossing M. guttatus and the selfing M. nasutus, we inspected nucleotide diversity within and between natural populations spanning the species' geographic ranges. High sequence similarity among populations of the selfing M. nasutus points to a single evolutionary origin for the species. Consistent with their distinct mating systems, all genetic variation in M. nasutus is distributed among populations, whereas M. guttatus exhibits appreciable levels of nucleotide diversity within populations. Silent genetic diversity is extensive in M. guttatus (mean theta(sil)/site = 0.077) and greatly exceeds the predicted twofold elevation in neutral variation for outcrossers relative to selfers. The finding of several M. guttatus sequences that share complete identity with sequences from M. nasutus suggests that recent asymmetric introgression may have occurred. We argue that exceptionally high nucleotide diversity in M. guttatus is consistent with a long-term history of directional introgression from M. nasutus to M. guttatus throughout the divergence of these two species.     

Genetic evidence on the demography of speciation in allopatric dolphin species

Under a neutral model, the stochastic lineage sorting that leads to gene monophyly proceeds slowly in large populations. Therefore, in many recent species with large population size, the genome will have mixed support for monophyly unless historical bottlenecks have accelerated coalescence. We use genealogical patterns in mitochondrial DNA and in introns of four nuclear loci to test for historical bottlenecks during the speciation and divergence of two temperate Lagenorhynchus dolphin species isolated by tropical Pacific waters (an antitropical distribution). Despite distinct morphologies, foraging behaviors, and mitochondrial DNAs, these dolphin species are polyphyletic at all four nuclear loci. The abundance of shared polymorphisms between these sister taxa is most consistent with the maintenance of large effective population sizes (5.09 x 10(4) to 10.9 x 10(4)) during 0.74-1.05 million years of divergence. A variety of population size histories are possible, however. We used gene tree coalescent probabilities to explore the rejection region for historical bottlenecks of different intensity given best estimates of effective population size under a strict isolation model of divergence. In L. obliquidens the data are incompatible with a colonization propagule of an effective size of 10 or fewer individuals. Although the ability to reject less extreme historical bottlenecks will require data from additional loci, the intermixed genealogical patterns observed between these dolphin sister species are highly probable only under an extended history of large population size. If similar demographic histories are inferred for other marine antitropical taxa, a parsimonious model for the Pleistocene origin of these distributions would not involve rare breaches of a constant dispersal barrier by small colonization propagules. Instead, a history of large population size in L. obliquidens and L. obscurus contributes to growing biological and environmental evidence that the equatorial barrier became permeable during glacial/interglacial cycles, leading to vicariant isolation of antitropical populations.     

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.     

Genealogical footprints of speciation processes in wild tomatoes: demography and evidence for historical gene flow

Multilocus studies assessing patterns of nucleotide polymorphism within and among closely related species provide access to genealogical information bearing on demographic and geographic aspects of their speciation history. However, the technical difficulties in obtaining sufficient sequence data have severely limited this approach thus far, especially in outbred plant taxa. We employ the analytical framework of divergence population genetics in testing the isolation model of speciation in three self-incompatible species of wild tomatoes (clade Lycopersicon), in particular the assumption of divergence without gene flow. Based on DNA sequence data for 13 nuclear loci, average levels of silent polymorphism vary more than three-fold among species. We estimate a large effective population size for the ancestral species, quite similar to that of the highly polymorphic L. peruvianum. The other two species, however, exhibit concordant signatures of population-size reduction. These demographic inferences are biologically plausible and consistent with results obtained from standard neutrality tests. While the isolation model cannot be rejected by goodness-of-fit criteria, patterns of intragenic linkage disequilibrium in L. peruvianum are indicative of historical introgression at least in some regions of the genome. Considered jointly with the geographic pattern of postzygotic reproductive isolation, our results suggest that speciation occurred under residual gene flow, implying natural selection as one of the evolutionary forces driving the diversification of tomato lineages.     

Small body size increases the regional differentiation of populations of tropical mantellid frogs (Anura: Mantellidae)

The processes affecting species diversification may also exert an influence on patterns of genetic variability within species. We evaluated the contributions of five variables potentially influencing clade diversification (body size, reproductive mode, range size, microhabitat and skin texture) on mtDNA divergence and polymorphism among populations of 40 species of frogs (Mantellidae) from two rainforest communities in Madagascar. We report an inverse association between body size and nucleotide divergence between populations but find no influence of other variables on genetic variation. Body size explained ca. 11% of the variation in nucleotide divergence between populations and was coupled with high F_ST levels and an absence of haplotype sharing in small‐bodied and medium‐sized frogs. Low dispersal ability is likely the proximate mechanism producing higher population differentiation in small mantellids. The lack of genetic cohesion among populations establishes regional genetic fragmentation which in turn has the potential to accelerate rates of allopatric speciation in small frogs relative to large species. However, there is little evidence of increased speciation rates in these or other small‐bodied organisms. We reconcile these contradictory observations by suggesting that lower dispersal ability also curbs colonization of new areas, decelerating diversification in weak dispersers. Our results imply that the intermediate dispersal model also applies to amphibians and may explain inconsistent previous results on the correlation of body size and speciation rate.     

DNA variation in a conifer,Cryptomeria japonica (Cupressaceae sensu lato)

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.     

Nucleotide variation in Quercus crispula Blume

Quercus is attractive for evolutionary studies, primarily for developing the concepts of the species, speciation and adaptation; however, remarkably little is known about levels of nucleotide polymorphism in the nuclear functional genes of this genus. This article provides the first characterization of levels of nucleotide polymorphism in 11 gene fragments in natural populations of a Quercus species, Quercus crispula Blume. Results show that the level of nucleotide variation in this oak is generally higher than that in conifers, as high as that in a European oak, but lower than that in an aspen. The level of population recombination is relatively high. Within-population inbreeding is negligible and between-population differentiation is modest. The decay of linkage disequilibrium is significantly faster in the species-wide samples and the three northernmost populations than in the other populations. Statistical tests support the hypothesis of a recent bottleneck for several populations in the southern part of Japan. The amounts and patterns of nucleotide variation, recombination and linkage disequilibrium, and genetic differentiation observed among populations of this species are contradictory to our expectations, given the recent colonization history of the northern Japan populations.     

Multilocus analysis of nucleotide variation and speciation in Oryza officinalis and its close relatives

Nucleotide variation in 10 unlinked nuclear genes was investigated in species-wide samples of Oryza officinalis and its close relatives (Oryza eichingeri and Oryza rhizomatis). Average estimates of nucleotide diversity were the lowest in O. rhizomatis ((sil) = 0.0038) and the highest in O. eichingeri ((sil) = 0.0057) that is disjunctly distributed in Africa and Sri Lanka. These wild rice species appeared to harbor relatively low levels of nucleotide variation relative to other plant species because the diversity level of O. eichingeri is only 23-46% of those in Zea species and 35% of that in Arabidopsis thaliana. The lower nucleotide diversity in these Oryza species could be best explained by their smaller historic effective population sizes. The speciation model test indicated that O. officinalis and its close relatives might have undergone a process of population contraction since divergence from their ancestor. Incongruent topologies among 10 gene trees, particularly regarding the positions of O. eichingeri and O. rhizomatis accessions might be attributed to lineage sorting arising from ancient polymorphism and hybridization/introgression between the Sri Lankan O. eichingeri and O. rhizomatis. However, the null hypothesis of the isolation model was not rejected for any contrast between taxa, which suggested that no subsequent gene flow shaped the present patterns of nucleotide variation since their divergence and that introgression was not pervasive in this group of species. Our molecular dating provides an approximate divergence time of 0.37 Myr between 2 geographical races of O. eichingeri, much more recent compared with the times of other speciation events in this group (0.63-0.68 Myr). A long-distance dispersal from West Africa to Sri Lanka was more likely to play a role in the disjunct distribution of O. eichingeri.     

Cryptic speciation on the high seas; global phylogenetics of the copepod family Eucalanidae

Few genetic data are currently available to assess patterns of population differentiation and speciation in planktonic taxa that inhabit the open ocean. A phylogenetic study of the oceanic copepod family Eucalanidae was undertaken to develop a model zooplankton taxon in which speciation events can be confidently identified. A global survey of 20 described species (526 individuals) sampled from 88 locations worldwide found high levels of cryptic diversity at the species level. Mitochondrial (16S rRNA, CO1) and nuclear (ITS2) DNA sequence data support 12 new genetic lineages as highly distinct from other populations with which they are currently considered conspecific. Out of these 12, at least four are new species. The circumglobal, boundary current species Rhincalanus nasutus was found to be a cryptic species complex, with genetic divergence between populations unrelated to geographic distance. 'Conspecific' populations of seven species exhibited varying levels of genetic differentiation between Atlantic and Pacific basins, suggesting that continental landmasses form barriers to dispersal for a subset of circumglobal species. A molecular phylogeny of the family based on both mitochondrial (16S rRNA) and nuclear (ITS2, 18S rRNA) gene loci supports monophyly of the family Eucalanidae, all four eucalanid genera and the 'pileatus' and 'subtenuis' species groups.     

Population genetics of speciation in two closely related wild tomatoes (Solanum section Lycopersicon)

We present a multilocus sequencing study to assess patterns of polymorphism and divergence in the closely related wild tomato species, Solanum peruvianum and S. chilense (Solanum section Lycopersicon, Solanaceae). The data set comprises seven mapped nuclear loci (approximately 9.3 kb of analyzed sequence across loci) and four local population samples per species that cover much of the species' range (between 80 and 88 sequenced alleles across both species). We employ the analytical framework of divergence population genetics (DPG) in evaluating the utility of the "isolation" model of speciation to explain observed patterns of polymorphism and divergence. Whereas the isolation model is not rejected by goodness-of-fit criteria established via coalescent simulations, patterns of intragenic linkage disequilibrium provide evidence for postdivergence gene flow at two of the seven loci. These results suggest that speciation occurred under residual gene flow, implying that natural selection is one of the evolutionary forces driving the divergence of these tomato species. This inference is fully consistent with their recent divergence, conservatively estimated to be 相似文献   

Nucleotide Variation, Linkage Disequilibrium and Founder-Facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (K_ST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.     

Latitudinal variation in genetic divergence of populations and the potential for future speciation

The increase in biological diversity with decreasing latitude is widely appreciated but the cause of the pattern is unknown. This pattern reflects latitudinal variation in both the origin of new species (cladogenesis) and the number of species that coexist. Here we address latitudinal variation in species origination, by examining population genetic processes that influence speciation. Previous data suggest a greater number of speciation events at lower latitudes. If speciation events occur more frequently at lower latitudes, we predicted that genetic divergence among populations within species, an important component of cladogenesis, should be greater among lower latitude populations. We tested this prediction using within-species patterns of mtDNA variation across 60 vertebrate species that collectively spanned six continents, two oceans, and 119 degrees latitude. We found greater genetic divergence of populations, controlling for geographic distance, at lower latitudes within species. This pattern remained statistically significant after removing populations that occur in localities previously covered by continental glaciers during the last glaciation. Results suggest that lower latitude populations within species exhibit greater evolutionary independence, increasing the likelihood that mutation, recombination, selection, and/or drift will lead to divergence of traits important for reproductive isolation and speciation. Results are consistent with a greater influence of seasonality, reduced energy, and/or glacial (Milankovitch) cycles acting on higher latitude populations, and represent one of the few tests of predictions of latitudinal variation in speciation rates using population genetic data.     

Polyploid speciation did not confer instant reproductive isolation in Capsella (Brassicaceae)

Polyploid formation is a major mode of sympatric speciation in flowering plants. Unlike other speciation processes, polyploidization is often assumed to confer instant reproductive isolation. Shared polymorphism across ploidy levels has therefore often been attributed to multiple polyploid origins, whereas the alternative hypothesis of introgressive hybridization has rarely been rigorously tested. Here, we sequence 12 nuclear loci representing 6 genes duplicated by polyploidy in 92 accessions of the tetraploid Capsella bursa-pastoris together with the corresponding loci in 21 accessions of its close diploid relative Capsella rubella. In C. bursa-pastoris accessions from western Eurasia, where the 2 species occur in partial sympatry, we find higher levels of nucleotide diversity than in accessions from eastern Eurasia, where C. rubella does not grow. Furthermore, haplotypes are shared across ploidy levels at 4 loci in western but not in eastern Eurasia. We test whether haplotype sharing is due to retention of ancestral polymorphism or due to hybridization and introgression using a coalescent-based isolation-with-migration model. In western but not in eastern Eurasia, there is evidence for unidirectional gene flow from C. rubella to C. bursa-pastoris. An independent estimate of the timing of dispersal of C. bursa-pastoris to eastern Eurasia indicates that it probably predated introgression. Our results show that polyploid speciation need not result in immediate and complete reproductive isolation, that postpolyploidization hybridization and introgression can contribute significantly to genetic variation in a newly formed polyploid, and that divergence population genetic analysis constitutes a powerful way of testing hypotheses on polyploid speciation.     

Recent radiation of endemic Caribbean Drosophila of the dunni subgroup inferred from multilocus DNA sequence variation

Studies of island endemism provide a unique opportunity to elucidate fundamental mechanisms of speciation. Here we examine intra- and interspecific DNA sequence variation at four unlinked genetic loci among populations of the Drosophila dunni subgroup to provide a detailed genealogical portrait of the process of speciation among these island endemic species. Our data indicate two major rounds of diversification that have shaped the D. dunni subgroup. The first occurred 1.6-2.6 million years ago and separated three major lineages, one in Puerto Rico and the Virgin Islands, a second in the northern Lesser Antilles and Barbados, and a third in St. Vincent and Grenada. A second round of diversification occurred in the last 96,000 years in the northern Lesser Antilles and Barbados. The four distinct species that resulted from this recent round of diversification maintain relatively high amounts of genetic variation, similar to that of a closely related mainland species, and share extensive ancestral polymorphism. These data suggest a minimal role for population bottlenecks linked to founder events in the history of the D. dunni subgroup. Further, the recent divergence of these island populations highlights the extremely rapid development of reproductive isolation and distinct patterns of abdominal pigmentation that has occurred in these species.     

The Pleistocene glacial cycles shaped the historical demography and phylogeography of a pine fungal endophyte

The fungal endophyte Lophodermium nitens is an obligate symbiont of soft pines inhabiting only two pine species in Mexico with a broad distribution of geographically isolated populations. A previous study for the hosts indicated a main east–west subdivision with recurrent gene flow within these regions and demographic expansion of populations. We took these patterns as null hypotheses to test for the demography and phylogeographical patterns of the fungus, given the obligatory relationship of the endophyte to the host and its reduced capacity for long-distance dispersal. For this purpose, we employed two nuclear DNA loci, fragments of the actin and chitin synthase I genes. Both loci showed high genetic variation, consisting of private single-copy alleles, as well as few ones at high frequency that were shared among almost all populations. In order to distinguish between shared polymorphism due to incomplete lineage sorting and gene flow posterior to population divergence, we applied the coalescent-based Isolation–Migration (IM) model. We found patterns of gene flow and isolation similar to those of the hosts as well as signs of population expansion. Mean migration time and divergence time estimates fell within the Pleistocene, previous to Last Glacial Maximum. The results presented here for L. nitens emphasize the potential use of endophytic fungi to deepen the knowledge of historical patterns and processes of their host plants.     

Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.     

The effects of subdivision on the genetic divergence of populations and species

Abstract. An island model of migration is used to study the effects of subdivision within populations and species on sample genealogies and on between-population or between-species measures of genetic variation. The model assumes that the number of demes within each population or species is large. When populations (or species), connected either by gene flow or historical association, are themselves subdivided into demes, changes in the migration rate among demes alter both the structure of genealogies and the time scale of the coalescent process. The time scale of the coalescent is related to the effective size of the population, which depends on the migration rate among demes. When the migration rate among demes within populations is low, isolation (or speciation) events seem more recent and migration rates among populations seem higher because the effective size of each population is increased. This affects the probability of reciprocal monophyly of two samples, the chance that a gene tree of a sample matches the species tree, and relative likelihoods of different types of polymorphic sites. It can also have a profound effect on the estimation of divergence times.     

Demographic history of a common pioneer tree,Zanthoxylum ailanthoides,reconstructed using isolation-with-migration model

The Japanese prickly ash, Zanthoxylum ailanthoides Siebold & Zucc. (Rutaceae), is one of the most common pioneer tree species that grow in disturbed areas, such as canopy gaps, in Japanese warm-temperate evergreen oak forests. The strong genetic structure of its current population might suggest long-term isolation of subpopulations by demographic events, in addition to ecological features of this species. Analyses of genome-wide nucleotide variation using model-based approaches are necessary to achieve a good understanding of the demographic history of a species. In this study, we analyzed the nucleotide variation in natural populations of Z. ailanthoides using a computer program that applied the isolation-with-migration model to quantitatively infer the demographic history. Nucleotide variations at 10 or 26 nuclear loci in six populations from a wide range of the species distribution were analyzed. The maximum likelihood estimate of the divergence time between the current populations in the two hypothetical refugia of Z. ailanthoides was approximately 24,000 years. Unexpectedly, the estimated size of the ancestral population was larger than the sizes of the two current populations. The results suggested a relatively recent divergence of these two populations and rapid formation of strong genetic structure among subpopulations. The large ancestral population may indicate a more complex demographic history during or after the last glacial period than the simple isolation-with-migration model implies.     

Population divergence with or without admixture: selecting models using an ABC approach

Genetic data have been widely used to reconstruct the demographic history of populations, including the estimation of migration rates, divergence times and relative admixture contribution from different populations. Recently, increasing interest has been given to the ability of genetic data to distinguish alternative models. One of the issues that has plagued this kind of inference is that ancestral shared polymorphism is often difficult to separate from admixture or gene flow. Here, we applied an approximate Bayesian computation (ABC) approach to select the model that best fits microsatellite data among alternative splitting and admixture models. We performed a simulation study and showed that with reasonably large data sets (20 loci) it is possible to identify with a high level of accuracy the model that generated the data. This suggests that it is possible to distinguish genetic patterns due to past admixture events from those due to shared polymorphism (population split without admixture). We then apply this approach to microsatellite data from an endangered and endemic Iberian freshwater fish species, in which a clustering analysis suggested that one of the populations could be admixed. In contrast, our results suggest that the observed genetic patterns are better explained by a population split model without admixture.