Molecular mapping of two loci that confer resistance to Asian rust in soybean

Asian soybean rust (ASR) is caused by the fungal pathogen Phakopsora pachyrhizi Sydow & Sydow. It was first identified in Brazil in 2001 and quickly infected soybean areas in several countries in South America. Primary efforts to combat this disease must involve the development of resistant cultivars. Four distinct genes that confer resistance against ASR have been reported: Rpp1, Rpp2, Rpp3, and Rpp4. However, no cultivar carrying any of those resistance loci has been released. The main objective of this study was to genetically map Rpp2 and Rpp4 resistance genes. Two F(2:3) populations, derived from the crosses between the resistant lines PI 230970 (Rpp2), PI 459025 (Rpp4) and the susceptible cultivar BRS 184, were used in this study. The mapping populations and parental lines were inoculated with a field isolate of P. pachyrhizi and evaluated for lesion type as resistant (RB lesions) or susceptible (TAN lesions). The mapping populations were screened with SSR markers, using the bulk segregant analysis (BSA) to expedite the identification of linked markers. Both resistance genes showed an expected segregation ratio for a dominant trait. This study allowed mapping Rpp2 and Rpp4 loci on the linkage groups J and G, respectively. The associated markers will be of great value on marker assisted selection for this trait.     

Distinct biphasic mRNA changes in response to Asian soybean rust infection

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is now established in all major soybean-producing countries. Currently, there is little information about the molecular basis of ASR-soybean interactions, which will be needed to assist future efforts to develop effective resistance. Toward this end, abundance changes of soybean mRNAs were measured over a 7-day ASR infection time course in mock-inoculated and infected leaves of a soybean accession (PI230970) carrying the Rpp2 resistance gene and a susceptible genotype (Embrapa-48). The expression profiles of differentially expressed genes (ASR-infected compared with the mock-inoculated control) revealed a biphasic response to ASR in each genotype. Within the first 12 h after inoculation (hai), which corresponds to fungal germination and penetration of the epidermal cells, differential gene expression changes were evident in both genotypes. mRNA expression of these genes mostly returned to levels found in mock-inoculated plants by 24 hai. In the susceptible genotype, gene expression remained unaffected by rust infection until 96 hai, a time period when rapid fungal growth began. In contrast, gene expression in the resistant genotype diverged from the mock-inoculated control earlier, at 72 h, demonstrating that Rpp2-mediated defenses were initiated prior to this time. These data suggest that ASR initially induces a nonspecific response that is transient or is suppressed when early steps in colonization are completed in both soybean genotypes. The race-specific resistance phenotype of Rpp2 is manifested in massive gene expression changes after the initial response prior to the onset of rapid fungal growth that occurs in the susceptible genotype.     

The Soybean Rhg1 locus for resistance to the soybean cyst nematode Heterodera glycines regulates the expression of a large number of stress- and defense-related genes in degenerating feeding cells

To gain new insights into the mechanism of soybean (Glycine max) resistance to the soybean cyst nematode (Heterodera glycines), we compared gene expression profiles of developing syncytia in soybean near-isogenic lines differing at Rhg1 (for resistance to Heterodera glycines), a major quantitative trait locus for resistance, by coupling laser capture microdissection with microarray analysis. Gene expression profiling revealed that 1,447 genes were differentially expressed between the two lines. Of these, 241 (16.8%) were stress- and defense-related genes. Several stress-related genes were up-regulated in the resistant line, including those encoding homologs of enzymes that lead to increased levels of reactive oxygen species and proteins associated with the unfolded protein response. These results indicate that syncytia induced in the resistant line are undergoing severe oxidative stress and imbalanced endoplasmic reticulum homeostasis, both of which likely contribute to the resistance reaction. Defense-related genes up-regulated within syncytia of the resistant line included those predominantly involved in apoptotic cell death, the plant hypersensitive response, and salicylic acid-mediated defense signaling; many of these genes were either partially suppressed or not induced to the same level by a virulent soybean cyst nematode population for successful nematode reproduction and development on the resistant line. Our study demonstrates that a network of molecular events take place during Rhg1-mediated resistance, leading to a highly complex defense response against a root pathogen.     

Molecular mapping of soybean rust resistance in soybean accession PI 561356 and SNP haplotype analysis of the Rpp1 region in diverse germplasm

Soybean rust (SBR), caused by Phakopsora pachyrhizi Sydow, is one of the most economically important and destructive diseases of soybean [Glycine max (L.) Merr.] and the discovery of novel SBR resistance genes is needed because of virulence diversity in the pathogen. The objectives of this research were to map SBR resistance in plant introduction (PI) 561356 and to identify single nucleotide polymorphism (SNP) haplotypes within the region on soybean chromosome 18 where the SBR resistance gene Rpp1 maps. One-hundred F(2:3) lines derived from a cross between PI 561356 and the susceptible experimental line LD02-4485 were genotyped with genetic markers and phenotyped for resistance to P. pachyrhizi isolate ZM01-1. The segregation ratio of reddish brown versus tan lesion type in the population supported that resistance was controlled by a single dominant gene. The gene was mapped to a 1-cM region on soybean chromosome 18 corresponding to the same interval as Rpp1. A haplotype analysis of diverse germplasm across a 213-kb interval that included Rpp1 revealed 21 distinct haplotypes of which 4 were present among 5 SBR resistance sources that have a resistance gene in the Rpp1 region. Four major North American soybean ancestors belong to the same SNP haplotype as PI 561356 and seven belong to the same haplotype as PI 594538A, the Rpp1-b source. There were no North American soybean ancestors belonging to the SNP haplotypes found in PI 200492, the source of Rpp1, or PI 587886 and PI 587880A, additional sources with SBR resistance mapping to the Rpp1 region.     

Characterization of nonhost resistance of Arabidopsis to the Asian soybean rust

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is a devastating disease of soybean. We report the use of the nonhost plant Arabidopsis thaliana to identify the genetic basis of resistance to P. pachyrhizi. Upon attack by P. pachyrhizi, epidermal cells of wild-type Arabidopsis accumulated H2O2, which likely orchestrates the frequently observed epidermal cell death. However, even when epidermal cell death occurred, fungal hyphae grew on and infection was terminated at the mesophyll boundary. These events were associated with expression of PDF1.2, suggesting that P. pachyrhizi, an ostensible biotroph, mimics aspects of a necrotroph. Extensive colonization of the mesophyll occurred in Arabidopsis pen mutants with defective penetration resistance. Although haustoria were found occasionally in mesophyll cells, the successful establishment of biotrophy failed, as evidenced by the cessation of fungal growth. Double mutants affected in either jasmonic acid or salicylic acid signaling in the pen3-1 background revealed the involvement of both pathways in nonhost resistance (NHR) of Arabidopsis to P. pachyrhizi. Interestingly, expression of AtNHL10, a gene that is expressed in tissue undergoing the hypersensitive response, was only triggered in infected pen3-1 mutants. Thus, a suppression of P. pachyrhizi-derived effectors by PEN3 can be inferred. Our results demonstrate that Arabidopsis can be used to study mechanisms of NHR to ASR.     

Constituents of soybean cultivars differing in insect resistance

Leaf samples of two insect-resistant plant introductions (PI 227687 and PI 229358) and two susceptible (Ransom and Coker Hampton 266A) cultivars of soybeans, Glycine max, were analyzed at different growth stages for their contents of total nitrogen, carbohydrates, organic acids, and sterols. The two susceptible cultivars accumulated more total nitrogen and at a faster rate than did the two resistant plant introductions. At pod-filling, the two resistant cultivars had equivalent soluble carbohydrates and 33% more than the susceptible cultivars. The quantity of organic acids was essentially the same for the two susceptible cultivars. The resistant cultivars had distinctly different quantities from each other as well as from the susceptible cultivars. The quantity of total sterol of these soybean cultivars varied during the growth of the plant. The resistant cultivars accumulated sterol faster and by pod-filling contained from 20 to 50% more sterol than did the susceptible cultivars.     

Biphasic gene expression changes elicited by Phakopsora pachyrhizi in soybean correlate with fungal penetration and haustoria formation

Inoculation of soybean (Glycine max) plants with Phakopsora pachyrhizi, the causal organism of Asian soybean rust, elicits a biphasic response characterized by a burst of differential gene expression in the first 12 h. A quiescent period occurs from 24 to 48 h after inoculation, in which P. pachyrhizi continues to develop but does not elicit strong host responses, followed by a second phase of intense gene expression. To correlate soybean responses with P. pachyrhizi growth and development, we inoculated the soybean cultivar Ankur (accession PI462312), which carries the Rpp3 resistance gene, with avirulent and virulent isolates of P. pachyrhizi. The avirulent isolate Hawaii 94-1 elicits hypersensitive cell death that limits fungal growth on Ankur and results in an incompatible response, while the virulent isolate Taiwan 80-2 grows extensively, sporulates profusely, and produces a compatible reaction. Inoculated leaves were collected over a 288-h time course for microarray analysis of soybean gene expression and microscopic analysis of P. pachyrhizi growth and development. The first burst in gene expression correlated with appressorium formation and penetration of epidermal cells, while the second burst of gene expression changes followed the onset of haustoria formation in both compatible and incompatible interactions. The proliferation of haustoria coincided with the inhibition of P. pachyrhizi growth in the incompatible interaction or the beginning of accelerated growth in the compatible interaction. The temporal relationships between P. pachyrhizi growth and host responses provide an important context in which to view interacting gene networks that mediate the outcomes of their interactions.     

Development of simple sequence repeat markers for the soybean rust fungus, Phakopsora pachyrhizi

Twenty-four simple sequence repeat markers were developed for Phakopsora pachyrhizi, a fungal pathogen of soybean (Glycine max) and other legumes. All 24 of the loci were evaluated on 28 isolates of P. pachyrhizi. Twenty-one loci were polymorphic, with allelic diversity ranging from two to eight alleles, and null alleles were observed for eight of the 24 loci. A preliminary screen with the closely related species, P. meibomiae, indicated that these primer pairs are specific to P. pachyrhizi.     

A microarray analysis for differential gene expression in the soybean genome using Bioconductor and R

This article describes specific procedures for conducting quality assessment of Affymetrix GeneChip(R) soybean genome data and for performing analyses to determine differential gene expression using the open-source R programming environment in conjunction with the open-source Bioconductor software. We describe procedures for extracting those Affymetrix probe set IDs related specifically to the soybean genome on the Affymetrix soybean chip and demonstrate the use of exploratory plots including images of raw probe-level data, boxplots, density plots and M versus A plots. RNA degradation and recommended procedures from Affymetrix for quality control are discussed. An appropriate probe-level model provides an excellent quality assessment tool. To demonstrate this, we discuss and display chip pseudo-images of weights, residuals and signed residuals and additional probe-level modeling plots that may be used to identify aberrant chips. The Robust Multichip Averaging (RMA) procedure was used for background correction, normalization and summarization of the AffyBatch probe-level data to obtain expression level data and to discover differentially expressed genes. Examples of boxplots and MA plots are presented for the expression level data. Volcano plots and heatmaps are used to demonstrate the use of (log) fold changes in conjunction with ordinary and moderated t-statistics for determining interesting genes. We show, with real data, how implementation of functions in R and Bioconductor successfully identified differentially expressed genes that may play a role in soybean resistance to a fungal pathogen, Phakopsora pachyrhizi. Complete source code for performing all quality assessment and statistical procedures may be downloaded from our web source: http://css.ncifcrf.gov/services/download/MicroarraySoybean.zip.     

Determining the optimal timing of foliar insecticide applications for control of soybean aphid (Hemiptera: Aphididae) on soybean

Field experiments were performed over 3 yr to examine the impact of insecticide application timing to control soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), populations and to prevent soybean yield losses. Experiments were conducted in early and late-planted soybean, Glycine max (L.) Merr. Insecticide applications were made based on soybean growth stages. In 2001, applications were made at V1, V3, R2, and R3 growth stages; in 2002 and 2003, applications were made at R2, R3, and R4 stages. Additional treatments consisted of an unsprayed control and a multiple spray treatment that received insecticide applications at 7-10-d intervals. Soybean aphid densities were recorded throughout the growing season, and yields were measured. Soybean aphid populations varied considerably across years and planting dates. In general, late-planted soybean exhibited higher aphid pressure than early planted soybean, and experiments in 2002 had lower aphid numbers than those in 2001 and 2003. The multiple spray treatment significantly increased yield over the control in four of the six experiments, the exceptions being 2002 late planted and 2003 early planted. This suggests that soybean aphid populations were not large enough to cause yield losses in these two experiments. The R3 spray treatment increased yield in three of the six experiments (2001 late planting, 2002 early planting, and 2003 late planting), the R2 spray treatment increased yield in two of six experiments (2001 and 2003 late plantings), and the V1 application increased yield over the control in the 2001 late-planted experiment. Results suggest that when aphid populations are high insecticide applications made at R2 and R3 plant stages are most effective in preventing yield loss.