Divergent population structure and climate associations of a chromosomal inversion polymorphism across the Mimulus guttatus species complex

Chromosomal rearrangement polymorphisms are common and increasingly found to be associated with adaptive ecological divergence and speciation. Rearrangements, such as inversions, reduce recombination in heterozygous individuals and thus can protect favourable allelic combinations at linked loci, facilitating their spread in the presence of gene flow. Recently, we identified a chromosomal inversion polymorphism that contributes to ecological adaptation and reproductive isolation between annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus. Here we evaluate the population genetic structure of this inverted region in comparison with the collinear regions of the genome across the M. guttatus species complex. We tested whether annual and perennial M. guttatus exhibit different patterns of divergence for loci in the inverted and noninverted regions of the genome. We then evaluated whether there are contrasting climate associations with these genomic regions through redundancy analysis. We found that the inversion exhibits broadly different patterns of divergence among annual and perennial M. guttatus and is associated with environmental variation across population accessions. This study is the first widespread population genetic survey of the diversity of the M. guttatus species complex. Our findings contribute to a greater understanding of morphological, ecological, and genetic evolutionary divergence across this highly diverse group of closely related ecotypes and species. Finally, understanding species relationships among M. guttatus sp. has hitherto been stymied by accumulated evidence of substantial gene flow among populations as well as designated species. Nevertheless, our results shed light on these relationships and provide insight into adaptation in life history traits within the complex.     

Standing and flowing: the complex origins of adaptive variation

A population faced with a new selection pressure can only adapt if appropriate genetic variation is available. This genetic variation might come from new mutations or from gene exchange with other populations or species, or it might already segregate in the population as standing genetic variation (which might itself have arisen from either mutation or gene flow). Understanding the relative importance of these sources of adaptive variation is a fundamental issue in evolutionary genetics (Orr & Betancourt 2001 ; Barrett & Schluter 2008 ; Gladyshev et al. 2008 ) and has practical implications for conservation, plant and animal breeding, biological control and infectious disease prevention (e.g. Robertson 1960 ; Soulé & Wilcox 1980 ; Prentis et al. 2008 ; Pennings 2012 ). In this issue of Molecular Ecology, Roesti et al. ( 2014 ) make an important contribution to this longstanding debate.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

Opportunities and challenges of next‐generation sequencing applications in ecological epigenetics

Evolutionary theory posits that adaptation can result when populations harbour heritable phenotypic variation for traits that increase tolerance to local conditions. However, the actual mechanisms that underlie heritable phenotypic variation are not completely understood (Keller 2014 ). Recently, the potential role of epigenetic mechanisms in the process of adaptive evolution has been the subject of much debate (Pigliucci & Finkelman 2014 ). Studies of variation in DNA methylation in particular have shown that natural populations harbour high amounts of epigenetic variation, which can be inherited across generations and can cause heritable trait variation independently of genetic variation (Kilvitis et al. 2014 ). While we have made some progress addressing the importance of epigenetics in ecology and evolution using methylation‐sensitive AFLP (MS‐AFLP), this approach provides relatively few anonymous and dominant markers per individual. MS‐AFLP are difficult to link to functional genomic elements or phenotype and are difficult to compare directly to genetic variation, which has limited the insights drawn from studies of epigenetic variation in natural nonmodel populations (Schrey et al. 2013 ). In this issue, Platt et al. provide an example of a promising approach to address this problem by applying a reduced representation bisulphite sequencing (RRBS) approach based on next‐generation sequencing methods in an ecological context.     

Population genomic analyses reveal possible drivers of population divergence

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. 2011 ), Saccharomyces uvarum (Almeida et al. 2014 ), Coccidioides spp. (Neafsey et al. 2010 ) and Cryptococcus gatti (Engelthaler et al. 2014 ). In this issue of Molecular Ecology, Branco et al. ( 2015 ) apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.     

Leaf shape evolution has a similar genetic architecture in three edaphic specialists within the Mimulus guttatus species complex

Background and Aims The genetic basis of leaf shape has long interested botanists because leaf shape varies extensively across the plant kingdom and this variation is probably adaptive. However, knowledge of the genetic architecture of leaf shape variation in natural populations remains limited. This study examined the genetic architecture of leaf shape diversification among three edaphic specialists in the Mimulus guttatus species complex. Lobed and narrow leaves have evolved from the entire, round leaves of M. guttatus in M. laciniatus, M. nudatus and a polymorphic serpentine M. guttatus population (M2L).Methods Bulk segregant analysis and next-generation sequencing were used to map quantitative trait loci (QTLs) that underlie leaf shape in an M. laciniatus × M. guttatus F₂ population. To determine whether the same QTLs contribute to leaf shape variation in M. nudatus and M2L, F₂s from M. guttatus × M. nudatus and lobed M2L × unlobed M. guttatus crosses were genotyped at QTLs from the bulk segregant analysis.Key Results Narrow and lobed leaf shapes in M. laciniatus, M. nudatus and M. guttatus are controlled by overlapping genetic regions. Several promising leaf shape candidate genes were found under each QTL.Conclusions The evolution of divergent leaf shape has taken place multiple times in the M. guttatus species complex and is associated with the occupation of dry, rocky environments. The genetic architecture of elongated and lobed leaves is similar across three species in this group. This may indicate that parallel genetic evolution from standing variation or new mutations is responsible for the putatively adaptive leaf shape variation in Mimulus.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

Pervasive selection or is it…? why are FST outliers sometimes so frequent?

It is now common for population geneticists to estimate F_ST for a large number of loci across the genome, before testing for selected loci as being outliers to the F_ST distribution. One surprising result of such F_ST scans is the often high proportion (>1% and sometimes >10%) of outliers detected, and this is often interpreted as evidence for pervasive local adaptation. In this issue of Molecular Ecolog, Fourcade et al. ( 2013 ) observe that a particularly high rate of F_ST outliers has often been found in river organisms, such as fishes or damselflies, despite there being no obvious reason why selection should affect a larger proportion of the genomes of these organisms. Using computer simulations, Fourcade et al. ( 2013 ) show that the strong correlation in co‐ancestry produced in long one‐dimensional landscapes (such as rivers, valleys, peninsulas, oceanic ridges or coastlines) greatly increases the neutral variance in F_ST, especially when the landscape is further reticulated into fractal networks. As a consequence, outlier tests have a high rate of false positives, unless this correlation can be taken into account. Fourcade et al.'s study highlights an extreme case of the general problem, first noticed by Robertson ( 1975a , b ) and Nei & Maruyama ( 1975 ), that correlated co‐ancestry inflates the neutral variance in F_ST when compared to its expectation under an island model of population structure. Similar warnings about the validity of outlier tests have appeared regularly since then but have not been widely cited in the recent genomics literature. We further emphasize that F_ST outliers can arise in many different ways and that outlier tests are not designed for situations where the genetic architecture of local adaptation involves many loci.     

Genomics and museum specimens

Nearly 25 years ago, Allan Wilson and colleagues isolated DNA sequences from museum specimens of kangaroo rats (Dipodomys panamintinus) and compared these sequences with those from freshly collected animals (Thomas et al. 1990 ). The museum specimens had been collected up to 78 years earlier, so the two samples provided a direct temporal comparison of patterns of genetic variation. This was not the first time DNA sequences had been isolated from preserved material, but it was the first time it had been carried out with a population sample. Population geneticists often try to make inferences about the influence of historical processes such as selection, drift, mutation and migration on patterns of genetic variation in the present. The work of Wilson and colleagues was important in part because it suggested a way in which population geneticists could actually study genetic change in natural populations through time, much the same way that experimentalists can do with artificial populations in the laboratory. Indeed, the work of Thomas et al. ( 1990 ) spawned dozens of studies in which museum specimens were used to compare historical and present‐day genetic diversity (reviewed in Wandeler et al. 2007 ). All of these studies, however, were limited by the same fundamental problem: old DNA is degraded into short fragments. As a consequence, these studies mostly involved PCR amplification of short templates, usually short stretches of mitochondrial DNA or microsatellites. In this issue, Bi et al. ( 2013 ) report a breakthrough that should open the door to studies of genomic variation in museum specimens. They used target enrichment (exon capture) and next‐generation (Illumina) sequencing to compare patterns of genetic variation in historic and present‐day population samples of alpine chipmunks (Tamias alpinus) (Fig. 1). The historic samples came from specimens collected in 1915, so the temporal span of this comparison is nearly 100 years.     

Ontogenetics of QTL: the genetic architecture of trichome density over time in Arabidopsis thaliana

Although much is known about the molecular genetic basis of trichome development in Arabidopsis thaliana, less is known about the underlying genetic basis of continuous variation in a trait known to be of adaptive importance: trichome density. The density of leaf trichomes is known to be a major determinant of herbivore damage in natural populations of A. thaliana and herbivores are a significant selective force on genetic variation for trichome density. A number of developmental changes occur during ontogeny in A. thaliana, including changes in trichome density. I used multiple interval mapping (MIM) analysis to identify QTL responsible for trichome density on both juvenile leaves and adult leaves in replicate, independent trials and asked whether those QTL changed with ontogeny. In both juvenile and adult leaves, I detected a single major QTL on chromosome 2 that explained much of the genetic variance. Although additional QTL were detected, there were no consistent differences in the genetic architecture of trichome density measured on juvenile and adult leaves. The finding of a single QTL of major effect for a trait of known adaptive importance suggests that genes of major effect may play an important role in adaptation.     

Two are better than one: combining landscape genomics and common gardens for detecting local adaptation in forest trees

Predicting likely species responses to an alteration of their local environment is key to decision‐making in resource management, ecosystem restoration and biodiversity conservation practice in the face of global human‐induced habitat disturbance. This is especially true for forest trees which are a dominant life form on Earth and play a central role in supporting diverse communities and structuring a wide range of ecosystems. In Europe, it is expected that most forest tree species will not be able to migrate North fast enough to follow the estimated temperature isocline shift given current predictions for rapid climate warming. In this context, a topical question for forest genetics research is to quantify the ability for tree species to adapt locally to strongly altered environmental conditions (Kremer et al. 2012 ). Identifying environmental factors driving local adaptation is, however, a major challenge for evolutionary biology and ecology in general but is particularly difficult in trees given their large individual and population size and long generation time. Empirical evaluation of local adaptation in trees has traditionally relied on fastidious long‐term common garden experiments (provenance trials) now supplemented by reference genome sequence analysis for a handful of economically valuable species. However, such resources have been lacking for most tree species despite their ecological importance in supporting whole ecosystems. In this issue of Molecular Ecology, De Kort et al. ( 2014 ) provide original and convincing empirical evidence of local adaptation to temperature in black alder, Alnus glutinosa L. Gaertn, a surprisingly understudied keystone species supporting riparian ecosystems. Here, De Kort et al. ( 2014 ) use an innovative empirical approach complementing state‐of‐the‐art landscape genomics analysis of A. glutinosa populations sampled in natura across a regional climate gradient with phenotypic trait assessment in a common garden experiment (Fig. 1 ). By combining the two methods, De Kort et al. ( 2014 ) were able to detect unequivocal association between temperature and phenotypic traits such as leaf size as well as with genetic loci putatively under divergent selection for temperature. The research by De Kort et al. ( 2014 ) provides valuable insight into adaptive response to temperature variation for an ecologically important species and demonstrates the usefulness of an integrated approach for empirical evaluation of local adaptation in nonmodel species (Sork et al. 2013 ).     

Differential adaptation to a harsh granite outcrop habitat between sympatric Mimulus species

Understanding which environmental variables and traits underlie adaptation to harsh environments is difficult because many traits evolve simultaneously as populations or species diverge. Here, we investigate the ecological variables and traits that underlie Mimulus laciniatus’ adaptation to granite outcrops compared to its sympatric, mesic‐adapted progenitor, Mimulus guttatus. We use fine‐scale measurements of soil moisture and herbivory to examine differences in selective forces between the species’ habitats, and measure selection on flowering time, flower size, plant height, and leaf shape in a reciprocal transplant using M. laciniatus × M. guttatus F₄ hybrids. We find that differences in drought and herbivory drive survival differences between habitats, that M. laciniatus and M. guttatus are each better adapted to their native habitat, and differential habitat selection on flowering time, plant stature, and leaf shape. Although early flowering time, small stature, and lobed leaf shape underlie plant fitness in M. laciniatus’ seasonally dry environment, increased plant size is advantageous in a competitive mesic environment replete with herbivores like M. guttatus’. Given that we observed divergent selection between habitats in the direction of species differences, we conclude that adaptation to different microhabitats is an important component of reproductive isolation in this sympatric species pair.     

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.     

High‐throughput SNPs for all: genotyping‐in‐thousands

Understanding the genetic structure of species is essential for conservation. It is only with this information that managers, academics, user groups and land‐use planners can understand the spatial scale of migration and local adaptation, source‐sink dynamics and effective population size. Such information is essential for a multitude of applications including delineating management units, balancing management priorities, discovering cryptic species and implementing captive breeding programmes. Species can range from locally adapted by hundreds of metres (Pavey et al. 2010 ) to complete species panmixia (Côté et al. 2013 ). Even more remarkable is that this essential information can be obtained without fully sequenced or annotated genomes, but from mere (putatively) nonfunctional variants. First with allozymes, then microsatellites and now SNPs, this neutral genetic variation carries a wealth of information about migration and drift. For many of us, it may be somewhat difficult to remember our understanding of species conservation before the widespread usage of these useful tools. However most species on earth have yet to give us that ‘peek under the curtain’. With the current diversity on earth estimated to be nearly 9 million species (Mora et al. 2011 ), we have a long way to go for a comprehensive meta‐phylogeographic understanding. A method presented in this issue by Campbell and colleagues (Campbell et al. 2015 ) is a tool that will accelerate the pace in this area. Genotyping‐in‐thousands (GT‐seq) leverages recent advancements in sequencing technology to save many hours and dollars over previous methods to generate this important neutral genetic information.     

Is local adaptation in Mimulus guttatus caused by trade‐offs at individual loci?

Local adaptation is considered to be the result of fitness trade‐offs for particular phenotypes across different habitats. However, it is unclear whether such phenotypic trade‐offs exist at the level of individual genetic loci. Local adaptation could arise from trade‐offs of alternative alleles at individual loci or by complementary sets of loci with different fitness effects of alleles in one habitat but selective neutrality in the alternative habitat. To evaluate the genome‐wide basis of local adaptation, we performed a field‐based quantitative trait locus (QTL) mapping experiment on recombinant inbred lines (RILs) created from coastal perennial and inland annual races of the yellow monkeyflower (Mimulus guttatus) grown reciprocally in native parental habitats. Overall, we detected 19 QTLs affecting one or more of 16 traits measured in two environments, most of small effect. We identified 15 additional QTL effects at two previously identified candidate QTLs [DIV ERGENCE (DIV)]. Significant QTL by environment interactions were detected at the DIV loci, which was largely attributable to genotypic differences at a single field site. We found no detectable evidence for trade‐offs for any one component of fitness, although DIV2 showed a trade‐off involving different fitness traits between sites, suggesting that local adaptation is largely controlled by non‐overlapping loci. This is surprising for an outcrosser, implying that reduced gene flow prevents the evolution of individuals adapted to multiple environments. We also determined that native genotypes were not uniformly adaptive, possibly reflecting fixed mutational load in one of the populations.     

Admixture increases diversity in managed honey bees: Reply to De la Rúa et al. (2013)

De la Rúa et al. (2013) express some concerns about the conclusions of our recent study showing that management increases genetic diversity of honey bees (Apis mellifera) by promoting admixture (Harpur et al. 2012). We provide a brief review of the literature on the population genetics of A. mellifera and show that we utilized appropriate sampling methods to estimate genetic diversity in the focal populations. Our finding of higher genetic diversity in two managed A. mellifera populations on two different continents is expected to be the norm given the large number of studies documenting admixture in honey bees. Our study focused on elucidating how management affects genetic diversity in honey bees, not on how to best manage bee colonies. We do not endorse the intentional admixture of honey bee populations, and we agree with De la Rúa et al. (2013) that native honey bee subspecies should be conserved.     

Expanding the population genetic perspective of cnidarian‐Symbiodinium symbioses

The modern synthesis was a seminal period in the biological sciences, establishing many of the core principles of evolutionary biology that we know today. Significant catalysts were the contributions of R.A. Fisher, J.B.S. Haldane and Sewall Wright (and others) developing the theoretical underpinning of population genetics, thus demonstrating adaptive evolution resulted from the interplay of forces such as natural selection and mutation within groups of individuals occupying the same space and time (i.e. a population). Given its importance, it is surprising that detailed population genetic data remain lacking for numerous organisms vital to many ecosystems. For example, the coral reef ecosystem is well recognized for its high biodiversity and productivity, numerous ecological services and significant economic and societal values (Moberg & Folke 1999; Cinner 2014). Many coral reef invertebrates form symbiotic relationships with single‐celled dinoflagellates within the genus Symbiodinium Freudenthal (Taylor 1974), with hosts providing these (typically) intracellular symbionts with by‐products of metabolism and in turn receiving photosynthetically fixed carbon capable of meeting hosts’ respiratory demands (Falkowski et al. 1984; Muscatine et al. 1984). Unfortunately, the health and integrity of the coral reef ecosystem has been significantly and negatively impacted by onslaughts like anthropogenic eutrophication and disease in addition to global climate change, with increased incidences of ‘bleaching’ events (characterized as the loss of photosynthetic pigments from the algal cell or massive reduction of Symbiodinium density from hosts’ tissue) and host mortality leading to staggering declines in geographic coverage (Bruno & Selig 2007) that have raised questions on the viability of this ecosystem as we know it (Bellwood et al. 2004; Parmesan 2006). One avenue towards anticipating the future of the coral reef ecosystem is by developing a broader and deeper understanding of the current genotypic diversity encompassed within and between populations of their keystone species, the scleractinian corals and dinoflagellate symbionts, as they potentially possess functional variation (either singularly or in combination) that may come under selection due to the ongoing and rapid environmental changes they are experiencing. However, such studies, especially for members of the genus Symbiodinium, are sparse. In this issue, Baums et al. (2014) provide a significant contribution by documenting the range‐wide population genetics of Symbiodinium ‘fitti’ (Fig. 1 ) in the context of complementary data from its host, the endangered Caribbean elkhorn coral Acropora palmata (Fig. 1 ). Notable results of this study include a single S. ‘fitti’ genotype typically dominates an individual A. palmata colony both spatially and temporally, gene flow among coral host populations is a magnitude higher to that of its symbiont populations, and the partners possess disparate patterns of genetic differentiation across the Greater Caribbean. The implications of such findings are discussed herein.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.