RNA-seq data comparisons of wild soybean genotypes in response to soybean cyst nematode (Heterodera glycines)

Physiological responses and transcriptome analysis of soybean under gradual water deficit

Soybean is an important food and oil crop widely cultivated globally. However, water deficit can seriously affect the yield and quality of soybeans. In order to ensure the stability and increase of soybean yield and improve agricultural water use efficiency (WUE), research on improving drought tolerance and the efficiency of water utilization of soybeans under drought stress has become particularly important. This study utilized the drought-tolerant variety Heinong 44 (HN44) and the drought-sensitive variety Suinong 14 (SN14) to analyze physiological responses and transcriptome changes during the gradual water deficit at the early seed-filling stage. The results indicated that under drought conditions, HN44 had smaller stomata, higher stomatal density, and lower stomatal conductance (Gs) and transpiration rate as compared to SN14. Additionally, HN44 had a higher abscisic acid (ABA) content and faster changes in stomatal morphology and Gs to maintain a dynamic balance between net photosynthetic rate (Pn) and Gs. Additionally, drought-tolerant variety HN44 had high instantaneous WUE under water deficit. Further, HN44 retained a high level of superoxide dismutase (SOD) activity and proline content, mitigating malondialdehyde (MDA) accumulation and drought-induced damage. Comprehensive analysis of transcriptome data revealed that HN44 had fewer differentially expressed genes (DEGs) under light drought stress, reacting insensitivity to water deficit. At the initial stage of drought stress, both varieties had a large number of upregulated DEGs to cope with the drought stress. Under severe drought stress, HN44 had fewer downregulated genes enriched in the photosynthesis pathway than SN14, while it had more upregulated genes enriched in the ABA-mediated signaling and glutathione metabolism pathways than SN14. During gradual water deficit, HN44 demonstrated better drought-tolerant physiological characteristics and water use efficiency than SN14 through key DEGs such as GmbZIP4, LOC100810474, and LOC100819313 in the major pathways. Key transcription factors were screened and identified, providing further clarity on the molecular regulatory pathways responsible for the physiological differences in drought tolerance among these varieties. This study deepened the understanding of the drought resistance mechanisms in soybeans, providing valuable references for drought-resistant soybean breeding.

Evaluation Soybean Cultivars for Reaction to Heterodera glycines Populations HG Types 7 and 1.3.4.7 in Northeast China

Soybean cyst nematode Heterodera glycines (SCN) is a major threat to global soybean production. Effective management of this disease is dependent on the development of resistant cultivars. Two SCN HG Types, 7 and 1.3.4.7. were previously identified as prevalent H. glycines populations in Northeast China. In order to evaluate soybean cultivars resistant to local SCN populations, 110 domestic commercial soybeans from different regions of Northeast China were assessed in the greenhouse to determine their potential as novel sources of resistance. The results suggested that cultivars responded differently to the two HG types. Of the 110 soybean cultivars evaluated, 24 accessions were classified as resistant or moderately resistant to HG Type 7, and five cultivars were classified as resistant or moderately resistant to HG Type 1.3.4.7. Among the tested cultivars, Kangxian 12 and Qingdou 13 had resistance response to both HG types 7 and 1.3.4.7. Thus, these broad-based SCN cultivars will be the valuable materials in the SCN resistance breeding program.

Transcriptome analysis reveals key drought-stress-responsive genes in soybean

Drought is the most common environmental stress and has had dramatic impacts on soybean (Glycine max L.) growth and yield worldwide. Therefore, to investigate the response mechanism underlying soybean resistance to drought stress, the drought-sensitive cultivar "Liaodou 15" was exposed to 7 (mild drought stress, LD), 17 (moderate drought stress, MD) and 27 (severe drought stress, SD) days of drought stress at the flowering stage followed by rehydration until harvest. A total of 2214, 3684 and 2985 differentially expressed genes (DEGs) in LD/CK1, MD/CK2, and SD/CK3, respectively, were identified by RNA-seq. Weighted gene co-expression network analysis (WGCNA) revealed the drought-response TFs such as WRKY (Glyma.15G021900, Glyma.15G006800), MYB (Glyma.15G190100, Glyma.15G237900), and bZIP (Glyma.15G114800), which may be regulated soybean drought resistance. Second, Glyma.08G176300 (NCED1), Glyma.03G222600 (SDR), Glyma.02G048400 (F3H), Glyma.14G221200 (CAD), Glyma.14G205200 (C4H), Glyma.19G105100 (CHS), Glyma.07G266200 (VTC) and Glyma.15G251500 (GST), which are involved in ABA and flavonoid biosynthesis and ascorbic acid and glutathione metabolism, were identified, suggesting that these metabolic pathways play key roles in the soybean response to drought. Finally, the soybean yield after rehydration was reduced by 50% under severe drought stress. Collectively, our study deepens the understanding of soybean drought resistance mechanisms and provides a theoretical basis for the soybean drought resistance molecular breeding and effectively adjusts water-saving irrigation for soybean under field production.

Functional Characterization of Ubiquitination Genes in the Interaction of Soybean—Heterodera glycines

Ubiquitination is a kind of post-translational modification of proteins that plays an important role in plant response to biotic and abiotic stress. The response of soybean GmPUB genes to soybean cyst nematode (SCN, Heterodera glycines) infection is largely unknown. In this study, quantitative real-time PCR (qRT-PCR) was performed to detect the relative expression of 49 GmPUB genes in susceptible cultivar William 82 and resistant cultivar Huipizhi after SCN inoculation. The results show that GmPUB genes responded to cyst nematode infection at 1 day post-inoculation (dpi), 5 dpi, 10 dpi and 15 dpi. The expression levels of GmPUB16A, GmPUB20A, GmCHIPA, GmPUB33A, GmPUB23A and GmPUB24A were dramatically changed during SCN infection. Furthermore, functional analysis of these GmPUB genes by overexpression and RNAi showed that GmPUB20A, GmPUB33A and GmPUB24A negatively regulated soybean resistance under SCN stress. The results from our present study provide insights into the complicated molecular mechanism of the interaction between soybean and SCN.

Key Soybean Seedlings Drought-Responsive Genes and Pathways Revealed by Comparative Transcriptome Analyses of Two Cultivars

Seedling drought stress is one of the most important constraints affecting soybean yield and quality. To unravel the molecular mechanisms under soybean drought tolerance, we conducted comprehensive comparative transcriptome analyses of drought-tolerant genotype Jindou 21 (JD) and drought-sensitive genotype Tianlong No.1 (N1) seedlings that had been exposed to drought treatment. A total of 6038 and 4112 differentially expressed genes (DEGs) were identified in drought-tolerant JD and drought-sensitive N1, respectively. Subsequent KEGG pathway analyses showed that numerous DEGs in JD are predominately involved in signal transduction pathways, including plant hormone signaling pathway, calcium signaling pathway, and MAPK signaling pathway. Interestingly, JA and BR plant hormone signal transduction pathways were found specifically participating in drought-tolerant JD. Meanwhile, the differentially expressed CPKs, CIPKs, MAPKs, and MAP3Ks of calcium and MAPK signaling pathway were only identified in JD. The number of DEGs involved in transcription factors (TFs) is larger in JD than that of in N1. Moreover, some differently expressed transcriptional factor genes were only identified in drought-tolerant JD, including FAR1, RAV, LSD1, EIL, and HB-PHD. In addition, this study suggested that JD could respond to drought stress by regulating the cell wall remodeling and stress-related protein genes such as EXPs, CALSs, CBPs, BBXs, and RD22s. JD is more drought tolerant than N1 owing to more DEGs being involved in multiple signal transduction pathways (JA, BR, calcium, MAPK signaling pathway), stress-related TFs, and proteins. The above valuable genes and pathways will deepen the understanding of the molecular mechanisms under drought stress in soybean.

Kinase-dead mutation: A novel strategy for improving soybean resistance to soybean cyst nematode Heterodera glycines

Protein kinases phosphorylate proteins for functional changes and are involved in nearly all cellular processes, thereby regulating almost all aspects of plant growth and development, and responses to biotic and abiotic stresses. We generated two independent co-expression networks of soybean genes using control and stress response gene expression data and identified 392 differentially highly interconnected kinase hub genes among the two networks. Of these 392 kinases, 90 genes were identified as "syncytium highly connected hubs", potentially essential for activating kinase signalling pathways in the nematode feeding site. Overexpression of wild-type coding sequences of five syncytium highly connected kinase hub genes using transgenic soybean hairy roots enhanced plant susceptibility to soybean cyst nematode (SCN; Heterodera glycines) Hg Type 0 (race 3). In contrast, overexpression of kinase-dead variants of these five syncytium kinase hub genes significantly enhanced soybean resistance to SCN. Additionally, three of the five tested kinase hub genes enhanced soybean resistance to SCN Hg Type 1.2.5.7 (race 2), highlighting the potential of the kinase-dead approach to generate effective and durable resistance against a wide range of SCN Hg types. Subcellular localization analysis revealed that kinase-dead mutations do not alter protein cellular localization, confirming the structure-function of the kinase-inactive variants in producing loss-of-function phenotypes causing significant decrease in nematode susceptibility. Because many protein kinases are highly conserved and are involved in plant responses to various biotic and abiotic stresses, our approach of identifying kinase hub genes and their inactivation using kinase-dead mutation could be translated for biotic and abiotic stress tolerance.

A Broad Review of Soybean Research on the Ongoing Race to Overcome Soybean Cyst Nematode

Plant pathogens greatly impact food security of the ever-growing human population. Breeding resistant crops is one of the most sustainable strategies to overcome the negative effects of these biotic stressors. In order to efficiently breed for resistant plants, the specific plant-pathogen interactions should be understood. Soybean is a short-day legume that is a staple in human food and animal feed due to its high nutritional content. Soybean cyst nematode (SCN) is a major soybean stressor infecting soybean worldwide including in China, Brazil, Argentina, USA and Canada. There are many Quantitative Trait Loci (QTLs) conferring resistance to SCN that have been identified; however, only two are widely used: rhg1 and Rhg4. Overuse of cultivars containing these QTLs/genes can lead to SCN resistance breakdown, necessitating the use of additional strategies. In this manuscript, a literature review is conducted on research related to soybean resistance to SCN. The main goal is to provide a current understanding of the mechanisms of SCN resistance and list the areas of research that could be further explored.

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.

Genetic Architecture of Early Vigor Traits in Wild Soybean

A worldwide food shortage has been projected as a result of the current increase in global population and climate change. In order to provide sufficient food to feed more people, we must develop crops that can produce higher yields. Plant early vigor traits, early growth rate (EGR), early plant height (EPH), inter-node length, and node count are important traits that are related to crop yield. Glycine soja, the wild counterpart to cultivated soybean, Glycine max, harbors much higher genetic diversity and can grow in diverse environments. It can also cross easily with cultivated soybean. Thus, it holds a great potential in developing soybean cultivars with beneficial agronomic traits. In this study, we used 225 wild soybean accessions originally from diverse environments across its geographic distribution in East Asia. We quantified the natural variation of several early vigor traits, investigated the relationships among them, and dissected the genetic basis of these traits by applying a Genome-Wide Association Study (GWAS) with genome-wide single nucleotide polymorphism (SNP) data. Our results showed positive correlation between all early vigor traits studied. A total of 12 SNPs significantly associated with EPH were identified with 4 shared with EGR. We also identified two candidate genes, Glyma.07G055800.1 and Glyma.07G055900.1, playing important roles in influencing trait variation in both EGR and EPH in G. soja.

Early transcriptional responses to soybean cyst nematode HG Type 0 show genetic differences among resistant and susceptible soybeans

KEY MESSAGE

Root transcriptome profiling of three soybean cultivars and a wild relative infected with soybean cyst nematode at migratory phase revealed differential resistance pathway responses between resistant and susceptible genotypes. The soybean cyst nematode (SCN), Heterodera glycines, is the most serious pathogen of soybean production throughout the world. Using resistant cultivars is the primary management strategy against SCN infestation. To gain insight into the still obscure mechanisms of genetic resistance to nematodes in different soybean genotypes, RNA-Seq profiling of the roots of Glycine max cv. Peking, Fayette, Williams 82, and a wild relative (Glycine soja PI 468916) was performed during SCN infection at the migratory phase. The analysis showed statistically significant changes of expression beginning at eight hours after inoculation in genes associated with defense mechanisms and pathways, such as the phenylpropanoid biosynthesis pathway, plant innate immunity and hormone signaling. Our results indicate the importance of the early plant response to migratory phase nematodes in pathogenicity determination. The transcriptome changes occurring during early SCN infection included a number of genes and pathways specific to the different resistant genotypes. We observed the most extensive resistant transcriptome reaction in PI 468916, where the resistant response was qualitatively different from that of commonly used G. max varieties.

Analysis of Soybean Somatic Embryogenesis Using Chromosome Segment Substitution Lines and Transcriptome Sequencing

Soybean is an important cash crop that is widely used as a source of vegetable protein and edible oil. The regeneration ability of soybean directly affects the application of biotechnology. In this study, we used the exogenous hormone 2,4-D to treat immature embryos. Different levels of somatic incidence were selected from the chromosome segment substitution lines (CSSLs) constructed by SN14 and ZYD00006. Transcriptome sequencing of extreme materials was performed, and 2666 differentially expressed genes were obtained. At the same time, a difference table was generated by combining the data on CSSL rearrangement. In the extreme materials, a total of 93 differentially expressed genes were predicted and were then analyzed by cluster analysis and Gene Ontology (GO) annotation. After screening and annotating the target genes, three differentially expressed genes with hormone pathways were identified. The expression patterns of the target genes were verified by real-time quantitative PCR (qRT-PCR). Haplotype polymorphism detection and linkage disequilibrium analysis were performed on the candidate gene Glyma.09g248200. This study provided more information on the regulation network of soybean somatic embryogenesis and regeneration processes, and further identified important genes in the soybean regeneration process and provided a theoretical basis for accelerating the application of biotechnology to soybean for improving its breeding efficiency.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

Neglected treasures in the wild - legume wild relatives in food security and human health

The legume family (Fabaceae) is the third-largest flowering family with over 18 000 species worldwide that are rich in proteins, oils, and nutrients. However, the production potential of legume-derived food cannot meet increasing global demand. Wild legumes represent a large group of wild species adaptive to diverse habitats and harbor rich genetic diversity for the improvement of the agronomic, nutritional, and medicinal values of the domesticated legumes. Accumulating evidence suggests that the genetic variation retained in these under-exploited leguminous wild relatives can be used to improve crop yield, nutrient contents, and resistance/tolerance to environmental stresses via the integration of omics, genetics, and genome-editing technologies.

Genome-wide association and genomic prediction identifies soybean cyst nematode resistance in common bean including a syntenic region to soybean Rhg1 locus

A genome-wide association study (GWAS) was applied to detect single nucleotide polymorphisms (SNPs) significantly associated with resistance to Heterodera glycines (HG) also known as the soybean cyst nematode (SCN) in the core collection of common bean, Phaseolus vulgaris. There were 84,416 SNPs identified in 363 common bean accessions. GWAS identified SNPs on chromosome (Chr) 1 that were significantly associated with resistance to HG type 2.5.7. These SNPs were in linkage disequilibrium with a gene cluster orthologous to the three genes at the Rhg1 locus in soybean. A novel signal on Chr 7 was detected and associated with resistance to HG type 1.2.3.5.6.7. Genomic predictions (GPs) for resistance to these two SCN HG types in common bean achieved prediction accuracy of 0.52 and 0.41, respectively. Our study generated a high-quality SNP panel for 363 common bean accessions and demonstrated that both GWAS and GP were effective strategies to understand the genetic architecture of SCN resistance in common bean.

Genome-wide analysis of genes encoding core components of the ubiquitin system in soybean (Glycine max) reveals a potential role for ubiquitination in host immunity against soybean cyst nematode

BACKGROUND

Ubiquitination is a major post-translational protein modification that regulates essentially all cellular and physiological pathways in eukaryotes. The ubiquitination process typically involves three distinct classes of enzymes, ubiquitin-activating enzyme (E1), ubiquitin-conjugating enzyme (E2) and ubiquitin ligase (E3). To date, a comprehensive identification and analysis of core components comprising of the whole soybean (Glycine max) ubiquitin system (UBS) has not been reported.

RESULTS

We performed a systematic, genome-wide analysis of genes that encode core members of the soybean UBS in this study. A total of 1431 genes were identified with high confidence to encode putative soybean UBS components, including 4 genes encoding E1s, 71 genes that encode the E2s, and 1356 genes encoding the E3-related components. Among the E3-encoding genes, 760 encode RING-type E3s, 124 encode U-box domain-containing E3s, and 472 encode F-box proteins. To find out whether the identified soybean UBS genes encode active enzymes, a set of genes were randomly selected and the enzymatic activities of their recombinant proteins were tested. Thioester assays indicated proteins encoded by the soybean E1 gene GmUBA1 and the majority of selected E2 genes are active E1 or E2 enzymes, respectively. Meanwhile, most of the purified RING and U-box domain-containing proteins displayed E3 activity in the in vitro ubiquitination assay. In addition, 1034 of the identified soybean UBS genes were found to express in at least one of 14 soybean tissues examined and the transcript level of 338 soybean USB genes were significantly changed after abiotic or biotic (Fusarium oxysporum and Rhizobium strains) stress treatment. Finally, the expression level of a large number of the identified soybean UBS-related genes was found significantly altered after soybean cyst nematode (SCN) treatment, suggesting the soybean UBS potentially plays an important role in soybean immunity against SCN.

CONCLUSIONS

Our findings indicate the presence of a large and diverse number of core UBS proteins in the soybean genome, which suggests that target-specific modification by ubiquitin is a complex and important part of cellular and physiological regulation in soybean. We also revealed certain members of the soybean UBS may be involved in immunity against soybean cyst nematode (SCN). This study sets up an essential foundation for further functional characterization of the soybean UBS in various physiological processes, such as host immunity against SCN.

The Untapped Genetic Reservoir: The Past, Current, and Future Applications of the Wild Soybean (Glycine soja)

There is a considerable demand for crop improvement, especially considering the increasing growth of world population, continuing climatic fluctuations, and rapidly evolving plant pests and pathogens. Crop wild relatives hold great potential in providing beneficial alleles for crop improvement. Wild soybean, Glycine soja (Siebold & Zucc.), the wild ancestor to the domesticated soybean (Glycine max (L.) Merr.), harbors a high level of genetic variation. Research on G. soja has been largely devoted to understanding the domestication history of the soybean, while little effort has been made to explore its genetic diversity for crop improvement. High genomic diversity and expanded traits make G. soja populations an excellent source for soybean improvement. This review summarizes recent successful research examples of applying wild soybeans in dissecting the genetic basis of various traits, with a focus on abiotic/biotic stress tolerance and resistance. We also discuss the limitations of using G. soja. Perspective future research is proposed, including the application of advanced biotechnology and emerging genomic data to further utilize the wild soybean to counterbalance the rising demand for superior crops. We proposed there is an urgent need for international collaboration on germplasm collection, resource sharing, and conservation. We hope to use the wild soybean as an example to promote the exploration and use of wild resources for crop improvement in order to meet future food requirements.