Genome-wide association studies for soybean epicotyl length in two environments using 3VmrMLM.

Germination of soybean seed is the imminent vital process after sowing. The status of plumular axis and radicle determine whether soybean seed can emerge normally. Epicotyl, an organ between cotyledons and first functional leaves, is essential for soybean seed germination, seedling growth and early morphogenesis. Epicotyl length (EL) is a quantitative trait controlled by multiple genes/QTLs. Here, the present study analyzes the phenotypic diversity and genetic basis of EL using 951 soybean improved cultivars and landraces from Asia, America, Europe and Africa. 3VmrMLM was used to analyze the associations between EL in 2016 and 2020 and 1,639,846 SNPs for the identification of QTNs and QTN-by-environment interactions (QEIs)".A total of 180 QTNs and QEIs associated with EL were detected. Among them, 74 QTNs (ELS_Q) and 16 QEIs (ELS_QE) were identified to be associated with ELS (epicotyl length of single plant emergence), and 60 QTNs (ELT_Q) and 30 QEIs (ELT_QE) were identified to be associated with ELT (epicotyl length of three seedlings). Based on transcript abundance analysis, GO (Gene Ontology) enrichment and haplotype analysis, ten candidate genes were predicted within nine genic SNPs located in introns, upstream or downstream, which were supposed to be directly or indirectly involved in the process of seed germination and seedling development., Of 10 candidate genes, two of them (Glyma.04G122400 and Glyma.18G183600) could possibly affect epicotyl length elongation. These results indicate the genetic basis of EL and provides a valuable basis for specific functional studies of epicotyl traits.

Soybean genetic resources contributing to sustainable protein production.

KEY MESSAGE: Genetic resources contributes to the sustainable protein production in soybean. Soybean is an important crop for food, oil, and forage and is the main source of edible vegetable oil and vegetable protein. It plays an important role in maintaining balanced dietary nutrients for human health. The soybean protein content is a quantitative trait mainly controlled by gene additive effects and is usually negatively correlated with agronomic traits such as the oil content and yield. The selection of soybean varieties with high protein content and high yield to secure sustainable protein production is one of the difficulties in soybean breeding. The abundant genetic variation of soybean germplasm resources is the basis for overcoming the obstacles in breeding for soybean varieties with high yield and high protein content. Soybean has been cultivated for more than 5000 years and has spread from China to other parts of the world. The rich genetic resources play an important role in promoting the sustainable production of soybean protein worldwide. In this paper, the origin and spread of soybean and the current status of soybean production are reviewed; the genetic characteristics of soybean protein and the distribution of resources are expounded based on phenotypes; the discovery of soybean seed protein-related genes as well as transcriptomic, metabolomic, and proteomic studies in soybean are elaborated; the creation and utilization of high-protein germplasm resources are introduced; and the prospect of high-protein soybean breeding is described.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array.

KEY MESSAGE: We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Estimation of soybean yield parameters under lodging conditions using RGB information from unmanned aerial vehicles.

The estimation of yield parameters based on early data is helpful for agricultural policymakers and food security. Developments in unmanned aerial vehicle (UAV) platforms and sensor technology help to estimate yields efficiency. Previous studies have been based on less cultivars (<10) and ideal experimental environments, it is not available in practical production. Therefore, the objective of this study was to estimate the yield parameters of soybean (Glycine max (L.) Merr.) under lodging conditions using RGB information. In this study, 17 time point data throughout the soybean growing season in Nanchang, Jiangxi Province, China, were collected, and the vegetation index, texture information, canopy cover, and crop height were obtained by UAV-image processing. After that, partial least squares regression (PLSR), logistic regression (Logistic), random forest regression (RFR), support vector machine regression (SVM), and deep learning neural network (DNN) were used to estimate the yield parameters. The results can be summarized as follows: (1) The most suitable time point to estimate the yield was flowering stage (48 days), which was when most of the soybean cultivars flowered. (2) The multiple data fusion improved the accuracy of estimating the yield parameters, and the texture information has a high potential to contribute to the estimation of yields, and (3) The DNN model showed the best accuracy of training (R2=0.66 rRMSE=32.62%) and validation (R2=0.50, rRMSE=43.71%) datasets. In conclusion, these results provide insights into both best estimate period selection and early yield estimation under lodging condition when using remote sensing.

Transcriptome analysis reveals differing response and tolerance mechanism of EPSPS and GAT genes among transgenic soybeans.

BACKGROUND: Glyphosate is a broad-spectrum, non-selective systemic herbicide. Introduction of glyphosate tolerance genes such as EPSPS or detoxification genes such as GAT can confer glyphosate tolerance on plants. Our previous study revealed that co-expression of EPSPS and GAT genes conferred higher glyphosate tolerance without "yellow flashing". However, the plant response to glyphosate at the transcriptional level was not investigated. METHODS AND RESULTS: To investigate the glyphosate tolerance mechanism, RNA-seq was conducted using four soybean genotypes, including two non-transgenic (NT) soybeans, ZH10 and MD12, and two GM soybeans, HJ698 and ZH10-6. Differentially expressed genes (DEGs) were identified in these soybeans before and after glyphosate treatment. Similar response to glyphosate in the two NT soybeans and the different effects of glyphosate on the two GM soybeans were identified. As treatment time was prolonged, the expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was increased or declined continuously in NT soybeans, and altered slightly in HJ698. However, the expression level of some DEGs was altered in ZH10-6 at 12 hpt, while similar expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was observed in ZH10-6 at 0 hpt and 72 hpt. These observations likely explain the higher glyphosate tolerance in ZH10-6 than in HJ698 and NT soybeans. CONCLUSIONS: These results suggested that GAT and EPSPS genes together play a crucial role in response to glyphosate, the GAT gene may work at the early stage of glyphosate exposure, whereas the EPSPS gene may be activated after the uptake of glyphosate by plants. These findings will provide valuable insight for the molecular basis underlying glyphosate tolerance or glyphosate detoxication.

Creation of Early Flowering Germplasm of Soybean by CRISPR/Cas9 Technology.

Soybean is an important economic crop and a typical short-day crop, sensitive to photoperiod, and has narrow geographical adaptative region, which limit the creation of transgenic materials and reduce the breeding efficiency of new varieties. In addition, the genetic transformation efficiency of soybean is lower than that of many other crops, and the available receptor genotypes are limited. In this study, Agrobacterium-mediated transformation were used to introduce the CRISPR/Cas9 expression vector into soybean cultivar Jack and generated targeted mutants of E1 gene controlling soybean flowering. We obtained two novel types of mutations, 11 bp and 40 bp deletion at E1 coding region, respectively, and frameshift mutations produced premature translation termination codons and truncated E1 proteins, causing obvious early flowering under long day condition. In addition, no off-target effects were observed by predicting and analyzing the potential off-target sites of E1 targets. Significant decreased E1 gene expression of two novel mutants showed that the truncated E1 protein disinhibited GmFT2a/5a and increasing GmFT2a/5a gene expressions resulted obvious early flowering. Homozygous trans-clean mutants without T-DNA elements were also obtained and showed early flowering under long day condition. The photo-insensitive soybean transformation receptor we created laid a foundation for breeding excellent transgenic receptors suitable for high latitudes.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis.

KEY MESSAGE: GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Development and utilization of a new chemically-induced soybean library with a high mutation density .

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.

Identification of Genomic Insertion and Flanking Sequence of G2-EPSPS and GAT Transgenes in Soybean Using Whole Genome Sequencing Method.

Molecular characterization of sequence flanking exogenous fragment insertion is essential for safety assessment and labeling of genetically modified organism (GMO). In this study, the T-DNA insertion sites and flanking sequences were identified in two newly developed transgenic glyphosate-tolerant soybeans GE-J16 and ZH10-6 based on whole genome sequencing (WGS) method. More than 22.4 Gb sequence data (∼21 × coverage) for each line was generated on Illumina HiSeq 2500 platform. The junction reads mapped to boundaries of T-DNA and flanking sequences in these two events were identified by comparing all sequencing reads with soybean reference genome and sequence of transgenic vector. The putative insertion loci and flanking sequences were further confirmed by PCR amplification, Sanger sequencing, and co-segregation analysis. All these analyses supported that exogenous T-DNA fragments were integrated in positions of Chr19: 50543767-50543792 and Chr17: 7980527-7980541 in these two transgenic lines. Identification of genomic insertion sites of G2-EPSPS and GAT transgenes will facilitate the utilization of their glyphosate-tolerant traits in soybean breeding program. These results also demonstrated that WGS was a cost-effective and rapid method for identifying sites of T-DNA insertions and flanking sequences in soybean.

Co-expression of G2-EPSPS and glyphosate acetyltransferase GAT genes conferring high tolerance to glyphosate in soybean.

Glyphosate is a widely used non-selective herbicide with broad spectrum of weed control around the world. At present, most of the commercial glyphosate tolerant soybeans utilize glyphosate tolerant gene CP4-EPSPS or glyphosate acetyltransferase gene GAT separately. In this study, both glyphosate tolerant gene G2-EPSPS and glyphosate degraded gene GAT were co-transferred into soybean and transgenic plants showed high tolerance to glyphosate. Molecular analysis including PCR, Sothern blot, qRT-PCR, and Western blot revealed that target genes have been integrated into genome and expressed effectively at both mRNA and protein levels. Furthermore, the glyphosate tolerance analysis showed that no typical symptom was observed when compared with a glyphosate tolerant line HJ06-698 derived from GR1 transgenic soybean even at fourfold labeled rate of Roundup. Chlorophyll and shikimic acid content analysis of transgenic plant also revealed that these two indexes were not significantly altered after glyphosate application. These results indicated that co-expression of G2-EPSPS and GAT conferred high tolerance to the herbicide glyphosate in soybean. Therefore, combination of tolerant and degraded genes provides a new strategy for developing glyphosate tolerant transgenic crops.