Genetic architecture of flowering time differentiation between locally adapted populations of Arabidopsis thaliana

An overview of floral regulatory genes in annual and perennial plants

The floral initiation in angiosperms is a complex process influenced by endogenous and exogenous signals. With this approach, we aim to provide a comprehensive review to integrate this complex floral regulatory process and summarize the regulatory genes and their functions in annuals and perennials. Seven primary paths leading to flowering have been discovered in Arabidopsis under several growth condition that include; photoperiod, ambient temperature, vernalization, gibberellins, autonomous, aging and carbohydrates. These pathways involve a series of interlinked signaling pathways that respond to both internal and external signals, such as light, temperature, hormones, and developmental cues, to coordinate the expression of genes that are involved in flower development. Among them, the photoperiodic pathway was the most important and conserved as some of the fundamental loci and mechanisms are shared even by closely related plant species. The activation of floral regulatory genes such as FLC, FT, LFY, and SOC1 that determine floral meristem identity and the transition to the flowering stage result from the merging of these pathways. Recent studies confirmed that alternative splicing, antisense RNA and epigenetic modification play crucial roles by regulating the expression of genes related to blooming. In this review, we documented recent progress in the floral transition time in annuals and perennials, with emphasis on the specific regulatory mechanisms along with the application of various molecular approaches including overexpression studies, RNA interference and Virus-induced flowering. Furthermore, the similarities and differences between annual and perennial flowering will aid significant contributions to the field by elucidating the mechanisms of perennial plant development and floral initiation regulation.

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.

A Genome-Wide Association Study Coupled With a Transcriptomic Analysis Reveals the Genetic Loci and Candidate Genes Governing the Flowering Time in Alfalfa (Medicago sativa L.)

The transition to flowering at the right time is very important for adapting to local conditions and maximizing alfalfa yield. However, the understanding of the genetic basis of the alfalfa flowering time remains limited. There are few reliable genes or markers for selection, which hinders progress in genetic research and molecular breeding of this trait in alfalfa. We sequenced 220 alfalfa cultivars and conducted a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs). The phenotypic analysis showed that the breeding status and geographical origin strongly influenced the alfalfa flowering time. Our GWAS revealed 63 loci significantly related to the flowering time. Ninety-five candidate genes were detected at these SNP loci within 40 kb (20 kb up- and downstream). Thirty-six percent of the candidate genes are involved in development and pollen tube growth, indicating that these genes are key genetic mechanisms of alfalfa growth and development. The transcriptomic analysis showed that 1,924, 2,405, and 3,779 differentially expressed genes (DEGs) were upregulated across the three growth stages, while 1,651, 2,613, and 4,730 DEGs were downregulated across the stages. Combining the results of our GWAS and transcriptome analysis, in total, 38 candidate genes (7 differentially expressed during the bud stage, 13 differentially expressed during the initial flowering stage, and 18 differentially expressed during the full flowering stage) were identified. Two SNPs located in the upstream region of the Msa0888690 gene (which is involved in isop renoids) were significantly related to flowering. The two significant SNPs within the upstream region of Msa0888690 existed as four different haplotypes in this panel. The genes identified in this study represent a series of candidate targets for further research investigating the alfalfa flowering time and could be used for alfalfa molecular breeding.

Distinct Cold Acclimation of Productivity Traits in Arabidopsis thaliana Ecotypes

Improvement of crop climate resilience will require an understanding of whole-plant adaptation to specific local environments. This review places features of plant form and function related to photosynthetic productivity, as well as associated gene-expression patterns, into the context of the adaptation of Arabidopsis thaliana ecotypes to local environments with different climates in Sweden and Italy. The growth of plants under common cool conditions resulted in a proportionally greater emphasis on the maintenance of photosynthetic activity in the Swedish ecotype. This is compared to a greater emphasis on downregulation of light-harvesting antenna size and upregulation of a host of antioxidant enzymes in the Italian ecotype under these conditions. This differential response is discussed in the context of the climatic patterns of the ecotypes’ native habitats with substantial opportunity for photosynthetic productivity under mild temperatures in Italy but not in Sweden. The Swedish ecotype’s response is likened to pushing forward at full speed with productivity under low temperature versus the Italian ecotype’s response of staying safe from harm (maintaining redox homeostasis) while letting productivity decline when temperatures are transiently cold. It is concluded that either strategy can offer directions for the development of climate-resilient crops for specific locations of cultivation.

Detection of the Progression of Anthesis in Field-Grown Maize Tassels: A Case Study

The tassel of the maize plant is responsible for the production and dispersal of pollen for subsequent capture by the silk (stigma) and fertilization of the ovules. Both the amount and timing of pollen shed are physiological traits that impact the production of a hybrid seed. This study describes an automated end-to-end pipeline that combines deep learning and image processing approaches to extract tassel flowering patterns from time-lapse camera images of plants grown under field conditions. Inbred lines from the SAM and NAM diversity panels were grown at the Curtiss farm at Iowa State University, Ames, IA, during the summer of 2016. Using a set of around 500 pole-mounted cameras installed in the field, images of plants were captured every 10 minutes of daylight hours over a three-week period. Extracting data from imaging performed under field conditions is challenging due to variabilities in weather, illumination, and the morphological diversity of tassels. To address these issues, deep learning algorithms were used for tassel detection, classification, and segmentation. Image processing approaches were then used to crop the main spike of the tassel to track reproductive development. The results demonstrated that deep learning with well-labeled data is a powerful tool for detecting, classifying, and segmenting tassels. Our sequential workflow exhibited the following metrics: mAP for tassel detection was 0.91, F1 score obtained for tassel classification was 0.93, and accuracy of semantic segmentation in creating a binary image from the RGB tassel images was 0.95. This workflow was used to determine spatiotemporal variations in the thickness of the main spike-which serves as a proxy for anthesis progression.

QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus)

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

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.

Synchronization of senescence and desynchronization of flowering in Arabidopsis thaliana

In a recent publication, we proposed that adjusting lifespan in order to synchronize senescence is important for timing of reproduction, and we quantified the synchrony of reproductive timing relative to germination timing. Here, in a second sequential seeding experiment (SSE), the germination timing of Arabidopsis thaliana accessions was manipulated and plants were then grown under two different temperature regimes. Life stage traits of plants in each temperature regime were analysed and it was evaluated whether the cohorts were grouped according to age and/or environmental conditions. While flowering-related traits showed desynchrony among cohorts, striking synchrony in the timing of senescence among cohorts for each group was found. A quantitative trait locus (QTL) analysis using a genotyped population of 'Cvi/Ler' recombinant inbred lines (RILs) was then conducted. Novel and known loci were assigned to flowering and senescence timing. However, senescence synchrony resulted in low variation in senescence time and weak QTL detection for flowering termination. Overlapping flowering and senescence genes with loci affecting either of those traits were found and suggest a potential interdependency of reproductive traits.

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.

QTL mapping of flowering time and biomass yield in tetraploid alfalfa (Medicago sativa L.)

BACKGROUND

The genetic and genomic basis of flowering time and biomass yield in alfalfa (Medicago sativa L.) remains poorly understood mainly due to the autopolyploid nature of the species and the lack of adequate genomic resources. We constructed linkage maps using genotyping-by-sequencing (GBS) based single dose allele (SDA) SNP and mapped alfalfa timing of flowering (TOF), spring yield (SY), and cumulative summer biomass (CSB) in a pseudo-testcross F1 population derived from a fall dormant (3010) and a non-dormant (CW 1010) cultivars. We analyzed the quantitative trait loci (QTL) to identify conserved genomic regions and detected molecular markers and potential candidate genes associated with the traits to improve alfalfa and provide genomic resources for the future studies.

RESULTS

This study showed that both fall dormant and non-dormant alfalfa cultivars harbored QTL for early and late flowering, suggesting that flowering time in alfalfa is not an indicator of its fall dormancy (FD) levels. A weak phenotypic correlation between the flowering time and fall dormancy (FD) in F1 and checks also corroborated that alfalfa FD and TOF are not the predictors of one another. The relationship between flowering time and alfalfa biomass yield was not strong, but the non-dormant had relatively more SY than dormant. Therefore, selecting superior alfalfa cultivars that are non-dormant, winter-hardy, and early flowering would allow for an early spring harvest with enhanced biomass. In this study, we found 25 QTL for TOF, 17 for SY and six QTL for CSB. Three TOF related QTL were stable and four TOF QTL were detected in the corresponding genomic locations of the flowering QTL of M. truncatula, an indication of possible evolutionarily conserved regions. The potential candidate genes for the SNP sequences of QTL regions were identified for all three traits and these genes would be potential targets for further molecular studies.

CONCLUSIONS

This research showed that variation in alfalfa flowering time after spring green up has no association with dormancy levels. Here we reported QTL, markers, and potential candidate genes associated with spring flowering time and biomass yield of alfalfa, which constitute valuable genomic resources for improving these traits via marker-assisted selection (MAS).

Molecular markers based on sequence variation in BoFLC1.C9 for characterizing early- and late-flowering cabbage genotypes

BACKGROUND

Cabbage (Brassica oleracea var. capitata) is popular worldwide for consumption as a leafy vegetable. Premature flowering is triggered by low temperature, and deteriorates quality of cabbage as vegetable. In general, growers prefer late-flowering varieties to assure good quality compact head. Here, we report BoFLC1.C9 as a gene with clear sequence variation between cabbage lines with different flowering times, and proposed as molecular marker to characterize early- and late-flowering cabbage lines.

RESULTS

We identified sequence variation of 67 bp insertions in intron 2, which were contributed in flowering time variation between two inbred lines through rapid down-regulation of the BoFLC1.C9 gene in early-flowering line compared to late-flowering one upon vernalization. One set of primer 'F7R7' proposed as marker, of which was explained with 83 and 80% of flowering time variation in 141 F₂ individuals and 20 commercial lines, respectively.

CONCLUSIONS

This F7R7 marker could be used as genetic tools to characterize flowering time variation and to select as well to develop early- and late-flowering cabbage cultivars.

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.

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.

Array of MADS-Box Genes: Facilitator for Rapid Adaptation?

In a world of global warming, the question emerges whether all plants have suitable mechanisms to keep pace with the rapidly changing environment. Most previous studies have focused on either the ability of plants to rapidly acclimatize via physiological and developmental plasticity, or long-term adaptation over thousands of years. However, we wonder whether plants can also adapt to changes in the environment within only a few generations. We hypothesize that rapidly evolving clusters of tandemly duplicated developmental control genes represent a source for fast adaptation. Specifically, we propose that a tandem cluster of FLC-like MADS-box genes involved in the transition to flowering in Arabidopsis functions as a facilitator for rapid adaptation to changes in ambient temperature.

Natural variation in genes potentially involved in plant architecture and adaptation in switchgrass (Panicum virgatum L.)

Background

Advances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.

Results

A total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima’s test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.

Conclusions

Several non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1193-2) contains supplementary material, which is available to authorized users.

Optimization of Photosynthetic Productivity in Contrasting Environments by Regulons Controlling Plant Form and Function

We review the role of a family of transcription factors and their regulons in maintaining high photosynthetic performance across a range of challenging environments with a focus on extreme temperatures and water availability. Specifically, these transcription factors include CBFs (C-repeat binding factors) and DREBs (dehydration-responsive element-binding), with CBF/DREB1 primarily orchestrating cold adaptation and other DREBs serving in heat, drought, and salinity adaptation. The central role of these modulators in plant performance under challenging environments is based on (i) interweaving of these regulators with other key signaling networks (plant hormones and redox signals) as well as (ii) their function in integrating responses across the whole plant, from light-harvesting and sugar-production in the leaf to foliar sugar export and water import and on to the plant's sugar-consuming sinks (growth, storage, and reproduction). The example of Arabidopsisthaliana ecotypes from geographic origins with contrasting climates is used to describe the links between natural genetic variation in CBF transcription factors and the differential acclimation of plant anatomical and functional features needed to support superior photosynthetic performance in contrasting environments. Emphasis is placed on considering different temperature environments (hot versus cold) and light environments (limiting versus high light), on trade-offs between adaptations to contrasting environments, and on plant lines minimizing such trade-offs.

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.

The FLC Locus: A Platform for Discoveries in Epigenetics and Adaptation

Our understanding of the detailed molecular mechanisms underpinning adaptation is still poor. One example for which mechanistic understanding of regulation has converged with studies of life history variation is Arabidopsis thaliana FLOWERING LOCUS C (FLC). FLC determines the need for plants to overwinter and their ability to respond to prolonged cold in a process termed vernalization. This review highlights how molecular analysis of vernalization pathways has revealed important insight into antisense-mediated chromatin silencing mechanisms that regulate FLC. In turn, such insight has enabled molecular dissection of the diversity in vernalization across natural populations of A. thaliana. Changes in both cotranscriptional regulation and epigenetic silencing of FLC are caused by noncoding polymorphisms at FLC. The FLC locus is therefore providing important concepts for how noncoding transcription and chromatin regulation influence gene expression and how these mechanisms can vary to underpin adaptation in natural populations.

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.

Habitat Temperature and Precipitation of Arabidopsis thaliana Ecotypes Determine the Response of Foliar Vasculature, Photosynthesis, and Transpiration to Growth Temperature

Acclimatory adjustments of foliar vascular architecture, photosynthetic capacity, and transpiration rate in Arabidopsis thaliana ecotypes (Italian, Polish [Col-0], Swedish) were characterized in the context of habitat of origin. Temperatures of the habitat of origin decreased linearly with increasing habitat latitude, but habitat precipitation was greatest in Italy, lowest in Poland, and intermediate in Sweden. Plants of the three ecotypes raised under three different growth temperature regimes (low, moderate, and high) exhibited highest photosynthetic capacities, greatest leaf thickness, highest chlorophyll a/b ratio and levels of β-carotene, and greatest levels of wall ingrowths in phloem transfer cells, and, in the Col-0 and Swedish ecotypes, of phloem per minor vein in plants grown at the low temperature. In contrast, vein density and minor vein tracheary to sieve element ratio increased with increasing growth temperature - most strongly in Col-0 and least strongly in the Italian ecotype - and transpirational water loss correlated with vein density and number of tracheary elements per minor vein. Plotting of these vascular features as functions of climatic conditions in the habitat of origin suggested that temperatures during the evolutionary history of the ecotypes determined acclimatory responses of the foliar phloem and photosynthesis to temperature in this winter annual that upregulates photosynthesis in response to lower temperature, whereas the precipitation experienced during the evolutionary history of the ecotypes determined adjustment of foliar vein density, xylem, and transpiration to temperature. In particular, whereas photosynthetic capacity, leaf thickness, and foliar minor vein phloem features increased linearly with increasing latitude and decreasing temperature of the habitats of origin in response to experimental growth at low temperature, transpiration rate, foliar vein density, and minor vein tracheary element numbers and cross-sectional areas increased linearly with decreasing precipitation level in the habitats of origin in response to experimental growth at high temperature. This represents a situation where temperature acclimation of the apparent capacity for water flux through the xylem and transpiration rate in a winter annual responded differently from that of photosynthetic capacity, in contrast to previous reports of strong relationships between hydraulic conductance and photosynthesis in other studies.

Growth temperature impact on leaf form and function in Arabidopsis thaliana ecotypes from northern and southern Europe

The plasticity of leaf form and function in European lines of Arabidopsis thaliana was evaluated in ecotypes from Sweden and Italy grown under contrasting (cool versus hot) temperature regimes. Although both ecotypes exhibited acclimatory adjustments, the Swedish ecotype exhibited more pronounced responses to the two contrasting temperature regimes in several characterized features. These responses included thicker leaves with higher capacities for photosynthesis, likely facilitated by a greater number of phloem cells per minor vein for the active loading and export of sugars, when grown under cool temperature as opposed to leaves with a higher vein density and a greater number of tracheary elements per minor vein, likely facilitating higher rates of transpirational water loss (and thus evaporative cooling), when grown under hot temperature with high water availability. In addition, only the Swedish ecotype exhibited reduced rosette growth and greater levels of foliar tocopherols under the hot growth temperature. These responses, and the greater responsiveness of the Swedish ecotype compared with the Italian ecotype, are discussed in the context of redox signalling networks and transcription factors, and the greater range of environmental conditions experienced by the Swedish versus the Italian ecotype during the growing season in their native habitats.

Genome to Phenome Mapping in Apple Using Historical Data

Apple ( X Borkh.) is one of the world's most valuable fruit crops. Its large size and long juvenile phase make it a particularly promising candidate for marker-assisted selection (MAS). However, advances in MAS in apple have been limited by a lack of phenotype and genotype data from sufficiently large samples. To establish genotype-phenotype relationships and advance MAS in apple, we extracted over 24,000 phenotype scores from the USDA-Germplasm Resources Information Network (GRIN) database and linked them with over 8000 single nucleotide polymorphisms (SNPs) from 689 apple accessions from the USDA apple germplasm collection clonally preserved in Geneva, NY. We find significant genetic differentiation between Old World and New World cultivars and demonstrate that the genetic structure of the domesticated apple also reflects the time required for ripening. A genome-wide association study (GWAS) of 36 phenotypes confirms the association between fruit color and the MYB1 locus, and we also report a novel association between the transcription factor, NAC18.1, and harvest date and fruit firmness. We demonstrate that harvest time and fruit size can be predicted with relatively high accuracies ( > 0.46) using genomic prediction. Rapid decay of linkage disequilibrium (LD) in apples means millions of SNPs may be required for well-powered GWAS. However, rapid LD decay also promises to enable extremely high resolution mapping of causal variants, which holds great potential for advancing MAS.

Population structure of a vector-borne plant parasite

Parasites are among the most diverse groups of life on Earth, yet complex natural histories often preclude studies of their speciation processes. The biology of parasitic plants facilitates in situ collection of data on both genetic structure and the mechanisms responsible for that structure. Here, we studied the role of mating, dispersal and establishment in host race formation of a parasitic plant. We investigated the population genetics of a vector-borne desert mistletoe (Phoradendron californicum) across two legume host tree species (Senegalia greggii and Prosopis velutina) in the Sonoran desert using microsatellites. Consistent with host race formation, we found strong host-associated genetic structure in sympatry, little genetic variation due to geographic site and weak isolation by distance. We hypothesize that genetic differentiation results from differences in the timing of mistletoe flowering by host species, as we found initial flowering date of individual mistletoes correlated with genetic ancestry. Hybrids with intermediate ancestry were detected genetically. Individuals likely resulting from recent, successful establishment events following dispersal between the host species were detected at frequencies similar to hybrids between host races. Therefore, barriers to gene flow between the host races may have been stronger at mating than at dispersal. We also found higher inbreeding and within-host individual relatedness values for mistletoes on the more rare and isolated host species (S. greggii). Our study spanned spatial scales to address how interactions with both vectors and hosts influence parasitic plant structure with implications for parasite virulence evolution and speciation.

Differences in light-harvesting, acclimation to growth-light environment, and leaf structural development between Swedish and Italian ecotypes of Arabidopsis thaliana

Leaf morphological differences have an impact on light distribution within the leaf, photosynthesis, and photoprotection in Arabidopsis thaliana ecotypes from near the limits of this species' latitudinal distribution in Europe. Leaf morphology, photosynthesis, and photoprotection were characterized in two Arabidopsis ecotypes from near the limits of this species' latitudinal distribution in Europe (63°N and 42°N). The Swedish ecotype formed thicker leaves and upregulated photosynthesis more substantially than the Italian ecotype in high-light environments. Conversely, the smaller rosette formed, and lesser aboveground biomass accumulated, by the Swedish versus the Italian ecotype in low growth-light environments is consistent with a lesser shade tolerance of the Swedish ecotype. The response of the thinner leaves of the Italian ecotype to evenly spaced daily periods of higher light against a background of otherwise non-fluctuating low light was to perform the same rate of photosynthesis with less chlorophyll, rather than exhibiting greater rates of photosynthesis. In contrast, the thicker leaves of the Swedish ecotype showed elevated photosynthetic performance in response to daily supplemental higher light periods in a low-light growth environment. These findings suggest significant self-shading in the lower depths of leaves of the Swedish ecotype by the chloroplasts residing in the upper portions of the leaf, resulting in a requirement for higher incident light to trigger photosynthetic upregulation in the lower portions of its thicker leaves. Conversely, photoprotective responses in the Italian ecotype suggest that more excess light penetrated into the lower depths of this ecotype's leaves. It is speculated that light absorption and the degree of utilization of this absorbed light inform cellular signaling networks that orchestrate leaf structural development, which, in turn, affects light distribution and the level of absorbed versus photosynthetically utilized light in a leaf.

Genetic variation in the flowering and yield formation of timothy (Phleum pratense L.) accessions after different photoperiod and vernalization treatments

Timothy is a perennial forage grass grown commonly in Boreal regions. This study explored the effect of vernalization and photoperiod (PP) on flowering and growth characteristics and how this related to changes in expression of three flowering related genes in accessions from different geographic origin. Large variation was found in accessions in their vernalization and PP responses. In southern accessions vernalization response or requirement was not observed, the heading date remained unchanged, and plants flowered without vernalization. On the contrary, northern types had obligatory requirement for vernalization and long PP, but the tiller elongation did not require vernalization at 16-h PP. Longer vernalization or PP treatments reduced the genotypical differences in flowering. Moreover, the vernalization saturation progressed stepwise from main tiller to lateral tillers, and this process was more synchronized in southern accessions. The expression of PpVRN1 was associated with vernalization while PpVRN3 accumulated at long PP. A crucial role for PpVRN3 in the transition to flowering was supported as in southern accession the transcript accumulated in non-vernalized plants after transfer to 16-h PP, and the apices transformed to generative stage. Differences in vernalization requirements were associated with variation in expression levels of PpVRN1 and PpVRN3, with higher expression levels in southern type. Most divergent transcript accumulation of PpMADS10 was found under different vernalization conditions. These differences between accessions can be translated into agronomic traits, such as the tiller composition of canopy, which affects the forage yield. The southern types, with minimal vernalization response, have fast re-growth ability and rapidly decreasing nutritive value, whereas northern types grow slowly and have better quality. This information can be utilized in breeding for new cultivars for longer growing seasons at high latitudes.

Heterosis and outbreeding depression in crosses between natural populations of Arabidopsis thaliana

Understanding the causes and architecture of genetic differentiation between natural populations is of central importance in evolutionary biology. Crosses between natural populations can result in heterosis if recessive or nearly recessive deleterious mutations have become fixed within populations because of genetic drift. Divergence between populations can also result in outbreeding depression because of genetic incompatibilities. The net fitness consequences of between-population crosses will be a balance between heterosis and outbreeding depression. We estimated the magnitude of heterosis and outbreeding depression in the highly selfing model plant Arabidopsis thaliana, by crossing replicate line pairs from two sets of natural populations (C↔R, B↔S) separated by similar geographic distances (Italy↔Sweden). We examined the contribution of different modes of gene action to overall differences in estimates of lifetime fitness and fitness components using joint scaling tests with parental, reciprocal F1 and F2, and backcross lines. One of these population pairs (C↔R) was previously demonstrated to be locally adapted, but locally maladaptive quantitative trait loci were also found, suggesting a role for genetic drift in shaping adaptive variation. We found markedly different genetic architectures for fitness and fitness components in the two sets of populations. In one (C↔R), there were consistently positive effects of dominance, indicating the masking of recessive or nearly recessive deleterious mutations that had become fixed by genetic drift. The other set (B↔S) exhibited outbreeding depression because of negative dominance effects. Additional studies are needed to explore the molecular genetic basis of heterosis and outbreeding depression, and how their magnitudes vary across environments.

Natural variation in cross-talk between glucosinolates and onset of flowering in Arabidopsis

Naturally variable regulatory networks control different biological processes including reproduction and defense. This variation within regulatory networks enables plants to optimize defense and reproduction in different environments. In this study we investigate the ability of two enzyme-encoding genes in the glucosinolate pathway, AOP2 and AOP3, to affect glucosinolate accumulation and flowering time. We have introduced the two highly similar enzymes into two different AOP (null) accessions, Col-0 and Cph-0, and found that the genes differ in their ability to affect glucosinolate levels and flowering time across the accessions. This indicated that the different glucosinolates produced by AOP2 and AOP3 serve specific regulatory roles in controlling these phenotypes. While the changes in glucosinolate levels were similar in both accessions, the effect on flowering time was dependent on the genetic background pointing to natural variation in cross-talk between defense chemistry and onset of flowering. This variation likely reflects an adaptation to survival in different environments.

Direct and indirect selection on flowering time, water-use efficiency (WUE, δ (13)C), and WUE plasticity to drought in Arabidopsis thaliana

Flowering time and water-use efficiency (WUE) are two ecological traits that are important for plant drought response. To understand the evolutionary significance of natural genetic variation in flowering time, WUE, and WUE plasticity to drought in Arabidopsis thaliana, we addressed the following questions: (1) How are ecophysiological traits genetically correlated within and between different soil moisture environments? (2) Does terminal drought select for early flowering and drought escape? (3) Is WUE plasticity to drought adaptive and/or costly? We measured a suite of ecophysiological and reproductive traits on 234 spring flowering accessions of A. thaliana grown in well-watered and season-ending soil drying treatments, and quantified patterns of genetic variation, correlation, and selection within each treatment. WUE and flowering time were consistently positively genetically correlated. WUE was correlated with WUE plasticity, but the direction changed between treatments. Selection generally favored early flowering and low WUE, with drought favoring earlier flowering significantly more than well-watered conditions. Selection for lower WUE was marginally stronger under drought. There were no net fitness costs of WUE plasticity. WUE plasticity (per se) was globally neutral, but locally favored under drought. Strong genetic correlation between WUE and flowering time may facilitate the evolution of drought escape, or constrain independent evolution of these traits. Terminal drought favored drought escape in these spring flowering accessions of A. thaliana. WUE plasticity may be favored over completely fixed development in environments with periodic drought.

Identification of QTLs affecting scopolin and scopoletin biosynthesis in Arabidopsis thaliana

BACKGROUND

Scopoletin and its glucoside scopolin are important secondary metabolites synthesized in plants as a defense mechanism against various environmental stresses. They belong to coumarins, a class of phytochemicals with significant biological activities that is widely used in medical application and cosmetics industry. Although numerous studies showed that a variety of coumarins occurs naturally in several plant species, the details of coumarins biosynthesis and its regulation is not well understood. It was shown previously that coumarins (predominantly scopolin and scopoletin) occur in Arabidopsis thaliana (Arabidopsis) roots, but until now nothing is known about natural variation of their accumulation in this model plant. Therefore, the genetic architecture of coumarins biosynthesis in Arabidopsis has not been studied before.

RESULTS

Here, the variation in scopolin and scopoletin content was assessed by comparing seven Arabidopsis accessions. Subsequently, a quantitative trait locus (QTL) mapping was performed with an Advanced Intercross Recombinant Inbred Lines (AI-RILs) mapping population EstC (Est-1 × Col). In order to reveal the genetic basis of both scopolin and scopoletin biosynthesis, two sets of methanol extracts were made from Arabidopsis roots and one set was additionally subjected to enzymatic hydrolysis prior to quantification done by high-performance liquid chromatography (HPLC). We identified one QTL for scopolin and five QTLs for scopoletin accumulation. The identified QTLs explained 13.86% and 37.60% of the observed phenotypic variation in scopolin and scopoletin content, respectively. In silico analysis of genes located in the associated QTL intervals identified a number of possible candidate genes involved in coumarins biosynthesis.

CONCLUSIONS

Together, our results demonstrate for the first time that Arabidopsis is an excellent model for studying the genetic and molecular basis of natural variation in coumarins biosynthesis in plants. It additionally provides a basis for fine mapping and cloning of the genes involved in scopolin and scopoletin biosynthesis. Importantly, we have identified new loci for this biosynthetic process.

Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

Lagging adaptation to warming climate in Arabidopsis thaliana

If climate change outpaces the rate of adaptive evolution within a site, populations previously well adapted to local conditions may decline or disappear, and banked seeds from those populations will be unsuitable for restoring them. However, if such adaptational lag has occurred, immigrants from historically warmer climates will outperform natives and may provide genetic potential for evolutionary rescue. We tested for lagging adaptation to warming climate using banked seeds of the annual weed Arabidopsis thaliana in common garden experiments in four sites across the species' native European range: Valencia, Spain; Norwich, United Kingdom; Halle, Germany; and Oulu, Finland. Genotypes originating from geographic regions near the planting site had high relative fitness in each site, direct evidence for broad-scale geographic adaptation in this model species. However, genotypes originating in sites historically warmer than the planting site had higher average relative fitness than local genotypes in every site, especially at the northern range limit in Finland. This result suggests that local adaptive optima have shifted rapidly with recent warming across the species' native range. Climatic optima also differed among seasonal germination cohorts within the Norwich site, suggesting that populations occurring where summer germination is common may have greater evolutionary potential to persist under future warming. If adaptational lag has occurred over just a few decades in banked seeds of an annual species, it may be an important consideration for managing longer-lived species, as well as for attempts to conserve threatened populations through ex situ preservation.

Evolutionary genetics of plant adaptation: insights from new model systems

Flowering time and mating system divergence are two of the most common adaptive transitions in plants. We review recent progress toward understanding the genetic basis of these adaptations in new model plant species. For flowering time, we find that individual crosses often reveal a simple genetic basis, but that the loci involved almost always vary within species and across environments, indicating a more complex genetic basis species-wide. Similarly, the transition to self-fertilization is often genetically complex, but this seems to depend on the amount of standing variation and time since species divergence. Recent population genomic studies also raise doubts about the long-term adaptive potential of self-fertilization, providing evidence that purifying selection is less effective in highly selfing species.

Major quantitative trait loci control divergence in critical photoperiod for flowering between selfing and outcrossing species of monkeyflower (Mimulus)

• Divergence in flowering time is a key contributor to reproductive isolation between incipient species, as it enforces habitat specialization and causes assortative mating even in sympatry. Understanding the genetic basis of flowering time divergence illuminates the origins and maintenance of species barriers. • We investigated the genetics of divergence in critical photoperiod for flowering between yellow monkeyflowers Mimulus guttatus (outcrosser, summer flowering) and Mimulus nasutus (selfer, spring flowering). We used quantitative trait locus (QTL) mapping of F2 hybrids and fine-mapping in nearly isogenic lines to characterize the genomic regions underlying a > 2 h critical photoperiod difference between allopatric populations, and then tested whether the same QTLs control flowering time in sympatry. • We identified two major QTLs that almost completely explain M. nasutus's ability to flower in early spring; they are shared by allopatric and sympatric population pairs. The smaller QTL is coincident with one that differentiates ecotypes within M. guttatus, but the larger effect QTL appears unique to M. nasutus. • Unlike floral traits associated with mating system divergence, large interspecific differences in flowering phenology depend on only a few loci. Major critical photoperiod QTLs may be 'speciation genes' and also restrict interspecific gene flow in secondary sympatry.

Adaptation to spring heat and drought in northeastern Spanish Arabidopsis thaliana

The extent to which a species' environmental range reflects adaptive differentiation remains an open question. Environmental gradients can lead to adaptive divergence when differences in stressors among sites along the gradient place conflicting demands on the balance of stress responses. The extent to which this is accomplished through stress tolerance vs stress avoidance is also an open question. We present results from a controlled environment study of 48 lineages of Arabidopsis thaliana collected along a gradient in northeastern Spain across which temperatures increase and precipitation decreases with decreasing elevation. We tested the extent to which clinal adaptive divergence in heat and drought is explained through tolerance and avoidance traits by subjecting plants to a dynamic growth chamber cycle of increasing heat and drought stress analogous to low elevation spring in northeastern Spain. Lineages collected at low elevation were the most fit and fitness scaled with elevation of origin. Higher fitness was associated with earlier bolting, greater early allocation to increased numbers of inflorescences, reduction in rosette leaf photosynthesis and earlier fruit ripening. We propose that this is a syndrome of avoidance through early flowering accompanied by restructuring of the organism that adapts A. thaliana to low-elevation Mediterranean climates.

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana

Organisms inhabiting different environments are often locally adapted, and yet despite a considerable body of theory, the genetic basis of local adaptation is poorly understood. Unanswered questions include the number and effect sizes of adaptive loci, whether locally favored loci reduce fitness elsewhere (i.e., fitness tradeoffs), and whether a lack of genetic variation limits adaptation. To address these questions, we mapped quantitative trait loci (QTL) for total fitness in 398 recombinant inbred lines derived from a cross between locally adapted populations of the highly selfing plant Arabidopsis thaliana from Sweden and Italy and grown for 3 consecutive years at the parental sites (>40,000 plants monitored). We show that local adaptation is controlled by relatively few genomic regions of small to modest effect. A third of the 15 fitness QTL we detected showed evidence of tradeoffs, which contrasts with the minimal evidence for fitness tradeoffs found in previous studies. This difference may reflect the power of our multiyear study to distinguish conditionally neutral QTL from those that reflect fitness tradeoffs. In Sweden, but not in Italy, the local genotype underlying fitness QTL was often maladaptive, suggesting that adaptation there is constrained by a lack of adaptive genetic variation, attributable perhaps to genetic bottlenecks during postglacial colonization of Scandinavia or to recent changes in selection regime caused by climate change. Our results suggest that adaptation to markedly different environments can be achieved through changes in relatively few genomic regions, that fitness tradeoffs are common, and that lack of genetic variation can limit adaptation.