Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean

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.

Pathways to de novo domestication of crop wild relatives

Growing knowledge about crop domestication, combined with increasingly powerful gene-editing toolkits, sets the stage for the continual domestication of crop wild relatives and other lesser-known plant species.

FRUITFULL-like genes regulate flowering time and inflorescence architecture in tomato

The timing of flowering and the inflorescence architecture are critical for the reproductive success of tomato (Solanum lycopersicum), but the gene regulatory networks underlying these traits have not been fully explored. Here, we show that the tomato FRUITFULL-like (FUL-like) genes FUL2 and MADS-BOX PROTEIN 20 (MBP20) promote the vegetative-to-reproductive transition and repress inflorescence branching by inducing floral meristem (FM) maturation. FUL1 fulfils a less prominent role and appears to depend on FUL2 and MBP20 for its upregulation in the inflorescence- and floral meristems. MBP10, the fourth tomato FUL-like gene, has probably lost its function. The tomato FUL-like proteins cannot homodimerize in in vitro assays, but heterodimerize with various other MADS-domain proteins, potentially forming distinct complexes in the transition meristem and FM. Transcriptome analysis of the primary shoot meristems revealed various interesting downstream targets, including four repressors of cytokinin signaling that are upregulated during the floral transition in ful1 ful2 mbp10 mbp20 mutants. FUL2 and MBP20 can also bind in vitro to the upstream regions of these genes, thereby probably directly stimulating cell division in the meristem upon the transition to flowering. The control of inflorescence branching does not occur via the cytokinin oxidase/dehydrogenases (CKXs) but may be regulated by repression of transcription factors such as TOMATO MADS-box gene 3 (TM3) and APETALA 2b (AP2b).

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

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.

Developing functional markers for vitamin E biosynthesis in oil palm

Vitamin E is essential for human health and plays positive roles in anti-oxidation. Previously, we detected large variation in vitamin E content among 161 oil palm accessions. In this study, twenty oil palm accessions with distinct variation in vitamin E contents (171.30 to 1 258.50 ppm) were selected for genetic variation analysis and developing functional markers associated with vitamin E contents. Thirty-seven homologous genes in oil palm belonging to vitamin E biosynthesis pathway were identified via BLASTP analysis, the lengths of which ranged from 426 to 25 717 bp (average 7 089 bp). Multiplex PCR sequencing for the 37 genes found 1 703 SNPs and 85 indels among the 20 oil palm accessions, with 226 SNPs locating in the coding regions. Clustering analysis for these polymorphic loci showed that the 20 oil palm accessions could be divided into five groups. Among these groups, group I included eight oil palm accessions whose vitamin E content (mean value: 893.50 ppm) was far higher than other groups (mean value 256.29 to 532.94 ppm). Correlation analysis between the markers and vitamin E traits showed that 134 SNP and 7 indel markers were significantly (p < 0.05) related with total vitamin E content. Among these functional markers, the indel EgTMT-1-24 was highly correlated with variation in vitamin E content, especially tocotrienol content. Our study identified a number of candidate function associated markers and provided clues for further research into molecular breeding for high vitamin E content oil palm.

Soybean Yield Formation Physiology - A Foundation for Precision Breeding Based Improvement

The continued improvement of crop yield is a fundamental driver in agriculture and is the goal of both plant breeders and researchers. Plant breeders have been remarkably successful in improving crop yield, as demonstrated by the continued release of varieties with improved yield potential. This has largely been accomplished through performance-based selection, without specific knowledge of the molecular mechanisms underpinning these improvements. Insight into molecular mechanisms has been provided by plant molecular, genetic, and biochemical research through elucidation of the function of genes and pathways that underlie many of the physiological processes that contribute to yield potential. Despite this knowledge, the impact of most genes and pathways on yield components have not been tested in key crops or in a field environment for yield assessment. This gap is difficult to bridge, but field-based physiological knowledge offers a starting point for leveraging molecular targets to successfully apply precision breeding technologies such as genome editing. A better understanding of both the molecular mechanisms underlying crop yield physiology and yield limiting processes under field conditions is essential for elucidating which combinations of favorable alleles are required for yield improvement. Consequently, one goal in plant biology should be to more fully integrate crop physiology, breeding, genetics, and molecular knowledge to identify impactful precision breeding targets for relevant yield traits. The foundation for this is an understanding of yield formation physiology. Here, using soybean as an example, we provide a top-down review of yield physiology, starting with the fact that yield is derived from a population of plants growing together in a community. We review yield and yield-related components to provide a basic overview of yield physiology, synthesizing these concepts to highlight how such knowledge can be leveraged for soybean improvement. Using genome editing as an example, we discuss why multiple disciplines must be brought together to fully realize the promise of precision breeding-based crop improvement.

A giant NLR gene confers broad-spectrum resistance to Phytophthora sojae in soybean

Phytophthora root and stem rot caused by P. sojae is a destructive soybean soil-borne disease found worldwide. Discovery of genes conferring broad-spectrum resistance to the pathogen is a need to prevent the outbreak of the disease. Here, we show that soybean Rps11 is a 27.7-kb nucleotide-binding site-leucine-rich repeat (NBS-LRR or NLR) gene conferring broad-spectrum resistance to the pathogen. Rps11 is located in a genomic region harboring a cluster of large NLR genes of a single origin in soybean, and is derived from rounds of unequal recombination. Such events result in promoter fusion and LRR expansion that may contribute to the broad resistance spectrum. The NLR gene cluster exhibits drastic structural diversification among phylogenetically representative varieties, including gene copy number variation ranging from five to 23 copies, and absence of allelic copies of Rps11 in any of the non-Rps11-donor varieties examined, exemplifying innovative evolution of NLR genes and NLR gene clusters.

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.

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.

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.

Multi-omic dissection of the drought resistance traits of soybean landrace LX

With diverse genetic backgrounds, soybean landraces are valuable resource for breeding programs. Herein, we apply multi-omic approaches to extensively characterize the molecular basis of drought tolerance in the soybean landrace LX. Initial screens established that LX performed better with PEG6000 treatment than control cultivars. LX germinated better than William 82 under drought conditions and accumulated more anthocyanin and flavonoids. Untargeted mass spectrometry in combination with transcriptomic analyses revealed the chemical diversity and genetic basis underlying the overall performance of LX landrace. Under control and drought conditions, significant differences in the expression of a suite of secondary metabolism genes, particularly those involved in the general phenylpropanoid pathway and flavonoid but not lignin biosynthesis, were seen in LX and William 82. The expression of these genes correlated with the corresponding metabolites in LX plants. Further correlation analysis between metabolites and transcripts identified pathway structural genes and transcription factors likely are responsible for the LX agronomic traits. The activities of some key biosynthetic genes or regulators were confirmed through heterologous expression in transgenic Arabidopsis and hairy root transformation in soybean. We propose a regulatory mechanism based on flavonoid secondary metabolism and adaptive traits of this landrace which could be of relevance to cultivated soybean.

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.

GmCCD4 controls carotenoid content in soybeans

To better understand the mechanisms regulating plant carotenoid metabolism in staple crop, we report the map-based cloning and functional characterization of the Glycine max carotenoid cleavage dioxygenase 4 (GmCCD4) gene, which encodes a carotenoid cleavage dioxygenase enzyme involved in metabolizing carotenoids into volatile β-ionone. Loss of GmCCD4 protein function in four Glycine max increased carotenoid content (gmicc) mutants resulted in yellow flowers due to excessive accumulation of carotenoids in flower petals. The carotenoid contents also increase three times in gmicc1 seeds. A genome-wide association study indicated that the GmCCD4 locus was one major locus associated with carotenoid content in natural population. Further analysis indicated that the haplotype-1 of GmCCD4 gene was positively associated with higher carotenoid levels in soybean cultivars and accumulated more β-carotene in engineered E. coli with ectopic expression of different GmCCD4 haplotypes. These observations uncovered that GmCCD4 was a negative regulator of carotenoid content in soybean, and its various haplotypes provide useful resources for future soybean breeding practice.

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.

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.

VrLELP controls flowering time under short-day conditions in Arabidopsis

Flowering time has a critically important effect on the reproduction of plants, and many components involved in flowering-time regulation have been identified in multiple plant species. However, studies of the flowering-time genes in mungbean (Vigna radiata) have been limited. Here, we characterized a novel mungbean gene, VrLELP, involved in flowering-time regulation in transgenic Arabidopsis. Subcellular localization analysis revealed that VrLELP was localized in the membrane, cytoplasm and nucleus and the nucleus and membrane contained higher signal than cytoplasm, similar to the empty vector control. The expression of VrLELP was higher in leaves and pods and lower in nodule roots relative to other tissues. The expression of VrLELP varied during flower development. The expression of VrLELP also varied during the day, reaching a peak after 12 h of illumination under long-day conditions. In contrast, under short-day conditions, the abundance of VrLELP transcripts changed little throughout the day. In addition, VrLELP delayed flowering time in transgenic Arabidopsis plants by suppressing the expression of the flowering-time genes CO and FT under short-day conditions. However, VrLELP did not affect flowering time under long-day conditions in Arabidopsis. Our study provides essential information for future studies of the molecular mechanisms of the flowering-time regulation system in mungbean.

Characterization of Mungbean CONSTANS-LIKE Genes and Functional Analysis of CONSTANS-LIKE 2 in the Regulation of Flowering Time in Arabidopsis

CONSTANS-LIKE (COL) genes play important roles in the regulation of plant growth and development, and they have been analyzed in many plant species. However, few studies have examined COL genes in mungbean (Vigna radiata). In this study, we identified and characterized 31 mungbean genes whose proteins contained B-Box domains. Fourteen were designated as VrCOL genes and were distributed on 7 of the 11 mungbean chromosomes. Based on their phylogenetic relationships, VrCOLs were clustered into three groups (I, II, and III), which contained 4, 6, and 4 members, respectively. The gene structures and conserved motifs of the VrCOL genes were analyzed, and two duplicated gene pairs, VrCOL1/VrCOL2 and VrCOL8/VrCOL9, were identified. A total of 82 cis-acting elements were found in the VrCOL promoter regions, and the numbers and types of cis-acting elements in each VrCOL promoter region differed. As a result, the expression patterns of VrCOLs varied in different tissues and throughout the day under long-day and short-day conditions. Among these VrCOL genes, VrCOL2 showed a close phylogenetic relationship with Arabidopsis thaliana CO and displayed daily oscillations in expression under short-day conditions but not long-day conditions. In addition, overexpression of VrCOL2 accelerated flowering in Arabidopsis under short-day conditions by affecting the expression of the flowering time genes AtFT and AtTSF. Our study lays the foundation for further investigation of VrCOL gene functions.

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.

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.

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.

Transcriptomic and biochemical analyses of the accumulation of sucrose in mungbean (Vigna radiata (L.) Wilczek) leaves after pod removal

Synchronous pod maturity is critical for increasing grain yield. The candidate genes involved in synchronous pod maturity were identified through RNA-seq and HPLC. Mungbean (Vigna radiata [L.] Wilczek), an important source of carbohydrate and protein in Asia, is characterized by nonsynchronous pod maturity; consequently, harvesting is labor intensive. Because pod maturity is associated with synthesis and remobilization of sucrose, we examined changes in sucrose levels and transcriptome in leaf (source) tissues after pod (sink) removal using two genotypes, VC1973A and V2984; VC1973A had higher synchronicity in pod maturity than V2984. After pod removal, much higher number of pods were produced in V2984 than VC1973A. The sucrose content of leaf tissues significantly decreased in V2984 because it continued to utilize assimilates from leaves for producing new pods, but significantly increased in VC1973A because of the loss of sink. Transcriptome analysis revealed that the number of differentially expressed genes was approximately fourfold higher in VC1973A than in those of V2984 after pod removal. The expression of two paralogous genes (Vradi01g05010 and Vradi10g08240), encoding beta-glucosidase enzymes, significantly decreased in VC1973A after pod removal and was significantly lower in depodded VC1973A than depodded V2984, indicating these two genes may participate in sucrose utilization for seed development by regulating the level of glucose. The results of this study will help elucidate the genetic basis of synchronous pod maturity in mungbean and facilitate the development of new cultivars with synchronous pod maturity.

Identification of Regulatory SNPs Associated with Vicine and Convicine Content of Vicia faba Based on Genotyping by Sequencing Data Using Deep Learning

Faba bean (Vicia faba) is a grain legume, which is globally grown for both human consumption as well as feed for livestock. Despite its agro-ecological importance the usage of Vicia faba is severely hampered by its anti-nutritive seed-compounds vicine and convicine (V+C). The genes responsible for a low V+C content have not yet been identified. In this study, we aim to computationally identify regulatory SNPs (rSNPs), i.e., SNPs in promoter regions of genes that are deemed to govern the V+C content of Vicia faba. For this purpose we first trained a deep learning model with the gene annotations of seven related species of the Leguminosae family. Applying our model, we predicted putative promoters in a partial genome of Vicia faba that we assembled from genotyping-by-sequencing (GBS) data. Exploiting the synteny between Medicago truncatula and Vicia faba, we identified two rSNPs which are statistically significantly associated with V+C content. In particular, the allele substitutions regarding these rSNPs result in dramatic changes of the binding sites of the transcription factors (TFs) MYB4, MYB61, and SQUA. The knowledge about TFs and their rSNPs may enhance our understanding of the regulatory programs controlling V+C content of Vicia faba and could provide new hypotheses for future breeding programs.

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).

Closing the Phenotyping Gap: High Resolution UAV Time Series for Soybean Growth Analysis Provides Objective Data from Field Trials

Close remote sensing approaches can be used for high throughput on-field phenotyping in the context of plant breeding and biological research. Data on canopy cover (CC) and canopy height (CH) and their temporal changes throughout the growing season can yield information about crop growth and performance. In the present study, sigmoid models were fitted to multi-temporal CC and CH data obtained using RGB imagery captured with a drone for a broad set of soybean genotypes. The Gompertz and Beta functions were used to fit CC and CH data, respectively. Overall, 90.4% fits for CC and 99.4% fits for CH reached an adjusted R2 > 0.70, demonstrating good performance of the models chosen. Using these growth curves, parameters including maximum absolute growth rate, early vigor, maximum height, and senescence were calculated for a collection of soybean genotypes. This information was also used to estimate seed yield and maturity (R8 stage) (adjusted R2 = 0.51 and 0.82). Combinations of parameter values were tested to identify genotypes with interesting traits. An integrative approach of fitting a curve to a multi-temporal dataset resulted in biologically interpretable parameters that were informative for relevant traits.

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.

GsMAS1 Encoding a MADS-box Transcription Factor Enhances the Tolerance to Aluminum Stress in Arabidopsis thaliana

The MADS-box transcription factors (TFs) are essential in regulating plant growth and development, and conferring abiotic and metal stress resistance. This study aims to investigate GsMAS1 function in conferring tolerance to aluminum stress in Arabidopsis. The GsMAS1 from the wild soybean BW69 line encodes a MADS-box transcription factor in Glycine soja by bioinformatics analysis. The putative GsMAS1 protein was localized in the nucleus. The GsMAS1 gene was rich in soybean roots presenting a constitutive expression pattern and induced by aluminum stress with a concentration-time specific pattern. The analysis of phenotypic observation demonstrated that overexpression of GsMAS1 enhanced the tolerance of Arabidopsis plants to aluminum (Al) stress with larger values of relative root length and higher proline accumulation compared to those of wild type at the AlCl3 treatments. The genes and/or pathways regulated by GsMAS1 were further investigated under Al stress by qRT-PCR. The results indicated that six genes resistant to Al stress were upregulated, whereas AtALMT1 and STOP2 were significantly activated by Al stress and GsMAS1 overexpression. After treatment of 50 μM AlCl3, the RNA abundance of AtALMT1 and STOP2 went up to 17-fold and 37-fold than those in wild type, respectively. Whereas the RNA transcripts of AtALMT1 and STOP2 were much higher than those in wild type with over 82% and 67% of relative expression in GsMAS1 transgenic plants, respectively. In short, the results suggest that GsMAS1 may increase resistance to Al toxicity through certain pathways related to Al stress in Arabidopsis.

Quantitative Genomic Dissection of Soybean Yield Components

Soybean is a crop of major economic importance with low rates of genetic gains for grain yield compared to other field crops. A deeper understanding of the genetic architecture of yield components may enable better ways to tackle the breeding challenges. Key yield components include the total number of pods, nodes and the ratio pods per node. We evaluated the SoyNAM population, containing approximately 5600 lines from 40 biparental families that share a common parent, in 6 environments distributed across 3 years. The study indicates that the yield components under evaluation have low heritability, a reasonable amount of epistatic control, and partially oligogenic architecture: 18 quantitative trait loci were identified across the three yield components using multi-approach signal detection. Genetic correlation between yield and yield components was highly variable from family-to-family, ranging from -0.2 to 0.5. The genotype-by-environment correlation of yield components ranged from -0.1 to 0.4 within families. The number of pods can be utilized for indirect selection of yield. The selection of soybean for enhanced yield components can be successfully performed via genomic prediction, but the challenging data collections necessary to recalibrate models over time makes the introgression of QTL a potentially more feasible breeding strategy. The genomic prediction of yield components was relatively accurate across families, but less accurate predictions were obtained from within family predictions and predicting families not observed included in the calibration set.

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.

Deciphering the Genetic Architecture of Plant Height in Soybean Using Two RIL Populations Sharing a Common M8206 Parent

Plant height (PH) is an important agronomic trait that is closely related to soybean yield and quality. However, it is a complex quantitative trait governed by multiple genes and is influenced by environment. Unraveling the genetic mechanism involved in PH, and developing soybean cultivars with desirable PH is an imperative goal for soybean breeding. In this regard, the present study used high-density linkage maps of two related recombinant inbred line (RIL) populations viz., MT and ZM evaluated in three different environments to detect additive and epistatic effect quantitative trait loci (QTLs) as well as their interaction with environments for PH in Chinese summer planting soybean. A total of eight and 12 QTLs were detected by combining the composite interval mapping (CIM) and mixed-model based composite interval mapping (MCIM) methods in MT and ZM populations, respectively. Among these QTLs, nine QTLs viz., QPH-2, qPH-6-2_MT, QPH-6, qPH-9-1_ZM, qPH-10-1_ZM, qPH-13-1_ZM, qPH-16-1_MT, QPH-17 and QPH-19 were consistently identified in multiple environments or populations, hence were regarded as stable QTLs. Furthermore, Out of these QTLs, three QTLs viz., qPH-4-2_ZM, qPH-15-1_MT and QPH-17 were novel. In particular, QPH-17 could detect in both populations, which was also considered as a stable and major QTL in Chinese summer planting soybean. Moreover, eleven QTLs revealed significant additive effects in both populations, and out of them only six showed additive by environment interaction effects, and the environment-independent QTLs showed higher additive effects. Finally, six digenic epistatic QTLs pairs were identified and only four additive effect QTLs viz., qPH-6-2_MT, qPH-19-1_MT/QPH-19, qPH-5-1_ZM and qPH-17-1_ZM showed epistatic effects. These results indicate that environment and epistatic interaction effects have significant influence in determining genetic basis of PH in soybean. These results would not only increase our understanding of the genetic control of plant height in summer planting soybean but also provide support for implementing marker assisted selection (MAS) in developing cultivars with ideal plant height as well as gene cloning to elucidate the mechanisms of plant height.

Rhizobial tRNA-derived small RNAs are signal molecules regulating plant nodulation

Rhizobial infection and root nodule formation in legumes require recognition of signal molecules produced by the bacteria and their hosts. Here, we show that rhizobial transfer RNA (tRNA)-derived small RNA fragments (tRFs) are signal molecules that modulate host nodulation. Three families of rhizobial tRFs were confirmed to regulate host genes associated with nodule initiation and development through hijacking the host RNA-interference machinery that involves ARGONAUTE 1. Silencing individual tRFs with the use of short tandem target mimics or by overexpressing their targets represses root hair curling and nodule formation, whereas repressing these targets with artificial microRNAs identical to the respective tRFs or mutating these targets with CRISPR-Cas9 promotes nodulation. Our findings thus uncover a bacterial small RNA-mediated mechanism for prokaryote-eukaryote interaction and may pave the way for enhancing nodulation efficiency in legumes.

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

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

QTL architecture of vine growth habit and gibberellin oxidase gene diversity in wild soybean (Glycine soja)

Vine growth habit (VGH) is a beneficial phenotype in many wild plants, and is considered an important domesticated-related trait in soybean. However, its genetic basis remains largely unclear. Hence, in the present study we used an integrated strategy combining linkage mapping and population genome diversity analyses to reveal the genetics of VGH in soybean. In this regard, two recombinant inbred line (RIL) populations derived by crossing a common wild soybean genotype (PI342618B) with two cultivated lines viz., NN 86-4 and NN 493-1 were used to map quantitative trait loci (QTL) for VGH. Here, we identified seven and five QTLs at flowering stage (R1) and maturity stage (R8), respectively, and among them qVGH-18-1, qVGH-18-2, qVGH-19-3, qVGH-19-4 were identified as major loci (R2 > 10% and detection time ≥2). However, qVGH-18-2 was considered as a main QTL for VGH being consistently identified in both RIL populations as well as all growth stages and cropping years. Out of all the annotated genes within qVGH-18-2, Glyma18g06870 was identified as the candidate gene and named as VGH1, which was a gibberellin oxidase (GAox) belongs to 2-oxoglutarate-dependent dioxygenase (2- ODD). Interestingly, there was one member of 2-ODD/GAox in qVGH-18-1 and qVGH-19-4 named as VGH2 and VGH3, respectively. Moreover, from sequencing data analysis VGH1 and three other GAox genes were found significantly divergent between vine and erect soybean with FST value larger than 0.25. Hence, GAox was assumed to play a major role in governing inheritance of VGH in soybean. Therefore, elucidating the genetic mechanism of GAox is very useful for exploring VGH and other stem traits, as well as genetic improvement of plant type in soybean.

Turning Meristems into Fortresses

TERMINAL FLOWER1 (TFL1) was named from knockout Arabidopsis thaliana mutants in which the inflorescence abnormally terminates into a flower. In wild type plants, the expression of TFL1 in the center of the inflorescence meristem represses the flower meristem identity genes LEAFY (LFY) and APETALA1 (AP1) to maintain indeterminacy. LFY and AP1 are activated by flowering signals that antagonize TFL1. Its characterization in numerous species revealed that the TFL1-mediated regulation of meristem fate has broader impacts on plant development than originally depicted in A. thaliana. By blocking floral transition, TFL1 genes participate in the control of juvenility, shoot growth pattern, inflorescence architecture, and the establishment of life history strategies. Here, we contextualize the role of the TFL1-mediated protection of meristem indeterminacy throughout plant development.

CRISPR/Cas9-mediated targeted mutagenesis of GmSPL9 genes alters plant architecture in soybean

BACKGROUND

The plant architecture has significant effects on grain yield of various crops, including soybean (Glycine max), but the knowledge on optimization of plant architecture in order to increase yield potential is still limited. Recently, CRISPR/Cas9 system has revolutionized genome editing, and has been widely utilized to edit the genomes of a diverse range of crop plants.

RESULTS

In the present study, we employed the CRISPR/Cas9 system to mutate four genes encoding SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors of the SPL9 family in soybean. These four GmSPL9 genes are negatively regulated by GmmiR156b, a target for the improvement of soybean plant architecture and yields. The soybean Williams 82 was transformed with the binary CRISPR/Cas9 plasmid, assembled with four sgRNA expression cassettes driven by the Arabidopsis thaliana U3 or U6 promoter, targeting different sites of these four SPL9 genes via Agrobacterium tumefaciens-mediated transformation. A 1-bp deletion was detected in one target site of the GmSPL9a and one target site of the GmSPL9b, respectively, by DNA sequencing analysis of two T0-generation plants. T2-generation spl9a and spl9b homozygous single mutants exhibited no obvious phenotype changes; but the T2 double homozygous mutant spl9a/spl9b possessed shorter plastochron length. In T4 generation, higher-order mutant plants carrying various combinations of mutations showed increased node number on the main stem and branch number, consequently increased total node number per plants at different levels. In addition, the expression levels of the examined GmSPL9 genes were higher in the spl9b-1 single mutant than wild-type plants, which might suggest a feedback regulation on the expression of the investigated GmSPL9 genes in soybean.

CONCLUSIONS

Our results showed that CRISPR/Cas9-mediated targeted mutagenesis of four GmSPL9 genes in different combinations altered plant architecture in soybean. The findings demonstrated that GmSPL9a, GmSPL9b, GmSPL9c and GmSPL9 function as redundant transcription factors in regulating plant architecture in soybean.

Effect of change from a determinate to a semi-determinate growth habit on the yield and lodging resistance of soybeans in the northeast region of Japan

Although an indeterminate growth habit is attractive to develop high-yield soybean varieties with higher number of pods (Glycine max (L). Merr.), lodging in indeterminate varieties remains a problem in Japan. As the semi-determinate varieties have shorter main stem length than the indeterminate varieties, this trait can be useful to improve varieties with high yield and low lodging risk. We introduced the genes Dt1 and Dt2, which regulate stem growth habit, into three determinate varieties by backcrossing and evaluated the resulting effects on yield and lodging tendency under four different growing environments. The yield and lodging degree of the semi-determinate and indeterminate lines were higher and more severe than those of the determinate lines. Despite the lower overall lodging score, the semi-determinate lines had marginally lower overall yield than that of the indeterminate lines. However, the effect of introduction of semi-determinate traits on yield and lodging degree was different in the three backgrounds, with the yield of semi-determinate lines being the highest and the difference in lodging degree between the semi-determinate and determinate lines being under 1.0 in one background. Therefore, semi-determinate growth habit has potential to develop high yielding varieties with low lodging risk.

GmPGL1, a Thiamine Thiazole Synthase, Is Required for the Biosynthesis of Thiamine in Soybean

Thiamine is an essential cofactor in several enzymatic reactions for all living organisms. Animals cannot synthesize thiamine and depend on their diet. Enhancing the content of thiamine is one of the most important goals of plant breeding to solve the thiamine deficiency associated with the low-thiamin staple crops. In this study, a Glycine max pale green leaf 1 (Gmpgl1) mutant was isolated from the EMS mutagenized population of soybean cultivar, Williams 82. Map-based cloning of the GmPGL1 locus revealed a single nucleotide deletion at the 292th nucleotide residue of the first exon of Glyma.10g251500 gene in Gmpgl1 mutant plant, encoding a thiamine thiazole synthase. Total thiamine contents decreased in both seedlings and seeds of the Gmpgl1 mutant. Exogenous application of thiazole restored the pale green leaf phenotype of the mutant. The deficiency of thiamine in Gmpgl1 mutant led to reduced activities of the pyruvate dehydrogenase (PDH) and pyruvate decarboxylase (PDC), and decreased contents of six amino acids as compared to that in the wild type plants. These results revealed that GmPGL1 played an essential role in thiamine thiazole biosynthesis.

Parallel domestication with a broad mutational spectrum of determinate stem growth habit in leguminous crops

Stem growth habit is a key plant architecture trait determining yield potential in grain legumes, and the phenotypic change from the indeterminate stem growth habit of wild mungbeans (Vigna radiata) to the determinate stem growth habit of cultivated mungbeans is a critical domestication transition. Here we show that indeterminate stem growth in wild mungbean is modulated by a single gene, VrDet1, which encodes a signaling protein of shoot apical meristems. The transition from an indeterminate to a determinate stem growth habit was achieved by selection of two linked point mutations in two putative cis-regulatory elements, resulting in a significant reduction in gene expression. Both the wild-type nucleotides corresponding to the two point mutations were essential for VrDet1 function. In addition, two highly diverse haplotypes of Vrdet1 were found in cultivated mungbeans, suggesting dual domestication of Vrdet1. VrDet1 was orthologous to Dt1 in wild soybean and PvTFL1y in wild common bean, where multiple loss-of-function mutations altering the coding sequences of individual genes were selected to produce determinate stems in cultivated accessions. Interspecific comparison of these orthologs in the wild and cultivated accessions reveals the most conservative interspecific and intraspecific parallel domestication events with the broadest mutational spectrum of a domestication trait in leguminous crops. We also found that interspecifically and functionally conserved promoters possess cis-regulatory elements that are highly conserved in kind but greatly variable in number and order, demonstrating the evolutionary dynamics of regulatory sequences. This work provides insights into the origins of cultivated mungbean and exemplifies the conservativeness and plasticity of the domestication processes of related crops.

Transcriptomic comparison reveals genetic variation potentially underlying seed developmental evolution of soybeans

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. However, the genetic variations underlying soybean domestication are not well known. Comparative transcriptomics revealed that a small portion of the orthologous genes might have been fast evolving. In contrast, three gene expression clusters were identified as divergent by their expression patterns, which occupied 37.44% of the total genes, hinting at an essential role for gene expression alteration in soybean domestication. Moreover, the most divergent stage in gene expression between wild and cultivated soybeans occurred during seed development around the cotyledon stage (15 d after fertilization, G15). A module in which the co-expressed genes were significantly down-regulated at G15 of wild soybeans was identified. The divergent clusters and modules included substantial differentially expressed genes (DEGs) between wild and cultivated soybeans related to cell division, storage compound accumulation, hormone response, and seed maturation processes. Chromosomal-linked DEGs, quantitative trait loci controlling seed weight and oil content, and selection sweeps revealed candidate DEGs at G15 in the fruit-related divergence of G. max and G. soja. Our work establishes a transcriptomic selection mechanism for altering gene expression during soybean domestication, thus shedding light on the molecular networks underlying soybean seed development and breeding strategy.

White seed color in common bean (Phaseolus vulgaris) results from convergent evolution in the P (pigment) gene

The presence of seed color in common bean (Phaseolus vulgaris) requires the dominant-acting P (pigment) gene, and white seed is a recessive phenotype in all domesticated races of the species. P was classically associated with seed size, thus describing it as the first genetic marker for a quantitative trait. The molecular structure of P was characterized to understand the selection of white seeds during bean diversification and the relationship of P to seed weight. P was identified by homology searches, a genome-wide association study (GWAS) and gene remodeling, and confirmed by gene silencing. Allelic variation was assessed by a combination of resequencing and marker development, and the relationship between P and seed weight was assessed by a GWAS study. P is a member of clade B of subclass IIIf of plant basic helix-loop-helix (bHLH) proteins. Ten race-specific P alleles conditioned the white seed phenotype, and each causative mutation affected at least one bHLH domain required for color expression. GWAS analysis confirmed the classic association of P with seed weight. In common bean, white seeds are the result of convergent evolution and, among plant species, orthologous convergence on a single transcription factor gene was observed.

Lotus japonicus NOOT-BOP-COCH-LIKE1 is essential for nodule, nectary, leaf and flower development

The NOOT-BOP-COCH-LIKE (NBCL) genes are orthologs of Arabidopsis thaliana BLADE-ON-PETIOLE1/2. The NBCLs are developmental regulators essential for plant shaping, mainly through the regulation of organ boundaries, the promotion of lateral organ differentiation and the acquisition of organ identity. In addition to their roles in leaf, stipule and flower development, NBCLs are required for maintaining the identity of indeterminate nitrogen-fixing nodules with persistent meristems in legumes. In legumes forming determinate nodules, without persistent meristem, the roles of NBCL genes are not known. We thus investigated the role of Lotus japonicus NOOT-BOP-COCH-LIKE1 (LjNBCL1) in determinate nodule identity and studied its functions in aerial organ development using LORE1 insertional mutants and RNA interference-mediated silencing approaches. In Lotus, LjNBCL1 is involved in leaf patterning and participates in the regulation of axillary outgrowth. Wild-type Lotus leaves are composed of five leaflets and possess a pair of nectaries at the leaf axil. Legumes such as pea and Medicago have a pair of stipules, rather than nectaries, at the base of their leaves. In Ljnbcl1, nectary development is abolished, demonstrating that nectaries and stipules share a common evolutionary origin. In addition, ectopic roots arising from nodule vascular meristems and reorganization of the nodule vascular bundle vessels were observed on Ljnbcl1 nodules. This demonstrates that NBCL functions are conserved in both indeterminate and determinate nodules through the maintenance of nodule vascular bundle identity. In contrast to its role in floral patterning described in other plants, LjNBCL1 appears essential for the development of both secondary inflorescence meristem and floral meristem.

Gene-based SNP identification and validation in soybean using next-generation transcriptome sequencing

Gene-based molecular markers are increasingly used in crop breeding programs for marker-assisted selection. However, identification of genetic variants associated with important agronomic traits has remained a difficult task in soybean. RNA-Seq provides an efficient way, other than assessing global expression variations of coding genes, to discover gene-based SNPs at the whole genome level. In this study, RNA isolated from four soybean accessions each with three replications was subjected to high-throughput sequencing and a range of 44.2-65.9 million paired-end reads were generated for each library. A total of 75,209 SNPs were identified among different genotypes after combination of replications, 89.1% of which were located in expressed regions and 27.0% resulted in amino acid changes. GO enrichment analysis revealed that most significant enriched genes with nonsynonymous SNPs were involved in ribonucleotide binding or catalytic activity. Of 22 SNPs subjected to PCR amplification and Sanger sequencing, all of them were validated. To test the utility of identified SNPs, these validated SNPs were also assessed by genotyping a relative large population with 393 wild and cultivated soybean accessions. These SNPs identified by RNA-Seq provide a useful resource for genetic and genomic studies of soybean. Moreover, the collection of nonsynonymous SNPs annotated with their predicted functional effects also provides a valuable asset for further discovery of genes, identification of gene variants, and development of functional markers.

Combining QTL-seq and linkage mapping to fine map a wild soybean allele characteristic of greater plant height

BACKGROUND

Plant height (PH) is an important agronomic trait and is closely related to yield in soybean [Glycine max (L.) Merr.]. Previous studies have identified many QTLs for PH. Due to the complex genetic background of PH in soybean, there are few reports on its fine mapping.

RESULTS

In this study, we used a mapping population derived from a cross between a chromosome segment substitution line CSSL3228 (donor N24852 (G. Soja), a receptor NN1138-2 (G. max)) and NN1138-2 to fine map a wild soybean allele of greater PH by QTL-seq and linkage mapping. We identified a QTL for PH in a 1.73 Mb region on soybean chromosome 13 through QTL-seq, which was confirmed by SSR marker-based classical QTL mapping in the mapping population. The linkage analysis showed that the QTLs of PH were located between the SSR markers BARCSOYSSR_13_1417 and BARCSOYSSR_13_1421 on chromosome 13, and the physical distance was 69.3 kb. RT-PCR and sequence analysis of possible candidate genes showed that Glyma.13 g249400 revealed significantly higher expression in higher PH genotypes, and the gene existed 6 differences in the amino acids encoding between the two parents.

CONCLUSIONS

Data presented here provide support for Glyma.13 g249400 as a possible candidate genes for higher PH in wild soybean line N24852.

Soybean MADS-box gene GmAGL1 promotes flowering via the photoperiod pathway

BACKGROUND

The MADS-box transcription factors are an ancient family of genes that regulate numerous physiological and biochemical processes in plants and facilitate the development of floral organs. However, the functions of most of these transcription factors in soybean remain unknown.

RESULTS

In this work, a MADS-box gene, GmAGL1, was overexpressed in soybean. Phenotypic analysis showed that GmAGL1 overexpression not only resulted in early maturation but also promoted flowering and affected petal development. Furthermore, the GmAGL1 was much more effective at promoting flowering under long-day conditions than under short-day conditions. Transcriptome sequencing analysis showed that before flowering, the photoperiod pathway photoreceptor CRY2 and several circadian rhythm genes, such as SPA1, were significantly down-regulated, while some other flowering-promoting circadian genes, such as GI and LHY, and downstream genes related to flower development, such as FT, LEAFY, SEP1, SEP3, FUL, and AP1, were up-regulated compared with the control. Other genes related to the flowering pathway were not noticeably affected.

CONCLUSIONS

The findings reported herein indicate that GmAGL1 may promote flowering mainly through the photoperiod pathway. Interestingly, while overexpression of GmAGL1 promoted plant maturity, no reduction in seed production or oil and protein contents was observed.

Elevation of soybean seed oil content through selection for seed coat shininess

Many leguminous species have adapted their seed coat with a layer of powdery bloom that contains hazardous allergens and makes the seeds less visible, offering duel protection against potential predators 1 . Nevertheless, a shiny seed surface without bloom is desirable for human consumption and health, and is targeted for selection under domestication. Here we show that seed coat bloom in wild soybeans is mainly controlled by Bloom1 (B1), which encodes a transmembrane transporter-like protein for biosynthesis of the bloom in pod endocarp. The transition from the 'bloom' to 'no-bloom' phenotypes is associated with artificial selection of a nucleotide mutation that naturally occurred in the coding region of B1 during soybean domestication. Interestingly, this mutation not only 'shined' the seed surface, but also elevated seed oil content in domesticated soybeans. Such an elevation of oil content in seeds appears to be achieved through b1-modulated upregulation of oil biosynthesis in pods. This study shows pleiotropy as a mechanism underlying the domestication syndrome 2 , and may pave new strategies for development of soybean varieties with increased seed oil content and reduced seed dust.

Ectopic expression of IiFUL isolated from Isatis indigotica could change the reproductive growth of Arabidopsis thaliana

The coding sequence of IiFUL in Isatis indigotica was isolated and was used in transformation of Arabidopsis. IiFUL overexpressing Arabidopsis plants exhibited early flowering phenotype, accompanied with the reduction of flower number and the production of terminal flower on the top of the main stems. In development process, the flowers located on the top of the main stems generated a lot of variations in phenotype, including abnormal swelling of pistil, withering and numerical change of stamens and petals, appearance of stigmatoid tissues and naked ovules at the margin or inside of sepals. Besides, secondary flower could be formed within the flowers on the top of the main stems. These observations illustrated that IiFUL mainly affected the development of inflorescence meristems and pistils, but its ectopic expression could also disturb the normal growth of other floral organs. Moreover, the fertile siliques produced by the lateral inflorescences of IiFUL overexpressing Arabidopsis plants showed indehiscent phenotype, and the shape of the cauline leaves was changed significantly. The results of quantitative real-time PCR revealed that higher transcriptional levels of IiFUL could be detected in flowers and silicles of I. indigotica. In comprehensive consideration of the previous reports about the dehiscence phenotype of Arabidopsis siliques and the fact that the siliques of IiFUL overexpressing Arabidopsis plants were indehiscent in the present work, it can be speculated that high expression of IiFUL in pericarp is likely the reason why the silicles of I. indigotica possess an indehiscent phenotype.

Photoperiodism dynamics during the domestication and improvement of soybean

Soybean (Glycine max) is a facultative short-day plant with a sensitive photoperiod perception and reaction system, which allows it to adjust its physiological state and gene regulatory networks to seasonal and diurnal changes in environmental conditions. In the past few decades, soybean cultivation has spread from East Asia to areas throughout the world. Biologists and breeders must now confront the challenge of understanding the molecular mechanism of soybean photoperiodism and improving agronomic traits to enable this important crop to adapt to geographical and environmental changes. In this review, we summarize the genetic regulatory network underlying photoperiodic responses in soybean. Genomic and genetic studies have revealed that the circadian clock, in conjunction with the light perception pathways, regulates photoperiodic flowering. Here, we provide an annotated list of 844 candidate flowering genes in soybean, with their putative biological functions. Many photoperiod-related genes have been intensively selected during domestication and crop improvement. Finally, we describe recent progress in engineering photoperiod-responsive genes for improving agronomic traits to enhance geographic adaptation in soybean, as well as future prospects for research on soybean photoperiodic responses.

Genome-wide Association Mapping of Qualitatively Inherited Traits in a Germplasm Collection

Genome-wide association (GWA) has been used as a tool for dissecting the genetic architecture of quantitatively inherited traits. We demonstrate here that GWA can also be highly useful for detecting many major genes governing categorically defined phenotype variants that exist for qualitatively inherited traits in a germplasm collection. Genome-wide association mapping was applied to categorical phenotypic data available for 10 descriptive traits in a collection of ∼13,000 soybean [ (L.) Merr.] accessions that had been genotyped with a 50,000 single nucleotide polymorphism (SNP) chip. A GWA on a panel of accessions of this magnitude can offer substantial statistical power and mapping resolution, and we found that GWA mapping resulted in the identification of strong SNP signals for 24 classical genes as well as several heretofore unknown genes controlling the phenotypic variants in those traits. Because some of these genes had been cloned, we were able to show that the narrow GWA mapping SNP signal regions that we detected for the phenotypic variants had chromosomal bp spans that, with just one exception, overlapped the bp region of the cloned genes, despite local variation in SNP number and nonuniform SNP distribution in the chip set.