Multiplex CRISPR/Cas9-mediated mutagenesis of alfalfa FLOWERING LOCUS Ta1 (MsFTa1) leads to delayed flowering time with improved forage biomass yield and quality

Alfalfa (Medicago sativa L.) is a perennial flowering plant in the legume family that is widely cultivated as a forage crop for its high yield, forage quality and related agricultural and economic benefits. Alfalfa is a photoperiod sensitive long-day (LD) plant that can accomplish its vegetative and reproductive phases in a short period of time. However, rapid flowering can compromise forage biomass yield and quality. Here, we attempted to delay flowering in alfalfa using multiplex CRISPR/Cas9-mediated mutagenesis of FLOWERING LOCUS Ta1 (MsFTa1), a key floral integrator and activator gene. Four guide RNAs (gRNAs) were designed and clustered in a polycistronic tRNA–gRNA system and introduced into alfalfa by Agrobacterium-mediated transformation. Ninety-six putative mutant lines were identified by gene sequencing and characterized for delayed flowering time and related desirable agronomic traits. Phenotype assessment of flowering time under LD conditions identified 22 independent mutant lines with delayed flowering compared to the control. Six independent Msfta1 lines containing mutations in all four copies of MsFTa1 accumulated significantly higher forage biomass yield, with increases of up to 78% in fresh weight and 76% in dry weight compared to controls. Depending on the harvesting schemes, many of these lines also had reduced lignin, acid detergent fibre (ADF) and neutral detergent fibre (NDF) content and significantly higher crude protein (CP) and mineral contents compared to control plants, especially in the stems. These CRISPR/Cas9-edited Msfta1 mutants could be introduced in alfalfa breeding programmes to generate elite transgene-free alfalfa cultivars with improved forage biomass yield and quality.     

Identification of loci controlling forage yield and nutritive value in diploid alfalfa using GBS-GWAS

Key message

We attempted to identify genomic regions controlling forage yield and nutritive value in alfalfa. Several candidate genes and associated genetic markers were identified that could potentially be useful for alfalfa breeding to more efficiently develop improved cultivars.

Abstract

Alfalfa is one of the most widely cultivated forage legumes worldwide and improving alfalfa forage yield and nutritive value is a major global breeding goal. Genotyping-by-sequencing (GBS) provides cost-effective molecular marker genotyping for genome-wide association studies (GWAS). Using more than 15,000 genome-wide single nucleotide polymorphisms (SNP) identified from GBS, we conducted a GWAS to investigate forage yield and nutritive value-related traits. We have detected a number of associations for all the traits evaluated and a number of associations detected were located on the Medicago truncatula genome. The SNP in a coding region of a cell wall biosynthesis gene was associated with several cell wall-related traits, and we suggest that it may be the causative polymorphism. Two other SNPs residing in meristematic development and early growth genes were found to associate with the total biomass yield. None of the SNPs associated with regrowth after harvest or with spring regrowth were mapped to the M. truncatula genome, possibly reflecting the fact that M. truncatula is an annual species related to alfalfa that typically has limited ability to regrow. The alleles we identify with the major impact on forage yield and nutritive value can be rapidly incorporated into our breeding program.

     

CONSTANS delays Arabidopsis flowering under short days

Long days (LD) promote flowering of Arabidopsis thaliana compared with short days (SD) by activating the photoperiodic pathway. Here we show that growth under very‐SD (3 h) or darkness (on sucrose) also accelerates flowering on a biological scale, indicating that SD actively repress flowering compared with very‐SD. CONSTANS (CO) repressed flowering under SD, and the early flowering of co under SD required FLOWERING LOCUS T (FT). FT was expressed at a basal level in the leaves under SD, but these levels were not enhanced in co. This indicates that the action of CO in A. thaliana is not the mirror image of the action of its homologue in rice. In the apex, CO enhanced the expression of TERMINAL FLOWER 1 (TFL1) around the time when FT expression is important to promote flowering. Under SD, the tfl1 mutation was epistatic to co and in turn ft was epistatic to tfl1. These observations are consistent with the long‐standing but not demonstrated model where CO can inhibit FT induction of flowering by affecting TFL1 expression.     

INDUCER OF CBF EXPRESSION 1 integrates cold signals into FLOWERING LOCUS C‐mediated flowering pathways in Arabidopsis

Plants constantly monitor changes in photoperiod and temperature throughout the year to synchronize flowering with optimal environmental conditions. In the temperate zones, both photoperiod and temperature fluctuate in a somewhat predictable manner through the seasons, although a transient shift to low temperature is also encountered during changing seasons, such as early spring. Although low temperatures are known to delay flowering by inducing the floral repressor FLOWERING LOCUS C (FLC), it is not fully understood how temperature signals are coordinated with photoperiodic signals in the timing of seasonal flowering. Here, we show that the cold signaling activator INDUCER OF CBF EXPRESSION 1 (ICE1), FLC and the floral promoter SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) constitute an elaborate signaling network that integrates cold signals into flowering pathways. The cold‐activated ICE1 directly induces the gene encoding FLC, which represses SOC1 expression, resulting in delayed flowering. In contrast, under floral promotive conditions, SOC1 inhibits the binding of ICE1 to the promoters of the FLC gene, inducing flowering with a reduction of freezing tolerance. These observations indicate that the ICE1‐FLC‐SOC1 signaling network contributes to the fine‐tuning of flowering during changing seasons.     

Gibberellic acid signaling is required for ambient temperature‐mediated induction of flowering in Arabidopsis thaliana

Distinct molecular mechanisms integrate changes in ambient temperature into the genetic pathways that govern flowering time in Arabidopsis thaliana. Temperature‐dependent eviction of the histone variant H2A.Z from nucleosomes has been suggested to facilitate the expression of FT by PIF4 at elevated ambient temperatures. Here we show that, in addition to PIF4, PIF3 and PIF5, but not PIF1 and PIF6, can promote flowering when expressed specifically in phloem companion cells (PCC), where they can induce FT and its close paralog, TSF. However, despite their strong potential to promote flowering, genetic analyses suggest that the PIF genes seem to have only a minor role in adjusting flowering in response to photoperiod or high ambient temperature. In addition, loss of PIF function only partially suppressed the early flowering phenotype and FT expression of the arp6 mutant, which is defective in H2A.Z deposition. In contrast, the chemical inhibition of gibberellic acid (GA) biosynthesis resulted in a strong attenuation of early flowering and FT expression in arp6. Furthermore, GA was able to induce flowering at low temperature (15°C) independently of FT, TSF, and the PIF genes, probably directly at the shoot apical meristem. Together, our results suggest that the timing of the floral transition in response to ambient temperature is more complex than previously thought and that GA signaling might play a crucial role in this process.     

Effects of plant age on performance of the tropical perennial fodder grass,Cenchrus ciliaris L. subjected to elevated ultraviolet‐B radiation

• Abiotic stress, notably high ultraviolet‐B (eUV‐B), limit growth and productivity of many crop plants, but information on response of forage grasses to eUV‐B radiation is rather limited.
• The present study was therefore conducted to increase our understanding of differential age‐related responses on growth, metabolism and fodder quality of Cenchrus ciliaris‐3108 (Buffel grass) to elevated UV‐B (eUV‐B: 7.2 kJ·m^?2·day^?1).
• Plant growth at both growth stages was notably reduced in response to eUV‐B, except for the number of nodes and tillers at vegetative and reproductive stages. At anthesis, tillering increased due to the perennial habit of this plant, but leaf senescence reduced the number of leaves per tiller. Unlike ambient UV‐B, eUV‐B at the vegetative stage resulted in diversion of photosynthate for the formation of secondary metabolites (tannins and phenolics), providing dual protection from photooxidative damage and from herbivory.
• The forage biomass as well as quality showed a marked decline under eUV‐B and relative nutritive value was reduced at both growth stages.

     

Protein accumulation and rumen stability of wheat γ‐gliadin fusion proteins in tobacco and alfalfa

The nutritional value of various crops can be improved by engineering plants to produce high levels of proteins. For example, because methionine deficiency limits the protein quality of Medicago Sativa (alfalfa) forage, producing alfalfa plants that accumulate high levels of a methionine‐rich protein could increase the nutritional value of that crop. We used three strategies in designing methionine‐rich recombinant proteins that could accumulate to high levels in plants and thereby serve as candidates for improving the protein quality of alfalfa forage. In tobacco, two fusion proteins, γ‐gliadin‐δ‐zein and γ‐δ‐zein, as well as δ‐zein co‐expressed with β‐zein, all formed protein bodies. However, the γ‐gliadin‐δ‐zein fusion protein accumulated to the highest level, representing up to 1.5% of total soluble protein (TSP) in one transformant. In alfalfa, γ‐gliadin‐δ‐zein accumulated to 0.2% of TSP, and in an in vitro rumen digestion assay, γ‐gliadin‐δ‐zein was more resistant to microbial degradation than Rubisco. Additionally, although it did not form protein bodies, a γ‐gliadin‐GFP fusion protein accumulated to much higher levels, 7% of TSP, than a recombinant protein comprised of an ER localization signal fused to GFP in tobacco. Based on our results, we conclude that γ‐gliadin‐δ‐zein is a potential candidate protein to use for enhancing methionine levels in plants and for improving rumen stability of forage protein. γ‐gliadin fusion proteins may provide a general platform for increasing the accumulation of recombinant proteins in transgenic plants.     

Bombus terrestris in a mass‐flowering pollinator‐dependent crop: A mutualistic relationship?

Bumblebees (Bombus spp.) rely on an abundant and diverse selection of floral resources to meet their nutritional requirements. In farmed landscapes, mass‐flowering crops can provide an important forage resource for bumblebees, with increased visitation from bumblebees into mass‐flowering crops having an additional benefit to growers who require pollination services. This study explores the mutualistic relationship between Bombus terrestris L. (buff‐tailed bumblebee), a common species in European farmland, and the mass‐flowering crop courgette (Cucurbita pepo L.) to see how effective B. terrestris is at pollinating courgette and in return how courgette may affect B. terrestris colony dynamics. By combining empirical data on nectar and pollen availability with model simulations using the novel bumblebee model Bumble‐BEEHAVE, we were able to quantify and simulate for the first time, the importance of courgette as a mass‐flowering forage resource for bumblebees. Courgette provides vast quantities of nectar to ensure a high visitation rate, which combined with abundant pollen grains, enables B. terrestris to have a high pollination potential. While B. terrestris showed a strong fidelity to courgette flowers for nectar, courgette pollen was not found in any pollen loads from returning foragers. Nonetheless, model simulations showed that early season courgette (nectar) increased the number of hibernating queens, colonies, and adult workers in the modeled landscapes. Synthesis and applications. Courgette has the potential to improve bumblebee population dynamics; however, the lack of evidence of the bees collecting courgette pollen in this study suggests that bees can only benefit from this transient nectar source if alternative floral resources, particularly pollen, are also available to fulfill bees’ nutritional requirements in space and time. Therefore, providing additional forage resources could simultaneously improve pollination services and bumblebee populations.     

Nitrogen as a key regulator of flowering in Fagus crenata: understanding the physiological mechanism of masting by gene expression analysis

The role of resource availability in determining the incidence of masting has been widely studied, but how floral transition and initiation are regulated by the resource level is unclear. We tested the hypothesis that floral transition is stimulated by high resource availabiltiy in Fagus crenata based on a new technique, the expression analyses of flowering genes. We isolated F. crenata orthologues of FLOWERING LOCUS T, LEAFY and APETALA1, and confirmed their functions using transgenic Arabidopsis thaliana. We monitored the gene expression levels for 5 years and detected a cycle of on and off years, which was correlated with fluctuations of the shoot‐nitrogen concentration. Nitrogen fertilisation resulted in the significantly higher expression of flowering genes than the control, where all of the fertilised trees flowered, whereas the control did not. Our findings identified nitrogen as a key regulator of mast flowering, thereby providing new empirical evidence to support the resource budget model.     

Natural variation and CRISPR/Cas9‐mediated mutation in GmPRR37 affect photoperiodic flowering and contribute to regional adaptation of soybean

Flowering time is a critical determinant of the geographic distribution and regional adaptability of soybean (Glycine max) and is strongly regulated by photoperiod and temperature. In this study, quantitative trait locus (QTL) mapping and subsequent candidate gene analysis revealed that GmPRR37, encoding a pseudo‐response regulator protein, is responsible for the major QTL qFT12‐2, which was identified from a population of 308 recombinant inbred lines (RILs) derived from a cross between a very late‐flowering soybean cultivar, ‘Zigongdongdou (ZGDD)’, and an extremely early‐flowering cultivar, ‘Heihe27 (HH27)’, in multiple environments. Comparative analysis of parental sequencing data confirmed that HH27 contains a non‐sense mutation that causes the loss of the CCT domain in the GmPRR37 protein. CRISPR/Cas9‐induced Gmprr37‐ZGDD mutants in soybean exhibited early flowering under natural long‐day (NLD) conditions. Overexpression of GmPRR37 significantly delayed the flowering of transgenic soybean plants compared with wild‐type under long photoperiod conditions. In addition, both the knockout and overexpression of GmPRR37 in soybean showed no significant phenotypic alterations in flowering time under short‐day (SD) conditions. Furthermore, GmPRR37 down‐regulated the expression of the flowering‐promoting FT homologues GmFT2a and GmFT5a, and up‐regulated flowering‐inhibiting FT homologue GmFT1a expression under long‐day (LD) conditions. We analysed haplotypes of GmPRR37 among 180 cultivars collected across China and found natural Gmprr37 mutants flower earlier and enable soybean to be cultivated at higher latitudes. This study demonstrates that GmPRR37 controls soybean photoperiodic flowering and provides opportunities to breed optimized cultivars with adaptation to specific regions and farming systems.     

Modeling daily flowering probabilities: expected impact of climate change on Japanese cherry phenology

Understanding the drivers of phenological events is vital for forecasting species’ responses to climate change. We developed flexible Bayesian survival regression models to assess a 29‐year, individual‐level time series of flowering phenology from four taxa of Japanese cherry trees (Prunus spachiana, Prunus × yedoensis, Prunus jamasakura, and Prunus lannesiana), from the Tama Forest Cherry Preservation Garden in Hachioji, Japan. Our modeling framework used time‐varying (chill and heat units) and time‐invariant (slope, aspect, and elevation) factors. We found limited differences among taxa in sensitivity to chill, but earlier flowering taxa, such as P. spachiana, were more sensitive to heat than later flowering taxa, such as P. lannesiana. Using an ensemble of three downscaled regional climate models under the A1B emissions scenario, we projected shifts in flowering timing by 2100. Projections suggest that each taxa will flower about 30 days earlier on average by 2100 with 2–6 days greater uncertainty around the species mean flowering date. Dramatic shifts in the flowering times of cherry trees may have implications for economically important cultural festivals in Japan and East Asia. The survival models used here provide a mechanistic modeling approach and are broadly applicable to any time‐to‐event phenological data, such as plant leafing, bird arrival time, and insect emergence. The ability to explicitly quantify uncertainty, examine phenological responses on a fine time scale, and incorporate conditions leading up to an event may provide future insight into phenologically driven changes in carbon balance and ecological mismatches of plants and pollinators in natural populations and horticultural crops.     

Seed traits are pleiotropically regulated by the flowering time gene PERPETUAL FLOWERING 1 (PEP1) in the perennial Arabis alpina

The life cycles of plants are characterized by two major life history transitions—germination and the initiation of flowering—the timing of which are important determinants of fitness. Unlike annuals, which make the transition from the vegetative to reproductive phase only once, perennials iterate reproduction in successive years. The floral repressor PERPETUAL FLOWERING 1 (PEP1), an ortholog of FLOWERING LOCUS C, in the alpine perennial Arabis alpina ensures the continuation of vegetative growth after flowering and thereby restricts the duration of the flowering episode. We performed greenhouse and garden experiments to compare flowering phenology, fecundity and seed traits between A. alpina accessions that have a functional PEP1 allele and flower seasonally and pep1 mutants and accessions that carry lesions in PEP1 and flower perpetually. In the garden, perpetual genotypes flower asynchronously and show higher winter mortality than seasonal ones. PEP1 also pleiotropically regulates seed dormancy and longevity in a way that is functionally divergent from FLC. Seeds from perpetual genotypes have shallow dormancy and reduced longevity regardless of whether they after‐ripened in plants grown in the greenhouse or in the experimental garden. These results suggest that perpetual genotypes have higher mortality during winter but compensate by showing higher seedling establishment. Differences in seed traits between seasonal and perpetual genotypes are also coupled with differences in hormone sensitivity and expression of genes involved in hormonal pathways. Our study highlights the existence of pleiotropic regulation of seed traits by hub developmental regulators such as PEP1, suggesting that seed and flowering traits in perennial plants might be optimized in a coordinated fashion.     

The trxG family histone methyltransferase SET DOMAIN GROUP 26 promotes flowering via a distinctive genetic pathway

Histone methylation is a major component in numerous processes such as determination of flowering time, which is fine‐tuned by multiple genetic pathways that integrate both endogenous and environmental signals. Previous studies identified SET DOMAIN GROUP 26 (SDG26) as a histone methyltransferase involved in the activation of flowering, as loss of function of SDG26 caused a late‐flowering phenotype in Arabidopsis thaliana. However, the SDG26 function and the underlying molecular mechanism remain largely unknown. In this study, we undertook a genetic analysis by combining the sdg26 mutant with mutants of other histone methylation enzymes, including the methyltransferase mutants Arabidopsis trithorax1 (atx1), sdg25 and curly leaf (clf), as well as the demethylase double mutant lsd1‐like1 lsd1‐like2 (ldl1 ldl2). We found that the early‐flowering mutants sdg25, atx1 and clf interact antagonistically with the late‐flowering mutant sdg26, whereas the late‐flowering mutant ldl1 ldl2 interacts synergistically with sdg26. Based on microarray analysis, we observed weak overlaps in the genes that were differentially expressed between sdg26 and the other mutants. Our analyses of the chromatin of flowering genes revealed that the SDG26 protein binds at the key flowering integrator SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1/AGAMOUS‐LIKE 20 (SOC1/AGL20), and is required for histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 36 trimethylation (H3K36me3) at this locus. Together, our results indicate that SDG26 promotes flowering time through a distinctive genetic pathway, and that loss of function of SDG26 causes a decrease in H3K4me3 and H3K36me3 at its target gene SOC1, leading to repression of this gene and the late‐flowering phenotype.