Flowering time QTL in natural populations ofArabidopsis thalianaand implications for their adaptive value

Planting Genomes in the Wild: Arabidopsis from Genetics History to the Ecology and Evolutionary Genomics Era

The genetics model system Arabidopsis thaliana (L.) Heynh. lives across a vast geographic range with contrasting climates, in response to which it has evolved diverse life histories and phenotypic adaptations. In the last decade, the cataloging of worldwide populations, DNA sequencing of whole genomes, and conducting of outdoor field experiments have transformed it into a powerful evolutionary ecology system to understand the genomic basis of adaptation. Here, we summarize new insights on Arabidopsis following the coordinated efforts of the 1001 Genomes Project, the latest reconstruction of biogeographic and demographic history, and the systematic genomic mapping of trait natural variation through 15 years of genome-wide association studies. We then put this in the context of local adaptation across climates by summarizing insights from 73 Arabidopsis outdoor common garden experiments conducted to date. We conclude by highlighting how molecular and genomic knowledge of adaptation can help us to understand species' (mal)adaptation under ongoing climate change.

Dual-trait genomic analysis in highly stratified Arabidopsis thaliana populations using genome-wide association summary statistics

Genome-wide association study (GWAS) is a powerful tool to identify genomic loci underlying complex traits. However, the application in natural populations comes with challenges, especially power loss due to population stratification. Here, we introduce a bivariate analysis approach to a GWAS dataset of Arabidopsis thaliana. We demonstrate the efficiency of dual-phenotype analysis to uncover hidden genetic loci masked by population structure via a series of simulations. In real data analysis, a common allele, strongly confounded with population structure, is discovered to be associated with late flowering and slow maturation of the plant. The discovered genetic effect on flowering time is further replicated in independent datasets. Using Mendelian randomization analysis based on summary statistics from our GWAS and expression QTL scans, we predicted and replicated a candidate gene AT1G11560 that potentially causes this association. Further analysis indicates that this locus is co-selected with flowering-time-related genes. The discovered pleiotropic genotype-phenotype map provides new insights into understanding the genetic correlation of complex traits.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

Genome-wide association study identifies a novel BMI1A QTL allele that confers FLC expression diversity in Arabidopsis thaliana

Identification and understanding of the genetic basis of natural variations in plants are essential for comprehending their phenotypic adaptation. Here, we report a genome-wide association study (GWAS) of FLOWERING LOCUS C (FLC) expression in 727 Arabidopsis accessions. We identified B LYMPHOMA MOLONEY MURINE LEUKEMIA VIRUS INSERTION REGION 1 HOMOLOG 1A (BMI1A) as a causal gene for one of the FLC expression quantitative trait loci (QTLs). Loss of function in BMI1A increases FLC expression and delays flowering time at 16 °C significantly compared with the wild type (Col-0). BMI1A activity is required for histone H3 lysine 27 trimethylation (H3K27me3) accumulation at the FLC, MADS AFFECTING FLOWERING 4 (MAF4), and MAF5 loci at low ambient temperature. We further uncovered two BMI1A haplotypes associated with the natural variation in FLC expression and flowering time at 16 °C, and demonstrated that polymorphisms in the BMI1A promoter region are the main contributor. Different BMI1A haplotypes are strongly associated with geographical distribution, and the low ambient temperature-sensitive BMI1A variants are associated with a lower mean temperature of the driest quarter of their collection sites compared with the temperature-non-responsive variants, indicating that the natural variations in BMI1A have adaptive functions in FLC expression and flowering time regulation. Therefore, our results provide new insights into the natural variations in FLC expression and flowering time diversity in plants.

Molecular Evaluation of the Effects of FLC Homologs and Coordinating Regulators on the Flowering Responses to Vernalization in Cabbage (Brassica oleracea var. capitata) Genotypes

The flowering loci of cabbage must be understood to boost their productivity. In this study, to clarify the flowering mechanisms of cabbage, we examined the three flowering repressors BoFLC1, 2 and 3, and the flowering regulators BoGI, BoCOOLAIR, and BoVIN3 of early (CAB1), middle (CAB3), and late (CAB5) flowering cabbage genotypes. Analysis of allele-specifically amplified genomic DNA and various sequence alignments demonstrated that maximal insertions and deletions influenced cabbage flowering behavior, notably in CAB3 and CAB5. Phylogenetic studies showed that BoFLC1, 2, and 3 in the CAB1, 3, and 5 genotypes had the highest homologies to other Brassica species, with CAB3 and 5 the most similar. Although CAB3 and CAB5 have comparable genetic patterns, flowering repressors and flowering regulators were investigated individually with and without vernalization to determine their minor flowering differences. The expression investigation revealed that vernalized CAB5 downregulated all BoFLC genes compared to CAB3 and, in contrast, CAB3 exhibited upregulated BoCOOLAIR. We hypothesized that the CAB3 BoFLC locus' additional insertions may have led to BoCOOLAIR overexpression and BoFLC downregulation. This study sheds light on cabbage genotypes-particularly those of CAB1 and CAB5-and suggests that structural variations in BoFLC2 and 3 bind flowering regulators, such as COOLAIR, which may affect cabbage flowering time.

Responsiveness to long days for flowering is reduced in Arabidopsis by yearly variation in growing season temperatures

Conservative flowering behaviours, such as flowering during long days in summer or late flowering at a high leaf number, are often proposed to protect against variable winter and spring temperatures which lead to frost damage if premature flowering occurs. Yet, due the many factors in natural environments relative to the number of individuals compared, assessing which climate characteristics drive these flowering traits has been difficult. We applied a multidisciplinary approach to 10 winter-annual Arabidopsis thaliana populations from a wide climactic gradient in Norway. We used a variable reduction strategy to assess which of 100 climate descriptors from their home sites correlated most to their flowering behaviours when tested for responsiveness to photoperiod after saturation of vernalization; then, assessed sequence variation of 19 known environmental-response flowering genes. Photoperiod responsiveness inversely correlated with interannual variation in timing of growing season onset. Time to flowering appeared driven by growing season length, curtailed by cold fall temperatures. The distribution of FLM, TFL2 and HOS1 haplotypes, genes involved in ambient temperature response, correlated with growing-season climate. We show that long-day responsiveness and late flowering may be driven not by risk of spring frosts, but by growing season temperature and length, perhaps to opportunistically maximize growth.

Genetic Architecture of Flowering Time Differs Between Populations With Contrasting Demographic and Selective Histories

Understanding the evolutionary factors that impact the genetic architecture of traits is a central goal of evolutionary genetics. Here, we investigate how quantitative trait variation accumulated over time in populations that colonized a novel environment. We compare the genetic architecture of flowering time in Arabidopsis populations from the drought-prone Cape Verde Islands and their closest outgroup population from North Africa. We find that trait polygenicity is severely reduced in the island populations compared to the continental North African population. Further, trait architectures and reconstructed allelic histories best fit a model of strong directional selection in the islands in accord with a Fisher-Orr adaptive walk. Consistent with this, we find that large-effect variants that disrupt major flowering time genes (FRI and FLC) arose first, followed by smaller effect variants, including ATX2 L125F, which is associated with a 4-day reduction in flowering time. The most recently arising flowering time-associated loci are not known to be directly involved in flowering time, consistent with an omnigenic signature developing as the population approaches its trait optimum. Surprisingly, we find no effect in the natural population of EDI-Cvi-0 (CRY2 V367M), an allele for which an effect was previously validated by introgression into a Eurasian line. Instead, our results suggest the previously observed effect of the EDI-Cvi-0 allele on flowering time likely depends on genetic background, due to an epistatic interaction. Altogether, our results provide an empirical example of the effects demographic history and selection has on trait architecture.

Ecological genetics of local adaptation in Arabidopsis: An 8-year field experiment

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Genetic Mapping of Genotype-by-Ploidy Effects in Arabidopsis thaliana

Plants can express different phenotypic responses following polyploidization, but ploidy-dependent phenotypic variation has so far not been assigned to specific genetic factors. To map such effects, segregating populations at different ploidy levels are required. The availability of an efficient haploid inducer line in Arabidopsis thaliana allows for the rapid development of large populations of segregating haploid offspring. Because Arabidopsis haploids can be self-fertilised to give rise to homozygous doubled haploids, the same genotypes can be phenotyped at both the haploid and diploid ploidy level. Here, we compared the phenotypes of recombinant haploid and diploid offspring derived from a cross between two late flowering accessions to map genotype × ploidy (G × P) interactions. Ploidy-specific quantitative trait loci (QTLs) were detected at both ploidy levels. This implies that mapping power will increase when phenotypic measurements of monoploids are included in QTL analyses. A multi-trait analysis further revealed pleiotropic effects for a number of the ploidy-specific QTLs as well as opposite effects at different ploidy levels for general QTLs. Taken together, we provide evidence of genetic variation between different Arabidopsis accessions being causal for dissimilarities in phenotypic responses to altered ploidy levels, revealing a G × P effect. Additionally, by investigating a population derived from late flowering accessions, we revealed a major vernalisation-specific QTL for variation in flowering time, countering the historical bias of research in early flowering accessions.

Climate drives rhizosphere microbiome variation and divergent selection between geographically distant Arabidopsis populations

Disentangling the contribution of climatic and edaphic factors to microbiome variation and local adaptation in plants requires an experimental approach to uncouple their effects and test for causality. We used microbial inocula, soil matrices and plant genotypes derived from two natural Arabidopsis thaliana populations in northern and southern Europe in an experiment conducted in climatic chambers mimicking seasonal changes in temperature, day length and light intensity of the home sites of the two genotypes. The southern A. thaliana genotype outperformed the northern genotype in the southern climate chamber, whereas the opposite was true in the northern climate chamber. Recipient soil matrix, but not microbial composition, affected plant fitness, and effects did not differ between genotypes. Differences between chambers significantly affected rhizosphere microbiome assembly, although these effects were small in comparison with the shifts induced by physicochemical differences between soil matrices. The results suggest that differences in seasonal changes in temperature, day length and light intensity between northern and southern Europe have strongly influenced adaptive differentiation between the two A. thaliana populations, whereas effects of differences in soil factors have been weak. By contrast, below-ground differences in soil characteristics were more important than differences in climate for rhizosphere microbiome differentiation.

Effects of primary seed dormancy on lifetime fitness of Arabidopsis thaliana in the field

BACKGROUND AND AIMS

Seed dormancy determines the environmental niche of plants in seasonal environments, and has consequences for plant performance that potentially go far beyond the seed and seedling stages. In this study, we examined the cascading effects of seed dormancy on the expression of subsequent life-history traits and fitness in the annual herb Arabidopsis thaliana.

METHODS

We planted seeds of >200 recombinant inbred lines (RILs) derived from a cross between two locally adapted populations (Italy and Sweden), and both parental genotypes at the native site of the Swedish population in three consecutive years. We quantified the relationship between primary seed dormancy and the expression of subsequent life-history traits and fitness in the RIL population with path analysis. To examine the effects of differences in dormancy on the relative fitness of the two parental genotypes, we planted dormant seeds during the seed dispersal period and non-dormant seeds during the germination period of the local population.

KEY RESULTS

In the RIL population, strong primary dormancy was associated with high seedling survival, but with low adult survival and fecundity, and path analysis indicated that this could be explained by effects on germination timing, rosette size and flowering start. The relationship between primary seed dormancy and germination proportion varied among years, and this was associated with differences in seasonal changes in soil moisture. The planting of dormant and non-dormant seeds indicated that the lower primary dormancy of the local Swedish genotype contributed to its higher germination proportion in two years and to its higher fecundity in one year.

CONCLUSIONS

Our results show that seed dormancy affects trait expression and fitness components across the life cycle, and suggest that among-year variation in the incidence of drought during the germination period should be considered when predicting the consequences of climatic change for population growth and evolution.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Flowering Times of Wild Arabidopsis Accessions From Across Norway Correlate With Expression Levels of FT, CO, and FLC Genes

Temperate species often require or flower most rapidly in the long daylengths, or photoperiods, experienced in summer or after prolonged periods of cold temperatures, referred to as vernalization. Yet, even within species, plants vary in the degree of responsiveness to these cues. In Arabidopsis thaliana, CONSTANS (CO) and FLOWERING LOCUS C (FLC) genes are key to photoperiod and vernalization perception and antagonistically regulate FLOWERING LOCUS T (FT) to influence the flowering time of the plants. However, it is still an open question as to how these genes vary in their interactions among wild accessions with different flowering behaviors and adapted to different microclimates, yet this knowledge could improve our ability to predict plant responses in variable natural conditions. To assess the relationships among these genes and to flowering time, we exposed 10 winter-annual Arabidopsis accessions from throughout Norway, ranging from early to late flowering, along with two summer-annual accessions to 14 weeks of vernalization and either 8- or 19-h photoperiods to mimic Norwegian climate conditions, then assessed gene expression levels 3-, 5-, and 8-days post vernalization. CO and FLC explained both FT levels and flowering time (days) but not rosette leaf number at flowering. The correlation between FT and flowering time increased over time. Although vernalization suppresses FLC, FLC was high in the late-flowering accessions. Across accessions, FT was expressed only at low FLC levels and did not respond to CO in the late-flowering accessions. We proposed that FT may only be expressed below a threshold value of FLC and demonstrated that these three genes correlated to flowering times across genetically distinct accessions of Arabidopsis.

Genetic differences in the temporal and environmental stability of transgenerational environmental effects

Environments influence the expression of phenotypes of individuals, their progeny, and even their grandprogeny. The duration of environmental effects and how they are modified by subsequent environments are predicted to be targets of natural selection in variable environments. However, little is known about the genetic basis of the temporal persistence of environmental effects and their stability of expression across subsequent environments, or even the extent to which natural genotypes differ in these attributes of environmental effects. We factorially manipulated the thermal environment experienced in three successive generations, to quantify the temporal persistence and environmental stability of temperature effects in contrasting genotypes of Arabidopsis thaliana. We found that genotypes differed in the manner in which environmental effects dissipated across successive generations, the manner in which responses to ancestral environments were stably expressed in present environments, the manner in which ancestral environments altered responses to present environments, and in the manner in which ancestral environments altered fitness in present conditions. Genetic variation exists in nature for these trait-specific environmental responses, suggesting that the temporal persistence and stability of environmental effects in variable environments have the potential to evolve in response to natural selection imposed by different environments and sequences of environments.

Joint analysis of days to flowering reveals independent temperate adaptations in maize

Domesticates are an excellent model for understanding biological consequences of rapid climate change. Maize (Zea mays ssp. mays) was domesticated from a tropical grass yet is widespread across temperate regions today. We investigate the biological basis of temperate adaptation in diverse structured nested association mapping (NAM) populations from China, Europe (Dent and Flint) and the United States as well as in the Ames inbred diversity panel, using days to flowering as a proxy. Using cross-population prediction, where high prediction accuracy derives from overall genomic relatedness, shared genetic architecture, and sufficient diversity in the training population, we identify patterns in predictive ability across the five populations. To identify the source of temperate adapted alleles in these populations, we predict top associated genome-wide association study (GWAS) identified loci in a Random Forest Classifier using independent temperate-tropical North American populations based on lines selected from Hapmap3 as predictors. We find that North American populations are well predicted (AUC equals 0.89 and 0.85 for Ames and USNAM, respectively), European populations somewhat well predicted (AUC equals 0.59 and 0.67 for the Dent and Flint panels, respectively) and that the Chinese population is not predicted well at all (AUC is 0.47), suggesting an independent adaptation process for early flowering in China. Multiple adaptations for the complex trait days to flowering in maize provide hope for similar natural systems under climate change.

Life-history trade-offs and the genetic basis of fitness in Arabidopsis thaliana

Resources allocated to survival cannot be used to increase fecundity, but the extent to which this trade-off constrains adaptation depends on overall resource status. Adaptation to local environmental conditions may therefore entail the evolution of traits that increase the amount of resources available to individuals (their resource status or 'condition'). We examined the relative contribution of trade-offs and increased condition to adaptive evolution in a recombinant inbred line population of Arabidopsis thaliana planted at the native sites of the parental ecotypes in Italy and Sweden in 2 years. We estimated genetic correlations among fitness components based on genotypic means and explored their causes with QTL mapping. The local ecotype produced more seeds per fruit than did the non-local ecotype, reflected in stronger adaptive differentiation than was previously shown based on survival and fruit number only. Genetic correlations between survival and overall fecundity, and between number of fruits and number of seeds per fruit, were positive, and there was little evidence of a trade-off between seed size and number. Quantitative trait loci for these traits tended to map to the same regions of the genome and showed positive pleiotropic effects. The results indicate that adaptive differentiation between the two focal populations largely reflects the evolution of increased ability to acquire resources in the local environment, rather than shifts in the relative allocation to different life-history traits. Differentiation both in phenology and in tolerance to cold is likely to contribute to the advantage of the local genotype at the two sites.

Molecular variation in a functionally divergent homolog of FCA regulates flowering time in Arabidopsis thaliana

The identification and functional characterization of natural variants in plants are essential for understanding phenotypic adaptation. Here we identify a molecular variation in At2g47310 that contributes to the natural variation in flowering time in Arabidopsis thaliana accessions. This gene, which we term SISTER of FCA (SSF), functions in an antagonistic manner to its close homolog FCA. Genome-wide association analysis screens two major haplotypes of SSF associated with the natural variation in FLC expression, and a single polymorphism, SSF-N414D, is identified as a main contributor. The SSF414N protein variant interacts more strongly with CUL1, a component of the E3 ubiquitination complex, than the SSF414D form, mediating differences in SSF protein degradation and FLC expression. FCA and SSF appear to have arisen through gene duplication after dicot-monocot divergence, with the SSF-N414D polymorphism emerging relatively recently within A. thaliana. This work provides a good example for deciphering the functional importance of natural polymorphisms in different organisms.

Natural variation represents valuable source for gene discovery. Here, the authors show that a homolog of Flowering Control Locus A (FCA) functions in an antagonistic manner to FCA in regulating Arabidopsis flowering time through interacting with CUL1-E3 and modulating FLC expression.

A Nested Association Mapping Panel in Arabidopsis thaliana for Mapping and Characterizing Genetic Architecture

Linkage and association mapping populations are crucial public resources that facilitate the characterization of trait genetic architecture in natural and agricultural systems. We define a large nested association mapping panel (NAM) from 14 publicly available recombinant inbred line populations (RILs) of Arabidopsis thaliana, which share a common recurrent parent (Col-0). Using a genotype-by-sequencing approach (GBS), we identified single nucleotide polymorphisms (SNPs; range 563-1525 per population) and subsequently built updated linkage maps in each of the 14 RIL sets. Simulations in individual RIL populations indicate that our GBS markers have improved power to detect small effect QTL and enhanced resolution of QTL support intervals in comparison to original linkage maps. Using these robust linkage maps, we imputed a common set of publicly available parental SNPs into each RIL linkage map, generating overlapping markers across all populations. Though ultimately depending on allele frequencies at causal loci, simulations of the NAM panel suggest that surveying between 4 to 7 of the 14 RIL populations provides high resolution of the genetic architecture of complex traits, relative to a single mapping population.

Natural variation in autumn expression is the major adaptive determinant distinguishing Arabidopsis FLC haplotypes

In Arabidopsis thaliana, winter is registered during vernalization through the temperature-dependent repression and epigenetic silencing of floral repressor FLOWERING LOCUS C (FLC). Natural Arabidopsis accessions show considerable variation in vernalization. However, which aspect of the FLC repression mechanism is most important for adaptation to different environments is unclear. By analysing FLC dynamics in natural variants and mutants throughout winter in three field sites, we find that autumnal FLC expression, rather than epigenetic silencing, is the major variable conferred by the distinct Arabidopsis FLChaplotypes. This variation influences flowering responses of Arabidopsis accessions resulting in an interplay between promotion and delay of flowering in different climates to balance survival and, through a post-vernalization effect, reproductive output. These data reveal how expression variation through non-coding cis variation at FLC has enabled Arabidopsis accessions to adapt to different climatic conditions and year-on-year fluctuations.

Genotypic variation in the persistence of transgenerational responses to seasonal cues

Phenotypes respond to environments experienced directly by an individual, via phenotypic plasticity, or to the environment experienced by ancestors, via transgenerational environmental effects. The adaptive value of environmental effects depends not only on the strength and direction of the induced response but also on how long the response persists within and across generations, and how stably it is expressed across environments that are encountered subsequently. Little is known about the genetic basis of those distinct components, or even whether they exhibit genetic variation. We tested for genetic differences in the inducibility, temporal persistence, and environmental stability of transgenerational environmental effects in Arabidopsis thaliana. Genetic variation existed in the inducibility of transgenerational effects on traits expressed across the life cycle. Surprisingly, the persistence of transgenerational effects into the third generation was uncorrelated with their induction in the second generation. Although environmental effects for some traits in some genotypes weakened over successive generations, others were stronger or even in the opposite direction in more distant generations. Therefore, transgenerational effects in more distant generations are not merely caused by the retention or dissipation of those expressed in prior generations, but they may be genetically independent traits with the potential to evolve independently.

Genetic dynamics of earlier maturity group emergence in south-to-north extension of Northeast China soybeans

KEY MESSAGE

This population genetic study is characterized with direct comparisons of days to flowering QTL-allele matrices between newly evolved and originally old maturity groups of soybeans to explore its evolutionary dynamics using the RTM-GWAS procedure. The Northeast China (NEC) soybeans are the major germplasm source of modern soybean production in Americas (> 80% of the world total). NEC is a relatively new soybean area in China, expanded after its nomadic status in the seventeenth century. At nine sites of four ecoregions in NEC, 361 varieties were tested for their days to flowering (DTF), a geography-sensitive trait as an indicator for maturity groups (MGs). The DTF reduced obviously along with soybeans extended to higher latitudes, ranging in 41-83 days and MG 000-III. Using the RTM-GWAS (restricted two-stage multi-locus model genome-wide association study) procedure, 81 QTLs with 342 alleles were identified, accounting for 77.85% genetic contribution (R2 = 0.01-7.74%/locus), and other 20.75% (98.60-77.85%, h2 = 98.60%) genetic variation was due to a collective of unmapped QTLs. With soybeans northward, breeding effort made the original MG I-III evolved to MG 0-00-000. In direct comparisons of QTL-allele matrices among MGs, the genetic dynamics are identified with local exotic introduction/migration (58.48%) as the first and selection against/exclusion of positive alleles causing new recombination (40.64%) as the second, while only a few allele emergence/mutation happened (0.88%, limited in MG 0, not in MG 00-000). In new MG emergence, 24 QTLs with 19 candidate genes are the major sources. A genetic potential of further DTF shortening (13-21 days) is predicted for NEC population. The QTL detection in individual ecoregions showed various ecoregion-specific QTLs-alleles/genes after co-localization treatment (removing the random environment shifting ones).

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

Two reproductive traits show contrasting genetic architectures in Plantago lanceolata

In many species, temperature-sensitive phenotypic plasticity (i.e., an individual's phenotypic response to temperature) displays a positive correlation with latitude, a pattern presumed to reflect local adaptation. This geographical pattern raises two general questions: (a) Do a few large-effect genes contribute to latitudinal variation in a trait? (b) Is the thermal plasticity of different traits regulated pleiotropically? To address the questions, we crossed individuals of Plantago lanceolata derived from northern and southern European populations. Individuals naturally exhibited high and low thermal plasticity in floral reflectance and flowering time. We grew parents and offspring in controlled cool- and warm-temperature environments, mimicking what plants would encounter in nature. We obtained genetic markers via genotype-by-sequencing, produced the first recombination map for this ecologically important nonmodel species, and performed quantitative trait locus (QTL) mapping of thermal plasticity and single-environment values for both traits. We identified a large-effect QTL that largely explained the reflectance plasticity differences between northern and southern populations. We identified multiple smaller-effect QTLs affecting aspects of flowering time, one of which affected flowering time plasticity. The results indicate that the genetic architecture of thermal plasticity in flowering is more complex than for reflectance. One flowering time QTL showed strong cytonuclear interactions under cool temperatures. Reflectance and flowering plasticity QTLs did not colocalize, suggesting little pleiotropic genetic control and freedom for independent trait evolution. Such genetic information about the architecture of plasticity is environmentally important because it informs us about the potential for plasticity to offset negative effects of climate change.

The vernalisation regulator FLOWERING LOCUS C is differentially expressed in biennial and annual Brassica napus

Plants in temperate areas evolved vernalisation requirement to avoid pre-winter flowering. In Brassicaceae, a period of extended cold reduces the expression of the flowering inhibitor FLOWERING LOCUS C (FLC) and paves the way for the expression of downstream flowering regulators. As with all polyploid species of the Brassicaceae, the model allotetraploid Brassica napus (rapeseed, canola) is highly duplicated and carries 9 annotated copies of Bna.FLC. To investigate whether these multiple homeologs and paralogs have retained their original function in vernalisation or undergone subfunctionalisation, we compared the expression patterns of all 9 copies between vernalisation-dependent (biennial, winter type) and vernalisation-independent (annual, spring type) accessions, using RT-qPCR with copy-specific primers and RNAseq data from a diversity set. Our results show that only 3 copies – Bna.FLC.A03b, Bna.FLC.A10 and to some extent Bna.FLC.C02 – are differentially expressed between the two growth types, showing that expression of the other 6 copies does not correlate with growth type. One of those 6 copies, Bna.FLC.C03b, was not expressed at all, indicating a pseudogene, while three further copies, Bna.FLC.C03a and Bna.FLC.C09ab, did not respond to cold treatment. Sequence variation at the COOLAIR binding site of Bna.FLC.A10 was found to explain most of the variation in gene expression. However, we also found that Bna.FLC.A10 expression is not fully predictive of growth type.

A Polygenic Genetic Architecture of Flowering Time in the Worldwide Arabidopsis thaliana Population

Here, we report an empirical study of the polygenic basis underlying the evolution of complex traits. Flowering time variation measured at 10 and 16°C in the 1,001-genomes Arabidopsis thaliana collection of natural accessions were used as a model. The polygenic architecture of flowering time was defined as the 48 loci that were significantly associated with flowering time-at 10 and/or 16°C and/or their difference-in this population. Contributions from alleles at flowering time associated loci to global and local adaptation were explored by evaluating their distribution across genetically and geographically defined subpopulations across the native range of the species. The dynamics in the genetic architecture of flowering time in response to temperature was evaluated by estimating how the effects of these loci on flowering changed with growth temperature. Overall, the genetic basis of flowering time was stable-about 2/3 of the flowering time loci had similar effects at 10°C and 16°C-but many loci were involved in gene by temperature interactions. Globally present alleles, mostly of moderate effect, contributed to the differences in flowering times between the subpopulations via subtle changes in allele frequencies. More extreme local adaptations were, on several occasions, due to regional alleles with relatively large effects, and their linkage disequilibrium-patterns suggest coevolution of functionally connected alleles within local populations. Overall, these findings provide a significant contribution to our understanding about the possible modes of global and local evolution of a complex adaptive trait in A. thaliana.

The complex genetic architecture of shoot growth natural variation in Arabidopsis thaliana

One of the main outcomes of quantitative genetics approaches to natural variation is to reveal the genetic architecture underlying the phenotypic space. Complex genetic architectures are described as including numerous loci (or alleles) with small-effect and/or low-frequency in the populations, interactions with the genetic background, environment or age. Linkage or association mapping strategies will be more or less sensitive to this complexity, so that we still have an unclear picture of its extent. By combining high-throughput phenotyping under two environmental conditions with classical QTL mapping approaches in multiple Arabidopsis thaliana segregating populations as well as advanced near isogenic lines construction and survey, we have attempted to improve our understanding of quantitative phenotypic variation. Integrative traits such as those related to vegetative growth used in this work (highlighting either cumulative growth, growth rate or morphology) all showed complex and dynamic genetic architecture with respect to the segregating population and condition. The more resolutive our mapping approach, the more complexity we uncover, with several instances of QTLs visible in near isogenic lines but not detected with the initial QTL mapping, indicating that our phenotyping accuracy was less limiting than the mapping resolution with respect to the underlying genetic architecture. In an ultimate approach to resolve this complexity, we intensified our phenotyping effort to target specifically a 3Mb-region known to segregate for a major quantitative trait gene, using a series of selected lines recombined every 100kb. We discovered that at least 3 other independent QTLs had remained hidden in this region, some with trait- or condition-specific effects, or opposite allelic effects. If we were to extrapolate the figures obtained on this specific region in this particular cross to the genome- and species-scale, we would predict hundreds of causative loci of detectable phenotypic effect controlling these growth-related phenotypes.

Assortative mating by flowering time and its effect on correlated traits in variable environments

Reproductive timing is a key life-history trait that impacts the pool of available mates, the environment experienced during flowering, and the expression of other traits through genetic covariation. Selection on phenology, and its consequences on other life-history traits, has considerable implications in the context of ongoing climate change and shifting growing seasons. To test this, we grew field-collected seed from the wildflower Mimulus guttatus in a greenhouse to assess the standing genetic variation for flowering time and covariation with other traits. We then created full-sib families through phenological assortative mating and grew offspring in three photoperiod treatments representing seasonal variation in daylength. We find substantial quantitative genetic variation for the onset of flowering time, which covaried with vegetative traits. The assortatively-mated offspring varied in their critical photoperiod by over two hours, so that families differed in their probability of flowering across treatments Allocation to flowering and vegetative growth changed across the daylength treatments, with consistent direction and magnitude of covariation among flowering time and other traits. Our results suggest that future studies of flowering time evolution should consider the joint evolution of correlated traits and shifting seasonal selection to understand how environmental variation influences life histories.

Integrating Selection Mapping With Genetic Mapping and Functional Genomics

Genomic scans for signatures of selection allow us to, in principle, detect variants and genes that underlie recent adaptations. By combining selection mapping with genetic mapping of traits known to be relevant to adaptation, we can simultaneously investigate whether genes and variants show signals of recent selection and whether they impact traits that have likely been selected. There are three ways to integrate selection mapping with genetic mapping or functional genomics: (1) To use genetic mapping data from other populations as a form of genome annotation. (2) To perform experimental evolution or artificial selection to be able to study selected variants when they segregate, either by performing genetic mapping before selection or by crossing the selected individuals to some reference population. (3) To perform a comparative study of related populations facing different selection regimes. This short review discusses these different ways of integrating selection mapping with genetic mapping and functional genomics, with examples of how each has been done.

A Genomic Variation Map Provides Insights into the Genetic Basis of Spring Chinese Cabbage (Brassica rapa ssp. pekinensis) Selection

Chinese cabbage is the most consumed leafy crop in East Asian countries. However, premature bolting induced by continuous low temperatures severely decreases the yield and quality of the Chinese cabbage, and therefore restricts its planting season and geographic distribution. In the past 40 years, spring Chinese cabbage with strong winterness has been selected to meet the market demand. Here, we report a genome variation map of Chinese cabbage generated from the resequencing data of 194 geographically diverse accessions of three ecotypes. In-depth analyses of the selection sweeps and genome-wide patterns revealed that spring Chinese cabbage was selected from a specific population of autumn Chinese cabbage around the area of Shandong peninsula in northern China. We identified 23 genomic loci that underwent intensive selection, and further demonstrated by gene expression and haplotype analyses that the incorporation of elite alleles of VERNALISATION INSENTIVE 3.1 (BrVIN3.1) and FLOWER LOCUS C 1 (BrFLC1) is a determinant genetic source of variation during selection. Moreover, we showed that the quantitative response of BrVIN3.1 to cold due to the sequence variations in the cis elements of the BrVIN3.1 promoter significantly contributes to bolting-time variation in Chinese cabbage. Collectively, our study provides valuable insights into the genetic basis of spring Chinese cabbage selection and will facilitate the breeding of bolting-resistant varieties by molecular-marker-assisted selection, transgenic or gene editing approaches.

Transcriptomic analysis identifies differentially expressed genes (DEGs) associated with bolting and flowering in Saposhnikovia divaricata

Saposhnikovia divaricata is a valuable Chinese medicinal herb; the transformation from vegetative growth to reproductive growth may lead to the decrease of its pharmacological activities. Therefore, the study of bolting and flowering for Saposhnikovia divaricata is warranted. The present study aimed to reveal differentially expressed genes (DEGs) and regularity of expression during the bolting and flowering process, and the results of this study might provide a theoretical foundation for the suppression of early bolting for future research and practical application. Three sample groups, early flowering, flower bud differentiation, and late flowering (groups A, B, and C, respectively) were selected. Transcriptomic analysis identified 67, 010 annotated unigenes, among which 50, 165 were differentially expressed including 16, 108 in A vs B, and 17, 459 in B vs C, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional classification analysis were performed on these differentially expressed genes, and five important pathways were significantly impacted (P ≤ 0.01): plant circadian rhythm, other glycan degradation, oxidative phosphorylation, plant hormone signal transduction, and starch and sucrose metabolism. Plant hormone signal transduction might play an important role in the bolting and flowering process. The differentially expressed indole-3-acetic acid (IAA) gene showed significant down-regulation during bolting and flowering, while the transport inhibitor response 1 (TIR1) gene showed no significant change during the bolting process. The expression of flowering related genes FLC, LYF, and AP1 also showed a greater difference at different development stages. In conclusion, we speculate that the decrease in auxin concentration is not caused by the degrading effect of TIR1 but by an alternative mechanism.

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

To defend or to grow: lessons from Arabidopsis C24

The emergence of Arabidopsis as a model species and the availability of genetic and genomic resources have resulted in the identification and detailed characterization of abiotic stress signalling pathways. However, this has led only to limited success in engineering abiotic stress tolerance in crops. This is because there needs to be a deeper understanding of how to combine resistances to a range of stresses with growth and productivity. The natural variation and genomic resources of Arabidopsis thaliana (Arabidopsis) are a great asset to understand the mechanisms of multiple stress tolerances. One natural variant in Arabidopsis is the accession C24, and here we provide an overview of the increasing research interest in this accession. C24 is highlighted as a source of tolerance for multiple abiotic and biotic stresses, and a key accession to understand the basis of basal immunity to infection, high water use efficiency, and water productivity. Multiple biochemical, physiological, and phenological mechanisms have been attributed to these traits in C24, and none of them constrains productivity. Based on the uniqueness of C24, we postulate that the use of variation derived from natural selection in undomesticated species provides opportunities to better understand how complex environmental stress tolerances and resource use efficiency are co-ordinated.

Genetic basis of photosynthetic responses to cold in two locally adapted populations of Arabidopsis thaliana

Local adaptation is common, but the traits and genes involved are often unknown. Physiological responses to cold probably contribute to local adaptation in wide-ranging species, but the genetic basis underlying natural variation in these traits has rarely been studied. Using a recombinant inbred (495 lines) mapping population from locally adapted populations of Arabidopsis thaliana from Sweden and Italy, we grew plants at low temperature and mapped quantitative trait loci (QTLs) for traits related to photosynthesis: maximal quantum efficiency (Fv/Fm), rapidly reversible photoprotection (NPQfast), and photoinhibition of PSII (NPQslow) using high-throughput, whole-plant measures of chlorophyll fluorescence. In response to cold, the Swedish line had greater values for all traits, and for every trait, large effect QTLs contributed to parental differences. We found one major QTL affecting all traits, as well as unique major QTLs for each trait. Six trait QTLs overlapped with previously published locally adaptive QTLs based on fitness measured in the native environments over 3 years. Our results demonstrate that photosynthetic responses to cold can vary dramatically within a species, and may predominantly be caused by a few QTLs of large effect. Some photosynthesis traits and QTLs probably contribute to local adaptation in this system.

QTL Mapping of Seed Glucosinolate Content Responsible for Environment in Brassica napus

Glucosinolates (GSLs) are a major class of secondary metabolites. The content of seed GSL is largely regulated by environments in rapeseed (Brassica napus). However, the genetic control of seed GSL content responsible for environment in B. napus has been poorly understood. In the current study, a doubled haploid (DH) population from a cross between winter and semi-winter lines of rapeseed was grown in two distinct eco-environments, Germany and China, to evaluate the eco-environment effect and dissect the quantitative trait loci (QTL) responsible for environment for seed GSL in rapeseed. The deviation value of GSL content between eco-environments (GSLE) was calculated for each line in the DH population and the QTLs for GSLE were detected. GSLE ranged from -46.90 to 36.13 μmol g^-1 meal in the DH population, suggesting the prominent eco-environmental effects for seed GSL in rapeseed. Four QTLs for GSLE were identified on chromosomes A04, A06, and A09 explaining 4.70∼9.93% of the phenotypic variation. Comparison of QTLs of seed GSL content between different eco-environments found three QTLs for GSL on A02 from 37.6 to 45.4 cM, A04 from 0 to 17.2 cM, and A09 from 67.0 to 98.6 cM exhibited significant difference of QTL effect between the German and Chinese eco-environments (P < 0.01), indicating the environment sensibility of these loci on seed GSL content. Moreover, flowering time (FT), an important environment adaptation trait in plant, was also investigated in this study. Comparative QTL analysis among GSLE, GSL, and FT revealed that three regions on chromosomes A02, A04, and A09 not only exhibited significant differences in QTL effect between Germany and China, but also co-located with the QTL intervals of GSLE and FT. Our results revealed that most of the GSL loci can influence GSL accumulation under different eco-environments, whereas the three QTL intervals on A02, A04, and A09 might be sensitive to the eco-environments for seed GSL content.

Natural variation of life-history traits, water use, and drought responses in Arabidopsis

The ability of plants to acquire and use water is critical in determining life-history traits such as growth, flowering, and allocation of biomass into reproduction. In this context, a combination of functionally linked traits is essential for plants to respond to environmental changes in a coordinated fashion to maximize resource use efficiency. We analyzed different water-use traits in Arabidopsis ecotypes to identify functionally linked traits that determine water use and plant growth performance. Water-use traits measured were (i) leaf-level water-use efficiency (WUE i ) to evaluate the amount of CO 2 fixed relative to water loss per leaf area and (ii) short-term plant water use at the vegetative stage (VWU) as a measure of whole-plant transpiration. Previously observed phenotypic variance in VWU, WUE i and life-history parameters, highlighted C24 as a valuable ecotype that combined drought tolerance, preferential reproductive biomass allocation, high WUE i , and reduced water use. We therefore screened 35 Arabidopsis ecotypes for these parameters, in order to assess whether the phenotypic combinations observed in C24 existed more widely within Arabidopsis ecotypes. All parameters were measured on a short dehydration cycle. A segmented regression analysis was carried out to evaluate the plasticity of the drought response and identified the breakpoint as a reliable measure of drought sensitivity. VWU was largely dependent on rosette area, but importantly the drought sensitivity and plasticity measures were independent of the transpiring leaf surface. A breakpoint at high rSWC indicated a more drought-sensitive plant that closed stomata early during the dehydration cycle and consequently showed stronger plasticity in leaf-level WUE i parameters. None of the sensitivity, plasticity, or water-use measurements were able to predict the overall growth performance; however, there was a general trade-off between vegetative and reproductive biomass. PCA and hierarchical clustering revealed that C24 was unique among the 35 ecotypes in uniting all the beneficial water use and stress tolerance traits, while also maintaining above average plant growth. We propose that a short dehydration cycle, measuring drought sensitivity and VWU is a fast and reliable screen for plant water use and drought response strategies.

Fine-tuning timing: natural variation informs the mechanistic basis of the switch to flowering in Arabidopsis thaliana

The evolution of diverse life history strategies has allowed Arabidopsis thaliana to adapt to worldwide locations, spanning a range of latitudinal and environmental conditions. Arabidopsis thaliana accessions are either vernalization-requiring winter annuals or rapid cyclers, with extensive natural variation in vernalization requirement and response. Genetic and molecular analysis of this variation has enhanced our understanding of the mechanisms involved in life history determination, with translation to both natural and crop systems in the Brassicaceae and beyond.

Acclimation of Swedish and Italian ecotypes of Arabidopsis thaliana to light intensity

This study addressed whether ecotypes of Arabidopsis thaliana from Sweden and Italy exhibited differences in foliar acclimation to high versus low growth light intensity, and compared CO₂ uptake under growth conditions with light- and CO₂-saturated intrinsic photosynthetic capacity and leaf morphological and vascular features. Differential responses between ecotypes occurred mainly at the scale of leaf architecture, with thicker leaves with higher intrinsic photosynthetic capacities and chlorophyll contents per leaf area, but no difference in photosynthetic capacity on a chlorophyll basis, in high light-grown leaves of the Swedish versus the Italian ecotype. Greater intrinsic photosynthetic capacity per leaf area in the Swedish ecotype was accompanied by a greater capacity of vascular infrastructure for sugar and water transport, but this was not associated with greater CO₂ uptake rates under growth conditions. The Swedish ecotype with its thick leaves is thus constructed for high intrinsic photosynthetic and vascular flux capacity even under growth chamber conditions that may not permit full utilization of this potential. Conversely, the Swedish ecotype was less tolerant of low growth light intensity than the Italian ecotype, with smaller rosette areas and lesser aboveground biomass accumulation in low light-grown plants. Foliar vein density and stomatal density were both enhanced by high growth light intensity with no significant difference between ecotypes, and the ratio of water to sugar conduits was also similar between the two ecotypes during light acclimation. These findings add to the understanding of the foliar vasculature's role in plant photosynthetic acclimation and adaptation.

Between semelparity and iteroparity: Empirical evidence for a continuum of modes of parity

The number of times an organism reproduces (i.e., its mode of parity) is a fundamental life-history character, and evolutionary and ecological models that compare the relative fitnesses of different modes of parity are common in life-history theory and theoretical biology. Despite the success of mathematical models designed to compare intrinsic rates of increase (i.e., density-independent growth rates) between annual-semelparous and perennial-iteroparous reproductive schedules, there is widespread evidence that variation in reproductive allocation among semelparous and iteroparous organisms alike is continuous. This study reviews the ecological and molecular evidence for the continuity and plasticity of modes of parity-that is, the idea that annual-semelparous and perennial-iteroparous life histories are better understood as endpoints along a continuum of possible strategies. I conclude that parity should be understood as a continuum of different modes of parity, which differ by the degree to which they disperse or concentrate reproductive effort in time. I further argue that there are three main implications of this conclusion: (1) that seasonality should not be conflated with parity; (2) that mathematical models purporting to explain the general evolution of semelparous life histories from iteroparous ones (or vice versa) should not assume that organisms can only display either an annual-semelparous life history or a perennial-iteroparous one; and (3) that evolutionary ecologists should base explanations of how different life-history strategies evolve on the physiological or molecular basis of traits underlying different modes of parity.

Interacting effects of genetic variation for seed dormancy and flowering time on phenology, life history, and fitness of experimental Arabidopsis thaliana populations over multiple generations in the field

Major alleles for seed dormancy and flowering time are well studied, and can interact to influence seasonal timing and fitness within generations. However, little is known about how this interaction controls phenology, life history, and population fitness across multiple generations in natural seasonal environments. To examine how seed dormancy and flowering time shape annual plant life cycles over multiple generations, we established naturally dispersing populations of recombinant inbred lines of Arabidopsis thaliana segregating early and late alleles for seed dormancy and flowering time in a field experiment. We recorded seasonal phenology and fitness of each genotype over 2 yr and several generations. Strong seed dormancy suppressed mid-summer germination in both early- and late-flowering genetic backgrounds. Strong dormancy and late-flowering genotypes were both necessary to confer a winter annual life history; other genotypes were rapid-cycling. Strong dormancy increased within-season fecundity in an early-flowering background, but decreased it in a late-flowering background. However, there were no detectable differences among genotypes in population growth rates. Seasonal phenology, life history, and cohort fitness over multiple generations depend strongly upon interacting genetic variation for dormancy and flowering. However, similar population growth rates across generations suggest that different life cycle genotypes can coexist in natural populations.

Circadian rhythms vary over the growing season and correlate with fitness components

Circadian clocks have evolved independently in all three domains of life, suggesting that internal mechanisms of time-keeping are adaptive in contemporary populations. However, the performance consequences of either discrete or quantitative clock variation have rarely been tested in field settings. Clock sensitivity of diverse segregating lines to the environment remains uncharacterized as do the statistical genetic parameters that determine evolutionary potential. In field studies with Arabidopsis thaliana, we found that major perturbations to circadian cycle length (referred to as clock period) via mutation reduce both survival and fecundity. Subtler adjustments via genomic introgression of naturally occurring alleles indicated that clock periods slightly >24 hr were adaptive, consistent with prior models describing how well the timing of biological processes is adjusted within a diurnal cycle (referred to as phase). In segregating recombinant inbred lines (RILs), circadian phase varied up to 2 hr across months of the growing season, and both period and phase expressed significant genetic variances. Performance metrics including developmental rate, size and fruit set were described by principal components (PC) analyses and circadian parameters correlated with the first PC, such that period lengths slightly >24 hr were associated with improved performance in multiple RIL sets. These experiments translate functional analyses of clock behaviour performed in controlled settings to natural ones, demonstrating that quantitative variation in circadian phase is highly responsive to seasonally variable abiotic factors. The results expand upon prior studies in controlled settings, showing that discrete and quantitative variation in clock phenotypes correlates with performance in nature.

Transcriptome predictors of coral survival and growth in a highly variable environment

Concern over rapid environmental shifts associated with climate change has led to a search for molecular markers of environmental tolerance. Climate-associated gene expression profiles exist for a number of systems, but have rarely been tied to fitness outcomes, especially in nonmodel organisms. We reciprocally transplanted corals between two backreef locations with more and less variable temperature regimes to disentangle effects of recent and native environment on survival and growth. Coral growth over 12 months was largely determined by local environment. Survival, however, was impacted by native environment; corals from the more variable environment had 22% higher survivorship. By contrast, corals native to the less variable environment had more variable survival. This might represent a "selective sieve" where poor survivors are filtered from the more stressful environment. We also find a potential fitness trade-off-corals with high survival under stressful conditions grew less in the more benign environment. Transcriptome samples taken a year before transplantation were used to examine gene expression patterns that predicted transplant survival and growth. Two separate clusters of coexpressed genes were predictive of survival in the two locations. Genes from these clusters are candidate biomarkers for predicting persistence of corals under future climate change scenarios.

Association mapping of loci controlling genetic and environmental interaction of soybean flowering time under various photo-thermal conditions

BACKGROUND

Soybean (Glycine max (L.) Merr.) is a short day plant. Its flowering and maturity time are controlled by genetic and environmental factors, as well the interaction between the two factors. Previous studies have shown that both genetic and environmental factors, mainly photoperiod and temperature, control flowering time of soybean. Additionally, these studies have reported gene × gene and gene × environment interactions on flowering time. However, the effects of quantitative trait loci (QTL) in response to photoperiod and temperature have not been well evaluated. The objectives of the current study were to identify the effects of loci associated with flowering time under different photo-thermal conditions and to understand the effects of interaction between loci and environment on soybean flowering.

METHODS

Different photoperiod and temperature combinations were obtained by adjusting sowing dates (spring sowing and summer sowing) or day-length (12 h, 16 h). Association mapping was performed on 91 soybean cultivars from different maturity groups (MG000-VIII) using 172 SSR markers and 5107 SNPs from the Illumina SoySNP6K iSelectBeadChip. The effects of the interaction between QTL and environments on flowering time were also analysed using the QTXNetwork.

RESULTS

Large-effect loci were detected on Gm 11, Gm 16 and Gm 20 as in previous reports. Most loci associated with flowering time are sensitive to photo-thermal conditions. Number of loci associated with flowering time was more under the long day (LD) than under the short day (SD) condition. The variation of flowering time among the soybean cultivars mostly resulted from the epistasis × environment and additive × environment interactions. Among the three candidate loci, i.e. Gm04_4497001 (near GmCOL3a), Gm16_30766209 (near GmFT2a and GmFT2b) and Gm19_47514601 (E3 or GmPhyA3), the Gm04_4497001 may be the key locus interacting with other loci for controlling soybean flowering time.

CONCLUSION

The effects of loci associated with the flowering time of soybean were dependent upon the photo-thermal conditions. This study facilitates the understanding of the genetic mechanism of soybean flowering and molecular breeding for the improvement of soybean adaptability to specific and/or broad regions.

Post-polyploidisation morphotype diversification associates with gene copy number variation

Genetic models for polyploid crop adaptation provide important information relevant for future breeding prospects. A well-suited model is Brassica napus, a recent allopolyploid closely related to Arabidopsis thaliana. Flowering time is a major adaptation trait determining life cycle synchronization with the environment. Here we unravel natural genetic variation in B. napus flowering time regulators and investigate associations with evolutionary diversification into different life cycle morphotypes. Deep sequencing of 35 flowering regulators was performed in 280 diverse B. napus genotypes. High sequencing depth enabled high-quality calling of single-nucleotide polymorphisms (SNPs), insertion-deletions (InDels) and copy number variants (CNVs). By combining these data with genotyping data from the Brassica 60 K Illumina® Infinium SNP array, we performed a genome-wide marker distribution analysis across the 4 ecogeographical morphotypes. Twelve haplotypes, including Bna.FLC.A10, Bna.VIN3.A02 and the Bna.FT promoter on C02_random, were diagnostic for the diversification of winter and spring types. The subspecies split between oilseed/kale (B. napus ssp. napus) and swedes/rutabagas (B. napus ssp. napobrassica) was defined by 13 haplotypes, including genomic rearrangements encompassing copies of Bna.FLC, Bna.PHYA and Bna.GA3ox1. De novo variation in copies of important flowering-time genes in B. napus arose during allopolyploidisation, enabling sub-functionalisation that allowed different morphotypes to appropriately fine-tune their lifecycle.