Search

Soybean, one of the most important crops worldwide, is severely affected by abiotic stress. Drought and flooding are the major abiotic stresses impacting soybean yield. In this regard, understanding water uptake by plants, its utilization and transport has great importance. In plants, water transport is mainly governed by channel forming aquaporin proteins (AQPs). Tonoplast intrinsic proteins (TIPs) belong to the plant-specific AQP subfamily and are known to have a role in abiotic stress tolerance. In this study, 23 soybean TIP genes were identified based on the latest soybean genome annotation. TIPs were characterized based on conserved structural features and phylogenetic distribution. Expression analysis of soybean TIP genes in various tissues and under abiotic stress conditions demonstrated tissue/stress-response specific differential expression. The natural variations for TIP genes were analyzed using whole genome re-sequencing data available for a set of 106 diverse soybean genotypes including wild types, landraces and elite lines. Results revealed 81 SNPs and several large insertions/deletions in the coding region of TIPs. Among these, non-synonymous SNPs are most likely to have a greater impact on protein function and are candidates for molecular studies as well as for the development of functional markers to assist breeding. The solute transport function of two TIPs was further validated by expression in Xenopus laevis oocytes. GmTIP1;5 was shown to facilitate the rapid movement of water across the oocyte membrane, while GmTIP2;5 facilitated the movement of water and boric acid. The present study provides an initial insight into the possible roles of soybean TIP genes under abiotic stress conditions. Our results will facilitate elucidation of their precise functions during abiotic stress responses and plant development, and will provide potential breeding targets for modifying water movement in soybean.

Soybean is a major crop in the world and an essential source for minerals, including iron (Fe) and zin (Zn). Deficiency of Fe and Zn in soil and soybean plants result in yield loss and poor seed nutritional qualities. Information on genomic regions and candidate genes controlling seed Fe and Zn accumulation in soybean seed is limited. Therefore, The objective of this research was to identify genetic regions, known as quantitative trait loci (QTL), and candidate genes that control the accumulation of Fe and Zn in soybean mature seeds. A ‘Forrest’ by ‘Williams 82’ (F × W82) recombinant inbred line (RIL) population (n = 306) was used and genotyped using a total of 5405 single nucleotides polymorphism (SNP) markers using Infinium SNP6K BeadChips. A two-year experiment was conducted across two environments: North Carolina in 2018 (NC) and Illinois in 2020 (IL). Only QTL with LOD scores ≥ 2.5, as identified by the composite interval mapping (CIM) method, are reported here. In total, 6 QTL were identified for seed Fe; specifically, 3 QTL (qFe-01-[NC-2018], qFe-02-[NC-2018], and qFe-03-[NC-2018]) were located on chromosomes 1, 2, and 6, respectively, in the NC environment, and 3 QTL (qFe-01-[IL-2020], qFe-02-[IL-2020], and qFe-03-[IL-2020]) were positioned on chromosomes 1, 2, and 12, respectively, in the IL environment. A total of 6 QTL associated with seed Zn were also identified; 4 QTL (qZn-01-[NC-2018]; qZn-02-[NC-2018]; qZn-03-[NC-2018]; and qZn-04-[NC-2018]), respectively on Chr 2, 3, 7, and 19 in NC; and 2 QTL (qZn-01-[IL-2020] and qZn-02-[IL-2020]), respectively, on Chr 5 and 8 in IL. Several functional genes encode Fe- and Zn-proteins, transcription factors, proteins-zinc finger motifs (involved in DNA binding and transcriptional regulation; crosstalk between the regulatory pathways of Zn and Fe transporters) were identified and located within the QTL interval. To our knowledge, and based on the literature available, the QTL identified here on Chr 2 and Chr 6 are novel and were not previously identified. This current research provides a new knowledge of the genetic basis of seed Fe and Zn and the markers associated with QTL. The QTL identified here will contribute to efficient marker assisted selection for higher Fe and Zn content in soybean seeds. The candidate genes and metal-responsive transcription factors may coordinate the expression of both Zn and Fe transporters in response to changes in metal availability, providing new knowledge on minerals uptake and transport mechanisms, allowing for possible genetic engineering application.

Drought and flooding are two major causes of severe yield loss in soybean worldwide. A lack of knowledge of the molecular mechanisms involved in drought and flood stress has been a limiting factor for the effective management of soybeans; therefore, it is imperative to assess the expression of genes involved in response to flood and drought stress. In this study, differentially expressed genes under drought and flooding conditions were investigated using Illumina RNA-Seq transcriptome profiling. A total of 2,724 and 3,498 differentially expressed genes (DEGs) were identified under drought and flooding treatments, respectively. These genes comprise 289 Transcription Factors (TFs) representing Basic Helix-loop Helix (bHLH), Ethylene Response Factors (ERFs), myeloblastosis (MYB), No apical meristem (NAC), and WRKY amino acid motif (WRKY) type major families known to be involved in the mechanism of stress tolerance. The expression of photosynthesis and chlorophyll synthesis related genes were significantly reduced under both types of stresses, which limit the metabolic processes and thus help prolong survival under extreme conditions. However, cell wall synthesis related genes were up-regulated under drought stress and down-regulated under flooding stress. Transcript profiles involved in the starch and sugar metabolism pathways were also affected under both stress conditions. The changes in expression of genes involved in regulating the flux of cell wall precursors and starch/sugar content can serve as an adaptive mechanism for soybean survival under stress conditions. This study has revealed the involvement of TFs, transporters, and photosynthetic genes, and has also given a glimpse of hormonal cross talk under the extreme water regimes, which will aid as an important resource for soybean crop improvement.

The phytohormone auxin plays a critical role in regulation of plant growth and development as well as plant responses to abiotic stresses. This is mainly achieved through its uneven distribution in plants via a polar auxin transport process. Auxin transporters are major players in polar auxin transport. The AUXIN RESISTANT 1 ⁄ LIKE AUX1 (AUX⁄LAX) auxin influx carriers belong to the amino acid permease family of proton-driven transporters and function in the uptake of indole-3-acetic acid (IAA). In this study, genome-wide comprehensive analysis of the soybean AUX⁄LAX (GmLAX) gene family, including phylogenic relationships, chromosome localization, and gene structure, were carried out. A total of 15 GmLAX genes, including seven duplicated gene pairs, were identified in the soybean genome. They were distributed on 10 chromosomes. Despite their higher percentage identities at the protein level, GmLAXs exhibited versatile tissue-specific expression patterns, indicating coordinated functioning during plant growth and development. Most GmLAXs were responsive to drought and dehydration stresses and auxin and abscisic acid (ABA) stimuli, in a tissue- and/or time point- sensitive mode. Several GmLAX members were involved in responding to salt stress. Sequence analysis revealed that promoters of GmLAXs contained different combinations of stress-related cis-regulatory elements. These studies suggest that the soybean GmLAXs were under control of a very complex regulatory network, responding to various internal and external signals. This study helps to identity candidate GmLAXs for further analysis of their roles in soybean development and adaption to adverse environments.

Resistance to soybean rust (SBR), caused by Phakopsora pachyrhizi Syd. & Syd., has been identified in many soybean germplasm accessions and is conferred by either dominant or recessive genes that have been mapped to six independent loci (Rpp1 -Rpp6), but No U.S. cultivars are resistant to SBR. The cultivar DT 2000 (PI 635999) has resistance to P. pachyrhizi isolates and field populations from the United States as well as Vietnam. A F6:7 recombinant inbred line (RIL) population derived from Williams 82 × DT 2000 was used to identify genomic regions associated with resistance to SBR in the field in Ha Noi, Vietnam, and in Quincy, Florida, in 2008. Bulked segregant analysis (BSA) was conducted using the soybean single nucleotide polymorphism (SNP) USLP 1.0 panel along with simple sequence repeat (SSR) markers to detect regions of the genome associated with resistance. BSA identified four BARC_SNP markers near the Rpp3 locus on chromosome (Chr.) 6. Genetic analysis identified an additional genomic region around the Rpp4 locus on Chr. 18 that was significantly associated with variation in the area under disease progress curve (AUDPC) values and sporulation in Vietnam. Molecular markers tightly linked to the DT 2000 resistance alleles on Chrs. 6 and 18 will be useful for marker-assisted selection and backcrossing in order to pyramid these genes with other available SBR resistance genes to develop new varieties with enhanced and durable resistance to SBR.

Food resources of the modern world are strained due to the increasing population. There is an urgent need for innovative methods and approaches to augment food production. Legume seeds are major resources of human food and animal feed with their unique nutrient compositions including oil, protein, carbohydrates, and other beneficial nutrients. Recent advances in next-generation sequencing (NGS) together with ‘omics’ technologies have considerably strengthened soybean research. The availability of well annotated soybean genome sequence along with hundreds of identified quantitative trait loci (QTL) associated with different seed traits can be used for gene discovery and molecular marker development for breeding applications. Despite the remarkable progress in these technologies, the analysis and mining of existing seed genomics data are still challenging due to the complexity of genetic inheritance, metabolic partitioning, and developmental regulations. Integration of ‘omics tools’ is an effective strategy to discover key regulators of various seed traits. In this review, recent advances in ‘omics’ approaches and their use in soybean seed trait investigations are presented along with the available databases and technological platforms and their applicability in the improvement of soybean. This article also highlights the use of modern breeding approaches, such as genome-wide association studies (GWAS), genomic selection (GS), and marker-assisted recurrent selection (MARS) for developing superior cultivars. A catalog of available important resources for major seed composition traits, such as seed oil, protein, carbohydrates, and yield traits are provided to improve the knowledge base and future utilization of this information in the soybean crop improvement programs.

Waterlogged condition due to flooding is one of the major abiotic stresses that drastically affect the soybean growth and yield around the world. Because of this, many breeders have focused on the development of waterlogging tolerance in soybean varieties, and thus, several tolerant varieties were developed. However, the physiological mechanism of waterlogging tolerance is not yet fully understood. We particularly studied the endogenous hormones regulation during waterlogging in two contrasting soybean genotypes. According to our results, adventitious roots better developed in the waterlogging tolerant line (WTL) than in the waterlogging susceptible line (WSL). Endogenous hormones also showed significant differences between WTL and WSL. The ethylene production ratio was higher in WTL than in WSL, and methionine was higher in WTL than in WSL. Other endogenous abscisic acid (ABA) contents were lower in WTL than in WSL. Conversely, gibberellic acid (GA) showed a tendency to be high in WTL, especially the levels of the bioactive GA4. The ratio of total GA and ABA was significantly higher in WTL than in WSL. Anatomical study of the root revealed that aerenchyma cells in the stele were better developed in WTL than in WSL.

Climate change affects agricultural productivity worldwide. Increased prices of food commodities are the initial indication of drastic edible yield loss, which is expected to surge further due to global warming. This situation has compelled plant scientists to develop climate change-resilient crops, which can withstand broad-spectrum stresses such as drought, heat, cold, salinity, flood and submergence, and pests along with increased productivity. Genomics appears to be a promising tool for deciphering the stress responsiveness of crop species with adaptation traits or in wild relatives towards identifying underlying genes, alleles or quantitative trait loci. Molecular breeding approaches have been proven helpful in enhancing the stress adaptation of crop plants, and recent advancement in next-generation sequencing along with high-throughput sequencing and phenotyping platforms have transformed molecular breeding to genomics-assisted breeding (GAB). In view of this, the present review elaborates the progress and prospects of GAB in improving climate change resilience in crop plants towards circumventing global food insecurity.

Abstract Development of high yielding cowpea varieties coupled with good taste and rich in essential minerals can promote consumption and thus nutrition and profitability. The sweet taste of cowpea grain is determined by its sugar content, which comprises mainly sucrose and galacto-oligosaccharides (GOS) including raffinose and stachyose. However, GOS are indigestible and their fermentation in the colon can produce excess intestinal gas, causing undesirable bloating and flatulence. In this study, we aimed to examine variation in grain sugar and mineral concentrations, then map quantitative trait loci (QTLs) and estimate genomic-prediction (GP) accuracies for possible application in breeding. Grain samples were collected from a multi-parent advanced generation intercross (MAGIC) population grown in California during 2016–2017. Grain sugars were assayed using high-performance liquid chromatography. Grain minerals were determined by inductively coupled plasma–optical emission spectrometry and combustion. Considerable variation was observed for sucrose (0.6–6.9%) and stachyose (2.3–8.4%). Major QTLs for sucrose (QSuc.vu-1.1), stachyose (QSta.vu-7.1), copper (QCu.vu-1.1) and manganese (QMn.vu-5.1) were identified. Allelic effects of major sugar QTLs were validated using the MAGIC grain samples grown in West Africa in 2017. GP accuracies for minerals were moderate (0.4–0.58). These findings help guide future breeding efforts to develop mineral-rich cowpea varieties with desirable sugar content.

Cultivated soybean (Glycine max L.) cv. Dunbar (PI 552538) and wild G. soja (PI 326582A) exhibited significant differences in root architecture and root-related traits. In this study, phenotypic variability for root traits among 251 BC2F5 backcross inbred lines (BILs) developed from the cross Dunbar/PI 326582A were identified. The root systems of the parents and BILs were evaluated in controlled environmental conditions using a cone system at seedling stage. The G. max parent Dunbar contributed phenotypically favorable alleles at a major quantitative trait locus on chromosome 8 (Satt315-I locus) that governed root traits (tap root length and lateral root number) and shoot length. This QTL accounted for >10% of the phenotypic variation of both tap root and shoot length. This QTL region was found to control various shoot- and root-related traits across soybean genetic backgrounds. Within the confidence interval of this region, eleven transcription factors (TFs) were identified. Based on RNA sequencing and Affymetrix expression data, key TFs including MYB, AP2-EREBP and bZIP TFs were identified in this QTL interval with high expression in roots and nodules. The backcross inbred lines with different parental allelic combination showed different expression pattern for six transcription factors selected based on their expression pattern in root tissues. It appears that the marker interval Satt315-I locus on chromosome 8 contain an essential QTL contributing to early root and shoot growth in soybean.

Sweetpotato highly produces carotenoids in storage roots. In this study, a cDNA encoding geranylgeranyl phyrophosphate synthase (GGPS), named IbGGPS, was isolated from sweetpotato storage roots. Green fluorescent protein (GFP) was fused to the C-terminus of IbGGPS to obtain an IbGGPS-GFP fusion protein that was transiently expressed in both epidermal cells of onion and leaves of tobacco. Confocal microscopic analysis determined that the IbGGPS-GFP protein was localized to specific areas of the plasma membrane of onion and chloroplasts in tobacco leaves. The coding region of IbGGPS was cloned into a binary vector under the control of 35S promoter and then transformed into Arabidopsis thaliana to obtain transgenic plants. High performance liquid chromatography (HPLC) analysis showed a significant increase of total carotenoids in transgenic plants. The seeds of transgenic and wild-type plants were germinated on an agar medium supplemented with polyethylene glycol (PEG). Transgenic seedlings grew significantly longer roots than wild-type ones did. Further enzymatic analysis showed an increased activity of superoxide dismutase (SOD) in transgenic seedlings. In addition, the level of malondialdehyde (MDA) was reduced in transgenics. qRT-PCR analysis showed altered expressions of several genes involved in the carotenoid biosynthesis in transgenic plants. These data results indicate that IbGGPS is involved in the biosynthesis of carotenoids in sweetpotato storage roots and likely associated with tolerance to osmotic stress.

In Arabidopsis, NAC (NAM, ATAF and CUC) transcription factors have been found to promote lateral root number through the auxin signaling pathway. In the present study, the role of water stress-inducible soybean GmNAC003 and GmNAC004 genes in the enhancement of lateral root development under water deficit conditions was investigated. Both genes were highly expressed in roots, leaves and flowers of soybean and were strongly induced by water stress and moderately induced by a treatment with abscisic acid (ABA). They showed a slight response to treatment with 2,4-dichlorophenoxyacetic acid (2,4-D). The transgenic Arabidopsis plants overexpressing GmNAC004 showed an increase in lateral root number and length under non-stress conditions and maintained higher lateral root number and length under mild water stress conditions compared to the wild-type (WT), while the transgenic plants overexpressing GmNAC003 did not show any response. However, LR development of GmNAC004 transgenic Arabidopsis plants was not enhanced in the water-stressed compared to the well-watered treatment. In the treatment with ABA, LR density of the GmNAC004 transgenic Arabidopsis was less suppressed than that of the WT, suggesting that GmNAC004 counteracts ABA-induced inhibition of lateral root development. In the treatment with 2,4-D, lateral root density was enhanced in both GmNAC004 transgenic Arabidopsis and WT plants but the promotion was higher in the transgenic plants. Conversely, in the treatment with naphthylphthalamic acid (NPA), lateral root density was inhibited and there was no difference in the phenotype of the GmNAC004 transgenic Arabidopsis and WT plants, indicating that auxin is required for the action of GmNAC004. Transcript analysis for a number of known auxin and ABA related genes showed that GmNAC004's role may suppress ABA signaling but promote auxin signaling to increase lateral root development in the Arabidopsis heterologous system.

Quantitative RT-PCR can be a very sensitive and powerful technique for measuring differential gene expression. Changes in gene expression induced by abiotic stresses are complex and multifaceted, which make determining stably expressed genes for data normalization difficult. To identify the most suitable reference genes for abiotic stress studies in soybean, 13 candidate genes collected from literature were evaluated for stability of expression under dehydration, high salinity, cold and ABA (abscisic acid) treatments using delta CT and geNorm approaches. Validation of reference genes indicated that the best reference genes are tissue- and stress-dependent. With respect to dehydration treatment, the Fbox/ABC, Fbox/60s gene pairs were found to have the highest expression stability in the root and shoot tissues of soybean seedlings, respectively. Fbox and 60s genes are the most suitable reference genes across dehydrated root and shoot tissues. Under salt stress the ELF1b/IDE and Fbox/ELF1b are the most stably expressed gene pairs in roots and shoots, respectively, while 60s/Fbox is the best gene pair in both tissues. For studying cold stress in roots or shoots, IDE/60s and Fbox/Act27 are good reference gene pairs, respectively. With regard to gene expression analysis under ABA treatment in either roots, shoots or across these tissues, 60s/ELF1b, ELF1b/Fbox and 60s/ELF1b are the most suitable reference genes, respectively. The expression of ELF1b/60s, 60s/Fbox and 60s/Fbox genes was most stable in roots, shoots and both tissues, respectively, under various stresses studied. Among the genes tested, 60s was found to be the best reference gene in different tissues and under various stress conditions. The highly ranked reference genes identified from this study were proved to be capable of detecting subtle differences in expression rates that otherwise would be missed if a less stable reference gene was used.

BACKGROUND: Synthesis of polyunsaturated fatty acids (PUFAs) in the endoplasmic reticulum of plants typically involves the fatty acid desaturases FAD2 and FAD3, which use cytochrome b(5) (Cb5) as an electron donor. Higher plants are reported to have multiple isoforms of Cb5, in contrast to a single Cb5 in mammals and yeast. Despite the wealth of information available on the roles of FAD2 and FAD3 in PUFA synthesis, information regarding the contributions of various Cb5 isoforms in desaturase-mediated reactions is limited. RESULTS: The present functional characterization of Cb5 polypeptides revealed that all Arabidopsis Cb5 isoforms are not similarly efficient in ω-6 desaturation, as evidenced by significant variation in their product outcomes in yeast-based functional assays. On the other hand, characterization of Cb5 polypeptides of soybean (Glycine max) suggested that similar ω-6 desaturation efficiencies were shared by various isoforms. With regard to ω-3 desaturation, certain Cb5 genes of both Arabidopsis and soybean were shown to facilitate the accumulation of more desaturation products than others when co-expressed with their native FAD3. Additionally, similar trends of differential desaturation product accumulation were also observed with most Cb5 genes of both soybean and Arabidopsis even if co-expressed with non-native FAD3. CONCLUSIONS: The present study reports the first description of the differential nature of the Cb5 genes of higher plants in fatty acid desaturation and further suggests that ω-3/ω-6 desaturation product outcome is determined by the nature of both the Cb5 isoform and the fatty acid desaturases.

The discovery of nucleotide diversity captured as single feature polymorphism (SFP) by using the expression array is a high-throughput and effective method in detecting genome-wide polymorphism. The efficacy of such method was tested in rice, and the results presented in the paper indicate high sensitivity in predicting SFP. The sensitivity of polymorphism detection was further demonstrated by the fact that no biasness was observed in detecting SFP with either single or multiple nucleotide polymorphisms. The high density SFP data that can be generated quite effectively by the current method has promise for high resolution genetic mapping studies, as physical location of features are well-defined on rice genome.

Abstract In addition to the challenge of meeting global demand for food production, there are increasing concerns about food safety and the need to protect consumer health from the negative effects of foodborne allergies. Certain bio‐molecules (usually proteins) present in food can act as allergens that trigger unusual immunological reactions, with potentially life‐threatening consequences. The relentless working lifestyles of the modern era often incorporate poor eating habits that include readymade prepackaged and processed foods, which contain additives such as peanuts, tree nuts, wheat, and soy‐based products, rather than traditional home cooking. Of the predominant allergenic foods (soybean, wheat, fish, peanut, shellfish, tree nuts, eggs, and milk), peanuts (Arachis hypogaea) are the best characterized source of allergens, followed by tree nuts (Juglans regia, Prunus amygdalus, Corylus avellana, Carya illinoinensis, Anacardium occidentale, Pistacia vera, Bertholletia excels), wheat (Triticum aestivum), soybeans (Glycine max), and kidney beans (Phaseolus vulgaris). The prevalence of food allergies has risen significantly in recent years including chance of accidental exposure to such foods. In contrast, the standards of detection, diagnosis, and cure have not kept pace and unfortunately are often suboptimal. In this review, we mainly focus on the prevalence of allergies associated with peanut, tree nuts, wheat, soybean, and kidney bean, highlighting their physiological properties and functions as well as considering research directions for tailoring allergen gene expression. In particular, we discuss how recent advances in molecular breeding, genetic engineering, and genome editing can be used to develop potential low allergen food crops that protect consumer health.

Abstract Soybean [Glycine max (L.) Merr.] is a globally important crop due to its valuable seed composition, versatile feed, food, and industrial end‐uses, and consistent genetic gain. Successful genetic gain in soybean has led to widespread adaptation and increased value for producers, processors, and consumers. Specific focus on the nutritional quality of soybean seed composition for food and feed has further elucidated genetic knowledge and bolstered breeding progress. Seed components are historical and current targets for soybean breeders seeking to improve nutritional quality of soybean. This article reviews genetic and genomic foundations for improvement of nutritionally important traits, such as protein and amino acids, oil and fatty acids, carbohydrates, and specific food‐grade considerations; discusses the application of advanced breeding technology such as CRISPR/Cas9 in creating seed composition variations; and provides future directions and breeding recommendations regarding soybean seed composition traits.

Abstract Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05‐1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Abstract Complete, gapless telomere‐to‐telomere chromosome assemblies are a prerequisite for comprehensively investigating the architecture of complex regions, like centromeres or telomeres and removing uncertainties in the order, spacing, and orientation of genes. Using complementary genomics technologies and assembly algorithms, we developed highly contiguous, nearly gapless, genome assemblies for two economically important soybean [Glycine max (L.) Merr] cultivars (Williams 82 and Lee). The centromeres were distinctly annotated on all the chromosomes of both assemblies. We further found that the canonical telomeric repeats were present at the telomeres of all chromosomes of both Williams 82 and Lee genomes. A total of 10 chromosomes in Williams 82 and eight in Lee were entirely reconstructed in single contigs without any gap. Using the combination of ab initio prediction, protein homology, and transcriptome evidence, we identified 58,287 and 56,725 protein‐coding genes in Williams 82 and Lee, respectively. The genome assemblies and annotations will serve as a valuable resource for studying soybean genomics and genetics and accelerating soybean improvement.

The nutritional value of soybean [Glycine max (L.) Merr.] for animals is influenced by soluble carbohydrates, such as sucrose and stachyose. Although sucrose is nutritionally desirable, stachyose is an antinutrient causing diarrhea and flatulence in non-ruminant animals. We conducted a genome-wide association study of 220 soybean accessions using 21,317 single nucleotide polymorphisms (SNPs) from the SoySNP50K iSelect Beadchip data to identify significant SNPs associated with sucrose and stachyose content. Seven significant SNPs were identified for sucrose content across chromosomes (Chrs.) 2, 8, 12, 17, and 20, while thirteen significant SNPs were identified for stachyose content across Chrs. 2, 5, 8, 9, 10, 13, 14, and 15. Among those significant SNPs, three sucrose-related SNPs on Chrs. 8 and 17 were novel, while twelve stachyose-related SNPs on Chrs. 2, 5, 8, 9, 10, 13, 14, and 15 were novel. Based on Phytozome, STRING, and GO annotation, 17 and 24 candidate genes for sucrose and stachyose content, respectively, were highly associated with the carbohydrate metabolic pathway. Among these, the publicly available RNA-seq Atlas database highlighted four candidate genes associated with sucrose (Glyma.08g361200 and Glyma.17g258100) and stachyose (Glyma.05g025300 and Glyma.13g077900) content, which had higher gene expression levels in developing seed and multiple parts of the soybean plant. The results of this study will extend knowledge of the molecular mechanism and genetic basis underlying sucrose and stachyose content in soybean seed. Furthermore, the novel candidate genes and SNPs can be valuable genetic resources that soybean breeders may utilize to modify carbohydrate profiles for animal and human usage.

The adoption of dicamba-tolerant (DT) soybean in the United States resulted in extensive off-target dicamba damage to non-DT vegetation across soybean-producing states. Although soybeans are highly sensitive to dicamba, the intensity of observed symptoms and yield losses are affected by the genetic background of genotypes. Thus, the objective of this study was to detect novel marker-trait associations and expand on previously identified genomic regions related to soybean response to off-target dicamba. A total of 551 non-DT advanced breeding lines derived from 232 unique bi-parental populations were phenotyped for off-target dicamba across nine environments for three years. Breeding lines were genotyped using the Illumina Infinium BARCSoySNP6K BeadChip. Filtered SNPs were included as predictors in Random Forest (RF) and Support Vector Machine (SVM) models in a forward stepwise selection loop to identify the combination of SNPs yielding the highest classification accuracy. Both RF and SVM models yielded high classification accuracies (0.76 and 0.79, respectively) with minor extreme misclassifications (observed tolerant predicted as susceptible, and vice-versa). Eight genomic regions associated with off-target dicamba tolerance were identified on chromosomes 6 [Linkage Group (LG) C2], 8 (LG A2), 9 (LG K), 10 (LG O), and 19 (LG L). Although the genetic architecture of tolerance is complex, high classification accuracies were obtained when including the major effect SNP identified on chromosome 6 as the sole predictor. In addition, candidate genes with annotated functions associated with phases II (conjugation of hydroxylated herbicides to endogenous sugar molecules) and III (transportation of herbicide conjugates into the vacuole) of herbicide detoxification in plants were co-localized with significant markers within each genomic region. Genomic prediction models, as reported in this study, can greatly facilitate the identification of genotypes with superior tolerance to off-target dicamba.

Abstract Soybean (Glycine max) is a major crop that contributes more than half of global oilseed production. Much research has been directed towards improvement of the fatty acid profile of soybean seeds through marker assisted breeding. Recently published soybean pangenomes, based on thousands of soybean lines, provide an opportunity to identify new alleles that may be involved in fatty acid biosynthesis. In this study, we identify fatty acid biosynthesis genes in soybean pangenomes based on sequence identity with known genes and examine their sequence diversity across diverse soybean collections. We find three possible instances of a gene missing in wild soybean, including FAD8 and FAD2‐2D, which may be involved in oleic and linoleic acid desaturation, respectively, although we recommend follow‐up research to verify the absence of these genes. More than half of the 53 fatty acid biosynthesis genes identified contained missense variants, including one linked with a previously identified QTL for oil quality. These variants were present in multiple studies based on either short read mappings or alignment of reference grade genomes. Missense variants were found in previously characterized genes including FAD2‐1A and FAD2‐1B, both of which are involved in desaturation of oleic acid, as well as uncharacterized candidate fatty acid biosynthesis genes. We find that the frequency of missense alleles in fatty acid biosynthesis genes has been reduced significantly more than the global average frequency of missense mutations during domestication, and missense variation in some genes is near absent in modern cultivars. This could be due to the selection for fatty acid profiles in seed, though future work should be conducted towards understanding the phenotypic impacts of these variants.

Understanding the genetic basis of seed Ni and Mo is essential. Since soybean is a major crop in the world and a major source for nutrients, including Ni and Mo, the objective of the current research was to map genetic regions (quantitative trait loci, QTL) linked to Ni and Mo concentrations in soybean seed. A recombinant inbred line (RIL) population was derived from a cross between ‘Forrest’ and ‘Williams 82’ (F × W82). A total of 306 lines was used for genotyping using 5405 single nucleotides polymorphism (SNP) markers using Infinium SNP6K BeadChips. A two-year experiment was conducted and included the parents and the RIL population. One experiment was conducted in 2018 in North Carolina (NC), and the second experiment was conducted in Illinois in 2020 (IL). Logarithm of the odds (LOD) of ≥2.5 was set as a threshold to report identified QTL using the composite interval mapping (CIM) method. A wide range of Ni and Mo concentrations among RILs was observed. A total of four QTL (qNi-01, qNi-02, and qNi-03 on Chr 2, 8, and 9, respectively, in 2018, and qNi-01 on Chr 20 in 2020) was identified for seed Ni. All these QTL were significantly (LOD threshold > 2.5) associated with seed Ni, with LOD scores ranging between 2.71–3.44, and with phenotypic variance ranging from 4.48–6.97%. A total of three QTL for Mo (qMo-01, qMo-02, and qMo-03 on Chr 1, 3, 17, respectively) was identified in 2018, and four QTL (qMo-01, qMo-02, qMo-03, and qMo-04, on Chr 5, 11, 14, and 16, respectively) were identified in 2020. Some of the current QTL had high LOD and significantly contributed to the phenotypic variance for the trait. For example, in 2018, Mo QTL qMo-01 on Chr 1 had LOD of 7.8, explaining a phenotypic variance of 41.17%, and qMo-03 on Chr 17 had LOD of 5.33, with phenotypic variance explained of 41.49%. In addition, one Mo QTL (qMo-03 on Chr 14) had LOD of 9.77, explaining 51.57% of phenotypic variance related to the trait, and another Mo QTL (qMo-04 on Chr 16) had LOD of 7.62 and explained 49.95% of phenotypic variance. None of the QTL identified here were identified twice across locations/years. Based on a search of the available literature and of SoyBase, the four QTL for Ni, identified on Chr 2, 8, 9, and 20, and the five QTL associated with Mo, identified on Chr 1, 17, 11, 14, and 16, are novel and not previously reported. This research contributes new insights into the genetic mapping of Ni and Mo, and provides valuable QTL and molecular markers that can potentially assist in selecting Ni and Mo levels in soybean seeds.

Late embryogenesis abundant (LEA) proteins play important roles in regulating plant growth and responses to various abiotic stresses. In this research, a genome-wide survey was conducted to recognize the LEA genes in Glycine max. A total of 74 GmLEA was identified and classified into nine subfamilies based on their conserved domains and the phylogenetic analysis. Subcellular localization, the duplication of genes, gene structure, the conserved motif, and the prediction of cis-regulatory elements and tissue expression pattern were then conducted to characterize GmLEAs. The expression profile analysis indicated that the expression of several GmLEAs was a response to drought and salt stress. The co-expression-based gene network analysis suggested that soybean LEA proteins may exert regulatory effects through the metabolic pathways. We further explored GnLEA4_19 function in Arabidopsis and the results suggests that overexpressed GmLEA4_19 in Arabidopsis increased plant height under mild or serious drought stress. Moreover, the overexpressed GmLEA4_19 soybean also showed a drought tolerance phenotype. These results indicated that GmLEA4_19 plays an important role in the tolerance to drought and will contribute to the development of the soybean transgenic with enhanced drought tolerance and better yield. Taken together, this study provided insight for better understanding the biological roles of LEA genes in soybean.

Soybean seed sugars are among the most abundant beneficial compounds for human and animal consumption in soybean seeds. Higher seed sugars such as sucrose are desirable as they contribute to taste and flavor in soy-based food. Therefore, the objectives of this study were to use the ‘Forrest’ by ‘Williams 82’ (F × W82) recombinant inbred line (RIL) soybean population (n = 309) to identify quantitative trait loci (QTLs) and candidate genes that control seed sugar (sucrose, stachyose, and raffinose) contents in two environments (North Carolina and Illinois) over two years (2018 and 2020). A total of 26 QTLs that control seed sugar contents were identified and mapped on 16 soybean chromosomes (chrs.). Interestingly, five QTL regions were identified in both locations, Illinois and North Carolina, in this study on chrs. 2, 5, 13, 17, and 20. Amongst 57 candidate genes identified in this study, 16 were located within 10 Megabase (MB) of the identified QTLs. Amongst them, a cluster of four genes involved in the sugars’ pathway was collocated within 6 MB of two QTLs that were detected in this study on chr. 17. Further functional validation of the identified genes could be beneficial in breeding programs to produce soybean lines with high beneficial sucrose and low raffinose family oligosaccharides.

Nutrient-efficient root system architecture (RSA) is becoming an important breeding objective for generating crop varieties with improved nutrient and water acquisition efficiency. Genetic variants shaping soybean RSA is key in improving nutrient and water acquisition. Here, we report on the use of an improved 2-dimensional high-throughput root phenotyping platform that minimizes background noise by imaging pouch-grown root systems submerged in water. We also developed a background image cleaning Python pipeline that computationally removes images of small pieces of debris and filter paper fibers, which can be erroneously quantified as root tips. This platform was used to phenotype root traits in 286 soybean lines genotyped with 5.4 million single-nucleotide polymorphisms. There was a substantially higher correlation in manually counted number of root tips with computationally quantified root tips (95% correlation), when the background was cleaned of nonroot materials compared to root images without the background corrected (79%). Improvements in our RSA phenotyping pipeline significantly reduced overestimation of the root traits influenced by the number of root tips. Genome-wide association studies conducted on the root phenotypic data and quantitative gene expression analysis of candidate genes resulted in the identification of 3 putative positive regulators of root system depth, total root length and surface area, and root system volume and surface area of thicker roots (DOF1-like zinc finger transcription factor, protein of unknown function, and C2H2 zinc finger protein). We also identified a putative negative regulator (gibberellin 20 oxidase 3) of the total number of lateral roots.

Abstract Improving seed yield is one of the main targets of soybean [Glycine max (L.) Merr.] breeding. Identification of loci that influence productivity and understanding their genetic mechanism will help marker‐assisted trait introgression. The present study evaluated a diverse panel of 541 soybean genotypes consisting of three maturity groups (MGs III–V) in four environments in Kansas, U.S. Data on seed yield, seed weight, shattering resistance, days to maturity, and plant height showed significant genotype, environmental, and genotype × environment interaction variations. Seed yield and shattering had moderate broad‐sense heritability (90%). The SoySNP50K iSelect BeadChip dataset was used to identify significantly associated loci via genome‐wide association studies (GWAS). A total of 19 single‐nucleotide polymorphisms (SNPs) were significantly associated with seed yield. Particularly, two stable seed yield quantitative trait loci (QTL) on chromosomes 9 and 17 were consistently detected in at least three out of four environments. Candidate gene analysis surrounding seed yield QTL on chromosome 9 showed that Glyma.09G048900, an oxygen binding protein, was the closest to the QTL peak. Similarly, Glyma.17G090200 and Glyma.17G090400 were within 20‐kb region of the seed yield QTL on chromosome 17. The candidate genes warrant further analysis to determine their functional mechanisms and develop markers for seed yield improvement.

The widespread adoption of genetically modified (GM) dicamba-tolerant (DT) soybean was followed by numerous reports of off-target dicamba damage and yield losses across most soybean-producing states. In this study, a subset of the USDA Soybean Germplasm Collection consisting of 382 genetically diverse soybean accessions originating from 15 countries was used to identify genomic regions associated with soybean response to off-target dicamba exposure. Accessions were genotyped with the SoySNP50K BeadChip and visually screened for damage in environments with prolonged exposure to off-target dicamba. Two models were implemented to detect significant marker-trait associations: the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) and a model that allows the inclusion of population structure in interaction with the environment (G×E) to account for variable patterns of genotype responses in different environments. Most accessions (84%) showed a moderate response, either moderately tolerant or moderately susceptible, with approximately 8% showing tolerance and susceptibility. No differences in off-target dicamba damage were observed across maturity groups and centers of origin. Both models identified significant associations in regions of chromosomes 10 and 19. The BLINK model identified additional significant marker-trait associations on chromosomes 11, 14, and 18, while the G×E model identified another significant marker-trait association on chromosome 15. The significant SNPs identified by both models are located within candidate genes possessing annotated functions involving different phases of herbicide detoxification in plants. These results entertain the possibility of developing non-GM soybean cultivars with improved tolerance to off-target dicamba exposure and potentially other synthetic auxin herbicides. Identification of genetic sources of tolerance and genomic regions conferring higher tolerance to off-target dicamba may sustain and improve the production of other non-DT herbicide soybean production systems, including the growing niche markets of organic and conventional soybean.

The J-protein family comprises molecular chaperones involved in plant growth, development, and stress responses. Little is known about this gene family in soybean. Hence, we characterized J-protein genes in soybean, with the most highly expressed and responsive during flower and seed development. We also revealed their phylogeny, structure, motif analysis, chromosome location, and expression. Based on their evolutionary links, we divided the 111 potential soybean J-proteins into 12 main clades (I–XII). Gene-structure estimation revealed that each clade had an exon-intron structure resembling or comparable to others. Most soybean J-protein genes lacked introns in Clades I, III, and XII. Moreover, transcriptome data obtained from a publicly accessible soybean database and RT-qPCR were used to examine the differential expression of DnaJ genes in various soybean tissues and organs. The expression level of DnaJ genes indicated that, among 14 tissues, at least one tissue expressed the 91 soybean genes. The findings suggest that J-protein genes could be involved in the soybean growth period and offer a baseline for further functional research into J-proteins' role in soybean. One important application is the identification of J-proteins that are highly expressed and responsive during flower and seed development in soybean. These genes likely play crucial roles in these processes, and their identification can contribute to breeding programs to improve soybean yield and quality.