The soybean stem growth habit gene Dt1 is an ortholog of Arabidopsis TERMINAL FLOWER1

Identification of a major locus for flowering pattern sheds light on plant architecture diversification in cultivated peanut

A major gene controls flowering pattern in peanut, possibly encoding a TFL1-like. It was subjected to gain/loss events of a deletion and changes in mRNA expression levels, partly explaining the evolution of flowering pattern in Arachis. Flowering pattern (FP) is a major characteristic differentiating the two subspecies of cultivated peanut (Arachis hypogaea L.). Subsp. fastigiata possessing flowers on the mainstem (MSF) and a sequential FP, whereas subsp. hypogaea lacks MSF and exhibits an alternate FP. FP is considered the main contributor to plant adaptability, and evidence indicates that its diversification occurred during the several thousand years of domestication. However, the genetic mechanism that controls FP in peanut is unknown. We investigated the genetics of FP in a recombinant inbred population, derivatives of an A. hypogaea by A. fastigiata cross. Lines segregated 1:1 for FP, indicating a single gene effect. Using Axiom_Arachis2 SNP-array, FP was mapped to a small segment in chromosome B02, wherein a Terminal Flowering 1-like (AhTFL1) gene with a 1492 bp deletion was found in the fastigiata line, leading to a truncated protein. Remapping FP in the RIL population with the AhTFL1 indel as a marker increased the LOD score from 53.3 to 158.8 with no recombination in the RIL population. The same indel was found co-segregating with the phenotype in two independent EMS-mutagenized M2 families, suggesting a hotspot for gene conversion. Also, AhTFL1 was significantly less expressed in the fastigiata line compared to hypogaea and in flowering than non-flowering branches. Sequence analysis of the AhTFL1 in peanut world collections indicated significant conservation, supporting the putative role of AhTFL1 in peanut speciation during domestication and modern cultivation.

A functionally divergent SOC1 homolog improves soybean yield and latitudinal adaptation

Soybean (Glycine max) grows in a wide range of latitudes, but it is extremely sensitive to photoperiod, which reduces its yield and ability to adapt to different environments. Therefore, understanding of the genetic basis of soybean adaptation is of great significance for breeding and improvement. Here, we characterized Tof18 (SOC1a) that conditions early flowering and growth habit under both short-day and long-day conditions. Molecular analysis confirmed that the two SOC1 homologs present in soybeans (SOC1a and SOC1b) underwent evolutionary functional divergence, with SOC1a having stronger effects on flowering time and stem node number than SOC1b due to transcriptional differences. soc1a soc1b double mutants showed stronger functional effects than either of the single mutants, perhaps due to the formation of SOC1a and SOC1b homodimers or heterodimers. Additionally, Tof18/SOC1a improves the latitudinal adaptation of cultivated soybeans, highlighting the functional importance of SOC1a. The Tof18G allele facilitates adaptation to high latitudes, whereas Tof18A facilitates adaptation to low latitudes. We demonstrated that SOC1s contribute to floral induction in both leaves and shoot apex through inter-regulation with FTs. The SOC1a-SOC1b-Dt2 complex plays essential roles in stem growth habit by directly binding to the regulatory sequence of Dt1, making the genes encoding these proteins potential targets for genome editing to improve soybean yield via molecular breeding. Since the natural Tof18A allele increases node number, introgressing this allele into modern cultivars could improve yields, which would help optimize land use for food production in the face of population growth and global warming.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes

INTRODUCTION

Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery.

OBJECTIVES

Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association.

METHODS

We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes.

RESULTS

The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene.

CONCLUSION

The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

BACKGROUND

With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research.

RESULTS

To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection.

CONCLUSION

The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

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

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

Identification and Characterization of Key Genes Responsible for Weedy and Cultivar Growth Types in Soybean

In cultivated plants, shoot morphology is an important factor that influences crop economic value. However, the effects of gene expression patterns on shoot morphology are not clearly understood. In this study, the molecular mechanism behind shoot morphology (including leaf, stem, and node) was analyzed using RNA sequencing to compare weedy (creeper) and cultivar (stand) growth types obtained in F7 derived from a cross of wild and cultivated soybeans. A total of 12,513 (in leaves), 14,255 (in stems), and 11,850 (in nodes) differentially expressed genes were identified among weedy and cultivar soybeans. Comparative transcriptome and expression analyses revealed 22 phytohormone-responsive genes. We found that GIBBERELLIN 2-OXIDASE 8 (GA2ox), SPINDLY (SPY), FERONIA (FER), AUXIN RESPONSE FACTOR 8 (ARF8), CYTOKININ DEHYDROGENASE-1 (CKX1), and ARABIDOPSIS HISTIDINE KINASE-3 (AHK3), which are crucial phytohormone response genes, were mainly regulated in the shoot of weedy and cultivar types. These results indicate that interactions between phytohormone signaling genes regulate shoot morphology in weedy and cultivar growth type plants. Our study provides insights that are useful for breeding and improving crops to generate high-yield soybean varieties.

The Genetic and Hormonal Inducers of Continuous Flowering in Orchids: An Emerging View

Orchids are the flowers of magnetic beauty. Vivid and attractive flowers with magnificent shapes make them the king of the floriculture industry. However, the long-awaited flowering is a drawback to their market success, and therefore, flowering time regulation is the key to studies about orchid flower development. Although there are some rare orchids with a continuous flowering pattern, the molecular regulatory mechanisms are yet to be elucidated to find applicable solutions to other orchid species. Multiple regulatory pathways, such as photoperiod, vernalization, circadian clock, temperature and hormonal pathways are thought to signalize flower timing using a group of floral integrators. This mini review, thus, organizes the current knowledge of floral time regulators to suggest future perspectives on the continuous flowering mechanism that may help to plan functional studies to induce flowering revolution in precious orchid species.

Progress in soybean functional genomics over the past decade

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

CRISPR gene editing of major domestication traits accelerating breeding for Solanaceae crops improvement

Domestication traits particularly fruit size and plant architecture and flowering are critical in transforming a progenitor's wild stature into a super improved plant. The latest advancements in the CRISPR system, as well as its rapid adoption, are speeding up plant breeding. Solanaceae has a varied range of important crops, with a few model crops, such as tomato and, more recently, groundcherry, serving as a foundation for developing molecular techniques, genome editing tools, and establishing standards for other crops. Domestication traits in agricultural plants are quantified and widely adopted under modern plant breeding to improve small-fruited and bushy crop species like goji berry. The molecular mechanisms of the FW2.2, FW3.2, FW11.3, FAS/CLV3, LC/WUS, SP, SP5G, and CRISPR genome editing technology have been described in detail here. Furthermore, special focus has been placed on CRISPR gene editing achievements for revolutionizing Solanaceae breeding and changing the overall crop landscape. This review seeks to provide a thorough overview of the CRISPR technique's ongoing advancements, particularly in Solanaceae, in terms of domesticated features, future prospects, and regulatory risks. We believe that this vigorous discussion will lead to a broader understanding of CRISPR gene editing as a tool for achieving key breeding goals in other Solanaceae minor crops with significant industrial value.

Detecting QTL and Candidate Genes for Plant Height in Soybean via Linkage Analysis and GWAS

Soybean is an important global crop for edible protein and oil, and plant height is a main breeding goal which is closely related to its plant shape and yield. In this research, a high-density genetic linkage map was constructed by 1996 SNP-bin markers on the basis of a recombinant inbred line population derived from Dongnong L13 × Henong 60. A total of 33 QTL related to plant height were identified, of which five were repeatedly detected in multiple environments. In addition, a 455-germplasm population with 63,306 SNP markers was used for multi-locus association analysis. A total of 62 plant height QTN were detected, of which 26 were detected repeatedly under multiple methods. Two candidate genes, Glyma.02G133000 and Glyma.05G240600, involving in plant height were predicted by pathway analysis in the regions identified by multiple environments and backgrounds, and validated by qRT-PCR. These results enriched the soybean plant height regulatory network and contributed to molecular selection-assisted breeding.

Green revolution to grain revolution: Florigen in the frontiers

Burgeoning human population dents, globally, the brimming buffer stock as well as gain in food grain production. However, an imminent global starvation was averted through precise scientific intervention and pragmatic policy changes in the 1960s and was eulogized as the "Green Revolution". Miracle rice and wheat obtained through morphometric changes in the ideotype of these two crops yielded bumper harvest that nucleated in Asia and translated into Latin America. The altered agronomic traits in these two crops were the result of tinkering with the phyto-hormone "Gibberellin'. Recently, another plant hormone 'Cytokinin' has gained prominence for its involvement in the grain revolution in rice and other field crops. Suo moto homeostasis of CK by the cytokinin oxidase enzyme governs the cardinal shoot apical meristem that produces new flowering primordia thereby enhancing grain number. Similarly, the flowering hormone 'Florigen' impacts sympodia formation, flowering, and fruit production in tomato. The role of heterozygosity induced heterosis by florigen in revolutionizing tomato production and cellular homeostasis of CK by CK oxidising enzyme (CKX) in enhancing rice production has been path-breaking. This review highlights role of phytohormones in grain revolution and crop specific fine-tuning of gibberellins, cytokinins and florigen to accomplish maximum yield potential in field crops.

Genetic control and allele variation among soybean maturity groups 000 through IX

Soybean [Glycinemax (L.) Merr.] maturity determines the growing region of a given soybean variety and is a primary factor in yield and other agronomic traits. The objectives of this research were to identify the quantitative trait loci (QTL) associated with maturity groups (MGs) and determine the genetic control of soybean maturity in each MG. Using data from 16,879 soybean accessions, genome-wide association (GWA) analyses were conducted for each paired MG and across MGs 000 through IX. Genome-wide association analyses were also performed using 184 genotypes (MGs V-IX) with days to flowering (DTF) and maturity (DTM) collected in the field. A total of 58 QTL were identified to be significantly associated with MGs in individual GWAs, which included 12 reported maturity loci and two stem termination genes. Genome-wide associations across MGs 000-IX detected a total of 103 QTL and confirmed 54 QTL identified in the individual GWAs. Of significant loci identified, qMG-5.2 had effects on the highest number (9) of MGs, followed by E2, E3, Dt2, qMG-15.5, E1, qMG-13.1, qMG-7.1, and qMG-16.1, which affected five to seven MGs. A high number of genetic loci (8-25) that affected MGs 0-V were observed. Stem termination genes Dt1 and Dt2 mainly had significant allele variation in MGs II-V. Genome-wide associations for DTF, DTM, and reproductive period (RP) in the diversity panel confirmed 15 QTL, of which seven were observed in MGs V-IX. The results generated can help soybean breeders manipulate the maturity loci for genetic improvement of soybean yield.

Overcoming the genetic compensation response of soybean florigens to improve adaptation and yield at low latitudes

The classical soybean (Glycine max) trait long juvenile (LJ) is essentially a reduction in sensitivity to short-day (SD) conditions for induction and completion of flowering, and has been introduced into soybean cultivars to improve yield in tropical environments. However, only one locus, J, is known to confer LJ in low-latitude varieties. Here, we defined two quantitative trait loci contributing to the LJ trait, LJ16.1 and LJ16.2, and identified them as the florigen (FT) homologs FT2a and FT5a, respectively. The two selected florigen variations both delay flowering time under SD conditions by repressing the floral meristem identity gene GmAPETALA1. Single mutants have a relatively subtle effect on flowering time and displayed a substantial genetic compensation response, but this was absent in ft2a ft5a double mutants, which showed an enhanced LJ phenotype that translated to higher yields under SD conditions. A survey of sequence diversity suggests that FT2a and FT5a variants have diverse origins and have played distinct roles as soybean spread to lower latitudes. Our results show that integration of variants in the florigen genes offers a strategy for customizing flowering time to adjust adaptation and improve crop productivity in tropical regions.

Overexpression of four MiTFL1 genes from mango delays the flowering time in transgenic Arabidopsis

BACKGROUND

TERMINAL FLOWER 1 (TFL1) belongs to the phosphatidylethanolamine-binding protein (PEBP) family, which is involved in inflorescence meristem development and represses flowering in several plant species. In the present study, four TFL1 genes were cloned from the mango (Mangifera indica L.) variety 'SiJiMi' and named MiTFL1-1, MiTFL1-2, MiTFL1-3 and MiTFL1-4.

RESULTS

Sequence analysis showed that the encoded MiTFL1 proteins contained a conserved PEBP domain and belonged to the TFL1 group. Expression analysis showed that the MiTFL1 genes were expressed in not only vegetative organs but also reproductive organs and that the expression levels were related to floral development. Overexpression of the four MiTFL1 genes delayed flowering in transgenic Arabidopsis. Additionally, MiTFL1-1 and MiTFL1-3 changed the flower morphology in some transgenic plants. Yeast two-hybrid (Y2H) analysis showed that several stress-related proteins interacted with MiTFL1 proteins.

CONCLUSIONS

The four MiTFL1 genes exhibited a similar expression pattern, and overexpression in Arabidopsis resulted in delayed flowering. Additionally, MiTFL1-1 and MiTFL1-3 overexpression affected floral organ development. Furthermore, the MiTFL1 proteins could interact with bHLH and 14-3-3 proteins. These results indicate that the MiTFL1 genes may play an important role in the flowering process in mango.

Increased copy number of gibberellin 2-oxidase 8 genes reduced trailing growth and shoot length during soybean domestication

Copy number variations (CNVs) play important roles in crop domestication. However, there is only very limited information on the involvement of CNVs in soybean domestication. Trailing growth and long shoots are soybean adaptations for natural habitats but cause lodging that hampers yield in cultivation. Previous studies have focused on Dt1/2 affecting the indeterminate/determinate growth habit, whereas the possible role of the gibberellin pathway remained unclear. In the present study, quantitative trait locus (QTL) mapping of a recombinant inbred population of 460 lines revealed a trailing-growth-and-shoot-length QTL. A CNV region within this QTL was identified, featuring the apical bud-expressed gibberellin 2-oxidase 8A/B, the copy numbers of which were positively correlated with expression levels and negatively with trailing growth and shoot length, and their effects were demonstrated by transgenic soybean and Arabidopsis thaliana. Based on the fixation index, this CNV region underwent intense selection during the initial domestication process.

Meta-GWAS for quantitative trait loci identification in soybean

We report a meta-Genome Wide Association Study involving 73 published studies in soybean [Glycine max L. (Merr.)] covering 17,556 unique accessions, with improved statistical power for robust detection of loci associated with a broad range of traits. De novo GWAS and meta-analysis were conducted for composition traits including fatty acid and amino acid composition traits, disease resistance traits, and agronomic traits including seed yield, plant height, stem lodging, seed weight, seed mottling, seed quality, flowering timing, and pod shattering. To examine differences in detectability and test statistical power between single- and multi-environment GWAS, comparison of meta-GWAS results to those from the constituent experiments were performed. Using meta-GWAS analysis and the analysis of individual studies, we report 483 peaks at 393 unique loci. Using stringent criteria to detect significant marker-trait associations, 59 candidate genes were identified, including 17 agronomic traits loci, 19 for seed-related traits, and 33 for disease reaction traits. This study identified potentially valuable candidate genes that affect multiple traits. The success in narrowing down the genomic region for some loci through overlapping mapping results of multiple studies is a promising avenue for community-based studies and plant breeding applications.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

Rapid excavating a FLOWERING LOCUS T-regulator NF-YA using genotyping-by-sequencing

Soybean (Glycine max (L.) Merrill) is one of the most important crop plants in the world as an important source of protein for both human consumption and livestock fodder. As flowering time contributes to yield, finding new QTLs and further identifying candidate genes associated with various flowering time are fundamental to enhancing soybean yield. In this study, a set of 120 recombinant inbred lines (RILs) which was developed from a cross of two soybean cultivars, Suinong4 (SN4) and ZK168, were genotyped by genotyping-by-sequencing (GBS) approach and phenotyped to expand the cognitive of flowering time by quantitative trait loci (QTL) analysis. Eventually, three stable QTLs related to flowering time which were detected separately located on chromosome 14, 18, and 19 under long-day (LD) conditions. We predicted candidate genes for each QTL and carried out association analyses between the putative causal alleles and flowering time. Moreover, a transient transfection assay was performed and showed that NUCLEAR FACTOR YA 1b (GmNF-YA1b) as a strong candidate for the QTL on chromosome 19 might affect flowering time by suppressing the expression of FLOWERING LOCUS T (GmFT) genes in soybean. QTLs detected in this study would provide fundamental resources for finding candidate genes and clarify the mechanisms of flowering which would be helpful for breeding novel high-yielding soybean cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01237-w.

Natural variation of the Dt2 promoter controls plant height and node number in semi-determinant soybean

Soybean (Glycine max (L.) Merrill) is an important legume crop worldwide. Plant height (PH) is a quantitative trait that is closely related to node number (NN) and internode length (IL) on the main stem, which together affect soybean yield. To identify candidate genes controlling these three traits in soybean, we examined a recombinant inbred line (RIL) population derived from a cross between two soybean varieties with semi-determinate stems (Dt1Dt1Dt2Dt2), JKK378 and HXW. A quantitative trait locus (QTL) named qPH18 was identified that simultaneously controls PH, NN, and IL; this region harbors the semi-determinant gene Dt2. Sequencing of the Dt2 promoter from JKK378 identified three polymorphisms relative to HXW, including two single nucleotide polymorphism (SNPs) and an 18-bp insertion/deletion polymorphism (Indel). Dt2 expression was lower in the qPH18JKK378 group than in the qPH18HXW group, whereas the expression level of the downstream gene Dt1 showed the opposite tendency. A transient transfection assay confirmed that Dt2 promoter activity is lower in JKK378 compared to HXW. We propose that the polymorphisms in the dominant Dt2 promoter underlie the differences in Dt2 expression and its downstream gene Dt1 in the two parents, thereby affecting PH, NN, IL, and grain weight per plant without altering stem growth habit. Compared to the PH18HXW allele, the qPH18JKK378 allele suppresses Dt2 expression, which releases the inhibition of Dt1 expression, thus enhancing NN and grain yield. Our findings shed light on the mechanism underlying NN and PH in soybean and provide a molecular marker to facilitate breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01235-y.

FT5a interferes with the Dt1-AP1 feedback loop to control flowering time and shoot determinacy in soybean

Flowering time and stem growth habit determine inflorescence architecture in soybean, which in turn influences seed yield. Dt1, a homolog of Arabidopsis TERMINAL FLOWER 1 (TFL1), is a major controller of stem growth habit, but its underlying molecular mechanisms remain unclear. Here, we demonstrate that Dt1 affects node number and plant height, as well as flowering time, in soybean under long-day conditions. The bZIP transcription factor FDc1 physically interacts with Dt1, and the FDc1-Dt1 complex directly represses the expression of APETALA1 (AP1). We propose that FT5a inhibits Dt1 activity via a competitive interaction with FDc1 and directly upregulates AP1. Moreover, AP1 represses Dt1 expression by directly binding to the Dt1 promoter, suggesting that AP1 and Dt1 form a suppressive regulatory feedback loop to determine the fate of the shoot apical meristem. These findings provide novel insights into the roles of Dt1 and FT5a in controlling the stem growth habit and flowering time in soybean, which determine the adaptability and grain yield of this important crop.

Identification of QTNs and their candidate genes for flowering time and plant height in soybean using multi-locus genome-wide association studies

Flowering time (FT) and plant height (PH) are important agronomic traits in soybean. However, their genetic foundations are not fully understood. Thus, in this study, a total of 106,013 single nucleotide polymorphisms in 286 soybean accessions were used to associate with the first and full FT (FT1 and FT2) and PH in 4 environments and their BLUP values using 6 multi-locus genome-wide association study methods. As a result, 38, 43, and 27 stable quantitative trait nucleotides (QTNs) were identified, respectively, for FT1, FT2, and PH across at least 3 methods and/or environments. Among these QTNs for FT1, FT2, and PH, 31, 36, and 21 were found to have significant phenotype differences across 2 alleles; 22, 18, and 13 were consistent with the corresponding loci in previous studies; 13 and 8 genes, with more than average expression level, around 64 FT and 27 PH QTNs were predicted as their corresponding candidate genes. Among these candidate genes, GmPRR3b, and GmGIa for FT, and GmTFL1b for PH were known, while some were new, e.g., GmPHYA4, GmVRN5, GmFPA, and GmSPA1 for FT, and Glyma.02g300200, GmFPA, and Glyma.13g339800 for PH. All the validated QTNs were used to design the best cross-combinations in 2 FT directions. In each FT direction, the best 5 cross-combinations were predicted, such as Heihe 54 × Qincha 1 for early FT, and Yingdejiadou × Wuhuabayuehuang for late FT. This study provides solid foundations for genetic basis, molecular biology, and breeding by design of soybean FT and PH.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01230-3.

Diversification in Functions and Expressions of Soybean FLOWERING LOCUS T Genes Fine-Tunes Seasonal Flowering

The proper timing of flowering in response to environmental changes is critical for ensuring crop yields. FLOWERING LOCUS T (FT) homologs of the phosphatidylethanolamine-binding protein family play important roles as floral integrators in many crops. In soybean, we identified 17 genes of this family, and characterized biological functions in flowering for ten FT homologs. Overexpression of GmFT homologs in Arabidopsis revealed that a set of GmFT homologs, including GmFT2a/2b, GmFT3a/3b, and GmFT5a/5b, promoted flowering similar to FT; in contrast, GmFT1a/1b, GmFT4, and GmFT6 delayed flowering. Consistently, expressions of GmFT2a, GmFT2b, and GmFT5a were induced in soybean leaves in response to floral inductive short days, whereas expressions of GmFT1a and GmFT4 were induced in response to long days. Exon swapping analysis between floral activator GmFT2a and floral repressor GmFT4 revealed that the segment B region in the fourth exon is critical for their antagonistic functions. Finally, expression analysis of GmFT2a, GmFT5a, and GmFT4 in soybean accessions exhibiting various flowering times indicated that the mRNA levels of GmFT2a and GmFT5a were higher in early flowering accessions than in late-flowering accessions, while GmFT4 showed the opposite pattern. Moreover, the relative mRNA levels between GmFT2a/GmFT5a and GmFT4 was important in determining day length-dependent flowering in soybean accessions. Taken together, our results suggest that the functions of GmFT homologs have diversified into floral activators and floral repressors during soybean evolution, and the timing of flowering in response to changing day length is determined by modulating the activities of antagonistic GmFT homologs.

Overexpression of GmMYB14 improves high-density yield and drought tolerance of soybean through regulating plant architecture mediated by the brassinosteroid pathway

MYB transcription factors (TFs) have been reported to regulate the biosynthesis of secondary metabolites, as well as to mediate plant adaption to abiotic stresses, including drought. However, the roles of MYB TFs in regulating plant architecture and yield potential remain poorly understood. Here, we studied the roles of the dehydration-inducible GmMYB14 gene in regulating plant architecture, high-density yield and drought tolerance through the brassinosteroid (BR) pathway in soybean. GmMYB14 was shown to localize to nucleus and has a transactivation activity. Stable GmMYB14-overexpressing (GmMYB14-OX) transgenic soybean plants displayed a semi-dwarfism and compact plant architecture associated with decreased cell size, resulting in a decrease in plant height, internode length, leaf area, leaf petiole length and leaf petiole angle, and improved yield in high density under field conditions. Results of the transcriptome sequencing suggested the involvement of BRs in regulating GmMYB14-OX plant architecture. Indeed, GmMYB14-OX plants showed reduced endogenous BR contents, while exogenous application of brassinolide could partly rescue the phenotype of GmMYB14-OX plants. Furthermore, GmMYB14 was shown to directly bind to the promoter of GmBEN1 and up-regulate its expression, leading to reduced BR content in GmMYB14-OX plants. GmMYB14-OX plants also displayed improved drought tolerance under field conditions. GmBEN1 expression was also up-regulated in the leaves of GmMYB14-OX plants under polyethylene glycol treatment, indicating that the GmBEN1-mediated reduction in BR level under stress also contributed to drought/osmotic stress tolerance of the transgenic plants. Our findings provided a strategy for stably increasing high-density yield and drought tolerance in soybean using a single TF-encoding gene.

Identifying new variation at the J locus, previously identified as e6, in long juvenile 'Paranagoiana' soybean

KEY MESSAGE

A previously identified soybean maturity locus, E6, is discovered to be J, with the long juvenile allele in Paranagoiana now deemed j-x. Soybean grown at latitudes of ~20° or lower can produce lower grain yields due to the short days. This limitation can be overcome by using the long juvenile trait (LJ) which delays flowering under short day conditions. Two LJ loci have been mapped to the same location on Gm04, J and E6. The objective of this research was to investigate the e6 allele in 'Paranagoiana' and determine if E6 and J are the same locus or linked loci. KASP markers showed that e6 lines did not have the j-1 allele of LJ PI 159925. A population fixed for E1 but segregating for E6, with e6 introgressed from Paranagoiana, showed single gene control for flowering and maturity under short days. Sequencing Glyma.04G050200, the J gene, with long amplification Taq found that the e6 line 'Paranagoiana' contains a Ty1-copia retrotransposon of ~10,000 bp, inserted within exon 4. PCR amplification of the cDNA of Glyma.04G050200 also showed differences between the mRNA sequences (presence of insertion in j-x). Hence, we conclude that the loci E6 and J are one locus and deem this new variation found in Paranagoiana as j-x.

De novo domestication of wild species to create crops with increased resilience and nutritional value

Creating crops with resistance to drought, soil salinity and insect damage, that simultaneously have higher nutritional quality, is challenging to conventional breeding due to the complex and diffuse genetic basis of those traits. Recent advances in gene editing technology, such as base editors and prime-editing, coupled with a deeper understanding of the genetic basis of domestication delivered by the analysis of crop 'pangenomes', open the exciting prospect of creating novel crops via manipulation of domestication-related genes in wild species. A de novo domestication platform may allow rapid and precise conversion of crop wild relatives into crops, while retaining many of the valuable resilience and nutritional traits left behind during domestication and breeding. Using the Solanaceae family as case in point, we discuss how such a knowledge-driven pipeline could be exploited to contribute to food security over the coming decades.

Identification of a Locus Controlling Compound Raceme Inflorescence in Mungbean [Vigna radiata (L.) R. Wilczek]

Mungbean [Vigna radiata (L.) R. Wilczek] produces a compound raceme inflorescence that branches into secondary inflorescences, which produce flowers. This architecture results in the less-domesticated traits of asynchronous pod maturity and multiple harvest times. This study identified the genetic factors responsible for the compound raceme of mungbean, providing a unique biological opportunity to improve simultaneous flowering. Using a recombinant inbred line (RIL) population derived from VC1973A, an elite cultivar with a compound raceme type, and IT208075, a natural mutant with a simple raceme type, a single locus that determined the inflorescence type was identified based on 1:1 segregation ratio in the F8 generation, and designated Comraceme. Linkage map analysis showed Comraceme was located on chromosome 4 within a marker interval spanning 520 kb and containing 64 genes. RILs carrying heterozygous fragments around Comraceme produced compound racemes, indicating this form was dominant to the simple raceme type. Quantitative trait loci related to plant architecture and inflorescence have been identified in genomic regions of soybean syntenic to Comraceme. In IT208075, 15 genes were present as distinct variants not observed in other landrace varieties or wild mungbean. These genes included Vradi04g00002481, a development-related gene encoding a B3 transcriptional factor. The upstream region of Vradi04g00002481 differed between lines producing the simple and compound types of raceme. Expression of Vradi04g00002481 was significantly lower at the early vegetative stage and higher at the early reproductive stage, in IT208075 than in VC1973A. Vradi04g00002481 was therefore likely to determine inflorescence type in mungbean. Although further study is required to determine the functional mechanism, this finding provides valuable genetic information for understanding the architecture of the compound raceme in mungbean.

Evolutionary QTL-allele changes in main stem node number among geographic and seasonal subpopulations of Chinese cultivated soybeans

The main stem node number (MSN) is a trait related to geographic adaptation, plant architecture and yield potential of soybean. The QTL-allele constitution of the Chinese Cultivated Soybean Population (CCSP) was identified using the RTM-GWAS (restricted two-stage multi-locus genome-wide association study) procedure, from which a QTL-allele matrix was established and then separated into submatrices to explore the genetic structure, evolutionary differentiation, breeding potential and candidate genes of MSN in CCSP. The MSN of 821 accessions varied from 8.8 to 31.1, with an average of 16.3 in Nanjing, China (32.07° N, 118.62° E), where the MSNs of all the materials could be evaluated in a standardized manner. Among the six geo-seasonal subpopulations, the MSN varied from 21.7 in a southern summer-autumn-sowing subpopulation (SA-IV) down to 13.5 in a northeastern spring-sowing subpopulation (SP-I). The materials were genotyped with restriction site-associated DNA-sequencing. Totally 142 main-effect QTLs (73.24% new) with 560 alleles contributing 72.98% to the phenotypic variance were identified. The evolutionary QTL-allele changes in MSN from SA-IV through SP-I showed that inheritance (78.93% of alleles) was the primary factor influencing the evolution of this trait, followed by allele emergence (19.64% alleles), allele exclusion (1.43% alleles), and recombination among retained alleles. In the evolutionary changes, 70 QTLs, including 12 newly emerged QTLs, with 118 alleles were involved. An increase potential of 2-8 nodes was predicted and 112 candidate genes were annotated and preliminarily verified with χ2-tests in the CCSP. The RTM-GWAS showed powerful in detecting QTL-allele system, assessing evolution factors, predicting optimal crosses and identifying candidate genes in a germplasm population.

Allelic characterization and protein structure analysis reveals the involvement of splice site mutation for growth habit differences in Lablab purpureus (L.) Sweet

BACKGROUND

Interrelationship between growth habit and flowering played a key role in the domestication history of pulses; however, the actual genes responsible for these traits have not been identified in Indian bean. Determinate growth habit is desirable due to its early flowering, photo-insensitivity, synchronous pod maturity, ease in manual harvesting and short crop duration. The present study aimed to identify, characterize and validate the gene responsible for growth habit by using a candidate gene approach coupled with sequencing, multiple sequence alignment, protein structure prediction and binding pocket analysis.

RESULTS

Terminal flowering locus was amplified from GPKH 120 (indeterminate) and GNIB-21 (determinate) using the primers designed from PvTFL1y locus of common bean. Gene prediction revealed that the length of the third and fourth exons differed between the two alleles. Allelic sequence comparison indicated a transition from guanine to adenine at the end of the third exon in GNIB 21. This splice site single-nucleotide polymorphism (SNP) was validated in germplasm lines by sequencing. Protein structure analysis indicated involvement of two binding pockets for interaction of terminal flowering locus (TFL) protein with other proteins.

CONCLUSION

The splice site SNP present at the end of the third exon of TFL locus is responsible for the transformation of shoot apical meristem into a reproductive fate in the determinate genotype GNIB 21. The splice site SNP leads to absence of 14 amino acids in mutant TFL protein of GNIB 21, rendering the protein non-functional. This deletion disturbed previously reported anion-binding pocket and secondary binding pocket due to displacement of small β-sheet away from an external loop. This finding may enable the modulation of growth habit in Indian bean and other pulse crops through genome editing.

Environmental and genetic regulation of plant height in soybean

BACKGROUND

Shoot architecture is fundamentally crucial to crop growth and productivity. As a key component of shoot architecture, plant height is known to be controlled by both genetic and environmental factors, though specific details remain scarce.

RESULTS

In this study, 308 representative soybean lines from a core collection and 168 F9 soybean progeny were planted at distinct field sites. The results demonstrated the presence of significant genotype × environment interaction (G × E) effects on traits associated with plant height in a natural soybean population. In total, 19 loci containing 51 QTLs (quantitative trait locus) for plant height were identified across four environments, with 23, 13 and 15 being QTLs for SH (shoot height), SNN (stem node number) and AIL (average internode length), respectively. Significant LOD ranging from 2.50 to 16.46 explained 2.80-26.10% of phenotypic variation. Intriguingly, only two loci, Loc11 and Loc19-1, containing 20 QTLs, were simultaneously detected across all environments. Results from Pearson correlation analysis and PCA (principal component analysis) revealed that each of the five agro-meteorological factors and four soil properties significantly affected soybean plant height traits, and that the corresponding QTLs had additive effects. Among significant environmental factors, AD (average day-length), AMaT (average maximum temperature), pH, and AN (available nitrogen) had the largest impacts on soybean plant height. Therefore, in spite of uncontrollable agro-meteorological factors, soybean shoot architecture might be remolded through combined efforts to produce superior soybean genetic materials while also optimizing soil properties.

CONCLUSIONS

Overall, the comprehensive set of relationships outlined herein among environment factors, soybean genotypes and QTLs in effects on plant height opens new avenues to explore in work aiming to increase soybean yield through improvements in shoot architecture.

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

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

Analysis of Shoot Architecture Traits in Edamame Reveals Potential Strategies to Improve Harvest Efficiency

Edamame is a type of green, vegetable soybean and improving shoot architecture traits for edamame is important for breeding of high-yield varieties by decreasing potential loss due to harvesting. In this study, we use digital imaging technology and computer vision algorithms to characterize major traits of shoot architecture for edamame. Using a population of edamame PIs, we seek to identify underlying genetic control of different shoot architecture traits. We found significant variations in the shoot architecture of the edamame lines including long-skinny and candle stick-like structures. To quantify the similarity and differences of branching patterns between these edamame varieties, we applied a topological measurement called persistent homology. Persistent homology uses algebraic geometry algorithms to measure the structural similarities between complex shapes. We found intriguing relationships between the topological features of branching networks and pod numbers in our plant population, suggesting combination of multiple topological features contribute to the overall pod numbers on a plant. We also identified potential candidate genes including a lateral organ boundary gene family protein and a MADS-box gene that are associated with the pod numbers. This research provides insight into the genetic regulation of shoot architecture traits and can be used to further develop edamame varieties that are better adapted to mechanical harvesting.

A novel mutation in TFL1 homolog sustaining determinate growth in cucumber (Cucumis sativus L.)

BSA-seq combined with whole-genome resequencing map-based cloning delimited the cucumber det-novel locus into a 44.5 kb region in chromosome 6 harboring a putative candidate gene encoding a phosphatidylethanolamine-binding protein (CsCEN). Determinate and indeterminate growth habits of cucumber can affect plant architecture and crop yield. The TERMINAL FLOWER 1 (TFL1) controls determinate/indeterminate growth in Arabidopsis. In this study, a novel mutation in cucumber TFL1 homolog (CsCEN) has shown to regulate determinate growth and product of terminal flowers in cucumber (Cucumis sativus L.), which is similar to the function of CsTFL1 as previously reported. Genetic analysis in two determinate genotypes (D226 and D082) and indeterminate genotype (CCMC) revealed that a single recessive gene is responsible for this determinate growth trait. With the combination of BSA-seq and whole-genome resequencing, the locus of determinate-novel (det-novel) trait was mapped to a 44.5 kb genomic region in chromosome 6. Sequence alignment identified one non-synonymous SNP mutation (A to T) in the third exon of CsCEN, resulting in an amino acid substitution (Thr to Pro), suggesting that determinate growth might be controlled by a novel gene CsCEN (Csa6G152360) which differed from the reported CsTFL1 gene. The CsCEN expression level in shoot apexes and axillary buds was significantly lower in D226 compared to CCMC, suggesting its essential role in sustaining indeterminate growth habit. Identification and characterization of the CsCEN in the present study provide a new insight into plant architecture modification and development of cucumber cultivars suited to mechanized production system.

Soybean AP1 homologs control flowering time and plant height

Flowering time and plant height are key agronomic traits that directly affect soybean (Glycine max) yield. APETALA1 (AP1) functions as a class A gene in the ABCE model for floral organ development, helping to specify carpel, stamen, petal, and sepal identities. There are four AP1 homologs in soybean, all of which are mainly expressed in the shoot apex. Here, we used clustered regularly interspaced short palindromic repeats (CRISPR) - CRISPR-associated protein 9 technology to generate a homozygous quadruple mutant, gmap1, with loss-of-function mutations in all four GmAP1 genes. Under short-day (SD) conditions, the gmap1 quadruple mutant exhibited delayed flowering, changes in flower morphology, and increased node number and internode length, resulting in plants that were taller than the wild type. Conversely, overexpression of GmAP1a resulted in early flowering and reduced plant height compared to the wild type under SD conditions. The gmap1 mutant and the overexpression lines also exhibited altered expression of several genes related to flowering and gibberellic acid metabolism, thereby providing insight into the role of GmAP1 in the regulatory networks controlling flowering time and plant height in soybean. Increased node number is the trait with the most promise for enhancing soybean pod number and grain yield. Therefore, the mutant alleles of the four AP1 homologs described here will be invaluable for molecular breeding of improved soybean yield.

SINGLE FLOWER TRUSS and SELF-PRUNING signal developmental and metabolic networks to guide cotton architectures

Patterns of indeterminate and determinate growth specify plant architecture and influence crop productivity. In cotton (Gossypium hirsutum), SINGLE FLOWER TRUSS (SFT) stimulates the transition to flowering and determinate growth, while its closely related antagonist SELF-PRUNING (SP) maintains meristems in indeterminate states to favor vegetative growth. Overexpressing GhSFT while simultaneously silencing GhSP produces highly determinate cotton with reduced foliage and synchronous fruiting. These findings suggest that GhSFT, GhSP, and genes in these signaling networks hold promise for enhancing 'annualized' growth patterns and improving cotton productivity and management. To identify the molecular programs underlying cotton growth habits, we used comparative co-expression networks, differential gene expression, and phenotypic analyses in cotton varieties expressing altered levels of GhSFT or GhSP. Using multiple cotton and tomato datasets, we identified diverse genetic modules highly correlated with SFT or SP orthologs which shared related Gene Ontologies in different crop species. Notably, altering GhSFT or GhSP levels in cotton affected the expression of genes regulating meristem fate and metabolic pathways. Further phenotypic analyses of gene products involved in photosynthesis, secondary metabolism, and cell wall biosynthesis showed that early changes in GhSFT and GhSP levels profoundly impacted later development in distal tissues. Identifying the molecular underpinnings of GhSFT and GhSP activities emphasizes their broad actions in regulating cotton architecture.

Prospects for the accelerated improvement of the resilient crop quinoa

Crops tolerant to drought and salt stress may be developed by two approaches. First, major crops may be improved by introducing genes from tolerant plants. For example, many major crops have wild relatives that are more tolerant to drought and high salinity than the cultivated crops, and, once deciphered, the underlying resilience mechanisms could be genetically manipulated to produce crops with improved tolerance. Secondly, some minor (orphan) crops cultivated in marginal areas are already drought and salt tolerant. Improving the agronomic performance of these crops may be an effective way to increase crop and food diversity, and an alternative to engineering tolerance in major crops. Quinoa (Chenopodium quinoa Willd.), a nutritious minor crop that tolerates drought and salinity better than most other crops, is an ideal candidate for both of these approaches. Although quinoa has yet to reach its potential as a fully domesticated crop, breeding efforts to improve the plant have been limited. Molecular and genetic techniques combined with traditional breeding are likely to change this picture. Here we analyse protein-coding sequences in the quinoa genome that are orthologous to domestication genes in established crops. Mutating only a limited number of such genes by targeted mutagenesis appears to be a promising route for accelerating the improvement of quinoa and generating a nutritious high-yielding crop that can meet the future demand for food production in a changing climate.

New Horizons for Dissecting Epistasis in Crop Quantitative Trait Variation

Uncovering the genes, variants, and interactions underlying crop diversity is a frontier in plant genetics. Phenotypic variation often does not reflect the cumulative effect of individual gene mutations. This deviation is due to epistasis, in which interactions between alleles are often unpredictable and quantitative in effect. Recent advances in genomics and genome-editing technologies are elevating the study of epistasis in crops. Using the traits and developmental pathways that were major targets in domestication and breeding, we highlight how epistasis is central in guiding the behavior of the genetic variation that shapes quantitative trait variation. We outline new strategies that illuminate how quantitative epistasis from modified gene dosage defines background dependencies. Advancing our understanding of epistasis in crops can reveal new principles and approaches to engineering targeted improvements in agriculture.

Combining focused identification of germplasm and core collection strategies to identify genebank accessions for central European soybean breeding

Environmental adaptation of crops is essential for reliable agricultural production and an important breeding objective. Genebanks provide genetic variation for the improvement of modern varieties, but the selection of suitable germplasm is frequently impeded by incomplete phenotypic data. We address this bottleneck by combining a Focused Identification of Germplasm Strategy (FIGS) with core collection methodology to select soybean (Glycine max) germplasm for Central European breeding from a collection of >17,000 accessions. By focussing on adaptation to high-latitude cold regions, we selected an "environmental precore" of 3,663 accessions using environmental data and compared the Donor opulation of Environments (DPE) in Asia and the Target Population of Environments (TPE) in Central Europe in the present and 2070. Using single nucleotide polymorphisms, we reduced the precore into two diverse core collections of 183 and 366 accessions to serve as diversity panels for evaluation in the TPE. Genetic differentiation between precore and non-precore accessions revealed genomic regions that control maturity, and novel candidate loci for environmental adaptation, demonstrating the potential of diversity panels for studying adaptation. Objective-driven core collections have the potential to increase germplasm utilization for abiotic adaptation by breeding for a rapidly changing climate, or de novo adaptation of crops to expand cultivation ranges.

Genetic dynamics of earlier maturity group emergence in south-to-north extension of Northeast China soybeans

KEY MESSAGE

This population genetic study is characterized with direct comparisons of days to flowering QTL-allele matrices between newly evolved and originally old maturity groups of soybeans to explore its evolutionary dynamics using the RTM-GWAS procedure. The Northeast China (NEC) soybeans are the major germplasm source of modern soybean production in Americas (> 80% of the world total). NEC is a relatively new soybean area in China, expanded after its nomadic status in the seventeenth century. At nine sites of four ecoregions in NEC, 361 varieties were tested for their days to flowering (DTF), a geography-sensitive trait as an indicator for maturity groups (MGs). The DTF reduced obviously along with soybeans extended to higher latitudes, ranging in 41-83 days and MG 000-III. Using the RTM-GWAS (restricted two-stage multi-locus model genome-wide association study) procedure, 81 QTLs with 342 alleles were identified, accounting for 77.85% genetic contribution (R2 = 0.01-7.74%/locus), and other 20.75% (98.60-77.85%, h2 = 98.60%) genetic variation was due to a collective of unmapped QTLs. With soybeans northward, breeding effort made the original MG I-III evolved to MG 0-00-000. In direct comparisons of QTL-allele matrices among MGs, the genetic dynamics are identified with local exotic introduction/migration (58.48%) as the first and selection against/exclusion of positive alleles causing new recombination (40.64%) as the second, while only a few allele emergence/mutation happened (0.88%, limited in MG 0, not in MG 00-000). In new MG emergence, 24 QTLs with 19 candidate genes are the major sources. A genetic potential of further DTF shortening (13-21 days) is predicted for NEC population. The QTL detection in individual ecoregions showed various ecoregion-specific QTLs-alleles/genes after co-localization treatment (removing the random environment shifting ones).

Knock-out of TERMINAL FLOWER 1 genes altered flowering time and plant architecture in Brassica napus

Background

TERMINAL FLOWER 1 (TFL1) is a member of phosphatidylethanolamine-binding protein (PEBP) family, which plays an important role in the determination of floral meristem identity and regulates flowering time in higher plants.

Results

Five BnaTFL1 gene copies were identified in the genome of Brassica napus. The phylogenetic analysis indicated that all five BnaTFL1 gene copies were clustered with their corresponding homologous copies in the ancestral species, B. rapa and B. oleracea. The expression of the BnaTFL1s were confined to flower buds, flowers, seeds, siliques and stem tissues and displayed distinct expression profiles. Knockout mutants of BnaC03.TFL1 generated by CRISPR/Cas9 exhibited early flowering phenotype, while the knockout mutants of the other gene copies had similar flowering time as the wild type. Furthermore, knock-out mutants of individual BnaTFL1 gene copy displayed altered plant architecture. The plant height, branch initiation height, branch number, silique number, number of seeds per silique and number of siliques on the main inflorescence were significantly reduced in the BnaTFL1 mutants.

Conclusions

Our results indicated that BnaC03.TFL1 negatively regulates flowering time in B. napus. BnaC03.TFL1 together with the other BnaTFL1 paralogues are essential for controlling the plant architecture.

Toward a "Green Revolution" for Soybean

Soybean (Glycine max), as an economically important food and oilseedcrop, is a major source of plant proteins and oils. Although considerable progress has been made in increasing the yields of rice, wheat, and maize through the "Green Revolution", little improvements have been made for soybean. With the increasing demand of soybean production and the rapid development of crop breeding technologies, time has come for this important crop to undergo a Green Revolution. Here, we briefly summarize the history of crop breeding and Green Revolution in other crops. We then discuss the possible directions and potential approaches toward achieving a Green Revolution for soybean. We provide our views and perspectives on how to breed new soybean varieties with improved yield.

Impacts of genomic research on soybean improvement in East Asia

It has been commonly accepted that soybean domestication originated in East Asia. Although East Asia has the historical merit in soybean production, the USA has become the top soybean producer in the world since 1950s. Following that, Brazil and Argentina have been the major soybean producers since 1970s and 1990s, respectively. China has once been the exporter of soybean to Japan before 1990s, yet she became a net soybean importer as Japan and the Republic of Korea do. Furthermore, the soybean yield per unit area in East Asia has stagnated during the past decade. To improve soybean production and enhance food security in these East Asian countries, much investment has been made, especially in the breeding of better performing soybean germplasms. As a result, China, Japan, and the Republic of Korea have become three important centers for soybean genomic research. With new technologies, the rate and precision of the identification of important genomic loci associated with desired traits from germplasm collections or mutants have increased significantly. Genome editing on soybean is also becoming more established. The year 2019 marked a new era for crop genome editing in the commercialization of the first genome-edited plant product, which is a high-oleic-acid soybean oil. In this review, we have summarized the latest developments in soybean breeding technologies and the remarkable progress in soybean breeding-related research in China, Japan, and the Republic of Korea.

Identification of QTL for perenniality and floral scent in cowpea (Vigna unguiculata [L.] Walp.)

Perennial habit and floral scent are major traits that distinguish domesticated cowpeas from their wild relatives. However, the genetic basis of these two important traits remains largely unknown in cowpea. Plant longevity, a perenniality-related trait, and floral scent, an outcrossing trait, were investigated using a RIL population derived from a cross between a domesticated and a wild cowpea. QTL analysis revealed three significant loci, one on chromosome 8 associated with plant longevity and two, on chromosomes 1 and 11, for floral scent. Genes within the QTL regions were identified. Genes encoding an F-box protein (Vigun08g215300) and two kinases (Vigun08g217000, Vigun08g217800), and involved in physiological processes including regulation of flowering time and plant longevity, were identified within the perenniality QTL region. A cluster of O-methyltransferase genes (Vigun11g096800, Vigun11g096900, Vigun11g097000, Vigun11g097600, and Vigun11g097800) was identified within the floral scent QTL region. These O-methyltransferase cowpea genes are orthologs of the Arabidopsis N-acetylserotonin O-methyltransferase (ASMT) gene, which is involved in the biosynthesis of melatonin. Melatonin is an indole derivative, which is an essential molecule for plant interactions with pollinators. These findings lay the foundation for further exploration of the genetic mechanisms of perenniality and floral scent in cowpea. Knowledge from this study can help in the development of new extended-growth cycle lines with increased yield or lines with increased outcrossing for population breeding.

Application of cucumber mosaic virus to efficient induction and long-term maintenance of virus-induced gene silencing in spinach

Virus-induced gene silencing (VIGS) is a useful tool for functional genomics in plants. In this study, we tried to apply cucumber mosaic virus (CMV) to efficient induction of VIGS in spinach. Although VIGS for spinach had been previously developed based on two viruses (beet curly top virus and tobacco rattle virus), they still have some problems with systemic movement and long-term maintenance of VIGS in spinach. Although ordinary CMV strains infect spinach inducing distinct mosaic symptoms, using a CMV pseudorecombinant, we can modify the viral pathogenicity to attenuate viral symptoms that may mask the silencing phenotype. We here successfully demonstrated the viral ability to silence the phytoene desaturase (PDS) and the dihydroflavonol 4-reductase (DFR) genes in spinach. Because CMV could quickly induce VIGS even at 7-10 days postinoculation and the virus did not disappear even at the flowering stage, this CMV-based VIGS system would contribute to functional genomics in spinach and especially to the elucidation of molecular mechanisms for some properties unique to spinach such as plasticity of sex expression; the CMV-induced VIGS can last until the flowering stage after the virus was inoculated onto the seedling.

The effects and interaction of soybean maturity gene alleles controlling flowering time, maturity, and adaptation in tropical environments

BACKGROUND

Soybean is native to the temperate zones of East Asia. Poor yields of soybean in West African countries may be partially attributed to inadequate adaptation of soybean to tropical environments. Adaptation will require knowledge of the effects of allelic combinations of major maturity genes (E1, E2, and E3) and stem architecture. The long juvenile trait (J) influences soybean flowering time in short, ~ 12 h days, which characterize tropical latitudes. Soybean plant architecture includes determinate or indeterminate stem phenotypes controlled by the Dt1 gene. Understanding the influence of these genetic components on plant development and adaptation is key to optimize phenology and improve soybean yield potential in tropical environments.

RESULTS

Soybean lines from five recombinant inbred populations were developed that varied in their combinations of targeted genes. The soybean lines were field tested in multiple environments and characterized for days to flowering (DTF), days to maturity (DTM), and plant height in locations throughout northern Ghana, and allelic combinations were determined for each line for associating genotype with phenotype. The results revealed significant differences based on genotype for DTF and DTM and allowed the comparison of different variant alleles of those genes. The mutant alleles of J and E1 had significant impact on DTF and DTM, and alleles of those genes interacted with each other for DTF but not DTM. The Dt1 gene significantly influenced plant height but not DTF or DTM.

CONCLUSIONS

This research identified major and minor effect alleles of soybean genes that can be combined to control DTF, DTM, and plant height in short day tropical environments in Ghana. These phenotypes contribute to adaptation to a low latitude environment that can be optimized in a soybean breeding program with targeted selection of desired allele combinations. The knowledge of the genetic control of these traits will enhance molecular breeding to produce optimally adapted soybean varieties targeted to tropical environments.

Genome-wide association and epistatic interactions of flowering time in soybean cultivar

Genome-wide association studies (GWAS) have enabled the discovery of candidate markers that play significant roles in various complex traits in plants. Recently, with increased interest in the search for candidate markers, studies on epistatic interactions between single nucleotide polymorphism (SNP) markers have also increased, thus enabling the identification of more candidate markers along with GWAS on single-variant-additive-effect. Here, we focused on the identification of candidate markers associated with flowering time in soybean (Glycine max). A large population of 2,662 cultivated soybean accessions was genotyped using the 180k Axiom® SoyaSNP array, and the genomic architecture of these accessions was investigated to confirm the population structure. Then, GWAS was conducted to evaluate the association between SNP markers and flowering time. A total of 93 significant SNP markers were detected within 59 significant genes, including E1 and E3, which are the main determinants of flowering time. Based on the GWAS results, multilocus epistatic interactions were examined between the significant and non-significant SNP markers. Two significant and 16 non-significant SNP markers were discovered as candidate markers affecting flowering time via interactions with each other. These 18 candidate SNP markers mapped to 18 candidate genes including E1 and E3, and the 18 candidate genes were involved in six major flowering pathways. Although further biological validation is needed, our results provide additional information on the existing flowering time markers and present another option to marker-assisted breeding programs for regulating flowering time of soybean.

Molecular and transcriptional characterization of phosphatidyl ethanolamine-binding proteins in wild peanuts Arachis duranensis and Arachis ipaensis

BACKGROUND

Phosphatidyl ethanolamine-binding proteins (PEBPs) are involved in the regulation of plant architecture and flowering time. The functions of PEBP genes have been studied in many plant species. However, little is known about the characteristics and expression profiles of PEBP genes in wild peanut species, Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanuts.

RESULTS

In this study, genome-wide identification methods were used to identify and characterize a total of 32 peanut PEBP genes, 16 from each of the two wild peanut species, A. duranensis and A. ipaensis. These PEBP genes were classified into 3 groups (TERMINAL FLOWER1-like, FLOWERING LOCUS T-like, and MOTHER OF FT AND TFL1-like) based on their phylogenetic relationships. The gene structures, motifs, and chromosomal locations for each of these PEBPs were analyzed. In addition, 4 interchromosomal duplications and 1 tandem duplication were identified in A. duranensis, and 2 interchromosomal paralogs and 1 tandem paralog were identified in A. ipaensis. Ninety-five different cis-acting elements were identified in the PEBP gene promoter regions and most genes had different numbers and types of cis-elements. As a result, the transcription patterns of these PEBP genes varied in different tissues and under long day and short day conditions during different growth phases, indicating the functional diversities of PEBPs in different tissues and their potential functions in plant photoperiod dependent developmental pathways. Moreover, our analysis revealed that AraduF950M/AraduWY2NX in A. duranensis, and Araip344D4/Araip4V81G in A. ipaensis are good candidates for regulating plant architecture, and that Aradu80YRY, AraduYY72S, and AraduEHZ9Y in A. duranensis and AraipVEP8T in A. ipaensis may be key factors regulating flowering time.

CONCLUSION

Sixteen PEBP genes were identified and characterized from each of the two diploid wild peanut genomes, A. duranensis and A. ipaensis. Genetic characterization and spatio-temporal expression analysis support their importance in plant growth and development. These findings further our understanding of PEBP gene functions in plant species.

A Post-domestication Mutation, Dt2, Triggers Systemic Modification of Divergent and Convergent Pathways Modulating Multiple Agronomic Traits in Soybean

The semi-determinate stem growth habit in leguminous crops, similar to the "green revolution" semi-dwarf trait in cereals, is a key plant architecture trait that affects several other traits determining grain yield. In soybean semi-determinacy is modulated by a post-domestication gain-of-function mutation in the gene, Dt2, which encodes an MADS-box transcription factor. However, its role in systemic modification of stem growth and other traits is unknown. In this study, we show that Dt2 functions not only as a direct repressor of Dt1, which prevents terminal flowering, but also as a direct activator of putative floral integrator/identity genes including GmSOC1, GmAP1, and GmFUL, which likely promote flowering. We also demonstrate that Dt2 functions as a direct repressor of the putative drought-responsive transcription factor gene GmDREB1D, and as a direct activator of GmSPCH and GmGRP7, which are potentially associated with asymmetric division of young epidermal cells and stomatal opening, respectively, and may affect the plant's water-use efficiency (WUE). Intriguingly, Dt2 was found to be a direct activator or repressor of the precursors of eight microRNAs targeting genes potentially associated with meristem maintenance, flowering time, stomatal density, WUE, and/or stress responses. This study thus reveals the molecular basis of pleiotropy associated with plant productivity, adaptability, and environmental resilience.

Indeterminate growth of the umbel inflorescence and bulb is associated with increased expression of the TFL1 homologue, AcTFL1, in onion

TERMINAL FLOWER1 (TFL1) is a key gene for maintenance of vegetative and inflorescence indeterminacy and architecture. In onion, flowering and bulbing are two distinct developmental phases, each under complex environmental regulatory control. We have identified two CEN/TFL1-like genes from onion designated as AcTFL1 and AcCEN1. AcTFL1 is expressed during bulbing and inflorescence development with expression increasing with indeterminate growth of the umbel and the bulb suggesting possible conservation of function. Increase in AcTFL1 expression during umbel growth is associated with a simultaneous reduction in expression of AcLFY. Expression of AcTFL1 within the bulb is lowest in the outermost layers and highest in the innermost (youngest) layers. Bulb storage at room temperature or in cold leads to a gradual reduction in AcTFL1 levels in the meristem-containing tissues, the decrease being faster in the variety not requiring vernalization. Constitutive expression of AcTFL1, but not AcCEN1 complements the Arabidopsis tfl1-14 mutant and delays flowering in wild type suggesting conservation of the AcTFL1 function even in the distantly related Arabidopsis. Taken together, AcTFL1 appears to be the functional counterpart of TFL1 and regulates indeterminate growth of the umbel inflorescence as well as bulb development in onion.