Fine mapping of the Asian soybean rust resistance gene Rpp2 from soybean PI 230970

Comparative Analysis of Transcriptomes Reveals Pathways and Verifies Candidate Genes for Clubroot Resistance in Brassica oleracea

Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.

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.

An pair of an atypical NLR encoding genes confer Asian soybean rust resistance in soybean

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is a devastating disease that is present in all major soybean-producing regions. The limited availability of resistant germplasm has resulted in a scarcity of commercial soybean cultivars that are resistant to the disease. To date, only the Chinese soybean landrace SX6907 has demonstrated an immune response to ASR. In this study, we present the isolation and characterization of Rpp6907-7 and Rpp6907-4, a gene pair that confer broad-spectrum resistance to ASR. Rpp6907-7 and Rpp6907-4 encode atypic nucleotide-binding leucine-rich repeat (NLR) proteins that are found to be required for NLR-mediated immunity. Genetic analysis shows that only Rpp6907-7 confers resistance, while Rpp6907-4 regulates Rpp6907-7 signaling activity by acting as a repressor in the absence of recognized effectors. Our work highlights the potential value of using Rpp6907 in developing resistant soybean cultivars.

Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

RNA-Seq and Comparative Transcriptomic Analyses of Asian Soybean Rust Resistant and Susceptible Soybean Genotypes Provide Insights into Identifying Disease Resistance Genes

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is one of the most destructive foliar diseases that affect soybeans. Developing resistant cultivars is the most cost-effective, environmentally friendly, and easy strategy for controlling the disease. However, the current understanding of the mechanisms underlying soybean resistance to P. pachyrhizi remains limited, which poses a significant challenge in devising effective control strategies. In this study, comparative transcriptomic profiling using one resistant genotype and one susceptible genotype was performed under infected and control conditions to understand the regulatory network operating between soybean and P. pachyrhizi. RNA-Seq analysis identified a total of 6540 differentially expressed genes (DEGs), which were shared by all four genotypes. The DEGs are involved in defense responses, stress responses, stimulus responses, flavonoid metabolism, and biosynthesis after infection with P. pachyrhizi. A total of 25,377 genes were divided into 33 modules using weighted gene co-expression network analysis (WGCNA). Two modules were significantly associated with pathogen defense. The DEGs were mainly enriched in RNA processing, plant-type hypersensitive response, negative regulation of cell growth, and a programmed cell death process. In conclusion, these results will provide an important resource for mining resistant genes to P. pachyrhizi infection and valuable resources to potentially pyramid quantitative resistance loci for improving soybean germplasm.

Genetic Mapping of Seven Kinds of Locus for Resistance to Asian Soybean Rust

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is one of the most serious soybean (Glycine max) diseases in tropical and subtropical regions. To facilitate the development of resistant varieties using gene pyramiding, DNA markers closely linked to seven resistance genes, namely, Rpp1, Rpp1-b, Rpp2, Rpp3, Rpp4, Rpp5, and Rpp6, were identified. Linkage analysis of resistance-related traits and marker genotypes using 13 segregating populations of ASR resistance, including eight previously published by our group and five newly developed populations, identified the resistance loci with markers at intervals of less than 2.0 cM for all seven resistance genes. Inoculation was conducted of the same population with two P. pachyrhizi isolates of different virulence, and two resistant varieties, 'Kinoshita' and 'Shiranui,' previously thought to only harbor Rpp5, was found to also harbor Rpp3. Markers closely linked to the resistance loci identified in this study will be used for ASR-resistance breeding and the identification of the genes responsible for resistance.

Mapping QTLs Controlling Soybean Rust Disease Resistance in Chiang Mai 5, an Induced Mutant Cultivar

Soybean rust (SBR) caused by the fungus Phakopsora pachyrhizi is an important folia disease of soybean (Glycine max). In this study, we identified QTLs controlling SBR in Chiang Mai 5 (CM5), an SBR-resistant cultivar developed by induced mutation breeding. A recombinant inbred line (RIL) population of 108 lines developed from a cross between Sukhothai 2 (SKT2, a susceptible cultivar) and CM5 was evaluated for SBR resistance under field conditions in Thailand. QTL analysis for the resistance in the RIL population identified a single QTL, qSBR18.1, for resistance. qSBR18.1 was mapped to a 212-kb region on chromosome 18 between simple sequence repeat markers Satt288 and sc21_3420 and accounted for 21.31-35.09% depending on the traits evaluated for resistance. The qSBR18.1 interval overlapped with genomic regions containing resistance to P. pachyrhizi 4 (Rpp4), a locus for SBR resistance. Three tightly linked genes, Glyma.18G226250, Glyma.18G226300, and Glyma.18G226500, each encoding leucine-rich repeat-containing protein, were identified as candidate genes for SBR resistance at the qSRB18.1. The qSBR18.1 would be useful for breeding of SBR resistance.

Genomic regions associated with resistance to soybean rust (Phakopsora pachyrhizi) under field conditions in soybean germplasm accessions from Japan, Indonesia and Vietnam

Eight soybean genomic regions, including six never before reported, were found to be associated with resistance to soybean rust (Phakopsora pachyrhizi) in the southeastern USA. Soybean rust caused by Phakopsora pachyrhizi is one of the most important foliar diseases of soybean [Glycine max (L.) Merr.]. Although seven Rpp resistance gene loci have been reported, extensive pathotype variation in and among fungal populations increases the importance of identifying additional genes and loci associated with rust resistance. One hundred and ninety-one soybean plant introductions from Japan, Indonesia and Vietnam, and 65 plant introductions from other countries were screened for resistance to P. pachyrhizi under field conditions in the southeastern USA between 2008 and 2015. The results indicated that 84, 69, and 49% of the accessions from southern Japan, Vietnam or central Indonesia, respectively, had negative BLUP values, indicating less disease than the panel mean. A genome-wide association analysis using SoySNP50K Infinium BeadChip data identified eight genomic regions on seven chromosomes associated with SBR resistance, including previously unreported regions of Chromosomes 1, 4, 6, 9, 13, and 15, in addition to the locations of the Rpp3 and Rpp6 loci. The six unreported genomic regions might contain novel Rpp loci. The identification of additional sources of rust resistance and associated genomic regions will further efforts to develop soybean cultivars with broad and durable resistance to soybean rust in the southern USA.

Identification and Characterization of Novel Sources of Resistance to Rust Caused by Uromyces pisi in Pisum spp

Pea rust is a major disease worldwide caused by Uromyces pisi in temperate climates. Only moderate levels of partial resistance against U. pisi have been identified so far in pea, urging for enlarging the levels of resistance available for breeding. Herein, we describe the responses to U. pisi of 320 Pisum spp. accessions, including cultivated pea and wild relatives, both under field and controlled conditions. Large variations for U. pisi infection response for most traits were observed between pea accessions under both field and controlled conditions, allowing the detection of genotypes with partial resistance. Simultaneous multi-trait indexes were applied to the datasets allowing the identification of partial resistance, particularly in accessions JI224, BGE004710, JI198, JI199, CGN10205, and CGN10206. Macroscopic observations were complemented with histological observations on the nine most resistant accessions and compared with three intermediates and three susceptible ones. This study confirmed that the reduced infection of resistant accessions was associated with smaller rust colonies due to a reduction in the number of haustoria and hyphal tips per colony. Additionally, a late acting hypersensitive response was identified for the first time in a pea accession (PI273209). These findings demonstrate that screening pea collections continues to be a necessary method in the search for complete resistance against U. pisi. In addition, the large phenotypic diversity contained in the studied collection will be useful for further association analysis and breeding perspectives.

Two antimicrobial genes from Aegilops tauschii Cosson identified by the Bacillus subtilis expression system

Antimicrobial genes play an important role as a primary defense mechanism in all multicellular organisms. We chose Bacillus subtilis as a target pathogen indicator and transferred the Aegilops tauschii Cosson cDNA library into B. subtilis cells. Expression of the candidate antimicrobial gene can inhibit B. subtilis cell growth. Using this strategy, we screened six genes that have an internal effect on the indicator bacteria. Then, the secreted proteins were extracted and tested; two genes, AtR100 and AtR472, were found to have strong external antimicrobial activities with broad-spectrum resistance against Xanthomonas oryzae pv. oryzicola, Clavibacter fangii, and Botrytis cinerea. Additionally, thermal stability tests indicated that the antimicrobial activities of both proteins were thermostable. Furthermore, these two proteins exhibited no significant hemolytic activities. To test the feasibility of application at the industrial level, liquid fermentation and spray drying of these two proteins were conducted. Powder dilutions were shown to have significant inhibitory effects on B. cinerea. Fluorescence microscopy and flow cytometry results showed that the purified protein impaired and targeted the cell membranes. This study revealed that these two antimicrobial peptides could potentially be used for replacing antibiotics, which would provide the chance to reduce the emergence of drug resistance.

Near-isogenic soybean lines carrying Asian soybean rust resistance genes for practical pathogenicity validation

Asian soybean rust caused by the fungal pathogen Phakopsora pachyrhizi is the most devastating disease of soybean. The host cultivar specificity of the pathogen shows considerable differentiation depending on the area and season of its emergence. Although resistance genes for P. pachyrhizi (Rpp) have been reported in several soybean varieties, the genetic background of these varieties is highly differentiated. Furthermore, some of the varieties harbor unknown genetic factors in addition to Rpp that could influence resistance reactions against the pathogen. In order to gain a comprehensive understanding of Rpp-P. pachyrhizi interactions, homogenous plant material harboring Rpp genes is necessary. In this study, we bred Rpp-near isogenic lines (Rpp-NILs), which retained identical plant characters originating from a single genetic background, and accordingly showed low-variant compatible/incompatible reactions against the pathogen. These Rpp-NILs can be used as genetic resources for studying P. pachyrhizi epidemiology and elucidating resistance mechanisms. Compatible/incompatible relationships between the soybean rust resistance gene Rpp and isolates of the pathogen P. pachyrhizi are clearly distinguishable using the Rpp-NILs bred in this study.

Assembly and annotation of a draft genome sequence for Glycine latifolia, a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939-Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked-reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome-scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91-bp centromere-specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92-bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein-coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine-specific orthologous gene families. A total of 304 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR-NBS-LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR-receptor-like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost-effectiveness of the application of Chromium linked-reads in diploid plant genome de novo assembly.

Discovery of a seventh Rpp soybean rust resistance locus in soybean accession PI 605823

KEY MESSAGE

A novel Rpp gene from PI 605823 for resistance to Phakopsora pachyrhizi was mapped on chromosome 19. Soybean rust, caused by the obligate biotrophic fungal pathogen Phakopsora pachyrhizi Syd. & P. Syd, is a disease threat to soybean production in regions of the world with mild winters. Host plant resistance conditioned by resistance to P. pachyrhizi (Rpp) genes has been found in numerous soybean accessions, and at least 10 Rpp genes or alleles have been mapped to six genetic loci. Identifying additional disease-resistance genes will facilitate development of soybean cultivars with durable resistance. PI 605823, a plant introduction from Vietnam, was previously identified as resistant to US populations of P. pachyrhizi in greenhouse and field trials. In this study, bulked segregant analysis using an F2 population derived from 'Williams 82' × PI 605823 identified a genomic region associated with resistance to P. pachyrhizi isolate GA12, which had been collected in the US State of Georgia in 2012. To further map the resistance locus, linkage mapping was carried out using single-nucleotide polymorphism markers and phenotypic data from greenhouse assays with an F2:3 population derived from Williams 82 × PI 605823 and an F4:5 population derived from '5601T' × PI 605823. A novel resistance gene, Rpp7, was mapped to a 154-kb interval (Gm19: 39,462,291-39,616,643 Glyma.Wm82.a2) on chromosome 19 that is different from the genomic locations of any previously reported Rpp genes. This new gene could be incorporated into elite breeding lines to help provide more durable resistance to soybean rust.

Fine mapping of the soybean aphid-resistance genes Rag6 and Rag3c from Glycine soja 85-32

KEY MESSAGE

Rag6 and Rag3c were delimited to a 49-kb interval on chromosome 8 and a 150-kb interval on chromosome 16, respectively. Structural variants in the exons of candidate genes were identified. The soybean aphid, an invasive species, has significantly threatened soybean production in North America since 2000. Host-plant resistance is known as an ideal management strategy for aphids. Two novel aphid-resistance loci, Rag6 and Rag3c, from Glycine soja 85-32, were previously detected in a 10.5-cM interval on chromosome 8 and a 7.5-cM interval on chromosome 16, respectively. Defining the exact genomic position of these two genes is critical for improving the effectiveness of marker-assisted selection for aphid resistance and for identification of the functional genes. To pinpoint the locations of Rag6 and Rag3c, four populations segregating for Rag6 and Rag3c were used to fine map these two genes. The availability of the Illumina Infinium SoySNP50K/8K iSelect BeadChip, combined with single-nucleotide polymorphism (SNP) markers discovered through the whole-genome re-sequencing of E12901, facilitated the fine mapping process. Rag6 was refined to a 49-kb interval on chromosome 8 with four candidate genes, including three clustered nucleotide-binding site leucine-rich repeat (NBS-LRR) genes and an amine oxidase encoding gene. Rag3c was refined to a 150-kb interval on chromosome 16 with 11 candidate genes, two of which are a LRR gene and a lipase gene. Moreover, by sequencing the whole-genome exome-capture of the resistant source (E12901), structural variants were identified in the exons of the candidate genes of Rag6 and Rag3c. The closely linked SNP markers and the candidate gene information presented in this study will be significant resources for integrating Rag6 and Rag3c into elite cultivars and for future functional genetics studies.

Molecular Soybean-Pathogen Interactions

Soybean hosts a wide variety of pathogens that cause significant yield losses. The importance of soybean as a major oilseed crop has led to research focused on its interactions with pathogens, such as Soybean mosaic virus, Pseudomonas syringae, Phytophthora sojae, Phakopsora pachyrhizi, and Heterodera glycines. Pioneering work on soybean's interactions with these organisms, which represent the five major pathogen groups (viruses, bacteria, oomycetes, fungi, and nematodes), has contributed to our understanding of the molecular mechanisms underlying virulence and immunity. These mechanisms involve conserved and unique features that validate the need for research in both soybean and homologous model systems. In this review, we discuss identification of effectors and their functions as well as resistance gene-mediated recognition and signaling. We also point out areas in which model systems and recent advances in resources and tools have provided opportunities to gain deeper insights into soybean-pathogen interactions.

A pigeonpea gene confers resistance to Asian soybean rust in soybean

Asian soybean rust (ASR), caused by the fungus Phakopsora pachyrhizi, is one of the most economically important crop diseases, but is only treatable with fungicides, which are becoming less effective owing to the emergence of fungicide resistance. There are no commercial soybean cultivars with durable resistance to P. pachyrhizi, and although soybean resistance loci have been mapped, no resistance genes have been cloned. We report the cloning of a P. pachyrhizi resistance gene CcRpp1 (Cajanus cajan Resistance against Phakopsora pachyrhizi 1) from pigeonpea (Cajanus cajan) and show that CcRpp1 confers full resistance to P. pachyrhizi in soybean. Our findings show that legume species related to soybean such as pigeonpea, cowpea, common bean and others could provide a valuable and diverse pool of resistance traits for crop improvement.

A novel Phakopsora pachyrhizi resistance allele (Rpp) contributed by PI 567068A

KEY MESSAGE

The Rpp6 locus of PI 567102B was mapped from 5,953,237 to 5,998,461 bp (chromosome 18); and a novel allele at the Rpp6 locus or tightly linked gene Rpp[PI567068A] of PI 567068A was mapped from 5,998,461 to 6,160,481 bp. Soybean rust (SBR), caused by the obligate, fungal pathogen Phakopsora pachyrhizi is an economic threat to soybean production, especially in the Americas. Host plant resistance is an important management strategy for SBR. The most recently described resistance to P. pachyrhizi (Rpp) gene is Rpp6 contributed by PI 567102B. Rpp6 was previously mapped to an interval of over four million base pairs on chromosome 18. PI 567068A was recently demonstrated to possess a resistance gene near the Rpp6 locus, yet PI 567068A gave a differential isolate reaction to several international isolates of P. pachyrhizi. The goals of this research were to fine map the Rpp6 locus of PI 567102B and PI 567068A and determine whether or not PI 567068A harbors a novel Rpp6 allele or another allele at a tightly linked resistance locus. Linkage mapping in this study mapped Rpp6 from 5,953,237 to 5,998,461 bp (LOD score of 58.3) and the resistance from PI 567068A from 5,998,461 to 6,160,481 bp (LOD score of 4.4) (Wm82.a1 genome sequence). QTL peaks were 139,033 bp apart from one another as determined by the most significant SNPs in QTL mapping. The results of haplotype analysis demonstrated that PI 567102B and PI 567068A share the same haplotype in the resistance locus containing both Rpp alleles, which was designated as the Rpp6/Rpp[PI567068A] haplotype. The Rpp6/Rpp[PI567068A] haplotype identified in this study can be used as a tool to rapidly screen other genotypes that possess a Rpp gene(s) and detect resistance at the Rpp6 locus in diverse germplasm.

Gaining insight into soybean defense responses using functional genomics approaches

Soybean pathogens significantly impact yield, resulting in over $4 billion dollars in lost revenue annually in the United States. Despite the deployment of improved soybean cultivars, pathogens continue to evolve to evade plant defense responses. Thus, there is an urgent need to identify and characterize gene networks controlling defense responses to harmful pathogens. In this review, we focus on major advances that have been made in identifying the genes and gene networks regulating defense responses with an emphasis on soybean-pathogen interactions that have been amenable to gene function analyses using gene silencing technologies. Further we describe new research striving to identify genes involved in durable broad-spectrum resistance. Finally, we consider future prospects for functional genomic studies in soybean and demonstrate that understanding soybean disease and stress tolerance will be expedited at an unprecedented pace.