Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Alternative splicing variants of IiSEP3 in Isatis indigotica are involved in floral transition and flower development

The SEPALLATA3 genes regulate several aspects of plant development. This study identified four distinct splicing isoforms of the SEPALLATA3 gene in Isatis indigotica (I. indigotica). IiSEP3-1 and IiSEP3-2 have eight exons and were named as IiSEP3-2/1. However, IiSEP3-3 and IiSEP3-4 with the missing sixth exon were labeled IiSEP3ΔK3. Furthermore, the IiSEP3-1 and IiSEP3-4 amino acids sequences lack the V90. IiSEP3 splicing variants were primarily expressed in floral organs, with petals showing the highest expression. Ectopic expression of IiSEP3-2 or IiSEP3-3 may cause early flowering and reduce the number of sepals, petals, and stamens. The ectopic expression of IiSEP3-2 resulted in cauline leaves and sepals converting to carpelloid structures. In contrast, the four floral whorls prematurely wilted, and the entire flower displayed an abortive state when IiSEP3-3 was expressed ectopically. Silencing the IiSEP3 gene of I. indigotica employing VIGS (tobacco rattle virus-mediated virus-induced gene silencing) technology using the TRV-IiSEP3-2/1 vector delayed flowering time and reduced the number of petals and stamens. Plants silenced with TRV-IiSEP3ΔK3 also exhibited similar phenotypes, including fewer sepals. The transcriptome analysis of silenced plants (TRV-IiSEP3-2/1 treatment group) indicated significant alterations in 1861 genes, with 1035 upregulated and 826 downregulated. TRV-IiSEP3ΔK3 treatment altered the expression of 2063 genes in plants, with 1289 genes upregulated and 774 genes transcription inhibited. Y2H and BIFC experiments revealed that IiSEP3-2 and IiSEP3-3 had distinct interacting proteins. Thus, we can conclude that IiSEP3-2 and IiSEP3-3 interact with different proteins, affecting floral transition and organ development in I. indigotica.

Prediction of plant complex traits via integration of multi-omics data

The formation of complex traits is the consequence of genotype and activities at multiple molecular levels. However, connecting genotypes and these activities to complex traits remains challenging. Here, we investigate whether integrating genomic, transcriptomic, and methylomic data can improve prediction for six Arabidopsis traits. We find that transcriptome- and methylome-based models have performances comparable to those of genome-based models. However, models built for flowering time using different omics data identify different benchmark genes. Nine additional genes identified as important for flowering time from our models are experimentally validated as regulating flowering. Gene contributions to flowering time prediction are accession-dependent and distinct genes contribute to trait prediction in different genotypes. Models integrating multi-omics data perform best and reveal known and additional gene interactions, extending knowledge about existing regulatory networks underlying flowering time determination. These results demonstrate the feasibility of revealing molecular mechanisms underlying complex traits through multi-omics data integration.

Evolutionary Dynamics of FLC-like MADS-Box Genes in Brassicaceae

MADS-box genes encode transcription factors that play important roles in the development and evolution of plants. There are more than a dozen clades of MADS-box genes in angiosperms, of which those with functions in the specification of floral organ identity are especially well-known. From what has been elucidated in the model plant Arabidopsis thaliana, the clade of FLC-like MADS-box genes, comprising FLC-like genes sensu strictu and MAF-like genes, are somewhat special among the MADS-box genes of plants since FLC-like genes, especially MAF-like genes, show unusual evolutionary dynamics, in that they generate clusters of tandemly duplicated genes. Here, we make use of the latest genomic data of Brassicaceae to study this remarkable feature of the FLC-like genes in a phylogenetic context. We have identified all FLC-like genes in the genomes of 29 species of Brassicaceae and reconstructed the phylogeny of these genes employing a Maximum Likelihood method. In addition, we conducted selection analyses using PAML. Our results reveal that there are three major clades of FLC-like genes in Brassicaceae that all evolve under purifying selection but with remarkably different strengths. We confirm that the tandem arrangement of MAF-like genes in the genomes of Brassicaceae resulted in a high rate of duplications and losses. Interestingly, MAF-like genes also seem to be prone to transposition. Considering the role of FLC-like genes sensu lato (s.l.) in the timing of floral transition, we hypothesize that this rapid evolution of the MAF-like genes was a main contributor to the successful adaptation of Brassicaceae to different environments.

Structural variation of GL1 gene determines the trichome formation in Brassica juncea

A 448 kb region on chromosome B02 was delimited to be associated with trichome trait in Brassica juncea, in which the BjuVB02G54610 gene with a structural variation of 3 kb structure variation (SV) encoding a MYB transcription factor was predicted as the possible candidate gene. Mustards (Brassica juncea) are allopolyploid crops in the worldwide, and trichomes are essential quality attributes that significantly influence its taste and palpability in vegetable-use cultivars. As important accessory tissues from specialized epidermal cells, trichomes also play an important role in mitigating biotic and abiotic stresses. In this study, we constructed a F₂ segregating population using YJ27 with intensive trichome leaves and 03B0307 with glabrous leaves as parents. By bulked segregant analysis (BSA-seq), we obtained a 2.1 Mb candidate region on B02 chromosome associated with the trichome or glabrous trait formation. Then, we used 13 Kompetitive Allele Specific PCR (KASP) markers for fine mapping and finally narrowed down the candidate region to about 448 kb in length. Interestingly, among the region, there was a 3 kb sequence deletion that located on the BjuVB02G54610 gene in the F₂ individuals with trichome leaves. Genotyping results of F₂ populations confirmed this deletion (R² = 81.44%) as a major QTL. Natural population re-sequencing analysis and genotyping results further validated the key role of the 3 kb structure variation (SV) of insertion/deletion type in trichome development in B. juncea. Our findings provide important information on the formation of trichomes and potential target gene for breeding vegetable mustards.

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean (Glycine max) nested association mapping population

Photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500-2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

Transcriptome analysis and identification of genes associated with floral transition and fruit development in rabbiteye blueberry (Vaccinium ashei)

Flowering and fruit set are important traits affecting fruit quality and yield in rabbiteye blueberry (Vaccinium ashei). Intense efforts have been made to elucidate the influence of vernalization and phytohormones on flowering, but the molecular mechanisms of flowering and fruit set remain unclear. To unravel these mechanisms, we performed transcriptome analysis to explore blueberry transcripts from flowering to early fruit stage. We divided flowering and fruit set into flower bud (S2), initial flower (S3), bloom flower (S4), pad fruit (S5), and cup fruit (S6) based on phenotype and identified 1,344, 69, 658, and 189 unique differentially expressed genes (DEGs) in comparisons of S3/S2, S4/S3, S5/S4, and S6/S5, respectively. There were obviously more DEGs in S3/S2 and S5/S4 than in S4/S3, and S6/S5, suggesting that S3/S2 and S5/S4 represent major transitions from buds to fruit in blueberry. GO and KEGG enrichment analysis indicated these DEGs were mostly enriched in phytohormone biosynthesis and signaling, transporter proteins, photosynthesis, anthocyanins biosynthesis, disease resistance protein and transcription factor categories, in addition, transcript levels of phytohormones and transporters changed greatly throughout the flowering and fruit set process. Gibberellic acid and jasmonic acid mainly acted on the early stage of flowering development like expression of the florigen gene FT, while the expression of auxin response factor genes increased almost throughout the process from bud to fruit development. Transporter proteins were mainly associated with minerals during the early flowering development stage and sugars during the early fruit stage. At the early fruit stage, anthocyanins started to accumulate, and the fruit was susceptible to diseases such as fungal infection. Expression of the transcription factor MYB86 was up-regulated during initial fruit development, which may promote anthocyanin accumulation. These results will aid future studies exploring the molecular mechanism underlying flowering and fruit set of rabbiteye blueberry.

Transposition and duplication of MADS-domain transcription factor genes in annual and perennial Arabis species modulates flowering

Annual and perennial species differ in their timing and intensity of flowering, but the underlying mechanisms are poorly understood. We hybridized closely related annual and perennial plants and used genetics, transgenesis, and genomics to characterize differences in the activity and function of their flowering-time genes. We identify a gene encoding a transcription factor that moved between chromosomes and is retained in the annual but absent from the perennial. This gene strongly delays flowering, and we propose that it has been retained in the annual to compensate for reduced activity of closely related genes. This study highlights the value of using direct hybridization between closely related plant species to characterize functional differences in fast-evolving reproductive traits.

The timing of reproduction is an adaptive trait in many organisms. In plants, the timing, duration, and intensity of flowering differ between annual and perennial species. To identify interspecies variation in these traits, we studied introgression lines derived from hybridization of annual and perennial species, Arabis montbretiana and Arabis alpina, respectively. Recombination mapping identified two tandem A. montbretiana genes encoding MADS-domain transcription factors that confer extreme late flowering on A. alpina. These genes are related to the MADS AFFECTING FLOWERING (MAF) cluster of floral repressors of other Brassicaceae species and were named A. montbretiana (Am) MAF-RELATED (MAR) genes. AmMAR1 but not AmMAR2 prevented floral induction at the shoot apex of A. alpina, strongly enhancing the effect of the MAF cluster, and MAR1 is absent from the genomes of all A. alpina accessions analyzed. Exposure of plants to cold (vernalization) represses AmMAR1 transcription and overcomes its inhibition of flowering. Assembly of the tandem arrays of MAR and MAF genes of six A. alpina accessions and three related species using PacBio long-sequence reads demonstrated that the MARs arose within the Arabis genus by interchromosomal transposition of a MAF1-like gene followed by tandem duplication. Time-resolved comparative RNA-sequencing (RNA-seq) suggested that AmMAR1 may be retained in A. montbretiana to enhance the effect of the AmMAF cluster and extend the duration of vernalization required for flowering. Our results demonstrate that MAF genes transposed independently in different Brassicaceae lineages and suggest that they were retained to modulate adaptive flowering responses that differ even among closely related species.

Temperature-Dependent Alternative Splicing of Precursor mRNAs and Its Biological Significance: A Review Focused on Post-Transcriptional Regulation of a Cold Shock Protein Gene in Hibernating Mammals

Multiple mRNA isoforms are often generated during processing such as alternative splicing of precursor mRNAs (pre-mRNA), resulting in a diversity of generated proteins. Alternative splicing is an essential mechanism for the functional complexity of eukaryotes. Temperature, which is involved in all life activities at various levels, is one of regulatory factors for controlling patterns of alternative splicing. Temperature-dependent alternative splicing is associated with various phenotypes such as flowering and circadian clock in plants and sex determination in poikilothermic animals. In some specific situations, temperature-dependent alternative splicing can be evoked even in homothermal animals. For example, the splicing pattern of mRNA for a cold shock protein, cold-inducible RNA-binding protein (CIRP or CIRBP), is changed in response to a marked drop in body temperature during hibernation of hamsters. In this review, we describe the current knowledge about mechanisms and functions of temperature-dependent alternative splicing in plants and animals. Then we discuss the physiological significance of hypothermia-induced alternative splicing of a cold shock protein gene in hibernating and non-hibernating animals.

Evolutionary processes from the perspective of flowering time diversity

Although it is well appreciated that genetic studies of flowering time regulation have led to fundamental advances in the fields of molecular and developmental biology, the ways in which genetic studies of flowering time diversity have enriched the field of evolutionary biology have received less attention despite often being equally profound. Because flowering time is a complex, environmentally responsive trait that has critical impacts on plant fitness, crop yield, and reproductive isolation, research into the genetic architecture and molecular basis of its evolution continues to yield novel insights into our understanding of domestication, adaptation, and speciation. For instance, recent studies of flowering time variation have reconstructed how, when, and where polygenic evolution of phenotypic plasticity proceeded from standing variation and de novo mutations; shown how antagonistic pleiotropy and temporally varying selection maintain polymorphisms in natural populations; and provided important case studies of how assortative mating can evolve and facilitate speciation with gene flow. In addition, functional studies have built detailed regulatory networks for this trait in diverse taxa, leading to new knowledge about how and why developmental pathways are rewired and elaborated through evolutionary time.

Enhanced Maps of Transcription Factor Binding Sites Improve Regulatory Networks Learned from Accessible Chromatin Data

Determining where transcription factors (TFs) bind in genomes provides insight into which transcriptional programs are active across organs, tissue types, and environmental conditions. Recent advances in high-throughput profiling of regulatory DNA have yielded large amounts of information about chromatin accessibility. Interpreting the functional significance of these data sets requires knowledge of which regulators are likely to bind these regions. This can be achieved by using information about TF-binding preferences, or motifs, to identify TF-binding events that are likely to be functional. Although different approaches exist to map motifs to DNA sequences, a systematic evaluation of these tools in plants is missing. Here, we compare four motif-mapping tools widely used in the Arabidopsis (Arabidopsis thaliana) research community and evaluate their performance using chromatin immunoprecipitation data sets for 40 TFs. Downstream gene regulatory network (GRN) reconstruction was found to be sensitive to the motif mapper used. We further show that the low recall of Find Individual Motif Occurrences, one of the most frequently used motif-mapping tools, can be overcome by using an Ensemble approach, which combines results from different mapping tools. Several examples are provided demonstrating how the Ensemble approach extends our view on transcriptional control for TFs active in different biological processes. Finally, a protocol is presented to effectively derive more complete cell type-specific GRNs through the integrative analysis of open chromatin regions, known binding site information, and expression data sets. This approach will pave the way to increase our understanding of GRNs in different cellular conditions.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

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

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.

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.

Population genomic scans suggest novel genes underlie convergent flowering time evolution in the introduced range of Arabidopsis thaliana

A long-standing question in evolutionary biology is whether the evolution of convergent phenotypes results from selection on the same heritable genetic components. Using whole-genome sequencing and genome scans, we tested whether the evolution of parallel longitudinal flowering time clines in the native and introduced ranges of Arabidopsis thaliana has a similar genetic basis. We found that common variants of large effect on flowering time in the native range do not appear to have been under recent strong selection in the introduced range. We identified a set of 38 new candidate genes that are putatively linked to the evolution of flowering time. A high degree of conditional neutrality of flowering time variants between the native and introduced range may preclude parallel evolution at the level of genes. Overall, neither gene pleiotropy nor available standing genetic variation appears to have restricted the evolution of flowering time to high-frequency variants from the native range or to known flowering time pathway genes.

Genome-Wide Association Mapping and Genomic Prediction Elucidate the Genetic Architecture of Morphological Traits in Arabidopsis

Quantitative traits in plants are controlled by a large number of genes and their interaction with the environment. To disentangle the genetic architecture of such traits, natural variation within species can be explored by studying genotype-phenotype relationships. Genome-wide association studies that link phenotypes to thousands of single nucleotide polymorphism markers are nowadays common practice for such analyses. In many cases, however, the identified individual loci cannot fully explain the heritability estimates, suggesting missing heritability. We analyzed 349 Arabidopsis accessions and found extensive variation and high heritabilities for different morphological traits. The number of significant genome-wide associations was, however, very low. The application of genomic prediction models that take into account the effects of all individual loci may greatly enhance the elucidation of the genetic architecture of quantitative traits in plants. Here, genomic prediction models revealed different genetic architectures for the morphological traits. Integrating genomic prediction and association mapping enabled the assignment of many plausible candidate genes explaining the observed variation. These genes were analyzed for functional and sequence diversity, and good indications that natural allelic variation in many of these genes contributes to phenotypic variation were obtained. For ACS11, an ethylene biosynthesis gene, haplotype differences explaining variation in the ratio of petiole and leaf length could be identified.

Environmental and genetic interactions reveal FLOWERING LOCUS C as a modulator of the natural variation for the plasticity of flowering in Arabidopsis

The timing of flowering initiation depends strongly on the environment, a property termed as the plasticity of flowering. Such plasticity determines the adaptive potential of plants because it provides phenotypic buffer against environmental changes, and its natural variation contributes to evolutionary adaptation. We addressed the genetic mechanisms of the natural variation for this plasticity in Arabidopsis thaliana by analysing a population of recombinant inbred lines derived from Don-0 and Ler accessions collected from distinct climates. Quantitative trait locus (QTL) mapping in four environmental conditions differing in photoperiod, vernalization treatment and ambient temperature detected the folllowing: (i) FLOWERING LOCUS C (FLC) as a large effect QTL affecting flowering time differentially in all environments; (ii) numerous QTL displaying smaller effects specifically in some conditions; and (iii) significant genetic interactions between FLC and other loci. Hence, the variation for the plasticity of flowering is determined by a combination of environmentally sensitive and specific QTL, and epistasis. Analysis of FLC from Don identified a new and more active allele likely caused by a cis-regulatory deletion covering the non-coding RNA COLDAIR. Further characterization of four FLC natural alleles showed different environmental and genetic interactions. Thus, FLC appears as a major modulator of the natural variation for the plasticity of flowering to multiple environmental factors.

Role of alternative pre-mRNA splicing in temperature signaling

Developmental plasticity enables plants to respond rapidly to changing environmental conditions, such as temperature fluctuations. Understanding how plants measure temperature and integrate this information into developmental programs at the molecular level will be essential to breed thermo-tolerant crop varieties. Recent studies identified alternative splicing (AS) as a possible 'molecular thermometer', allowing plants to quickly adjust the abundance of functional transcripts to environmental perturbations. In this review, recent advances regarding the effects of temperature-responsive AS on plant development will be discussed, with emphasis on the circadian clock and flowering time control. The challenge for the near future will be to understand the molecular mechanisms by which temperature can influence AS regulation.

MAF2 Is Regulated by Temperature-Dependent Splicing and Represses Flowering at Low Temperatures in Parallel with FLM

Plants enter their reproductive phase when the environmental conditions are favourable for the successful production of progeny. The transition from vegetative to reproductive phase is influenced by several environmental factors including ambient temperature. In the model plant Arabidopsis thaliana, SHORT VEGETATIVE PHASE (SVP) is critical for this pathway; svp mutants cannot modify their flowering time in response to ambient temperature. SVP encodes a MADS-box transcription factor that directly represses genes that promote flowering. SVP binds DNA in complexes with other MADS-box transcription factors, including FLOWERING LOCUS M (FLM), which acts with SVP to repress the floral transition at low temperatures. Small temperature changes post-transcriptionally regulate FLM through temperature-dependent alternative splicing (TD-AS). As ambient temperature increases, the predominant FLM splice isoform shifts to encode a protein incapable of exerting a repressive effect on flowering. Here we characterize a closely related MADS-box transcription factor, MADS AFFECTING FLOWERING2 (MAF2), which has independently evolved TD-AS. At low temperatures the most abundant MAF2 splice variant encodes a protein that interacts with SVP to repress flowering. At increased temperature the relative abundance of splice isoforms shifts in favour of an intron-retaining variant that introduces a premature termination codon. We show that this isoform encodes a protein that cannot interact with SVP or repress flowering. At lower temperatures MAF2 and SVP repress flowering in parallel with FLM and SVP, providing an additional input to sense ambient temperature for the control of flowering.

Genetic architecture of variation in heading date among Asian rice accessions

BACKGROUND

Heading date, a crucial factor determining regional and seasonal adaptation in rice (Oryza sativa L.), has been a major selection target in breeding programs. Although considerable progress has been made in our understanding of the molecular regulation of heading date in rice during last two decades, the previously isolated genes and identified quantitative trait loci (QTLs) cannot fully explain the natural variation for heading date in diverse rice accessions.

RESULTS

To genetically dissect naturally occurring variation in rice heading date, we collected QTLs in advanced-backcross populations derived from multiple crosses of the japonica rice accession Koshihikari (as a common parental line) with 11 diverse rice accessions (5 indica, 3 aus, and 3 japonica) that originate from various regions of Asia. QTL analyses of over 14,000 backcrossed individuals revealed 255 QTLs distributed widely across the rice genome. Among the detected QTLs, 128 QTLs corresponded to genomic positions of heading date genes identified by previous studies, such as Hd1, Hd6, Hd3a, Ghd7, DTH8, and RFT1. The other 127 QTLs were detected in different chromosomal regions than heading date genes.

CONCLUSIONS

Our results indicate that advanced-backcross progeny allowed us to detect and confirm QTLs with relatively small additive effects, and the natural variation in rice heading date could result from combinations of large- and small-effect QTLs. We also found differences in the genetic architecture of heading date (flowering time) among maize, Arabidopsis, and rice.

Natural variation for anthocyanin accumulation under high-light and low-temperature stress is attributable to the ENHANCER OF AG-4 2 (HUA2) locus in combination with PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2

Growing conditions combining high light intensities and low temperatures lead to anthocyanin accumulation in plants. This response was contrasted between two Arabidopsis thaliana accessions, which were used to decipher the genetic and molecular bases underlying the variation of this response. Quantitative trait loci (QTLs) for flowering time (FT) and anthocyanin accumulation under a high-light and low-temperature scenario versus a control environment were mapped. Major QTLs were confirmed using near-isogenic lines. Candidate genes were examined using mutants and gene expression studies as well as transgenic complementation. Several QTLs were found for FT and for anthocyanin content, of which one QTL co-located at the ENHANCER OF AG-4 2 (HUA2) locus. That HUA2 is a regulator of both pathways was confirmed by the analysis of loss-of-function mutants. For a strong expression of anthocyanin, additional allelic variation was detected for the PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2 genes which control the anthocyanin pathway. The genetic control of variation for anthocyanin content was dissected in A. thaliana and shown to be affected by a common regulator of flowering and anthocyanin biosynthesis together with anthocyanin-specific regulators.

Gene regulatory variation mediates flowering responses to vernalization along an altitudinal gradient in Arabidopsis

Steep environmental gradients provide ideal settings for studies of potentially adaptive phenotypic and genetic variation in plants. The accurate timing of flowering is crucial for reproductive success and is regulated by several pathways, including the vernalization pathway. Among the numerous genes known to enable flowering in response to vernalization, the most prominent is FLOWERING LOCUS C (FLC). FLC and other genes of the vernalization pathway vary extensively among natural populations and are thus candidates for the adaptation of flowering time to environmental gradients such as altitude. We used 15 natural Arabidopsis (Arabidopsis thaliana) genotypes originating from an altitudinal gradient (800-2,700 m above sea level) in the Swiss Alps to test whether flowering time correlated with altitude under different vernalization scenarios. Additionally, we measured the expression of 12 genes of the vernalization pathway and its downstream targets. Flowering time correlated with altitude in a nonlinear manner for vernalized plants. Flowering time could be explained by the expression and regulation of the vernalization pathway, most notably by AGAMOUS LIKE19 (AGL19), FLOWERING LOCUS T (FT), and FLC. The expression of AGL19, FT, and VERNALIZATION INSENSITIVE3 was associated with altitude, and the regulation of MADS AFFECTING FLOWERING2 (MAF2) and MAF3 differed between low- and high-altitude genotypes. In conclusion, we found clinal variation across an altitudinal gradient both in flowering time and the expression and regulation of genes in the flowering time control network, often independent of FLC, suggesting that the timing of flowering may contribute to altitudinal adaptation.

QTL meta-analysis in Arabidopsis reveals an interaction between leaf senescence and resource allocation to seeds

Sequential and monocarpic senescence are observed at vegetative and reproductive stages, respectively. Both facilitate nitrogen (N) remobilization and control the duration of carbon (C) fixation. Genetic and environmental factors control N and C resource allocation to seeds. Studies of natural variation in Arabidopsis thaliana revealed differences between accessions for leaf senescence phenotypes, seed N and C contents, and N remobilization efficiency-related traits. Here, a quantitative genetics approach was used to gain a better understanding of seed filling regulation in relation to leaf senescence and resource allocation. For that purpose, three Arabidopsis recombinant inbred line populations (Ct-1×Col-0, Cvi-0×Col-0, Bur-0×Col-0) were used to map QTL (quantitative trait loci) for ten traits related to senescence, resource allocation, and seed filling. The use of common markers across the three different maps allowed direct comparisons of the positions of the detected QTL in a single consensus map. QTL meta-analysis was then used to identify interesting regions (metaQTL) where QTL for several traits co-localized. MetaQTL were compared with positions of candidate genes known to be involved in senescence processes and flowering time. Finally, investigation of the correlation between yield and seed N concentration in the three populations suggests that leaf senescence disrupts the negative correlation generally observed between these two traits.

Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesis

Deciphering the genetic and molecular bases of quantitative variation is a long-standing challenge in plant biology because it is essential for understanding evolution and for accelerating plant breeding. Recent multi-trait analyses at different phenotypic levels are uncovering the pleiotropy and the genetic regulation underlying high-level complex traits. Thus, the number of known causal loci, genes and nucleotide polymorphisms is expanding. Current plant causal catalogs contain ∼400 genes and natural polymorphisms revealing several dysfunctional allelic series that involve multiple mutations. In addition, repeated evolution of quantitative traits mediated by large effect alleles is found across plant phylogeny. Finally, systematic analyses of genetic and environmental interactions are beginning to elucidate the molecular mechanisms of relevant interactions.

Structural variation and genome complexity: is dispensable really dispensable?

Structural variants (SVs) such as copy number variants (CNVs) and presence/absence variants (PAVs) substantially contribute to genetic variation and have an important effect on phenotypic diversity. Since unbalanced SVs are by definition sequences present only in some individuals, they have therefore been referred to as dispensable genome and are not necessary for survival, even though they may provide an important contribution to phenotypic diversity within the species. However, some multi-copy sequences of the dispensable genomes (e.g., multigene families) may be needed in a given proportion by each individual, thus belonging to a conditionally dispensable portion of the pan-genome. Another interesting aspect reported by recent studies is that the rate at which SVs are formed might be influenced by the mating system and by common environmental stresses. In conclusion the dispensable genome plays an important role in genome evolution and in the complex interplay between the genome and the environment.

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.

Acceleration of flowering in Arabidopsis thaliana by Cape Verde Islands alleles of FLOWERING H is dependent on the floral promoter FD

Flowering time in the model plant Arabidopsis thaliana is regulated by both external environmental signals and internal developmental pathways. Natural variation at the FLOWERING H (FLH) locus has previously been described, with alleles present in the Cape Verde Islands accession causing early flowering, particularly after vernalization. The mechanism of FLH-induced early flowering is not understood. Here, the integration of FLH activity into the known flowering time pathways is described using molecular and genetic approaches. The identification of molecular markers that co-segregated with the FLH locus allowed the generation of multiple combinations of FLH alleles with mutations in flowering time genes in different flowering pathways. Combining an early flowering FLH allele with mutations in vernalization pathway genes that regulate FLC expression revealed that FLH appears to act in parallel to FLC. Surprisingly, the early flowering allele of FLH requires the floral integrator FD, but not FT, to accelerate flowering. This suggests a model in which some alleles of FLH are able to affect the FD-dependent activity of the floral activator complex.

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

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

The flowering repressor SVP underlies a novel Arabidopsis thaliana QTL interacting with the genetic background

The timing of flowering initiation is a fundamental trait for the adaptation of annual plants to different environments. Large amounts of intraspecific quantitative variation have been described for it among natural accessions of many species, but the molecular and evolutionary mechanisms underlying this genetic variation are mainly being determined in the model plant Arabidopsis thaliana. To find novel A. thaliana flowering QTL, we developed introgression lines from the Japanese accession Fuk, which was selected based on the substantial transgression observed in an F₂ population with the reference strain Ler. Analysis of an early flowering line carrying a single Fuk introgression identified Flowering Arabidopsis QTL1 (FAQ1). We fine-mapped FAQ1 in an 11 kb genomic region containing the MADS transcription factor gene SHORT VEGETATIVE PHASE (SVP). Complementation of the early flowering phenotype of FAQ1-Fuk with a SVP-Ler transgen demonstrated that FAQ1 is SVP. We further proved by directed mutagenesis and transgenesis that a single amino acid substitution in SVP causes the loss-of-function and early flowering of Fuk allele. Analysis of a worldwide collection of accessions detected FAQ1/SVP-Fuk allele only in Asia, with the highest frequency appearing in Japan, where we could also detect a potential ancestral genotype of FAQ1/SVP-Fuk. In addition, we evaluated allelic and epistatic interactions of SVP natural alleles by analysing more than one hundred transgenic lines carrying Ler or Fuk SVP alleles in five genetic backgrounds. Quantitative analyses of these lines showed that FAQ1/SVP effects vary from large to small depending on the genetic background. These results support that the flowering repressor SVP has been recently selected in A. thaliana as a target for early flowering, and evidence the relevance of genetic interactions for the intraspecific evolution of FAQ1/SVP and flowering time.

In many plant species, the timing of flowering initiation shows abundant quantitative variation among natural varieties, which reflects the importance of this trait for adaptation to different environments. Currently, a major goal in plant biology is to determine the molecular and evolutionary bases of this natural genetic variation. In this study we demonstrate that the central flowering regulator SHORT VEGETATIVE PHASE (SVP), encoding a MADS transcription factor, is involved in the flowering natural variation of the model organism Arabidopsis thaliana. In particular, we prove that a structural change caused by a single amino acid substitution generates a SVP early flowering allele that is distributed only in Asia. Furthermore, genetic interactions have been shown to be a component of the natural variation for many important adaptive traits. However, very few studies, either in animals or plants, have systematically addressed the extent of genetic interactions among specific alleles responsible for the natural variation of complex traits. Our study shows that the flowering effects of SVP natural alleles depend significantly on the genetic background; and, subsequently, we demonstrate the relevance of epistasis for the evolution of this crucial transcription factor and flowering time.

Functional analysis of splice variant expression of MADS AFFECTING FLOWERING 2 of Arabidopsis thaliana

The MADS-AFFECTING FLOWERING 2 (MAF2) gene of Arabidopsis thaliana has been characterized as a repressor of flowering. The molecular basis of MAF2 gene function and role of alternative MAF2 transcripts in flowering time modulation is not understood. MAF2 splice variant expression was quantified in cold-acclimated plants by quantitative RT-PCR. Cold influenced the abundance of splice variants and prompted a functional study of splice forms. Individual variants were overexpressed in the Col background and were assayed for their ability to delay flowering. Overexpression of MAF2 variants 2 and 4 had limited effect on flowering time. Overexpression of MAF2 splice variant 1 resulted in early flowering and affected the expression of the endogenous MAF2 gene and its paralogues, confounding functional assessment. In the Ll-2 Arabidopsis accession, a MAF2, MAF3, MAF4 and FLC null line, MAF2 var1 was consistent in its effect on reproductive delay under ambient and reduced temperatures, indicating that it acts as a repressor of flowering.

Divergence of flowering genes in soybean

Soybean genome sequences were blasted with Arabidopsis thaliana regulatory genes involved in photoperioddependent flowering. This approach enabled the identification of 118 genes involved in the flowering pathway. Two genome sequences of cultivated (Williams 82) and wild (IT182932) soybeans were employed to survey functional DNA variations in the flowering-related homologs. Forty genes exhibiting nonsynonymous substitutions between G. max and G. soja were catalogued. In addition, 22 genes were found to co-localize with QTLs for six traits including flowering time, first flower, pod maturity, beginning of pod, reproductive period, and seed filling period. Among the genes overlapping the QTL regions, two LHY/CCA1 genes, GI and SFR6 contained amino acid changes. The recently duplicated sequence regions of the soybean genome were used as additional criteria for the speculation of the putative function of the homologs. Two duplicated regions showed redundancy of both flowering-related genes and QTLs. ID 12398025, which contains the homeologous regions between chr 7 and chr 16, was redundant for the LHY/CCA1 and SPA1 homologs and the QTLs. Retaining of the CRY1 gene and the pod maturity QTLs were observed in the duplicated region of ID 23546507 on chr 4 and chr 6. Functional DNA variation of the LHY/CCA1 gene (Glyma07g05410) was present in a counterpart of the duplicated region on chr 7, while the gene (Glyma16g01980) present in the other portion of the duplicated region on chr 16 did not show a functional sequence change. The gene list catalogued in this study provides primary insight for understanding the regulation of flowering time and maturity in soybean.

Genetic architecture of flowering-time variation in Arabidopsis thaliana

The onset of flowering is an important adaptive trait in plants. The small ephemeral species Arabidopsis thaliana grows under a wide range of temperature and day-length conditions across much of the Northern hemisphere, and a number of flowering-time loci that vary between different accessions have been identified before. However, only few studies have addressed the species-wide genetic architecture of flowering-time control. We have taken advantage of a set of 18 distinct accessions that present much of the common genetic diversity of A. thaliana and mapped quantitative trait loci (QTL) for flowering time in 17 F₂ populations derived from these parents. We found that the majority of flowering-time QTL cluster in as few as five genomic regions, which include the locations of the entire FLC/MAF clade of transcription factor genes. By comparing effects across shared parents, we conclude that in several cases there might be an allelic series caused by rare alleles. While this finding parallels results obtained for maize, in contrast to maize much of the variation in flowering time in A. thaliana appears to be due to large-effect alleles.