Search

'S17-2243C' (Reg. no. CV-557, PI 700003) is a semi-determinate, maturity group IV (relative maturity 4.9), non-genetically modified (non-GM) soybean [Glycine max (L.) Merr.] cultivar developed and released by the University of Missouri-Fisher Delta Research, Extension, and Education Center. S17-2243C was developed to meet the growing demands for new non-GM soybean cultivars with high yield and elevated seed oil content. S17-2243C is resistant to stem canker and charcoal rot and has tolerance to salinity conditions. Seed of S17-2243C has averaged 232 g kg(-1) of oil concentration on a dry weight basis, which was significantly higher than all check cultivars in the 2020 USDA Uniform Soybean Tests, Southern States. S17-2243C was tested against high-yielding private and public soybean cultivars from 2018 to 2021 in 80 locations across 12 states, including Alabama, Arkansas, Illinois, Kentucky, Louisiana, Mississippi, Missouri, North Carolina, Ohio, South Carolina, Tennessee, and Virginia. With high yield potential, broad adaptability, early maturity, elevated seed oil content, and non-GM traits, S17-2243C is an excellent cultivar choice for soybean growers adopting alternative growing systems and benefiting from premium prices offered for non-GM soybean products. Missouri Soybean Merchandising Council [6066-21220-014-000D]; USDA-ARS Project; [301] Published version Missouri Soybean Merchandising Council, Grant/Award Number: 301; USDA-ARS Project, Grant/Award Number: 6066-21220-014-000D Public domain – authored by a U.S. government employee

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+, we identified an allelic series of mutations in five soybean acyl-ACP thioesterase genes, including GmFATA1A, GmFATB1A, GmFATB1B, GmFATB2A, and GmFATB2B. Additionally, we discovered mutations at GmFATA1A to confer high oleic acid (up to 34.5%) content, while mutations at GmFATB presented low palmitic acid (as low as 5.6%) and high oleic acid (up to 36.5%) phenotypes. The obtained soybean mutants with altered fatty acid content can be used in soybean breeding program for improving soybean oil composition traits.

Genetic analysis identified a unique combination of major QTL for resistance to important soybean nematodes concurrently present in a single soybean accession, which has not been reported earlier. An exotic soybean [Glycine max (L.) Merr.] accession, PI 567305, was reported to be highly resistant to three important nematode species, soybean cyst (SCN), root-knot (RKN), and reniform (RN) nematodes. However, genetic basis controlling broad-spectrum resistance in this germplasm has not been investigated. We report results of genetic analysis to identify genomic loci conferring resistance to these nematode species. A bi-parental population consisting of 242 F8-derived recombinant inbred lines (RILs) was developed from a cross of a nematode susceptible cultivar, Magellan, and resistant accession, PI 567305. The RILs were phenotyped for nematode resistance to three SCN HG types. They were genotyped using the Infinium SoySNP6K BeadChips and genotype-by-sequencing (GBS) methods in an attempt to evaluate the cost-effectiveness and efficiency of these two genotyping platforms. Genetic analysis confirmed the major QTL on chromosomes (Chrs) 10 and 18 with broad-spectrum resistance to the three nematodes present in this germplasm. Haplotype and copy number variation analyses of SCN resistance QTL indicated that PI 567305 has a different haplotype, which is associated with likely a unique SCN resistance mechanism different from Peking- or PI 88788-type resistance. The evaluations of both Infinium Beadchip- and GBS-based genotyping technologies provided comprehensive insights for researchers to choose a cost-effective and efficient platform for QTL mapping and for other genomic studies in soybeans.

Phosphorus (P) is an essential macronutrient for all organisms. Phosphate (Pi) deficiency reduces grain yield and quality in wheat. Understanding how wheat responds to Pi deficiency at the global transcriptional level remains limited. We revisited the available RNA-seq transcriptome from Pi-starved wheat roots and shoots subjected to Pi starvation. Genome-wide transcriptome resetting was observed under Pi starvation, with a total of 917 and 2338 genes being differentially expressed in roots and shoots, respectively. Chromosomal distribution analysis of the gene triplets and differentially expressed genes (DEGs) revealed that the D genome displayed genome induction bias and, specifically, the chromosome 2D might be a key contributor to Pi-limiting triggered gene expression response. Alterations in multiple metabolic pathways pertaining to secondary metabolites, transcription factors and Pi uptake-related genes were evidenced. This study provides genomic insight and the dynamic landscape of the transcriptional changes contributing to the hexaploid wheat during Pi starvation. The outcomes of this study and the follow-up experiments have the potential to assist the development of Pi-efficient wheat cultivars.

There is a decrease in the land available for rice cultivation due to the rapid conversion to urban uses. Subsequently, acid soil could be an alternative land cultivating rice, but will require the use of aluminum (Al)-tolerant rice varieties. This Al tolerance trait is genetically controlled, and there is a need to discover more genes needed to develop Al-tolerant rice. Therefore, the objective of this study was to clone and characterize a novel Al tolerance gene isolated from a local cultivar of Indonesian rice. The gene cloning was conducted based on the rye/rice microsynteny relationship. In addition, the root growth and gene expression analyses were performed to verify the role of the gene on Al tolerance in gene-silenced rice and in overexpressed transgenic tobacco. The results showed an Al tolerance candidate gene, OsGERLP, was successfully cloned from rice cv. Hawara Bunar, with its gene encoding a protein similar to a bacterial ribosomal L32 protein. Additionally, the analysis showed that low gene expression caused the gene-silenced rice to be sensitive to Al, while high expression induced the Al tolerance in transgenic tobacco. Furthermore, it was discovered that the gene expression level in both plants was in line with the lower expression of the OsFRDL4 gene in the silenced rice and the high expression of the MATE gene in transgenic tobacco also with the higher citrate secretion from transgenic tobacco roots. In conclusion, the OsGERLP gene could act as a regulator for other Al tolerance genes, with the potential to develop Al-tolerant rice varieties.

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.

Epigenetics is defined as changes in gene expression that are not associated with changes in DNA sequence but due to the result of methylation of DNA and post-translational modifications to the histones. These epigenetic modifications are known to regulate gene expression by bringing changes in the chromatin state, which underlies plant development and shapes phenotypic plasticity in responses to the environment and internal cues. This review articulates the role of histone modifications and DNA methylation in modulating biotic and abiotic stresses, as well as crop improvement. It also highlights the possibility of engineering epigenomes and epigenome-based predictive models for improving agronomic traits.

Key message This review provides a comprehensive atlas of QTLs, genes, and alleles conferring resistance to 28 important diseases in all major soybean production regions in the world. Abstract Breeding disease-resistant soybean [Glycine max (L.) Merr.] varieties is a common goal for soybean breeding programs to ensure the sustainability and growth of soybean production worldwide. However, due to global climate change, soybean breeders are facing strong challenges to defeat diseases. Marker-assisted selection and genomic selection have been demonstrated to be successful methods in quickly integrating vertical resistance or horizontal resistance into improved soybean varieties, where vertical resistance refers to R genes and major effect QTLs, and horizontal resistance is a combination of major and minor effect genes or QTLs. This review summarized more than 800 resistant loci/alleles and their tightly linked markers for 28 soybean diseases worldwide, caused by nematodes, oomycetes, fungi, bacteria, and viruses. The major breakthroughs in the discovery of disease resistance gene atlas of soybean were also emphasized which include: (1) identification and characterization of vertical resistance genes reside rhg1 and Rhg4 for soybean cyst nematode, and exploration of the underlying regulation mechanisms through copy number variation and (2) map-based cloning and characterization of Rps11 conferring resistance to 80% isolates of Phytophthora sojae across the USA. In this review, we also highlight the validated QTLs in overlapping genomic regions from at least two studies and applied a consistent naming nomenclature for these QTLs. Our review provides a comprehensive summary of important resistant genes/QTLs and can be used as a toolbox for soybean improvement. Finally, the summarized genetic knowledge sheds light on future directions of accelerated soybean breeding and translational genomics studies.

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

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.

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306), map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020), identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

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.

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.

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

Soybean is sensitive to flooding stress that may result in poor seed quality and significant yield reduction. Soybean production under flooding could be sustained by developing flood-tolerant cultivars through breeding programs. Conventionally, soybean tolerance to flooding in field conditions is evaluated by visually rating the shoot injury/damage due to flooding stress, which is labor-intensive and subjective to human error. Recent developments of field high-throughput phenotyping technology have shown great potential in measuring crop traits and detecting crop responses to abiotic and biotic stresses. The goal of this study was to investigate the potential in estimating flood-induced soybean injuries using UAV-based image features collected at different flight heights. The flooding injury score (FIS) of 724 soybean breeding plots was taken visually by breeders when soybean showed obvious injury symptoms. Aerial images were taken on the same day using a five-band multispectral and an infrared (IR) thermal camera at 20, 50, and 80 m above ground. Five image features, i.e., canopy temperature, normalized difference vegetation index, canopy area, width, and length, were extracted from the images at three flight heights. A deep learning model was used to classify the soybean breeding plots to five FIS ratings based on the extracted image features. Results show that the image features were significantly different at three flight heights. The best classification performance was obtained by the model developed using image features at 20 m with 0.9 for the five-level FIS. The results indicate that the proposed method is very promising in estimating FIS for soybean breeding.

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

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

Flooding threatens soybean production and limits soybean yields worldwide. The most effective and economic approach to decrease loss of yield due to flooding is to develop flood-tolerant soybean cultivars. The objective of this study was to identify genetic loci and candidate genes associated with flooding tolerance. A panel of 384 soybean plant introductions (PIs) was evaluated for flooding tolerance in two consecutive years in the field. The plant foliar damage score was used to index soybean response to flooding stress. A total of 42,291 SNP markers were obtained from the Illumina Infinium SoySNP50K BeadChip database. After filtration for quality control, 31,125 SNPs were used for genome-wide association mapping utilizing four different models (regression linear model (GLM), mixed linear model (MLM), compressed mixed linear model (CMLM), and enriched compressed mixed linear model (ECMLM)). Fourteen SNPs were identified to be associated with flooding tolerance across all environments and models at a significance level of −Log10 (P) ≥ 2.5. Five SNPs were located within the coding regions of five candidate genes. Several PIs with lower best linear unbiased prediction (BLUPs) of the breeding values and a large number of favorable flood-tolerant alleles were found as new genetic sources for use in soybean breeding programs.

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.

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.