GmRmd1 encodes a TIR-NBS-BSP protein and confers resistance to powdery mildew in soybean

Identification of Quantitative Trait Loci for Node Number, Pod Number, and Seed Number in Soybean

Optimizing soybean yield remains a crucial challenge in meeting global food security demands. In this study, we report a comprehensive genetic analysis of yield-related traits in soybeans using a recombinant inbred line (RIL) population derived from crosses between 'Qihuang 34' (GH34) and 'Dongsheng 16' (DS16). Phenotypic analysis across two years (2023-2024) revealed significant variations between parental lines. Through high-density genetic mapping with 6297 SLAF markers spanning 2945.26 cM across 20 chromosomes, we constructed a genetic map with an average marker distance of 0.47 cM and 99.17% of gaps under 5 cM. QTL analysis identified ten significant loci across both years: in 2023, we detected six QTLs, including a major main stem node number (MSNN) QTL on chromosome 19 (LOD = 22.59, PVE = 24.57%), two seed number (SN) QTLs on chromosomes 14 and 18 (LOD = 2.52-2.85, PVE = 7.35% combined), and one pod number (PN) QTL on chromosome 20 (LOD = 4.68, PVE = 5.85%). The 2024 analysis revealed four major QTLs, notably a cluster on chromosome 19 harboring significant loci for MSNN (LOD = 37.92, PVE = 43.59%), PN (LOD = 18.16, PVE = 23.02%), and SN (LOD = 15.24, PVE = 19.59%). Within the stable chromosome 19 region, we identified seventeen candidate genes involved in crucial developmental processes. Gene expression analysis revealed distinct temporal patterns between parental lines during vegetative and reproductive stages, with GH34 showing dramatically higher expression of key reproductive genes Glyma.19G201300 and Glyma.19G201400 during the R1 stage. Our findings provide new insights into the genetic architecture of soybean stem node development and yield components, offering multiple promising targets for molecular breeding programs aimed at crop improvement.

Identification of quantitative trait loci for lodging and related agronomic traits in soybean (Glycine max [L.] Merr.)

BACKGROUND

Lodging, a crucial agronomic trait linked to soybean yield, poses a significant challenge in soybean production. Nevertheless, there has been less research on soybean lodging compared to other important agronomic traits, hindering progress in breeding high-yield soybeans. Our goals were to investigate lodging, pinpoint quantitative trait loci (QTL) linked to lodging, and forecast potential candidate genes linked to this trait. To achieve this, we employed a recombinant inbred line (RIL) population derived from a cross between Guizao 1 and B13 (GB) across various environments.

RESULTS

The lodging score of the RIL population was found to be significantly positively correlated with flowering time, maturity time, plant height, number of main stem nodes, stem diameter, and internode length, with correlation coefficients ranging from 0.457 to 0.783. A total of 84 QTLs associated with soybean lodging and related traits were identified using the GB population. The contribution of phenotypic variance ranged from 1.26 to 66.87%, with LOD scores ranging from 2.52 to 69.22. Additionally, within these QTLs, a stable major QTL associated with lodging was newly discovered in the GB population. Out of the ten major QTLs associated with other related traits, nine of them were situated within the qLD-4-1 interval of the major lodging score locus, displaying phenotypic variations ranging from 12.10 to 66.87%. Specific alterations in gene expression were revealed through the analysis of resequencing data from the two parental lines, potentially indicating their significant roles in lodging. Subsequently, it was determined through qRT-PCR that four genes are likely to be the major genes controlling soybean lodging.

CONCLUSIONS

This study's findings offer valuable insights into the genetic underpinnings of soybean lodging resistance traits. By comprehending the potential genetic factors associated with lodging, this research lays the groundwork for breeding high-yield soybeans with improved lodging resistance.

Allelic variability in the Rpp1 locus conferring resistance to Asian soybean rust revealed by genome-wide association

Soybean is a crucial crop for the Brazilian economy, but it faces challenges from the biotrophic fungus Phakopsora pachyrhizi, which causes Asian Soybean Rust (ASR). In this study, we aimed to identify SNPs associated with resistance within the Rpp1 locus, which is effective against Brazilian ASR populations. We employed GWAS and re-sequencing analyzes to pinpoint SNP markers capable of differentiating between soybean accessions harboring the Rpp1, Rpp1-b and other alternative alleles in the Rpp1 locus and from susceptible soybean cultivars. Seven SNP markers were found to be associated with ASR resistance through GWAS, with three of them defining haplotypes that efficiently distinguished the accessions based on their ASR resistance and source of the Rpp gene. These haplotypes were subsequently validated using a bi-parental population and a diverse set of Rpp sources, demonstrating that the GWAS markers co-segregate with ASR resistance. We then examined the presence of these haplotypes in a diverse set of soybean genomes worldwide, finding a few new potential sources of Rpp1/Rpp1-b. Further genomic sequence analysis revealed nucleotide differences within the genes present in the Rpp1 locus, including the ULP1-NBS-LRR genes, which are potential R gene candidates. These results provide valuable insights into ASR resistance in soybean, thus helping the development of resistant soybean varieties through genetic breeding programs.

A comprehensive review of soybean RNL and TIR domain proteins

Both prokaryotic and eukaryotic organisms use the nucleotide-binding domain/leucine-rich repeat (NBD/LRR)-triggered immunity (NLR-triggered immunity) signaling pathway to defend against pathogens. Plant NLRs are intracellular immune receptors that can bind to effector proteins secreted by pathogens. Dicotyledonous plants express a type of NLR known as TIR domain-containing NLRs (TNLs). TIR domains are enzymes that catalyze the production of small molecules that are essential for immune signaling and lead to plant cell death. The activation of downstream TNL signaling components, such as enhanced disease susceptibility 1 (EDS1), phytoalexin deficient 4 (PAD4), and senescence-associated gene 101 (SAG101), is facilitated by these small molecules. Helper NLRs (hNLRs) and the EDS1-PAD4/SAG101 complex associate after activation, causing the hNLRs to oligomerize, translocate to the plasma membrane (PM), and produce cation-selective channels. According to a recent theory, cations enter cells through pores created by oligomeric hNLRs and trigger cell death. Occasionally, TNLs can self-associate to create higher-order oligomers. Here, we categorized soybean TNLs based on the protein domains that they possess. We believe that TNLs may help soybean plants effectively fight pathogens by acting as a source of genetic resistance. In summary, the purpose of this review is to elucidate the range of TNLs that are expressed in soybean.

Genome-wide association study and haplotype analysis reveal novel candidate genes for resistance to powdery mildew in soybean

Powdery mildew disease (PMD) is caused by the obligate biotrophic fungus Microsphaera diffusa Cooke & Peck (M. diffusa) and results in significant yield losses in soybean (Glycine max (L.) Merr.) crops. By identifying disease-resistant genes and breeding soybean accessions with enhanced resistance, we can effectively mitigate the detrimental impact of PMD on soybeans. We analyzed PMD resistance in a diversity panel of 315 soybean accessions in two locations over 3 years, and candidate genes associated with PMD resistance were identified through genome-wide association studies (GWAS), haplotype analysis, qRT-PCR, and EMS mutant analysis. Based on the GWAS approach, we identified a region on chromosome 16 (Chr16) in which 21 genes form a gene cluster that is highly correlated with PMD resistance. In order to validate and refine these findings, we conducted haplotype analysis of 21 candidate genes and indicated there are single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) variations of Glyma.16G214000, Glyma.16G214200, Glyma.16G215100 and Glyma.16G215300 within the coding and promoter regions that exhibit a strong association with resistance against PMD. Subsequent structural analysis of candidate genes within this cluster revealed that in 315 accessions, the majority of accessions exhibited resistance to PMD when Glyma.16G214300, Glyma.16G214800 and Glyma.16G215000 were complete; however, they demonstrated susceptibility to PMD when these genes were incomplete. Quantitative real-time PCR assays (qRT-PCR) of possible candidate genes showed that 14 candidate genes (Glyma.16G213700, Glyma.16G213800, Glyma.16G213900, Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214500, Glyma.16G214585, Glyma.16G214669, Glyma.16G214700, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300) were involved in PMD resistance. Finally, we evaluated the PMD resistance of mutant lines from the Williams 82 EMS mutations library, which revealed that mutants of Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300, exhibited sensitivity to PMD. Combined with the analysis results of GWAS, haplotypes, qRT-PCR and mutants, the genes Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300 were identified as highly correlated with PMD resistance. The candidate genes identified above are all NLR family genes, and these discoveries deepen our understanding of the molecular basis of PMD resistance in soybeans and will be useful for guiding breeding strategies.

Application of bulk segregant RNA-Seq (BSR-Seq) and allele-specific primers to study soybean powdery mildew resistance

BACKGROUND

Powdery mildew (PM) is one of the important soybean diseases, and host resistance could practically contribute to soybean PM management. To date, only the Rmd locus on chromosome (Chr) 16 was identified through traditional QTL mapping and GWAS, and it remains unclear if the bulk segregant RNA-Seq (BSR-Seq) methodology is feasible to explore additional PM resistance that might exist in other varieties.

RESULTS

BSR-Seq was applied to contrast genotypes and gene expressions between the resistant bulk (R bulk) and the susceptible bulk (S bulk), as well as the parents. The ∆(SNP-index) and G' value identified several QTL and significant SNPs/Indels on Chr06, Chr15, and Chr16. Differentially expressed genes (DEGs) located within these QTL were identified using HISAT2 and Kallisto, and allele-specific primers (AS-primers) were designed to validate the accuracy of phenotypic prediction. While the AS-primers on Chr06 or Chr15 cannot distinguish the resistant and susceptible phenotypes, AS-primers on Chr16 exhibited 82% accuracy prediction with an additive effect, similar to the SSR marker Satt431.

CONCLUSIONS

Evaluation of additional AS-primers in the linkage disequilibrium (LD) block on Chr16 further confirmed the resistant locus, derived from the resistant parental variety 'Kaohsiung 11' ('KS11'), not only overlaps with the Rmd locus with unique up-regulated LRR genes (Glyma.16G213700 and Glyma.16G215100), but also harbors a down-regulated MLO gene (Glyma.16G145600). Accordingly, this study exemplified the feasibility of BSR-Seq in studying biotrophic disease resistance in soybean, and showed the genetic makeup of soybean variety 'KS11' comprising the Rmd locus and one MLO gene.

Enhancing powdery mildew resistance in soybean by targeted mutation of MLO genes using the CRISPR/Cas9 system

BACKGROUND

Powdery mildew is a major disease that causes great losses in soybean yield and seed quality. Disease-resistant varieties, which are generated by reducing the impact of susceptibility genes through mutation in host plants, would be an effective approach to protect crops from this disease. The Mildew Locus O (MLO) genes are well-known susceptibility genes for powdery mildew in plant. In this study, we utilized the CRISPR/Cas9 system to induce targeted mutations in the soybean GmMLO genes to improve powdery mildew resistance.

RESULTS

A dual-sgRNA CRISPR/Cas9 construct was designed and successfully transferred into the Vietnamese soybean cultivar DT26 through Agrobacterium tumefaciens-mediated transformation. Various mutant forms of the GmMLO genes including biallelic, chimeric and homozygous were found at the T0 generation. The inheritance and segregation of CRISPR/Cas9-induced mutations were confirmed and validated at the T1 and T2 generations. Out of six GmMLO genes in the soybean genome, we obtained the Gmmlo02/Gmmlo19/Gmmlo23 triple and Gmmlo02/Gmmlo19/Gmmlo20/Gmmlo23 quadruple knockout mutants at the T2 generation. When challenged with Erysiphe diffusa, a fungus that causes soybean powdery mildew, all mutant plants showed enhanced resistance to the pathogen, especially the quadruple mutant. The powdery mildew severity in the mutant soybeans was reduced by up to 36.4% compared to wild-type plants. In addition, no pleiotropic effect on soybean growth and development under net-house conditions was observed in the CRISPR/Cas9 mutants.

CONCLUSIONS

Our results indicate the involvement of GmMLO02, GmMLO19, GmMLO20 and GmMLO23 genes in powdery mildew susceptibility in soybean. Further research should be conducted to investigate the roles of individual tested genes and the involvement of other GmMLO genes in this disease infection mechanism. Importantly, utilizing the CRISPR/Cas9 system successfully created the Gmmlo transgene-free homozygous mutant lines with enhanced resistance to powdery mildew, which could be potential materials for soybean breeding programs.

Identification of major QTLs for soybean seed size and seed weight traits using a RIL population in different environments

INTRODUCTION

The seed weight of soybean [Glycine max (L.) Merr.] is one of the major traits that determine soybean yield and is closely related to seed size. However, the genetic basis of the synergistic regulation of traits related to soybean yield is unclear.

METHODS

To understand the molecular genetic basis for the formation of soybean yield traits, the present study focused on QTLs mapping for seed size and weight traits in different environments and target genes mining.

RESULTS

A total of 85 QTLs associated with seed size and weight traits were identified using a recombinant inbred line (RIL) population developed from Guizao1×B13 (GB13). We also detected 18 environmentally stable QTLs. Of these, qSL-3-1 was a novel QTL with a stable main effect associated with seed length. It was detected in all environments, three of which explained more than 10% of phenotypic variance (PV), with a maximum of 15.91%. In addition, qSW-20-3 was a novel QTL with a stable main effect associated with seed width, which was identified in four environments. And the amount of phenotypic variance explained (PVE) varied from 9.22 to 21.93%. Five QTL clusters associated with both seed size and seed weight were summarized by QTL cluster identification. Fifteen candidate genes that may be involved in regulating soybean seed size and weight were also screened based on gene function annotation and GO enrichment analysis.

DISCUSSION

The results provide a biologically basic reference for understanding the formation of soybean seed size and weight traits.