Speciation and introgression between Mimulus nasutus and Mimulus guttatus

A carnivorous plant genetic map: pitcher/insect-capture QTL on a genetic linkage map of Sarracenia

This study presents the first genetic map for a carnivorous plant, mapping 64 QTLs in Sarracenia to provide the genetic basis for differences between the pitfall and lobster-trap insect-capture strategies.

The study of carnivorous plants can afford insight into their unique evolutionary adaptations and their interactions with prokaryotic and eukaryotic species. For Sarracenia (pitcher plants), we identified 64 quantitative trait loci (QTL) for insect-capture traits of the pitchers, providing the genetic basis for differences between the pitfall and lobster-trap strategies of insect capture. The linkage map developed here is based upon the F2 of a cross between Sarracenia rosea and Sarracenia psittacina; we mapped 437 single nucleotide polymorphism and simple sequence repeat markers. We measured pitcher traits which differ between S. rosea and S. psittacina, mapping 64 QTL for 17 pitcher traits; there are hot-spot locations where multiple QTL map near each other. There are epistatic interactions in many cases where there are multiple loci for a trait. The QTL map uncovered the genetic basis for the differences between pitfall- and lobster-traps, and the changes that occurred during the divergence of these species. The longevity and clonability of Sarracenia plants make the F2 mapping population a resource for mapping more traits and for phenotype-to-genotype studies.

Local PCA Shows How the Effect of Population Structure Differs Along the Genome

Population structure leads to systematic patterns in measures of mean relatedness between individuals in large genomic data sets, which are often discovered and visualized using dimension reduction techniques such as principal component analysis (PCA). Mean relatedness is an average of the relationships across locus-specific genealogical trees, which can be strongly affected on intermediate genomic scales by linked selection and other factors. We show how to use local PCA to describe this intermediate-scale heterogeneity in patterns of relatedness, and apply the method to genomic data from three species, finding in each that the effect of population structure can vary substantially across only a few megabases. In a global human data set, localized heterogeneity is likely explained by polymorphic chromosomal inversions. In a range-wide data set of Medicago truncatula, factors that produce heterogeneity are shared between chromosomes, correlate with local gene density, and may be caused by linked selection, such as background selection or local adaptation. In a data set of primarily African Drosophila melanogaster, large-scale heterogeneity across each chromosome arm is explained by known chromosomal inversions thought to be under recent selection and, after removing samples carrying inversions, remaining heterogeneity is correlated with recombination rate and gene density, again suggesting a role for linked selection. The visualization method provides a flexible new way to discover biological drivers of genetic variation, and its application to data highlights the strong effects that linked selection and chromosomal inversions can have on observed patterns of genetic variation.

The role of gene flow in rapid and repeated evolution of cave-related traits in Mexican tetra, Astyanax mexicanus

Understanding the molecular basis of repeatedly evolved phenotypes can yield key insights into the evolutionary process. Quantifying gene flow between populations is especially important in interpreting mechanisms of repeated phenotypic evolution, and genomic analyses have revealed that admixture occurs more frequently between diverging lineages than previously thought. In this study, we resequenced 47 whole genomes of the Mexican tetra from three cave populations, two surface populations and outgroup samples. We confirmed that cave populations are polyphyletic and two Astyanax mexicanus lineages are present in our data set. The two lineages likely diverged much more recently than previous mitochondrial estimates of 5-7 mya. Divergence of cave populations from their phylogenetically closest surface population likely occurred between ~161 and 191 k generations ago. The favoured demographic model for most population pairs accounts for divergence with secondary contact and heterogeneous gene flow across the genome, and we rigorously identified gene flow among all lineages sampled. Therefore, the evolution of cave-related traits occurred more rapidly than previously thought, and trogolomorphic traits are maintained despite gene flow with surface populations. The recency of these estimated divergence events suggests that selection may drive the evolution of cave-derived traits, as opposed to disuse and drift. Finally, we show that a key trogolomorphic phenotype QTL is enriched for genomic regions with low divergence between caves, suggesting that regions important for cave phenotypes may be transferred between caves via gene flow. Our study shows that gene flow must be considered in studies of independent, repeated trait evolution.

Gene regulatory divergence between locally adapted ecotypes in their native habitats

Local adaptation is a key driver of ecological specialization and the formation of new species. Despite its importance, the evolution of gene regulatory divergence among locally adapted populations is poorly understood, especially how that divergence manifests in nature. Here, we evaluate gene expression divergence and allele-specific gene expression responses for locally adapted coastal perennial and inland annual accessions of the yellow monkeyflower, Mimulus guttatus, in a field reciprocal transplant experiment. Overall, 6765 (73%) of surveyed genes were differentially expressed between coastal and inland habitats, while 7213 (77%) were differentially expressed between the coastal perennial and inland annual accessions. Cis-regulatory variation was pervasive, affecting 79% (5532) of differentially expressed genes. We detected trans effects for 52% (3611) of differentially expressed genes. Expression plasticity of alleles across habitats (G × E interactions) appears to be relatively common (affecting 18% of transcripts) and could minimize fitness trade-offs at loci that contribute to local adaptation. We also found evidence that at least one chromosomal inversion may act as supergene by holding together haplotypes of differentially expressed genes, but this pattern depends on habitat context. Our results highlight multiple key patterns regarding the relationship between gene expression and the evolution of locally adapted populations.

Digest: Hybrid incompatibilities and introgression in wild monkeyflowers

How are alleles that are detrimental to fitness maintained in natural populations? Zuellig and Sweigart (2018a) find that alleles from a two-locus hybrid incompatibility system segregate at considerable frequencies in two species of monkeyflowers, suggesting that despite providing a fitness cost, these alleles remain polymorphic as a consequence of gene flow between the two species. The system provides the potential to understand the evolutionary trajectory of hybrid incompatibilities and their role in speciation.

Higher Rates of Protein Evolution in the Self-Fertilizing Plant Arabidopsis thaliana than in the Out-Crossers Arabidopsis lyrata and Arabidopsis halleri

The common transition from out-crossing to self-fertilization in plants decreases effective population size. This is expected to result in a reduced efficacy of natural selection and in increased rates of protein evolution in selfing plants compared with their outcrossing congeners. Prior analyses, based on a very limited number of genes, detected no differences between the rates of protein evolution in the selfing Arabidopsis thaliana compared with the out-crosser Arabidopsis lyrata. Here, we reevaluate this trend using the complete genomes of A. thaliana, A. lyrata, Arabidopsis halleri, and the outgroups Capsella rubella and Thellungiella parvula. Our analyses indicate slightly but measurably higher nonsynonymous divergences (dN), synonymous divergences (dS) and dN/dS ratios in A. thaliana compared with the other Arabidopsis species, indicating that purifying selection is indeed less efficacious in A. thaliana.

Model-based detection and analysis of introgressed Neanderthal ancestry in modern humans

Genetic evidence has revealed that the ancestors of modern human populations outside Africa and their hominin sister groups, notably Neanderthals, exchanged genetic material in the past. The distribution of these introgressed sequence tracts along modern-day human genomes provides insight into the selective forces acting on them and the role of introgression in the evolutionary history of hominins. Studying introgression patterns on the X-chromosome is of particular interest, as sex chromosomes are thought to play a special role in speciation. Recent studies have developed methods to localize introgressed ancestries, reporting long regions that are depleted of Neanderthal introgression and enriched in genes, suggesting negative selection against the Neanderthal variants. On the other hand, enriched Neanderthal ancestry in hair- and skin-related genes suggests that some introgressed variants facilitated adaptation to new environments. Here, we present a model-based introgression detection method called dical-admix. We demonstrate its efficiency and accuracy through extensive simulations and apply it to detect tracts of Neanderthal introgression in modern human individuals from the 1000 Genomes Project. Our findings are largely concordant with previous studies, consistent with weak selection against Neanderthal ancestry. We find evidence that selection against Neanderthal ancestry was due to higher genetic load in Neanderthals resulting from small effective population size, rather than widespread Dobzhansky-Müller incompatibilities (DMIs) that could contribute to reproductive isolation. Moreover, we confirm the previously reported low level of introgression on the X-chromosome, but find little evidence that DMIs contributed to this pattern.

A two-locus hybrid incompatibility is widespread, polymorphic, and active in natural populations of Mimulus

Reproductive isolation, which is essential for the maintenance of species in sympatry, is often incomplete between closely related species. In these taxa, reproductive barriers must evolve within species, without being degraded by ongoing gene flow. To better understand this dynamic, we investigated the frequency and geographic distribution of alleles underlying a two-locus, hybrid lethality system between naturally hybridizing species of monkeyflower (Mimulus guttatus and M. nasutus). We found that M. guttatus typically carries hybrid lethality alleles at one locus (hl13) and M. nasutus typically carries hybrid lethality alleles at the other locus (hl14). As a result, natural hybrids carry incompatible alleles at both loci, and express hybrid lethality in later generations. We also discovered considerable polymorphism at both hl13 and hl14 within both species. For M. guttatus, polymorphism at both loci occurs within populations, meaning that incompatible allele pairings likely arise through intraspecific gene flow. Genetic variation at markers linked to hl13 and hl14 suggest that introgression from M. nasutus is the primary driver of this polymorphism within M. guttatus. Additionally, patterns of introgression at the two hybrid lethality loci suggest that natural selection eliminates incompatible allele pairings, suggesting that even weak reproductive barriers might promote genomic divergence between species.

Loter: A Software Package to Infer Local Ancestry for a Wide Range of Species

Admixture between populations provides opportunity to study biological adaptation and phenotypic variation. Admixture studies rely on local ancestry inference for admixed individuals, which consists of computing at each locus the number of copies that originate from ancestral source populations. Existing software packages for local ancestry inference are tuned to provide accurate results on human data and recent admixture events. Here, we introduce Loter, an open-source software package that does not require any biological parameter besides haplotype data in order to make local ancestry inference available for a wide range of species. Using simulations, we compare the performance of Loter to HAPMIX, LAMP-LD, and RFMix. HAPMIX is the only software severely impacted by imperfect haplotype reconstruction. Loter is the less impacted software by increasing admixture time when considering simulated and admixed human genotypes. For simulations of admixed Populus genotypes, Loter and LAMP-LD are robust to increasing admixture times by contrast to RFMix. When comparing length of reconstructed and true ancestry tracts, Loter and LAMP-LD provide results whose accuracy is again more robust than RFMix to increasing admixture times. We apply Loter to individuals resulting from admixture between Populus trichocarpa and Populus balsamifera and lengths of ancestry tracts indicate that admixture took place ∼100 generations ago. We expect that providing a rapid and parameter-free software for local ancestry inference will make more accessible genomic studies about admixture processes.

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less-fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co-occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre-date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.

The genomic impact of historical hybridization with massive mitochondrial DNA introgression

BACKGROUND

The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA.

RESULTS

Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions.

CONCLUSIONS

General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.

Divergent strategies in pre- and postzygotic reproductive isolation between two closely related Dianthus species

Quantifying the relative contribution of multiple isolation barriers to gene flow between recently diverged species is essential for understanding speciation processes. In parapatric populations, local adaptation is thought to be a major contributor to the evolution of reproductive isolation. However, extrinsic postzygotic barriers assessed in reciprocal transplant experiments are often neglected in empirical assessments of multiple isolation barriers. We analyzed multiple isolation barriers between two closely related species of the plant genus Dianthus, a genus characterized by the most rapid species diversification in plants reported so far. Although D. carthusianorum L. and D. sylvestris Wulf. can easily be hybridized in crossing experiments, natural hybrids are rare. We found that in parapatry, pollinator-mediated prezygotic reproductive isolation barriers are important for both D. carthusianorum (0.761) and D. sylvestris (0.468). In contrast to D. carthusianorum, high hybrid viability in D. sylvestris (-0.491) was counteracted by strong extrinsic postzygotic isolation (0.900). Our study highlights the importance of including reciprocal transplant experiments for documenting extrinsic postzygotic isolation and demonstrates clearly divergent strategies and hence asymmetric pre- and postzygotic reproductive isolation between closely related species. It also suggests that pollinator-mediated and ecological isolation could have interacted in synergistic ways, further stimulating rapid speciation in Dianthus.

Variation in recombination frequency and distribution across eukaryotes: patterns and processes

Recombination, the exchange of DNA between maternal and paternal chromosomes during meiosis, is an essential feature of sexual reproduction in nearly all multicellular organisms. While the role of recombination in the evolution of sex has received theoretical and empirical attention, less is known about how recombination rate itself evolves and what influence this has on evolutionary processes within sexually reproducing organisms. Here, we explore the patterns of, and processes governing recombination in eukaryotes. We summarize patterns of variation, integrating current knowledge with an analysis of linkage map data in 353 organisms. We then discuss proximate and ultimate processes governing recombination rate variation and consider how these influence evolutionary processes. Genome-wide recombination rates (cM/Mb) can vary more than tenfold across eukaryotes, and there is large variation in the distribution of recombination events across closely related taxa, populations and individuals. We discuss how variation in rate and distribution relates to genome architecture, genetic and epigenetic mechanisms, sex, environmental perturbations and variable selective pressures. There has been great progress in determining the molecular mechanisms governing recombination, and with the continued development of new modelling and empirical approaches, there is now also great opportunity to further our understanding of how and why recombination rate varies.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Natural selection interacts with recombination to shape the evolution of hybrid genomes

To investigate the consequences of hybridization between species, we studied three replicate hybrid populations that formed naturally between two swordtail fish species, estimating their fine-scale genetic map and inferring ancestry along the genomes of 690 individuals. In all three populations, ancestry from the "minor" parental species is more common in regions of high recombination and where there is linkage to fewer putative targets of selection. The same patterns are apparent in a reanalysis of human and archaic admixture. These results support models in which ancestry from the minor parental species is more likely to persist when rapidly uncoupled from alleles that are deleterious in hybrids. Our analyses further indicate that selection on swordtail hybrids stems predominantly from deleterious combinations of epistatically interacting alleles.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.

Gene duplicates cause hybrid lethality between sympatric species of Mimulus

Hybrid incompatibilities play a critical role in the evolution and maintenance of species. We have discovered a simple genetic incompatibility that causes lethality in hybrids between two closely related species of yellow monkeyflower (Mimulus guttatus and M. nasutus). This hybrid incompatibility, which causes one sixteenth of F₂ hybrid seedlings to lack chlorophyll and die shortly after germination, occurs between sympatric populations that are connected by ongoing interspecific gene flow. Using complimentary genetic mapping and gene expression analyses, we show that lethality occurs in hybrids that lack a functional copy of the critical photosynthetic gene pTAC14. In M. guttatus, this gene was duplicated, but the ancestral copy is no longer expressed. In M. nasutus, the duplication is missing altogether. As a result, hybrids die when they are homozygous for the nonfunctional M. guttatus copy and missing the duplicate from M. nasutus, apparently due to misregulated transcription of key photosynthetic genes. Our study indicates that neutral evolutionary processes may play an important role in the evolution of hybrid incompatibilities and opens the door to direct investigations of their contribution to reproductive isolation among naturally hybridizing species.

The evolution of hybrid incompatibilities (gene interactions that cause hybrids to be sterile or inviable) is a common outcome of genomic divergence between lineages. However, evaluating the importance of hybrid incompatibilities for speciation requires that we identify the causal genes and evolutionary forces in recently diverged, wild species. We discovered that hybrid seedling lethality between two closely related sister species of yellow monkeyflower is caused by duplicate copies of a gene critical for chloroplast development. Because each lineage carries its one functional gene copy in a distinct genomic location, some hybrids inherit only inactive (or missing) alleles. We infer that hybrid lethality in this young species pair has arisen through divergent resolution of gene duplicates by degenerative mutations and (likely) genetic drift. These findings are an important step toward understanding the evolutionary dynamics of hybrid incompatibility genes in nature, as well as the role of such genes in species divergence.

Chloroplastic and nuclear diversity of wild beets at a large geographical scale: Insights into the evolutionary history of the Beta section

Historical demographic processes and mating systems are believed to be major factors in the shaping of the intraspecies genetic diversity of plants. Among Caryophyllales, the Beta section of the genus Beta, within the Amaranthaceae/Chenopodiaceae alliance, is an interesting study model with species and subspecies (Beta macrocarpa, Beta patula, Beta vulgaris maritima and B.v. adanensis) differing in geographical distribution and mating system. In addition, one of the species, B. macrocarpa, mainly diploid, varies in its level of ploidy with a tetraploid cytotype described in the Canary Islands and in Portugal. In this study, we analyzed the nucleotide diversity of chloroplastic and nuclear sequences on a representative sampling of species and subspecies of the Beta section (except B. patula). Our objectives were (1) to assess their genetic relationships through phylogenetic and multivariate analyses, (2) relate their genetic diversity to their mating system, and (3) reconsider the ploidy status and the origin of the Canarian Beta macrocarpa.

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2-3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster- and species-specific marker analyses revealed that, apart from four early-generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1-3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Gene flow analysis method, the D-statistic, is robust in a wide parameter space

Background

We evaluated the sensitivity of the D-statistic, a parsimony-like method widely used to detect gene flow between closely related species. This method has been applied to a variety of taxa with a wide range of divergence times. However, its parameter space and thus its applicability to a wide taxonomic range has not been systematically studied. Divergence time, population size, time of gene flow, distance of outgroup and number of loci were examined in a sensitivity analysis.

Result

The sensitivity study shows that the primary determinant of the D-statistic is the relative population size, i.e. the population size scaled by the number of generations since divergence. This is consistent with the fact that the main confounding factor in gene flow detection is incomplete lineage sorting by diluting the signal. The sensitivity of the D-statistic is also affected by the direction of gene flow, size and number of loci. In addition, we examined the ability of the f-statistics, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\widehat{f}}_G $$\end{document}f^G and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\widehat{f}}_{hom} $$\end{document}f^hom, to estimate the fraction of a genome affected by gene flow; while these statistics are difficult to implement to practical questions in biology due to lack of knowledge of when the gene flow happened, they can be used to compare datasets with identical or similar demographic background.

Conclusions

The D-statistic, as a method to detect gene flow, is robust against a wide range of genetic distances (divergence times) but it is sensitive to population size. The D-statistic should only be applied with critical reservation to taxa where population sizes are large relative to branch lengths in generations.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-2002-4) contains supplementary material, which is available to authorized users.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

Mechanisms of Transmission Ratio Distortion at Hybrid Sterility Loci Within and Between Mimulus Species

Hybrid incompatibilities are a common correlate of genomic divergence and a potentially important contributor to reproductive isolation. However, we do not yet have a detailed understanding of how hybrid incompatibility loci function and evolve within their native species, or why they are dysfunctional in hybrids. Here, we explore these issues for a well-studied, two-locus hybrid incompatibility between hybrid male sterility 1 (hms1) and hybrid male sterility 2 (hms2) in the closely related yellow monkeyflower species Mimulus guttatus and M. nasutus By performing reciprocal backcrosses with introgression lines (ILs), we find evidence for gametic expression of the hms1-hms2 incompatibility. Surprisingly, however, hybrid transmission ratios at hms1 do not reflect this incompatibility, suggesting that additional mechanisms counteract the effects of gametic sterility. Indeed, our backcross experiment shows hybrid transmission bias toward M. guttatus through both pollen and ovules, an effect that is particularly strong when hms2 is homozygous for M. nasutus alleles. In contrast, we find little evidence for hms1 transmission bias in crosses within M. guttatus, providing no indication of selfish evolution at this locus. Although we do not yet have sufficient genetic resolution to determine if hybrid sterility and transmission ratio distortion (TRD) map to the same loci, our preliminary fine-mapping uncovers a genetically independent hybrid lethality system involving at least two loci linked to hms1 This fine-scale dissection of TRD at hms1 and hms2 provides insight into genomic differentiation between closely related Mimulus species and reveals multiple mechanisms of hybrid dysfunction.

Habitat Predicts Levels of Genetic Admixture in Saccharomyces cerevisiae

Genetic admixture can provide material for populations to adapt to local environments, and this process has played a crucial role in the domestication of plants and animals. The model yeast, Saccharomyces cerevisiae, has been domesticated multiple times for the production of wine, sake, beer, and bread, but the high rate of admixture between yeast lineages has so far been treated as a complication for population genomic analysis. Here, we make use of the low recombination rate at centromeres to investigate admixture in yeast using a classic Bayesian approach and a locus-by-locus phylogenetic approach. Using both approaches, we find that S. cerevisiae from stable oak woodland habitats are less likely to show recent genetic admixture compared with those isolated from transient habitats such as fruits, wine, or human infections. When woodland yeast strains do show recent genetic admixture, the degree of admixture is lower than in strains from other habitats. Furthermore, S. cerevisiae populations from oak woodlands are genetically isolated from each other, with only occasional migration between woodlands and local fruit habitats. Application of the phylogenetic approach suggests that there is a previously undetected population in North Africa that is the closest outgroup to the European S. cerevisiae, including the domesticated Wine population. Careful testing for admixture in S. cerevisiae leads to a better understanding of the underlying population structure of the species and will be important for understanding the selective processes underlying domestication in this economically important species.

Population-genomic inference of the strength and timing of selection against gene flow

The interplay of divergent selection and gene flow is key to understanding how populations adapt to local environments and how new species form. Here, we use DNA polymorphism data and genome-wide variation in recombination rate to jointly infer the strength and timing of selection, as well as the baseline level of gene flow under various demographic scenarios. We model how divergent selection leads to a genome-wide negative correlation between recombination rate and genetic differentiation among populations. Our theory shows that the selection density (i.e., the selection coefficient per base pair) is a key parameter underlying this relationship. We then develop a procedure for parameter estimation that accounts for the confounding effect of background selection. Applying this method to two datasets from Mimulus guttatus, we infer a strong signal of adaptive divergence in the face of gene flow between populations growing on and off phytotoxic serpentine soils. However, the genome-wide intensity of this selection is not exceptional compared with what M. guttatus populations may typically experience when adapting to local conditions. We also find that selection against genome-wide introgression from the selfing sister species M. nasutus has acted to maintain a barrier between these two species over at least the last 250 ky. Our study provides a theoretical framework for linking genome-wide patterns of divergence and recombination with the underlying evolutionary mechanisms that drive this differentiation.

African genomes illuminate the early history and transition to selfing in Arabidopsis thaliana

Over the past 20 y, many studies have examined the history of the plant ecological and molecular model, Arabidopsis thaliana, in Europe and North America. Although these studies informed us about the recent history of the species, the early history has remained elusive. In a large-scale genomic analysis of African A. thaliana, we sequenced the genomes of 78 modern and herbarium samples from Africa and analyzed these together with over 1,000 previously sequenced Eurasian samples. In striking contrast to expectations, we find that all African individuals sampled are native to this continent, including those from sub-Saharan Africa. Moreover, we show that Africa harbors the greatest variation and represents the deepest history in the A. thaliana lineage. Our results also reveal evidence that selfing, a major defining characteristic of the species, evolved in a single geographic region, best represented today within Africa. Demographic inference supports a model in which the ancestral A. thaliana population began to split by 120-90 kya, during the last interglacial and Abbassia pluvial, and Eurasian populations subsequently separated from one another at around 40 kya. This bears striking similarities to the patterns observed for diverse species, including humans, implying a key role for climatic events during interglacial and pluvial periods in shaping the histories and current distributions of a wide range of species.

A Hidden Markov Model Approach for Simultaneously Estimating Local Ancestry and Admixture Time Using Next Generation Sequence Data in Samples of Arbitrary Ploidy

Admixture-the mixing of genomes from divergent populations-is increasingly appreciated as a central process in evolution. To characterize and quantify patterns of admixture across the genome, a number of methods have been developed for local ancestry inference. However, existing approaches have a number of shortcomings. First, all local ancestry inference methods require some prior assumption about the expected ancestry tract lengths. Second, existing methods generally require genotypes, which is not feasible to obtain for many next-generation sequencing projects. Third, many methods assume samples are diploid, however a wide variety of sequencing applications will fail to meet this assumption. To address these issues, we introduce a novel hidden Markov model for estimating local ancestry that models the read pileup data, rather than genotypes, is generalized to arbitrary ploidy, and can estimate the time since admixture during local ancestry inference. We demonstrate that our method can simultaneously estimate the time since admixture and local ancestry with good accuracy, and that it performs well on samples of high ploidy-i.e. 100 or more chromosomes. As this method is very general, we expect it will be useful for local ancestry inference in a wider variety of populations than what previously has been possible. We then applied our method to pooled sequencing data derived from populations of Drosophila melanogaster on an ancestry cline on the east coast of North America. We find that regions of local recombination rates are negatively correlated with the proportion of African ancestry, suggesting that selection against foreign ancestry is the least efficient in low recombination regions. Finally we show that clinal outlier loci are enriched for genes associated with gene regulatory functions, consistent with a role of regulatory evolution in ecological adaptation of admixed D. melanogaster populations. Our results illustrate the potential of local ancestry inference for elucidating fundamental evolutionary processes.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

Population-level genetic variation and climate change in a biodiversity hotspot

INTRODUCTION

Estimated future climate scenarios can be used to predict where hotspots of endemism may occur over the next century, but life history, ecological and genetic traits will be important in informing the varying responses within myriad taxa. Essential to predicting the consequences of climate change to individual species will be an understanding of the factors that drive genetic structure within and among populations. Here, I review the factors that influence the genetic structure of plant species in California, but are applicable elsewhere; existing levels of genetic variation, life history and ecological characteristics will affect the ability of an individual taxon to persist in the presence of anthropogenic change.

FACTORS INFLUENCING THE DISTRIBUTION OF GENETIC VARIATION

Persistence in the face of climate change is likely determined by life history characteristics: dispersal ability, generation time, reproductive ability, degree of habitat specialization, plant-insect interactions, existing genetic diversity and availability of habitat or migration corridors. Existing levels of genetic diversity in plant populations vary based on a number of evolutionary scenarios that include endemism, expansion since the last glacial maximum, breeding system and current range sizes.

REGIONAL PRIORITIES AND EXAMPLES

A number of well-documented examples are provided from the California Floristic Province. Some predictions can be made for the responses of plant taxa to rapid environmental changes based on geographic position, evolutionary history, existing genetic variation, and ecological amplitude.

CONCLUSIONS, SOLUTIONS AND RECOMMENDATIONS

The prediction of how species will respond to climate change will require a synthesis drawing from population genetics, geography, palaeontology and ecology. The important integration of the historical factors that have shaped the distribution and existing genetic structure of California's plant taxa will enable us to predict and prioritize the conservation of species and areas most likely to be impacted by rapid climate change, human disturbance and invasive species.

Finding the Genomic Basis of Local Adaptation: Pitfalls, Practical Solutions, and Future Directions

Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.

Hybridization and hybrid speciation under global change

Contents 1170 I. 1170 II. 1172 III. 1175 IV. 1180 V. 1183 1184 References 1184 SUMMARY: An unintended consequence of global change is an increase in opportunities for hybridization among previously isolated lineages. Here we illustrate how global change can facilitate the breakdown of reproductive barriers and the formation of hybrids, drawing on the flora of the British Isles for insight. Although global change may ameliorate some of the barriers preventing hybrid establishment, for example by providing new ecological niches for hybrids, it will have limited effects on environment-independent post-zygotic barriers. For example, genic incompatibilities and differences in chromosome numbers and structure within hybrid genomes are unlikely to be affected by global change. We thus speculate that global change will have a larger effect on eroding pre-zygotic barriers (eco-geographical isolation and phenology) than post-zygotic barriers, shifting the relative importance of these two classes of reproductive barriers from what is usually seen in naturally produced hybrids where pre-zygotic barriers are the largest contributors to reproductive isolation. Although the long-term fate of neo-hybrids is still to be determined, the massive impact of global change on the dynamics and distribution of biodiversity generates an unprecedented opportunity to study large numbers of unpredicted, and often replicated, hybridization 'experiments', allowing us to peer into the birth and death of evolutionary lineages.

The genetic architecture of local adaptation and reproductive isolation in sympatry within the Mimulus guttatus species complex

The genetic architecture of local adaptation has been of central interest to evolutionary biologists since the modern synthesis. In addition to classic theory on the effect size of adaptive mutations by Fisher, Kimura and Orr, recent theory addresses the genetic architecture of local adaptation in the face of ongoing gene flow. This theory predicts that with substantial gene flow between populations local adaptation should proceed primarily through mutations of large effect or tightly linked clusters of smaller effect loci. In this study, we investigate the genetic architecture of divergence in flowering time, mating system-related traits, and leaf shape between Mimulus laciniatus and a sympatric population of its close relative M. guttatus. These three traits are probably involved in M. laciniatus' adaptation to a dry, exposed granite outcrop environment. Flowering time and mating system differences are also reproductive isolating barriers making them 'magic traits'. Phenotypic hybrids in this population provide evidence of recent gene flow. Using next-generation sequencing, we generate dense SNP markers across the genome and map quantitative trait loci (QTLs) involved in flowering time, flower size and leaf shape. We find that interspecific divergence in all three traits is due to few QTL of large effect including a highly pleiotropic QTL on chromosome 8. This QTL region contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in other Mimulus species. Our results are consistent with theory, indicating that local adaptation and reproductive isolation with gene flow should be due to few loci with large and pleiotropic effects.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Genome sequencing of Metrosideros polymorpha (Myrtaceae), a dominant species in various habitats in the Hawaiian Islands with remarkable phenotypic variations

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

Genetic loci with parent-of-origin effects cause hybrid seed lethality in crosses between Mimulus species

In flowering plants, F1 hybrid seed lethality is a common outcome of crosses between closely related diploid species, but the genetic basis of this early-acting and potentially widespread form of postzygotic reproductive isolation is largely unknown. We intercrossed two closely related species of monkeyflower, Mimulus guttatus and Mimulus tilingii, to characterize the mechanisms and strength of postzygotic reproductive isolation. Then, using a reciprocal backcross design, we performed high-resolution genetic mapping to determine the genetic architecture of hybrid seed lethality and directly test for loci with parent-of-origin effects. We found that F1 hybrid seed lethality is an exceptionally strong isolating barrier between Mimulus species, with reciprocal crosses producing < 1% viable seeds. This form of postzygotic reproductive isolation appears to be highly polygenic, indicating that multiple incompatibility loci have accumulated rapidly between these closely related Mimulus species. It is also primarily caused by genetic loci with parent-of-origin effects, suggesting a possible role for imprinted genes in the evolution of Mimulus hybrid seed lethality. Our findings suggest that divergence in loci with parent-of-origin effects, which is probably driven by genomic coevolution within lineages, might be an important source of hybrid incompatibilities between flowering plant species.

Absence of postmating barriers between a selfing vs. outcrossing Chilean Mimulus species pair

PREMISE OF THE STUDY

Reproductive isolation between sympatric species pairs may be maintained by both pre- and postmating barriers. Here we evaluate potential barriers to mating between the outcrossing Mimulus luteus and its more highly selfing sympatric congener, M. cupreus, two members of the South American luteus complex of Mimulus.

METHODS

Seed set was compared following autonomous self-pollination, manual pollination, conspecific outcrossing, and sympatric and allopatric hybridization, for laboratory-maintained inbred lines and wild-collected accessions. Survival and reproductive fitness of hybrids relative to parental species were examined across environments that differed with respect to temperature and soil nutrients, two factors that vary across the ranges of M. luteus and M. cupreus.

KEY RESULTS

Mimulus luteus was minimally capable of autonomous self-fertilization, consistent with reliance on an animal pollinator, whereas M. cupreus was a successful selfer across all tested accessions. Postmating barriers to hybridization are negligible, in both low- and high-stress environments, across multiple sympatric and allopatric populations.

CONCLUSION

As in the North American M. guttatus-M. nasutus species pair, postmating barriers contribute little to isolation between M. luteus and M. cupreus. This result reinforces the importance of premating barriers, specifically species differences in reliance on, and accessibility to, animal pollinators. A unique aspect of the M. luteus-M. cupreus pair is the recent gain of red floral anthocyanin pigmentation in M. cupreus. On the basis of species differences in vegetative anthocyanin production, a facultative stress-protective response, we propose a potential stress-protective role for the constitutive floral anthocyanins of M. cupreus.

A genomic perspective on hybridization and speciation

Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.

Recent selection for self-compatibility in a population of Leavenworthia alabamica

The evolution of self-compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self-incompatibility (SI) to self-fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S-locus. Despite permitting inbreeding, mutations at the S-locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self- and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S-locus, we collected sequence data from a ∼690 Kb region surrounding the S-locus, as well as from regions not linked to the S-locus. To test for recent positive selection acting at the S-locus, we examined polymorphism and the site-frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.

Comparative Transcriptomics Indicates a Role for SHORT VEGETATIVE PHASE (SVP) Genes in Mimulus guttatus Vernalization Response

The timing of reproduction in response to variable environmental conditions is critical to plant fitness, and is a major driver of taxon differentiation. In the yellow monkey flower, Mimulus guttatus, geographically distinct North American populations vary in their photoperiod and chilling (vernalization) requirements for flowering, suggesting strong local adaptation to their surroundings. Previous analyses revealed quantitative trait loci (QTL) underlying short-day mediated vernalization responsiveness using two annual M. guttatus populations that differed in their vernalization response. To narrow down candidate genes responsible for this variation, and to reveal potential downstream genes, we conducted comparative transcriptomics and quantitative PCR (qPCR) in shoot apices of parental vernalization responsive IM62, and unresponsive LMC24 inbred lines grown under different photoperiods and temperatures. Our study identified several metabolic, hormone signaling, photosynthetic, stress response, and flowering time genes that are differentially expressed between treatments, suggesting a role for their protein products in short-day-mediated vernalization responsiveness. Only a small subset of these genes intersected with candidate genes from the previous QTL study, and, of the main candidates tested with qPCR under nonpermissive conditions, only SHORT VEGETATIVE PHASE (SVP) gene expression met predictions for a population-specific short-day-repressor of flowering that is repressed by cold.

Powerful methods for detecting introgressed regions from population genomic data

Understanding the types and functions of genes that are able to cross species boundaries-and those that are not-is an important step in understanding the forces maintaining species as largely independent lineages across the remainder of the genome. With large next-generation sequencing data sets we are now able to ask whether introgression has occurred across the genome, and multiple methods have been proposed to detect the signature of such events. Here, we introduce a new summary statistic that can be used to test for introgression, RNDmin , that makes use of the minimum pairwise sequence distance between two population samples relative to divergence to an outgroup. We find that our method offers a modest increase in power over other, related tests, but that all such tests have high power to detect introgressed loci when migration is recent and strong. RNDmin is robust to variation in the mutation rate, and remains reliable even when estimates of the divergence time between sister species are inaccurate. We apply RNDmin to population genomic data from the African mosquitoes Anopheles quadriannulatus and A. arabiensis, identifying three novel candidate regions for introgression. Interestingly, one of the introgressed loci is on the X chromosome, but outside of an inversion separating these two species. Our results suggest that significant, but rare, sharing of alleles is occurring between species that diverged more than 1 million years ago, and that application of these methods to additional systems are likely to reveal similar results.

Pollinator-mediated assemblage processes in California wildflowers

Community assembly is the result of multiple ecological and evolutionary forces that influence species coexistence. For flowering plants, pollinators are often essential for plant reproduction and establishment, and pollinator-mediated interactions may influence plant community composition. Here, we use null models and community phylogenetic analyses of co-occurrence patterns to determine the role of pollinator-mediated processes in structuring plant communities dominated by congeners. We surveyed three species-rich genera (Limnanthes, Mimulus and Clarkia) with centres of diversity in the Sierra Nevada of California. Each genus contains species that co-flower and share pollinators, and each has a robust phylogeny. Within each genus, we surveyed 44-48 communities at three spatial scales, measured floral and vegetative traits and tested for segregation or aggregation of: (i) species, (ii) floral traits (which are likely to be influenced by pollinators), and (iii) vegetative traits (which are likely affected by other environmental factors). We detected both aggregation and segregation of floral traits that were uncorrelated with vegetative trait patterns; we infer that pollinators have shaped the community assembly although the mechanisms may be varied (competition, facilitation, or filtering). We also found that mating system differences may play an important role in allowing species co-occurrence. Together, it appears that pollinators influence community assemblage in these three clades.

The Limits of Natural Selection in a Nonequilibrium World

Evolutionary theory predicts that factors such as a small population size or low recombination rate can limit the action of natural selection. The emerging field of comparative population genomics offers an opportunity to evaluate these hypotheses. However, classical theoretical predictions assume that populations are at demographic equilibrium. This assumption is likely to be violated in the very populations researchers use to evaluate selection's limits: populations that have experienced a recent shift in population size and/or effective recombination rates. Here we highlight theory and data analyses concerning limitations on the action of natural selection in nonequilibrial populations and argue that substantial care is needed to appropriately test whether species and populations show meaningful differences in selection efficacy. A move toward model-based inferences that explicitly incorporate nonequilibrium dynamics provides a promising approach to more accurately contrast selection efficacy across populations and interpret its significance.