Analysis of expressed sequence tags from roots of resistant soybean infected by the soybean cyst nematode

Glycine max polygalacturonase inhibiting protein 11 (GmPGIP11) functions in the root to suppress Heterodera glycines parasitism

Pathogen-secreted polygalacturonases (PGs) alter plant cell wall structure by cleaving the α-(1 → 4) linkages between D-galacturonic acid residues in homogalacturonan (HG), macerating the cell wall, facilitating infection. Plant PG inhibiting proteins (PGIPs) disengage pathogen PGs, impairing infection. The soybean cyst nematode, Heterodera glycines, obligate root parasite produces secretions, generating a multinucleate nurse cell called a syncytium, a byproduct of the merged cytoplasm of 200-250 root cells, occurring through cell wall maceration. The common cytoplasmic pool, surrounded by an intact plasma membrane, provides a source from which H. glycines derives nourishment but without killing the parasitized cell during a susceptible reaction. The syncytium is also the site of a naturally-occurring defense response that happens in specific G. max genotypes. Transcriptomic analyses of RNA isolated from the syncytium undergoing the process of defense have identified that one of the 11 G. max PGIPs, GmPGIP11, is expressed during defense. Functional transgenic analyses show roots undergoing GmPGIP11 overexpression (OE) experience an increase in its relative transcript abundance (RTA) as compared to the ribosomal protein 21 (GmRPS21) control, leading to a decrease in H. glycines parasitism as compared to the overexpression control. The GmPGIP11 undergoing RNAi experiences a decrease in its RTA as compared to the GmRPS21 control with transgenic roots experiencing an increase in H. glycines parasitism as compared to the RNAi control. Pathogen associated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) components are shown to influence GmPGIP11 expression while numerous agricultural crops are shown to have homologs.

Identification of Candidate Genes Controlling Soybean Cyst Nematode Resistance in "Handou 10" Based on Genome and Transcriptome Analyzes

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is a highly destructive pathogen for soybean production worldwide. The use of resistant varieties is the most effective way of preventing yield loss. Handou 10 is a commercial soybean variety with desirable agronomic traits and SCN resistance, however genes underlying the SCN resistance in the variety are unknown. An F_2:8 recombinant inbred line (RIL) population derived from a cross between Zheng 9525 (susceptible) and Handou 10 was developed and its resistance to SCN HG type 2.5.7 (race 1) and 1.2.5.7 (race 2) was identified. We identified seven quantitative trait loci (QTLs) with additive effects. Among these, three QTLs on Chromosomes 7, 8, and 18 were resistant to both races. These QTLs could explain 1.91-7.73% of the phenotypic variation of SCN's female index. The QTLs on chromosomes 8 and 18 have already been reported and were most likely overlapped with rhg1 and Rhg4 loci, respectively. However, the QTL on chromosome 7 was novel. Candidate genes for the three QTLs were predicted through genes functional analysis and transcriptome analysis of infected roots of Handou 10 vs. Zheng 9525. Transcriptome analysis performed also indicated that the plant-pathogen interaction played an important role in the SCN resistance for Handou 10. The information will facilitate SCN-resistant gene cloning, and the novel resistant gene will be a source for improving soybeans' resistance to SCN.

Early transcriptional responses to soybean cyst nematode HG Type 0 show genetic differences among resistant and susceptible soybeans

KEY MESSAGE

Root transcriptome profiling of three soybean cultivars and a wild relative infected with soybean cyst nematode at migratory phase revealed differential resistance pathway responses between resistant and susceptible genotypes. The soybean cyst nematode (SCN), Heterodera glycines, is the most serious pathogen of soybean production throughout the world. Using resistant cultivars is the primary management strategy against SCN infestation. To gain insight into the still obscure mechanisms of genetic resistance to nematodes in different soybean genotypes, RNA-Seq profiling of the roots of Glycine max cv. Peking, Fayette, Williams 82, and a wild relative (Glycine soja PI 468916) was performed during SCN infection at the migratory phase. The analysis showed statistically significant changes of expression beginning at eight hours after inoculation in genes associated with defense mechanisms and pathways, such as the phenylpropanoid biosynthesis pathway, plant innate immunity and hormone signaling. Our results indicate the importance of the early plant response to migratory phase nematodes in pathogenicity determination. The transcriptome changes occurring during early SCN infection included a number of genes and pathways specific to the different resistant genotypes. We observed the most extensive resistant transcriptome reaction in PI 468916, where the resistant response was qualitatively different from that of commonly used G. max varieties.

Genomics of Plant Disease Resistance in Legumes

The constant interactions between plants and pathogens in the environment and the resulting outcomes are of significant importance for agriculture and agricultural scientists. Disease resistance genes in plant cultivars can break down in the field due to the evolution of pathogens under high selection pressure. Thus, the protection of crop plants against pathogens is a continuous arms race. Like any other type of crop plant, legumes are susceptible to many pathogens. The dawn of the genomic era, in which high-throughput and cost-effective genomic tools have become available, has revolutionized our understanding of the complex interactions between legumes and pathogens. Genomic tools have enabled a global view of transcriptome changes during these interactions, from which several key players in both the resistant and susceptible interactions have been identified. This review summarizes some of the large-scale genomic studies that have clarified the host transcriptional changes during interactions between legumes and their plant pathogens while highlighting some of the molecular breeding tools that are available to introgress the traits into breeding programs. These studies provide valuable insights into the molecular basis of different levels of host defenses in resistant and susceptible interactions.

Regulatory interplay between soybean root and soybean cyst nematode during a resistant and susceptible reaction

BACKGROUND

Plant-parasitic nematodes (PPNs) are obligate parasites that feed on the roots of living host plants. Often, these nematodes can lay hundreds of eggs, each capable of surviving without a host for as long as 12 years. When it comes to wreaking havoc on agricultural yield, few nematodes can compare to the soybean cyst nematode (SCN). Quantifying soybean (Glycine max) transcription factor binding sites (TFBSs) during a late-stage SCN resistant and susceptible reaction can shed light onto the systematic interplay between host and pathogen, thereby elucidating underlying cis-regulatory mechanisms.

RESULTS

We sequenced the soybean root transcriptome at 6 and 8 days upon independent inoculation with a virulent and avirulent SCN population. Genes such as β-1,4 glucanase, chalcone synthase, superoxide dismutase and various heat shock proteins (HSPs) exhibited reaction-specific expression profiles. Several likely defense-response genes candidates were also identified which are believed to confer SCN resistance. To explore magnitude of TFBS representation during SCN pathogenesis, a multivariate statistical software identified 46 over-represented TFBSs which capture soybean regulatory dynamics across both reactions.

CONCLUSIONS

Our results reveal a set of soybean TFBSs which are over-represented solely throughout a resistant and susceptible SCN reaction. This set furthers our understanding of soybean cis-regulatory dynamics by providing reaction-specific levels of over-representation at 6 and 8 days after inoculation (dai) with SCN.

Engineered resistance and hypersusceptibility through functional metabolic studies of 100 genes in soybean to its major pathogen, the soybean cyst nematode

During pathogen attack, the host plant induces genes to ward off the pathogen while the pathogen often produces effector proteins to increase susceptibility of the host. Gene expression studies of syncytia formed in soybean root by soybean cyst nematode (Heterodera glycines) identified many genes altered in expression in resistant and susceptible roots. However, it is difficult to assess the role and impact of these genes on resistance using gene expression patterns alone. We selected 100 soybean genes from published microarray studies and individually overexpressed them in soybean roots to determine their impact on cyst nematode development. Nine genes reduced the number of mature females by more than 50 % when overexpressed, including genes encoding ascorbate peroxidase, β-1,4-endoglucanase, short chain dehydrogenase, lipase, DREPP membrane protein, calmodulin, and three proteins of unknown function. One gene encoding a serine hydroxymethyltransferase decreased the number of mature cyst nematode females by 45 % and is located at the Rhg4 locus. Four genes increased the number of mature cyst nematode females by more than 200 %, while thirteen others increased the number of mature cyst nematode females by more than 150 %. Our data support a role for auxin and ethylene in susceptibility of soybean to cyst nematodes. These studies highlight the contrasting gene sets induced by host and nematode during infection and provide new insights into the interactions between host and pathogen at the molecular level. Overexpression of some of these genes result in a greater decrease in the number of cysts formed than recognized soybean cyst nematode resistance loci.

Analysis of response mechanism in soybean under low oxygen and flooding stresses using gel-base proteomics technique

A proteomics approach was used to analyze the response mechanism in soybean seedlings under low oxygen and flooding stresses. Three-day-old soybean seedlings were subjected to low oxygen and flooding stresses. Growth of root was suppressed in both stresses with more extent of suppression in flooded seedlings at 3 and 6 days following the treatments. Proteins were extracted from roots and separated by two-dimensional polyacrylamide gel electrophoresis. Of total 1,233 protein spots, 27 protein spots were commonly changed under low oxygen and flooding stresses; while the differential change in 4 and 18 protein spots was specific to low oxygen and flooding stresses, respectively. Proteins related to metabolism and energy were increased; while protein destination/storage related proteins were decreased commonly under low oxygen and flooding stresses. Protein specie, TCP domain class transcription factor was decreased specifically under low oxygen stress; while decrease of nine proteins related to metabolism, protein destination/storage and disease/defense was specific in flooded seedlings. The decrease in majority of the proteins related to protein destination/storage specifically in flooded seedlings implies the misfolding of proteins resulting in flooded injuries in an independent way of oxygen deprivation. These results suggest that decrease in proteins related to protein destination/storage and disease/defense causes more growth suppression in soybean seedlings under flooding stress compared to low oxygen stress.

Transcriptome analysis of resistant soybean roots infected by Meloidogyne javanica

Soybean is an important crop for Brazilian agribusiness. However, many factors can limit its production, especially root-knot nematode infection. Studies on the mechanisms employed by the resistant soybean genotypes to prevent infection by these nematodes are of great interest for breeders. For these reasons, the aim of this work is to characterize the transcriptome of soybean line PI 595099-Meloidogyne javanica interaction through expression analysis. Two cDNA libraries were obtained using a pool of RNA from PI 595099 uninfected and M. javanica (J(2)) infected roots, collected at 6, 12, 24, 48, 96, 144 and 192 h after inoculation. Around 800 ESTs (Expressed Sequence Tags) were sequenced and clustered into 195 clusters. In silico subtraction analysis identified eleven differentially expressed genes encoding putative proteins sharing amino acid sequence similarities by using BlastX: metallothionein, SLAH4 (SLAC1 Homologue 4), SLAH1 (SLAC1 Homologue 1), zinc-finger proteins, AN1-type proteins, auxin-repressed proteins, thioredoxin and nuclear transport factor 2 (NTF-2). Other genes were also found exclusively in nematode stressed soybean roots, such as NAC domain-containing proteins, MADS-box proteins, SOC1 (suppressor of overexpression of constans 1) proteins, thioredoxin-like protein 4-Coumarate-CoA ligase and the transcription factor (TF) MYBZ2. Among the genes identified in non-stressed roots only were Ser/Thr protein kinases, wound-induced basic protein, ethylene-responsive family protein, metallothionein-like protein cysteine proteinase inhibitor (cystatin) and Putative Kunitz trypsin protease inhibitor. An understanding of the roles of these differentially expressed genes will provide insights into the resistance mechanisms and candidate genes involved in soybean-M. javanica interaction and contribute to more effective control of this pathogen.

Analysis of ascorbate peroxidase genes expressed in resistant and susceptible wheat lines infected by the cereal cyst nematode, Heterodera avenae

Changes in ascorbate peroxidase (APX) enzyme activity in response to nematode (Heterodera avenae) attack were studied in roots of three hexaploid wheat lines carrying Cre2, Cre5, or Cre7 nematode resistance genes and the susceptible Triticum aestivum cv. Anza. A spectrophotometric analysis was carried out with root extracts of infected plants 4, 7, 11, and 14 days after nematode inoculation using uninfected plant as control. APX induction in infected resistant genotypes was similar and higher than in the susceptible control. The introgression wheat/Aegilops ventricosa H-93-8 line, carrying the Cre2 gene, and its parental line H-10-15 as susceptible control were used to analyze whether this increase of activity was correlated with the induction of APX gene expression. Genes encoding cytosolic forms of APX were induced in roots of both lines in response to nematode infection. This induction took place both earlier and with greater intensity in the resistant line than in the susceptible one, and it was also higher in the root area at the site of nematode attachment.

Microarray Detection Call Methodology as a Means to Identify and Compare Transcripts Expressed within Syncytial Cells from Soybean (Glycine max) Roots Undergoing Resistant and Susceptible Reactions to the Soybean Cyst Nematode (Heterodera glycines)

Background. A comparative microarray investigation was done using detection call methodology (DCM) and differential expression analyses. The goal was to identify genes found in specific cell populations that were eliminated by differential expression analysis due to the nature of differential expression methods. Laser capture microdissection (LCM) was used to isolate nearly homogeneous populations of plant root cells. Results. The analyses identified the presence of 13,291 transcripts between the 4 different sample types. The transcripts filtered down into a total of 6,267 that were detected as being present in one or more sample types. A comparative analysis of DCM and differential expression methods showed a group of genes that were not differentially expressed, but were expressed at detectable amounts within specific cell types. Conclusion. The DCM has identified patterns of gene expression not shown by differential expression analyses. DCM has identified genes that are possibly cell-type specific and/or involved in important aspects of plant nematode interactions during the resistance response, revealing the uniqueness of a particular cell population at a particular point during its differentiation process.

The first set of EST resource for gene discovery and marker development in pigeonpea (Cajanus cajan L.)

BACKGROUND

Pigeonpea (Cajanus cajan (L.) Millsp) is one of the major grain legume crops of the tropics and subtropics, but biotic stresses [Fusarium wilt (FW), sterility mosaic disease (SMD), etc.] are serious challenges for sustainable crop production. Modern genomic tools such as molecular markers and candidate genes associated with resistance to these stresses offer the possibility of facilitating pigeonpea breeding for improving biotic stress resistance. Availability of limited genomic resources, however, is a serious bottleneck to undertake molecular breeding in pigeonpea to develop superior genotypes with enhanced resistance to above mentioned biotic stresses. With an objective of enhancing genomic resources in pigeonpea, this study reports generation and analysis of comprehensive resource of FW- and SMD- responsive expressed sequence tags (ESTs).

RESULTS

A total of 16 cDNA libraries were constructed from four pigeonpea genotypes that are resistant and susceptible to FW ('ICPL 20102' and 'ICP 2376') and SMD ('ICP 7035' and 'TTB 7') and a total of 9,888 (9,468 high quality) ESTs were generated and deposited in dbEST of GenBank under accession numbers GR463974 to GR473857 and GR958228 to GR958231. Clustering and assembly analyses of these ESTs resulted into 4,557 unique sequences (unigenes) including 697 contigs and 3,860 singletons. BLASTN analysis of 4,557 unigenes showed a significant identity with ESTs of different legumes (23.2-60.3%), rice (28.3%), Arabidopsis (33.7%) and poplar (35.4%). As expected, pigeonpea ESTs are more closely related to soybean (60.3%) and cowpea ESTs (43.6%) than other plant ESTs. Similarly, BLASTX similarity results showed that only 1,603 (35.1%) out of 4,557 total unigenes correspond to known proteins in the UniProt database (or= 5 sequences detected 102 single nucleotide polymorphisms (SNPs) in 37 contigs. As an example, a set of 10 contigs were used for confirming in silico predicted SNPs in a set of four genotypes using wet lab experiments. Occurrence of SNPs were confirmed for all the 6 contigs for which scorable and sequenceable amplicons were generated. PCR amplicons were not obtained in case of 4 contigs. Recognition sites for restriction enzymes were identified for 102 SNPs in 37 contigs that indicates possibility of assaying SNPs in 37 genes using cleaved amplified polymorphic sequences (CAPS) assay.

CONCLUSION

The pigeonpea EST dataset generated here provides a transcriptomic resource for gene discovery and development of functional markers associated with biotic stress resistance. Sequence analyses of this dataset have showed conservation of a considerable number of pigeonpea transcripts across legume and model plant species analysed as well as some putative pigeonpea specific genes. Validation of identified biotic stress responsive genes should provide candidate genes for allele mining as well as candidate markers for molecular breeding.

Syncytium gene expression in Glycine max([PI 88788]) roots undergoing a resistant reaction to the parasitic nematode Heterodera glycines

The plant parasitic nematode, Heterodera glycines is the major pathogen of Glycine max (soybean). H. glycines accomplish parasitism by creating a nurse cell known as the syncytium from which it feeds. The syncytium undergoes two developmental phases. The first is a parasitism phase where feeding sites are selected, initiating the development of the syncytium. During this earlier phase (1-4 days post infection), syncytia undergoing resistant and susceptible reactions appear the same. The second phase is when the resistance response becomes evident (between 4 and 6dpi) and is completed by 9dpi. Analysis of the resistant reaction of G. max genotype PI 88788 (G. max([PI 88788])) to H. glycines population NL1-RHg/HG-type 7 (H. glycines([NL1-RHg/HG-type 7])) is accomplished by laser microdissection of syncytia at 3, 6 and 9dpi. Comparative analyses are made to pericycle and their neighboring cells isolated from mock-inoculated roots. These analyses reveal induced levels of the jasmonic acid biosynthesis and 13-lipoxygenase pathways. Direct comparative analyses were also made of syncytia at 6 days post infection to those at 3dpi (base line). The comparative analyses were done to identify localized gene expression that characterizes the resistance phase of the resistant reaction. The most highly induced pathways include components of jasmonic acid biosynthesis, 13-lipoxygenase pathway, S-adenosyl methionine pathway, phenylpropanoid biosynthesis, suberin biosynthesis, adenosylmethionine biosynthesis, ethylene biosynthesis from methionine, flavonoid biosynthesis and the methionine salvage pathway. In comparative analyses of 9dpi to 6dpi (base line), these pathways, along with coumarin biosynthesis, cellulose biosynthesis and homogalacturonan degradation are induced. The experiments presented here strongly implicate the jasmonic acid defense pathway as a factor involved in the localized resistant reaction of G. max([PI 88788]) to H. glycines([NL1-RHg/HG-type 7]).

A gene expression analysis of syncytia laser microdissected from the roots of the Glycine max (soybean) genotype PI 548402 (Peking) undergoing a resistant reaction after infection by Heterodera glycines (soybean cyst nematode)

The syncytium is a nurse cell formed within the roots of Glycine max by the plant parasitic nematode Heterodera glycines. Its development and maintenance are essential for nematode survival. The syncytium appears to undergo two developmental phases during its maturation into a functional nurse cell. The first phase is a parasitism phase where the nematode establishes the molecular circuitry that during the second phase ensures a compatible interaction with the plant cell. The cytological features of syncytia undergoing susceptible or resistant reactions appear the same during the parasitism phase. Depending on the outcome of any defense response, the second phase is a period of syncytium maintenance (susceptible reaction) or failure (resistant reaction). In the analyses presented here, the localized gene expression occurring at the syncytium during the resistant reaction was studied. This was accomplished by isolating syncytial cells from Glycine max genotype Peking (PI 548402) by laser capture microdissection. Microarray analyses using the Affymetrix soybean GeneChip directly compared Peking syncytia undergoing a resistant reaction to those undergoing a susceptible reaction during the parasitism phase of the resistant reaction. Those analyses revealed lipoxygenase-9 and lipoxygenase-4 as the most highly induced genes in the resistant reaction. The analysis also identified induced levels of components of the phenylpropanoid pathway. These genes included phenylalanine ammonia lyase, chalcone isomerase, isoflavone reductase, cinnamoyl-CoA reductase and caffeic acid O-methyltransferase. The presence of induced levels of these genes implies the importance of jasmonic acid and phenylpropanoid signaling pathways locally at the site of the syncytium during the resistance phase of the resistant reaction. The analysis also identified highly induced levels of four S-adenosylmethionine synthetase genes, the EARLY-RESPONSIVE TO DEHYDRATION 2 gene and the 14-3-3 gene known as GENERAL REGULATORY FACTOR 2. Subsequent analyses studied microdissected syncytial cells at 3, 6 and 9 days post infection (dpi) during the course of the resistant reaction, resulting in the identification of signature gene expression profiles at each time point in a single G. max genotype, Peking.

The nematode resistance allele at the rhg1 locus alters the proteome and primary metabolism of soybean roots

Heterodera glycines, the soybean cyst nematode (SCN), causes the most damaging chronic disease of soybean (Glycine max). Host resistance requires the resistance allele at rhg1. Resistance destroys the giant cells created in the plant's roots by the nematodes about 24 to 48 h after commencement of feeding. In addition, 4 to 8 d later, a systemic acquired resistance develops that discourages later infestations. The molecular mechanisms that control the rhg1-mediated resistance response appear to be multigenic and complex, as judged by transcript abundance changes, even in near isogenic lines (NILs). This study aimed to focus on key posttranscriptional changes by identifying proteins and metabolites that were increased in abundance in both resistant and susceptible NILs. Comparisons were made among NILs 10 d after SCN infestation and without SCN infestation. Two-dimensional gel electrophoresis resolved more than 1,000 protein spots on each gel. Only 30 protein spots with a significant (P < 0.05) difference in abundance of 1.5-fold or more were found among the four treatments. The proteins in these spots were picked, trypsin digested, and analyzed using quadrupole time-of-flight tandem mass spectrometry. Protein identifications could be made for 24 of the 30 spots. Four spots contained two proteins, so that 28 distinct proteins were identified. The proteins were grouped into six functional categories. Metabolite analysis by gas chromatography-mass spectrometry identified 131 metabolites, among which 58 were altered by one or more treatment; 28 were involved in primary metabolism. Taken together, the data showed that 17 pathways were altered by the rhg1 alleles. Pathways altered were associated with systemic acquired resistance-like responses, including xenobiotic, phytoalexin, ascorbate, and inositol metabolism, as well as primary metabolisms like amino acid synthesis and glycolysis. The pathways impacted by the rhg1 allelic state and SCN infestation agreed with transcript abundance analyses but identified a smaller set of key proteins. Six of the proteins lay within the same small region of the interactome identifying a key set of 159 interacting proteins involved in transcriptional control, nuclear localization, and protein degradation. Finally, two proteins (glucose-6-phosphate isomerase [EC 5.3.1.9] and isoflavone reductase [EC 1.3.1.45]) and two metabolites (maltose and an unknown) differed in resistant and susceptible NILs without SCN infestation and may form the basis of a new assay for the selection of resistance to SCN in soybean.

Population-specific gene expression in the plant pathogenic nematode Heterodera glycines exists prior to infection and during the onset of a resistant or susceptible reaction in the roots of the Glycine max genotype Peking

BACKGROUND

A single Glycine max (soybean) genotype (Peking) reacts differently to two different populations of Heterodera glycines (soybean cyst nematode) within the first twelve hours of infection during resistant (R) and susceptible (S) reactions. This suggested that H. glycines has population-specific gene expression signatures. A microarray analysis of 7539 probe sets representing 7431 transcripts on the Affymetrix soybean GeneChip were used to identify population-specific gene expression signatures in pre-infective second stage larva (pi-L2) prior to their infection of Peking. Other analyses focused on the infective L2 at 12 hours post infection (i-L2(12h)), and the infective sedentary stages at 3 days post infection (i-L2(3d)) and 8 days post infection (i-L2/L3(8d)).

RESULTS

Differential expression and false discovery rate (FDR) analyses comparing populations of pi-L2 (i.e., incompatible population, NL1-RHg to compatible population, TN8) identified 71 genes that were induced in NL1-RHg as compared to TN8. These genes included putative gland protein G23G12, putative esophageal gland protein Hgg-20 and arginine kinase. The comparative analysis of pi-L2 identified 44 genes that were suppressed in NL1-RHg as compared to TN8. These genes included a different Hgg-20 gene, an EXPB1 protein and a cuticular collagen. By 12 h, there were 7 induced genes and 0 suppressed genes in NL1-RHg. By 3d, there were 9 induced and 10 suppressed genes in NL1-RHg. Substantial changes in gene expression became evident subsequently. At 8d there were 13 induced genes in NL1-RHg. This included putative gland protein G20E03, ubiquitin extension protein, putative gland protein G30C02 and beta-1,4 endoglucanase. However, 1668 genes were found to be suppressed in NL1-RHg. These genes included steroid alpha reductase, serine proteinase and a collagen protein.

CONCLUSION

These analyses identify a genetic expression signature for these two populations both prior to and subsequently as they undergo an R or S reaction. The identification of genes like steroid alpha reductase and serine proteinase that are involved in feeding and nutritional uptake as being highly suppressed during the R response at 8d may indicate genes that the plant is targeting. The analyses also identified numerous putative parasitism genes that are differentially expressed. The 1668 genes that are suppressed in NL1-RHg, and hence induced in TN8 may represent genes that are important during the parasitic stages of H. glycines development. The potential for different arrays of putative parasitism genes to be expressed in different nematode populations may indicate how H. glycines evolve mechanisms to overcome resistance.

Expression patterns in soybean resistant to Phakopsora pachyrhizi reveal the importance of peroxidases and lipoxygenases

Soybean rust caused by Phakopsora pachyrhizi Sydow is a devastating foliar disease that has spread to most soybean growing regions throughout the world, including the USA. Four independent rust resistance genes, Rpp1-Rpp4, have been identified in soybean that recognize specific isolates of P. pachyrhizi. A suppressive subtraction hybridization (SSH) complementary DNA (cDNA) library was constructed from the soybean accession PI200492, which contains Rpp1, after inoculation with two different isolates of P. pachyrhizi that result in susceptible or immune reactions. Both forward and reverse SSH were performed using cDNA from messenger RNA pooled from 1, 6, 12, 24, and 48 h post-inoculation. A total of 1,728 SSH clones were sequenced and compared to sequences in GenBank for similarity. Microarray analyses were conducted on a custom 7883 soybean-cDNA clone array encompassing all of the soybean-rust SSH clones and expressed sequence tags from four other soybean cDNA libraries. Results of the microarray revealed 558 cDNA clones differentially expressed in the immune reaction. The majority of the upregulated cDNA clones fell into the functional category of defense. In particular, cDNA clones with similarity to peroxidases and lipoxygenases were prevalent. Downregulated cDNA clones included those with similarity to cell-wall-associated protein, such as extensins, proline-rich proteins, and xyloglucan endotransglycosylases.

A decline in transcript abundance for Heterodera glycines homologs of Caenorhabditis elegans uncoordinated genes accompanies its sedentary parasitic phase

Background

Heterodera glycines (soybean cyst nematode [SCN]), the major pathogen of Glycine max (soybean), undergoes muscle degradation (sarcopenia) as it becomes sedentary inside the root. Many genes encoding muscular and neuromuscular components belong to the uncoordinated (unc) family of genes originally identified in Caenorhabditis elegans. Previously, we reported a substantial decrease in transcript abundance for Hg-unc-87, the H. glycines homolog of unc-87 (calponin) during the adult sedentary phase of SCN. These observations implied that changes in the expression of specific muscle genes occurred during sarcopenia.

Results

We developed a bioinformatics database that compares expressed sequence tag (est) and genomic data of C. elegans and H. glycines (CeHg database). We identify H. glycines homologs of C. elegans unc genes whose protein products are involved in muscle composition and regulation. RT-PCR reveals the transcript abundance of H. glycines unc homologs at mobile and sedentary stages of its lifecycle. A prominent reduction in transcript abundance occurs in samples from sedentary nematodes for homologs of actin, unc-60B (cofilin), unc-89, unc-15 (paromyosin), unc-27 (troponin I), unc-54 (myosin), and the potassium channel unc-110 (twk-18). Less reduction is observed for the focal adhesion complex gene Hg-unc-97.

Conclusion

The CeHg bioinformatics database is shown to be useful in identifying homologs of genes whose protein products perform roles in specific aspects of H. glycines muscle biology. Our bioinformatics comparison of C. elegans and H. glycines genomic data and our Hg-unc-87 expression experiments demonstrate that the transcript abundance of specific H. glycines homologs of muscle gene decreases as the nematode becomes sedentary inside the root during its parasitic feeding stages.

Identification of Heterodera glycines (soybean cyst nematode [SCN]) cDNA sequences with high identity to those of Caenorhabditis elegans having lethal mutant or RNAi phenotypes

The soybean cyst nematode (SCN; Heterodera glycines) is a devastating obligate parasite of Glycine max (soybean) causing one billion dollars in losses to the US economy per year and over ten billion dollars in losses worldwide. While much is understood about the pathology of H. glycines, its genome sequence is not well characterized or fully sequenced. We sought to create bioinformatic tools to mine the H. glycines nucleotide database. One way is to use a comparative genomics approach by anchoring our analysis with an organism, like the free-living nematode Caenorhabditis elegans. Unlike H. glycines, the C. elegans genome is fully sequenced and is well characterized with a number of lethal genes identified through experimental methods. We compared an EST database of H. glycines with the C. elegans genome. Our goal was identifying genes that may be essential for H. glycines survival and would serve as an automated pipeline for RNAi studies to both study and control H. glycines. Our analysis yielded a total of nearly 8334 conserved genes between H. glycines and C. elegans. Of these, 1508 have lethal phenotypes/phenocopies in C. elegans. RNAi of a conserved ribosomal gene from H. glycines (Hg-rps-23) yielded dead and dying worms as shown by positive Sytox fluorescence. Endogenous Hg-rps-23 exhibited typical RNA silencing as shown by RT-PCR. However, an unrelated gene Hg-unc-87 did not exhibit RNA silencing in the Hg-rps-23 dsRNA-treated worms, demonstrating the specificity of the silencing.

Divergent evolution of arrested development in the dauer stage of Caenorhabditis elegans and the infective stage of Heterodera glycines

BACKGROUND

The soybean cyst nematode Heterodera glycines is the most important parasite in soybean production worldwide. A comprehensive analysis of large-scale gene expression changes throughout the development of plant-parasitic nematodes has been lacking to date.

RESULTS

We report an extensive genomic analysis of H. glycines, beginning with the generation of 20,100 expressed sequence tags (ESTs). In-depth analysis of these ESTs plus approximately 1,900 previously published sequences predicted 6,860 unique H. glycines genes and allowed a classification by function using InterProScan. Expression profiling of all 6,860 genes throughout the H. glycines life cycle was undertaken using the Affymetrix Soybean Genome Array GeneChip. Our data sets and results represent a comprehensive resource for molecular studies of H. glycines. Demonstrating the power of this resource, we were able to address whether arrested development in the Caenorhabditis elegans dauer larva and the H. glycines infective second-stage juvenile (J2) exhibits shared gene expression profiles. We determined that the gene expression profiles associated with the C. elegans dauer pathway are not uniformly conserved in H. glycines and that the expression profiles of genes for metabolic enzymes of C. elegans dauer larvae and H. glycines infective J2 are dissimilar.

CONCLUSION

Our results indicate that hallmark gene expression patterns and metabolism features are not shared in the developmentally arrested life stages of C. elegans and H. glycines, suggesting that developmental arrest in these two nematode species has undergone more divergent evolution than previously thought and pointing to the need for detailed genomic analyses of individual parasite species.

Timecourse microarray analyses reveal global changes in gene expression of susceptible Glycine max (soybean) roots during infection by Heterodera glycines (soybean cyst nematode)

Changes in gene expression within roots of Glycine max (soybean), cv. Kent, susceptible to infection by Heterodera glycines (the soybean cyst nematode [SCN]), at 6, 12, and 24 h, and 2, 4, 6, and 8 days post-inoculation were monitored using microarrays containing more than 6,000 cDNA inserts. Replicate, independent biological samples were examined at each time point. Gene expression was analyzed statistically using T-tests, ANOVA, clustering algorithms, and online analytical processing (OLAP). These analyses allow the user to query the data in several ways without importing the data into third-party software. RT-PCR confirmed that WRKY6 transcription factor, trehalose phosphate synthase, EIF4a, Skp1, and CLB1 were differentially induced across most time-points. Other genes induced across most timepoints included lipoxygenase, calmodulin, phospholipase C, metallothionein-like protein, and chalcone reductase. RT-PCR demonstrated enhanced expression during the first 12 h of infection for Kunitz trypsin inhibitor and sucrose synthase. The stress-related gene, SAM-22, phospholipase D and 12-oxophytodienoate reductase were also induced at the early time-points. At 6 and 8 dpi there was an abundance of transcripts expressed that encoded genes involved in transcription and protein synthesis. Some of those genes included ribosomal proteins, and initiation and elongation factors. Several genes involved in carbon metabolism and transport were also more abundant. Those genes included glyceraldehyde 3-phosphate dehydrogenase, fructose-bisphosphate aldolase and sucrose synthase. These results identified specific changes in gene transcript levels triggered by infection of susceptible soybean roots by SCN.

Laser capture microdissection (LCM) and expression analyses of Glycine max (soybean) syncytium containing root regions formed by the plant pathogen Heterodera glycines (soybean cyst nematode)

Roots of soybean, Glycine max cv. Kent L. Merr., plants susceptible to the soybean cyst nematode (SCN), Heterodera glycines Ichinohe, were inoculated and allowed to develop feeding sites (syncytia) for 8 days. Root samples enriched in syncytial cells were collected using laser capture microdissection (LCM). RNA was extracted and used to make a cDNA library and expressed sequence tags (ESTs) were produced and used for a Gene Ontology (GO) analysis. RT-PCR results indicated enhanced expression of an aquaporin (GmPIP2,2), alpha-tubulin (GmTubA1), beta-tubulin (GmTubB4) and several other genes in syncytium-enriched samples as compared to samples extracted from whole roots. While RT-PCR data showed increased transcript levels of GmPIP2,2 from LCM tissue enriched in syncytial cells, in situ hybridization showed prominent GmPIP2,2 hybridization to RNA in the parenchymal cells tightly juxtaposed to the syncytium. Immunolocalization indicated stronger alpha-tubulin signal within the syncytium as compared to surrounding tissue. However, alpha-tubulin labeling appeared diffuse or clumped. Thus, LCM allowed for the isolation of tissue enriched for syncytial cells, providing material suitable for a variety of molecular analyses.

SGMD: the Soybean Genomics and Microarray Database

The Soybean Genomics and Microarray Database (SGMD) attempts to provide an integrated view of the interaction of soybean with the soybean cyst nematode and contains genomic, EST and microarray data with embedded analytical tools allowing correlation of soybean ESTs with their gene expression profiles. SGMD provides analytical tools to mine the microarray data quickly by integrating many analysis methods within the database itself. The expression profiles of genes at time intervals during the first 8 days of nematode invasion is searchable by gene name or GenBank accession number. Recent developments include the addition of a searchable database for soybean cyst nematode ESTs and photographs of the invasion process at time points examined using microarrays. SGMD is completely accessible from the web at: http://psi081.ba.ars.usda.gov/SGMD/default.htm.