skip to main content


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

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 withMimulus guttatus(yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a majorQTLcausing dense leaf trichomes in thermally adapted plants to a <50‐kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independentM. guttatuspopulations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferredQTLpeak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3MYBtranscription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r2 = 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 ofMimulus guttatusadaptation 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.

 
more » « less
NSF-PAR ID:
10246885
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
25
Issue:
22
ISSN:
0962-1083
Page Range / eLocation ID:
p. 5647-5662
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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 betweenMimulus laciniatusand a sympatric population of its close relativeM. guttatus. These three traits are probably involved inM. 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 denseSNPmarkers 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 fewQTLof large effect including a highly pleiotropicQTLon chromosome 8. ThisQTLregion contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in otherMimulusspecies. 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.

     
    more » « less
  2. Abstract

    Floral attraction traits can significantly affect pollinator visitation patterns, but adaptive evolution of these traits may be constrained by correlations with other traits. In some cases, molecular pathways contributing to floral attraction are well characterized, offering the opportunity to explore loci potentially underlying variation among individuals. Here, we quantify the range of variation in floralUVpatterning (i.e.UV‘bulls‐eye nectar guides) among crop and wild accessions ofBrassica rapa. We then use experimental crosses to examine the genetic architecture, candidate loci and biochemical underpinnings of this patterning as well as phenotypic manipulations to test the ecological impact. We find qualitative variation inUVpatterning between wild (commonly lackingUVpatterns) and crop (commonly exhibitingUVpatterns) accessions. Similar to the majority of crops, recombinant inbred lines (RILs) derived from an oilseed crop × WIfast‐plant®cross exhibitUVpatterns, the size of which varies extensively among genotypes. InRILs, we further observe strong statistical‐genetic andQTLcorrelations within petal morphological traits and within measurements of petalUVpatterning; however, correlations between morphology andUVpatterning are weak or nonsignificant, suggesting thatUVpatterning is regulated and may evolve independently of overall petal size.HPLCanalyses reveal a high concentration of sinapoyl glucose inUV‐absorbing petal regions, which, in concert with physical locations ofUV‐traitQTLs, suggest a regulatory and structural gene as candidates underlying observed quantitative variation. Finally, insects prefer flowers withUVbulls‐eye patterns over those that lack patterns, validating the importance ofUVpatterning in pollen‐limited populations ofB. rapa.

     
    more » « less
  3. Abstract

    Incompletely isolated species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain distinct species in the face of ongoing gene flow. Here, we use field surveys and reduced representation sequencing to characterize the patterns of reproductive isolation, admixture and genomic divergence between populations of the outcrossing wildflowerMimulus guttatusand selfingM. nasutus. Focusing on a single site where these two species have come into secondary contact, we find that phenological isolation is strong, although incomplete, and is likely driven by divergence in response to photoperiod. In contrast to previous field studies, which have suggested that F1‐hybrid formation might be rare, we discover patterns of genomic variation consistent with ongoing introgression. Strikingly, admixed individuals vary continuously from highly admixed to nearly pureM. guttatus, demonstrating ongoing hybridization and asymmetric introgression fromM. nasutusintoM. guttatus. Patterns of admixture and divergence across the genome show that levels of introgression are more variable than expected by chance. Some genomic regions show a reduced introgression, including one region that overlaps a critical photoperiodQTL, whereas other regions show elevated levels of interspecific gene flow. In addition, we observe a genome‐wide negative relationship between absolute divergence and the local recombination rate, potentially indicating natural selection againstM. nasutusancestry inM. guttatusgenetic backgrounds. Together, our results suggest thatMimulusspeciation is both ongoing and dynamic and that a combination of divergence in phenology and mating system, as well as selection against interspecific alleles, likely maintains these sympatric species.

     
    more » « less
  4. Abstract

    Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on atRNAligase gene (Migut.N02091;RLG1a) exhibiting unprecedented, and fitness‐relevant,CNVwithin an annual population of the yellow monkeyflowerMimulus guttatus.RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate‐frequency three‐copy variants ofRLG1a (trip+;5/35 = 14%), andtrip+lines exhibited elevatedRLG1a expression under multiple conditions.trip+carriers, in addition to being over‐represented in late‐flowering and large‐flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rareRLG1a variant (high+) that carries 250–300 copies ofRLG1a totalling ~5.7 Mb (20–40% of a chromosome). In the progeny of ahigh+carrier, Mendelian segregation of diagnostic alleles andqPCR‐based copy counts indicate thathigh+is a single tandem array unlinked to the single‐copyRLG1a locus. In the wild,high+carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; bothp < 0.01), while single‐copy individuals were twice as fecund as eitherCNVtype in a lush year (2016:p < 0.005). Our results demonstrate fluctuating selection onCNVs affecting phenological traits in a wild population, suggest that planttRNAligases mediate stress‐responsive life‐history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

     
    more » « less
  5. Abstract

    Phenotypic plasticity allows organisms to change their phenotype in response to shifts in the environment. While a central topic in current discussions of evolutionary potential, a comprehensive understanding of the genetic underpinnings of plasticity is lacking in systems undergoing adaptive diversification. Here, we investigate the genetic basis of phenotypic plasticity in a textbook adaptive radiation, Lake Malawi cichlid fishes. Specifically, we crossed two divergent species to generate an F3hybrid mapping population. At early juvenile stages, hybrid families were split and reared in alternate foraging environments that mimicked benthic/scraping or limnetic/sucking modes of feeding. These alternate treatments produced a variation in morphology that was broadly similar to the major axis of divergence among Malawi cichlids, providing support for the flexible stem theory of adaptive radiation. Next, we found that the genetic architecture of several morphological traits was highly sensitive to the environment. In particular, of 22 significant quantitative trait loci (QTL), only one was shared between the environments. In addition, we identifiedQTLacting across environments with alternate alleles being differentially sensitive to the environment. Thus, our data suggest that while plasticity is largely determined by loci specific to a given environment, it may also be influenced by loci operating across environments. Finally, our mapping data provide evidence for the evolution of plasticity via genetic assimilation at an important regulatory locus,ptch1. In all, our data address long‐standing discussions about the genetic basis and evolution of plasticity. They also underscore the importance of the environment in affecting developmental outcomes, genetic architectures, morphological diversity and evolutionary potential.

     
    more » « less