Transcriptome-wide expression landscape and starch synthesis pathway co-expression network in sorghum.

The gene expression landscape across different tissues and developmental stages reflects their biological functions and evolutionary patterns. Integrative and comprehensive analyses of all transcriptomic data in an organism are instrumental to obtaining a comprehensive picture of gene expression landscape. Such studies are still very limited in sorghum, which limits the discovery of the genetic basis underlying complex agricultural traits in sorghum. We characterized the genome-wide expression landscape for sorghum using 873 RNA-sequencing (RNA-seq) datasets representing 19 tissues. Our integrative analysis of these RNA-seq data provides the most comprehensive transcriptomic atlas for sorghum, which will be valuable for the sorghum research community for functional characterizations of sorghum genes. Based on the transcriptome atlas, we identified 595 housekeeping genes (HKGs) and 2080 tissue-specific expression genes (TEGs) for the 19 tissues. We identified different gene features between HKGs and TEGs, and we found that HKGs have experienced stronger selective constraints than TEGs. Furthermore, we built a transcriptome-wide co-expression network (TW-CEN) comprising 35 modules with each module enriched in specific Gene Ontology terms. High-connectivity genes in TW-CEN tend to express at high levels while undergoing intensive selective pressure. We also built global and seed-preferential co-expression networks of starch synthesis pathways, which indicated that photosynthesis and microtubule-based movement play important roles in starch synthesis. The global transcriptome atlas of sorghum generated by this study provides an important functional genomics resource for trait discovery and insight into starch synthesis regulation in sorghum.

Genetic basis and selection of glyceollin elicitation in wild soybean.

Glyceollins, a family of phytoalexins elicited in legume species, play crucial roles in environmental stress response (e.g., defending against pathogens) and human health. However, little is known about the genetic basis of glyceollin elicitation. In the present study, we employed a metabolite-based genome-wide association (mGWA) approach to identify candidate genes involved in glyceollin elicitation in genetically diverse and understudied wild soybeans subjected to soybean cyst nematode. In total, eight SNPs on chromosomes 3, 9, 13, 15, and 20 showed significant associations with glyceollin elicitation. Six genes fell into two gene clusters that encode glycosyltransferases in the phenylpropanoid pathway and were physically close to one of the significant SNPs (ss715603454) on chromosome 9. Additionally, transcription factors (TFs) genes such as MYB and WRKY were also found as promising candidate genes within close linkage to significant SNPs on chromosome 9. Notably, four significant SNPs on chromosome 9 show epistasis and a strong signal for selection. The findings describe the genetic foundation of glyceollin biosynthesis in wild soybeans; the identified genes are predicted to play a significant role in glyceollin elicitation regulation in wild soybeans. Additionally, how the epistatic interactions and selection influence glyceollin variation in natural populations deserves further investigation to elucidate the molecular mechanism of glyceollin biosynthesis.

Genome-Wide Association Studies Prioritize Genes Controlling Seed Size and Reproductive Period Length in Soybean.

Hundred-seed weight (HSW) and reproductive period length (RPL) are two major agronomic traits critical for soybean production and adaptation. However, both traits are quantitatively controlled by multiple genes that have yet to be comprehensively elucidated due to the lack of major genes; thereby, the genetic basis is largely unknown. In the present study, we conducted comprehensive genome-wide association analyses (GWAS) of HSW and RPL with multiple sets of accessions that were phenotyped across different environments. The large-scale analysis led to the identification of sixty-one and seventy-four significant QTLs for HSW and RPL, respectively. An ortholog-based search analysis prioritized the most promising candidate genes for the QTLs, including nine genes (TTG2, BZR1, BRI1, ANT, KLU, EOD1/BB, GPA1, ABA2, and ABI5) for HSW QTLs and nine genes (such as AGL8, AGL9, TOC1, and COL4) and six known soybean flowering time genes (E2, E3, E4, Tof11, Tof12, and FT2b) for RPL QTLs. We also demonstrated that some QTLs were targeted during domestication to drive the artificial selection of both traits towards human-favored traits. Local adaptation likely contributes to the increased genomic diversity of the QTLs underlying RPL. The results provide additional insight into the genetic basis of HSW and RPL and prioritize a valuable resource of candidate genes that merits further investigation to reveal the complex molecular mechanism and facilitate soybean improvement.

The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean.

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

Genome-wide identification and characterization of greenbug-inducible NAC transcription factors in sorghum.

BACKGROUND: Sorghum (Sorghum bicolor) is an important cereal crop grown worldwide because of its multipurpose uses such as food, forage, and bioenergy feedstock and its wide range of adaption even in marginal environments. Greenbug can cause severe damage to sorghum plants and yield loss. Plant NAC transcription factors (TFs) have been reported to have diverse functions in plant development and plant defense but has not been studied in sorghum yet. METHODS AND RESULTS: In this study, a comprehensive analysis of the sorghum NAC (SbNAC) gene family was conducted through genome-wide analysis. A total of 112 NAC genes has been identified in the sorghum genome. These SbNAC genes are phylogenetically clustered into 15 distinct subfamilies and unevenly distribute in clusters at the telomeric ends of each chromosome. Twelve pairs of SbNAC genes are possibly involved in the segmental duplication among nine chromosomes except chromosome 10. Structure analysis showed the diverse structures with a highly variable number of exons in the SbNAC genes. Furthermore, most of the SbNAC genes showed specific temporal and spatial expression patterns according to the results of RNA-seq analysis, suggesting their diverse functions during sorghum growth and development. We have also identified nine greenbug-inducible SbNAC genes by comparing the expression profiles between two sorghum genotypes (susceptible BTx623 and resistant PI607900) in response to greenbug infestation. CONCLUSIONS: Our systematic analysis of the NAC gene expression profiles provides both a preliminary survey into their roles in plant defense against insect pests and a useful reference for in-depth characterization of the SbNAC genes and the regulatory network that contributes genetic resistance to aphids.

Soybean type-B response regulator GmRR1 mediates phosphorus uptake and yield by modifying root architecture.

Phosphorus (P) plays a pivotal role in plant growth and development. Low P stress can greatly hamper plant growth. Here, we identified a QTL (named QPH-9-1), which is associated with P efficiency across multiple environments through linkage analysis and genome-wide association study. Furthermore, we successfully cloned the underlying soybean (Glycine max) gene GmRR1 (a soybean type-B Response Regulator 1) that encodes a type-B response regulator protein. Knockout of GmRR1 resulted in a substantial increase in plant height, biomass, P uptake efficiency, and yield-related traits due to the modification of root structure. In contrast, overexpression of GmRR1 in plants resulted in a decrease in these phenotypes. Further analysis revealed that knockout of GmRR1 substantially increased the levels of auxin and ethylene in roots, thereby promoting root hair formation and growth by promoting the formation of root hair primordium and lengthening the root apical meristem. Yeast two-hybrid, bimolecular fluorescence complementation, and dual-luciferase assays demonstrated an interaction between GmRR1 and Histidine-containing Phosphotransmitter protein 1. Expression analysis suggested that these proteins coparticipated in response to low P stress. Analysis of genomic sequences showed that GmRR1 underwent a selection during soybean domestication. Taken together, this study provides further insights into how plants respond to low P stress by modifying root architecture through phytohormone pathways.

A novel soybean malectin-like receptor kinase-encoding gene, GmMLRK1, provides resistance to soybean mosaic virus.

Soybean mosaic virus (SMV) severely damages soybean [Glycine max (L.) Merr.] yield and seed quality. Moreover, the underlying genetic determinants of resistance to SMV remain largely unknown. Here, we performed a genome-wide association study (GWAS) of SMV resistance in a panel of 219 diverse soybean accessions across four environments and identified a new resistance-related gene, GmMLRK1, at the major resistance locus Rsv4 on chromosome 2. GmMLRK1 encodes a malectin-like receptor kinase (RK) that was induced earlier and to a greater degree in leaves of the SMV-resistant cultivar Kefeng No. 1 than in those of the susceptible cultivar Nannong 1138-2 after inoculation. We demonstrated that soybean plants overexpressing GmMLRK1 show broad-spectrum resistance to both strains SC7 and SC3 on the basis of reduced viral accumulation, increased reactive oxygen species production, and local cell death associated with the hypersensitive response. In contrast, GmMLRK1 knockout mutants were more susceptible to both pathotypes. Haplotype analysis revealed the presence of five haplotypes (H1-H5) within the soybean population, and only H1 provided SMV resistance, which was independent of its tightly linked SMV resistance gene RNase-H at the same locus. These results report a novel gene that adds new understanding of SMV resistance and can be used for breeding resistant soybean accessions.

GmFtsH25 overexpression increases soybean seed yield by enhancing photosynthesis and photosynthates.

Increasing plant photosynthetic capacity is a promising approach to boost yields, but it is particularly challenging in C3 crops, such as soybean (Glycine max (L.) Merr.). Here, we identified GmFtsH25, encoding a member of the filamentation temperature-sensitive protein H protease family, as a major gene involved in soybean photosynthesis, using linkage mapping and a genome-wide association study. Overexpressing GmFtsH25 resulted in more grana thylakoid stacks in chloroplasts and increased photosynthetic efficiency and starch content, while knocking out GmFtsH25 produced the opposite phenotypes. GmFtsH25 interacted with photosystem I light harvesting complex 2 (GmLHCa2), and this interaction may contribute to the observed enhanced photosynthesis. GmFtsH25 overexpression lines had superior yield traits, such as yield per plant, compared to the wild type and knockout lines. Additionally, we identified an elite haplotype of GmFtsH25, generated by natural mutations, which appears to have been selected during soybean domestication. Our study sheds light on the molecular mechanism by which GmFtsH25 modulates photosynthesis and provides a promising strategy for improving the yields of soybean and other crops.

GmEIL4 enhances soybean (Glycine max) phosphorus efficiency by improving root system development.

Phosphorus (P) deficiency seriously affects plant growth and development and ultimately limits the quality and yield of crops. Here, a new P efficiency-related major quantitative trait locus gene, GmEIL4 (encoding an ethylene-insensitive 3-like 1 protein), was cloned at qP2, which was identified by linkage analysis and genome-wide association study across four environments. Overexpressing GmEIL4 significantly improved the P uptake efficiency by increasing the number, length and surface area of lateral roots of hairy roots in transgenic soybeans, while interfering with GmEIL4 resulted in poor root phenotypic characteristics compared with the control plants under low P conditions. Interestingly, we found that GmEIL4 interacted with EIN3-binding F box protein 1 (GmEBF1), which may regulate the root response to low P stress. We conclude that the expression of GmEIL4 was induced by low-P stress and that overexpressing GmEIL4 improved P accumulation by regulating root elongation and architecture. Analysis of allele variation of GmEIL4 in 894 soybean accessions suggested that GmEIL4 is undergoing artificial selection during soybean evolution, which will benefit soybean production. Together, this study further elucidates how plants respond to low P stress by modifying root structure and provides insight into the great potential of GmEIL4 in crop P-efficient breeding.

Whole-genome resequencing reveals signature of local adaptation and divergence in wild soybean.

Global climate change has threatened world crop production and food security. Decoding the adaptive genetic basis of wild relatives provides an invaluable genomic resource for climate-smart crop breedinG. Here, we performed whole-genome sequencing of 185 diverse wild soybean (Glycine soja) accessions collected from three major agro-ecological zones in China to parse the genomic basis of local adaptation in wild soybean. The population genomic diversity pattern exhibited clear agro-ecological zone-based population structure, and multiple environmental factors were observed to contribute to the genetic divergence. Demographic analysis shows that wild soybeans from the three ecological zones diverged about 1 × 105 years ago, and then the effective population sizes have undergone different degrees of expansions. Genome-environment association identified multiple genes involved in the local adaptation, such as flowering time and temperature-related genes. A locus containing two adjacent MADS-box transcription factors on chromosome 19 was identified for multiple environmental factors, and it experienced positive selection that enables the adaptation to high-latitude environment. This study provides insights into the genetic mechanism of ecological adaptation in wild soybean that may facilitate climate-resilient soybean breeding.

POWR1 is a domestication gene pleiotropically regulating seed quality and yield in soybean.

Seed protein, oil content and yield are highly correlated agronomically important traits that essentially account for the economic value of soybean. The underlying molecular mechanisms and selection of these correlated seed traits during soybean domestication are, however, less known. Here, we demonstrate that a CCT gene, POWR1, underlies a large-effect protein/oil QTL. A causative TE insertion truncates its CCT domain and substantially increases seed oil content, weight, and yield while decreasing protein content. POWR1 pleiotropically controls these traits likely through regulating seed nutrient transport and lipid metabolism genes. POWR1 is also a domestication gene. We hypothesize that the TE insertion allele is exclusively fixed in cultivated soybean due to selection for larger seeds during domestication, which significantly contributes to shaping soybean with increased yield/seed weight/oil but reduced protein content. This study provides insights into soybean domestication and is significant in improving seed quality and yield in soybean and other crop species.

The purple acid phosphatase GmPAP17 predominantly enhances phosphorus use efficiency in soybean.

Purple acid phosphatase (PAP) is an important plant acid phosphatase, which can secrete to the rhizosphere to decompose organophosphorus, promote phosphorus use efficiency, plant growth and development. However, little is known about the functions of intracellular PAP in plants, especially for soybean. Our previous study integrating QTL mapping and transcriptome analysis identified an promising low phosphorus (LP)-induced gene GmPAP17. Here, we determined that GmPAP17 was mainly expressed in roots and had a strong response to LP stress. Furthermore, and the relative expression in the root of LP tolerant genotypes NN94-156 was significantly greater than that of LP sensitive genotype Bogao after LP stress treatment. The overexpression of GmPAP17 significantly enhanced both acid phosphatase activity and growth performance of hairy roots under LP stress condition, it was vice versa for RNAi interference of GmPAP17, indicating that GmPAP17 plays an important role in P use efficiency. Moreover, yeast two-hybrid and bimolecular fluorescence complementation analysis showed that GmRAP2.2 was involved in the regulation network of GmPAP17. Taken together, our results suggest that GmPAP17 is a novel plant PAP that functions in the adaptation of soybean to LP stress, possibly through its involvement in P recycling in plants.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes.

BACKGROUND: With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research. RESULTS: To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection. CONCLUSION: The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

A polygalacturonase gene PG031 regulates seed coat permeability with a pleiotropic effect on seed weight in soybean.

KEY MESSAGE: A QTL gene PG031 regulates the seed coat permeability and seed weight. The critical SNP that can explain the variation of permeability in soybean population can be used for seed improvement. Seed coat permeability is a critical trait for soybean and is tightly associated with seed storage longevity, germination, soy-food processing, and other commercially important traits. However, the molecular mechanism of such an important trait in soybean is largely unclear. In the present study, we uncovered a polygalacturonase (PG) gene, PG031, which controls seed coat permeability in soybean. PG031 exhibited tissue expression specificity in flowers while it was strongly induced in the seed coat and radical upon imbibition. Subcellular localization localized PG031 to the cell wall, suggesting its role specific to the cell wall of the seed coat. Natural variation analysis reveals three haplotypes (PG031289H, PG031289Y, and PG031Hap3) and the single nucleotide polymorphism (SNP) variation for H289Y may explain the variation in permeability in cultivated soybean population. Overexpression of impermeable allele PG031289H significantly reduced the seed coat permeability and 100-seed weight in transgenic seeds through decreasing intracellular spaces of the osteosclereid layer and parenchyma of the seed coat to decline water accessing the seed. PG031 was also located within a quantitative trait locus (QTL) explaining ~ 15% of total phenotypic variation in permeability, nominating it the QTL gene controlling permeability. PG031289Y allele associated with high permeability and high seed weight is experiencing ongoing artificial selection. The results provide insight into the genetic mechanism of seed coat permeability and indicate its potential for the improvement of permeability-associated seed traits in soybean.

GmWRKY46, a WRKY transcription factor, negatively regulates phosphorus tolerance primarily through modifying root morphology in soybean.

Phosphorus (P) deficiency affects soybean growth and development, resulting in significant reduction of yields. However, the regulatory mechanism of P deficiency tolerance in soybean is still largely unclear. WRKY transcription factors are a family of regulators involved in a variety of abiotic stresses in plants while rarely reported in P deficiency. Here, we demonstrated that a soybean GmWRKY46 gene, belonging to group III of WRKY TF family, was involved in the regulation of P deficiency tolerance in soybean. The expression of GmWRKY46 in low P sensitive soybean varieties was significantly higher than that in tolerant soybean varieties. It was primarily expressed in roots and strongly induced by P deprivation. GmWRKY46 was localized in the nucleus. Compared with the control expressing the empty vector, overexpression of GmWRKY46 in soybean hairy roots exhibited more sensitive phenotypes to low P stress, while the RNA interfered GmWRKY46 significantly enhanced P deficiency tolerance by increasing the proliferation, elongation and P absorption efficiency of hairy roots. Expression patterns of a number of P-responsive genes (GmPht1;1, GmPht1;4, GmPTF1, GmACP1, GmPAP21 and GmExpansin-A7) were altered in both overexpression and gene silenced plants. The results provided a novel insight into how soybean responds to low P stress and new gene that may be used to improve soybean low P tolerance through gene editing approach.

Novel target sites for soybean yield enhancement by photosynthesis.

Photosynthesis plays an important role in plant growth and development. Increasing photosynthetic rate is a main objective of improving crop productivity. Chlorophyll fluorescence is an effective method for quickly evaluating photosynthesis. In this study, four representative chlorophyll fluorescence parameters, that is, maximum quantum efficiency of photosystem II, quantum efficiency of PSII, photochemical quenching, and non-photochemical quenching, of 219 diverse soybean accessions were measured across three environments. The underlying genetic architecture was analyzed by genome-wide association study. Forty-eight SNPs were detected to associate with the four traits and explained 10.43-20.41% of the phenotypic variation. Nine candidate genes in the stable QTLs were predicted. Great differences in the expression levels of the candidate genes existed between the high photosynthetic efficiency accessions and low photosynthetic efficiency accessions. In all, we uncover 17 QTLs associated with photosynthesis-related traits and nine genes that may participate in the regulation of photosynthesis, which can provide references for revealing the genetic mechanism of photosynthesis. These QTLs and candidate genes will provide new targets for crop yield improvement through increasing photosynthesis.

Functional study and elite haplotype identification of soybean phosphate starvation response transcription factors GmPHR14 and GmPHR32.

Phosphorus (P) is one of the important mineral elements required for plant growth and development. However, because of the low mobility in soil, P deficiency has been an important factor limiting soybean production. Here, we identified 14 PHR (phosphate starvation response) genes in soybean genome and verified that two previously unreported GmPHR members, GmPHR14 and GmPHR32, were involved in low-P stress tolerance in soybean. GmPHR14 and GmPHR32 were present in two diverged branches of the phylogenic tree. Both genes were highly expressed in roots and root nodules and were induced by P deficiency. GmPHR14 and GmPHR32 both were expressed in the nucleus. The 211 amino acids in the N terminus of GmPHR32 were found to be required for the transcriptional activity. Overexpressing GmPHR14 or GmPHR32 in soybean hairy roots significantly increased roots and shoots dry weight under low-P condition, and overexpressing GmPHR14 additionally significantly increased roots P concentration under low-P condition. GmPHR14 and GmPHR32 were polymorphic in soybean population and the elite haplotype2 (Hap2) for both genes was preferentially present in improved cultivars and showed significantly higher shoots dry weight under low-P condition than the other two haplotypes. These results suggested GmPHR14 and GmPHR32 both positively regulated low-P responses in soybean, and would shed light on the molecular mechanism of low-P stress tolerance. Furthermore, the identified elite haplotypes would be useful in P-efficient soybean breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01301-z.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development.

BACKGROUND: Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement. RESULTS: In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules. CONCLUSIONS: Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Novel resistance strategies to soybean cyst nematode (SCN) in wild soybean.

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is the most damaging soybean pest worldwide and management of SCN remains challenging. The current SCN resistant soybean cultivars, mainly developed from the cultivated soybean gene pool, are losing resistance due to SCN race shifts. The domestication process and modern breeding practices of soybean cultivars often involve strong selection for desired agronomic traits, and thus, decreased genetic variation in modern cultivars, which consequently resulted in limited sources of SCN resistance. Wild soybean (Glycine soja) is the wild ancestor of cultivated soybean (Glycine max) and it's gene pool is indisputably more diverse than G. max. Our aim is to identify novel resistant genetic resources from wild soybean for the development of new SCN resistant cultivars. In this study, resistance response to HG type 2.5.7 (race 5) of SCN was investigated in a newly identified SCN resistant ecotype, NRS100. To understand the resistance mechanism in this ecotype, we compared RNA seq-based transcriptomes of NRS100 with two SCN-susceptible accessions of G. soja and G. max, as well as an extensively studied SCN resistant cultivar, Peking, under both control and nematode J2-treated conditions. The proposed mechanisms of resistance in NRS100 includes the suppression of the jasmonic acid (JA) signaling pathway in order to allow for salicylic acid (SA) signaling-activated resistance response and polyamine synthesis to promote structural integrity of root cell walls. Our study identifies a set of novel candidate genes and associated pathways involved in SCN resistance and the finding provides insight into the mechanism of SCN resistance in wild soybean, advancing the understanding of resistance and the use of wild soybean-sourced resistance for soybean improvement.

Selection of GmSWEET39 for oil and protein improvement in soybean.

Soybean [Glycine max (L.) Merr.] was domesticated from wild soybean (G. soja Sieb. and Zucc.) and has been further improved as a dual-use seed crop to provide highly valuable oil and protein for food, feed, and industrial applications. However, the underlying genetic and molecular basis remains less understood. Having combined high-confidence bi-parental linkage mapping with high-resolution association analysis based on 631 whole sequenced genomes, we mapped major soybean protein and oil QTLs on chromosome15 to a sugar transporter gene (GmSWEET39). A two-nucleotide CC deletion truncating C-terminus of GmSWEET39 was strongly associated with high seed oil and low seed protein, suggesting its pleiotropic effect on protein and oil content. GmSWEET39 was predominantly expressed in parenchyma and integument of the seed coat, and likely regulates oil and protein accumulation by affecting sugar delivery from maternal seed coat to the filial embryo. We demonstrated that GmSWEET39 has a dual function for both oil and protein improvement and undergoes two different paths of artificial selection. A CC deletion (CC-) haplotype H1 has been intensively selected during domestication and extensively used in soybean improvement worldwide. H1 is fixed in North American soybean cultivars. The protein-favored (CC+) haplotype H3 still undergoes ongoing selection, reflecting its sustainable role for soybean protein improvement. The comprehensive knowledge on the molecular basis underlying the major QTL and GmSWEET39 haplotypes associated with soybean improvement would be valuable to design new strategies for soybean seed quality improvement using molecular breeding and biotechnological approaches.