Genetic architecture of life history traits and environment-specific trade-offs

Local climate and vernalization sensitivity predict the latitudinal patterns of flowering onset in the crop wild relative Linum bienne Mill

BACKGROUND AND AIMS

The timing of flowering onset is often correlated with latitude, indicative of climatic gradients. Flowering onset in temperate species commonly requires exposure to cold temperatures, known as vernalization. Hence, population differentiation of flowering onset with latitude might reflect adaptation to the local climatic conditions experienced by populations.

METHODS

Within its western range, seeds from Linum bienne populations (the wild relative of cultivated Linum usitatissimum) were used to describe the latitudinal differentiation of flowering onset to determine its association with the local climate of the population. A vernalization experiment including different crop cultivars was used to determine how vernalization accelerates flowering onset, in addition to the vernalization sensitivity response among populations and cultivars. Additionally, genetic differentiation of L. bienne populations along the latitudinal range was scrutinized using microsatellite markers.

KEY RESULTS

Flowering onset varied with latitude of origin, with southern populations flowering earlier than their northern counterparts. Vernalization reduced the number of days to flowering onset, but vernalization sensitivity was greater in northern populations compared with southern ones. Conversely, vernalization delayed flowering onset in the crop, exhibiting less variation in sensitivity. In L. bienne, both flowering onset and vernalization sensitivity were better predicted by the local climate of the population than by latitude itself. Microsatellite data unveiled genetic differentiation of populations, forming two groups geographically partitioned along latitude.

CONCLUSIONS

The consistent finding of latitudinal variation across experiments suggests that both flowering onset and vernalization sensitivity in L. bienne populations are under genetic regulation and might depend on climatic cues at the place of origin. The association with climatic gradients along latitude suggests that the climate experienced locally drives population differentiation of the flowering onset and vernalization sensitivity patterns. The genetic population structure suggests that past population history could have influenced the flowering initiation patterns detected, which deserves further work.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

The genetic architecture and evolution of life-history divergence among perennials in the Mimulus guttatus species complex

Ecological divergence is a fundamental source of phenotypic diversity between closely related species, yet the genetic architecture of most ecologically relevant traits is poorly understood. Differences in elevation can impose substantial divergent selection on both complex, correlated suites of traits (such as life-history), as well as novel adaptations. We use the Mimulus guttatus species complex to assess if the divergence in elevation is accompanied by trait divergence in a group of closely related perennials and determine the genetic architecture of this divergence. We find that divergence in elevation is associated with differences in life-history, as well as a unique trait, the production of rhizomes. The divergence between two perennials is largely explained by few mid-to-large effect quantitative trait loci (QTLs). However, the presence of QTLs with correlated, but opposing effects on multiple traits leads to some hybrids with transgressive trait combinations. Lastly, we find that the genetic architecture of the ability to produce rhizomes changes through development, wherein most hybrids produce rhizomes, but only later in development. Our results suggest that elevational differences may shape life-history divergence between perennials, but aspects of the genetic architecture of divergence may have implications for hybrid fitness in nature.

Experimental evolution of adaptive divergence under varying degrees of gene flow

Adaptive divergence is the key evolutionary process generating biodiversity by means of natural selection. Yet, the conditions under which it can arise in the presence of gene flow remain contentious. To address this question, we subjected 132 sexually reproducing fission yeast populations, sourced from two independent genetic backgrounds, to disruptive ecological selection and manipulated the level of migration between environments. Contrary to theoretical expectations, adaptive divergence was most pronounced when migration was either absent (allopatry) or maximal (sympatry), but was much reduced at intermediate rates (parapatry and local mating). This effect was apparent across central life-history components (survival, asexual growth and mating) but differed in magnitude between ancestral genetic backgrounds. The evolution of some fitness components was constrained by pervasive negative correlations (trade-off between asexual growth and mating), while others changed direction under the influence of migration (for example, survival and mating). In allopatry, adaptive divergence was mainly conferred by standing genetic variation and resulted in ecological specialization. In sympatry, divergence was mainly mediated by novel mutations enriched in a subset of genes and was characterized by the repeated emergence of two strategies: an ecological generalist and an asexual growth specialist. Multiple loci showed consistent evidence for antagonistic pleiotropy across migration treatments providing a conceptual link between adaptation and divergence. This evolve-and-resequence experiment shows that rapid ecological differentiation can arise even under high rates of gene flow. It further highlights that adaptive trajectories are governed by complex interactions of gene flow, ancestral variation and genetic correlations.

Asexual reproduction and growth rate: independent and plastic life history traits in Neurospora crassa

Trade-offs among traits influencing fitness are predicted by life history theory because resources allocated to one function are unavailable to another. Here we examine the relationship between two such traits, asexual reproduction and growth rate, in the filamentous fungus Neurospora crassa, where shared genetic and physiological factors and a source-sink energetic relationship between growth and reproduction may constrain the evolution of these traits. To test growth-reproduction relationships in this species, we independently selected on mycelial growth rate or asexual spore production in a heterogeneous lab-derived population and evaluated the response of the non-selected traits. Combined with phenotypes for the 20 wild strains used to produce the heterogeneous population and the genome-wide genotypes of 468 strains, these data show that growth and reproduction are highly plastic in N. crassa and do not trade off either among wild strains or after laboratory selection in two environments. Rather, we find no predictable growth-reproduction relationship in the environments tested, indicating an effective absence of genetic constraint between these traits. Our results suggest that growth rate and asexual reproduction may not respond predictably to environmental change and suggest that reliance on a single trait as a proxy for fitness in fungal studies may be inadvisable.

Inflorescence photosynthetic contribution to fitness releases Arabidopsis thaliana plants from trade-off constraints on early flowering

Leaves are thought to be the primary carbon source for reproduction in plants, so a positive relationship between vegetative size and reproductive output is expected, establishing a trade-off between time to reproduction and reproductive output. A common response to higher temperatures due to climate changes is the induction of earlier transition into reproduction. Thus, in annual plants, earlier transition into flowering can potentially constrain plant size and reduce seed production. However, trade-offs between early reproduction and fitness are not always observed, suggesting mechanisms to escape the constraints of early flowering do exist. Here, we test whether inflorescence photosynthesis contribution to the reproductive output of Arabidopsis thaliana can offset the cost of early reproduction. We followed the development, growth rate and fitness of 15 accessions, and removed all rosette leaves at flowering (prior to the completion of inflorescence development or any fruit production) in half of the plants to determine the ability of inflorescences to maintain fitness in the absence of leaves. Although leaf removal significantly reduced fruit number, seed weight and plant height, even the most severely impacted accessions maintained 35% of their fitness with the inflorescence as the sole photosynthetic organ; and some accessions experienced no reduction in fitness. Differences between accessions in their ability to maintain fitness after leaf removal is best explained by earlier flowering time and the ability to maintain as many or more branches after leaf removal as in the control treatment. Although earlier flowering does constrain plant vegetative size, we found that inflorescence photosynthesis can significantly contribute to seed production, explaining why early flowering plants can maintain high fitness despite a reduction in vegetative size. Thus, plants can be released from the usually assumed trade-offs associated with earlier reproduction, and selection on inflorescence traits can mediate the impact of climate change on phenology.

Allocation to male vs female floral function varies by currency and responds differentially to density and moisture stress

Allocation of finite resources to separate reproductive functions is predicted to vary across environments and affect fitness. Biomass is the most commonly measured allocation currency; however, in comparison with nutrients it may be less limited and express different environmental and evolutionary responses. Here, we measured carbon, nitrogen, phosphorus, and biomass allocation among floral whorls in recombinant inbred lines of Brassica rapa in multiple environments to characterize the genetic architecture of floral allocation, including its sensitivity to environmental heterogeneity and to choice of currency. Mass, carbon, and nitrogen allocation to female whorls (pistils and sepals) decreased under high density, whereas nitrogen allocation to male organs (stamens) decreased under drought. Phosphorus allocation decreased by half in pistils under drought, while stamen phosphorus was unaffected by environment. While the contents of each currency were positively correlated among whorls, selection to improve fitness through female (or male) function typically favored increased allocation to pistils (or stamens) but decreased allocation to other whorls. Finally, genomic regions underlying correlations among allocation metrics were mapped, and loci related to nitrogen uptake and floral organ development were located within mapped quantitative trait loci. Our candidate gene identification suggests that nutrient uptake may be a limiting step in maintaining male allocation. Taken together, allocation to male vs female function is sensitive to distinct environmental stresses, and the choice of currency affects the interpretation of floral allocation responses to the environment. Further, genetic correlations may counter the evolution of allocation patterns that optimize fitness through female or male function.

Natural variation in life history strategy of Arabidopsis thaliana determines stress responses to drought and insects of different feeding guilds

Plants are sessile organisms and, consequently, are exposed to a plethora of stresses in their local habitat. As a result, different populations of a species are subject to different selection pressures leading to adaptation to local conditions and intraspecific divergence. The annual brassicaceous plant Arabidopsis thaliana is an attractive model for ecologists and evolutionary biologists due to the availability of a large collection of resequenced natural accessions. Accessions of A. thaliana display one of two different life cycle strategies: summer and winter annuals. We exposed a collection of 308 European Arabidopsis accessions, that have been genotyped for 250K SNPs, to a range of stresses: one abiotic stress (drought), four biotic stresses (Pieris rapae caterpillars, Plutella xylostella caterpillars, Frankliniella occidentalis thrips and Myzus persicae aphids) and two combined stresses (drought plus P. rapae and Botrytis cinerea fungus plus P. rapae). We identified heritable genetic variation for responses to the different stresses, estimated by narrow-sense heritability. We found that accessions displaying different life cycle strategies differ in their response to stresses. Winter annuals are more resistant to drought, aphids and thrips and summer annuals are more resistant to P. rapae and P. xylostella caterpillars. Summer annuals are also more resistant to the combined stresses of drought plus P. rapae and infection by the fungus Botryris cinerea plus herbivory by P. rapae. Adaptation to drought displayed a longitudinal gradient. Finally, trade-offs were recorded between the response to drought and responses to herbivory by caterpillars of the specialist herbivore P. rapae.

Environmental effects on the structure of the G-matrix

Genetic correlations between traits determine the multivariate response to selection in the short term, and thereby play a causal role in evolutionary change. Although individual studies have documented environmentally induced changes in genetic correlations, the nature and extent of environmental effects on multivariate genetic architecture across species and environments remain largely uncharacterized. We reviewed the literature for estimates of the genetic variance-covariance (G) matrix in multiple environments, and compared differences in G between environments to the divergence in G between conspecific populations (measured in a common garden). We found that the predicted evolutionary trajectory differed as strongly between environments as it did between populations. Between-environment differences in the underlying structure of G (total genetic variance and the relative magnitude and orientation of genetic correlations) were equal to or greater than between-population differences. Neither environmental novelty, nor the difference in mean phenotype predicted these differences in G. Our results suggest that environmental effects on multivariate genetic architecture may be comparable to the divergence that accumulates over dozens or hundreds of generations between populations. We outline avenues of future research to address the limitations of existing data and characterize the extent to which lability in genetic correlations shapes evolution in changing environments.

Genotypic variation in biomass allocation in response to field drought has a greater affect on yield than gas exchange or phenology

BACKGROUND

Plant performance in agricultural and natural settings varies with moisture availability, and understanding the range of potential drought responses and the underlying genetic architecture is important for understanding how plants will respond to both natural and artificial selection in various water regimes. Here, we raised genotypes of Brassica rapa under well-watered and drought treatments in the field. Our primary goal was to understand the genetic architecture and yield effects of different drought-escape and dehydration-avoidance strategies.

RESULTS

Drought treatments reduced soil moisture by 62 % of field capacity. Drought decreased biomass accumulation and fruit production by as much as 48 %, whereas instantaneous water-use efficiency and root:shoot ratio increased. Genotypes differed in the mean value of all traits and in the sensitivity of biomass accumulation, root:shoot ratio, and fruit production to drought. Bivariate correlations involving gas-exchange and phenology were largely constant across environments, whereas those involving root:shoot varied across treatments. Although root:shoot was typically unrelated to gas-exchange or yield under well-watered conditions, genotypes with low to moderate increases in root:shoot allocation in response to drought survived the growing season, maintained maximum photosynthesis levels, and produced more fruit than genotypes with the greatest root allocation under drought. QTL for gas-exchange and yield components (total biomass or fruit production) had common effects across environments while those for root:shoot were often environment-specific.

CONCLUSIONS

Increases in root allocation beyond those needed to survive and maintain favorable water relations came at the cost of fruit production. The environment-specific effects of root:shoot ratio on yield and the differential expression of QTL for this trait across water regimes have important implications for efforts to improve crops for drought resistance.

Paralogous SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes differentially regulate leaf initiation and reproductive phase change in petunia

MAIN CONCLUSION

Duplicated petunia clade-VI SPL genes differentially promote the timing of inflorescence and flower development, and leaf initiation rate. The timing of plant reproduction relative to favorable environmental conditions is a critical component of plant fitness, and is often associated with variation in plant architecture and habit. Recent studies have shown that overexpression of the microRNA miR156 in distantly related annual species results in plants with perennial characteristics, including late flowering, weak apical dominance, and abundant leaf production. These phenotypes are largely mediated through the negative regulation of a subset of genes belonging to the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) family of transcription factors. In order to determine how and to what extent paralogous SPL genes have partitioned their roles in plant growth and development, we functionally characterized petunia clade-VI SPL genes under different environmental conditions. Our results demonstrate that PhSBP1and PhSBP2 differentially promote discrete stages of the reproductive transition, and that PhSBP1, and possibly PhCNR, accelerates leaf initiation rate. In contrast to the closest homologs in annual Arabidopsis thaliana and Mimulus guttatus, PhSBP1 and PhSBP2 transcription is not mediated by the gibberellic acid pathway, but is positively correlated with photoperiod and developmental age. The developmental functions of clade-VI SPL genes have, thus, evolved following both gene duplication and speciation within the core eudicots, likely through differential regulation and incomplete sub-functionalization.

Timing of shoot development transitions affects degree of perenniality in Arabidopsis lyrata (Brassicaceae)

BACKGROUND

Perenniality is best understood in quantitative terms, involving the relationship between production vs. turnover of meristems, biomass, or energy reserves. Previous quantitative trait locus (QTL) studies using divergent populations of the perennial rock cress Arabidopsis lyrata have shown that trade-offs in vegetative growth vs. reproduction are due to cascading effects of differences in early vegetative development, which contribute to local adaptation. However, details of the developmental differences and how they affect perenniality remained unclear. In this study, we investigated in detail the developmental differences in perenniality between populations. A. lyrata from Norway and North Carolina populations, representing contrasting environments and degrees of perenniality, were grown under controlled conditions, and data were collected on plant phenology and shoot-level development. We tested hypotheses that differences in perenniality involve strict allocation of lateral meristems to vegetative vs. reproductive fates, or alternatively quantitative effects of pre-reproductive vegetative development.

RESULTS

The two populations showed large differences in the degree of vegetative development on individual shoots prior to reproductive transitions. The number of leaves produced on shoots prior to bolting, and not strict meristem allocation or variation in apical dominance, was able to explain variation in the number of inflorescences on individual plants. These results suggested that allocation of time to shoot vegetative vs. reproductive development could be a major factor in resource allocation differences between the populations.

CONCLUSIONS

Based on these results and those of previous QTL studies, we propose a model in which the degree of shoot vegetative development shapes the developmental context for reproduction and subsequent vegetative growth in different environments. Climate-specific effects of shoot development patterns on reproductive output and survival may result in divergent evolutionary trajectories along a perenniality continuum, which may have broader relevance for plant life history evolution.

Increased plant productivity and decreased microbial respiratory C loss by plant growth-promoting rhizobacteria under elevated CO₂

Increased plant productivity and decreased microbial respiratory C loss can potentially mitigate increasing atmospheric CO₂, but we currently lack effective means to achieve these goals. Soil microbes may play critical roles in mediating plant productivity and soil C/N dynamics under future climate scenarios of elevated CO₂ (eCO₂) through optimizing functioning of the root-soil interface. By using a labeling technique with (13)C and (15)N, we examined the effects of plant growth-promoting Pseudomonas fluorescens on C and N cycling in the rhizosphere of a common grass species under eCO₂. These microbial inoculants were shown to increase plant productivity. Although strong competition for N between the plant and soil microbes was observed, the plant can increase its capacity to store more biomass C per unit of N under P. fluorescens addition. Unlike eCO₂ effects, P. fluorescens inoculants did not change mass-specific microbial respiration and accelerate soil decomposition related to N cycling, suggesting these microbial inoculants mitigated positive feedbacks of soil microbial decomposition to eCO₂. The potential to mitigate climate change by optimizing soil microbial functioning by plant growth-promoting Pseudomonas fluorescens is a prospect for ecosystem management.

Multivariate genetic analysis of plant responses to water deficit and high temperature revealed contrasting adaptive strategies

How genetic factors control plant performance under stressful environmental conditions is a central question in ecology and for crop breeding. A multivariate framework was developed to examine the genetic architecture of performance-related traits in response to interacting environmental stresses. Ecophysiological and life history traits were quantified in the Arabidopsis thaliana Ler × Cvi mapping population exposed to constant soil water deficit and high air temperature. The plasticity of the genetic variance-covariance matrix (G-matrix) was examined using mixed-effects models after regression into principal components. Quantitative trait locus (QTL) analysis was performed on the predictors of genotype effects and genotype by environment interactions (G × E). Three QTLs previously identified for flowering time had antagonistic G × E effects on carbon acquisition and the other traits (phenology, growth, leaf morphology, and transpiration). This resulted in a size-dependent response of water use efficiency (WUE) to high temperature but not soil water deficit, indicating that most of the plasticity of carbon acquisition and WUE to temperature is controlled by the loci that control variation of development, size, growth, and transpiration. A fourth QTL, MSAT2.22, controlled the response of carbon acquisition to specific combinations of watering and temperature irrespective of plant size and development, growth, and transpiration rate, which resulted in size-independent plasticity of WUE. These findings highlight how the strategies to optimize plant performance may differ in response to water deficit and high temperature (or their combination), and how different G × E effects could be targeted to improve plant tolerance to these stresses.

Quantitative trait loci × environment interactions for plant morphology vary over ontogeny in Brassica rapa

Growth in plants occurs via the addition of repeating modules, suggesting that the genetic architecture of similar subunits may vary between earlier- and later-developing modules. These complex environment × ontogeny interactions are not well elucidated, as studies examining quantitative trait loci (QTLs) expression over ontogeny have not included multiple environments. Here, we characterized the genetic architecture of vegetative traits and onset of reproduction over ontogeny in recombinant inbred lines of Brassica rapa in the field and glasshouse. The magnitude of genetic variation in plasticity of seedling internodes was greater than in those produced later in ontogeny. We correspondingly detected that QTLs for seedling internode length were environment-specific, whereas later in ontogeny the majority of QTLs affected internode lengths in all treatments. The relationship between internode traits and onset of reproduction varied with environment and ontogenetic stage. This relationship was observed only in the glasshouse environment and was largely attributable to one environment-specific QTL. Our results provide the first evidence of a QTL × environment × ontogeny interaction, and provide QTL resolution for differences between early- and later-stage plasticity for stem elongation. These results also suggest potential constraints on morphological evolution in early vs later modules as a result of associations with reproductive timing.

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment-specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2-5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.

Quantitative variation in water-use efficiency across water regimes and its relationship with circadian, vegetative, reproductive, and leaf gas-exchange traits

Drought limits light harvesting, resulting in lower plant growth and reproduction. One trait important for plant drought response is water-use efficiency (WUE). We investigated (1) how the joint genetic architecture of WUE, reproductive characters, and vegetative traits changed across drought and well-watered conditions, (2) whether traits with distinct developmental bases (e.g. leaf gas exchange versus reproduction) differed in the environmental sensitivity of their genetic architecture, and (3) whether quantitative variation in circadian period was related to drought response in Brassica rapa. Overall, WUE increased in drought, primarily because stomatal conductance, and thus water loss, declined more than carbon fixation. Genotypes with the highest WUE in drought expressed the lowest WUE in well-watered conditions, and had the largest vegetative and floral organs in both treatments. Thus, large changes in WUE enabled some genotypes to approach vegetative and reproductive trait optima across environments. The genetic architecture differed for gas-exchange and vegetative traits across drought and well-watered conditions, but not for floral traits. Correlations between circadian and leaf gas-exchange traits were significant but did not vary across treatments, indicating that circadian period affects physiological function regardless of water availability. These results suggest that WUE is important for drought tolerance in Brassica rapa and that artificial selection for increased WUE in drought will not result in maladaptive expression of other traits that are correlated with WUE.