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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

Strategies for breeding crops for future environments

The green revolution successfully increased agricultural output in the early 1960s by relying primarily on three pillars: plant breeding, irrigation, and chemical fertilization. Today, the need to reduce the use of chemical fertilizers, water scarcity, and future environmental changes, together with a growing population, requires innovative strategies to adapt to a new context and prevent food shortages. Therefore, scientists from around the world are directing their efforts to breed crops for future environments to sustainably produce more nutritious food. Herein, we propose scientific avenues to be reinforced in selecting varieties, including crop wild relatives, either for monoculture or mixed cropping systems, taking advantage of plant-microbial interactions, while considering the diversity of organisms associated with crops and unlocking combinatorial nutritional stresses.

Genetic structure and ecological niche space of lentil's closest wild relative, Lens orientalis (Boiss.) Schmalh

Crops arose from wild ancestors and to understand their domestication it is essential to compare the cultivated species with their crop wild relatives. These represent an important source of further crop improvement, in particular in relation to climate change. Although there are about 58,000 Lens accessions held in genebanks, only 1% are wild. We examined the geographic distribution and genetic diversity of the lentil's immediate progenitor L. orientalis. We used Genotyping by Sequencing (GBS) to identify and characterize differentiation among accessions held at germplasm collections. We then determined whether genetically distinct clusters of accessions had been collected from climatically distinct locations. Of the 195 genotyped accessions, 124 were genuine L. orientalis with four identified genetic groups. Although an environmental distance matrix was significantly correlated with geographic distance in a Mantel test, the four identified genetic clusters were not found to occupy significantly different environmental space. Maxent modelling gave a distinct predicted distribution pattern centred in the Fertile Crescent, with intermediate probabilities of occurrence in parts of Turkey, Greece, Cyprus, Morocco, and the south of the Iberian Peninsula with NW Africa. Future projections did not show any dramatic alterations in the distribution according to the climate change scenarios tested. We have found considerable diversity in L. orientalis, some of which track climatic variability. The results of the study showed the genetic diversity of wild lentil and indicate the importance of ongoing collections and in situ conservation for our future capacity to harness the genetic variation of the lentil progenitor.

The Roles of the PSEUDO-RESPONSE REGULATORs in Circadian Clock and Flowering Time in Medicago truncatula

PSEUDO-RESPONSE REGULATORs (PRRs) play key roles in the circadian rhythms and flowering in plants. Here, we identified the four members of the PRR family in Medicago truncatula, including MtPRR9a, MtPRR9b, MtPRR7 and MtPRR5, and isolated their Tnt1 retrotransposon-tagged mutants. They were expressed in different organs and were nuclear-localized. The four MtPRRs genes played important roles in normal clock rhythmicity maintenance by negatively regulating the expression of MtGI and MtLHY. Surprisingly, the four MtPRRs functioned redundantly in regulating flowering time under long-day conditions, and the quadruple mutant flowered earlier. Moreover, MtPRR can recruit the MtTPL/MtTPR corepressors and the other MtPRRs to form heterodimers to constitute the core mechanism of the circadian oscillator.

Genetic diversity and local adaption of alfalfa populations (Medicago sativa L.) under long-term grazing

Genomic information on alfalfa adaptation to long-term grazing is useful for alfalfa genetic improvement. In this study, 14 alfalfa populations were collected from long-term grazing sites (> 25 years) across four soil zones in western Canada. Alfalfa cultivars released between 1926 and 1980 were used to compare degree of genetic variation of the 14 populations. Six agro-morphological and three nutritive value traits were evaluated from 2018 to 2020. The genotyping-by-sequencing (GBS) data of the alfalfa populations and environmental data were used for genotype-environment association (GEA). Both STRUCTURE and UPGMA based on 19,853 SNPs showed that the 14 alfalfa populations from long-term grazing sites had varying levels of parentages from alfalfa sub-species Medicago sativa and M. falcata. The linear regression of STRUCTURE membership probability on phenotypic data indicated genetic variations of forage dry matter yield, spring vigor and plant height were low, but genetic variations of regrowth, fall plant height, days to flower and crude protein were still high for the 14 alfalfa populations from long-term grazing sites. The GEA identified 31 SNPs associated with 13 candidate genes that were mainly associated with six environmental factors of. Candidate genes underlying environmental factors were associated with a variety of proteins, which were involved in plant responses to abiotic stresses, i.e., drought, cold and salinity-alkali stresses.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift

Resurrection studies are a useful tool to measure how phenotypic traits have changed in populations through time. If these trait modifications correlate with the environmental changes that occurred during the time period, it suggests that the phenotypic changes could be a response to selection. Selfing, through its reduction of effective size, could challenge the ability of a population to adapt to environmental changes. Here, we used a resurrection study to test for adaptation in a selfing population of Medicago truncatula, by comparing the genetic composition and flowering times across 22 generations. We found evidence for evolution toward earlier flowering times by about two days and a peculiar genetic structure, typical of highly selfing populations, where some multilocus genotypes (MLGs) are persistent through time. We used the change in frequency of the MLGs through time as a multilocus fitness measure and built a selection gradient that suggests evolution toward earlier flowering times. Yet, a simulation model revealed that the observed change in flowering time could be explained by drift alone, provided the effective size of the population is small enough (<150). These analyses suffer from the difficulty to estimate the effective size in a highly selfing population, where effective recombination is severely reduced.

Annual and perennial Medicago show signatures of parallel adaptation to climate and soil in highly conserved genes

Human induced environmental change may require rapid adaptation of plant populations and crops, but the genomic basis of environmental adaptation remain poorly understood. We analysed polymorphic loci from the perennial crop Medicago sativa (alfalfa or lucerne) and the annual legume model species M. truncatula to search for a common set of candidate genes that might contribute to adaptation to abiotic stress in both annual and perennial Medicago species. We identified a set of candidate genes of adaptation associated with environmental gradients along the distribution of the two Medicago species. Candidate genes for each species were detected in homologous genomic linkage blocks using genome-environment (GEA) and genome-phenotype association analyses. Hundreds of GEA candidate genes were species-specific, of these, 13.4% (M. sativa) and 24% (M. truncatula) were also significantly associated with phenotypic traits. A set of 168 GEA candidates were shared by both species, which was 25.4% more than expected by chance. When combined, they explained a high proportion of variance for certain phenotypic traits associated with adaptation. Genes with highly conserved functions dominated among the shared candidates and were enriched in gene ontology terms that have shown to play a central role in drought avoidance and tolerance mechanisms by means of cellular shape modifications and other functions associated with cell homeostasis. Our results point to the existence of a molecular basis of adaptation to abiotic stress in Medicago determined by highly conserved genes and gene functions. We discuss these results in light of the recently proposed omnigenic model of complex traits.

Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing

Elymus sibiricus L. is an important cold-season grass with excellent cold and drought tolerance, good palatability, and nutrition. Flowering time is a key trait that affects forage and seed yield. Development of EST-SSR (expressed sequence tag simple sequence repeat) markers based on flowering genes contributes to the improvement of flowering traits. In the study, we detected 155 candidate genes related to flowering traits from 10,591 unigenes via transcriptome sequencing in early- and late-flowering genotypes. These candidate genes were mainly involved in the photoperiodic pathway, vernalization pathway, central integrator, and gibberellin pathway. A total of 125 candidate gene-based EST-SSRs were developed. Further, 15 polymorphic EST-SSRs closely associated to 13 candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions, including two contrasting panels (early-flowering and late-flowering). Among them, primer 28366, designed from heading date 3a (HD3a), effectively distinguished early- and late-flowering genotypes using a specifically amplified band of 175 bp. The polymorphic information content (PIC) value ranged from 0.12 to 0.48, with an average of 0.25. The unweighted pair group method analysis (UPGMA) cluster and structure analysis showed that the 20 E. sibiricus genotypes with similar flowering times tended to group together. These newly developed EST-SSR markers have the potential to be used for molecular markers assisted selection and germplasm evaluation of flowering traits in E. sibiricus.

A linkage disequilibrium-based statistical test for Genome-Wide Epistatic Selection Scans in structured populations

The quest for signatures of selection using single nucleotide polymorphism (SNP) data has proven efficient to uncover genes involved in conserved and/or adaptive molecular functions, but none of the statistical methods were designed to identify interacting alleles as targets of selective processes. Here, we propose a statistical test aimed at detecting epistatic selection, based on a linkage disequilibrium (LD) measure accounting for population structure and heterogeneous relatedness between individuals. SNP-based ([Formula: see text]) and window-based ([Formula: see text]) statistics fit a Student distribution, allowing to test the significance of correlation coefficients. As a proof of concept, we use SNP data from the Medicago truncatula symbiotic legume plant and uncover a previously unknown gene coadaptation between the MtSUNN (Super Numeric Nodule) receptor and the MtCLE02 (CLAVATA3-Like) signaling peptide. We also provide experimental evidence supporting a MtSUNN-dependent negative role of MtCLE02 in symbiotic root nodulation. Using human HGDP-CEPH SNP data, our new statistical test uncovers strong LD between SLC24A5 (skin pigmentation) and EDAR (hairs, teeth, sweat glands development) world-wide, which persists after correction for population structure and relatedness in Central South Asian populations. This result suggests that epistatic selection or coselection could have contributed to the phenotypic make-up in some human populations. Applying this approach to genome-wide SNP data will facilitate the identification of coadapted gene networks in model or non-model organisms.

Physical Dormancy Release in Medicago truncatula Seeds Is Related to Environmental Variations

Seed dormancy and timing of its release is an important developmental transition determining the survival of individuals, populations, and species in variable environments. Medicago truncatula was used as a model to study physical seed dormancy at the ecological and genetics level. The effect of alternating temperatures, as one of the causes releasing physical seed dormancy, was tested in 178 M. truncatula accessions over three years. Several coefficients of dormancy release were related to environmental variables. Dormancy varied greatly (4-100%) across accessions as well as year of experiment. We observed overall higher physical dormancy release under more alternating temperatures (35/15 °C) in comparison with less alternating ones (25/15 °C). Accessions from more arid climates released dormancy under higher experimental temperature alternations more than accessions originating from less arid environments. The plasticity of physical dormancy can probably distribute the germination through the year and act as a bet-hedging strategy in arid environments. On the other hand, a slight increase in physical dormancy was observed in accessions from environments with higher among-season temperature variation. Genome-wide association analysis identified 136 candidate genes related to secondary metabolite synthesis, hormone regulation, and modification of the cell wall. The activity of these genes might mediate seed coat permeability and, ultimately, imbibition and germination.

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS® . Based on 99 085 high-quality SNPs, landraces were classified into three sub-populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour-joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait-associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub-populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.

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.

Temperature- versus precipitation-limitation shape local temperature tolerance in a Holarctic freshwater crustacean

Species with wide geographical distributions are often adapted locally to the prevailing temperatures. To understand how they respond to ongoing climatic change, we must appreciate the interplay between temperature, seasonality and the organism's life cycle. The temperature experienced by many organisms results from an often-overlooked combination of climate and phenology. Summer-active (high latitude) populations are expected to adapt to local summer temperatures, but this is not expected for populations that outlive the summer in their dormant stage (low latitude, precipitation-limited). We recorded reproduction and survival in genotypes from 123 Holarctic populations of Daphnia magna during a multi-generation thermal ramp experiment. Genotypes from summer-active populations showed a positive relationship between heat tolerance and local summer temperature, whereas winter-active populations did not. These findings are consistent with the hypothesis that D. magna adapts to the local temperatures the animals experience during their planktonic phase. We conclude that predicting local temperature adaptation, in particular in the light of climate change, needs to consider the phenology of geographically wide-ranging species.

A local score approach improves GWAS resolution and detects minor QTL: application to Medicago truncatula quantitative disease resistance to multiple Aphanomyces euteiches isolates

Quantitative trait loci (QTL) with small effects, which are pervasive in quantitative phenotypic variation, are difficult to detect in genome-wide association studies (GWAS). To improve their detection, we propose to use a local score approach that accounts for the surrounding signal due to linkage disequilibrium, by accumulating association signals from contiguous single markers. Simulations revealed that, in a GWAS context with high marker density, the local score approach outperforms single SNP p-value-based tests for detecting minor QTL (heritability of 5-10%) and is competitive with regard to alternative methods, which also aggregate p-values. Using more than five million SNPs, this approach was applied to identify loci involved in Quantitative Disease Resistance (QDR) to different isolates of the plant root rot pathogen Aphanomyces euteiches, from a GWAS performed on a collection of 174 accessions of the model legume Medicago truncatula. We refined the position of a previously reported major locus, underlying MYB/NB-ARC/tyrosine kinase candidate genes conferring resistance to two closely related A. euteiches isolates belonging to pea pathotype I. We also discovered a diversity of minor resistance QTL, not detected using p-value-based tests, some of which being putatively shared in response to pea (pathotype I and III) and/or alfalfa (race 1 and 2) isolates. Candidate genes underlying these QTL suggest pathogen effector recognition and plant proteasome as key functions associated with M. truncatula resistance to A. euteiches. GWAS on any organism can benefit from the local score approach to uncover many weak-effect QTL.

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

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

Functional Genomics and Flowering Time in Medicago truncatula: An Overview

Flowering time is an important trait that influences adaptation and yield in many crop legumes. Both the inherent earliness of flowering and the degree to which it is responsive to environmental factors determine both the eco-geographic range across which crops can be successfully grown and the seasonal cycles most suitable for production. This chapter will provide a brief review of studies investigating the genetic control of flowering time in Medicago truncatula.

Soil environment is a key driver of adaptation in Medicago truncatula: new insights from landscape genomics

Spatial differences in environmental selective pressures interact with the genomes of organisms, ultimately leading to local adaptation. Landscape genomics is an emergent research area that uncovers genome-environment associations, thus allowing researchers to identify candidate loci for adaptation to specific environmental variables. In the present study, we used latent factor mixed models (LFMMs) and Moran spectral outlier detection/randomization (MSOD-MSR) to identify candidate loci for adaptation to 10 environmental variables (climatic, soil and atmospheric) among 43 515 single nucleotide polymorphisms (SNPs) from 202 accessions of the model legume Medicago truncatula. Soil variables were associated with a large number of candidate loci identified through both LFMMs and MSOD-MSR. Genes tagged by candidate loci associated with drought and salinity are involved in the response to biotic and abiotic stresses, while those tagged by candidates associated with soil nitrogen and atmospheric nitrogen, participate in the legume-rhizobia symbiosis. Candidate SNPs identified through both LFMMs and MSOD-MSR explained up to 56% of variance in flowering traits. Our findings highlight the importance of soil in driving adaptation in the system and elucidate the basis of evolutionary potential of M. truncatula to respond to global climate change and anthropogenic disruption of the nitrogen cycle.

Medicago truncatula: Genetic and Genomic Resources

Medicago truncatula was chosen by the legume community, along with Lotus japonicus, as a model plant to study legume biology. Since then, numerous resources and tools have been developed for M. truncatula. These include, for example, its genome sequence, core ecotype collections, transformation/regeneration methods, extensive mutant collections, and a gene expression atlas. This review aims to describe the different genetic and genomic tools and resources currently available for M. truncatula. We also describe how these resources were generated and provide all the information necessary to access these resources and use them from a practical point of view. © 2017 by John Wiley & Sons, Inc.

Plant hormone signaling in flowering: An epigenetic point of view

Reproduction is one of the most important phases in an organism's lifecycle. In the case of angiosperm plants, flowering provides the major developmental transition from the vegetative to the reproductive stage, and requires genetic and epigenetic reprogramming to ensure the success of seed production. Flowering is regulated by a complex network of genes that integrate multiple environmental cues and endogenous signals so that flowering occurs at the right time; hormone regulation, signaling and homeostasis are very important in this process. Working alone or in combination, hormones are able to promote flowering by epigenetic regulation. Some plant hormones, such as gibberellins, jasmonic acid, abscisic acid and auxins, have important effects on chromatin compaction mediated by DNA methylation and histone posttranslational modifications, which hints at the role that epigenetic regulation may play in flowering through hormone action. miRNAs have been viewed as acting independently from DNA methylation and histone modification, ignoring their potential to interact with hormone signaling - including the signaling of auxins, gibberellins, ethylene, jasmonic acid, salicylic acid and others - to regulate flowering. Therefore, in this review we examine new findings about interactions between epigenetic mechanisms and key players in hormone signaling to coordinate flowering.