Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation.

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

Genetics and the Evolution of Prezygotic Isolation.

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Ancient hybridization leads to the repeated evolution of red flowers across a monkeyflower radiation.

The reuse of old genetic variation can promote rapid diversification in evolutionary radiations, but in most cases, the historical events underlying this divergence are not known. For example, ancient hybridization can generate new combinations of alleles that sort into descendant lineages, potentially providing the raw material to initiate divergence. In the Mimulus aurantiacus species complex, there is evidence for widespread gene flow among members of this radiation. In addition, allelic variation in the MaMyb2 gene is responsible for differences in flower color between the closely related ecotypes of subspecies puniceus, contributing to reproductive isolation by pollinators. Previous work suggested that MaMyb2 was introgressed into the red-flowered ecotype of puniceus. However, additional taxa within the radiation have independently evolved red flowers from their yellow-flowered ancestors, raising the possibility that this introgression had a more ancient origin. In this study, we used repeated tests of admixture from whole-genome sequence data across this diverse radiation to demonstrate that there has been both ancient and recurrent hybridization in this group. However, most of the signal of this ancient introgression has been removed due to selection, suggesting that widespread barriers to gene flow are in place between taxa. Yet, a roughly 30 kb region that contains the MaMyb2 gene is currently shared only among the red-flowered taxa. Patterns of admixture, sequence divergence, and extended haplotype homozygosity across this region confirm a history of ancient hybridization, where functional variants have been preserved due to positive selection in red-flowered taxa but lost in their yellow-flowered counterparts. The results of this study reveal that selection against gene flow can reduce genomic signatures of ancient hybridization, but that historical introgression can provide essential genetic variation that facilitates the repeated evolution of phenotypic traits between lineages.

Integrating top-down and bottom-up approaches to understand the genetic architecture of speciation across a monkeyflower hybrid zone.

Understanding the phenotypic and genetic architecture of reproductive isolation is a long-standing goal of speciation research. In several systems, large-effect loci contributing to barrier phenotypes have been characterized, but such causal connections are rarely known for more complex genetic architectures. In this study, we combine "top-down" and "bottom-up" approaches with demographic modelling toward an integrated understanding of speciation across a monkeyflower hybrid zone. Previous work suggests that pollinator visitation acts as a primary barrier to gene flow between two divergent red- and yellow-flowered ecotypes of Mimulus aurantiacus. Several candidate isolating traits and anonymous single nucleotide polymorphism loci under divergent selection have been identified, but their genomic positions remain unknown. Here, we report findings from demographic analyses that indicate this hybrid zone formed by secondary contact, but that subsequent gene flow was restricted by widespread barrier loci across the genome. Using a novel, geographic cline-based genome scan, we demonstrate that candidate barrier loci are broadly distributed across the genome, rather than mapping to one or a few "islands of speciation." Quantitative trait locus (QTL) mapping reveals that most floral traits are highly polygenic, with little evidence that QTL colocalize, indicating that most traits are genetically independent. Finally, we find little evidence that QTL and candidate barrier loci overlap, suggesting that some loci contribute to other forms of reproductive isolation. Our findings highlight the challenges of understanding the genetic architecture of reproductive isolation and reveal that barriers to gene flow other than pollinator isolation may play an important role in this system.

Multiple evolutionary origins of glyphosate resistance in Lolium multiflorum.

The multitude of herbicide resistance patterns that have evolved in different weed species is a remarkable example of the rapid adaptation to anthropogenic-driven disturbance. Recently, resistance to glyphosate was identified in multiple populations of Lolium multiflorum in Oregon. We used phenotypic approaches, as well as population genomic and gene expression analyses, to determine whether known mechanisms were responsible for glyphosate resistance and whether resistance phenotypes evolved independently in different populations, and to identify potential loci contributing to resistance. We found no evidence of genetic alterations or expression changes at known target and non-target-site resistance mechanisms of glyphosate. Population genomic analyses indicated that resistant populations tended to have largely distinct ancestry from one another, suggesting that glyphosate resistance did not spread among populations by gene flow. Rather, resistance appears to have evolved independently on different genetic backgrounds. We also detected potential loci associated with the resistance phenotype, some of which encode proteins with potential effects on herbicide metabolism. Our results suggest that Oregon populations of L. multiflorum evolved resistance to glyphosate due to a novel mechanism. Future studies that characterize the gene or genes involved in resistance will be necessary to confirm this conclusion.

Population genomic and historical analysis suggests a global invasion by bridgehead processes in Mimulus guttatus.

Imperfect historical records and complex demographic histories present challenges for reconstructing the history of biological invasions. Here, we combine historical records, extensive worldwide and genome-wide sampling, and demographic analyses to investigate the global invasion of Mimulus guttatus from North America to Europe and the Southwest Pacific. By sampling 521 plants from 158 native and introduced populations genotyped at >44,000 loci, we determined that invasive M. guttatus was first likely introduced to the British Isles from the Aleutian Islands (Alaska), followed by admixture from multiple parts of the native range. We hypothesise that populations in the British Isles then served as a bridgehead for vanguard invasions worldwide. Our results emphasise the highly admixed nature of introduced M. guttatus and demonstrate the potential of introduced populations to serve as sources of secondary admixture, producing novel hybrids. Unravelling the history of biological invasions provides a starting point to understand how invasive populations adapt to novel environments.

Widespread selection and gene flow shape the genomic landscape during a radiation of monkeyflowers.

Speciation genomic studies aim to interpret patterns of genome-wide variation in light of the processes that give rise to new species. However, interpreting the genomic "landscape" of speciation is difficult, because many evolutionary processes can impact levels of variation. Facilitated by the first chromosome-level assembly for the group, we use whole-genome sequencing and simulations to shed light on the processes that have shaped the genomic landscape during a radiation of monkeyflowers. After inferring the phylogenetic relationships among the 9 taxa in this radiation, we show that highly similar diversity (π) and differentiation (FST) landscapes have emerged across the group. Variation in these landscapes was strongly predicted by the local density of functional elements and the recombination rate, suggesting that the landscapes have been shaped by widespread natural selection. Using the varying divergence times between pairs of taxa, we show that the correlations between FST and genome features arose almost immediately after a population split and have become stronger over time. Simulations of genomic landscape evolution suggest that background selection (BGS; i.e., selection against deleterious mutations) alone is too subtle to generate the observed patterns, but scenarios that involve positive selection and genetic incompatibilities are plausible alternative explanations. Finally, tests for introgression among these taxa reveal widespread evidence of heterogeneous selection against gene flow during this radiation. Combined with previous evidence for adaptation in this system, we conclude that the correlation in FST among these taxa informs us about the processes contributing to adaptation and speciation during a rapid radiation.

Variation in ecophysiological traits might contribute to ecogeographic isolation and divergence between parapatric ecotypes of Mimulus aurantiacus.

Many forms of reproductive isolation contribute to speciation, and early-acting barriers may be especially important, because they have the first opportunity to limit gene flow. Ecogeographic isolation occurs when intrinsic traits of taxa contribute to disjunct geographic distributions, reducing the frequency of intertaxon mating. Characterizing this form of isolation requires knowledge of both the geographic arrangement of suitable habitats in nature and the identification of phenotypes involved in shaping geographic distributions. In Mimulus aurantiacus, red- and yellow-flowered ecotypes are incompletely isolated by divergent selection exerted by different pollinators. However, these emerging taxa are largely isolated spatially, with a hybrid zone occurring along a narrow region of contact. In order to assess whether responses to abiotic conditions contribute to the parapatric distribution of ecotypes, we measured a series of ecophysiological traits from populations along a transect, including drought sensitivity, leaf area and the concentrations of vegetative flavonoids. In contrast to the abrupt transitions in floral phenotypes, we found that ecophysiological traits exhibited a continuous geographic transition that largely mirrors variation in climatological variables. These traits may impede gene flow across a continuous environmental gradient, but they would be unlikely to result in ecotypic divergence alone. Nevertheless, we found a genetic correlation between vegetative and floral traits, providing a potential link between the two forms of isolation. Although neither barrier appears sufficient to cause divergence on its own, the combined impacts of local adaptation to abiotic conditions and regional adaptation to pollinators may interact to drive discontinuous variation in the face of gene flow in this system.

Geographic cline analysis as a tool for studying genome-wide variation: a case study of pollinator-mediated divergence in a monkeyflower.

A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome-wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome-wide patterns of divergence. Aided by a first-draft genome assembly, we combined an FCT scan and geographic cline analysis to characterize patterns of genome-wide divergence between divergent pollination ecotypes of Mimulus aurantiacus. FCT analysis of 58 872 SNPs generated via RAD-seq revealed little ecotypic differentiation (mean FCT  = 0.041), although a small number of loci were moderately-to-highly diverged. Consistent with our previous results from the gene MaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under-utilized framework for studying genome-wide patterns of divergence.

Genomewide variation provides insight into evolutionary relationships in a monkeyflower species complex (Mimulus sect. Diplacus).

PREMISE OF THE STUDY: Evolutionary radiations provide excellent opportunities to study the origins of biodiversity, but rapid divergence and ongoing gene flow make inferring evolutionary relationships among taxa difficult. Consequently, combining morphological and genomic analyses will be necessary to clarify the evolutionary history of radiations. We used an integrative approach to shed light on relationships within a diverse radiation of monkeyflowers (Mimulus section Diplacus) with a controversial taxonomic history. METHODS: We used genomewide single nucleotide polymorphism data and a combination of phylogenetic and population genomic analyses to infer the evolutionary relationships within the group. Tests for hybridization were performed to reveal sources of shared variation, and multivariate analyses of floral trait data were conducted to examine the correspondence between phenotypic and phylogenetic data. KEY RESULTS: We identified four primary clades with evidence for some shared variation among them. We also detected evidence for recent gene flow between closely related subclades and populations. Strong discordance between floral trait and molecular data provides evidence for divergent and convergent phenotypic evolution. CONCLUSIONS: Mimulus section Diplacus has all the hallmarks of a rapid radiation, including diverse taxa that are at different stages of divergence, extensive shared variation among taxa, and complex patterns of phenotypic evolution. Our findings will direct future evolutionary research and have important taxonomic implications that highlight the need for a new revision of section Diplacus.