Genetic redundancy in soybean photoresponses associated with duplication of the phytochrome A gene

The MADS-box transcription factor GmFULc promotes GmZTL4 gene transcription to modulate maturity in soybean

Flowering time and maturity are crucial agronomic traits that affect the regional adaptability of soybean plants. The development of soybean cultivars with early maturity adapted to longer days and colder climates of high latitudes is very important for ensuring normal ripening before frost begins. FUL belongs to the MADS-box transcription factor family and has several duplicated members in soybeans. In this study, we observed that overexpression of GmFULc in the Dongnong 50 cultivar promoted soybean maturity, while GmFULc knockout mutants exhibited late maturity. Chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) revealed that GmFULc could bind to the CArG, bHLH and homeobox motifs. Further investigation revealed that GmFULc could directly bind to the CArG motif in the promoters of the GmZTL3 and GmZTL4 genes. Overexpression of GmZTL4 promoted soybean maturity, whereas the ztl4 mutants exhibited delayed maturity. Moreover, we found that the cis element box 4 motif of the GmZTL4 promoter, a motif of light response elements, played an important role in controlling the growth period. Deletion of this motif shortened the growth period by increasing the expression levels of GmZTL4. Functional investigations revealed that short-day treatment promoted the binding of GmFULc to the promoter of GmZTL4 and inhibited the expression of E1 and E1Lb, ultimately resulting in the promotion of flowering and early maturation. Taken together, these findings suggest a novel photoperiod regulatory pathway in which GmFULc directly activates GmZTL4 to promote earlier maturity in soybean.

Subfunctionalisation and self-repression of duplicated E1 homologues finetunes soybean flowering and adaptation

Soybean is a photoperiod-sensitive staple crop. Its photoperiodic flowering has major consequences for latitudinal adaptation and grain yield. Here, we identify and characterise a flowering locus named Time of flower 4b (Tof4b), which encodes E1-Like b (E1Lb), a homologue of the key soybean floral repressor E1. Tof4b protein physically associates with the promoters of two FLOWERING LOCUS T (FT) genes to repress their transcription and delay flowering to impart soybean adaptation to high latitudes. Three E1 homologues undergo subfunctionalisation and show differential subcellular localisation. Moreover, they all possess self-repression capability and each suppresses the two homologous counterparts. Subfunctionalisation and the transcriptional regulation of E1 genes collectively finetune flowering time and high-latitude adaptation in soybean. We propose a model for the functional fate of the three E1 genes after the soybean whole-genome duplication events, refine the molecular mechanisms underlying high-latitude adaption, and provide a potential molecular-breeding resource.

A critical suppression feedback loop determines soybean photoperiod sensitivity

Photoperiod sensitivity is crucial for soybean flowering, adaptation, and yield. In soybean, photoperiod sensitivity centers around the evening complex (EC) that regulates the transcriptional level of the core transcription factor E1, thereby regulating flowering. However, little is known about the regulation of the activity of EC. Our study identifies how E2/GIGANTEA (GI) and its homologs modulate photoperiod sensitivity through interactions with the EC. During long days, E2 interacts with the blue-light receptor flavin-binding, kelch repeat, F box 1 (FKF1), leading to the degradation of J/ELF3, an EC component. EC also suppresses E2 expression by binding to its promoter. This interplay forms a photoperiod regulatory loop, maintaining sensitivity to photoperiod. Disruption of this loop leads to losing sensitivity, affecting soybean's adaptability and yield. Understanding this loop's dynamics is vital for molecular breeding to reduce soybean's photoperiod sensitivity and develop cultivars with better adaptability and higher yields, potentially leading to the creation of photoperiod-insensitive varieties for broader agricultural applications.

Molecular selection of soybean towards adaptation to Central European agroclimatic conditions

Europe is highly dependent on soybean meal imports and anticipates an increase of domestic plant protein production. Ongoing climate change resulted in northward shift of plant hardiness zones, enabling spring-sowing of freezing-sensitive crops, including soybean. However, it requires efficient reselection of germplasm adapted to relatively short growing season and long-day photoperiod. In the present study, a PCR array has been implemented, targeting early maturity (E1-E4, E7, E9, and E10), pod shattering (qPHD1), and growth determination (Dt1) genes. This array was optimized for routine screening of soybean diversity panel (204 accessions), subjected to the 2018-2020 survey of phenology, morphology, and yield-related traits in a potential cultivation region in Poland. High broad-sense heritability (0.84-0.88) was observed for plant height, thousand grain weight, maturity date, and the first pod height. Significant positive correlations were identified between the number of seeds and pods per plant, between these two traits and seed yield per plant as well as between flowering, maturity, plant height, and first pod height. PCR array genotyping revealed high genetic diversity, yielding 98 allelic combinations. The most remarkable correlations were identified between flowering and E7 or E1, between maturity and E4 or E7 and between plant height and Dt1 or E4. The study demonstrated high applicability of this PCR array for molecular selection of soybean towards adaptation to Central Europe, designating recessive qPHD1 and dominant Dt1, E3, and E4 alleles as major targets to align soybean growth season requirements with the length of the frost-free period, improve plant performance, and increase yield.

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

Environment-induced heritable variations are common in Arabidopsis thaliana

Parental or ancestral environments can induce heritable phenotypic changes, but whether such environment-induced heritable changes are a common phenomenon remains unexplored. Here, we subject 14 genotypes of Arabidopsis thaliana to 10 different environmental treatments and observe phenotypic and genome-wide gene expression changes over four successive generations. We find that all treatments caused heritable phenotypic and gene expression changes, with a substantial proportion stably transmitted over all observed generations. Intriguingly, the susceptibility of a genotype to environmental inductions could be predicted based on the transposon abundance in the genome. Our study thus challenges the classic view that the environment only participates in the selection of heritable variation and suggests that the environment can play a significant role in generating of heritable variations.

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.

Floral-promoting GmFT homologs trigger photoperiodic after-effects: An important mechanism for early-maturing soybean varieties to regulate reproductive development and adapt to high latitudes

Soybean (Glycine max) is a typical short-day plant, but has been widely cultivated in high-latitude long-day (LD) regions because of the development of early-maturing genotypes which are photoperiod-insensitive. However, some early-maturing varieties exhibit significant responses to maturity under different daylengths but not for flowering, depicting an evident photoperiodic after-effect, a poorly understood mechanism. In this study, we investigated the postflowering responses of 11 early-maturing soybean varieties to various preflowering photoperiodic treatments. We confirmed that preflowering SD conditions greatly promoted maturity and other postflowering developmental stages. Soybean homologs of FLOWERING LOCUS T (FT), including GmFT2a, GmFT3a, GmFT3b and GmFT5a, were highly accumulated in leaves under preflowering SD treatment. More importantly, they maintained a high expression level after flowering even under LD conditions. E1 RNAi and GmFT2a overexpression lines showed extremely early maturity regardless of preflowering SD and LD treatments due to constitutively high levels of floral-promoting GmFT homolog expression throughout their life cycle. Collectively, our data indicate that high and stable expression of floral-promoting GmFT homologs play key roles in the maintenance of photoperiodic induction to promote postflowering reproductive development, which confers early-maturing varieties with appropriate vegetative growth and shortened reproductive growth periods for adaptation to high latitudes.

Mechanisms underlying key agronomic traits and implications for molecular breeding in soybean

Soybean (Glycine max [L.] Merr.) is an important crop that provides protein and vegetable oil for human consumption. As soybean is a photoperiod-sensitive crop, its cultivation and yield are limited by the photoperiodic conditions in the field. In contrast to other major crops, soybean has a special plant architecture and a special symbiotic nitrogen fixation system, representing two unique breeding directions. Thus, flowering time, plant architecture, and symbiotic nitrogen fixation are three critical or unique yield-determining factors. This review summarizes the progress made in our understanding of these three critical yield-determining factors in soybean. Meanwhile, we propose potential research directions to increase soybean production, discuss the application of genomics and genomic-assisted breeding, and explore research directions to address future challenges, particularly those posed by global climate changes.

Dissection of the E8 locus in two early maturing Canadian soybean populations

Soybean [Glycine max (L.) Merr.] is a short-day crop for which breeders want to expand the cultivation range to more northern agro-environments by introgressing alleles involved in early reproductive traits. To do so, we investigated quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) regions comprised within the E8 locus, a large undeciphered region (~7.0 Mbp to 44.5 Mbp) associated with early maturity located on chromosome GM04. We used a combination of two mapping algorithms, (i) inclusive composite interval mapping (ICIM) and (ii) genome-wide composite interval mapping (GCIM), to identify major and minor regions in two soybean populations (QS15524F2:F3 and QS15544RIL) having fixed E1, E2, E3, and E4 alleles. Using this approach, we identified three main QTL regions with high logarithm of the odds (LODs), phenotypic variation explained (PVE), and additive effects for maturity and pod-filling within the E8 region: GM04:16,974,874-17,152,230 (E8-r1); GM04:35,168,111-37,664,017 (E8-r2); and GM04:41,808,599-42,376,237 (E8-r3). Using a five-step variant analysis pipeline, we identified Protein far-red elongated hypocotyl 3 (Glyma.04G124300; E8-r1), E1-like-a (Glyma.04G156400; E8-r2), Light-harvesting chlorophyll-protein complex I subunit A4 (Glyma.04G167900; E8-r3), and Cycling dof factor 3 (Glyma.04G168300; E8-r3) as the most promising candidate genes for these regions. A combinatorial eQTL mapping approach identified significant regulatory interactions for 13 expression traits (e-traits), including Glyma.04G050200 (Early flowering 3/E6 locus), with the E8-r3 region. Four other important QTL regions close to or encompassing major flowering genes were also detected on chromosomes GM07, GM08, and GM16. In GM07:5,256,305-5,404,971, a missense polymorphism was detected in the candidate gene Glyma.07G058200 (Protein suppressor of PHYA-105). These findings demonstrate that the locus known as E8 is regulated by at least three distinct genomic regions, all of which comprise major flowering genes.

Genome-wide association study reveals GmFulb as candidate gene for maturity time and reproductive length in soybeans (Glycine max)

The ability of soybean [Glycine max (L.) Merr.] to adapt to different latitudes is attributed to genetic variation in major E genes and quantitative trait loci (QTLs) determining flowering time (R1), maturity (R8), and reproductive length (RL). Fully revealing the genetic basis of R1, R8, and RL in soybeans is necessary to enhance genetic gains in soybean yield improvement. Here, we performed a genome-wide association analysis (GWA) with 31,689 single nucleotide polymorphisms (SNPs) to detect novel loci for R1, R8, and RL using a soybean panel of 329 accessions with the same genotype for three major E genes (e1-as/E2/E3). The studied accessions were grown in nine environments and observed for R1, R8 and RL in all environments. This study identified two stable peaks on Chr 4, simultaneously controlling R8 and RL. In addition, we identified a third peak on Chr 10 controlling R1. Association peaks overlap with previously reported QTLs for R1, R8, and RL. Considering the alternative alleles, significant SNPs caused RL to be two days shorter, R1 two days later and R8 two days earlier, respectively. We identified association peaks acting independently over R1 and R8, suggesting that trait-specific minor effect loci are also involved in controlling R1 and R8. From the 111 genes highly associated with the three peaks detected in this study, we selected six candidate genes as the most likely cause of R1, R8, and RL variation. High correspondence was observed between a modifying variant SNP at position 04:39294836 in GmFulb and an association peak on Chr 4. Further studies using map-based cloning and fine mapping are necessary to elucidate the role of the candidates we identified for soybean maturity and adaptation to different latitudes and to be effectively used in the marker-assisted breeding of cultivars with optimal yield-related traits.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

The FLOWERING LOCUS T 5b positively regulates photoperiodic flowering and improves the geographical adaptation of soybean

Plants can sense the photoperiod to flower at the right time. As a sensitive short-day crop, soybean (Glycine max) flowering varies greatly depending on photoperiods, affecting yields. Adaptive changes in soybeans rely on variable genetic loci such as E1 and FLOWERING LOCUS T orthologs. However, the precise coordination and control of these molecular components remain largely unknown. In this study, we demonstrate that GmFT5b functions as a crucial factor for soybean flowering. Overexpressed or mutated GmFT5b resulted in significantly early or later flowering, altering expression profiles for several downstream flowering-related genes under a long-day photoperiod. GmFT5b interacts with the transcription factor GmFDL15, suggesting transcriptional tuning of flowering time regulatory genes via the GmFT5b/GmFDL15 complex. Notably, GmFT5a partially compensated for GmFT5b function, as ft5a ft5b double mutants exhibited an enhanced late-flowering phenotype. Association mapping revealed that GmFT5b was associated with flowering time, maturity, and geographical distribution of soybean accessions, all associated with the E1 locus. Therefore, GmFT5b is a valuable target for enhancing regional adaptability. Natural variants or multiple mutants in this region can be utilized to generate optimized soybean varieties with precise flowering times.

Introgression of a Danbaekkong high-protein allele across different genetic backgrounds in soybean

Soybean meal is a major component of livestock feed due to its high content and quality of protein. Understanding the genetic control of protein is essential to develop new cultivars with improved meal protein. Previously, a genomic region on chromosome 20 significantly associated with elevated protein content was identified in the cultivar Danbaekkong. The present research aimed to introgress the Danbaekkong high-protein allele into elite lines with different genetic backgrounds by developing and deploying robust DNA markers. A multiparent population consisting of 10 F5-derived populations with a total of 1,115 recombinant inbred lines (RILs) was developed using "Benning HP" as the donor parent of the Danbaekkong high-protein allele. A new functional marker targeting the 321-bp insertion in the gene Glyma.20g085100 was developed and used to track the Danbaekkong high-protein allele across the different populations and enable assessment of its effect and stability. Across all populations, the high-protein allele consistently increased the content, with an increase of 3.3% in seed protein. A total of 103 RILs were selected from the multiparent population for yield testing in five environments to assess the impact of the high-protein allele on yield and to enable the selection of new breeding lines with high protein and high yield. The results indicated that the high-protein allele impacts yield negatively in general; however, it is possible to select high-yielding lines with high protein content. An analysis of inheritance of the Chr 20 high-protein allele in Danbaekkong indicated that it originated from a Glycine soja line (PI 163453) and is the same as other G. soja lines studied. A survey of the distribution of the allele across 79 G. soja accessions and 35 Glycine max ancestors of North American soybean cultivars showed that the high-protein allele is present in all G. soja lines evaluated but not in any of the 35 North American soybean ancestors. These results demonstrate that G. soja accessions are a valuable source of favorable alleles for improvement of protein composition.

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.

PH13 improves soybean shade traits and enhances yield for high-density planting at high latitudes

Shading in combination with extended photoperiods can cause exaggerated stem elongation (ESE) in soybean, leading to lodging and reduced yields when planted at high-density in high-latitude regions. However, the genetic basis of plant height in adaptation to these regions remains unclear. Here, through a genome-wide association study, we identify a plant height regulating gene on chromosome 13 (PH13) encoding a WD40 protein with three main haplotypes in natural populations. We find that an insertion of a Ty1/Copia-like retrotransposon in the haplotype 3 leads to a truncated PH13H3 with reduced interaction with GmCOP1s, resulting in accumulation of STF1/2, and reduced plant height. In addition, PH13H3 allele has been strongly selected for genetic improvement at high latitudes. Deletion of both PH13 and its paralogue PHP can prevent shade-induced ESE and allow high-density planting. This study provides insights into the mechanism of shade-resistance and offers potential solutions for breeding high-yielding soybean cultivar for high-latitude regions.

Identification of genetic loci for early maturity in spring bread wheat using the association analysis and gene dissection

Background

Early maturity in spring bread wheat is highly desirable in the regions where it enables the plants to evade high temperatures and plant pathogens at the end of the growing season.

Methods

To reveal the genetic loci responsible for the maturity time association analysis was carried out based on phenotyping for an 11-year period and high-throughput SNP genotyping of a panel of the varieties contrasting for this trait. The expression of candidate genes was verified using qPCR. The association between the SNP markers and the trait was validated using the biparental F2:3 population.

Results

Our data showed that under long-day conditions, the period from seedling to maturity is mostly influenced by the time from heading to maturity, rather than the heading time. The QTLs associated with the trait were located on 2A, 3B, 4A, 5B, 7A and 7B chromosomes with the 7BL locus being the most significant and promising for its SNPs accelerated the maturity time by about 9 days. Gene dissection in this locus detected a number of candidates, the best being TraesCS7B02G391800 (bZIP9) and TraesCS7B02G412200 (photosystem II reaction center). The two genes are predominantly expressed in the flag leaf while flowering. The effect of the SNPs was verified in F2:3 population and confirmed the association of the 4A, 5B and 7BL loci with the maturity time.

Comparative Analysis of Circadian Transcriptomes Reveals Circadian Characteristics between Arabidopsis and Soybean

The circadian clock, an endogenous timing system, exists in nearly all organisms on Earth. The plant circadian clock has been found to be intricately linked with various essential biological activities. Extensive studies of the plant circadian clock have yielded valuable applications. However, the distinctions of circadian clocks in two important plant species, Arabidopsis thaliana and Glycine max (soybean), remain largely unexplored. This study endeavors to address this gap by conducting a comprehensive comparison of the circadian transcriptome profiles of Arabidopsis and soybean to uncover their distinct circadian characteristics. Utilizing non-linear regression fitting (COS) integrated with weights, we identified circadian rhythmic genes within both organisms. Through an in-depth exploration of circadian parameters, we unveiled notable differences between Arabidopsis and soybean. Furthermore, our analysis of core circadian clock genes shed light on the distinctions in central oscillators between these two species. Additionally, we observed that the homologous genes of Arabidopsis circadian clock genes in soybean exert a significant influence on the regulation of flowering and maturity of soybean. This phenomenon appears to stem from shifts in circadian parameters within soybean genes. These findings highlight contrasting biological activities under circadian regulation in Arabidopsis and soybean. This study not only underscores the distinctive attributes of these species, but also offers valuable insights for further scrutiny into the soybean circadian clock and its potential applications.

Identification of E1-E4 allele combinations and ecological adaptability of soybean varieties from different geographical origins in China

The duration of soybean growth and development is regulated by E1-E4 allele genes, which form the basis for ecological adaptation related to biomass accumulation, flowering and pod formation, maturation, and yield. To elucidate the effects of different combinations of E1-E4 allele genes on soybean ecological adaptation, this study conducted competitive allele-specific PCR (KASP) analysis and photoperiod gene typing on 101 main soybean cultivars from different latitudes in China. The ecological adaptation of these cultivars in Sichuan was also investigated. The results showed that within a certain range (60-95 days), soybean varieties with a genotype combination of E1/e2-ns/E3/E4 exhibited a longer growth period and demonstrated higher biomass and yield, displaying overall better performance. These varieties showed strong ecological adaptation in the Chengdu Plain region and are suitable for introduction in similar low to mid-latitude areas like the Chengdu Plain (30°N~32°N). Conversely, soybean varieties carrying a higher number of recessive alleles of E1-E4 are not suitable for introduction in this region.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

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.

Marker-Assisted Selection for Early Maturing E Loci in Soybean Yielded Prospective Breeding Lines for High Latitudes of Northern Kazakhstan

The photoperiodic sensitivity of soybean (Glycine max L.) is one of the limiting factors affecting plant growth and yield. At higher latitudes, early flowering and maturity with neutral reaction to photoperiods are required for adaptation of soybean plants to long-day conditions. Currently, the production and distribution of new varieties of soybeans adapted to widespread agricultural regions in northern Kazakhstan is in strong demand. Eleven soybean hybrid populations were obtained from crosses between 17 parents with four maturity groups, MG 000, 00, 0, and I. Marker-assisted selection (MAS) was assessed for suitable SSR markers and successfully applied for genes E1, E3, E4, and E7, targeting homozygous genotypes with recessive alleles. The identified and selected genotypes were propagated and tested in the conditions of 53° N latitude in the Kostanay region of northern Kazakhstan. Finally, 20 early maturing F₄ breeding lines were identified and developed with genotypes e1 e3 E4 e7, e1 E3 E4 e7, and e1 E3 e4 e7, all completing their growth period within 92-102 days. These breeding lines were developed by MAS and should provide very prospective superior varieties of soybean for northern Kazakhstan through a strategy that may be very helpful to other countries with high latitudes.

Genome-Wide Association Analysis for Resistance to Coniothyrium glycines Causing Red Leaf Blotch Disease in Soybean

Soybean is a high oil and protein-rich legume with several production constraints. Globally, several fungi, viruses, nematodes, and bacteria cause significant yield losses in soybean. Coniothyrium glycines (CG), the causal pathogen for red leaf blotch disease, is the least researched and causes severe damage to soybean. The identification of resistant soybean genotypes and mapping of genomic regions associated with resistance to CG is critical for developing improved cultivars for sustainable soybean production. This study used single nucleotide polymorphism (SNP) markers generated from a Diversity Arrays Technology (DArT) platform to conduct a genome-wide association (GWAS) analysis of resistance to CG using 279 soybean genotypes grown in three environments. A total of 6395 SNPs was used to perform the GWAS applying a multilocus model Fixed and random model Circulating Probability Unification (FarmCPU) with correction of the population structure and a statistical test p-value threshold of 5%. A total of 19 significant marker-trait associations for resistance to CG were identified on chromosomes 1, 5, 6, 9, 10, 12, 13, 15, 16, 17, 19, and 20. Approximately 113 putative genes associated with significant markers for resistance to red leaf blotch disease were identified across soybean genome. Positional candidate genes associated with significant SNP loci-encoding proteins involved in plant defense responses and that could be associated with soybean defenses against CG infection were identified. The results of this study provide valuable insight for further dissection of the genetic architecture of resistance to CG in soybean. They also highlight SNP variants and genes useful for genomics-informed selection decisions in the breeding process for improving resistance traits in soybean.

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.

Haplotype Analysis of GmSGF14 Gene Family Reveals Its Roles in Photoperiodic Flowering and Regional Adaptation of Soybean

Flowering time and photoperiod sensitivity are fundamental traits that determine soybean adaptation to a given region or a wide range of geographic environments. The General Regulatory Factors (GRFs), also known as 14-3-3 family, are involved in protein-protein interactions in a phosphorylation-dependent manner, thus regulating ubiquitous biological processes, such as photoperiodic flowering, plant immunity and stress response. In this study, 20 soybean GmSGF14 genes were identified and divided into two categories according to phylogenetic relationships and structural characteristics. Real-time quantitative PCR analysis revealed that GmSGF14g, GmSGF14i, GmSGF14j, GmSGF14k, GmSGF14m and GmSGF14s were highly expressed in all tissues compared to other GmSGF14 genes. In addition, we found that the transcript levels of GmSGF14 family genes in leaves varied significantly under different photoperiodic conditions, indicating that their expression responds to photoperiod. To explore the role of GmSGF14 in the regulation of soybean flowering, the geographical distribution of major haplotypes and their association with flowering time in six environments among 207 soybean germplasms were studied. Haplotype analysis confirmed that the GmSGF14m^H4 harboring a frameshift mutation in the 14-3-3 domain was associated with later flowering. Geographical distribution analysis demonstrated that the haplotypes related to early flowering were frequently found in high-latitude regions, while the haplotypes associated with late flowering were mostly distributed in low-latitude regions of China. Taken together, our results reveal that the GmSGF14 family genes play essential roles in photoperiodic flowering and geographical adaptation of soybean, providing theoretical support for further exploring the function of specific genes in this family and varietal improvement for wide adaptability.

The genetic architecture of soybean photothermal adaptation to high latitudes

Soybean is a major plant protein source for both human food and animal feed, but to meet global demands as well as a trend towards regional production, soybean cultivation needs to be expanded to higher latitudes. In this study, we developed a large diversity panel consisting of 1503 early-maturing soybean lines and used genome-wide association mapping to dissect the genetic architecture underlying two crucial adaptation traits, flowering time and maturity. This revealed several known maturity loci, E1, E2, E3, and E4, and the growth habit locus Dt2 as causal candidate loci, and also a novel putative causal locus, GmFRL1, encoding a homolog of the vernalization pathway gene FRIGIDA-like 1. In addition, the scan for quantitative trait locus (QTL)-by-environment interactions identified GmAPETALA1d as a candidate gene for a QTL with environment-dependent reversed allelic effects. The polymorphisms of these candidate genes were identified using whole-genome resequencing data of 338 soybeans, which also revealed a novel E4 variant, e4-par, carried by 11 lines, with nine of them originating from Central Europe. Collectively, our results illustrate how combinations of QTL and their interactions with the environment facilitate the photothermal adaptation of soybean to regions far beyond its center of origin.

Diverse flowering responses subjecting to ambient high temperature in soybean under short-day conditions

Flowering time is one of important agronomic traits determining the crop yield and affected by high temperature. When facing high ambient temperature, plants often initiate early flowering as an adaptive strategy to escape the stress and ensure successful reproduction. However, here we find opposing ways in the short-day crop soybean to respond to different levels of high temperatures, in which flowering accelerates when temperature changes from 25 to 30 °C, but delays when temperature reaches 35 °C under short day. phyA-E1, possibly photoperiodic pathway, is crucial for 35 °C-mediated late flowering, however, does not contribute to promoting flowering at 30 °C. 30 °C-induced up-regulation of FT2a and FT5a leads to early flowering, independent of E1. Therefore, distinct responsive mechanisms are adopted by soybean when facing different levels of high temperatures for successful flowering and reproduction.

Shade avoidance syndrome in soybean and ideotype toward shade tolerance

The shade avoidance syndrome (SAS) in soybean can have destructive effects on yield, as essential carbon resources reserved for yield are diverted to the petiole and stem for exaggerated elongation, resulting in lodging and susceptibility to disease. Despite numerous attempts to reduce the unfavorable impacts of SAS for the development of cultivars suitable for high-density planting or intercropping, the genetic bases and fundamental mechanisms of SAS remain largely unclear. The extensive research conducted in the model plant Arabidopsis provides a framework for understanding the SAS in soybean. Nevertheless, recent investigations suggest that the knowledge obtained from model Arabidopsis may not be applicable to all processes in soybean. Consequently, further efforts are required to identify the genetic regulators of SAS in soybean for molecular breeding of high-yield cultivars suitable for density farming. In this review, we present an overview of the recent developments in SAS studies in soybean and suggest an ideal planting architecture for shade-tolerant soybean intended for high-yield breeding.

Genetic Diversity of Korean Wild Soybean Core Collections and Genome-Wide Association Study for Days to Flowering

The utilization of wild soybean germplasms in breeding programs increases genetic diversity, and they contain the rare alleles of traits of interest. Understanding the genetic diversity of wild germplasms is essential for determining effective strategies that can improve the economic traits of soybeans. Undesirable traits make it challenging to cultivate wild soybeans. This study aimed to construct a core subset of 1467 wild soybean accessions of the total population and analyze their genetic diversity to understand their genetic variations. Genome-wild association studies were conducted to detect the genetic loci underlying the time to flowering for a core subset collection, and they revealed the allelic variation in E genes for predicting maturity using the available resequencing data of wild soybean. Based on principal component and cluster analyses, 408 wild soybean accessions in the core collection covered the total population and were explained by 3 clusters representing the collection regions, namely, Korea, China, and Japan. Most of the wild soybean collections in this study had the E1e2E3 genotype according to association mapping and a resequencing analysis. Korean wild soybean core collections can provide helpful genetic resources to identify new flowering and maturity genes near the E gene loci and genetic materials for developing new cultivars, facilitating the introgression of genes of interest from wild soybean.

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.

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Altered regulation of flowering expands growth ranges and maximizes yields in major crops

Flowering time influences reproductive success in plants and has a significant impact on yield in grain crops. Flowering time is regulated by a variety of environmental factors, with daylength often playing an important role. Crops can be categorized into different types according to their photoperiod requirements for flowering. For instance, long-day crops include wheat (Triticum aestivum), barley (Hordeum vulgare), and pea (Pisum sativum), while short-day crops include rice (Oryza sativa), soybean (Glycine max), and maize (Zea mays). Understanding the molecular regulation of flowering and genotypic variation therein is important for molecular breeding and crop improvement. This paper reviews the regulation of flowering in different crop species with a particular focus on how photoperiod-related genes facilitate adaptation to local environments.

Establishment of a novel experimental system for studying the photoperiodic response of short-day dicots using soybean 'cotyledon-only plant' as material

Soybean is an important model crop for photoperiodic response studies in plants and contributes significantly to the study of plant development and physiology in the past century. Because soybean plant is much bigger in size and longer in life cycle than Arabidopsis, it needs much more space for growth and time for investigation, which significantly hamper the efficiency of research. In the current study, we tested the photoperiodic response of a distinctive artificially-made cotyledon-only plant (COP) using a photoperiod-sensitive soybean variety Zigongdongdou (ZGDD) and other varieties with diverse sensitivity to photoperiod. ZGDD COPs flowered 39.4 ± 2.5 d after emergence under short-day conditions but maintained vegetative growth under long-day and night break conditions, which is similar to the case in the intact ZGDD plants. The COPs of early-maturing and medium-maturing soybean varieties also grew and flowered normally under natural day-length conditions. At the molecular level, the key genes in the photoperiodic pathway such as E1, GmFT1a, GmFT2a, and GmFT5a in the COPs also showed the same photoperiod sensitivity as in the intact plants. In addition, a simpler material of COP with only one cotyledon and root was generated and found to be sensitive to photoperiod as well. Notably, the COPs are only one-fifth the height of intact plants and one-third the maximum diameter of the intact plants grown in chambers 30 d after emergence. Based on COPs, we established a novel experimental system characterized by an entire photoperiodic response and longer longevity of cotyledons in addition to small plant size, ensuring the consistency, reliability, and stability of plant materials. COPs have the potential to be a novel model material for studies of the developmental biology of soybean and other dicots.

CRISPR/Cas9-mediated targeted mutation of the E1 decreases photoperiod sensitivity, alters stem growth habits, and decreases branch number in soybean

The distribution of elite soybean (Glycine max) cultivars is limited due to their highly sensitive to photoperiod, which affects the flowering time and plant architecture. The recent emergence of CRISPR/Cas9 technology has uncovered new opportunities for genetic manipulation of soybean. The major maturity gene E1 of soybean plays a critical role in soybean photoperiod response. Here, we performed CRISPR/Cas9-mediated targeted mutation of E1 gene in soybean cultivar Tianlong1 carrying the dominant E1 to investigate its precise function in photoperiod regulation, especially in plant architecture regulation. Four types of mutations in the E1 coding region were generated. No off-target effects were observed, and homozygous trans-clean mutants without T-DNA were obtained. The photoperiod sensitivity of e1 mutants decreased relative to the wild type plants; however, e1 mutants still responded to photoperiod. Further analysis revealed that the homologs of E1, E1-La, and E1-Lb, were up-regulated in the e1 mutants, indicating a genetic compensation response of E1 and its homologs. The e1 mutants exhibited significant changes in the architecture, including initiation of terminal flowering, formation of determinate stems, and decreased branch numbers. To identify E1-regulated genes related to plant architecture, transcriptome deep sequencing (RNA-seq) was used to compare the gene expression profiles in the stem tip of the wild-type soybean cultivar and the e1 mutants. The expression of shoot identity gene Dt1 was significantly decreased, while Dt2 was significantly upregulated. Also, a set of MADS-box genes was up-regulated in the stem tip of e1 mutants which might contribute to the determinate stem growth habit.

QNE1 is a key flowering regulator determining the length of the vegetative period in soybean cultivars

The soybean E1 gene is a major regulator that plays an important role in flowering time and maturity. However, it remains unclear how cultivars carrying the dominant E1 allele adapt to the higher latitudinal areas of northern China. We mapped the novel quantitative trait locus QNE1 (QTL near E1) for flowering time to the region proximal to E1 on chromosome 6 in two mapping populations. Positional cloning revealed Glyma.06G204300, encoding a TCP-type transcription factor, as a strong candidate gene for QNE1. Association analysis further confirmed that functional single nucleotide polymorphisms (SNPs) at nucleotides 686 and 1,063 in the coding region of Glyma.06G204300 were significantly associated with flowering time. The protein encoded by the candidate gene is localized primarily to the nucleus. Furthermore, soybean and Brassica napus plants overexpressing Glyma.06G204300 exhibited early flowering. We conclude that despite their similar effects on flowering time, QNE1 and E4 may control flowering time through different regulatory mechanisms, based on expression studies and weighted gene co-expression network analysis of flowering time-related genes. Deciphering the molecular basis of QNE1 control of flowering time enriches our knowledge of flowering gene networks in soybean and will facilitate breeding soybean cultivars with broader latitudinal adaptation.

QTL analyses of soybean root system architecture revealed genetic relationships with shoot-related traits

The genetic basis of soybean root system architecture (RSA) and the genetic relationship between shoot and RSA were revealed by integrating data from recombinant inbred population grafting and QTL mapping. Variations in root system architecture (RSA) affect the functions of roots and thus play vital roles in plant adaptations and agricultural productivity. The aim of this study was to unravel the genetic relationship between RSA traits and shoot-related traits in soybean. This study characterized RSA variability at seedling stage in a recombinant inbred population, derived from a cross between cultivated soybean C08 and wild soybean W05, and performed high-resolution quantitative trait locus (QTL) mapping. In total, 34 and 41 QTLs were detected for RSA-related and shoot-related traits, respectively, constituting eight QTL clusters. Significant QTL correspondence was found between shoot biomass and RSA-related traits, consistent with significant correlations between these phenotypes. RSA-related QTLs also overlapped with selection regions in the genome, suggesting the cultivar RSA could be a partial consequence of domestication. Using reciprocal grafting, we confirmed that shoot-derived signals affected root development and the effects were controlled by multiple loci. Meanwhile, RSA-related QTLs were found to co-localize with four soybean flowering-time loci. Consistent with the phenotypes of the parental lines of our RI population, diminishing the function of flowering controlling E1 family through RNA interference (RNAi) led to reduced root growth. This implies that the flowering time-related genes within the RSA-related QTLs are actually contributing to RSA. To conclude, this study identified the QTLs that determine RSA through controlling root growth indirectly via regulating shoot functions, and discovered superior alleles from wild soybean that could be used to improve the root structure in existing soybean cultivars.

Identification of two quantitative genes controlling soybean flowering using bulked-segregant analysis and genetic mapping

Photoperiod responsiveness is important to soybean production potential and adaptation to local environments. Varieties from temperate regions generally mature early and exhibit extremely low yield when grown under inductive short-day (SD) conditions. The long-juvenile (LJ) trait is essentially a reduction and has been introduced into soybean cultivars to improve yield in tropical environments. In this study, we used next-generation sequencing (NGS)-based bulked segregant analysis (BSA) to simultaneously map qualitative genes controlling the LJ trait in soybean. We identified two genomic regions on scaffold_32 and chromosome 18 harboring loci LJ32 and LJ18, respectively. Further, we identified LJ32 on the 228.7-kb scaffold_32 as the soybean pseudo-response-regulator gene Tof11 and LJ18 on a 301-kb region of chromosome 18 as a novel PROTEIN FLOWERING LOCUS T-RELATED gene, Glyma.18G298800. Natural variants of both genes contribute to LJ trait regulation in tropical regions. The molecular identification and functional characterization of Tof11 and LJ18 will enhance understanding of the molecular mechanisms underlying the LJ trait and provide useful genetic resources for soybean molecular breeding in tropical regions.

Soybean GmHY2a encodes a phytochromobilin synthase that regulates internode length and flowering time

Plant height and flowering time are important agronomic traits that directly affect soybean [Glycine max (L.) Merr.] adaptability and yield. Here, the Glycine max long internode 1 (Gmlin1) mutant was selected from an ethyl methyl sulfonate (EMS)-mutated Williams 82 population due to its long internodes and early flowering. Using bulked segregant analysis (BSA), the Gmlin1 locus was mapped to Glyma.02G304700, a homologue of the Arabidopsis HY2 gene, which encodes a phytochromobilin (PΦB) synthase involved in phytochrome chromophore synthesis. Mutation of GmHY2a results in failure of the de-etiolation response under both red and far-red light. The Gmlin1 mutant exhibits a constitutive shade avoidance response under normal light, and the mutations influence the auxin and gibberellin pathways to promote internode elongation. The Gmlin1 mutant also exhibits decreased photoperiod sensitivity. In addition, the soybean photoperiod repressor gene E1 is down-regulated in the Gmlin1 mutant, resulting in accelerated flowering. The nuclear import of phytochrome A (GmphyA) and GmphyB following light treatment is decreased in Gmlin1 protoplasts, indicating that the weak light response of the Gmlin1 mutant is caused by a decrease in functional phytochrome. Together, these results indicate that GmHY2a plays an important role in soybean phytochrome biosynthesis and provide insights into the adaptability of the soybean plant.

Novel and multifaceted regulations of photoperiodic flowering by phytochrome A in soybean

Photoperiod is an important environmental cue. Plants can distinguish the seasons and flower at the right time through sensing the photoperiod. Soybean is a sensitive short-day crop, and the timing of flowering varies greatly at different latitudes, thus affecting yields. Soybean cultivars in high latitudes adapt to the long day by the impairment of two phytochrome genes, PHYA3 and PHYA2, and the legume-specific flowering suppressor, E1. However, the regulating mechanism underlying phyA and E1 in soybean remains largely unknown. Here, we classified the regulation of the E1 family by phyA2 and phyA3 at the transcriptional and posttranscriptional levels, revealing that phyA2 and phyA3 regulate E1 by directly binding to LUX proteins, the critical component of the evening complex, to regulate the stability of LUX proteins. In addition, phyA2 and phyA3 can also directly associate with E1 and its homologs to stabilize the E1 proteins. Therefore, phyA homologs control the core flowering suppressor E1 at both the transcriptional and posttranscriptional levels, to double ensure the E1 activity. Thus, our results disclose a photoperiod flowering mechanism in plants by which the phytochrome A regulates LUX and E1 activity.

Plant clock modifications for adapting flowering time to local environments

During and after the domestication of crops from ancestral wild plants, humans selected cultivars that could change their flowering time in response to seasonal daylength. Continuous selection of this trait eventually allowed the introduction of crops into higher or lower latitudes and different climates from the original regions where domestication initiated. In the past two decades, numerous studies have found the causal genes or alleles that change flowering time and have assisted in adapting crop species such as barley (Hordeum vulgare), wheat (Triticum aestivum L.), rice (Oryza sativa L.), pea (Pisum sativum L.), maize (Zea mays spp. mays), and soybean (Glycine max (L.) Merr.) to new environments. This updated review summarizes the genes or alleles that contributed to crop adaptation in different climatic areas. Many of these genes are putative orthologs of Arabidopsis (Arabidopsis thaliana) core clock genes. We also discuss how knowledge of the clock's molecular functioning can facilitate molecular breeding in the future.

Genomic Dissection and Diurnal Expression Analysis Reveal the Essential Roles of the PRR Gene Family in Geographical Adaptation of Soybean

Pseudo-response regulator (PRR) family members serve as key components of the core clock of the circadian clock, and play important roles in photoperiodic flowering, stress tolerance, growth, and the development of plants. In this study, 14 soybean PRR genes were identified, and classified into three groups according to phylogenetic analysis and structural characteristics. Real-time quantitative PCR analysis revealed that 13 GmPRRs exhibited obvious rhythmic expression under long-day (LD) and short-day (SD) conditions, and the expression of 12 GmPRRs was higher under LD in leaves. To evaluate the effects of natural variations in GmPRR alleles on soybean adaptation, we examined the sequences of GmPRRs among 207 varieties collected across China and the US, investigated the flowering phenotypes in six environments, and analyzed the geographical distributions of the major haplotypes. The results showed that a majority of non-synonymous mutations in the coding region were associated with flowering time, and we found that the nonsense mutations resulting in deletion of the CCT domain were related to early flowering. Haplotype analysis demonstrated that the haplotypes associated with early flowering were mostly distributed in Northeast China, while the haplotypes associated with late flowering were mostly cultivated in the lower latitudes of China. Our study of PRR family genes in soybean provides not only an important guide for characterizing the circadian clock-controlled flowering pathway but also a theoretical basis and opportunities to breed varieties with adaptation to specific regions and farming systems.

Modified photoperiod response of CsFT promotes day neutrality and early flowering in cultivated cucumber

Map-based cloning and photoperiod response detection suggested that CsFT is the critical gene for cucumber photoperiod domestication. Photoperiod sensitivity is important for sensing seasonal changes and local adaptation. However, day-length sensitivity limits crop geographical adaptation and it should be modified during domestication. Cucumber was domesticated in southern Asia and is currently cultivated worldwide across a wide range of latitudes, but its photoperiod sensitivity and its change during cucumber domestication are unknown. Here, we confirmed wild cucumber (Hardwickii) was a short-day plant, and its flowering depends on short-day (SD) conditions, while the cultivated cucumber (9930) is a day-neutral plant that flowers independently of day length. A photoperiod sensitivity locus (ps-1) was identified by the 9930 × Hardwickii F₂ segregating populations, which span a ~ 970 kb region and contain 60 predicted genes. RNA-seq analysis showed that the critical photoperiod pathway gene FLOWERING LOCUS T (CsFT) within the ps-1 locus exhibits differential expression between 9930 and Hardwickii, which was confirmed by qRT-PCR detection. CsFT in Hardwickii was sensitive to day length and could be significantly induced by SD conditions, whereas CsFT was highly expressed in 9930 and was insensitive to day length. Moreover, the role of CsFT in promoting flowering was verified by overexpression of CsFT in Arabidopsis. We also identified the genetic variations existing in the promoter of CsFT among the different geographic cucumbers and suggest they have possible roles in photoperiod domestication. The results of this study suggest that a variation in photoperiod sensitivity of CsFT is associated with day neutrality and early flowering in cultivated cucumber and could contribute to cucumber cultivation in diverse regions throughout the world.

Development of Ethyl Methanesulfonate Mutant Edamame Soybean (Glycine max (L.) Merr.) Populations and Forward and Reverse Genetic Screening for Early-Flowering Mutants

Induced mutation is a viable breeding strategy that is widely utilized in the development of elite plant varieties. We aimed to improve a variety of edamame by constructing novel mutant populations using the ethyl methanesulfonate in soybeans (Glycine max (L.) Merr.). In the M₂ population, the flowering stage showed a considerable standard deviation compared to the wild type, confirming that the mutant populations had the expected DNA mutations. To identify the DNA mutations in the mutant populations, we used the targeting induced local lesions in genomes (TILLING) method, which is a reverse genetic method, to search for soybean flowering-related gene mutants. A total of 30 mutants from E1, E3, E4, and PhyA1 genes, which are known to be highly effective genes, or their homologous gene for flowering and maturation found in soybean quantitative trait locus analyses were isolated from our TILLING screening. Among these mutants, there were eleven nonsynonymous substitution mutants, one nonsense mutant, and two single nucleotide deletion mutants that could be expected to reduce or eliminate gene function. The e1, e3, and e4 mutants obtained in this study flowered considerably earlier than the wild type. In particular, the e1 mutant with a nonsynonymous substitution flowered approximately 1 month after sowing regardless of the sowing date, and its harvest date was approximately 1 month earlier than that of the wild type. Mutations identified using the TILLING method could not only be used as gel-based DNA markers with the same manipulation method, but the mutations could also be detected as DNA markers by the high-resolution melting method. These results indicate that mutations achieved without chromosome modification by crossbreeding are effective for early and practical improvement of superior varieties and that efficient selection of mutants by reverse genetics is an effective method for the identification of genetic modifications. The edamame mutant populations developed in this study are believed to possess various useful alleles which may be applicable in the search for mutations that lead to improved edamame yield and eating quality beyond the flowering stage.

Transcriptome Analysis of Short-Day Photoperiod Inducement in Adzuki Bean (Vigna angularis L.) Based on RNA-Seq

The flowering characteristics of adzuki bean are influenced by several environmental factors. Light is an important ecological factor that induces flowering in adzuki bean, but to date, there have been few reports on the transcriptomic features of photoperiodic regulation of adzuki bean flowering. This study is based on RNA sequencing (RNA-seq) techniques to elucidate the expression of light-related regulatory genes under short-day photoperiod inducement of adzuki bean flowering, providing an important theoretical basis for its accelerated breeding. Short-day photoperiod inducement of 10 h was conducted for 5 day, 10 day, and 15 day periods on "Tang shan hong xiao dou" varieties, which are more sensitive to short-day photoperiod conditions than the other varieties. Plants grown under natural light (14.5 h) for 5 days, 10 days, and 15 days were used as controls to compare the progress of flower bud differentiation and flowering characteristics. The topmost unfolded functional leaves were selected for transcriptome sequencing and bioinformatics analysis. The short-day photoperiod inducement promoted flower bud differentiation and advanced flowering time in adzuki bean. Transcriptomic analysis revealed 5,608 differentially expressed genes (DEGs) for the combination of CK-5d vs. SD-5d, CK-10d vs. SD-10d, and CK-15d vs. SD-15d. The three groups of the DEGs were analyzed using the Gene Ontology (GO) and the Kyoto Encyclopedia of Genomes and Genomes (KEGG) databases; the DEGs were associated with flowering, photosystem, and the circadian rhythm and were mainly concentrated in the hormone signaling and metabolism, circadian rhythm, and antenna protein pathways; So, 13 light-related genes across the three pathways were screened for differential and expression characteristics. Through the functional annotations of orthologs, these genes were related to flowering, which were supposed to be good candidate genes in adzuki bean. The findings provide a deep understanding of the molecular mechanisms of adzuki bean flowering in response to short-day photoperiod inducement, which laid a foundation for the functional verification of genes in the next step, and provide an important reference for the molecular breeding of adzuki bean.

GmMDE genes bridge the maturity gene E1 and florigens in photoperiodic regulation of flowering in soybean

Soybean (Glycine max) is highly sensitive to photoperiod, which affects flowering time and plant architecture and thus limits the distribution range of elite soybean cultivars. The major maturity gene E1 confers the most prominent effect on photoperiod sensitivity, but its downstream signaling pathway remains largely unknown. Here, we confirm that the encoded E1 protein is a transcriptional repressor. The expression of seven GmMDE genes (Glycine max MADS-box genes downregulated by E1) was suppressed when E1 was overexpressed and promoted when E1 was knocked out through clustered regularly-interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9)-mediated mutagenesis. These GmMDEs exhibited similar tissue specificity and expression patterns, including in response to photoperiod, E1 expression, and E1 genotype. E1 repressed GmMDE promoter activity. Results for two GmMDEs showed that E1 epigenetically silences their expression by directly binding to their promoters to increase H3K27me3 levels. The overexpression of GmMDE06 promoted flowering and post-flowering termination of stem growth. The late flowering phenotype of E1-overexpressing soybean lines was reversed by the overexpression of GmMDE06, placing GmMDE06 downstream of E1. The overexpression of GmMDE06 increased the expression of the soybean FLOWERING LOCUS T orthologs GmFT2a and GmFT5a, leading to feedback upregulation of GmMDE, indicating that GmMDE and GmFT2a/GmFT5a form a positive regulatory feedback loop promoting flowering. GmMDE06 also promoted post-flowering termination of stem growth by repressing the expression of the shoot identity gene Dt1. The E1-GmMDEs-GmFT2a/5a-Dt1 signaling pathway illustrates how soybean responds to photoperiod by modulating flowering time and post-flowering stem termination.

Current overview on the genetic basis of key genes involved in soybean domestication

Modern crops were created through the domestication and genetic introgression of wild relatives and adaptive differentiation in new environments. Identifying the domestication-related genes and unveiling their molecular diversity provide clues for understanding how the domesticated variants were selected by ancient people, elucidating how and where these crops were domesticated. Molecular genetics and genomics have explored some domestication-related genes in soybean (Glycine max). Here, we summarize recent studies about the quantitative trait locus (QTL) and genes involved in the domestication traits, introduce the functions of these genes, clarify which alleles of domesticated genes were selected during domestication. A deeper understanding of soybean domestication could help to break the bottleneck of modern breeding by highlighting unused genetic diversity not selected in the original domestication process, as well as highlighting promising new avenues for the identification and research of important agronomic traits among different crop species.

Expanding the gene pool for soybean improvement with its wild relatives

Genetic diversity is a cornerstone of crop improvement, However, cultivated soybean (Glycine max) has undergone several genetic bottlenecks, including domestication in China, the introduction of landraces to other areas of the world and, latterly, selective breeding, leading to low genetic diversity the poses a major obstacle to soybean improvement. By contrast, there remains a relatively high level of genetic diversity in soybean's wild relatives, especially the perennial soybeans (Glycine subgenus Glycine), which could serve as potential gene pools for improving soybean cultivars. Wild soybeans are phylogenetically diversified and adapted to various habitats, harboring resistance to various biotic and abiotic stresses. Advances in genome and transcriptome sequencing enable alleles associated with desirable traits that were lost during domestication of soybean to be discovered in wild soybean. The collection and conservation of soybean wild relatives and the dissection of their genomic features will accelerate soybean breeding and facilitate sustainable agriculture and food production.

Identification of the accessible chromatin regions in six tissues in the soybean

Accessible chromatin regions (ACRs) are tightly associated with gene expressions in the genome. Conserved non-coding cis-regulatory elements, such as transcription factor binding motifs, are usually found in ACRs, indicating an essential regulatory role of ACRs in the plant genome architecture. However, there have been few studies on soybean ACRs, especially those focusing on specific tissues. Hence, in this study, with the convenient ATAC-seq, we identified the ACRs in six soybean tissues, including root, leaf bud, flower, flower bud, developing seed, and pod. In total, the ACRs occupied about 3.3% of the entire soybean genome. By integrating the results from RNA-seq and transcription factor (TF) ChIP-seq, ACRs were found to be tightly associated with gene expressions and TF binding capacities in soybean. Together, these data provide a comprehensive understanding of the genomic features of ACRs in soybean. As a collection of essential genomic resources, these processed data are made available at datahub.wildsoydb.org.

Candidate Genes Modulating Reproductive Timing in Elite US Soybean Lines Identified in Soybean Alleles of Arabidopsis Flowering Orthologs With Divergent Latitude Distribution

Adaptation of soybean cultivars to the photoperiod in which they are grown is critical for optimizing plant yield. However, despite its importance, only the major loci conferring variation in flowering time and maturity of US soybean have been isolated. By contrast, over 200 genes contributing to floral induction in the model organism Arabidopsis thaliana have been described. In this work, putative alleles of a library of soybean orthologs of these Arabidopsis flowering genes were tested for their latitudinal distribution among elite US soybean lines developed in the United States. Furthermore, variants comprising the alleles of genes with significant differences in latitudinal distribution were assessed for amino acid conservation across disparate genera to infer their impact on gene function. From these efforts, several candidate genes from various biological pathways were identified that are likely being exploited toward adaptation of US soybean to various maturity groups.