Next Article in Journal
Evaluating the Prebiotic Properties of Agar Oligosaccharides Obtained from the Red Alga Gracilaria fisheri via Enzymatic Hydrolysis
Next Article in Special Issue
Characterization of the Regulatory Network under Waterlogging Stress in Soybean Roots via Transcriptome Analysis
Previous Article in Journal
Phytochemical Profiling and Antioxidant Capacity of Traditional Plants, Northern Thailand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soybean Variety Saedanbaek Confers a New Resistance Allele to Phytophthora sojae

1
Department of Crop Science, College of Agriculture and Life Sciences, Chungnam National University, Daejeon 34134, Republic of Korea
2
Division of Crop Cultivation and Environment Research, Department of Central Area Crop Science, National Institute of Crop Science, Suwon 16613, Republic of Korea
3
Department of Crop Science and Biotechnology, College of Bioresource Science, Dankook University, Cheonan 31116, Republic of Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2023, 12(23), 3957; https://doi.org/10.3390/plants12233957
Submission received: 4 October 2023 / Revised: 13 November 2023 / Accepted: 21 November 2023 / Published: 24 November 2023

Abstract

:
Phytophthora root and stem rot (PRSR) disease results in substantial losses in soybean production worldwide. The occurrence of PRSR caused by Phytophthora sojae Kaufmann & Gerdemann has become increasingly important for soybean production in the Republic of Korea, but domestic soybean–P. sojae interaction has been less studied. The disease has been managed by developing varieties harboring resistance to the Phytophthora sojae (Rps) gene. The present study aimed to identify a major gene locus conferring resistance to new P. sojae isolate 2858 in the recombinant inbred line population derived from a cross between parental lines ‘Daepung’ (susceptible) and ‘Saedanbaek’ (resistant). Seventy-three recombination inbred lines (RILs) were evaluated for resistance to P. sojae isolate 2858. A resistance locus was identified in the approximate 3.3–4.3 megabase pair region on chromosome 3 using both single-marker and linkage analyses. The Rps of Saedanbaek (RpsSDB) was located on the well-known Rps gene/allele cluster region, which also partially overlapped with a locus previously identified in the Korean soybean variety, ‘Daewon’, resistant to another P. sojae isolate 2457 (RpsDW). Approximately 402 kilobase pairs of the interval region overlapped, including six nucleotide-binding site-leucine-rich repeat (NBS-LRR)-coding genes. Additional phenotypic assays revealed that Saedanbaek was susceptible to isolate 2457 and that Daewon was susceptible to isolate 2858, indicating that RpsSDB and RpsDW are different genes or alleles that confer race-specific resistance to the two P. sojae isolates. These results provide information that will be helpful for breeders developing P. sojae-resistant cultivars.

1. Introduction

Soybean [Glycine max (L.) Merr.] is one of the most economically valuable crops and an important source of protein and oil worldwide. Soybean production is affected by biotic and abiotic stressors. Phytophthora root and stem rot (PRSR), caused by the soil-borne pathogen, Phytophthora sojae Kaufmann & Gerdemann, is one of the most prevalent diseases leading to yield loss in soybean production [1]. Approximately 11 million tons of annual yield loss of soybean due to diseases was estimated in the United States and Canada from 2010 to 2019 [1,2]. Soybean is the primary host of the causal agent P. sojae, known to have a narrow host range [3,4]. It is a self-fertile oomycete that usually overwinters as oospores. Sporangia develop when soils are saturated, in which many zoospores will form. Once zoospores are released in soil water, they are attracted to soybean roots by chemotaxis [5]. It eventually encysts, germinates, and invades the root [4]. The symptoms of PRSR in the susceptible soybean include seed decay and seeding damping-off during the early growth stage, and brown stem lesion and collapsing tissue at the later season.
Soybean–P. sojae interaction has been most extensively studied in the U.S. in the past several decades [3]. Genetic diversity in the P. sojae population is distinguishable among multiple states of the U.S. and the frequency changes of diverse pathotypes were discovered through periodical state-wide studies over decades [3]. Similarly, the diversity and shift of pathotypes of the P. sojae population vary by country and by geographical region within a country based on the phenotypic assays using differential varieties and molecular analyses in China, Japan, and Brazil [6,7,8]. The development of resistant cultivars is the most effective way of preventing PRSR [7]. PRSR is mainly controlled by a single qualitative gene, the so-called resistance to the P. sojae (Rps) gene, and quantitative resistance has been reported as an alternative regulatory mechanism [3]. Rps-mediated resistance shows race-specific features, which are determined by the gene-for-gene relationship between the Rps of soybean and avirulence (Avr) gene of P. sojae [9,10]. Thus, identification of Rps loci and race-specific resistance alleles is essential for developing resistant varieties of P. sojae. Particular Rps genes were preferred, selected, and employed in newly developed resistant soybean varieties based on the genetic diversity of the local P. sojae population [7,11].
To date, more than 30 Rps alleles have been identified on the 10 chromosomes of soybean using bi-parental genetic mapping [12,13]. These Rps loci are not evenly distributed in the genome but tend to cluster in a few genomic regions [12,14]. Among these loci, approximately 20 Rps alleles have been mapped to the R-gene-rich region on chromosome 3 (Rps1a, 1b, 1c, 1d, 1k, 7, 9, UN1, Yu25, YD25, YD29, HN, Q, WY, HC18, X, GZ, unnamed in cv. Waseshiroge and Daewon, T1/T2/T3, and 14) [15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. These alleles were mapped to a range of regions harboring many nucleotide-binding site-leucine-rich repeats (NBS-LRRs or LRR) and serine/threonine protein kinase (STK) genes, which are primarily known as disease resistance genes in plant species [32,33]. Other loci were mapped on chromosome 2 (RpsZS18), 7 (Rps11), 10 (RpsSu), 13 (Rps3a, 3b, 3c, 8, SN10, and CD), 16 (Rps2 and UN2), 17 (Rps10), 18 (Rps4, 5, 6, JS, 12, and 13), and 19 (RpsYB30) [16,18,34,35,36,37,38,39,40,41,42,43,44,45,46].
In the Republic of Korea, PRSR was initially reported in the Chungchengnam-do in the late 1990s [47]; however, only a few studies were conducted in recent years. Twenty soybean cultivars widely planted in the Republic of Korea were screened to identify Rps-resistant cultivars against four isolates (P-9662, P-98145, 2457, and 3444-1) [48]. Among the cultivars, ten showed susceptible reactions to each isolate, and isolate-specific resistance reactions were observed in 11 varieties. ‘Daewon’ shows Rps-mediated resistance to isolate 2457 and the resistance gene was detected on chromosome 3 using linkage analysis of a recombination inbred line (RIL) population [24]. Another resistant variety, Socheong2, displayed resistance to P. sojae isolates 40412 and 2457 [49]. Linkage mapping revealed three Rps loci on chromosomes 3 and 18, and two at distant genomic locations on chromosome 18 [49].
Soil-borne pathogens such as P. sojae and Fusarium spp. occur more frequently under wet soil conditions [3,4]. PRSR is more frequently observed during the growing season because of increased soybean cultivation in converted paddy fields, where soil moisture levels tend to be high, owing to poor drainage. Consequently, PRSR has become a serious problem for soybean production in the Republic of Korea [48]. A preliminary screening showed that the cultivar ‘Saedanbaek’ is resistant to P. sojae isolate 2858, which has not yet been genetically characterized for the resistance. Thus, this study aimed to identify a genomic region associated with R-gene-mediated resistance to P. sojae isolate 2858 using an RIL population derived from a cross between ‘Daepung’ and Saedanbaek.

2. Results

2.1. Phenotypic Assay for Resistance to P. sojae

The two parents, Daepung and Saedanbaek, showed typical susceptibility and resistance reactions to P. sojae isolate 2858 in repeated experiments (Figure 1). Of the RILs, 25 and 43 exhibited resistance and susceptibility, respectively (Table 1). Intermediate reactions were observed in four RILs, with 45–60% dead seedlings. A goodness-of-fit test demonstrated that the segregation ratio agreed marginally with the expected ratio in the F6 generation (34.9:2.2:34.9 for resistance: intermediate: susceptibility), suggesting that the trait was mediated by a single gene (χ2 = 6.05, p = 0.05) (Table 1).

2.2. Single-Marker Analysis of Variance (ANOVA)

Singe-marker ANOVA was performed to detect significant marker-trait associations for 26,130 SNP. In total, 28 SNP markers were highly significant and located in a genomic region of 3.3–4.3 megabase pairs (Mbp) on chromosome 3, for which p-values ranged from 3.7 × 10−16 to 4.4 × 10−23 (Table 2). The estimated R2 ranged from 71% to 82%, indicating that the resistance was controlled by a single major gene at this locus (Table 2). No SNP was significant on the other chromosomes based on the threshold of p < 3.8 × 10−7.

2.3. Identification of a Resistance Locus Using Linkage Analysis

In total, 963 SNPs evenly distributed over 20 chromosomes were retained after binning 26,130 SNPs, which were subsequently used to construct genetic maps. The number of SNPs that were integrated into each chromosome ranged from 33 to 63, with an average of 48 (Table 3). The total length of the genetic map was approximately 1586 cM, with an average marker interval of 1.7 cM (Table 3). In the genetic bin map, the genetic positions of the mapped SNPs were in good agreement with the physical positions based on the reference, Glyma.Wm82.a2 (Figure S1). Linkage analysis identified the same genomic region (3.3–4.3 Mbp) with an extremely high logarithm of odds (LOD) of 38.7, explaining 93% of the phenotypic variance (Table 4, Figure 2). The peak of LOD was located at the AX-90404495 on the genetic bin map, which represented all the 28 significant SNPs detected by ANOVA. The additive effect of the locus was 45.5, indicating that the Saedanbaek allele conferred resistance to P. sojae (Table 4).

2.4. Comparison of the Rps Locus of Saedanbaek (RpsSDB) with That of Daewon (RpsDW)

The resistant allele of Saedanbaek identified in the present study was named RpsSDB to distinguish it from other colocalized Rps alleles (Figure 3). RpsSDB partially overlapped with RpsDW detected in ‘Daewon’, which was previously mapped on chromosome 3 [24]. To distinguish between the two Rps loci, recombinants within the target interval were screened using RIL populations. As RILs with recombination breaks within the RpsSDB were not detected in the Saedanbaek population, RpsSDB and RpsDW overlapped in an approximately 402 kbp interval (3,893,390–4,295,128 bp) (Figure 3). According to the soybean reference genome (Glyma.Wm82.a2.v1), 53 and 31 genes, in total, were annotated in the interval of RpsSDB and RpsDW intervals, including seven and eight copies of NBS-LRR- or STK-coding genes, the two main categories of disease resistance genes, respectively (Table S1). The overlapped interval included six LRR-coding genes, i.e., Glyma.03G034200, Glyma.03G034400, Glyma.03G034500, Glyma.03G034800, Glyma.03G034900, and Glyma.03G035300 (Figure 3, Table S1).

2.5. Distinguished Race-Specificity of RpsSDB and RpsDW

Saedanbaek and Daewon were inoculated with isolates 2858 and 2457, and their race specificity was confirmed. Saedanbaek was resistant to isolate 2858 but susceptible to isolate 2457 (Figure 4). Daewon showed susceptibility and resistance to isolates 2858 and 2457, respectively (Figure 4). To confirm the reciprocal race-specific reactions of RpsSDB and RpsDW, 25 RILs harboring the Saedanbaek-resistant allele on chromosome 3 and 35 RILs that included the Daewon-resistant allele on chromosome 3 were inoculated with isolates 2858 and 2457, respectively. Daepung was susceptible to these two isolates. Saedanbaek and 25 RILs that were resistant to isolate 2858 were susceptible to isolate 2457 (Table 5). Similarly, Daewon and 35 RILs that were resistant to isolate 2457 were confirmed to be susceptible to isolate 2858 (Table 5). As all the tested RILs and the two resistant parents exhibited susceptibility to the tested isolates with 100% seedling death, race-specific resistance of the respective RpsSDB and RpsDW was confirmed, and either one could effectively protect against only one isolate.

3. Discussion

Characterization of valuable genetic sources for a target trait is the most important initial step in plant breeding. Previously, few Korean varieties have exhibited resistance to some selective P. sojae isolates, and two varieties, Daewon and Socheong2, have been genetically characterized for their inheritance of resistance to P. sojae isolates 2457 and 40412, respectively [24,48,49]. In the present study, the Korean high-protein variety, Saedanbaek, was used to analyze its Rps-mediated resistance to P. sojae isolate 2858 using the Daepung × Saedanbaek RIL population. Soybean resistance to isolate 2858 was investigated for the first time in this study. The majority of RILs were clearly segregated into either resistant (R) or sensitive (S) phenotypes following hypocotyl inoculation, and a single resistance gene locus (RpsSDB) was mapped to the 3.3–4.3 Mbp region of chromosome 3 via both single-marker ANOVA and linkage analysis (Figure 2, Table 4). This locus explained approximately 82% and 93% of the phenotypic variation with ANOVA and composite interval mapping, respectively, indicating that the resistance was controlled by a single major gene in a qualitative manner (Figure 2, Table 4). These high R2 or PVE values of the identified SNP strongly support that distinct phenotypes are clearly segregated in the mapping population and a strong association between phenotypic and genotypic data is quite evident. Also, the size of the mapping population (i.e., N = 72) is sufficient to identify a major locus controlling resistance to P. sojae.
The distribution of Rps genes appears to be uneven over 20 chromosomes; more than 20 Rps genes have been identified on chromosome 3 from different resistance sources [12,13,14]. This is comparable with another ‘hot spot’ of Rps, the approximately 53−56 Mbp region on chromosome 18, where Rps4, 5, 6, 12, JS, and unnamed in cv. CheonAl partially overlapped or were physically separated based on their approximate positions [12,36,38,39,40,49,50]. Some of the Rps alleles in these clusters were roughly mapped by genetic analyses with small numbers of molecular markers and/or small sizes of mapping populations without further fine-scale mapping. Thus, it is unclear whether they all are indeed allelic to each other or belong to adjacent different loci. In the present study, Saedanbaek and Daewon exhibited reciprocal race-specific reactions with isolates 2858 and 2457. However, the two presumably-adjacent Rps genes were not fully defined for their physical positions due to lack of a few RILs with recombination breakpoints within the target interval, which facilitate a narrowing-down of the target interval as used in the fine mapping of RpsWY and RpsGZ [15,28]. Using the 180 K SNP data, haplotypes of this overlap were compared among the three parental varieties, though any feature determining the allelic relationship was not discovered. Given this results, genetic allelism tests should be performed using the F2:3 Daewon × Saedanbaek population and the two P. sojae isolates to confirm the allelic relationship between RpsSDB and RpsDW. Deep-depth sequencing analysis should assist to elucidate the relationship of the two Rps genes for further study. More accurate determination of these genes/alleles will aid in the marker-assisted selection for multiple Rps-stacked lines.
The NBS-LRR genes are the extremely large family of plant disease resistance genes and their clustered/tandemly-repeated nature or the evolutionary complexity has been well documented [32,33]. According to the gene annotation (Glym.Wm82.a2.v1), multiple copies of NBS-LRR- and STK-coding genes were commonly distributed in the two hotspot regions mentioned above. The intervals in the respective RpsDW and unnamed Rps of CheonAl included 14 and 9 genes [24,50]. Another R-gene RpsUN1 was fine-mapped within a 151 kbp interval that harbored five candidate genes, including three predicted R-gene homologs [18,51]. RpsHN was also mapped to a 278 kbp region, where eight genes were annotated, including one STK gene [19]. The overlapping 403 kbp region between RpsSDB and RpsDW includes six NBS-LRR genes (Figure 3). Of the six genes, Glyma.03G034900 has been previously highlighted as a candidate gene for P. sojae resistance in RpsHN and RpsYD25 [19,52]. The expression of this gene increased 12 h after P. sojae inoculation in Meng8206 (resistant), whereas it did not significantly change in LinMeng6-46 (susceptible) [19]. Nearby, slightly outside the overlapping region, Glyma.03G034600 was also identified as a candidate gene for RpsUN1 based on fine mapping and gene expression analysis [51]. These Rps clusters will be interesting resources for studying the dynamics of the R-genes and understanding them will provide insights regarding the mechanism of P. sojae resistance.
The most significant marker, AX-90404495, was mapped approximately 25 and 20 cM from two flanking markers. As noted above, this is not a single SNP, but is representative of 28 co-segregating SNPs located over 3.3–4.3 Mbp. Considering the genomic distances to the flanking SNPs based on the reference genome (Glyma.Wm82.a2), the distances between the markers were notably large, implying that the recombination frequency might be high around this R-gene locus (Figure S1). Interestingly, high recombination frequencies were observed in the similar genomic region of chromosome 3 from Shuurei × Tosan-231 population [30]. The possible explanations for this result are as follows: (i) the small size of the mapping population may lead to a relatively biased recombination frequency in this region [53], (ii) chromosomal rearrangement or structurally variable regions by genotype may exist between the representative U.S. variety (i.e., Williams82) and non-U.S. variety [54], or (iii) unequal crossing over, a driving force of R-gene evolution, may increase the variability of this region [55]. These factors can result in an overestimated recombinant frequency, leading to prolonged marker intervals.
Genomic rearrangements or recombination, such as gene duplication, unequal exchange, abnormal chromosomal recombination, and gene conversion, have been reported as important mechanisms that affect the copy number of disease-resistant genes [55]. In soybeans, Rps11, an NBS-LRR gene that was subjected to map-based cloning, was derived from rounds of unequal crossover events, resulting in promoter fusion and LRR expansion [56]. Consequently, this change led to broad-spectrum resistance to 127 P. sojae isolates widely distributed across the USA. Copy number variations in the Rhg1-b allele of soybean have also been associated with resistance to the soybean cyst nematode (SCN) [57]. The susceptible varieties for SCN have one copy of the 31 kbp segment, which encodes resistant proteins for SCN, whereas the rhg1-b haplotype harbors 10 tandem copies with increased resistance to SCN. The maize Rp1 locus is characterized as a complex region related to resistance to the fungus Puccinia sorghi, and some Rp1 alleles are meiotically unstable, showing unequal exchange in this region [58]. Unequal exchange in the Rp1 locus leads to variations in the repeat number of Rp1, resulting in altered disease resistance. These reports suggest that recombination events in repeated LRR genes may affect the number or structure of LRR genes, leading to changes in resistance and race-specific reactions.
The development of broad-spectrum disease-resistant varieties is an ideal goal of soybean breeding programs. Some loci have been associated with broad-spectrum resistance to several P. sojae isolates. RpsYD25 identified in the Chinese soybean cultivar Yudou25 exhibited broad-spectrum resistance to P. sojae isolates with different pathotypes in China [59]. Among the 26 P. sojae pathotypes, RpsYD25 showed resistance to 21 isolates with differing pathotypes. Another Rps locus, Rps11, in the soybean landrace PI 594527, showed broad-spectrum resistance to 16 P. sojae isolates [56]. Tightly linked or pyramided multiple Rps genes also facilitate to protect soybeans against a broad-range of P. sojae isolates [30,60]. As most of the tested Korean soybean varieties were susceptible to four P. sojae isolates, and some were resistant to only a few isolates [48], the pyramiding of multiple Rps alleles from different resistance sources is a desirable approach in developing new resistant varieties. For this, the identification of germplasms with resistance to diverse pathotypes should be prioritized.
In summary, a locus, named RpsSDB, resistant to P. sojae isolate 2858, was mapped in the 921 kbp region, which is a well-known Rps gene cluster harboring enriched NBS-LRR annotated genes, on chromosome 3 in the Daepung × Saedanbaek population. The physical position of RpsSDB partially overlapped with the previously detected RpsDW in the 402 kbp interval, which harbors six annotated NBS-LRR genes. Interestingly, RpsSDB and RpsDW reciprocally showed race-specific resistance to either isolate 2457 or 2858, implying that they may be two allelic forms of one gene or two different genes. The current genetic mapping and haplotype comparison did not differentiate the chromosomal positions of the two Rps from those different varieties. In future studies, high-throughput sequencing, relative expression analysis, or genetic tests for allelism using a cross of Saedanbaek × Daewon will be able to clarify whether RpsSDB and RpsDW are two different genes adjacently located in the genome. The results of this study will assist soybean breeders in developing soybean lines with broad spectrum resistance to P. sojae infection.

4. Materials and Methods

4.1. Plant Materials

Daepung, Saedanbaek, and 72 RILs derived from crosses between Daepung and Saedanbaek were used in this study. Daepung is an elite Korean cultivar with lodging tolerance [61], and Saedanbaek is widely used in making tofu because of its high protein content and resistance to bacterial pustules [62]. Both parental lines exhibited white flowers, yellow seeds, and determinate growth habits. In our preliminary experiments, two soybean varieties, Daepung and Saedanbaek, were classified as susceptible or resistant to P. sojae isolate 2858 (Figure 1). In 2019, the RIL population was obtained from Dankook University in Cheonan, Republic of Korea. To propagate the seeds, the F6-derived RIL population was grown in a research field at the Chungnam National University, Daejeon, Republic of Korea. In this study, two parents and 72 RILs were used in subsequent experiments.

4.2. Evaluation of Resistance to P. sojae Isolate 2858

The P. sojae isolate 2858 was originally collected in Gyeonggi-do and provided by the Department of Southern Area Crop Science, National Institute of Crop Science, Rural Development Administration, Miryang, Gyeongsangbuk-do, Republic of Korea. A single oospore of the isolate 2858 was re-isolated from the original isolate for uniformity of results and then it was stocked for the following phenotypic assays. The inoculum was prepared as previously described [63]. Briefly, P. sojae was incubated in V8 juice medium for 7 days. The V8 medium, which was fully covered by P. sojae mycelia, was macerated using a 50 mL syringe and transferred into a 10 mL syringe before inoculation. Ten to twelve seedlings of each genotype were grown in plastic pots and inoculated using the hypocotyl inoculation technique when they were 7 days old after sowing [63]. A short slit (about 1 cm) was made using a scalpel on the upper part of the hypocotyl, and 0.2–0.4 mL of a mycelial slurry was placed inside the slit. The inoculated plants were maintained under dark and humid conditions (>90% relative humidity, 25 °C) for 16 h, and moved into a plant growth chamber (14/10 h of light/dark condition at 25 °C with >70% relative humidity). Seven days after inoculation, the percentage of dead seedlings was calculated and (2)the reaction of each RIL was determined in a categorical manner based on the percentage of dead seedlings; resistant (R) for ≤20%, intermediate (I) for 20–80%, or susceptible (S) for ≥80%. The phenotypic assay was repeated thrice.

4.3. DNA Extraction and Single Nucleotide Polymorphism (SNP) Genotyping

The genomic DNA of parents and RILs was extracted using cetyltrimethylammonium bromide (CTAB) [64] and genotyping was performed using Axiom® 180K SoyaSNP arrays (Affymetrix, Santa Clara, CA, USA) [65]. The SNP data were preprocessed as follows: in total, 169,028 high-quality SNPs were obtained from the SNP array and 142,036 undefined or non-polymorphic SNPs were removed. Additionally, 862 SNPs were eliminated because of ≥10% missing. In total, 26,130 SNPs were included in the final dataset for subsequent genetic analyses.

4.4. Goodness-of-Fit Test and Single-Marker Analysis of Variance

The segregation ratio for P. sojae-resistant reactions was investigated in the RIL population and the Chi-square (χ2) test was used. Marker-trait associations were analyzed using single-marker ANOVA using the lme4 package 1.1.31 [66] in R-4.3.2 [67]. Bonferroni correction was used to determine a conservative threshold for p-value (i.e., p < 3.8 × 10−7) using the following equation: α/N, where α is 0.01 and N is the number of SNPs.

4.5. Construction of Genetic Map and Linkage Analysis

A genetic map of the Daepung × Saedanbaek RIL population was constructed using IciMapping version 4.2 following the users’ instructions [68]. Of the 26,130 SNPs, redundant SNPs were filtered using the BIN function. Filtered SNPs were grouped and ordered using the MAP function. Linkage analysis was conducted using the BIP function, where the average values of the percentage of dead seedlings over replicates were used as the phenotype. With mapping steps of 1 cM, a 1000-permutation test was performed to determine the LOD threshold at p < 0.05 [69]. Graphical presentations of the genetic map and QTL plots of LOD scores were produced using MapChart 2.32 [70].

4.6. Cross-Validation of the Allelic Relationship between Saedanbaek and Daewon

As described in Section 4.2, additional P. sojae inoculation assays were performed using the hypocotyl inoculation technique [63] to determine whether two different genes were located adjacently and conferred resistance to the respective isolate in each resistant variety. Saedanbaek, and 25 RILs derived from a cross between Daepung and Saedanbaek, which showed resistance to isolate 2858, were inoculated with isolate 2457. Similarly, Daewon and 35 RILs derived from a cross between the Daepung and Daewon populations that exhibited resistance to isolate 2457 were also tested using isolate 2858.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12233957/s1, Figure S1: Collinearity between genetic (cM) and physical (Mbp) positions of the mapped SNPs; Table S1. List of candidate genes in the RpsSDB and RpsDW region.

Author Contributions

Conceptualization, S.L.; investigation, H.J.Y. and S.L.; resources, S.K., J.-M.K. and I.-J.K.; data analysis and curation, H.J.Y. and S.L.; writing and original draft preparation, H.J.Y., K.-C.S. and S.L.; writing, review, and editing, K.-C.S., S.K., J.-M.K., I.-J.K. and S.L.; project administration, S.L.; and funding acquisition, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was in part funded by the “Cooperative Research Program for Agriculture Science and Technology” (Project No. PJ01416804 or RS-2019-RD008897) of the Rural Development Administration, Republic of Korea.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

The authors thank Yun-Woo Jang in the department of Southern Area Crop Science, National Institute of Crop Science, Rural Development Administration, Miryang, Republic of Korea, for providing the P. sojae isolate 2858.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allen, T.W.; Bradley, C.A.; Sisson, A.J.; Byamukama, E.; Chilvers, M.I.; Coker, C.M.; Collins, A.A.; Damicone, J.P.; Dorrance, A.E.; Dufault, N.S. Soybean yield loss estimates due to diseases in the United States and Ontario, Canada, from 2010 to 2014. Plant Health Prog. 2017, 18, 19–27. [Google Scholar] [CrossRef]
  2. Bradley, C.A.; Allen, T.W.; Sisson, A.J.; Bergstrom, G.C.; Bissonnette, K.M.; Bond, J.; Byamukama, E.; Chilvers, M.I.; Collins, A.A.; Damicone, J.P.; et al. Soybean yield loss estimates due to diseases in the United States and Ontario, Canada, from 2015 to 2019. Plant Health Prog. 2021, 22, 483–495. [Google Scholar] [CrossRef]
  3. Dorrance, A.E. Management of Phytophthora sojae of soybean: A review and future perspectives. Can. J. Plant Pathol. 2018, 40, 210–219. [Google Scholar] [CrossRef]
  4. Schmitthenner, A.F. Problems and progress in control of Phytophthora root rot of soybean. Plant Dis. 1985, 69, 362–368. [Google Scholar] [CrossRef]
  5. Morris, P.F.; Ward, E. Chemoattraction of zoospores of the soybean pathogen, Phytophthora sojae, by isoflavones. Physiol. Mol. Plant Pathol. 1992, 40, 17–22. [Google Scholar] [CrossRef]
  6. Zhang, Z.; Zhao, Y.; Yu, H.; Chen, Y.; Gu, X.; Wen, J. Pathotypes of Phytophthora sojae and their distribution in Jilin, China. J. Plant Pathol. 2021, 103, 241–248. [Google Scholar] [CrossRef]
  7. Sugimoto, T.; Kato, M.; Yoshida, S.; Matsumoto, I.; Kobayashi, T.; Kaga, A.; Hajika, M.; Yamamoto, R.; Watanabe, K.; Aino, M.; et al. Pathogenic diversity of Phytophthora sojae and breeding strategies to develop Phytophthora-resistant soybeans. Breed. Sci. 2012, 61, 511–522. [Google Scholar] [CrossRef]
  8. Batista, I.C.A.; Silva, M.P.C.; Silva Junior, A.L.; Arrigada, M.P.G.; de Camargo, M.P.; Figueiredo, A.; Júnior, B.T.H.; Mizubuti, E.S.G. A shift in pathotype diversity and complexity of Phytophthora sojae in Brazil. Plant Dis 2023, 107, 1968–1972. [Google Scholar] [CrossRef]
  9. Flor, H.H. Current status of the gene-for-gene concept. Annu. Rev. Phytopathol. 1971, 9, 275–296. [Google Scholar] [CrossRef]
  10. Dodds, P.N. From gene-for-gene to resistosomes: Flor’s enduring legacy. Mol. Plant-Microbe Interact. 2023, 36, 461–467. [Google Scholar] [CrossRef]
  11. Tremblay, V.; McLaren, D.L.; Kim, Y.M.; Strelkov, S.E.; Conner, R.L.; Wally, O.; Bélanger, R.R. Molecular assessment of pathotype diversity of Phytophthora sojae in Canada highlights declining sources of resistance in soybean. Plant Dis. 2021, 105, 4006–4013. [Google Scholar] [CrossRef] [PubMed]
  12. Jang, I.-H.; Lee, S. A review and perspective on soybean (Glycine max L.) breeding for the resistance to Phytophthora sojae in Korea. Plant Breed. Biotechnol. 2020, 8, 114–130. [Google Scholar] [CrossRef]
  13. Chandra, S.; Choudhary, M.; Bagaria, P.K.; Nataraj, V.; Kumawat, G.; Choudhary, J.R.; Sonah, H.; Gupta, S.; Wani, S.H.; Ratnaparkhe, M.B. Progress and prospectus in genetics and genomics of Phytophthora root and stem rot resistance in soybean (Glycine max L.). Front. Genet. 2022, 13, 939182. [Google Scholar] [CrossRef] [PubMed]
  14. Lin, F.; Chhapekar, S.S.; Vieira, C.C.; Da Silva, M.P.; Rojas, A.; Lee, D.; Liu, N.; Pardo, E.M.; Lee, Y.-C.; Dong, Z.; et al. Breeding for disease resistance in soybean: A global perspective. Theor. Appl. Genet. 2022, 135, 3773–3872. [Google Scholar] [CrossRef] [PubMed]
  15. Cheng, Y.; Ma, Q.; Ren, H.; Xia, Q.; Song, E.; Tan, Z.; Li, S.; Zhang, G.; Nian, H. Fine mapping of a Phytophthora-resistance gene RpsWY in soybean (Glycine max L.) by high-throughput genome-wide sequencing. Theor. Appl. Genet. 2017, 130, 1041–1051. [Google Scholar] [CrossRef] [PubMed]
  16. Demirbas, A.; Rector, B.; Lohnes, D.; Fioritto, R.; Graef, G.; Cregan, P.; Shoemaker, R.; Specht, J. Simple sequence repeat markers linked to the soybean Rps genes for Phytophthora resistance. Crop Sci. 2001, 41, 1220–1227. [Google Scholar] [CrossRef]
  17. Li, Y.; Sun, S.; Zhong, C.; Wang, X.; Wu, X.; Zhu, Z. Genetic mapping and development of co-segregating markers of RpsQ, which provides resistance to Phytophthora sojae in soybean. Theor. Appl. Genet. 2017, 130, 1223–1233. [Google Scholar] [CrossRef]
  18. Lin, F.; Zhao, M.; Ping, J.; Johnson, A.; Zhang, B.; Abney, T.S.; Hughes, T.J.; Ma, J. Molecular mapping of two genes conferring resistance to Phytophthora sojae in a soybean landrace PI 567139B. Theor. Appl. Genet. 2013, 126, 2177–2185. [Google Scholar] [CrossRef]
  19. Niu, J.; Guo, N.; Sun, J.; Li, L.; Cao, Y.; Li, S.; Huang, J.; Zhao, J.; Zhao, T.; Xing, H. Fine mapping of a resistance gene RpsHN that controls Phytophthora sojae using recombinant inbred lines and secondary populations. Front. Plant Sci. 2017, 8, 538. [Google Scholar] [CrossRef]
  20. Sugimoto, T.; Yoshida, S.; Watanabe, K.; Aino, M.; Kanto, T.; Maekawa, K.; Irie, K. Identification of SSR markers linked to the Phytophthora resistance gene Rps1-d in soybean. Plant Breed. 2007, 127, 154–159. [Google Scholar] [CrossRef]
  21. Sun, S.; Wu, X.; Zhao, J.; Wang, Y.; Tang, Q.; Yu, D.; Gai, J.; Xing, H. Characterization and mapping of RpsYu25, a novel resistance gene to Phytophthora sojae. Plant Breed. 2011, 130, 139–143. [Google Scholar] [CrossRef]
  22. Weng, C.; Yu, K.; Anderson, T.; Poysa, V. Mapping genes conferring resistance to Phytophthora root rot of soybean, Rps1a and Rps7. J. Hered. 2001, 92, 442–446. [Google Scholar] [CrossRef]
  23. Wu, X.-L.; Zhang, B.-Q.; Sun, S.; Zhao, J.-M.; Yang, F.; Guo, N.; Gai, J.-Y.; Xing, H. Identification, genetic analysis and mapping of resistance to Phytophthora sojae of Pm28 in Soybean. Agric. Sci. China 2011, 10, 1506–1511. [Google Scholar] [CrossRef]
  24. Jang, I.-H.; Kang, I.J.; Kim, J.-M.; Kang, S.-T.; Jang, Y.E.; Lee, S. Genetic mapping of a resistance locus to Phytophthora sojae in the Korean soybean cultivar Daewon. Plant Pathol. J. 2020, 36, 591–599. [Google Scholar] [CrossRef]
  25. Zhang, J.; Xia, C.; Wang, X.; Duan, C.; Sun, S.; Wu, X.; Zhu, Z. Genetic characterization and fine mapping of the novel Phytophthora resistance gene in a Chinese soybean cultivar. Theor. Appl. Genet. 2013, 126, 1555–1561. [Google Scholar] [CrossRef] [PubMed]
  26. Zhong, C.; Sun, S.; Li, Y.; Duan, C.; Zhu, Z. Next-generation sequencing to identify candidate genes and develop diagnostic markers for a novel Phytophthora resistance gene, RpsHC18, in soybean. Theor. Appl. Genet. 2017, 131, 525–538. [Google Scholar] [CrossRef] [PubMed]
  27. Zhong, C.; Li, Y.; Sun, S.; Duan, C.; Zhu, Z. Genetic mapping and molecular characterization of a broad-spectrum Phytophthora sojae resistance gene in Chinese soybean. Int. J. Mol. Sci. 2019, 20, 1809. [Google Scholar] [CrossRef]
  28. Jiang, B.; Cheng, Y.; Cai, Z.; Li, M.; Jiang, Z.; Ma, R.; Yuan, Y.; Xia, Q.; Nian, H. Fine mapping of a Phytophthora-resistance locus RpsGZ in soybean using genotyping-by-sequencing. BMC Genet. 2020, 21, 280. [Google Scholar] [CrossRef]
  29. Sugimoto, T.; Yoshida, S.; Kaga, A.; Hajika, M.; Watanabe, K.; Aino, M.; Tatsuda, K.; Yamamoto, R.; Matoh, T.; Walker, D. Genetic analysis and identification of DNA markers linked to a novel Phytophthora sojae resistance gene in the Japanese soybean cultivar Waseshiroge. Euphytica 2011, 182, 133. [Google Scholar] [CrossRef]
  30. Matsuoka, J.I.; Takahashi, M.; Yamada, T.; Kono, Y.; Yamada, N.; Takahashi, K.; Moriwaki, J.; Akamatsu, H. Identification of three closely linked loci conferring broad-spectrum Phytophthora sojae resistance in soybean variety Tosan-231. Theor. Appl. Genet. 2021, 134, 2151–2165. [Google Scholar] [CrossRef]
  31. Chen, L.; Wang, W.; Ping, J.; Fitzgerald, J.C.; Cai, G.; Clark, C.B.; Aggarwal, R.; Ma, J. Identification and molecular mapping of Rps14, a gene conferring broad-spectrum resistance to Phytophthora sojae in soybean. Theor. Appl. Genet. 2021, 134, 3863–3872. [Google Scholar] [CrossRef] [PubMed]
  32. Van Wersch, S.; Li, X. Stronger when together: Clustering of plant NLR disease resistance genes. Trends Plant Sci. 2019, 24, 688–699. [Google Scholar] [CrossRef] [PubMed]
  33. Tamborski, J.; Krasileva, K.V. Evolution of plant NLRs: From natural history to precise modifications. Annu. Rev. Plant Biol. 2020, 71, 355–378. [Google Scholar] [CrossRef]
  34. Burnham, K.; Dorrance, A.; Francis, D.; Fioritto, R.; St Martin, S. Rps 8, a new locus in soybean for resistance to Phytophthora sojae. Crop Sci. 2003, 43, 101–105. [Google Scholar]
  35. Gordon, S.G.; Martin, S.K.S.; Dorrance, A.E. Rps 8 maps to a resistance gene rich region on soybean molecular linkage group F. Crop Sci. 2006, 46, 168–173. [Google Scholar] [CrossRef]
  36. Gordon, S.; Kowitwanich, K.; Pipatpongpinyo, W.; Martin, S.S.; Dorrance, A. Molecular marker analysis of soybean plant introductions with resistance to Phytophthora sojae. Phytopathology 2007, 97, 113–118. [Google Scholar] [CrossRef] [PubMed]
  37. Yu, A.; Xu, P.; Wang, J.; Zhang, S.; Wu, J.; Li, W.; Chen, W.; Li, N.; Fan, S.; Wang, X. Genetic analysis and SSR mapping of gene resistance to Phytophthora sojae race 1 in soybean cv Suinong 10. Chin. J. Oil Crop Sci. 2010, 32, 462–466. [Google Scholar]
  38. Sahoo, D.K.; Abeysekara, N.S.; Cianzio, S.R.; Robertson, A.E.; Bhattacharyya, M.K. A novel Phytophthora sojae resistance Rps12 gene mapped to a genomic region that contains several Rps genes. PLoS ONE 2017, 12, e0169950. [Google Scholar] [CrossRef]
  39. Sandhu, D.; Gao, H.; Cianzio, S.; Bhattacharyya, M.K. Deletion of a disease resistance nucleotide-binding-site leucine-rich-repeat-like sequence is associated with the loss of the Phytophthora resistance gene Rps4 in soybean. Genetics 2004, 168, 2157–2167. [Google Scholar] [CrossRef]
  40. Sun, J.; Li, L.; Zhao, J.; Huang, J.; Yan, Q.; Xing, H.; Guo, N. Genetic analysis and fine mapping of RpsJS, a novel resistance gene to Phytophthora sojae in soybean [Glycine max (L.) Merr.]. Theor. Appl. Genet. 2014, 127, 913–919. [Google Scholar] [CrossRef]
  41. Huang, J.; Guo, N.; Li, Y.; Sun, J.; Hu, G.; Zhang, H.; Li, Y.; Zhang, X.; Zhao, J.; Xing, H.; et al. Phenotypic evaluation and genetic dissection of resistance to Phytophthora sojae in the Chinese soybean mini core collection. BMC Genet. 2016, 17, 85. [Google Scholar] [CrossRef] [PubMed]
  42. Yao, H.-Y.; Xiaoming, W.; Wu, X.-F.; Xiao, Y.-N.; Zhu, Z.-D. Molecular mapping of Phytophthora resistance gene in Soybean cultivar Zaoshu18. J. Plant Genet. Resour. 2010, 11, 213–217. [Google Scholar]
  43. Ping, J.; Fitzgerald, J.C.; Zhang, C.; Lin, F.; Bai, Y.; Wang, D.; Aggarwal, R.; Rehman, M.; Crasta, O.; Ma, J. Identification and molecular mapping of Rps11, a novel gene conferring resistance to Phytophthora sojae in soybean. Theor. Appl. Genet. 2016, 129, 445–451. [Google Scholar] [CrossRef] [PubMed]
  44. Wu, X.; Zhou, B.; Sun, S.; Zhao, J.; Chen, S.; Gai, J.; Xing, H. Genetic analysis and mapping of resistance to Phytophthora sojae of Pm14 in soybean. Sci. Agric. Sin. 2011, 44, 456–460. [Google Scholar]
  45. Zhang, J.; Xia, C.; Duan, C.; Sun, S.; Wang, X.; Wu, X.; Zhu, Z. Identification and candidate gene analysis of a novel Phytophthora resistance gene Rps10 in a Chinese soybean cultivar. PLoS ONE 2013, 8, e69799. [Google Scholar] [CrossRef]
  46. Zhu, Z.; Huo, Y.; Wang, X.; Huang, J.; Wu, X. Molecular identification of a novel Phytophthora resistance gene in soybean. Acta Agron. Sin. 2007, 33, 154–157. [Google Scholar]
  47. Jee, H.; Kim, W.; Cho, W. Occurrence of Phytophthora root rot on soybean (Glycine max) and identification of the causal fungus. RDA J. Crop Prot. 1998, 40, 16–22. [Google Scholar]
  48. Kang, I.J.; Kang, S.; Jang, I.H.; Jang, Y.W.; Shim, H.K.; Heu, S.; Lee, S. Identification of new isolates of Phytophthora sojae and the reactions of Korean soybean cultivars following hypocotyl inoculation. Plant Pathol. J. 2019, 35, 698. [Google Scholar] [CrossRef]
  49. Jang, Y.E.; Jang, I.H.; Kang, I.J.; Kim, J.-M.; Kang, S.-T.; Lee, S. Two isolate-specific resistance loci for Phytophthora sojae in the soybean Socheong2. Korean J. Breed. Sci. 2020, 52, 398–407. [Google Scholar] [CrossRef]
  50. You, H.J.; Kang, E.J.; Kang, I.J.; Kim, J.-M.; Kang, S.-T.; Lee, S. Identification of a locus associated with resistance to phytophthora sojae in the soybean elite line ‘CheonAl’. Korean J. Crop Sci. 2023, 68, 134–146. [Google Scholar] [CrossRef]
  51. Li, L.; Lin, F.; Wang, W.; Ping, J.; Fitzgerald, J.C.; Zhao, M.; Li, S.; Sun, L.; Cai, C.; Ma, J. Fine mapping and candidate gene analysis of two loci conferring resistance to Phytophthora sojae in soybean. Theor. Appl. Genet. 2016, 129, 2379–2386. [Google Scholar] [CrossRef]
  52. Zhong, C.; Sun, S.; Zhang, X.; Duan, C.; Zhu, Z. Fine mapping, candidate gene identification and co-segregating marker development for the Phytophthora root rot resistance gene RpsYD25. Front. Genet. 2020, 11, 799. [Google Scholar] [CrossRef]
  53. Wang, Y.; Sun, S.; Liu, B.; Wang, H.; Deng, J.; Liao, Y.; Wang, Q.; Cheng, F.; Wang, X.; Wu, J. A sequence-based genetic linkage map as a reference for Brassica rapa pseudochromosome assembly. BMC Genet. 2011, 12, 239. [Google Scholar] [CrossRef] [PubMed]
  54. Million, C.R.; Wijeratne, S.; Cassone, B.J.; Lee, S.; Rouf Mian, M.A.; McHale, L.K.; Dorrance, A.E. Hybrid genome assembly of a major quantitative disease resistance locus in soybean toward Fusarium graminearum. Plant Genome 2019, 12, 180102. [Google Scholar] [CrossRef] [PubMed]
  55. Michelmore, R.W.; Meyers, B.C. Clusters of resistance genes in plants evolve by divergent selection and a birth-and-death process. Genome Res. 1998, 8, 1113–1130. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, W.; Chen, L.; Fengler, K.; Bolar, J.; Llaca, V.; Wang, X.; Clark, C.B.; Fleury, T.J.; Myrvold, J.; Oneal, D.; et al. A giant NLR gene confers broad-spectrum resistance to Phytophthora sojae in soybean. Nat. Commun. 2021, 12, 6263. [Google Scholar] [CrossRef]
  57. Cook, D.E.; Lee, T.G.; Guo, X.; Melito, S.; Wang, K.; Bayless, A.M.; Wang, J.; Hughes, T.J.; Willis, D.K.; Clemente, T.E.; et al. Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean. Science 2012, 338, 1206–1209. [Google Scholar] [CrossRef]
  58. Sudupak, M.A.; Bennetzen, J.L.; Hulbert, S.H. Unequal exchange and meiotic instability of disease-resistance genes in the Rp1 region of maize. Genetics 1993, 133, 119–125. [Google Scholar] [CrossRef]
  59. Zhang, J.; Sun, S.; Wang, G.; Duan, C.; Wang, X.; Wu, X.; Zhu, Z. Characterization of Phytophthora resistance in soybean cultivars/lines bred in Henan province. Euphytica 2014, 196, 375–384. [Google Scholar] [CrossRef]
  60. Sahoo, D.K.; Das, A.; Huang, X.; Cianzio, S.; Bhattacharyya, M.K. Tightly linked Rps12 and Rps13 genes provide broad-spectrum Phytophthora resistance in soybean. Sci. Rep. 2021, 11, 16907. [Google Scholar] [CrossRef]
  61. Park, K.-Y.; Moon, J.-K.; Yun, H.-T.; Lee, Y.-H.; Kim, S.-L.; Ryu, Y.-H.; Kim, Y.-H.; Ku, J.-H.; Roh, J.-H.; Lee, E.-S.; et al. A new soybean cultivar for fermented soyfood and tofu with high yield, “Daepung”. Korean. J. Breed. Sci. 2005, 37, 111–112. [Google Scholar]
  62. Kim, H.T.; Ko, J.M.; Baek, I.Y.; Jeon, M.K.; Han, W.Y.; Park, K.Y.; Lee, B.W.; Lee, Y.H.; Jung, C.S.; Oh, K.W.; et al. Soybean cultivar for tofu, ‘Saedanbaek’ with disease resistance, and high protein content. Korean. J. Breed. Sci. 2014, 46, 295–301. [Google Scholar] [CrossRef]
  63. Dorrance, A.; Jia, H.; Abney, T. Evaluation of soybean differentials for their interaction with Phytophthora sojae. Plant Health Prog. 2004, 5, 9. [Google Scholar] [CrossRef]
  64. Doyle, J.J.; Doyle, J.L. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 1987, 19, 11–15. [Google Scholar]
  65. Lee, Y.G.; Jeong, N.; Kim, J.H.; Lee, K.; Kim, K.H.; Pirani, A.; Ha, B.K.; Kang, S.T.; Park, B.S.; Moon, J.K.; et al. Development, validation and genetic analysis of a large soybean SNP genotyping array. Plant J. 2015, 81, 625–636. [Google Scholar] [CrossRef]
  66. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  67. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  68. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015, 3, 269–283. [Google Scholar] [CrossRef]
  69. Churchill, G.A.; Doerge, R.W. Empirical threshold values for quantitative trait mapping. Genetics 1994, 138, 963–971. [Google Scholar] [CrossRef] [PubMed]
  70. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phenotypic reactions of Daepung and Saedanbaek following inoculation with Phytophthora sojae isolate 2858.
Figure 1. Phenotypic reactions of Daepung and Saedanbaek following inoculation with Phytophthora sojae isolate 2858.
Plants 12 03957 g001
Figure 2. The genomic region for resistance to Phytophthora sojae isolate 2858 identified on chromosome 3 in the Daepung × Saedanbaek population. The bands on the chromosome indicate that SNP marker loci integrated into the genetic map. The LOD graph to the right of the chromosome show the most likely position of the resistance gene locus. The hatched line on the LOD graph indicates the LOD threshold. The 1- and 2-LOD intervals are displayed as bar and loid lines, respectively.
Figure 2. The genomic region for resistance to Phytophthora sojae isolate 2858 identified on chromosome 3 in the Daepung × Saedanbaek population. The bands on the chromosome indicate that SNP marker loci integrated into the genetic map. The LOD graph to the right of the chromosome show the most likely position of the resistance gene locus. The hatched line on the LOD graph indicates the LOD threshold. The 1- and 2-LOD intervals are displayed as bar and loid lines, respectively.
Plants 12 03957 g002
Figure 3. Comparison of the identified Rps loci between Saedanbaek (RpsSDB) and Daewon (RpsDW). Several nucleotide-binding site-leucine-rich repeat (NBS-LRR)- and serine/threonine protein kinase (STK)-coding genes were annotated in this interval.
Figure 3. Comparison of the identified Rps loci between Saedanbaek (RpsSDB) and Daewon (RpsDW). Several nucleotide-binding site-leucine-rich repeat (NBS-LRR)- and serine/threonine protein kinase (STK)-coding genes were annotated in this interval.
Plants 12 03957 g003
Figure 4. Susceptible reactions of Saedanbaek (left) and Daewon (right) following inoculation with isolates 2457 and 2858 of Phytophthora sojae, respectively.
Figure 4. Susceptible reactions of Saedanbaek (left) and Daewon (right) following inoculation with isolates 2457 and 2858 of Phytophthora sojae, respectively.
Plants 12 03957 g004
Table 1. Goodness-of-fit test for segregation ratios in the 72 recombinant inbred lines (RILs) of the Daepung × Saedanbaek population following inoculation with Phytophthora sojae isolate 2858.
Table 1. Goodness-of-fit test for segregation ratios in the 72 recombinant inbred lines (RILs) of the Daepung × Saedanbaek population following inoculation with Phytophthora sojae isolate 2858.
Parents and RILsObserved 1Expected 1Goodness-of-Fit
RISRISχ2p
Daepung (P1) 50
Saedanbaek (P2)50
Daepung × Saedanbaek (RILs)2544334.92.2534.96.050.05
1 R, resistant reaction; S, susceptible reaction; I, intermediate reaction.
Table 2. Twenty-eight single nucleotide polymorphisms (SNPs) significantly associated with resistance to Phytophthora sojae isolate 2858.
Table 2. Twenty-eight single nucleotide polymorphisms (SNPs) significantly associated with resistance to Phytophthora sojae isolate 2858.
Chr 1Position (bp) 2SNP IDDaepung Allele (S)Saedanbaek Allele (R)AdjustedR2 4
GenotypeFrequencyGenotypeFrequencyp-Value 3
33,373,644AX-90419199AA0.64GG0.331.54 × 10−220.81
33,395,315AX-90404495TT0.64CC0.364.41 × 10−230.82
33,403,744AX-90448600CC0.64TT0.364.41 × 10−230.82
33,403,812AX-90372396AA0.64GG0.364.41 × 10−230.82
33,417,978AX-90354028CC0.64TT0.364.41 × 10−230.82
33,423,161AX-90380038TT0.64CC0.364.41 × 10−230.82
33,458,537AX-90524133CC0.64TT0.364.41 × 10−230.82
33,464,766AX-90414569TT0.64GG0.364.41 × 10−230.82
33,464,955AX-90387750AA0.64GG0.364.41 × 10−230.82
33,484,755AX-90378648TT0.63GG0.364.41 × 10−230.82
33,485,682AX-90470694TT0.64CC0.339.62 × 10−190.76
33,517,886AX-90449575AA0.64CC0.364.41 × 10−230.82
33,554,469AX-90388258TT0.63AA0.363.68 × 10−160.71
33,695,759AX-90317052GG0.63AA0.374.41 × 10−230.82
33,709,588AX-90394026CC0.63TT0.364.41 × 10−230.82
33,759,276AX-90440228CC0.64TT0.364.41 × 10−230.82
33,847,841AX-90375748AA0.63GG0.374.41 × 10−230.82
33,849,443AX-90499181GG0.64AA0.364.41 × 10−230.82
33,855,656AX-90467840CC0.63TT0.374.41 × 10−230.82
33,875,061AX-90449384TT0.64CC0.342.18 × 10−220.81
34,272,521AX-90417885GG0.63AA0.374.41 × 10−230.82
34,272,894AX-90347629TT0.63CC0.374.41 × 10−230.82
34,277,380AX-90310078TT0.63CC0.374.41 × 10−230.82
34,283,885AX-90465452TT0.64AA0.364.41 × 10−230.82
34,284,091AX-90328472CC0.64TT0.303.98 × 10−200.78
34,291,232AX-90331552TT0.63CC0.374.41 × 10−230.82
34,291,566AX-90307867AA0.63TT0.374.41 × 10−230.82
34,295,128AX-90317436AA0.63GG0.374.41 × 10−230.82
1 Chromosome. 2 Physical positions are based on soybean genome Glyma2 (http://soybase.org). 3 Adjusted p-value was obtained using Bonferroni correction (p < 1.91 × 10−6). 4 Phenotypic variance explained by each marker.
Table 3. Information regarding the genetic map of the Daepung × Saedanbaek population.
Table 3. Information regarding the genetic map of the Daepung × Saedanbaek population.
Chr. 1Total Length
(cM) (a)
No. of Total
SNP Markers
No. of Unique
Loci (b) 2
Avg. Marker
Interval (cM) (a/b) 3
187.41485581.5
2100.41335631.6
3127.01165482.6
462.91338461.4
568.61358511.3
679.21583611.3
770.21094461.5
899.21350521.9
954.21198331.6
1074.41515521.4
1173.9498401.8
1249.4987361.4
1392.21704432.1
1466.71092471.4
1595.41936571.7
1659.71281391.5
1792.3709412.3
1899.61822591.7
1966.21594441.5
2067.61086471.4
Total158626,130963
1 Chromosome. 2 Actual number of loci where markers are positioned in the genetic map. 3 Average interval between loci, which was calculated by dividing total length (cM) by no. of unique loci.
Table 4. A genomic region identified for resistance to Phytophthora sojae isolate 2858 in the Daepung × Saedanbaek population.
Table 4. A genomic region identified for resistance to Phytophthora sojae isolate 2858 in the Daepung × Saedanbaek population.
Chr. 1Position (bp) 2Flanking MarkersLOD 3PVE (%) 4
33,373,644…4,295,128AX-90419199…AX-9031743638.793.0
1 Chromosome. 2 Physical position in bp based on the genome version Glyma2. 3 Logarithm of odds (LOD). 4 Phenotypic variance (%) explained by the QTL.
Table 5. Response of Saedanbaek, Daewon, and RILs of two populations following inoculation with Phytophthora sojae isolates 2858 and 2457.
Table 5. Response of Saedanbaek, Daewon, and RILs of two populations following inoculation with Phytophthora sojae isolates 2858 and 2457.
Parents and RILsReaction to P. sojae Isolates Following Inoculation
28582457
DaepungSS
SaedanbaekRS
RILs 1 (N = 25)All RAll S
DaepungSS
DaewonSR
RILs 2 (N = 35)All SAll R
1 The RILs with resistance to the isolate 2858 from the Daepung × Saedanbaek population. 2 The RILs with resistance to the isolate 2457 from the Daepung × Daewon population.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

You, H.J.; Shim, K.-C.; Kang, I.-J.; Kim, J.-M.; Kang, S.; Lee, S. Soybean Variety Saedanbaek Confers a New Resistance Allele to Phytophthora sojae. Plants 2023, 12, 3957. https://doi.org/10.3390/plants12233957

AMA Style

You HJ, Shim K-C, Kang I-J, Kim J-M, Kang S, Lee S. Soybean Variety Saedanbaek Confers a New Resistance Allele to Phytophthora sojae. Plants. 2023; 12(23):3957. https://doi.org/10.3390/plants12233957

Chicago/Turabian Style

You, Hee Jin, Kyu-Chan Shim, In-Jeong Kang, Ji-Min Kim, Sungtaeg Kang, and Sungwoo Lee. 2023. "Soybean Variety Saedanbaek Confers a New Resistance Allele to Phytophthora sojae" Plants 12, no. 23: 3957. https://doi.org/10.3390/plants12233957

APA Style

You, H. J., Shim, K. -C., Kang, I. -J., Kim, J. -M., Kang, S., & Lee, S. (2023). Soybean Variety Saedanbaek Confers a New Resistance Allele to Phytophthora sojae. Plants, 12(23), 3957. https://doi.org/10.3390/plants12233957

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop