The Evolutionary Forces Shaping Cis- and Trans-Regulation of Gene Expression within a Population of Outcrossing Plants.

Understanding the persistence of genetic variation within populations has long been a goal of evolutionary biology. One promising route toward achieving this goal is using population genetic approaches to describe how selection acts on the loci associated with trait variation. Gene expression provides a model trait for addressing the challenge of the maintenance of variation because it can be measured genome-wide without information about how gene expression affects traits. Previous work has shown that loci affecting the expression of nearby genes (local or cis-eQTLs) are under negative selection, but we lack a clear understanding of the selective forces acting on variants that affect the expression of genes in trans. Here, we identify loci that affect gene expression in trans using genomic and transcriptomic data from one population of the obligately outcrossing plant, Capsella grandiflora. The allele frequencies of trans-eQTLs are consistent with stronger negative selection acting on trans-eQTLs than cis-eQTLs, and stronger negative selection acting on trans-eQTLs associated with the expression of multiple genes. However, despite this general pattern, we still observe the presence of a trans-eQTL at intermediate frequency that affects the expression of a large number of genes in the same coexpression module. Overall, our work highlights the different selective pressures shaping variation in cis- and trans-regulation.

Standing genetic variation in a tissue-specific enhancer underlies selfing-syndrome evolution in Capsella.

Mating system shifts recurrently drive specific changes in organ dimensions. The shift in mating system from out-breeding to selfing is one of the most frequent evolutionary transitions in flowering plants and is often associated with an organ-specific reduction in flower size. However, the evolutionary paths along which polygenic traits, such as size, evolve are poorly understood. In particular, it is unclear how natural selection can specifically modulate the size of one organ despite the pleiotropic action of most known growth regulators. Here, we demonstrate that allelic variation in the intron of a general growth regulator contributed to the specific reduction of petal size after the transition to selfing in the genus Capsella Variation within this intron affects an organ-specific enhancer that regulates the level of STERILE APETALA (SAP) protein in the developing petals. The resulting decrease in SAP activity leads to a shortening of the cell proliferation period and reduced number of petal cells. The absence of private polymorphisms at the causal region in the selfing species suggests that the small-petal allele was captured from standing genetic variation in the ancestral out-crossing population. Petal-size variation in the current out-crossing population indicates that several small-effect mutations have contributed to reduce petal-size. These data demonstrate how tissue-specific regulatory elements in pleiotropic genes contribute to organ-specific evolution. In addition, they provide a plausible evolutionary explanation for the rapid evolution of flower size after the out-breeding-to-selfing transition based on additive effects of segregating alleles.

Association mapping reveals the role of purifying selection in the maintenance of genomic variation in gene expression.

The evolutionary forces that maintain genetic variation in quantitative traits within populations remain poorly understood. One hypothesis suggests that variation is under purifying selection, resulting in an excess of low-frequency variants and a negative correlation between minor allele frequency and selection coefficients. Here, we test these predictions using the genetic loci associated with total expression variation (eQTLs) and allele-specific expression variation (aseQTLs) mapped within a single population of the plant Capsella grandiflora. In addition to finding eQTLs and aseQTLs for a large fraction of genes, we show that alleles at these loci are rarer than expected and exhibit a negative correlation between phenotypic effect size and frequency. Overall, our results show that the distribution of frequencies and effect sizes of the loci responsible for local expression variation within a single outcrossing population are consistent with the effects of purifying selection.

Spatially and temporally varying selection on intrapopulation quantitative trait loci for a life history trade-off in Mimulus guttatus.

Why do populations remain genetically variable despite strong continuous natural selection? Mutation reconstitutes variation eliminated by selection and genetic drift, but theoretical and experimental studies each suggest that mutation-selection balance insufficient to explain extant genetic variation in most complex traits. The alternative hypothesis of balancing selection, wherein selection maintains genetic variation, is an aggregate of multiple mechanisms (spatial and temporal heterogeneity in selection, frequency-dependent selection, antagonistic pleiotropy, etc.). Most of these mechanisms have been demonstrated for Mendelian traits, but there is little comparable data for loci affecting quantitative characters. Here, we report a 3-year field study of selection on intrapopulation quantitative trait loci (QTL) of flower size, a highly polygenic trait in Mimulus guttatus. The QTL exhibit antagonistic pleiotropy: alleles that increase flower size, reduce viability, but increase fecundity. The magnitude and direction of selection fluctuates yearly and on a spatial scale of metres. This study provides direct evidence of balancing selection mechanisms on QTL of an ecologically relevant trait.

Explaining the heritability of an ecologically significant trait in terms of individual quantitative trait loci.

Most natural populations display substantial genetic variation in behaviour, morphology, physiology, life history and the susceptibility to disease. A major challenge is to determine the contributions of individual loci to variation in complex traits. Quantitative trait locus (QTL) mapping has identified genomic regions affecting ecologically significant traits of many species. In nearly all cases, however, the importance of these QTLs to population variation remains unclear. In this paper, we apply a novel experimental method to parse the genetic variance of floral traits of the annual plant Mimulus guttatus into contributions of individual QTLs. We first use QTL-mapping to identify nine loci and then conduct a population-based breeding experiment to estimate V(Q), the genetic variance attributable to each QTL. We find that three QTLs with moderate effects explain up to one-third of the genetic variance in the natural population. Variation at these loci is probably maintained by some form of balancing selection. Notably, the largest effect QTLs were relatively minor in their contribution to heritability.

Contribution of chromosomal polymorphisms to the G-matrix of Mimulus guttatus.

Evolution of genetic (co)variances (the G-matrix) fundamentally influences multitrait divergence. Here, we isolated the contribution of two chromosomal quantitative trait loci (QTLs), a meiotic drive locus and a polymorphic inversion, to the overall G-matrix for a suite of floral, phenological and male fitness traits in a population of Mimulus guttatus. This allowed us to predict the evolution of trait means and genetic (co)variances as a function of allele frequencies, and to evaluate theories about the maintenance of genetic variation in fitness. Individuals generated using a replicated F(2) breeding design were grown under common conditions, genotyped and measured for trait values. Significant additive genetic variance existed for all traits, and most genetic covariances were significantly nonzero. Both QTLs contribute to the additive genetic (co)variances of multiple traits. Pleiotropy was not generally consistent, either between QTLs or with the genetic background. Shifts in allele frequencies at either QTL are predicted to result in substantial changes in the G-matrix. Both QTLs contribute substantially to the genetic variation in pollen viability. The Drive QTL, and perhaps also the inversion, demonstrates the contribution of balancing selection to the maintenance of genetic variation in fitness.

Mutational disruption of the ABCC2 gene in fall armyworm, Spodoptera frugiperda, confers resistance to the Cry1Fa and Cry1A.105 insecticidal proteins.

The use of Bt proteins in crops has revolutionized insect pest management by offering effective season-long control. However, field-evolved resistance to Bt proteins threatens their utility and durability. A recent example is field-evolved resistance to Cry1Fa and Cry1A.105 in fall armyworm (Spodoptera frugiperda). This resistance has been detected in Puerto Rico, mainland USA, and Brazil. A S. frugiperda population with suspected resistance to Cry1Fa was sampled from a maize field in Puerto Rico and used to develop a resistant lab colony. The colony demonstrated resistance to Cry1Fa and partial cross-resistance to Cry1A.105 in diet bioassays. Using genetic crosses and proteomics, we show that this resistance is due to loss-of-function mutations in the ABCC2 gene. We characterize two novel mutant alleles from Puerto Rico. We also find that these alleles are absent in a broad screen of partially resistant Brazilian populations. These findings confirm that ABCC2 is a receptor for Cry1Fa and Cry1A.105 in S. frugiperda, and lay the groundwork for genetically enabled resistance management in this species, with the caution that there may be several distinct ABCC2 resistances alleles in nature.

A Segregating Inversion Generates Fitness Variation in Yellow Monkeyflower (Mimulus guttatus).

Polymorphic chromosomal rearrangements can bind hundreds of genes into single genetic loci with diverse effects. Rearrangements are often associated with local adaptation and speciation and may also be an important component of genetic variation within populations. We genetically and phenotypically characterize a segregating inversion (inv6) in the Iron Mountain (IM) population of Mimulus guttatus (yellow monkeyflower). We initially mapped inv6 as a region of recombination suppression in three F2 populations resulting from crosses among IM plants. In each case, the F1 parent was heterozygous for a derived haplotype, homogenous across markers spanning over 5 Mb of chromsome 6. In the three F2 populations, inv6 reduced male and female fitness components. In addition,i nv6 carriers suffered an ∼30% loss of pollen viability in the field. Despite these costs, inv6 exists at moderate frequency (∼8%) in the natural population, suggesting counterbalancing fitness benefits that maintain the polymorphism. Across 4 years of monitoring in the field, inv6 had an overall significant positive effect on seed production (lifetime female fitness) of carriers. This benefit was particularly strong in harsh years and may be mediated (in part) by strong positive effects on flower production. These data suggest that opposing fitness effects maintain an intermediate frequency, and as a consequence, inv6 generates inbreeding depression and high genetic variance. We discuss these findings in relation to the theory of inbreeding depression and the maintenance of fitness variation.

Divergent sorting of a balanced ancestral polymorphism underlies the establishment of gene-flow barriers in Capsella.

In the Bateson-Dobzhansky-Muller model of genetic incompatibilities post-zygotic gene-flow barriers arise by fixation of novel alleles at interacting loci in separated populations. Many such incompatibilities are polymorphic in plants, implying an important role for genetic drift or balancing selection in their origin and evolution. Here we show that NPR1 and RPP5 loci cause a genetic incompatibility between the incipient species Capsella grandiflora and C. rubella, and the more distantly related C. rubella and C. orientalis. The incompatible RPP5 allele results from a mutation in C. rubella, while the incompatible NPR1 allele is frequent in the ancestral C. grandiflora. Compatible and incompatible NPR1 haplotypes are maintained by balancing selection in C. grandiflora, and were divergently sorted into the derived C. rubella and C. orientalis. Thus, by maintaining differentiated alleles at high frequencies, balancing selection on ancestral polymorphisms can facilitate establishing gene-flow barriers between derived populations through lineage sorting of the alternative alleles.

Identifying the genes underlying quantitative traits: a rationale for the QTN programme.

The goal of identifying the genes or even nucleotides underlying quantitative and adaptive traits has been characterized as the 'QTN programme' and has recently come under severe criticism. Part of the reason for this criticism is that much of the QTN programme has asserted that finding the genes and nucleotides for adaptive and quantitative traits is a fundamental goal, without explaining why it is such a hallowed goal. Here we outline motivations for the QTN programme that offer general insight, regardless of whether QTNs are of large or small effect, and that aid our understanding of the mechanistic dynamics of adaptive evolution. We focus on five areas: (i) vertical integration of insight across different levels of biological organization, (ii) genetic parallelism and the role of pleiotropy in shaping evolutionary dynamics, (iii) understanding the forces maintaining genetic variation in populations, (iv) distinguishing between adaptation from standing variation and new mutation, and (v) the role of genomic architecture in facilitating adaptation. We argue that rather than abandoning the QTN programme, we should refocus our efforts on topics where molecular data will be the most effective for testing hypotheses about phenotypic evolution.

Repeated evolutionary changes of leaf morphology caused by mutations to a homeobox gene.

Elucidating the genetic basis of morphological changes in evolution remains a major challenge in biology. Repeated independent trait changes are of particular interest because they can indicate adaptation in different lineages or genetic and developmental constraints on generating morphological variation. In animals, changes to "hot spot" genes with minimal pleiotropy and large phenotypic effects underlie many cases of repeated morphological transitions. By contrast, only few such genes have been identified from plants, limiting cross-kingdom comparisons of the principles of morphological evolution. Here, we demonstrate that the REDUCED COMPLEXITY (RCO) locus underlies more than one naturally evolved change in leaf shape in the Brassicaceae. We show that the difference in leaf margin dissection between the sister species Capsella rubella and Capsella grandiflora is caused by cis-regulatory variation in the homeobox gene RCO-A, which alters its activity in the developing lobes of the leaf. Population genetic analyses in the ancestral C. grandiflora indicate that the more-active C. rubella haplotype is derived from a now rare or lost C. grandiflora haplotype via additional mutations. In Arabidopsis thaliana, the deletion of the RCO-A and RCO-B genes has contributed to its evolutionarily derived smooth leaf margin, suggesting the RCO locus as a candidate for an evolutionary hot spot. We also find that temperature-responsive expression of RCO-A can explain the phenotypic plasticity of leaf shape to ambient temperature in Capsella, suggesting a molecular basis for the well-known negative correlation between temperature and leaf margin dissection.

A case of paramyotonia congenita without periodic paralysis: electrophysiological and molecular genetic studies.

Paramyotonia congenita (PC), first described in 1886 by Eulenberg, is characterized by cold and exercise-induced muscle stiffness and intermittent flaccid paresis not necessarily related to cold or myotonia. Several authors segregated a pure form of PC, which has no periodic paralysis, even after cold exposure. The existence of this phenotype has been debated in the literature. We describe electrophysiological and molecular genetic features of a patient with PC who had no history of periodic paralysis. Immersion in cold water or potassium load could not induce clinical paralysis. However, repetitive nerve stimulation and exercise test demonstrated a drop in compound muscle action potential amplitude. Genetic analysis revealed the substitution of valine for glycine on the human skeletal muscle sodium channel (SCN4A) gene. The G1306V mutation is rare in the classic form of PC, and moreover might be the first in pure paramyotonia.

Distinct light and clock modulation of cytosolic free Ca2+ oscillations and rhythmic CHLOROPHYLL A/B BINDING PROTEIN2 promoter activity in Arabidopsis.

Plants have circadian oscillations in the concentration of cytosolic free calcium ([Ca(2+)](cyt)). To dissect the circadian Ca(2+)-signaling network, we monitored circadian [Ca(2+)](cyt) oscillations under various light/dark conditions (including different spectra) in Arabidopsis thaliana wild type and photoreceptor and circadian clock mutants. Both red and blue light regulate circadian oscillations of [Ca(2+)](cyt). Red light signaling is mediated by PHYTOCHROME B (PHYB). Blue light signaling occurs through the redundant action of CRYPTOCHROME1 (CRY1) and CRY2. Blue light also increases the basal level of [Ca(2+)](cyt), and this response requires PHYB, CRY1, and CRY2. Light input into the oscillator controlling [Ca(2+)](cyt) rhythms is gated by EARLY FLOWERING3. Signals generated in the dark also regulate the circadian behavior of [Ca(2+)](cyt). Oscillations of [Ca(2+)](cyt) and CHLOROPHYLL A/B BINDING PROTEIN2 (CAB2) promoter activity are dependent on the rhythmic expression of LATE ELONGATED HYPOCOTYL and CIRCADIAN CLOCK-ASSOCIATED1, but [Ca(2+)](cyt) and CAB2 promoter activity are uncoupled in the timing of cab1 (toc1-1) mutant but not in toc1-2. We suggest that the circadian oscillations of [Ca(2+)](cyt) and CAB2 promoter activity are regulated by distinct oscillators with similar components that are used in a different manner and that these oscillators may be located in different cell types in Arabidopsis.