Extensive Analysis of GmFTL and GmCOL Expression in Northern Soybean Cultivars in Field Conditions

Editing the nuclear localization signals of E1 and E1Lb enables the production of tropical soybean in temperate growing regions

Soybean is a typical short-day crop, and most commercial soybean cultivars are restricted to a relatively narrow range of latitudes due to photoperiod sensitivity. Photoperiod sensitivity hinders the utilization of soybean germplasms across geographical regions. When grown in temperate regions, tropical soybean responds to prolonged day length by increasing the vegetative growth phase and delaying flowering and maturity, which often pushes the harvest window past the first frost date. In this study, we used CRISPR/LbCas12a to edit a North American subtropical soybean cultivar named 06KG218440 that belongs to maturity group 5.5. By designing one gRNA to edit the nuclear localization signal (NLS) regions of both E1 and E1Lb, we created a series of new germplasms with shortened flowering time and time to maturity and determined their favourable latitudinal zone for cultivation. The novel partial function alleles successfully achieve yield and early maturity trade-offs and exhibit good agronomic traits and high yields in temperate regions. This work offers a straightforward editing strategy to modify subtropical and tropical soybean cultivars for temperate growing regions, a strategy that could be used to enrich genetic diversity in temperate breeding programmes and facilitate the introduction of important crop traits such as disease tolerance or high yield.

Molecular breeding for improvement of photothermal adaptability in soybean

Soybean (Glycine max (L.) Merr.) is a typical short-day and temperate crop that is sensitive to photoperiod and temperature. Responses of soybean to photothermal conditions determine plant growth and development, which affect its architecture, yield formation, and capacity for geographic adaptation. Flowering time, maturity, and other traits associated with photothermal adaptability are controlled by multiple major-effect and minor-effect genes and genotype-by-environment interactions. Genetic studies have identified at least 11 loci (E1-E4, E6-E11, and J) that participate in photoperiodic regulation of flowering time and maturity in soybean. Molecular cloning and characterization of major-effect flowering genes have clarified the photoperiod-dependent flowering pathway, in which the photoreceptor gene phytochrome A, circadian evening complex (EC) components, central flowering repressor E1, and FLOWERING LOCUS T family genes play key roles in regulation of flowering time, maturity, and adaptability to photothermal conditions. Here, we provide an overview of recent progress in genetic and molecular analysis of traits associated with photothermal adaptability, summarizing advances in molecular breeding practices and tools for improving these traits. Furthermore, we discuss methods for breeding soybean varieties with better adaptability to specific ecological regions, with emphasis on a novel strategy, the Potalaization model, which allows breeding of widely adapted soybean varieties through the use of multiple molecular tools in existing elite widely adapted varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01406-z.

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.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Comparative Genomic and Expression Analysis Insight into Evolutionary Characteristics of PEBP Genes in Cultivated Peanuts and Their Roles in Floral Induction

Phosphatidyl ethanolamine-binding proteins (PEBPs) are involved in regulating flowering time and various developmental processes. Functions and expression patterns in cultivated peanuts (Arachis hypogaea L.) remain unknown. In this study, 33 PEBP genes in cultivated peanuts were identified and divided into four subgroups: FT, TFL, MFT and FT-like. Gene structure analysis showed that orthologs from A and B genomes in cultivated peanuts had highly similar structures, but some orthologous genes have subgenomic dominance. Gene collinearity and phylogenetic analysis explain that some PEBP genes play key roles in evolution. Cis-element analysis revealed that PEBP genes are mainly regulated by hormones, light signals and stress-related pathways. Multiple PEPB genes had different expression patterns between early and late-flowering genotypes. Further detection of its response to temperature and photoperiod revealed that PEBPs ArahyM2THPA, ArahyEM6VH3, Arahy4GAQ4U, ArahyIZ8FG5, ArahyG6F3P2, ArahyLUT2QN, ArahyDYRS20 and ArahyBBG51B were the key genes controlling the flowering response to different flowering time genotypes, photoperiods and temperature. This study laid the foundation for the functional study of the PEBP gene in cultivated peanuts and the adaptation of peanuts to different environments.

Natural allelic variation of GmST05 controlling seed size and quality in soybean

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.

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.

A near-complete genome sequence of mungbean (Vigna radiata L.) provides key insights into the modern breeding program

Mungbean (Vigna radiata L.), a fast-growing legume species, is an important source of carbohydrates and proteins in developing countries of Asia. Here, we constructed a near-complete genome sequence of mungbean with a scaffold N50 value of 5.2 Mb and only a 0.4% gap, with a total scaffold size of 475 Mb. We identified several misassembled pseudomolecules (Chr03, Chr04, Chr05, and Chr08) in the previous draft assembly; Chr03, Chr04, and Chr08 were assembled into one chromosome, and Chr05 was broken into two chromosomes in the improved reference genome assembly, thus providing more accurate linkage information to breeders. Additionally, using an ultra-high-resolution linkage map constructed based on resequencing data, we identified several quantitative trait loci (QTLs) and the underlying candidate genes affecting synchronous pod maturity (SPM). Mungbean homologs of two soybean ([Glycine max (L.) Merr.] flowering genes, E3 (phytochrome A) and J (early flowering 3), were identified as candidate genes for the QTLs, and the candidate genes for plant height, node number, and SPM showed critical nucleotide substitutions between the reference cultivar and other genotypes (landraces and wild accessions). Based on the analysis of genetic diversity among 276 accessions collected from 23 countries, we identified 36 selective sweep regions and observed that the overall genetic diversity of cultivars decreased to 30% of that in wild accessions postdomestication. The near-complete genome sequence of mungbean represents an important resource for genome-assisted improvement in the mungbean breeding program.

Fine-Tuning Florigen Increases Field Yield Through Improving Photosynthesis in Soybean

Crop yield has been maintaining its attraction for researchers because of the demand of global population growth. Mutation of flowering activators, such as florigen, increases plant biomass at the expense of later flowering, which prevents crop maturity in the field. As a result, it is difficult to apply flowering activators in agriculture production. Here, we developed a strategy to utilize florigen to significantly improve soybean yield in the field. Through the screening of transgenic lines of RNAi-silenced florigen homologs in soybean (Glycine-max-Flowering Locus T Like, GmFTL), we identified a line, GmFTL-RNAi#1, with minor changes in both GmFTL expression and flowering time but with notable increase in soybean yield. As expected, GmFTL-RNAi#1 matured normally in the field and exhibited markedly high yield over multiple locations and years, indicating that it is possible to reach a trade-off between flowering time and high yield through the fine-tuning expression of flowering activators. Further studies uncovered an unknown mechanism by which GmFTL negatively regulates photosynthesis, a substantial source of crop yield, demonstrating a novel function of florigen. Thus, because of the highly conserved functions of florigen in plants and the classical RNAi approach, the findings provide a promising strategy to harness early flowering genes to improve crop yield.

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.

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.

Positional Cloning of the Flowering Time QTL qFT12-1 Reveals the Link Between the Clock Related PRR Homolog With Photoperiodic Response in Soybeans

Flowering time and maturity are important agronomic traits for soybean cultivars to adapt to different latitudes and achieve maximal yield. Genetic studies on genes and quantitative trait loci (QTL) that control flowering time and maturity are extensive. In particular, the molecular bases of E1-E4, E6, E9, E10, and J have been deciphered. For a better understanding of regulation of flowering time gene networks, we need to understand if more molecular factors carrying different biological functions are also involved in the regulation of flowering time in soybeans. We developed a population derived from a cross between a landrace Jilincailihua (male) and a Chinese cultivar Chongnong16 (female). Both parents carry the same genotypes of E1e2E3HaE4 at E1, E2, E3, and E4 loci. Nighty-six individuals of the F₂ population were genotyped with Illumina SoySNP8k iSelect BeadChip. A total of 2,407 polymorphic single nucleotide polymorphism (SNP) markers were used to construct a genetic linkage map. One major QTL, qFT12-1, was mapped to an approximately 567-kB region on chromosome 12. Genotyping and phenotyping of recombinant plant whose recombination events were occurring within the QTL region allowed us to narrow down the QTL region to 56.4 kB, in which four genes were annotated. Allelism and association analysis indicated Glyma.12G073900, a PRR7 homolog, is the strongest candidate gene for qFT12-1. The findings of this study disclosed the possible involvement of circadian clock gene in flowering time regulation of soybeans.

Molecular mechanisms of flowering under long days and stem growth habit in soybean

Precise timing of flowering is critical to crop adaptation and productivity in a given environment. A number of classical E genes controlling flowering time and maturity have been identified in soybean [Glycine max (L.) Merr.]. The public availability of the soybean genome sequence has accelerated the identification of orthologues of Arabidopsis flowering genes and their functional analysis, and has allowed notable progress towards understanding the molecular mechanisms of flowering in soybean. Great progress has been made particularly in identifying genes and modules that inhibit flowering in long-day conditions, because a reduced or absent inhibition of flowering by long daylengths is an essential trait for soybean, a short-day (SD) plant, to expand its adaptability toward higher latitude environments. In contrast, the molecular mechanism of floral induction by SDs remains elusive in soybean. Here we present an update on recent work on molecular mechanisms of flowering under long days and of stem growth habit, outlining the progress in the identification of genes responsible, the interplay between photoperiod and age-dependent miRNA-mediated modules, and the conservation and divergence in photoperiodic flowering and stem growth habit in soybean relative to other legumes, Arabidopsis, and rice (Oryza sativa L.).

A soybean quantitative trait locus that promotes flowering under long days is identified as FT5a, a FLOWERING LOCUS T ortholog

FLOWERING LOCUS T (FT) is an important floral integrator whose functions are conserved across plant species. In soybean, two orthologs, FT2a and FT5a, play a major role in initiating flowering. Their expression in response to different photoperiods is controlled by allelic combinations at the maturity loci E1 to E4, generating variation in flowering time among cultivars. We determined the molecular basis of a quantitative trait locus (QTL) for flowering time in linkage group J (Chromosome 16). Fine-mapping delimited the QTL to a genomic region of 107kb that harbors FT5a We detected 15 DNA polymorphisms between parents with the early-flowering (ef) and late-flowering (lf) alleles in the promoter region, an intron, and the 3' untranslated region of FT5a, although the FT5a coding regions were identical. Transcript abundance of FT5a was higher in near-isogenic lines for ef than in those for lf, suggesting that different transcriptional activities or mRNA stability caused the flowering time difference. Single-nucleotide polymorphism (SNP) calling from re-sequencing data for 439 cultivated and wild soybean accessions indicated that ef is a rare haplotype that is distinct from common haplotypes including lf The ef allele at FT5a may play an adaptive role at latitudes where early flowering is desirable.