Figures
Abstract
Resistance to Phytophthora sojae isolate PsMC1 was evaluated in 102 F2∶3 families derived from a cross between the resistant soybean cultivar Wandou 15 and the susceptible cultivar Williams and genotyped using simple sequence repeat (SSR) markers. The segregation ratio of resistant, segregating, and susceptible phenotypes in the population suggested that the resistance in Wandou 15 was dominant and monogenic. Twenty-six polymorphic SSR markers were identified on soybean chromosome 17 (Molecular linkage group D2; MLG D2), which were linked to the resistance gene based on bulked segregation analysis (BSA). Markers Sattwd15-24/25 and Sattwd15-47 flanked the resistance gene at a distance of 0.5 cM and 0.8 cM, respectively. Two cosegregating markers, Sattwd15-28 and Sattwd15-32, were also screened in this region. This is the first Rps resistance gene mapped on chromosome 17, which is designated as Rps10. Eight putative genes were found in the mapped region between markers Sattwd15-24/25 and Sattwd15-47. Among them, two candidate genes encoding serine/threonine (Ser/Thr) protein kinases in Wandou 15 and Williams were identified and sequenced. And the differences in genomic sequence and the putative amino acid sequence, respectively, were identified within each candidate gene between Wandou 15 and Williams. This novel gene Rps10 and the linked markers should be useful in developing soybean cultivars with durable resistance to P. sojae.
Citation: Zhang J, Xia C, Duan C, Sun S, Wang X, Wu X, et al. (2013) Identification and Candidate Gene Analysis of a Novel Phytophthora Resistance Gene Rps10 in a Chinese Soybean Cultivar. PLoS ONE 8(7): e69799. https://doi.org/10.1371/journal.pone.0069799
Editor: Zhengguang Zhang, Nanjing Agricultural University, China
Received: January 27, 2013; Accepted: June 12, 2013; Published: July 25, 2013
Copyright: © 2013 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was mainly supported by the Special Fund for Agro-scientific Research in the Public Interest (201303018) and the Program of Protection of Crop Germplasm Resources (NB2010-2130135-25-14) from the Ministry of Agriculture of the People's Republic of China. The funders had no role in this study design, data collection and analysis, decision publish, or preparation of the manuscript.
Competing interests: The Institute of Crop Science, Chinese Academy of Agricultural Sciences has applied to the Patent Office, China, for a patent on this gene. The authors confirm that this not alter their adherence to all the PLOS ONE policies on sharing data and materials. And the authors have declared that no competing interests exist.
Introduction
Phytophthora root rot (PRR), caused by the Phytophthora sojae Kaufmann & Gerdemann, is one of the most destructive diseases of soybean (Glycine max (L.) Merr.) [1]. Since it was first reported in the United States in 1948, PRR has been observed in other soybean-producing areas, including Australia, Europe, Asia, and Africa [1], [2]. In China, PRR was first recorded in 1991 in Heilongjiang Province [3], and since then, this disease has occurred in almost all soybean production-areas of the province [4], [5]. P. sojae can infect soybean plants at any stage of growth. It not only kills seedlings and young plants, it also causes wilting and death of adult plants in later stages [6]. PRR reduces soybean yield by 10–40%, and severe infection can result in a complete yield loss during severe epiphytotic outbreaks [7], [8].
The most effective method to control PRR is to plant resistant cultivars, and genes conferring resistance to P. sojae (Rps) have been widely used in commercial soybean cultivars [1], [9]. To date, twenty Rps genes/alleles, involved 14 loci distributed on six different soybean chromosomes, have been reported, including Rps1a, Rps1b, Rps1c, Rps1d, Rps1k, Rps2, Rps3a, Rps3b, Rps3c, Rps4, Rps5, Rps6, Rps7, Rps8, Rps9, RpsYB30, RpsYD25, RpsZS18, RpsSN10, and the Rps gene in Waseshiroge [10]–[16]. Following deployment of each single Rps gene, races of P. sojae were subsequently identified that were virulent to plants carrying the Rps gene. Single Rps genes have been effective for 8 to 15 years, depending on inoculum density and environmental conditions [17]. Chen et al. [18] and Zhu et al. [19] reported that most of the Rps genes (except Rps1c and Rps1k) were not effective against the Chinese populations of P. sojae.
P. sojae interacts with soybean in a gene-for-gene model, and 55 physiological races of P. sojae, and many more pathotypes, have been identified in the United States and other countries [20]–[22]. The populations of P. sojae in China have highly diverse virulence, as in the United States [4], [5], [19]. Zhang et al. [5] found 12 races of P. sojae present in the Heilongjiang Province of China and identified 14 intermediate reaction types on a set of eight soybean differentials. Cui et al. [4] identified 30 virulence pathotypes from the isolates of Heilongjiang Province that were not characterized to any race based on the published race definitions. As the complexity of the virulence pathotypes of P. sojae continues to increase in soybean fields, identification and incorporation of new sources with stable and durable resistance into commercial cultivars using marker-assisted selection (MAS) in soybean breeding programs will be essential to control this disease effectively.
The soybean genome has been sequenced and the reference sequence could be used for developing new markers and identifying candidate genes [23], [24]. Using the soybean genome sequence, Sugimoto et al. [24] localized an Rps gene to a 2.5 cM region of chromosome 3. Saghai Maroof et al. [23] mapped Rsv4 to a physical interval of less than 100 kb on chromosome 2. Sequence analysis revealed that this region contained several predicted transcription factors and unknown protein products, indicating that Rsv4 likely belonged to a new class of resistance genes that interfere with viral infection and cell-to-cell movement, and delay vascular movement, which differed from the Rsv1 and Rsv3 genes, which were characterized as nucleotide binding site-leucine rich repeat (NBS-LRR) type proteins.
Soybean cultivar Wandou 15 was registered in Anhui province in 1996 and is fatty, rich in protein, matures early and is resistant to soybean mosaic virus and downy mildew; therefore, it has been widely planted in Huaihe and the middle-lower Yangtze areas in China [25]. Zhu et al. [26] reported that Wandou 15 was resistant to two new races of P. sojae. Chen et al. [18] postulated that the cultivar carried gene combinations of Rps1c or Rps1k and Rps4. Recently, Xia et al. [27] suggested that Wandou 15 might carry a novel Rps gene.
Here, we report (1) the inheritance of Phytophthora resistance in Wandou 15 by investigating phenotypic data; (2) the fine mapping of the Rps gene with new SSR markers in the genomic regions associated with PRR resistance; and (3) the prediction and cloning of the candidate Rps gene(s) for PRR resistance in Wandou 15.
Results
Phenotype analysis for mapping population
The resistant parent Wandou 15 plants showed no symptoms in response to isolate PsMC1 at 6 days post inoculation (DPI), while the susceptible parent Williams plants showed severe rot at the inoculated location at 6 DPI and all plants ultimately died (Figure 1). Among the 102 F2∶3 families of the mapping population, 31 families were identified as homozygous-resistant (R), 29 families were homozygous-susceptible (S) and 42 families were segregating (Rs) for isolate PsMC1, according to mortality (Table 1). A segregation ratio of 31∶42∶29 in the F2∶3 population fitted well with the genetic model of 1∶2∶1 ratio (χ2 = 3.25, p = 0.20), indicating that Phytophthora resistance in Wandou 15 is controlled by a single dominant locus.
Genetic mapping of the Rps gene in Wandou 15
To map the Rps gene in Wandou 15, the SSR markers described in Soybase (http://soybase.org) was screened using the bulk segregation analysis (BSA). Nine SSR markers (Sat_222, Sat_292, Satt514, Satt461, Satt528, Satt574, Satt543, Satt615, and Satt301) on chromosome 17 were screened and showed polymorphisms between Wandou 15 and Williams, as well as the resistant bulk and the susceptible bulk. The linkage analysis further revealed that the Rps gene in Wandou 15 was linked to these nine SSR markers and was located between Satt543 and Satt615 at a genetic distance of 8.2 cM (Figure 2A). This was a novel Phytophthora resistance locus that differed from the 14 previously reported Rps loci. It was designated as Rps10.
(A) Genetic map of the Phytophthora resistance gene Rps10 was deduced from the segregation analysis of 102 F2∶3 families derived from the cross of Wandou 15 and Williams. The map was generated in Joinmap v. 4.0 using Kosambi's mapping function (Marker names and distances are on the both side of Linkage map). (B) Physical map and gene annotations between Sattwd15-24/25 and Sattwd15-47 on chromosome 17 of Williams 82. Positions of the tightly linked and cosegregating markers are indicated on the right of the chromosome line, and locations of the putative gene models (filled triangles) and genes with no functional annotation for the locus (open triangles) are indicated on the left of the line.
To further fine map Rps10 on chromosome 17, new SSR markers in this region were designed and SSR hunter was used to detect SSRs. BLAST analysis of SSR markers Satt543 and Satt615 detected a 4.29 Mb fragment in the Glycine max v1.0 release of the soybean genome sequence in this region. Four hundred and fourteen out of 671 detected SSR loci with more than 10 bp repeat motifs were selected to design primers to survey polymorphisms between the parental cultivars. Among them, 17 new SSR markers showed polymorphisms between Wandou 15 and Williams, and were used to analyze the F2∶3 population (Table 2).
Based on the genotype data of the 26 polymorphic SSR markers, Rps10 was further localized between markers Sattwd15-24/25 and Sattwd15-47 with a genetic distance of 1.3 cM (Figure 2A). Furthermore, two markers, Sattwd15-28 and Sattwd15-32, cosegregated with Rps10 in this population (Figure 2A).
Candidate gene prediction and analysis
To locate the Rps10 region in the soybean genome sequence, BLASTN searches for the sequences of SSR markers tightly linked to Rps10 were performed initially against the soybean genome sequence, release Glycine max v1.0 (http://www.Phytozome.net). Based on the public soybean physical map [28], the mapped length of the region flanked by markers Sattwd15-24/25 and Sattwd15-47 on chromosome 17 is approximately 311 kb. By inspection of the soybean gene annotation database, eight candidate genes were detected in the mapped region (Figure 2B; Table S1). Four out of the eight genes were annotated with known protein function: one for Prokaryotic DNA topoisomerase, one for Glycosyl transferase family 2, and the other two genes, Glyma17g28950.1 and Glyma17g28970.1, were annotated as serine/threonine (Ser/Thr) protein kinases with plant-type (Figure 2B; Table S1). Furthermore, the two SSR markers Sattwd15-28 and Sattwd15-32 cosegregated with Rps10 were found in Glyma17g28950.1 and Glyma17g28970.1, respectively. The sequences of Wandou 15 and Williams amplified by both markers shared 92.40–96.64% and 96.05–100% identity with Glyma17g28950.1 and Glyma17g28970.1 genes, respectively (Figure S1A, S1B). Thus, Glyma17g28950.1 and Glyma17g28970.1 were also co-segregated with Rps10 phenotype, and most likely the referenced candidate genes of Rps10 (Figure 3A).
(A) The structures of the gene models Glyma17g28950.1, Glyma17g28970.1 from 309×105 to 310×105 bp on chromosome 17 of Williams 82. (B) Predicated structure of the candidate gene Rps10-1. (C) Predicted structure of the putative protein Rps10-1. (D) Predicated structure of the candidate gene Rps10-2. (E) Predicted structure of the putative protein Rps10-2. Filled rectangular arrows with orientations (from 5′ to 3′) indicate the predicated genes. Filled rectangles indicate the exons. Open rectangular arrows with orientations (N-C) indicate the predicated proteins. The filled rectangle indicates the Ser/Thr protein kinase domain.
Candidate gene analysis with allelic sequence comparison
To obtain the fragment of genomic DNA of the two candidate genes in Wandou 15, two pairs of primers were designed according to the regions of the 5′-UTR and 3′-UTR of the gene models Glyma17g28950.1 and Glyma17g28970.1, respectively. These two candidate genes of Rps10 in Wandou 15 were named as Rps10-1 and Rps10-2. The full lengths of Rps10-1 (GenBank No. JX682937) and Rps10-2 (GenBank No. JX682938) were 4294 bp and 4963 bp, respectively (Figure 3B, 3D). Compared with the known genomic sequence, Rps10-1 and Rps10-2 shared 98.72% and 99.62% identity with the sequences of Glyma17g28950.1 and Glyma17g28970.1 in Williams 82. Rps10-1 was predicted to have three exons and two introns (Figure 3B), and encoded a putative 507-aa protein (Figure 3C). Rps10-2 was predicted to have seven exons and six introns (Figure 3D), and encoded a putative 602-aa protein (Figure 3E).
To analyze the susceptible alleles, these two genes were also completely sequenced and investigated for the presence of the characteristic motifs. The allelic sequences rpswm-1 and rpswm-2 in Williams were 4272 bp, 4966 bp in length, and shared 98.35% and 98.23% identity with Rps10-1 and Rps10-2, respectively. The rpswm-1 was predicted to contain three exons and four introns, and was nearly identical in size and sequence with Rps10-1 (Figure S2A). And rpswm-2 was predicted to contain five exons and six introns, which was significantly different from the structure of the Rps10-2 (Figure S2B). Comparison of the predicted protein sequences revealed the polymorphisms in the amino acid sequence present between Rps10-1 and rpswm-1, and between Rps10-2 and rpswm-2 (Figure 4). These results indicated that the sequence variations either between Rps10-1 and rpswm-1, or between Rps10-2 and rpswm-2, may account for the resistance to PRR in Wandou 15, or for susceptibility to PRR in Williams.
(A) The alignment of amino acid sequences of the predicted proteins encoded by Rps10-1 and rpswm-1 genes. (B) The alignment of amino acid sequences of the predicted proteins encoded by Rps10-2 and rpswm-2 genes. The black boxes indicate the conserved domains within Rps10-1 and rpswm-1, Rps10-2 and rpswm-2. Open boxes show the amino acids that differ between Rps10-1 and rpswm-1 and Rps10-2 and rpswm-2.
Discussion
The results of the present study indicated that resistance to P. sojae in Wandou 15 is controlled by a single dominant locus, Rps10, located on soybean chromosome 17. Thus, Rps10 is a novel Rps gene given that it is the first Phytophthora resistance gene identified on chromosome 17. Two markers Sattwd15-28 and Sattwd15-32 were identified that cosegregated with Rps10. Markers Sattwd15-24/25 and Sattwd15-47 were found to flank the gene at a distance of 0.5 and 0.8 cM, respectively. It was theoretically estimated the MAS accuracies would be 100% and 99.2% using both cosegregating markers and flanking markers, respectively.
Several other disease-resistant quantitative trait loci (QTLs) or recessive genes have been mapped at chromosome 17. To the best of our knowledge, three QTLs conferring sudden death syndrome, three QTLs for resistance to soybean cyst nematode, seven QTLs conferring sclerotinia stem rot resistance, rxp gene and several QTLs conferring bacterial leaf pustule resistance are also located at different regions of chromosome 17 [29]–[33]. Rps10 is the first single dominant disease resistance gene identified on chromosome 17 that confers resistance to P. sojae. Kang et al. [34] reported that only nine NBS-LRR genes were located on this chromosome. These data indicate that Rps10 was located in a disease resistance gene-poor region.
To date, 19 Rps genes/alleles have been reported, but only three Rps genes, Rps1k, Rps2 and Rps4, have been cloned and characterized as NBS-LRR genes, which involved recognizing the presence of pathogens and ultimately confer resistance [35]–[37]. In this study, eight putative genes were found in the mapped region and the four un-annotated genes out of the eight genes were analyzed by BLAST searching the nucleotide sequence at the NCBI and the putative protein sequence at the EMBL-EBI databases. Glyma17g28740.1 might function as a DNA/RNA polymerase or in nucleic acid binding. Glyma17g28980.1 shared 100% identity with a gene (GenBank No. XM_003551110) encoding the predicted Glycine max DNA topoisomerase 1-like. Glyma17g29050.1 might encode the 60S ribosomal protein L30-like predicted by the gene (GenBank No. XM_003552594). However, Glyma17g29040.1 only shared 82% identity in nucleotide sequence with Glycine max clone GM_WBc0124O13 (GenBank No. AC236224). There was no similar level of correspondence among the proteins in the EMBL-EBI database. We suggested that Glyma17g28950.1 and Glyma17g28970.1, encoding plant-type Ser/Thr protein kinases, were the candidate genes of Rps10. The two candidate genes, Rps10-1 and Rps10-2, about 14 kb apart, were identified at this region. Using the InterProScan program analysis, their predicted domain structure was: 1) a hydrophobic N-terminal putative signal peptide, and 2) a C-terminal Ser/Thr protein kinase catalytic domain. Additionally, Rps10-1 contained a putative transmembrane domain followed by a membrane transfer-stop signal. These analyses indicated that Rps10-1 and Rps10-2 might be the Ser/Thr protein kinases. The difference in structures between the Ser/Thr genes representing Rps10 and the NBS-LRR genes Rps1k, Rps2, and Rps4 indicated that the mechanism of resistance in the soybean-P. sojae interaction was likely to be different.
Some resistance genes belonging to the RLK families have been characterized as an intracellular Ser/Thr protein kinase, such as Pto, Xa21, Xa26, ZmPto and NgRLK1 [38]–[42]. Based on the structure of the putative extracellular domain, the plant RLKs were divided into seven classes: the S-locus glycoprotein like (SLG-like), Leucine-Rich Repeats (LRRs), epidermal growth factor-like (EGF-like) repeats, putative carbohydrate-binding lectin (Lectin-like), the tumor necrosis factor receptor (TNFR), pathogenesis-related (PR5) and N-glycosylation site type of RLKs [43]. Interestingly, the N-terminal 116 aa of predicted Rps10-1, and 302 aa of predicted Rps10-2 have no similarities to the extracellular regions of the known RLKs. The kinase domain of Rps10-1 and Rps10-2 showed 22–41% sequence similarity to other plant RLKs (data not shown). In addition, Rps10-1 and Rps10-2 showed a low level of sequence identity (4–5%) with the calcium-dependent protein kinase SK5-like proteins in soybean (data not shown). The above analyses suggested that Rps10-1 and Rps10-2 might encode a novel type of plant RLK.
RLKs are unusual membrane-associated plant protein kinases, some of which have important roles in pathogen resistance. For example, the Ser/Thr protein kinase Pto, which confers resistance to Pseudomonas syringae pv. tomato, is autophosphorylated at the Thr38 and Ser198 sites after monitoring the signal of the interaction between effector AvrPto and tomato Prf. Prf is further phosphorylated by another Ser/Thr kinase Pti1 to induce a series of mitogen-activated protein kinase cascades, resulting in a hypersensitive response (HR) [44]–[47]. The rice Xa21 kinase domain interacts with the protein XB3 via its ankyrin repeat domain and trans-phosphorylation occurs. The physical interaction with XB3 is thought to stabilize the Xa21 protein, thereby maintaining R protein levels and hence the capacity to fully activate the defense response.
The differences in the predicted amino acid sequences between Rps10-1 and rpswm-1, Rps10-2 and rpswm-2 might explain the resistance/susceptibility of the parental cultivars. However, the resistance mechanism of Rps10 is still ambiguous. The previous revealed that phosphorylation played an important role in regulating components of signaling pathways involved in the HR. Analysis in the KinasePhos database revealed that seven out 44 Ser-phosphorylation sites, and one out 22 Thr- phosphorylation sites, were detected in Rps10-1.Seven out 43 Ser-phosphorylation sites, and one out 15 Thr-phosphorylation sites, were detected in the Rps10-2 (data not shown). These predicted phosphorylation sites are very important for further investigation of the resistance mechanism of Rps10. The functional significance of the proposed novel Phytophthora resistance mechanism remains to be elucidated.
Materials and Methods
Plant materials
Wandou 15 is a PRR-resistant soybean cultivar, and Williams is a susceptible cultivar. The F1 plants from the cross between Wandou 15 and Williams were self-pollinated to produce an F2 population. Each F2 plant was threshed individually to yield seeds of 102 F2∶3 families. The 102 F2∶3 families were used to investigate both phenotypic and genetic data to map the resistance gene(s) for PRR.
Phytophthora sojae inoculation and PRR evaluation
The P. sojae isolate PsMC1 (pathotype 1a, 1c, 1k, 2, 3b, 3c, 4, 5, 6, 7, 8, ZS18) was used as the inoculum to evaluate the phenotype of parental materials and the F2∶3 families. Fourteen-day-old seedlings of Wandou 15 (20 individuals), Williams (20 individuals), and each F2∶3 family (30 individual seedlings) were inoculated with PsMC1 using the hypocotyl-inoculation technique as described by Haas and Buzzell [48]. The negative control for both parental cultivars was performed in the same way with pure carrot agar slurry instead of PsMC1 inoculum slurry. After inoculation, the plants were placed in a mist room with relative humidity 100% and an average temperature of 25°C for 2 d. Then they were moved to a greenhouse with average temperature 25°C.
The reactions were recorded by the mortality of inoculated plants in each F2∶3 family and parental cultivars 6 DPI. A family with 0–20% seedling death was scored as homozygous-resistant (R), a family with 80–100% seedling death was considered homozygous-susceptible (S), while a family with 21–79% was scored as segregating (Rs) [49].
DNA preparation and pooling for bulk segregation analysis
Equivalent amounts of leaf tissues from 30 seedlings of each family were bulked and placed in liquid nitrogen, and ground into a powder. Genomics DNA was extracted by the Cetyl trimethylammonium bromide method [50]. Resistant and susceptible bulks of DNA for the BSA [51] were prepared by pooling an equal amount of DNA (1 µg) from each of 10 selected homozygous-resistant and 10 susceptible F2∶3 families, respectively. The finial concentration of the DNA bulks was adjusted to 40 ng/µl for PCR.
SSR markers design and screening
The Rps gene in Wandou 15 was first mapped using the SSR markers described in Soybase (http://soybase.org). Subsequently, new SSR markers were designed according to the sequence between flanking markers Satt543 and Satt615, downloaded from Phytozome (http://www.phytozome.net/soybean), in which the SSR was searched by SSR hunter 1.3 (http://en.bio-soft.net/dna/SSRHunter). The primers were designed using Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA) with default parameters, based on the simple repeat sequence.
The PCR reaction was performed in a final volume of 15 µl in a TGRADIENT Thermocycler (Biometra, Goettingen, Germany), and comprised 40 ng genomic DNA, 2.0 µl 10×PCR reaction buffer (2.0 mM MgCl2) (TIANGEN, Beijing, China), 0.2 mM of each dNTP (TIANGEN), 1.0 U of Taq DNA polymerase (TIANGEN) and 0.2 µM of each primer. The PCR amplification consisted of an initial denaturation at 95°C for 5 min; 35 cycles of denaturation at 95°C for 40 s, annealing at temperatures 45–52°C for 40 s and extension at 72°C for 40 s; with a final extension at 72°C for 10 min. PCR products were mixed with 4 µl of 6×loading buffer (0.25% bromophenol blue, 0.25% xylene cyanol FF and 40% sucrose) separated on an 8% polyacrylamide gel.
The polymorphic SSR markers were screened according to the polymorphic DNA fragment between the parents, resistant and susceptible bulks, and further tested in the entire F2∶3 mapping population. SSR markers were scored as AA (homozygous for the Wandou 15 allele), BB (homozygous for the Williams allele) or AB (heterozygous) in the 102 F2∶3 families.
Data analysis and genetic linkage map construction
The segregation patterns of phenotypes and SSR markers selected in the mapping population were tested for the goodness-of-fit to Mendelian segregation ratio using Chi-squared (χ2) analysis. A genetic linkage map of Rps10 was constructed with Joinmap 4.0 linkage analysis software [52].
Candidate gene(s) prediction and sequence analysis
Sequences of candidate gene models within the closely linked SSR markers were obtained from the Phytozome database (http://www.Phytozome.net). And two pairs of primers were designed according to the regions of the 5′-UTR and 3′-UTR of the gene models Glyma17g28950.1 and Glyma17g28970.1, respectively. PCR was performed on genomic DNA using the standard protocol with Taq DNA polymerase (TaKaRa, Dalian, China). The allelic variations of the candidate gene(s) in these two cultivars were compared. Sequences alignment was performed with BioEdit software [53]. Gene prediction was performed using the program GENESCAN (http://genes.mit.edu/GENSCAN.html) and verified by Gene Finder (http://rulai.cshl.org/tools/genefinder/) and GenomeScan (http://genes.mit.edu/genomescan.html). InterProScan (http://www.ebi.ac.uk/InterProScan) was used for protein annotation.
Supporting Information
Figure S1.
Alignment of sequences corresponding to the locus in Wandou 15, Williams amplified using the markers Sattwd15-28 and Sattwd15-32, and the referred region in Glyma17g28950.1 and Glyma17g28970.1. (A) Alignment of the sequences in Wandou 15, Williams amplified using the marker Sattwd15-28, and the referred region in Glyma17g28950.1 (B) Alignment of the sequences in Wandou 15, Williams amplified using the marker Sattwd15-32, and the referred region in Glyma17g28970.1 (EPS).
https://doi.org/10.1371/journal.pone.0069799.s001
(TIF)
Figure S2.
Structures of the gene rpswm-1 and rpswm-2, the allelic of Rps10-1 and Rps10-2 in susceptible cultivar Williams. (A) Predicated structure of the gene rpswm-1. (B) Predicated structure of the gene rpswm-2. Filled rectangular arrows with orientations (from 5′ to 3′) indicate the predicated genes. Filled rectangles indicate the exons (EPS).
https://doi.org/10.1371/journal.pone.0069799.s002
(TIF)
Table S1.
Annotation of genes between markers Sattwd15-24/25 and Sattwd15-47 with the physical location from 30,796,875 to 31,108.028 on chromosome 17 (DOC).
https://doi.org/10.1371/journal.pone.0069799.s003
(DOC)
Acknowledgments
We thank Dr. Schmitthenner at the Ohio State University, USA and Dr. Lijuan Qiu at the Institute of Crop Science, CAAS for providing the cultivars/lines carrying a single Rps at the locus from Rps1 to Rps8 and the susceptible control cultivar Williams.
Author Contributions
Conceived and designed the experiments: ZZ JZ. Performed the experiments: JZ CX. Analyzed the data: CD SS. Contributed reagents/materials/analysis tools: X. Wang X. Wu. Wrote the paper: JZ ZZ.
References
- 1.
Schmitthenner AF (1999) Phytophthora rot of soybean. In: Hartman GL, Sinclair JB and Rupe JC, editors. Compendium of soybean diseases. 4th edn, The American Phytopathological Society Press, St. Paul, Minnesota 39–42
- 2. Kaufmann MJ, Gerdemann JW (1958) Root and stem rot of soybean caused by Phytophthora sojae n. sp. Phytopathology 48: 201–208.
- 3. Shen CY, Su YC (1991) Discovery and preliminary studies of Phytophthora megasperma on soybean in China. Acta Phytopathol Sinica 21: 298.
- 4. Cui KL, Yin WX, Tang QH, Dong SM, Zheng XB, et al. (2010) Distribution, pathotypes, and metalaxyl sensitivity of Phytophthora sojae from Heilongjiang and Fujian provinces in China. Plant Dis 94: 881–884.
- 5. Zhang SZ, Xu PF, Wu JJ, Xue AG, Zhang JX, et al. (2010) Races of Phytophthora sojae and their virulences on soybean cultivars in Heilongjiang, China. Plant Dis 94: 87–91.
- 6. Dorrance AE, McClure SA, St. Martin SK (2003) Effect of partial resistance on Phytophthora stem rot incidence and yield of soybean in Ohio. Plant Dis 87: 308–312.
- 7.
Anderson TR, Tenuta A (2003) Phytophthora rot. pages 155–156 in: diseases of field crops in Canada . Bailey KL, Gossen BD, Gugel RK, Morrall RAA, editors. Canadian Phytopathological Society. University Extension Press, University of Saskatchewan, Saskatoon, Canada
- 8. Li BY, Ma SM (1999) The occurrence of soybean Phytophthora root rot and its control. Chin J Oil Crop Sci 21: 47–50.
- 9. Dorrance AE, Schmitthenner AF (2000) New sources of resistance to Phytophthora sojae in the soybean plant introductions. Plant Dis 84: 1303–1308.
- 10. Fan AY, Wang XM, Fang XP, Wu XF, Zhu ZD (2009) Molecular identification of Phytophthora resistance gene in soybean cultivar Yudou 25. Acta Agron Sinica 35: 1844–1850.
- 11. Sugimoto T, Kato M, Yoshida S, Matsumoto I, Kobayashi T, et al. (2012) Pathogenic diversity of Phytophthora sojae and breeding strategies to develop Phytophthora–resistant soybeans. Breed Sci 61: 511–522.
- 12. Sun S, Wu XL, Zhao JM, Wang YC, Tang QH, et al. (2011) Characterization and mapping of RpsYu25, a novel resistance gene to Phytophthora sojae. Plant Breed 130: 139–143.
- 13. Wu XL, Zhang BQ, Sun S, Zhao JM, Yang F, et al. (2011) Identification, genetic analysis and mapping of resistance to Phytophthora sojae of Pm28 in soybean. Agric Sci in China 10: 1506–1511.
- 14. Yao HY, Wang XM, Wu XF, Xiao YN, Zhu ZD (2010) Molecular mapping of Phytophthora resistance gene in soybean cultivar zaoshu18. J Plant Genet Resour 11: 213–217.
- 15. Yu AL, Xu PF, Wang JS, Zhang SZ, Wu JJ, et al. (2010) Genetic analysis and SSR mapping of gene resistance to Phytophthora sojae race 1 in soybean cv. Suinong 10. Chin J Oil Crop Sci 32: 462–466.
- 16. Zhu ZD, Huo YL, Wang XM, Huang JB, Wu XF (2007) Molecular identification of a novel Phytophthora resistance gene in soybean. Acta Agron Sinica 33: 154–157.
- 17. Schmitthenner AF (1985) Problems and progress in control of Phytophthora root rot of soybean. Plant Dis 69: 362–368.
- 18. Chen XL, Zhu ZD, Wang XM, Xiao YN, Wu XF (2008) Postulation of Phytophthora resistance genes in soybean cultivars or lines. Sci Agric Sinica 41: 1227–1234.
- 19. Zhu ZD, Wang HB, Wang XM, Chang RZ, Wu XF (2003) Distribution and virulence diversity of Phytophthora sojae in China. Sci Agric Sinica 36: 793–799.
- 20.
Grau CR, Dorrance AE, Bond J, Russin J (2004) Fungal diseases. In: Boerma HR, Specht JE (eds) Soybeans: improvement, production and uses. 3rd edn. Agronomy Monogr. American Soc Agron Madison WI 679–763.
- 21. Nelson BD, Mallik I, McEwen D, Christianson T (2008) Pathotypes, distribution, and metalaxyl sensitivity of Phytophthora sojae from North Dakota. Plant Dis 92: 1062–1066.
- 22. Ryley MJ, Obst NR, Irwin JAG, Drenth A (1998) Changes in the racial composition of Phytophthora sojae in Australia between 1979 and 1996. Plant Dis 82: 1048–1054.
- 23. Saghai Maroof MA, Tucker DM, Skoneczka JA, Bowman BC, Tripathy S, et al. (2010) Fine mapping and candidate gene discovery of the soybean mosaic virus resistance gene, Rsv4. The Plant Genome 3: 14–22.
- 24. Sugimoto T, Yoshida S, Kaga A, Hajika M, Watanabe K, et al. (2011) Genetic analysis and Identification of DNA markers linked to a novel Phytophthora sojae resistance gene in the Japanese soybean cultivar Waseshiroge. Euphytica 182: 133–145.
- 25. Li JS (2007) The breeding and cultivation technique of soybean cultivar Wandou 15. Anhui Agri Sci Bull 13: 108.
- 26. Zhu ZD, Wang XM, Wang HB, Wu XF, Zhang L (2001) Identification and race of Phytophthora sojae isolates collected in Mengcheng, Anhui Province. Acta Phytopathol Sinica 1: 236–240.
- 27. Xia CJ, Zhang JQ, Wang XM, Wu XF, Liu ZX, et al. (2011) Analyses of resistance genes to Phytophthora root rot in soybean germplasm. Chinese J Oil Crop Sci 33: 396–401.
- 28. Schmutz J, Cannon SB, Schlueter J, Ma JX, Mitros T, et al. (2010) Genome sequence of the palaeopolyploid soybean. Nature 463: 178–183.
- 29. Concibido VC, Diers BW, Arelli PR (2004) A decade of QTL mapping for cyst nematode resistance in soybean. Crop Sci 44: 1121–1131.
- 30. Farias-Neto AF, Hashmi R, Schmidt ME, Carlson SR, Hart–man GL, et al. (2007) Mapping and conformation of a sudden death syndrome resistance QTL on linkage group D2 from the soybean genotypes ‘PI 567374’ and ‘Ripley’. Mol Breeding 20: 53–62.
- 31. Kim YT, Oh J, Kim KH, Uhm JY, Lee BM (2010) Isolation and characterization of NgRLK1, a receptor–like kinase of Nicotiana glutinosa that interacts with the elicitin of Phytophthora capsici. Mol Biol Rep 37: 717–727.
- 32. Venancio S, Arahana, Graef GL, Specht JE, Steadman JR, et al. (2001) Identification of QTLs for resistance to Sclerotinia sclerotiorum in Soybean. Crop Sci 41: 180–188.
- 33. Yuan JZ, Bashir R, Salas G, Sharma H, Srour A, et al. (2012) New approaches to selecting resistance or tolerance to SDS and fusarium root rot. J Plant Genome Sci 1: 10–17.
- 34. Kang YJ, Kim KH, Shim S, Yoon MY, Sun SL, et al. (2012) Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean. BMC Plant Biol 2012 12: 139 biomedcentral. com/1471-2229/12/139.
- 35. Gao H, Bhattacharyya MK (2008) The soybean–Phytophthora resistance locus Rps1–k encompasses coiled coil–nucleotide binding–leucine rich repeat–like genes and repetitive sequences. BMC Plant Biol 8: 1–14.
- 36. Graham MA, Marek LF, Shoemaker RC (2002) Organization, expression and evolution of a disease resistance gene cluster in soybean. Genetics 162: 1961–1977.
- 37. Sandhu D, Gao H, Cianzio S, Bhattacharyya MK (2004) 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. Genet 168: 2157–2167.
- 38. Kim DH, Kim KH, Van K, Kim MY, Lee SH (2010) Fine mapping of a resistance gene to bacterial leaf pustule. Theor Appl Genet 120: 1443–1450.
- 39. Martin GB, Brommonschenkel SH, Chunwongse J, Frary A, Ganal MW, et al. (1993) Map-based cloning of a protein kinase gene conferring disease resistance in tomato. Science 262: 1432–1436.
- 40. Song WY, Wang GL, Chen LL, Kim HS, Pi LY, et al. (1995) A receptor kinase–like protein encoded by the rice disease resistance gene Xa21. Science 270: 1804–1806.
- 41. Sun X, Cao Y, Yang Z, Xu C, Li X, et al. (2004) Xa26, a gene conferring resistance to Xanthomonas oryzae pv. oryzae in rice, encodes an LRR receptor kinase–like protein. Plant J 37: 517–527.
- 42. Zou HW, Song ZJ, Wu ZY, Zhang XH, Liu HF, et al. (2011) Isolation and analysis of ZmPto from maize, a homologue to Pto. POJ 4: 53–59.
- 43. Satterlee JS, Sussman MR (1998) Unusual membrane–associated protein kinases in higher plants. Membrane Biol 164: 205–213.
- 44. Salmero JM, Oldroyd GED, Rommens CMT, Scofield SR, Kim HS, et al. (1996) Tomato Prf is a member of the leucine-rich repeat class of plant disease resistance genes and lies embedded within the Pto kinase gene cluster. Cell 86: 123–133.
- 45. Sessa GD, Ascenzo M, Martin GB (2000) Thr38 and Ser198 are Pto autophosphorylation sites required for the AvrPto/Pto-mediated hypersensitive response. EMBO J 19: 2257–2269.
- 46. Zhou JM, Loh YT, Bressan RA, Martin GB (1995) The tomato gene Pti1 encodes a serine/threonine kinase that is phosphorylated by Pto and is involved in the hypersensitive response. Cell 83: 925–935.
- 47. Zhou JM, Tang XY, Martin GB (1997) The Pto kinase conferring resistance to tomato bacterial speck disease interacts with proteins that bind a cis-element of pathogenesis-related genes. EMBO J 16: 3207–3218.
- 48. Haas JH, Buzzell RI (1976) New races 5 and 6 of Phytophthora megasperma var. sojae and differential reactions of soybean cultivars for races 1 and 6. Phytopathology 66: 1361–1362.
- 49. Gordon SG, St Martin SK, Dorrance AE (2006) Rps8 maps to a resistance gene rich region on soybean molecular linkage group F. Crop Sci 46: 168–173.
- 50. Allen GC, Flores–Vergara MA, Krasynanski S, Kumar S, Thompson WF (2006) A modified protocol for rapid DNA isolation from plant tissues using cetyl trimethylammonium bromide. Nat Protoc 1: 2320–2325.
- 51. Michelmore RW, Paran I, Kesseli RV (1991) Identification of markers linked to disease–resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations. Proc Natl Acad Sci 88: 9828–9832.
- 52.
Van Ooijen JW (2006) JoinMap 4 Software for the calculation of genetic linkage maps in experimental populations. Kyazma BV, Wageningen, Netherlands.
- 53. Hall TA (1999) BioEdit: A user–friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 41: 95–98.