Evolutionary trajectories of duplicated FT homologues and their roles in soybean domestication

Genome editing, an opportunity to revive soybean cultivation in Europe

Soybean (Glycine max Merr.) is the world's most important oilseed crop and its ability to fix atmospheric nitrogen makes it a cornerstone of sustainable agriculture. Despite its importance, Europe relies heavily on imports, leading to environmental and economic vulnerabilities. To address these challenges, the European Union has implemented policies to boost local soybean production, emphasizing sustainable practices and reduced dependency on imports. However, conventional breeding methods are time-consuming and may not keep pace with the rapid environmental and consumer habit changes. Genome-editing technologies, such as CRISPR-Cas, offer precise and efficient tools for developing soybean varieties tailored to European conditions. These technologies can enhance traits related to precocity, stress responses, yield and quality that are essential for adapting to climate change and promoting ecological sustainability. This review explores the integration of genome editing (GE) in soybean breeding, highlighting its potential in advancing the agroecological transition in Europe. By having a clear regulation and enhancing breeding efforts, GE can significantly contribute to developing resilient and sustainable soybean varieties, fostering a competitive and environmentally friendly European agriculture.

Identification of genes associated with sex expression and sex determination in hemp (Cannabis sativa L.)

Dioecy in flowering plants has evolved independently many times, and thus the genetic mechanisms underlying sex determination are diverse. In hemp (Cannabis sativa), sex is controlled by a pair of sex chromosomes (XX for females and XY for males). In an attempt to understand the molecular mechanism responsible for sex expression in hemp plants, we carried out RNA sequencing of male and female plants at different developmental stages. Using a pipeline involving differential gene expression analysis and weighted gene co-expression network analysis, we identified genes important for male and female flower development. We also demonstrate that sex-biased expression is already established at very early vegetative stages, before the onset of reproductive development, and identify several genes encoding transcription factors of the REM, bZIP, and MADS families as candidate sex-determination genes in hemp. Our findings demonstrate that the gene regulatory networks governing male and female development in hemp diverge at a very early stage, leading to profound morphological differences between male and female hemp plants.

Evolution of FLOWERING LOCUS T-like genes in angiosperms: a core Lamiales-specific diversification

Plant life history is determined by two transitions, germination and flowering time, in which the phosphatidylethanolamine-binding proteins (PEBPs) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) play key regulatory roles. Compared with the highly conserved TFL1-like genes, FT-like genes vary significantly in copy numbers in gymnosperms, and monocots within the angiosperms, while sporadic duplications can be observed in eudicots. Here, via a systematic analysis of the PEBPs in angiosperms with a special focus on 12 representative species featuring high-quality genomes in the order Lamiales, we identified a successive lineage-specific but systematic expansion of FT-like genes in the families of core Lamiales. The first expansion event generated FT1-like genes mainly via a core Lamiales-specific whole-genome duplication (cL-WGD), while a likely random duplication produced the FT2-like genes in the lineages containing Scrophulariaceae and the rest of the core Lamiales. Both FT1- and FT2-like genes were further amplified tandemly in some families. These expanded FT-like genes featured highly diverged expression patterns and structural variation, indicating functional diversification. Intriguingly, some core Lamiales contained the relict MOTHER OF FT AND TFL1 like 2 (MFT2) that probably expanded in the common ancestor of angiosperms. Our data showcase the highly dynamic lineage-specific expansion of the FT-like genes, and thus provide important and fresh evolutionary insights into the gene regulatory network underpinning flowering time diversity in Lamiales and, more generally, in angiosperms.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

GmNF-YC4 delays soybean flowering and maturation by directly repressing GmFT2a and GmFT5a expression

Flowering time and growth period are key agronomic traits which directly affect soybean (Glycine max (L.) Merr.) adaptation to diverse latitudes and farming systems. The FLOWERING LOCUS T (FT) homologs GmFT2a and GmFT5a integrate multiple flowering regulation pathways and significantly advance flowering and maturity in soybean. Pinpointing the genes responsible for regulating GmFT2a and GmFT5a will improve our understanding of the molecular mechanisms governing growth period in soybean. In this study, we identified the Nuclear Factor Y-C (NFY-C) protein GmNF-YC4 as a novel flowering suppressor in soybean under long-day (LD) conditions. GmNF-YC4 delays flowering and maturation by directly repressing the expression of GmFT2a and GmFT5a. In addition, we found that a strong selective sweep event occurred in the chromosomal region harboring the GmNF-YC4 gene during soybean domestication. The GmNF-YC4Hap3 allele was mainly found in wild soybean (Glycine soja Siebold & Zucc.) and has been eliminated from G. max landraces and improved cultivars, which predominantly contain the GmNF-YC4Hap1 allele. Furthermore, the Gmnf-yc4 mutants displayed notably accelerated flowering and maturation under LD conditions. These alleles may prove to be valuable genetic resources for enhancing soybean adaptability to higher latitudes.

Light recovery after maize harvesting promotes soybean flowering in a maize-soybean relay strip intercropping system

Moving from sole cropping to intercropping is a transformative change in agriculture, contributing to yield. Soybeans adapt to light conditions in intercropping by adjusting the onset of reproduction and the inflorescence architecture to optimize reproductive success. Maize-soybean strip intercropping (MS), maize-soybean relay strip intercropping (IS), and sole soybean (SS) systems are typical soybean planting systems with significant differences in light environments during growth periods. To elucidate the effect of changes in the light environment on soybean flowering processes and provide a theoretical basis for selecting suitable varieties in various planting systems to improve yields, field experiments combining planting systems (IS, MS, and SS) and soybean varieties (GQ8, GX7, ND25, and NN996) were conducted in 2021 and 2022. Results showed that growth recovery in the IS resulted in a balance in the expression of TERMINAL FLOWER 1 (TFL1) and FLOWERING LOCUS T (FT) in the meristematic tissues of soybeans, which promoted the formation of new branches or flowers. IS prolonged the flowering time (2-7 days) and increased the number of forming flowers compared with SS (93.0 and 169%) and MS (67.3 and 103.3%) at the later soybean flowering stage. The higher carbon and nitrogen content in the middle and bottom canopies of soybean contributed to decreased flower abscission by 26.7 and 30.2%, respectively, compared with SS. Canopy light environment recovery promoted branch and flower formation and transformation of flowers into pods with lower flower-pod abscission, which contributed to elevating soybean yields in late-maturing and multibranching varieties (ND25) in IS.

Circadian and photoperiodic regulation of the vegetative to reproductive transition in plants

As sessile organisms, plants must respond constantly to ever-changing environments to complete their life cycle; this includes the transition from vegetative growth to reproductive development. This process is mediated by photoperiodic response to sensing the length of night or day through circadian regulation of light-signaling molecules, such as phytochromes, to measure the length of night to initiate flowering. Flowering time is the most important trait to optimize crop performance in adaptive regions. In this review, we focus on interplays between circadian and light signaling pathways that allow plants to optimize timing for flowering and seed production in Arabidopsis, rice, soybean, and cotton. Many crops are polyploids and domesticated under natural selection and breeding. In response to adaptation and polyploidization, circadian and flowering pathway genes are epigenetically reprogrammed. Understanding the genetic and epigenetic bases for photoperiodic flowering will help improve crop yield and resilience in response to climate change.

Soybean WRINKLED1 protein GmWRI1a promotes flowering under long-day conditions via regulating expressions of flowering-related genes

Flowering time is an important agronomic character that influences the adaptability and yield of soybean [Glycine max (L.) Merrill]. WRINKLED 1 (WRI1) plays an important regulatory role in plant growth and development. In this study, we found that the expression of GmWIR1a could be induced by long days. Compared with the wild type, transgenic soybean overexpressing GmWRI1a showed earlier flowering and maturity under long days but no significant changes under short days. Overexpression of GmWRI1a led to up-regulated expression of genes involved in the regulation of flowering time. The GmWRI1a protein was able to directly bind to the promoter regions of GmAP1, GmFUL1a, GmFUL2 and up-regulated their expression. GmCOL3 was identified by yeast one-hybrid library screening using the GmWRI1a promoter as bait. GmCOL3 was revealed to be a nucleus-localized protein that represses the transcription of GmWRI1a. Expression of GmCOL3 was induced by short days. Taken together, the results show that overexpression of GmWRI1a promotes flowering under long days by promoting the transcriptional activity of flowering-related genes in soybean, and that GmCOL3 binds to the GmWRI1a promoter and directly down-regulates its transcription. This discovery reveals a new function for GmWRI1a, which regulates flowering and maturity in soybean.

An overview of floral regulatory genes in annual and perennial plants

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

An Improved Genome-Wide Association Procedure Explores Gene-Allele Constitutions and Evolutionary Drives of Growth Period Traits in the Global Soybean Germplasm Population

In soybeans (Glycine max (L.) Merr.), their growth periods, DSF (days of sowing-to-flowering), and DFM (days of flowering-to-maturity) are determined by their required accumulative day-length (ADL) and active temperature (AAT). A sample of 354 soybean varieties from five world eco-regions was tested in four seasons in Nanjing, China. The ADL and AAT of DSF and DFM were calculated from daily day-lengths and temperatures provided by the Nanjing Meteorological Bureau. The improved restricted two-stage multi-locus genome-wide association study using gene-allele sequences as markers (coded GASM-RTM-GWAS) was performed. (i) For DSF and its related ADL_DSF and AAT_DSF, 130-141 genes with 384-406 alleles were explored, and for DFM and its related ADL_DFM and AAT_DFM, 124-135 genes with 362-384 alleles were explored, in a total of six gene-allele systems. DSF shared more ADL and AAT contributions than DFM. (ii) Comparisons between the eco-region gene-allele submatrices indicated that the genetic adaptation from the origin to the geographic sub-regions was characterized by allele emergence (mutation), while genetic expansion from primary maturity group (MG)-sets to early/late MG-sets featured allele exclusion (selection) without allele emergence in addition to inheritance (migration). (iii) Optimal crosses with transgressive segregations in both directions were predicted and recommended for breeding purposes, indicating that allele recombination in soybean is an important evolutionary drive. (iv) Genes of the six traits were mostly trait-specific involved in four categories of 10 groups of biological functions. GASM-RTM-GWAS showed potential in detecting directly causal genes with their alleles, identifying differential trait evolutionary drives, predicting recombination breeding potentials, and revealing population gene networks.

FLOWERING LOCUS T indel variants confer vernalization-independent and photoperiod-insensitive flowering of yellow lupin (Lupinus luteus L.)

Ongoing climate change has considerably reduced the seasonal window for crop vernalization, concurrently expanding cultivation area into northern latitudes with long-day photoperiod. To address these changes, cool season legume breeders need to understand molecular control of vernalization and photoperiod. A key floral transition gene integrating signals from these pathways is the Flowering locus T (FT). Here, a recently domesticated grain legume, yellow lupin (Lupinus luteus L.), was explored for potential involvement of FT homologues in abolition of vernalization and photoperiod requirements. Two FTa (LlutFTa1a and LlutFTa1b) and FTc (LlutFTc1 and LlutFTc2) homologues were identified and sequenced for two contrasting parents of a reference recombinant inbred line (RIL) population, an early-flowering cultivar Wodjil and a late-flowering wild-type P28213. Large deletions were detected in the 5' promoter regions of three FT homologues. Quantitative trait loci were identified for flowering time and vernalization response in the RIL population and in a diverse panel of wild and domesticated accessions. A 2227 bp deletion found in the LlutFTc1 promoter was linked with early phenology and vernalization independence, whereas LlutFTa1a and LlutFTc2 indels with photoperiod responsiveness. Comparative mapping highlighted convergence of FTc1 indel evolution in two Old World lupin species, addressing both artificial selection during domestication and natural adaptation to short season environmental conditions. We concluded that rapid flowering in yellow lupin is associated with the de-repression of the LlutFTc1 homologue from the juvenile phase, putatively due to the elimination of all binding sites in the promoter region for the AGAMOUS-like 15 transcription factor.

Putting the pea in photoPEAriod

This article comments on:

Williams O, Vander Schoor JK, Butler JB, Ridge S, Sussmilch FC, Hecht VFG, Weller JL. 2022. The genetic architecture of flowering time changes in pea from wild to crop. Journal of Experimental Botany 73,3978–3990.

Comparative analysis of buds transcriptome and identification of two florigen gene AkFTs in Amorphophallus konjac

Leaves and flowers of Amorphophallus konjac do not develop simultaneously thus unique features can be elucidated through study of flowering transformation in A. konjac. In this study, transcriptome libraries of A. konjac leaf buds (LB) and flower buds (FB) were constructed followed by high-throughput sequencing. A total of 68,906 unigenes with an average length of 920 bp were obtained after library assembly. Out of these genes, 24,622 unigenes had annotation information. A total of 6859 differentially expressed genes (DEGs) were identified through differential expression analysis using LB as control. Notably, 2415 DEGs were upregulated whereas 4444 DEGs were downregulated in the two transcriptomes. Go and KEGG analysis showed that the DEGs belonged to 44 functional categories and were implicated in 98 metabolic pathways and 38 DEGs involved in plant hormone signal transduction. Several genes were mined that may be involved in A. konjac flower bud differentiation and flower organ development. Eight DEGs were selected for verification of RNA-seq results using qRT-PCR analysis. Two FLOWERING LOCUS T (FT) genes named AkFT1 and AkFT2 were identified though homologous analysis may be the florigen gene implicated in modulation of A. konjac flowering. These genes were significantly upregulated in flower buds compared with the expression levels on leaf buds. Overexpression of AkFT genes though heterologous expression in Arabidopsis showed that the transgenics flowered at a very early stage relative to wild type plants. These findings indicate that AkFT1 and AkFT2 function as regulation genes in A. konjac flowering development and the two genes may present similar functions during flowering transition.

Genetic Diversity and Selection Footprints in the Genome of Brazilian Soybean Cultivars

Although Brazil is currently the largest soybean producer in the world, only a small number of studies have analyzed the genetic diversity of Brazilian soybean. These studies have shown the existence of a narrow genetic base. The objectives of this work were to analyze the population structure and genetic diversity, and to identify selection signatures in the genome of soybean germplasms from different companies in Brazil. A panel consisting of 343 soybean lines from Brazil, North America, and Asia was genotyped using genotyping by sequencing (GBS). Population structure was assessed by Bayesian and multivariate approaches. Genetic diversity was analyzed using metrics such as the fixation index, nucleotide diversity, genetic dissimilarity, and linkage disequilibrium. The software BayeScan was used to detect selection signatures between Brazilian and Asian accessions as well as among Brazilian germplasms. Region of origin, company of origin, and relative maturity group (RMG) all had a significant influence on population structure. Varieties belonging to the same company and especially to the same RMG exhibited a high level of genetic similarity. This result was exacerbated among early maturing accessions. Brazilian soybean showed significantly lower genetic diversity when compared to Asian accessions. This was expected, because the crop's region of origin is its main genetic diversity reserve. We identified 7 genomic regions under selection between the Brazilian and Asian accessions, and 27 among Brazilian varieties developed by different companies. Associated with these genomic regions, we found 96 quantitative trait loci (QTLs) for important soybean breeding traits such as flowering, maturity, plant architecture, productivity components, pathogen resistance, and seed composition. Some of the QTLs associated with the markers under selection have genes of great importance to soybean's regional adaptation. The results reported herein allowed to expand the knowledge about the organization of the genetic variability of the Brazilian soybean germplasm. Furthermore, it was possible to identify genomic regions under selection possibly associated with the adaptation of soybean to Brazilian environments.

Sequencing the USDA core soybean collection reveals gene loss during domestication and breeding

The gene content of plants varies between individuals of the same species due to gene presence/absence variation, and selection can alter the frequency of specific genes in a population. Selection during domestication and breeding will modify the genomic landscape, though the nature of these modifications is only understood for specific genes or on a more general level (e.g., by a loss of genetic diversity). Here we have assembled and analyzed a soybean (Glycine spp.) pangenome representing more than 1,000 soybean accessions derived from the USDA Soybean Germplasm Collection, including both wild and cultivated lineages, to assess genomewide changes in gene and allele frequency during domestication and breeding. We identified 3,765 genes that are absent from the Lee reference genome assembly and assessed the presence/absence of all genes across this population. In addition to a loss of genetic diversity, we found a significant reduction in the average number of protein-coding genes per individual during domestication and subsequent breeding, though with some genes and allelic variants increasing in frequency associated with selection for agronomic traits. This analysis provides a genomic perspective of domestication and breeding in this important oilseed crop.

Functional Redundancy of FLOWERING LOCUS T 3b in Soybean Flowering Time Regulation

Photoperiodic flowering is an important agronomic trait that determines adaptability and yield in soybean and is strongly influenced by FLOWERING LOCUS T (FT) genes. Due to the presence of multiple FT homologs in the genome, their functions in soybean are not fully understood. Here, we show that GmFT3b exhibits functional redundancy in regulating soybean photoperiodic flowering. Bioinformatic analysis revealed that GmFT3b is a typical floral inducer FT homolog and that the protein is localized to the nucleus. Moreover, GmFT3b expression was induced by photoperiod and circadian rhythm and was more responsive to long-day (LD) conditions. We generated a homozygous ft3b knockout and three GmFT3b-overexpressing soybean lines for evaluation under different photoperiods. There were no significant differences in flowering time between the wild-type, the GmFT3b overexpressors, and the ft3b knockouts under natural long-day, short-day, or LD conditions. Although the downstream flowering-related genes GmFUL1 (a, b), GmAP1d, and GmLFY1 were slightly down-regulated in ft3b plants, the floral inducers GmFT5a and GmFT5b were highly expressed, indicating potential compensation for the loss of GmFT3b. We suggest that GmFT3b acts redundantly in flowering time regulation and may be compensated by other FT homologs in soybean.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Towards understanding the biological foundations of perenniality

Perennial life cycles enable plants to have remarkably long lifespans, as exemplified by trees that can live for thousands of years. For this, they require sophisticated regulatory networks that sense environmental changes and initiate adaptive responses in their growth patterns. Recent research has gradually elucidated fundamental mechanisms underlying the perennial life cycle. Intriguingly, several conserved components of the floral transition pathway in annuals such as Arabidopsis thaliana also participate in these regulatory mechanisms underpinning perenniality. Here, we provide an overview of perennials' physiological features and summarise their recently discovered molecular foundations. We also highlight the importance of deepening our understanding of perenniality in the development of perennial grain crops, which are promising elements of future sustainable agriculture.

Molecular Genetic Understanding of Photoperiodic Regulation of Flowering Time in Arabidopsis and Soybean

The developmental switch from a vegetative phase to reproduction (flowering) is essential for reproduction success in flowering plants, and the timing of the floral transition is regulated by various environmental factors, among which seasonal day-length changes play a critical role to induce flowering at a season favorable for seed production. The photoperiod pathways are well known to regulate flowering time in diverse plants. Here, we summarize recent progresses on molecular mechanisms underlying the photoperiod control of flowering in the long-day plant Arabidopsis as well as the short-day plant soybean; furthermore, the conservation and diversification of photoperiodic regulation of flowering in these two species are discussed.

Development of a marker panel for genotyping of domestic soybean cultivars for genes controlling the duration of vegetation and response to photoperiod

Soybean, Glycine max L., is one of the most important agricultural crops grown in a wide range of latitude. In this regard, in soybean breeding, it is necessary to pay attention to the set of genes that control the transition to the f lowering stage, which will make it possible to adapt genotypes to local growing conditions as accurately as possible. The possibilities of soybean breeding for this trait have now signif icantly expanded due to identif ication of the main genes (E1–E4, GmFT2a, GmFT5a) that control the processes of f lowering and maturation in soybean, depending on the day length. The aim of this work was to develop a panel of markers for these genes, which could be used for a rapid and eff icient genotyping of domestic soybean cultivars and selection of plant material based on sensitivity to photoperiod and the duration of vegetation. Combinations of 10 primers, both previously developed and our own, were tested to identify different alleles of the E1–E4, GmFT2a, and GmFT5a genes using 10 soybean cultivars from different maturity groups. As a result, 5 combinations of dominant and recessive alleles for the E1–E4 genes were identif ied: (1) e1-nl(e1-as)/ e2-ns/e3-tr(e3-fs)/e4; (2) e1-as/e2-ns/e3-tr/E4; (3) e1-as/e2-ns/E3-Ha/e4; (4) E1/e2-ns/e3-tr/E4; (5) e1-nl/e2-ns/E3-Ha/E4. The studied cultivars contained the most common alleles of the GmFT2a and GmFT5a genes, with the exception of the ‘Cassidi’ cultivar having a rare dominant allele GmFT5a-H4. The degree of earliness of cultivars positively correlated with the number of recessive genes E1–E4, which is consistent with the data of foreign authors on different sets of cultivars from Japan and North China. Thus, the developed panel of markers can be successfully used in the selection of soybean for earliness and sensitivity to photoperiod.

Transcriptomal dissection of soybean circadian rhythmicity in two geographically, phenotypically and genetically distinct cultivars

BACKGROUND

In soybean, some circadian clock genes have been identified as loci for maturity traits. However, the effects of these genes on soybean circadian rhythmicity and their impacts on maturity are unclear.

RESULTS

We used two geographically, phenotypically and genetically distinct cultivars, conventional juvenile Zhonghuang 24 (with functional J/GmELF3a, a homolog of the circadian clock indispensable component EARLY FLOWERING 3) and long juvenile Huaxia 3 (with dysfunctional j/Gmelf3a) to dissect the soybean circadian clock with time-series transcriptomal RNA-Seq analysis of unifoliate leaves on a day scale. The results showed that several known circadian clock components, including RVE1, GI, LUX and TOC1, phase differently in soybean than in Arabidopsis, demonstrating that the soybean circadian clock is obviously different from the canonical model in Arabidopsis. In contrast to the observation that ELF3 dysfunction results in clock arrhythmia in Arabidopsis, the circadian clock is conserved in soybean regardless of the functional status of J/GmELF3a. Soybean exhibits a circadian rhythmicity in both gene expression and alternative splicing. Genes can be grouped into six clusters, C1-C6, with different expression profiles. Many more genes are grouped into the night clusters (C4-C6) than in the day cluster (C2), showing that night is essential for gene expression and regulation. Moreover, soybean chromosomes are activated with a circadian rhythmicity, indicating that high-order chromosome structure might impact circadian rhythmicity. Interestingly, night time points were clustered in one group, while day time points were separated into two groups, morning and afternoon, demonstrating that morning and afternoon are representative of different environments for soybean growth and development. However, no genes were consistently differentially expressed over different time-points, indicating that it is necessary to perform a circadian rhythmicity analysis to more thoroughly dissect the function of a gene. Moreover, the analysis of the circadian rhythmicity of the GmFT family showed that GmELF3a might phase- and amplitude-modulate the GmFT family to regulate the juvenility and maturity traits of soybean.

CONCLUSIONS

These results and the resultant RNA-seq data should be helpful in understanding the soybean circadian clock and elucidating the connection between the circadian clock and soybean maturity.

Molecular mechanisms for the photoperiodic regulation of flowering in soybean

Photoperiodic flowering is one of the most important factors affecting regional adaptation and yield in soybean (Glycine max). Plant adaptation to long-day conditions at higher latitudes requires early flowering and a reduction or loss of photoperiod sensitivity; adaptation to short-day conditions at lower latitudes involves delayed flowering, which prolongs vegetative growth for maximum yield potential. Due to the influence of numerous major loci and quantitative trait loci (QTLs), soybean has broad adaptability across latitudes. Forward genetic approaches have uncovered the molecular basis for several of these major maturity genes and QTLs. Moreover, the molecular characterization of orthologs of Arabidopsis thaliana flowering genes has enriched our understanding of the photoperiodic flowering pathway in soybean. Building on early insights into the importance of the photoreceptor phytochrome A, several circadian clock components have been integrated into the genetic network controlling flowering in soybean: E1, a repressor of FLOWERING LOCUS T orthologs, plays a central role in this network. Here, we provide an overview of recent progress in elucidating photoperiodic flowering in soybean, how it contributes to our fundamental understanding of flowering time control, and how this information could be used for molecular design and breeding of high-yielding soybean cultivars.

Natural variation and artificial selection of photoperiodic flowering genes and their applications in crop adaptation

Flowering links vegetative growth and reproductive growth and involves the coordination of local environmental cues and plant genetic information. Appropriate timing of floral initiation and maturation in both wild and cultivated plants is important to their fitness and productivity in a given growth environment. The domestication of plants into crops, and later crop expansion and improvement, has often involved selection for early flowering. In this review, we analyze the basic rules for photoperiodic adaptation in several economically important and/or well-researched crop species. The ancestors of rice (Oryza sativa), maize (Zea mays), soybean (Glycine max), and tomato (Solanum lycopersicum) are short-day plants whose photosensitivity was reduced or lost during domestication and expansion to high-latitude areas. Wheat (Triticum aestivum) and barley (Hordeum vulgare) are long-day crops whose photosensitivity is influenced by both latitude and vernalization type. Here, we summarize recent studies about where these crops were domesticated, how they adapted to photoperiodic conditions as their growing area expanded from domestication locations to modern cultivating regions, and how allelic variants of photoperiodic flowering genes were selected during this process. A deeper understanding of photoperiodic flowering in each crop will enable better molecular design and breeding of high-yielding cultivars suited to particular local environments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00039-0.

The Synchronized Efforts to Decipher the Molecular Basis for Soybean Maturity Loci E1, E2, and E3 That Regulate Flowering and Maturity

The general concept of photoperiodism, i.e., the photoperiodic induction of flowering, was established by Garner and Allard (1920). The genetic factor controlling flowering time, maturity, or photoperiodic responses was observed in soybean soon after the discovery of the photoperiodism. E1, E2, and E3 were named in 1971 and, thereafter, genetically characterized. At the centennial celebration of the discovery of photoperiodism in soybean, we recount our endeavors to successfully decipher the molecular bases for the major maturity loci E1, E2, and E3 in soybean. Through systematic efforts, we successfully cloned the E3 gene in 2009, the E2 gene in 2011, and the E1 gene in 2012. Recently, successful identification of several circadian-related genes such as PRR3a, LUX, and J has enriched the known major E1-FTs pathway. Further research progresses on the identification of new flowering and maturity-related genes as well as coordinated regulation between flowering genes will enable us to understand profoundly flowering gene network and determinants of latitudinal adaptation in soybean.

Molecular characterization of a soybean FT homologue, GmFT7

Soybean (Glycine max) is a vital oilseed legume crop that provides protein and oil for humans and feedstock for animals. Flowering is a prerequisite for seed production. Floral transition, from vegetative to reproductive stage, in a plant, is regulated by environmental (light, temperature) and endogenous factors. In Arabidopsis, Flowering Locus T (FT) protein is shown to be a mobile signal that moves from leaf to shoot apical meristem to induce flowering. However, FTs role in soybean is not fully resolved due to the presence of multiple (ten) homologs in the genome. Two of the ten FT homologs (GmFT2a and GmFT5a) have a role in the floral transition while GmFT1a and GmFT4 suppress soybean flowering. Recent deep sequencing data revealed that six FT homologs are expressed in shoot apical meristem and leaves during floral transition. One FT homolog, GmFT7 showed strong expression during soybean floral transition. Though bioinformatic analyses revealed that GmFT7 had high similarity with GmFT2a, ectopic GmFT7 expression in Arabidopsis could not promote flowering or rescue the late-flowering phenotype of Arabidopsis ft-10 mutant.

Cotranscriptional and Posttranscriptional Features of the Transcriptome in Soybean Shoot Apex and Leaf

Transcription is the first step of central dogma, in which the genetic information stored in DNA is copied into RNA. In addition to mature RNA sequencing (RNA-seq), high-throughput nascent RNA assays have been established and applied to provide detailed transcriptional information. Here, we present the profiling of nascent RNA from trifoliate leaves and shoot apices of soybean. In combination with nascent RNA (chromatin-bound RNA, CB RNA) and RNA-seq, we found that introns were largely spliced cotranscriptionally. Although alternative splicing (AS) was mainly determined at nascent RNA biogenesis, differential AS between the leaf and shoot apex at the mature RNA level did not correlate well with cotranscriptional differential AS. Overall, RNA abundance was moderately correlated between nascent RNA and mature RNA within each tissue, but the fold changes between the leaf and shoot apex were highly correlated. Thousands of novel transcripts (mainly non-coding RNA) were detected by CB RNA-seq, including the overlap of natural antisense RNA with two important genes controlling soybean reproductive development, FT2a and Dt1. Taken together, we demonstrated the adoption of CB RNA-seq in soybean, which may shed light on gene expression regulation of important agronomic traits in leguminous crops.

CaFT-LIKE is a flowering promoter in pepper and functions as florigen in tomato

We identified a pepper late-flowering mutant that is disrupted in the sequence of CaFT-LIKE, the ortholog of tomato SINGLE FLOWER TRUSS (SFT). Heterologous expression in tomato indicated that CaFT-LIKE has a conserved function as a flowering promoter and can rescue the wild-type phenotype of the tomato sft mutant. CaFT-LIKE confers a graft-transmissible signal for flowering initiation in tomato, implicating its function as a florigen. To test the relationship between CaFT-LIKE and FASCICULATE (FA), the ortholog of tomato SELF PRUNING (SP), we constructed the double mutant Caft-like fa. The phenotype of Caft-like fa resembled that of Caft-like, indicating epistasis of Caft-like over fa in controlling flowering time and sympodial shoot structure. To examine the association between the expression pattern of flowering genes and natural variation in flowering time, the expression levels of CaFT-LIKE and the flowering repressor CaAP2 were determined in a panel of early-flowering cultivars and late-flowering landraces and wild accessions. Strong positive and negative correlations between flowering time and expression levels of CaAP2 and CaFT-LIKE, respectively, were observed, indicating that high-expression alleles of CaFT-LIKE and low-expression alleles of CaAP2 were selected for early flowering during pepper domestication and breeding.

Evolutionary Insights into the Nature of Plant Domestication

Domestication is a co-evolutionary process that occurs when wild plants are brought into cultivation by humans, leading to origin of new species and/or differentiated populations that are critical for human survival. Darwin used domesticated species as early models for evolution, highlighting their variation and the key role of selection in species differentiation. Over the last two decades, a growing synthesis of plant genetics, genomics, and archaeobotany has led to challenges to old orthodoxies and the advent of fresh perspectives on how crop domestication and diversification proceed. I discuss four new insights into plant domestication - that in general domestication is a protracted process, that unconscious (natural) selection plays a prominent role, that interspecific hybridization may be an important mechanism for crop species diversification and range expansion, and that similar genes across multiple species underlies parallel/convergent phenotypic evolution between domesticated taxa. Insights into the evolutionary origin and diversification of crop species can help us in developing new varieties (and possibly even new species) to deal with current and future environmental challenges in a sustainable manner.

A Domestication-Associated Gene GmPRR3b Regulates the Circadian Clock and Flowering Time in Soybean

Improved soybean cultivars have been adapted to grow at a wide range of latitudes, enabling expansion of cultivation worldwide. However, the genetic basis of this broad adaptation is still not clear. Here, we report the identification of GmPRR3b as a major flowering time regulatory gene that has been selected during domestication and genetic improvement for geographic expansion. Through a genome-wide association study of a diverse soybean landrace panel consisting of 279 accessions, we identified 16 candidate quantitative loci associated with flowering time and maturity time. The strongest signal resides in the known flowering gene E2, verifying the effectiveness of our approach. We detected strong signals associated with both flowering and maturity time in a genomic region containing GmPRR3b. Haplotype analysis revealed that GmPRR3bH6 is the major form of GmPRR3b that has been utilized during recent breeding of modern cultivars. mRNA profiling analysis showed that GmPRR3bH6 displays rhythmic and photoperiod-dependent expression and is preferentially induced under long-day conditions. Overexpression of GmPRR3bH6 increased main stem node number and yield, while knockout of GmPRR3bH6 using CRISPR/Cas9 technology delayed growth and the floral transition. GmPRR3bH6 appears to act as a transcriptional repressor of multiple predicted circadian clock genes, including GmCCA1a, which directly upregulates J/GmELF3a to modulate flowering time. The causal SNP (Chr12:5520945) likely endows GmPRR3bH6 a moderate but appropriate level of activity, leading to early flowering and vigorous growth traits preferentially selected during broad adaptation of landraces and improvement of cultivars.

Analysis of the quinoa genome reveals conservation and divergence of the flowering pathways

Quinoa (Chenopodium quinoa Willd.) is a grain crop grown in the Andes renowned as a highly nutritious plant exhibiting tolerance to abiotic stress such as drought, cold and high salinity. Quinoa grows across a range of latitudes corresponding to differing day lengths, suggesting regional adaptations of flowering regulation. Improved understanding and subsequent modification of the flowering process, including flowering time, ensuring high yields, is one of the key factors behind expansion of cultivation zones and goals of the crop improvement programs worldwide. However, our understanding of the molecular basis of flower initiation and development in quinoa is limited. Here, we use a computational approach to perform genome-wide identification and analysis of 611 orthologues of the Arabidopsis thaliana flowering genes. Conservation of the genes belonging to the photoperiod, gibberellin and autonomous pathways was observed, while orthologues of the key genes found in the vernalisation pathway (FRI, FLC) were absent from the quinoa genome. Our analysis indicated that on average each Arabidopsis flowering gene has two orthologous copies in quinoa. Several genes including orthologues of MIF1, FT and TSF were identified as homologue-rich genes in quinoa. We also identified 459 quinoa-specific genes uniquely expressed in the flower and/or meristem, with no known orthologues in other species. The genes identified provide a resource and framework for further studies of flowering in quinoa and related species. It will serve as valuable resource for plant biologists, crop physiologists and breeders to facilitate further research and establishment of modern breeding programs for quinoa.

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.

B-box Proteins in Arachis duranensis: Genome-Wide Characterization and Expression Profiles Analysis

B-box (BBX) proteins are important factors involved in plant growth and developmental regulation, and they have been identified in many species. However, information on the characteristics and transcription patterns of BBX genes in wild peanut are limited. In this study, we identified and characterized 24 BBX genes from a wild peanut, Arachis duranensis. Many characteristics were analyzed, including chromosomal locations, phylogenetic relationships, and gene structures. Arachis duranensis B-box (AdBBX) proteins were grouped into five classes based on the diversity of their conserved domains: I (3 genes), II (4 genes), III (4 genes), IV (9 genes), and V (4 genes). Fifteen distinct motifs were found in the 24 AdBBX proteins. Duplication analysis revealed the presence of two interchromosomal duplicated gene pairs, from group II and IV. In addition, 95 kinds of cis-acting elements were found in the genes’ promoter regions, 53 of which received putative functional predictions. The numbers and types of cis-acting elements varied among different AdBBX promoters, and, as a result, AdBBX genes exhibited distinct expression patterns in different tissues. Transcriptional profiling combined with synteny analysis suggests that AdBBX8 may be a key factor involved in flowering time regulation. Our study will provide essential information for further functional investigation of AdBBX genes.

Functional divergence between soybean FLOWERING LOCUS T orthologues FT2a and FT5a in post-flowering stem growth

Genes in the FLOWERING LOCUS T (FT) family integrate external and internal signals to control various aspects of plant development. In soybean (Glycine max), FT2a and FT5a play a major role in floral induction, but their roles in post-flowering reproductive development remain undetermined. Ectopic overexpression analyses revealed that FT2a and FT5a similarly induced flowering, but FT5a was markedly more effective than FT2a for the post-flowering termination of stem growth. The down-regulation of Dt1, a soybean orthologue of Arabidopsis TERMINAL FLOWER1, in shoot apices in early growing stages of FT5a-overexpressing plants was concomitant with highly up-regulated expression of APETALA1 orthologues. The Dt2 gene, a repressor of Dt1, was up-regulated similarly by the overexpression of FT2a and FT5a, suggesting that it was not involved in the control of stem termination by FT5a. In addition to the previously reported interaction with FDL19, a homologue of the Arabidopsis bZIP protein FD, both FT2a and FT5a interacted with FDL12, but only FT5a interacted with FDL06. Our results suggest that FT2a and FT5a have different functions in the control of post-flowering stem growth. A specific interaction of FT5a with FDL06 may play a key role in determining post-flowering stem growth in soybean.

Natural variations of FT family genes in soybean varieties covering a wide range of maturity groups

BACKGROUND

Flowering time and maturity are among the most important adaptive traits in soybean (Glycine max (L.) Merill). Flowering Locus T (FT) family genes function as key flowering integrators, with flowering-promoting members GmFT2a/GmFT5a and flowering-inhibiting members GmFT4/GmFT1a antagonistically regulating vegetative and reproductive growth. However, to date, the relations between natural variations of FT family genes and the diversity of flowering time and maturity in soybean are not clear. Therefore, we conducted this study to discover natural variations in FT family genes in association with flowering time and maturity.

RESULTS

Ten FT family genes, GmFT1a, GmFT1b, GmFT2a, GmFT2b, GmFT3a, GmFT3b, GmFT4, GmFT5a, GmFT5b and GmFT6, were cloned and sequenced in the 127 varieties evenly covering all 14 known maturity groups (MG0000-MGX). They were diversified at the genome sequence polymorphism level. GmFT3b and GmFT5b might have experienced breeding selection in the process of soybean domestication and breeding. Haplotype analysis showed that a total of 17 haplotypes had correlative relationships with flowering time and maturity among the 10 FT genes, namely, 1a-H3, 1b-H1, 1b-H6, 1b-H7, 2a-H1, 2a-H3, 2a-H4, 2a-H9, 2b-H3, 2b-H4, 2b-H6, 2b-H7, 3b-H4, 5a-H1, 5a-H2, 5a-H4 and 5b-H1. Based on the association analysis, 38 polymorphic sites had a significant association with flowering time at the level of p < 0.01.

CONCLUSIONS

Some natural variations exist within the 10 FT family genes, which might be involved in soybean adaptation to different environments and have an influence on diverse flowering time and maturity. This study will facilitate the understanding of the roles of FTs in flowering and maturity.

Parallel origins of photoperiod adaptation following dual domestications of common bean

Common bean (Phaseolus vulgaris L.) is an important grain legume domesticated independently in Mexico and Andean South America approximately 8000 years ago. Wild forms are obligate short-day plants, and relaxation of photoperiod sensitivity was important for expansion to higher latitudes and subsequent global spread. To better understand the nature and origin of this key adaptation, we examined its genetic control in progeny of a wide cross between a wild accession and a photoperiod-insensitive cultivar. We found that photoperiod sensitivity is under oligogenic control, and confirm a major effect of the Ppd locus on chromosome 1. The red/far-red photoreceptor gene PHYTOCHROME A3 (PHYA3) was identified as a strong positional candidate for Ppd, and sequencing revealed distinct deleterious PHYA3 mutations in photoperiod-insensitive Andean and Mesoamerican accessions. These results reveal the independent origins of photoperiod insensitivity within the two major common bean gene pools and demonstrate the conserved importance of PHYA genes in photoperiod adaptation of short-day legume species.

Transcriptome-Enabled Network Inference Revealed the GmCOL1 Feed-Forward Loop and Its Roles in Photoperiodic Flowering of Soybean

Photoperiodic flowering, a plant response to seasonal photoperiod changes in the control of reproductive transition, is an important agronomic trait that has been a central target of crop domestication and modern breeding programs. However, our understanding about the molecular mechanisms of photoperiodic flowering regulation in crop species is lagging behind. To better understand the regulatory gene networks controlling photoperiodic flowering of soybeans, we elucidated global gene expression patterns under different photoperiod regimes using the near isogenic lines (NILs) of maturity loci (E loci). Transcriptome signatures identified the unique roles of the E loci in photoperiodic flowering and a set of genes controlled by these loci. To elucidate the regulatory gene networks underlying photoperiodic flowering regulation, we developed the network inference algorithmic package CausNet that integrates sparse linear regression and Granger causality heuristics, with Gaussian approximation of bootstrapping to provide reliability scores for predicted regulatory interactions. Using the transcriptome data, CausNet inferred regulatory interactions among soybean flowering genes. Published reports in the literature provided empirical verification for several of CausNet's inferred regulatory interactions. We further confirmed the inferred regulatory roles of the flowering suppressors GmCOL1a and GmCOL1b using GmCOL1 RNAi transgenic soybean plants. Combinations of the alleles of GmCOL1 and the major maturity locus E1 demonstrated positive interaction between these genes, leading to enhanced suppression of flowering transition. Our work provides novel insights and testable hypotheses in the complex molecular mechanisms of photoperiodic flowering control in soybean and lays a framework for de novo prediction of biological networks controlling important agronomic traits in crops.

Abundant Small Genetic Alterations after Upland Cotton Domestication

Domestication has long been recognized as the most direct and effective way to intentionally influence morphological and physiological phenotypes in plants and animals. Consequently, understanding how small genetic alterations contribute to domestication is of considerable importance. In this study, we resequenced the genome of the wild upland cotton variety Gossypium hirsutum var. yucatanense, the putative wild ancestor of cultivated upland cotton, and then compared single nucleotide polymorphism (SNP) and short insertion and deletion (InDel) variations of the genome with the cultivated accession (TM-1) of G. hirsutum. We found approximately 6.6 million SNPs and 0.7 million InDels between the two genomes. Most of the small genetic variations were anchored in the noncoding regions. With regard to potential coding genes, we found 24,035 genes with nonsynonymous SNPs. Interestingly, 2603 genes in domesticated cotton are found that have changed the positions of stop codons or shifted reading frames from that in G. hirsutum var. yucatanense. This suggests that domestication may have been selected for mutations that restored gene function or that wild cotton has undergone a number of gene inactivation events since its divergence from cultivated cotton. The former scenario seems most likely due to the intense selective pressure applied during the domestication process. These results demonstrate that, within a relatively short period of time, the cotton genome has been readjusted through small genetic changes. The current study provides useful clues for seeking interesting genes for cotton improvement.

Natural variation in Brassica FT homeologs influences multiple agronomic traits including flowering time, silique shape, oil profile, stomatal morphology and plant height in B. juncea

Natural structural variants of regulatory proteins causing quantitative phenotypic consequences have not been reported in plants. Herein, we show that 28 natural structural variants of FT homeologs, isolated from 6 species of Brassica, differ with respect to amino-acid substitutions in regions critical for interactions with FD and represent two evolutionarily distinct categories. Analysis of structural models of selected candidates from Brassica juncea (BjuFT_AAMF1) and Brassica napus (BnaFT_CCLF) predicted stronger binding between BjuFT and Arabidopsis thaliana FD. Over-expression of BjuFT and BnaFT in wild type and ft-10 mutant backgrounds of Arabidopsis validated higher potency of BjuFT in triggering floral transition. Analysis of gain-of-function and artificial miRNA mediated silenced lines of B. juncea implicated Brassica FT in multiple agronomic traits beyond flowering, consistent with a pleiotropic effect. Several dependent and independent traits such as lateral branching, silique shape, seed size, oil-profile, stomatal morphology and plant height were found altered in mutant lines. Enhanced FT levels caused early flowering, which in turn was positively correlated to a higher proportion of desirable fatty acids (PUFA). However, higher FT levels also resulted in altered silique shape and reduced seed size, suggesting trait trade-offs. Modulation of FT levels for achieving optimal balance of trait values and parsing pair-wise interactions among a reportoire of regulatory protein homeologs in polyploid genomes are indeed future areas of crop research.

A high-density consensus linkage map of white lupin highlights synteny with narrow-leafed lupin and provides markers tagging key agronomic traits

White lupin (Lupinus albus L.) is a valuable source of seed protein, carbohydrates and oil, but requires genetic improvement to attain its agronomic potential. This study aimed to (i) develop a new high-density consensus linkage map based on new, transcriptome-anchored markers; (ii) map four important agronomic traits, namely, vernalization requirement, seed alkaloid content, and resistance to anthracnose and Phomopsis stem blight; and, (iii) define regions of synteny between the L. albus and narrow-leafed lupin (L. angustifolius L.) genomes. Mapping of white lupin quantitative trait loci (QTLs) revealed polygenic control of vernalization responsiveness and anthracnose resistance, as well as a single locus regulating seed alkaloid content. We found high sequence collinearity between white and narrow-leafed lupin genomes. Interestingly, the white lupin QTLs did not correspond to previously mapped narrow-leafed lupin loci conferring vernalization independence, anthracnose resistance, low alkaloids and Phomopsis stem blight resistance, highlighting different genetic control of these traits. Our suite of allele-sequenced and PCR validated markers tagging these QTLs is immediately applicable for marker-assisted selection in white lupin breeding. The consensus map constitutes a platform for synteny-based gene cloning approaches and can support the forthcoming white lupin genome sequencing efforts.