Search

Introduction Legume crops are an important source of protein and oil for human health and in fixing atmospheric N2 for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement. Objectives Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution. Methods We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies. Results We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of “QTL-hotspot” identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, ‘Legumepedia’. Conclusion This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

Plant parasitic nematodes are a major yield-limiting factor of soybean in the United States and Canada. It has been indicated that soybean cyst nematode (SCN; Heterodera glycines Ichinohe) and reniform nematode (RN; Rotylenchulus reniformis Linford and Oliveira) resistance could be genetically related. For many years, fragmentary data have shown this relationship. This report evaluates RN reproduction on 418 plant introductions (PIs) selected from the U.S. Department of Agriculture Soybean Germplasm Collection with reported SCN resistance. The germplasm was divided into two tests of 214 PIs reported as resistant and 204 PIs reported as moderately resistant to SCN. The defining and reporting of RN resistance changed several times in the last 30 years, causing inconsistencies in RN resistance classification among multiple experiments. Comparison of four RN resistance classification methods was performed: (i) ≤10% as compared with the susceptible check, (ii) using normalized reproduction index (RI) values, and using (iii) transformed data log10(x), and (iv) transformed data log10(x + 1) in an optimal univariate k-means clustering analysis. The method of transformed data log10(x) was selected as the most accurate for classification of RN resistance. Among 418 PIs with reported SCN resistance, the log10(x) method grouped 59 PIs (15%) as resistant and 130 PIs (31%) as moderately resistant to RN. Genotyping of a subset of the most resistant PIs to both nematode species revealed their strong correlation with rhg1-a allele. This research identified genotypes with resistance to two nematode species and potential new sources of RN resistance that could be valuable to breeders in developing resistant cultivars.

The qSCN10 locus with broad-spectrum SCN resistance was fine-mapped to a 379-kb region on chromosome 10 in soybean accession PI 567516C. Candidate genes and potential application benefits of this locus were discussed. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pests of soybean, causing significant yield losses worldwide every year. Genetic resistance has been the major strategy to control this pest. However, the overuse of the same genetic resistance derived primarily from PI 88788 has led to the genetic shifts in nematode populations and resulted in the reduced effectiveness in soybean resistance to SCN. Therefore, novel genetic resistance resources, especially those with broad-spectrum resistance, are needed to develop new resistant cultivars to cope with the genetic shifts in nematode populations. In this study, a quantitative trait locus (QTL) qSCN10 previously identified from a soybean landrace PI 567516C was confirmed to confer resistance to multiple SCN HG Types. This QTL was further fine-mapped to a 379-kb region. There are 51 genes in this region. Four of them are defense-related and were regulated by SCN infection, suggesting their potential role in mediating resistance to SCN. The phylogenetic and haplotype analyses of qSCN10 as well as other information indicate that this locus is different from other reported resistance QTL or genes. There was no yield drag or other unfavorable traits associated with this QTL when near-isogenic lines with and without qSCN10 were tested in a SCN-free field. Therefore, our study not only provides further insight into the genetic basis of soybean resistance to SCN, but also identifies a novel genetic resistance resource for breeding soybean for durable, broad-spectrum resistance to this pest.

The qSCN18 QTL from PI 56756C was confirmed and fine-mapped to improve soybean resistance to the SCN population HG Type 2.5.7 using near-isogenic lines carrying recombination crossovers within the QTL region. The QTL underlying resistance was fine-mapped to a 166-Kbp region on chromosome 18, and the candidate genes were selected based on genomic analyses. Soybean cyst nematode (SCN, Heterodera glycines, Ichinohe) is the most devastating pathogen of soybean. Understanding the genetic basis of SCN resistance is crucial for managing this parasite in the field. Two major loci, rhg1 and Rhg4, were previously characterized as valuable resources for SCN resistance. However, their continuous use has caused shifts in the virulence of SCN populations, which can overcome the resistance conferred by these two major loci. Reduced effectiveness became a major concern in the soybean industry due to continuous use of rhg1 for decades. Thus, it is imperative to identify sources of SCN resistance for durable SCN management. A novel QTL qSCN18 was identified in PI567516C. To fine-map qSCN18 and identify resistance genes, a large backcross population was developed. Nineteen near-isogenic lines (NILs) carrying recombination crossovers within the QTL region were identified. The first phase of fine-mapping narrowed the QTL region to 549-Kbp, whereas the second phase confined the region to 166-Kbp containing 23 genes. Two flanking markers, MK-1 and MK-6, were developed and validated to detect the presence of the qSCN18 resistance allele. Haplotype analysis clustered the fine-mapped qSCN18 region from PI 567516C with the cqSCN-007 locus previously mapped in the wild soybean accession PI 468916. The NILs were developed to further characterize the causal gene(s) harbored in this QTL. This study also confirmed the previously identified qSCN18. The results will facilitate marker-assisted selection (MAS) introducing the qSCN18 locus from PI 567516C into high-yielding soybean cultivars with durable resistance to SCN.

Recent growth in crop genomic and trait data have opened opportunities for the application of novel approaches to accelerate crop improvement. Machine learning and deep learning are at the forefront of prediction-based data analysis. However, few approaches for genotype to phenotype prediction compare machine learning with deep learning and further interpret the models that support the predictions. This study uses genome wide molecular markers and traits across 1110 soybean individuals to develop accurate prediction models. For 13/14 sets of predictions, XGBoost or random forest outperformed deep learning models in prediction performance. Top ranked SNPs by F-score were identified from XGBoost, and with further investigation found overlap with significantly associated loci identified from GWAS and previous literature. Feature importance rankings were used to reduce marker input by up to 90%, and subsequent models maintained or improved their prediction performance. These findings support interpretable machine learning as an approach for genomic based prediction of traits in soybean and other crops.

A whole-genome resequencing-derived SNP dataset identified six quantitative trait loci (QTL) significantly associated with colonization of soybean by an arbuscular mycorrhizal fungus (Rhizophagus intraradices). Candidate genes identified in these QTL regions include homologs to known nodulin protein families and other symbiosis-specific genes. Arbuscular mycorrhizal fungi (AMF) form associations with over 80% of all terrestrial plant species and assist their host plants by increasing their nutrient uptake, drought tolerance, and resilience against pathogens and pests. Genotypic variation of crop plants to AMF colonization has been identified in crops, including soybean; however, the genetics controlling levels of AMF colonization in soybean are unknown. The overall goal of our study was to identify genomic regions associated with mycorrhizal colonization in soybean using genome-wide association analysis. A diverse panel of 350 exotic soybean genotypes inoculated with Rhizophagus intraradices were microscopically evaluated for root colonization using a modified gridline intersect method. Root colonization differed significantly (P

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

In the last decade, more than 70 quantitative trait loci (QTL) related to soybean [Glycine max (L.) Merr.] partial resistance (PR) against Phytophthora sojae have been identified by genome-wide association studies (GWAS). However, most of them have either a minor effect on the resistance level or are specific to a single phenotypic variable or one isolate, thereby limiting their use in breeding programs. In this study, we have used an analytical approach combining (a) the phenotypic characterization of a diverse panel of 357 soybean accessions for resistance to P. sojae captured through a single variable, corrected dry weight; (b) a new hydroponic assay allowing the inoculation of a combination of P. sojae isolates covering the spectrum of commercially relevant Rps genes; and (c) exhaustive genotyping through whole-genome resequencing (WGS). This led to the identification of a novel P. sojae resistance QTL with a relatively major effect compared with the previously reported QTL. The QTL interval, spanning ∼500 kb on chromosome (Chr) 15, does not colocalize with previously reported QTL for P. sojae resistance. Plants carrying the favorable allele at this QTL were 60% more resistant. Eight genes were found to reside in the linkage disequilibrium (LD) block containing the peak single-nucleotide polymorphism (SNP) including Glyma.15G217100, which encodes a major latex protein (MLP)-like protein, with a functional annotation related to pathogen resistance. Expression analysis of Glyma.15G217100 indicated that it was nearly eight times more highly expressed in a group of plant introductions (PIs) carrying the resistant (R) allele compared with those carrying the susceptible (S) allele within a short period after inoculation. These results offer new and valuable options to develop improved soybean cultivars with broad resistance to P. sojae through marker-assisted selection.

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

Efficient crop improvement depends on the application of accurate genetic information contained in diverse germplasm resources. Here we report a reference-grade genome of wild soybean accession W05, with a final assembled genome size of 1013.2 Mb and a contig N50 of 3.3 Mb. The analytical power of the W05 genome is demonstrated by several examples. First, we identify an inversion at the locus determining seed coat color during domestication. Second, a translocation event between chromosomes 11 and 13 of some genotypes is shown to interfere with the assignment of QTLs. Third, we find a region containing copy number variations of the Kunitz trypsin inhibitor (KTI) genes. Such findings illustrate the power of this assembly in the analysis of large structural variations in soybean germplasm collections. The wild soybean genome assembly has wide applications in comparative genomic and evolutionary studies, as well as in crop breeding and improvement programs., Wild relatives of crop plants are invaluable germplasm for genetic improvement. Here, Xie et al. report a reference-grade wild soybean genome and show that it can be used to identify structural variation and refine quantitative trait loci.

Summary Soybean cyst nematode (SCN) is the most devastating plant‐parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1‐b locus, while ‘Peking’ requires rhg1‐a and Rhg4 for SCN resistance. In the present study, whole genome re‐sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG‐Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and ‘Cloud’. At least 5.6 copies of the PI88788‐type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a ‘Peking’‐type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad‐based resistance to SCN and providing new insight into epistasis, haplotype‐compatibility, CNV, promoter variation and its impact on broad‐based disease resistance in plants.

Two novel QTL for resistance to Pythium ultimum var. ultimum were identified in soybean using an Illumina SNP Chip and whole genome re-sequencing. Pythium ultimum var. ultimum is one of numerous Pythium spp. that causes severe pre- and post-emergence damping-off of seedlings and root rot of soybean [Glycine max (L.) Merr.]. The objective of this research was to identify quantitative trait loci (QTL) for resistance to P. ultimum var. ultimum in a recombinant inbred line population derived from a cross of 'Magellan' (moderately resistant) and PI 438489B (susceptible). Two different mapping approaches were utilized: the universal soybean linkage panel (USLP 1.0) and the bin map constructed from whole genome re-sequencing (WGRS) technology. Two genomic regions associated with variation in three disease-related parameters were detected using both approaches, with the bin map providing higher resolution. Using WGRS, the first QTL were mapped within a 350-kbp region on Chr. 6 and explained 7.5-13.5% of the phenotypic variance. The second QTL were positioned in a 260-kbp confidence interval on Chr. 8 and explained 6.3-16.8% of the phenotypic variation. Candidate genes potentially associated with disease resistance were proposed. High-resolution genetic linkage maps with a number of significant SNP markers could benefit marker-assisted breeding and dissection of the molecular mechanisms underlying soybean resistance to Pythium damping-off in 'Magellan.' Additionally, the outputs of this study may encourage more screening of diverse soybean germplasm and utilization of genome-wide association studies to understand the genetic basis of quantitative disease resistance.

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

Soybean acyl-ACP thioesterase gene family have been characterized; GmFATA1A mutants were discovered to confer high oleic acid, while GmFATB mutants presented low palmitic and high oleic acid seed content. Soybean oil stability and quality are primarily determined by the relative proportions of saturated versus unsaturated fatty acids. Commodity soybean typically contains 11% palmitic acid, as the primary saturated fatty acids. Reducing palmitic acid content is the principal approach to minimize the levels of saturated fatty acids in soybean. Though high palmitic acid enhances oxidative stability of soybean oil, it is negatively correlated with oil and oleic acid content and can cause coronary heart diseases for humans. For plants, acyl-acyl carrier protein (ACP) thioesterases (TEs) are a group of enzymes to hydrolyze acyl group and release free fatty acid from plastid. Among them, GmFATB1A has become the main target to genetically reduce the palmitic acid content in soybean. However, the role of members in soybean acyl-ACP thioesterase gene family is largely unknown. In this study, we characterized two classes of TEs, GmFATA, and GmFATB in soybean. We also denominated two GmFATA members and discovered six additional members that belong to GmFATB gene family through phylogenetic, syntenic, and in silico analysis. Using TILLING-by-Sequencing

Reniform nematode (RN, Rotylenchulus reniformis Linford & Oliveira) has emerged as one of the most important plant parasitic nematodes of soybean [Glycine max (L.) Merr.]. Planting resistant varieties is the most effective strategy for nematode management. The objective of this study was to identify quantitative trait loci (QTL) for RN resistance in an exotic soybean line, PI 438489B, using two linkage maps constructed from the Universal Soybean Linkage Panel (USLP 1.0) and next‐generation whole‐genome resequencing (WGRS) technology. Two QTL controlling RN resistance were identified—the soybean cyst nematode (SCN, Heterodera glycines) resistance gene GmSNAP18 at the rhg1 locus and its paralog GmSNAP11. Strong association between resistant phenotype and haplotypes of the GmSNAP11 and GmSNAP18 was observed. The results indicated that GmSNAP11 possibly could have epistatic effect on GmSNAP18, or vice versa, with the presence of a significant correlation in RN resistance of rhg1‐a GmSNAP18 vs. rhg1‐b GmSNAP18. Most importantly, our preliminary data suggested that GmSNAP18 and GmSNAP11 proteins physically interact in planta, suggesting that they belong to the same pathway for resistance. Unlike GmSNAP18, no indication of GmSNAP11 copy number variation was found. Moreover, gene‐based single nucleotide polymorphism (SNP) markers were developed for rapid detection of RN or SCN resistance at these loci. Our analysis substantiates synergic interaction between GmSNAP11 and GmSNAP18 genes and confirms their roles in RN as well as SCN resistance. These results could contribute to a better understanding of evolution and subfunctionalization of genes conferring resistance to multiple nematode species and provide a framework for further investigations.

The gene content of plants varies between individuals of the same species due to gene presence/absence variation, and selection can alter the frequency of specific genes in a population. Selection during domestication and breeding will modify the genomic landscape, though the nature of these modifications is only understood for specific genes or on a more general level (e.g., by a loss of genetic diversity). Here we have assembled and analyzed a soybean (Glycine spp.) pangenome representing more than 1,000 soybean accessions derived from the USDA Soybean Germplasm Collection, including both wild and cultivated lineages, to assess genomewide changes in gene and allele frequency during domestication and breeding. We identified 3,765 genes that are absent from the Lee reference genome assembly and assessed the presence/absence of all genes across this population. In addition to a loss of genetic diversity, we found a significant reduction in the average number of protein-coding genes per individual during domestication and subsequent breeding, though with some genes and allelic variants increasing in frequency associated with selection for agronomic traits. This analysis provides a genomic perspective of domestication and breeding in this important oilseed crop.

We report characteristics of soybean genetic diversity and structure from the resequencing of 481 diverse soybean accessions, comprising 52 wild (Glycine soja) selections and 429 cultivated (Glycine max) varieties (landraces and elites). This data was used to identify 7.8 million SNPs, to predict SNP effects relative to genic regions, and to identify the genetic structure, relationships, and linkage disequilibrium. We found evidence of distinct, mostly independent selection of lineages by particular geographic location. Among cultivated varieties, we identified numerous highly conserved regions, suggesting selection during domestication. Comparisons of these accessions against the whole U.S. germplasm genotyped with the SoySNP50K iSelect BeadChip revealed that over 95% of the re-sequenced accessions have a high similarity to their SoySNP50K counterparts. Probable errors in seed source or genotype tracking were also identified in approximately 5% of the accessions., Measurement(s) genetic variation • SNP • Linkage Disequilibrium Technology Type(s) DNA sequencing • bioinformatics analysis • computational phylogenetic analysis Sample Characteristic - Organism Glycine soja • Glycine max Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13568552

Alternative splicing (AS) is a common post-transcriptional regulatory mechanism that modulates gene expression to increase proteome diversity. Increasing evidence indicates that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate drought responses in soybean remains poorly understood. In this study, we performed a genome-wide analysis of AS events in soybean (Glycine max) roots grown under various drought conditions using the high-throughput RNA-sequencing method, identifying 385, 989, 1429, and 465 AS events that were significantly differentially spliced under very mild drought stress, mild drought stress, severe drought stress, and recovery after severe drought conditions, respectively. Among them, alternative 3&prime, splice sites and skipped exons were the major types of AS. Overall, 2120 genes that experienced significant AS regulation were identified from these drought-treated root samples. Gene Ontology term analysis indicated that the AS regulation of binding activity has vital roles in the drought response of soybean root. Notably, the genes encoding splicing regulatory factors in the spliceosome pathway and mRNA surveillance pathway were enriched according to the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Splicing regulatory factor-related genes in soybean root also responded to drought stress and were alternatively spliced under drought conditions. Taken together, our data suggest that drought-responsive AS acts as a direct or indirect mode to regulate drought response of soybean roots. With further in-depth research of the function and mechanism of AS in the process of abiotic stress, these results will provide a new strategy for enhancing stress tolerance of plants.

Aluminium (Al) toxicity inhibits soybean root growth, leading to insufficient water and nutrient uptake. Two soybean lines (‘Magellan’ and PI 567731) were identified differing in Al tolerance, as determined by primary root length ratio, total root length ratio, and root tip number ratio under Al stress. Serious root necrosis was observed in PI 567731, but not in Magellan under Al stress. An F8 recombinant inbred line population derived from a cross between Magellan and PI 567731 was used to map the quantitative trait loci (QTL) for Al tolerance. Three QTL on chromosomes 3, 13, and 20, with tolerant alleles from Magellan, were identified. qAl_Gm13 and qAl_Gm20 explained large phenotypic variations (13–27%) and helped maintain root elongation and initiation under Al stress. In addition, qAl_Gm13 and qAl_Gm20 were confirmed in near-isogenic backgrounds and were identified to epistatically regulate Al tolerance via internal detoxification instead of Al3+ exclusion. Phylogenetic and pedigree analysis identified the tolerant alleles of both loci derived from the US ancestral line, A.K.[FC30761], originally from China. Our results provide novel genetic resources for breeding Al-tolerant soybean and suggest that internal detoxification contributes to soybean tolerance to excessive soil Al.

Transporters, a class of membrane proteins that facilitate exchange of solutes including diverse molecules and ions across the cellular membrane, are vital component for the survival of all organisms. Understanding plant transporters is important to get insight of the basic cellular processes, physiology, and molecular mechanisms including nutrient uptake, signaling, response to external stress, and many more. In this regard, extensive analysis of transporters predicted in soybean and other plant species was performed. In addition, an integrated database for soybean transporter protein, SoyTD, was developed that will facilitate the identification, classification, and extensive characterization of transporter proteins by integrating expression, gene ontology, conserved domain and motifs, gene structure organization, and chromosomal distribution features. A comprehensive analysis was performed to identify highly confident transporters by integrating various prediction tools. Initially, 7541 transmembrane (TM) proteins were predicted in the soybean genome; out of these, 3306 non-redundant transporter genes carrying two or more transmembrane domains were selected for further analysis. The identified transporter genes were classified according to a standard transporter classification (TC) system. Comparative analysis of transporter genes among 47 plant genomes provided insights into expansion and duplication of transporter genes in land plants. The whole genome resequencing (WGRS) and tissue-specific transcriptome datasets of soybean were integrated to investigate the natural variants and expression profile associated with transporter(s) of interest. Overall, SoyTD provides a comprehensive interface to study genetic and molecular function of soybean transporters. SoyTD is publicly available at http://artemis.cyverse.org/soykb_dev/SoyTD/.

Soybean [Glycine max (L.) Merr.] seed is a valuable source of protein and oil worldwide. Traditionally, the natural variations were heavily used in conventional soybean breeding programs to select desired traits. However, traditional plant breeding is encumbered with low frequencies of spontaneous mutations. In mutation breeding, genetic variations from induced mutations provide abundant sources of alterations in important soybean traits, this facilitated the development of soybean germplasm with modified seed composition traits to meet the different needs of end users. In this study, a total of 2366 &lsquo, Forrest&rsquo, derived M2 families were developed for both forward and reverse genetic studies. A subset of 881 M3 families was forward genetically screened to measure the contents of protein, oil, carbohydrates, and fatty acids. A total of 14 mutants were identified to have stable seed composition phenotypes observed in both M3 and M4 generations. Correlation analyses have been conducted among ten seed composition traits and compared to a collection of 103 soybean germplasms. Mainly, ethyl methanesulfonate (EMS) mutagenesis had a strong impact on the seed-composition correlation that was observed among the 103 soybean germplasms, which offers multiple benefits for the soybean farmers and industry to breed for desired multiple seed phenotypes.

Soybean is a widely grown grain legume and one of the most important economic crop species. Brassinosteroids play a crucial role in plant vegetative growth and reproductive development. However, it remains unclear how BRs regulate the developmental processes in soybean, and the molecular mechanism underlying soybean early development is largely unexplored. In this study, we first characterized how soybean early vegetative growth was specifically regulated by the BR biosynthesis inhibitor propiconazole; this characterization included shortened root and shoot lengths, reduced leaf area, and decreased chlorophyll content. In addition, the growth inhibition induced by Pcz could be rescued by exogenous brassinolide application. The RNA-seq technique was employed to investigate the BR regulatory networks during soybean early vegetative development. Identification and analysis of differentially expressed genes indicated that BRs orchestrate a wide range of cellular activities and biological processes in soybean under various BR concentrations. The regulatory networks between BRs and multiple hormones or stress-related pathways were investigated. The results provide a comprehensive view of the physiological functions of BRs and new insights into the molecular mechanisms at the transcriptional level of BR regulation of soybean early development.

Developing crops with better root systems is a promising strategy to ensure productivity in both optimum and stress environments. Root system architectural traits in 397 soybean accessions were characterized and a high-density single nucleotide polymorphisms (SNPs)-based genome-wide association study was performed to identify the underlying genes associated with root structure. SNPs associated with root architectural traits specific to landraces and elite germplasm pools were detected. Four loci were detected in landraces for lateral root number (LRN) and distribution of root thickness in diameter Class I with a major locus on chromosome 16. This major loci was detected in the coding region of unknown protein, and subsequent analyses demonstrated that root traits are affected with mutated haplotypes of the gene. In elite germplasm pool, 3 significant SNPs in alanine-glyoxalate aminotransferase, Leucine-Rich Repeat receptor/No apical meristem, and unknown functional genes were found to govern multiple traits including root surface area and volume. However, no major loci were detected for LRN in elite germplasm. Nucleotide diversity analysis found evidence of selective sweeps around the landraces LRN gene. Soybean accessions with minor and mutated allelic variants of LRN gene were found to perform better in both water-limited and optimal field conditions.

Summary The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high‐resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K‐SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K‐SNP markers. Importantly, the bin map derived from next‐generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication‐related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome‐wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL‐based haplotype using whole‐genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

The objective of this study was to determine the genetic relationship between the oleic acid and protein content. The genotypes having high oleic acid and elevated protein (HOEP) content were crossed with five elite lines having normal oleic acid and average protein (NOAP) content. The selected accessions were grown at six environments in three different locations and phenotyped for protein, oil, and fatty acid components. The mean protein content of parents, HOEP, and NOAP lines was 34.6%, 38%, and 34.9%, respectively. The oleic acid concentration of parents, HOEP, and NOAP lines was 21.7%, 80.5%, and 20.8%, respectively. The HOEP plants carried both FAD2-1A (S117N) and FAD2-1B (P137R) mutant alleles contributing to the high oleic acid phenotype. Comparative genome analysis using whole-genome resequencing data identified six genes having single nucleotide polymorphism (SNP) significantly associated with the traits analyzed. A single SNP in the putative gene Glyma.10G275800 was associated with the elevated protein content, and palmitic, oleic, and linoleic acids. The genes from the marker intervals of previously identified QTL did not carry SNPs associated with protein content and fatty acid composition in the lines used in this study, indicating that all the genes except Glyma.10G278000 may be the new genes associated with the respective traits.

A major novel quantitative disease resistance locus, qRfg_Gm06, for Fusarium graminearum was genetically mapped to chromosome 6. Genomic-assisted haplotype analysis within this region identified three putative candidate genes. Fusarium graminearum causes seed, root rot, and seedling damping-off in soybean which contributes to reduced stands and yield. A cultivar Magellan and PI 567516C were identified with low and high levels of partial resistance to F. graminearum, respectively. Quantitative disease resistance loci (QDRL) were mapped with 241 F7:8 recombinant inbred lines (RILs) derived from a cross of Magellan × PI 567516C. Phenotypic evaluation for resistance to F. graminearum used the rolled towel assay in a randomized incomplete block design. The genetic map was constructed from 927 polymorphic single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers. One major QDRL qRfg_Gm06 was detected and mapped to chromosome 6 with a LOD score of 20.3 explaining 40.2% of the total phenotypic variation. This QDRL was mapped to a ~400 kb genomic region of the Williams 82 reference genome. Genome mining of this region identified 14 putative candidate disease resistance genes. Haplotype analysis of this locus using whole genome re-sequencing (WGRS) of 106 diverse soybean lines narrowed the list to three genes. A SNP genotyping Kompetitive allele-specific PCR (KASP) assay was designed for one of the genes and was validated in a subset of the RILs and all 106 diverse lines.

Summary Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole‐genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single‐nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1‐Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene‐centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss‐of‐function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Flood has an important effect on plant growth by affecting their physiologic and biochemical properties. Soybean is one of the main cultivated crops in the world and the United States is one of the largest soybean producers. However, soybean plant is sensitive to flood stress that may cause slow growth, low yield, small crop production and result in significant economic loss. Therefore, it is critical to develop soybean cultivars that are tolerant to flood. One of the current bottlenecks in developing new crop cultivars is slow and inaccurate plant phenotyping that limits the genetic gain. This study aimed to develop a low-cost 3D imaging system to quantify the variation in the growth and biomass of soybean due to flood at its early growth stages. Two cultivars of soybeans, i.e. flood tolerant and flood sensitive, were planted in plant pots in a controlled greenhouse. A low-cost 3D imaging system was developed to take measurements of plant architecture including plant height, plant canopy width, petiole length, and petiole angle. It was found that the measurement error of the 3D imaging system was 5.8% in length and 5.0% in angle, which was sufficiently accurate and useful in plant phenotyping. Collected data were used to monitor the development of soybean after flood treatment. Dry biomass of soybean plant was measured at the end of the vegetative stage (two months after emergence). Results show that four groups had a significant difference in plant height, plant canopy width, petiole length, and petiole angle. Flood stress at early stages of soybean accelerated the growth of the flood-resistant plants in height and the petiole angle, however, restrained the development in plant canopy width and the petiole length of flood-sensitive plants. The dry biomass of flood-sensitive plants was near two to three times lower than that of resistant plants at the end of the vegetative stage. The results indicate that the developed low-cost 3D imaging system has the potential for accurate measurements in plant architecture and dry biomass that may be used to improve the accuracy of plant phenotyping.

Background Salinity is an abiotic stress that negatively affects soybean [Glycine max (L.) Merr.] seed yield. Although a major gene for salt tolerance was identified and consistently mapped to chromosome (Chr.) 3 by linkage mapping studies, it does not fully explain genetic variability for tolerance in soybean germplasm. In this study, a genome-wide association study (GWAS) was performed to map genomic regions for salt tolerance in a diverse panel of 305 soybean accessions using a single nucleotide polymorphism (SNP) dataset derived from the SoySNP50K iSelect BeadChip. A second GWAS was also conducted in a subset of 234 accessions using another 3.7 M SNP dataset derived from a whole-genome resequencing (WGRS) study. In addition, three gene-based markers (GBM) of the known gene, Glyma03g32900, on Chr. 3 were also integrated into the two datasets. Salt tolerance among soybean lines was evaluated by leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC). Results For both association studies, a major locus for salt tolerance on Chr. 3 was confirmed by a number of significant SNPs, of which three gene-based SNP markers, Salt-20, Salt14056 and Salt11655, had the highest association with all four traits studied. Also, additional genomic regions on Chrs. 1, 8, and 18 were found to be associated with various traits measured in the second GWAS using the WGRS-derived SNP dataset. Conclusions A region identified on Chr. 8 was identified to be associated with all four traits and predicted as a new minor locus for salt tolerance in soybean. The candidate genes harbored in this minor locus may help reveal the molecular mechanism involved in salt tolerance and to improve tolerance in soybean cultivars. The significant SNPs will be useful for marker-assisted selection for salt tolerance in soybean breeding programs.

Background Brassinosteroids (BRs) play a crucial role in plant vegetative growth and reproductive development. The transcription factors BZR1 and BES1/BZR2 are well characterized as downstream regulators of the BR signaling pathway in Arabidopsis and rice. Soybean contains four BZR1-like proteins (GmBZLs), and it was reported that GmBZL2 plays a conserved role in BR signaling regulation. However, the roles of other GmBZLs have not been thoroughly studied, and the targets of GmBZLs in soybean remain unclear. Results In this study, we first characterized GmBZL3 in soybean from gene expression patterns, conserved domains in coding sequences, and genomic replication times of four GmBZL orthologous. The results indicated that GmBZL3 might play conserved roles during soybean development. The overexpression of GmBZL3P219L in the Arabidopsis BR-insensitive mutant bri1–5 partially rescued the phenotypic defects including BR-insensitivity, which provides further evidence that GmBZL3 functions are conserved between soybean and the homologous Arabidopsis genes. In addition, the identification of the GmBZL3 target genes through ChIP-seq technology revealed that BR has broad roles in soybean and regulates multiple pathways, including other hormone signaling, disease-related, and immunity response pathways. Moreover, the BR-regulated GmBZL3 target genes were further identified, and the results demonstrate that GmBZL3 is a major transcription factor responsible for BR-regulated gene expression and soybean growth. A comparison of GmBZL3 and AtBZR1/BES1 targets demonstrated that GmBZL3 might play conserved as well as specific roles in the soybean BR signaling network. Finally, the identification of two natural soybean varieties of the GmBZL3 mutantion by SNP analysis could facilitate the understanding of gene function during soybean development in the future. Conclusions We illustrate here that GmBZL3 orchestrates a genome-wide transcriptional response that underlies BR-mediated soybean early vegetative growth, and our results support that BRs play crucial regulatory roles in soybean morphology and gene expression levels. Electronic supplementary material The online version of this article (10.1186/s12870-019-1677-2) contains supplementary material, which is available to authorized users.

In soybeans, eighteen members constitute the serine hydroxymethyltransferase (GmSHMT) gene family, of which the cytosolic-targeted GmSHMT08c member has been reported to mediate resistance to soybean cyst nematode (SCN). This work presents a comprehensive study of the SHMT gene family members, including synteny, phylogeny, subcellular localizations, haplotypes, protein homology modeling, mutational, and expression analyses. Phylogenetic analysis showed that SHMT genes are divided into four classes reflecting their subcellular distribution (cytosol, nucleus, mitochondrion, and chloroplast). Subcellular localization of selected GmSHMT members supports their in-silico predictions and phylogenetic distribution. Expression and functional analyses showed that GmSHMT genes display many overlapping, but some divergent responses during SCN infection. Furthermore, mutational analysis reveals that all isolated EMS mutants that lose their resistance to SCN carry missense and nonsense mutations at the GmSHMT08c, but none of the Gmshmt08c mutants carried mutations in the other GmSHMT genes. Haplotype clustering analysis using the whole genome resequencing data from a collection of 106 diverse soybean germplams (15X) was performed to identify allelic variants and haplotypes within the GmSHMT gene family. Interestingly, only the cytosolic-localized GmSHMT08c presented SNP clusters that were associated with SCN resistance, supporting our mutational analysis. Although eight GmSHMT members respond to the nematode infestation, functional and mutational analysis has shown the absence of functional redundancy in resistance to SCN. Structural analysis and protein homology modeling showed the presence of spontaneous mutations at important residues within the GmSHMT proteins, suggesting the presence of altered enzyme activities based on substrate affinities. Due to the accumulation of mutations during the evolution of the soybean genome, the other GmSHMT members have undergone neofunctionalization and subfunctionalization events.

Slow canopy wilting (SW) is a water conservation trait controlled by quantitative trait loci (QTLs) in late maturity group soybeans [Glycine max (L.) Merr.]. Recently, two exotic (landraces) plant introductions (PI 567690 and PI 567731) were identified as new SW lines in early maturity groups. Here, we show that the two PIs share the same water conservation strategy of limited maximum transpiration rates as PI 416937. However, in contrast to PI 416937, the transpiration rates of these PIs were sensitive to an aquaporin inhibitor, indicating an independence between limited maximum transpiration and the lack of silver-sensitive aquaporins. Yield tests of selected recombinant inbred lines from two elite/exotic crosses provide direct evidence to support the benefit of SW in drought tolerance. Four SW QTLs mapped in a Pana×PI 567690 cross at multiple environments were found to be co-located with previous reports. Moreover, two new SW QTLs were mapped on chromosomes 6 and 10 from a Magellan×PI 567731 cross. These two QTLs explain the observed relatively large contributions of 20-30% and were confirmed in a near-isogenic background. These findings demonstrate the importance of SW in yield protection under drought and provide genetic resources for improving drought tolerance in early maturity group soybeans.

The superior agronomic and human nutritional properties of grain legumes (pulses) make them an ideal foundation for future sustainable agriculture. Legume-based farming is particularly important in Africa, where small-scale agricultural systems dominate the food production landscape. Legumes provide an inexpensive source of protein and nutrients to African households as well as natural fertilization for the soil. Although the consumption of traditionally grown legumes has started to decline, the production of soybeans (Glycine max Merr.) is spreading fast, especially across southern Africa. Predictions of future land-use allocation and production show that the soybean is poised to dominate future production across Africa. Land use models project an expansion of harvest area, whereas crop models project possible yield increases. Moreover, a seed change in farming strategy is underway. This is being driven largely by the combined cash crop value of products such as oils and the high nutritional benefits of soybean as an animal feed. Intensification of soybean production has the potential to reduce the dependence of Africa on soybean imports. However, a successful "soybean bonanza" across Africa necessitates an intensive research, development, extension, and policy agenda to ensure that soybean genetic improvements and production technology meet future demands for sustainable production.

Xyloglucan endotransglycosylases/hydrolases (XTHs) are a class of enzymes involved in the construction and remodeling of cellulose/xyloglucan crosslinks and play an important role in regulating cell wall extensibility. However, little is known about this class of enzymes in soybean. Here, 61 soybean XTH genes (GmXTHs) were identified and classified into three subgroups through comparative phylogenetic analysis. Genome duplication greatly contributed to the expansion of GmXTH genes in soybean. A conserved amino acid motif responsible for the catalytic activity was identified in all GmXTHs. Further expression analysis revealed that most GmXTHs exhibited a distinct organ-specific expression pattern, and the expression level of many GmXTH genes was significantly associated with ethylene and flooding stress. To illustrate a possible role of XTH genes in regulating stress responses, the Arabidopsis AtXTH31 gene was overexpressed in soybean. The generated transgenic plants exhibited improved tolerance to flooding stress, with a higher germination rate and longer roots/hypocotyls during the seedling stage and vegetative growth stages. In summary, our combined bioinformatics and gene expression pattern analyses suggest that GmXTH genes play a role in regulating soybean stress responses. The enhanced soybean flooding tolerance resulting from the expression of an Arabidopsis XTH also supports the role of XTH genes in regulating plant flooding stress responses.

The objectives were to (a) quantify the effects of high daytime temperature (HDT) from gametogenesis to full bloom on photosynthesis and pod set in soybean (Glycine max L. Merril) genotypes and (b) assess the relationships among photosynthesis, cardinal temperatures for pollen germination, in vitro pollen germination percentage, canopy reflectance, and pod-set percentage. Three field experiments were conducted, and Experiment I had HDT between gametogenesis and full bloom (36.5°C to 38.6°C) compared with Experiments II and III (29.5°C to 31.6°C; optimum temperature). HDT decreased photosynthesis (22%) and pod-set percent (11%) compared with Experiment III. Cultivars had higher photosynthesis and pod-set percent than plant introduction (PI) lines. The cultivars (i.e., IA3023 and KS4694) and PI lines (i.e., PI393540 and PI588026A) were HDT tolerant and susceptible, respectively. The decreased pod-set percentage in susceptible genotypes (PI lines) was associated with pollen characteristics. Significant positive (r2 ≥ 0.67) association between photosynthesis, cardinal temperatures for pollen germination (Topt and Tmax ) with pod-set percentage was observed. However, a negative (r2 ≥ -0.43) association between photosynthesis and pod set with canopy reflectance at visible spectrum was observed. In vitro pollen germination and canopy reflectance at visible spectrum can be used as a high-throughput phenotypic tool for breeding HDT-tolerant genotypes.

Natural genetic variations in waterlogging tolerance are controlled by multiple genes mapped as quantitative trait loci (QTLs) in major crops, including soybean (Glycine max L.). In this research, 2 novel QTLs associated with waterlogging tolerance were mapped from an elite/exotic soybean cross. The subsequent research was focused on a major QTL (qWT_Gm03) with the tolerant allele from the exotic parent. This QTL was isolated into near-isogenic backgrounds, and its effects on waterlogging tolerance were validated in multiple environments. Fine mapping narrowed qWT_Gm03 into a genomic region of380 Kbp excluding Rps1 gene for Phytophthora sojae resistance. The tolerant allele of qWT_Gm03 promotes root growth under nonstress conditions and favourable root plasticity under waterlogging, resulting in improved waterlogging tolerance, yield, and drought tolerance-related traits, possibly through more efficient water/nutrient uptakes. Meanwhile, involvement of auxin pathways was also identified in the regulation of waterlogging tolerance, as the genotypic differences of qWT_Gm03 in waterlogging tolerance and formation of adventitious/aerial roots can be complemented by an exogenous auxin-biosynthesis inhibitor. These findings provided genetic resources to address the urgent demand of improving waterlogging tolerance in soybean and revealed the determinant roles of root architecture and plasticity in the plant adaptation to waterlogging.

Circadian clocks are a great evolutionary innovation and provide competitive advantage during the day/night cycle and under changing environmental conditions. The circadian clock mediates expression of a large proportion of genes in plants, achieving a harmonious relationship between energy metabolism, photosynthesis, and biotic and abiotic stress responses. Here it is shown that multiple paralogues of clock genes are present in soybean (Glycine max) and mediate flooding and drought responses. Differential expression of many clock and SUB1 genes was found under flooding and drought conditions. Furthermore, natural variation in the amplitude and phase shifts in PRR7 and TOC1 genes was also discovered under drought and flooding conditions, respectively. PRR3 exhibited flooding- and drought-specific splicing patterns and may work in concert with PRR7 and TOC1 to achieve energy homeostasis under flooding and drought conditions. Higher expression of TOC1 also coincides with elevated levels of abscisic acid (ABA) and variation in glucose levels in the morning and afternoon, indicating that this response to abiotic stress is mediated by ABA, endogenous sugar levels, and the circadian clock to fine-tune photosynthesis and energy utilization under stress conditions. It is proposed that the presence of multiple clock gene paralogues with variation in DNA sequence, phase, and period could be used to screen exotic germplasm to find sources for drought and flooding tolerance. Furthermore, fine tuning of multiple clock gene paralogues (via a genetic engineering approach) should also facilitate the development of flooding- and drought-tolerant soybean varieties.

Key message We identified and characterized a mutant of soybean stachyose synthase gene controlling reduced stachyose content which benefit the soybean seed composition breeding program in the future. Abstract It has been shown that in soybean, increased sucrose and reduced raffinose family oligosaccharides would have a positive impact on the world’s feed industry by improving digestibility and feed efficiency. We searched for new sources of modified oligosaccharide content in a subset of the USDA Soybean Germplasm Collection and then identified plant introduction (PI) 603176A as having ultra-low stachyose content (0.5 %). We identified a 33-bp deletion mutant in the putative stachyose synthase gene (STS gene, Glyma19g40550) of PI 603176A. A co-dominate indel marker was successfully developed from this 33-bp deletion area and was genetically mapped into two F2:3 populations and a F4:5 population, which associated with low stachyose content in the progeny lines. These observations provided strong evidence that the STS gene is responsible for stachyose biosynthesis in the soybean plant. Expression of the sts gene remained at the normal level, suggesting the loss of function in the gene is due to defective protein function. This gene-based perfect genetic marker for low stachyose content can be useful for marker-assisted selection in soybean molecular breeding programs. Electronic supplementary material The online version of this article (doi:10.1007/s00122-015-2575-0) contains supplementary material, which is available to authorized users.

Many sources of drought and flooding tolerance have been identified in soybean, however underlying molecular and physiological mechanisms are poorly understood. Therefore, it is important to illuminate different plant responses to these abiotic stresses and understand the mechanisms that confer tolerance. Towards this goal we used four contrasting soybean (Glycine max) genotypes (PI 567690 - drought tolerant, Pana - drought susceptible, PI 408105A - flooding tolerant, S99-2281 - flooding susceptible) grown under greenhouse conditions and compared genotypic responses to drought and flooding at the physiological, biochemical, and cellular level. We also quantified these variations and tried to infer their role in drought and flooding tolerance in soybean. Our results revealed that different mechanisms contribute to reduction in net photosynthesis under drought and flooding stress. Under drought stress, ABA and stomatal conductance are responsible for reduced photosynthetic rate; while under flooding stress, accumulation of starch granules played a major role. Drought tolerant genotypes PI 567690 and PI 408105A had higher plastoglobule numbers than the susceptible Pana and S99-2281. Drought stress increased the number and size of plastoglobules in most of the genotypes pointing to a possible role in stress tolerance. Interestingly, there were seven fibrillin proteins localized within the plastoglobules that were up-regulated in the drought and flooding tolerant genotypes PI 567690 and PI 408105A, respectively, but down-regulated in the drought susceptible genotype Pana. These results suggest a potential role of Fibrillin proteins, FBN1a, 1b and 7a in soybean response to drought and flooding stress.

Nuclear factor-Y (NF-Y), a heterotrimeric transcription factor, is composed of NF-YA, NF-YB and NF-YC proteins. In plants, there are usually more than 10 genes for each family and their members have been identified to be key regulators in many developmental and physiological processes controlling gametogenesis, embryogenesis, nodule development, seed development, abscisic acid (ABA) signaling, flowering time, primary root elongation, blue light responses, endoplasmic reticulum (ER) stress response and drought tolerance. Taking the advantages of the recent soybean genome draft and information on functional characterizations of nuclear factor Y (NF-Y) transcription factor family in plants, we identified 21 GmNF-YA, 32 GmNF-YB, and 15 GmNF-YC genes in the soybean (Glycine max) genome. Phylogenetic analyses show that soybean’s proteins share strong homology to Arabidopsis and many of them are closely related to functionally characterized NF-Y in plants. Expression analysis in various tissues of flower, leaf, root, seeds of different developmental stages, root hairs under rhizobium inoculation, and drought-treated roots and leaves revealed that certain groups of soybean NF-Y are likely involved in specific developmental and stress responses. This study provides extensive evaluation of the soybean NF-Y family and is particularly useful for further functional characterization of GmNF-Y proteins in seed development, nodulation and drought adaptation of soybean. Electronic supplementary material The online version of this article (doi:10.1007/s00438-014-0978-2) contains supplementary material, which is available to authorized users.

The confirmation of a major locus associated with salt tolerance and mapping of a new locus, which could be beneficial for improving salt tolerance in soybean. Breeding soybean for tolerance to high salt conditions is important in some regions of the USA and world. Soybean cultivar Fiskeby III (PI 438471) in maturity group 000 has been reported to be highly tolerant to multiple abiotic stress conditions, including salinity. In this study, a mapping population of 132 F2 families derived from a cross of cultivar Williams 82 (PI 518671, moderately salt sensitive) and Fiskeby III (salt tolerant) was analyzed to map salt tolerance genes. The evaluation for salt tolerance was performed by analyzing leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC) after treatment with 120 mM NaCl under greenhouse conditions. Genotypic data for the F2 population were obtained using the SoySNP6K Illumina Infinium BeadChip assay. A major allele from Fiskeby III was significantly associated with LSS, CCR, LSC, and LCC on chromosome (Chr.) 03 with LOD scores of 19.1, 11.0, 7.7 and 25.6, respectively. In addition, a second locus associated with salt tolerance for LSC was detected and mapped on Chr. 13 with an LOD score of 4.6 and an R 2 of 0.115. Three gene-based polymorphic molecular markers (Salt-20, Salt14056 and Salt11655) on Chr.03 showed a strong predictive association with phenotypic salt tolerance in the present mapping population. These molecular markers will be useful for marker-assisted selection to improve salt tolerance in soybean.

Key message Genetic improvement of soybean protein meal is a complex process because of negative correlation with oil, yield, and temperature. This review describes the progress in mapping and genomics, identifies knowledge gaps, and highlights the need of integrated approaches. Abstract Meal protein derived from soybean [Glycine max (L) Merr.] seed is the primary source of protein in poultry and livestock feed. Protein is a key factor that determines the nutritional and economical value of soybean. Genetic improvement of soybean seed protein content is highly desirable, and major quantitative trait loci (QTL) for soybean protein have been detected and repeatedly mapped on chromosomes (Chr.) 20 (LG-I), and 15 (LG-E). However, practical breeding progress is challenging because of seed protein content’s negative genetic correlation with seed yield, other seed components such as oil and sucrose, and interaction with environmental effects such as temperature during seed development. In this review, we discuss rate-limiting factors related to soybean protein content and nutritional quality, and potential control factors regulating seed storage protein. In addition, we describe advances in next-generation sequencing technologies for precise detection of natural variants and their integration with conventional and high-throughput genotyping technologies. A syntenic analysis of QTL on Chr. 15 and 20 was performed. Finally, we discuss comprehensive approaches for integrating protein and amino acid QTL, genome-wide association studies, whole-genome resequencing, and transcriptome data to accelerate identification of genomic hot spots for allele introgression and soybean meal protein improvement. Electronic supplementary material The online version of this article (doi:10.1007/s00122-017-2955-8) contains supplementary material, which is available to authorized users.

Background Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling. Methods Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips. Results Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others. Conclusions These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water. Electronic supplementary material The online version of this article (doi:10.1186/s12918-016-0374-0) contains supplementary material, which is available to authorized users.

Much research has been conducted on the changes in gene expression of the model plant Arabidopsis to low-oxygen stress. Flooding results in a low oxygen environment in the root zone. However, there is ample evidence that tolerance to soil flooding is more than tolerance to low oxygen alone. In this study, we investigated the physiological response and differential expression of root-related transcription factors (TFs) associated with the tolerance of soybean plants to soil flooding. Differential responses of PI408105A and S99-2281 plants to ten days of soil flooding were evaluated at physiological, morphological and anatomical levels. Gene expression underlying the tolerance response was investigated using qRT-PCR of root-related TFs, known anaerobic genes, and housekeeping genes. Biomass of flood-sensitive S99-2281 roots remained unchanged during the entire 10 days of flooding. Flood-tolerant PI408105A plants exhibited recovery of root growth after 3 days of flooding. Flooding induced the development of aerenchyma and adventitious roots more rapidly in the flood-tolerant than the flood-sensitive genotype. Roots of tolerant plants also contained more ATP than roots of sensitive plants at the 7th and 10th days of flooding. Quantitative transcript analysis identified 132 genes differentially expressed between the two genotypes at one or more time points of flooding. Expression of genes related to the ethylene biosynthesis pathway and formation of adventitious roots was induced earlier and to higher levels in roots of the flood-tolerant genotype. Three potential flood-tolerance TFs which were differentially expressed between the two genotypes during the entire 10-day flooding duration were identified. This study confirmed the expression of anaerobic genes in response to soil flooding. Additionally, the differential expression of TFs associated with soil flooding tolerance was not qualitative but quantitative and temporal. Functional analyses of these genes will be necessary to reveal their potential to enhance flooding tolerance of soybean cultivars.

The QTLs controlling alpha-linolenic acid concentration from wild soybean were mapped on nine soybean chromosomes with various phenotypic variations. New QTLs for alpha-linolenic acid were detected in wild soybean. Alpha-linolenic acid (ALA) is a polyunsaturated fatty acid desired in human and animal diets. Some wild soybean (Glycine soja) genotypes are high in ALA. The objective of this study was to identify quantitative trait loci (QTLs) controlling ALA concentration in a wild soybean accession, PI483463. In total, 188 recombinant inbred lines of F5:6, F5:7, and F5:8 generations derived from a cross of wild soybean PI483463 (~15 % ALA) and cultivar Hutcheson (~9 % ALA) were planted in four environments. Harvested seeds were used to measure fatty acid concentration. Single nucleotide polymorphism markers of the universal soybean linkage panel (USLP 1.0) and simple sequence repeat markers were used for molecular genotyping. Nine putative QTLs were identified that controlled ALA concentration by model-based composite interval mapping and mapped to different soybean chromosomes. The QTLs detected in four environments explained 2.4-7.9 % of the total phenotypic variation (PV). Five QTLs, qALA5_3, qALA6_1, qALA14_1, qALA15_1, and qALA17_1, located on chromosomes 5, 6, 14, 15, and 17 were identified by model-based composite interval mapping and composite interval mapping in two individual environments. Among them, qALA6_1 showed the highest contribution to the PV with 10.0-10.2 % in two environments. The total detected QTLs for additive and epistatic effects explained 52.4 % of the PV for ALA concentration. These findings will provide useful information for understanding genetic structure and marker-assisted breeding programs to increase ALA concentration in seeds derived from wild soybean PI483463.

The reaction catalyzed by squalene synthase (EC.2.5.1.21) that converts two molecules of farnesyl pyrophosphate to squalene represents a crucial branch point of the isoprenoid pathway in diverting carbon flux towards the biosynthesis of sterols. In the present study two soybean squalene synthase genes, GmSQS1 and GmSQS2, were identified in the soybean genome and functionally characterized for their roles in sterol biosynthesis. Both genes encode a deduced protein of 413 amino acids. Complementation assays showed that the two genes were able to convert yeast sterol auxotrophy erg9 mutant to sterol prototrophy. Expression of GmSQS1 and GmSQS2 was ubiquitous in roots, stem, leaves, flower and young seeds of soybean, however GmSQS1 transcript was preferential in roots while GmSQS2 transcript was more in leaves. Their expression was lower in response to dehydration treatments suggesting they might be negative regulators of water stress adaptation. Transgenic Arabidopsis plants overexpressing GmSQS1 driven by either constitutive or seed-specific promoters showed increases in the major end product sterols: campesterol, sitosterol and stigmasterol, which resulted in up to 50% increase in total sterol content in the seeds. The increase in the end product sterols by GmSQS1 overexpression was at the level achievable by previously reported overexpression of individual or combination of other key enzymes in the sterol pathway. Together the data demonstrate that soybean SQS genes play an important role in diverting carbon flux to the biosynthesis of the end product sterols in the seeds.

Soybean Knowledge Base (http://soykb.org) is a comprehensive web resource developed for bridging soybean translational genomics and molecular breeding research. It provides information for six entities including genes/proteins, microRNAs/sRNAs, metabolites, single nucleotide polymorphisms, plant introduction lines and traits. It also incorporates many multi-omics datasets including transcriptomics, proteomics, metabolomics and molecular breeding data, such as quantitative trait loci, traits and germplasm information. Soybean Knowledge Base has a new suite of tools such as In Silico Breeding Program for soybean breeding, which includes a graphical chromosome visualizer for ease of navigation. It integrates quantitative trait loci, traits and germplasm information along with genomic variation data, such as single nucleotide polymorphisms, insertions, deletions and genome-wide association studies data, from multiple soybean cultivars and Glycine soja.

The objective of this study was to use next-generation sequencing technologies to dissect quantitative trait loci (QTL) for southern root-knot nematode (RKN) resistance into individual genes in soybean. Two hundred forty-six recombinant inbred lines (RIL) derived from a cross between Magellan (susceptible) and PI 438489B (resistant) were evaluated for RKN resistance in a greenhouse and sequenced at an average of 0.19× depth. A sequence analysis pipeline was developed to identify and validate single-nucleotide polymorphisms (SNPs), infer the parental source of each SNP allele, and genotype the RIL population. Based on 109,273 phased SNPs, recombination events in RILs were identified, and a total of 3,509 bins and 3,489 recombination intervals were defined. About 50.8% of bins contain 1 to 10 genes. A linkage map was subsequently constructed by using bins as molecular markers. Three QTL for RKN resistance were identified. Of these, one major QTL was mapped to bin 10 of chromosome 10, which is 29.7 kb in size and harbors three true genes and two pseudogenes. Based on sequence variations and gene-expression analysis, the candidate genes underlying the major QTL for RKN resistance were pinpointed. They are Glyma10g02150 and Glyma10g02160, encoding a pectin methylesterase inhibitor and a pectin methylesterase inhibitor -pectin methylesterase, respectively. This QTL mapping approach not only combines SNP discovery, SNP validation, and genotyping, but also solves the issues caused by genome duplication and repetitive sequences. Hence, it can be widely used in crops with a reference genome to enhance QTL mapping accuracy.

Background With the advances in next-generation sequencing (NGS) technology and significant reductions in sequencing costs, it is now possible to sequence large collections of germplasm in crops for detecting genome-scale genetic variations and to apply the knowledge towards improvements in traits. To efficiently facilitate large-scale NGS resequencing data analysis of genomic variations, we have developed “PGen”, an integrated and optimized workflow using the Extreme Science and Engineering Discovery Environment (XSEDE) high-performance computing (HPC) virtual system, iPlant cloud data storage resources and Pegasus workflow management system (Pegasus-WMS). The workflow allows users to identify single nucleotide polymorphisms (SNPs) and insertion-deletions (indels), perform SNP annotations and conduct copy number variation analyses on multiple resequencing datasets in a user-friendly and seamless way. Results We have developed both a Linux version in GitHub (https://github.com/pegasus-isi/PGen-GenomicVariations-Workflow) and a web-based implementation of the PGen workflow integrated within the Soybean Knowledge Base (SoyKB), (http://soykb.org/Pegasus/index.php). Using PGen, we identified 10,218,140 single-nucleotide polymorphisms (SNPs) and 1,398,982 indels from analysis of 106 soybean lines sequenced at 15X coverage. 297,245 non-synonymous SNPs and 3330 copy number variation (CNV) regions were identified from this analysis. SNPs identified using PGen from additional soybean resequencing projects adding to 500+ soybean germplasm lines in total have been integrated. These SNPs are being utilized for trait improvement using genotype to phenotype prediction approaches developed in-house. In order to browse and access NGS data easily, we have also developed an NGS resequencing data browser (http://soykb.org/NGS_Resequence/NGS_index.php) within SoyKB to provide easy access to SNP and downstream analysis results for soybean researchers. Conclusion PGen workflow has been optimized for the most efficient analysis of soybean data using thorough testing and validation. This research serves as an example of best practices for development of genomics data analysis workflows by integrating remote HPC resources and efficient data management with ease of use for biological users. PGen workflow can also be easily customized for analysis of data in other species. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1227-y) contains supplementary material, which is available to authorized users.

Integration of genetic analysis, molecular biology, and genomic approaches drastically enhanced our understanding of genetic control of nematode resistance and provided effective breeding strategies in soybeans. Three nematode species, including soybean cyst (SCN, Heterodera glycine), root-knot (RKN, Meloidogyne incognita), and reniform (RN, Rotylenchulus reniformis), are the most destructive pests and have spread to soybean growing areas worldwide. Host plant resistance has played an important role in their control. This review focuses on genetic, genomic studies, and breeding efforts over the past two decades to identify and improve host resistance to these three nematode species. Advancements in genetics, genomics, and bioinformatics have improved our understanding of the molecular and genetic mechanisms of nematode resistance and enabled researchers to generate large-scale genomic resources and marker-trait associations. Whole-genome resequencing, genotyping-by-sequencing, genome-wide association studies, and haplotype analyses have been employed to map and dissect genomic locations for nematode resistance. Recently, two major SCN-resistant loci, Rhg1 and Rhg4, were cloned and other novel resistance quantitative trait loci (QTL) have been discovered. Based on these discoveries, gene-specific DNA markers have been developed for both Rhg1 and Rhg4 loci, which were useful for marker-assisted selection. With RKN resistance QTL being mapped, candidate genes responsible for RKN resistance were identified, leading to the development of functional single nucleotide polymorphism markers. So far, three resistances QTL have been genetically mapped for RN resistance. With nematode species overcoming the host plant resistance, continuous efforts in the identification and deployment of new resistance genes are required to support the development of soybean cultivars with multiple and durable resistance to these pests.