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)

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.

The heterologous expression of conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs suppresses Meloidogyne incognita parasitism in Gossypium hirsutum (upland cotton)

Two conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs function in defense to the parasitic soybean cyst nematode Heterodera glycines. Gene Ontology analyses of RNA seq data obtained from MAPK3-1-overexpressing (OE) and MAPK3-2-OE roots compared to their control, as well as MAPK3-1-RNA interference (RNAi) and MAPK3-2-RNAi compared to their control, hierarchically orders the induced and suppressed genes, strengthening the hypothesis that their heterologous expression in Gossypium hirsutum (upland cotton) would impair parasitism by the root knot nematode (RKN) Meloidogyne incognita. MAPK3-1 expression (E) in G. hirsutum suppresses the production of M. incognita root galls, egg masses, and second stage juveniles (J2s) by 80.32%, 82.37%, and 88.21%, respectfully. Unexpectedly, egg number increases by 28.99% but J2s are inviable. MAPK3-2-E effects are identical, statistically. MAPK3-1-E and MAPK3-2-E decreases root mass 1.49-fold and 1.55-fold, respectively, as compared to the pRAP15-ccdB-E control. The reproductive factor (RF) of M. incognita for G. hirsutum roots expressing MAPK3-1-E or MAPK3-2-E decreases 60.39% and 50.46%, respectively, compared to controls. The results are consistent with upstream pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) functioning in defense to H. glycines. The experiments showcase the feasibility of employing MAPK3, through heterologous expression, to combat M. incognita parasitism, possibly overcoming impediments otherwise making G. hirsutum's defense platform deficient. MAPK homologs are identified in other important crop species for future functional analyses.

Glycine max Homologs of DOESN'T MAKE INFECTIONS 1, 2, and 3 Function to Impair Heterodera glycines Parasitism While Also Regulating Mitogen Activated Protein Kinase Expression

Glycine max root cells developing into syncytia through the parasitic activities of the pathogenic nematode Heterodera glycines underwent isolation by laser microdissection (LM). Microarray analyses have identified the expression of a G. max DOESN'T MAKE INFECTIONS3 (DMI3) homolog in syncytia undergoing parasitism but during a defense response. DMI3 encodes part of the common symbiosis pathway (CSP) involving DMI1, DMI2, and other CSP genes. The identified DMI gene expression, and symbiosis role, suggests the possible existence of commonalities between symbiosis and defense. G. max has 3 DMI1, 12 DMI2, and 2 DMI3 paralogs. LM-assisted gene expression experiments of isolated syncytia under further examination here show G. max DMI1-3, DMI2-7, and DMI3-2 expression occurring during the defense response in the H. glycines-resistant genotypes G.max _{[Peking/PI548402]} and G.max _[PI88788] indicating a broad and consistent level of expression of the genes. Transgenic overexpression (OE) of G. max DMI1-3, DMI2-7, and DMI3-2 impairs H. glycines parasitism. RNA interference (RNAi) of G. max DMI1-3, DMI2-7, and DMI3-2 increases H. glycines parasitism. The combined opposite outcomes reveal a defense function for these genes. Prior functional transgenic analyses of the 32-member G. max mitogen activated protein kinase (MAPK) gene family has determined that 9 of them act in the defense response to H. glycines parasitism, referred to as defense MAPKs. RNA-seq analyses of root RNA isolated from the 9 G. max defense MAPKs undergoing OE or RNAi reveal they alter the relative transcript abundances (RTAs) of specific DMI1, DMI2, and DMI3 paralogs. In contrast, transgenically-manipulated DMI1-3, DMI2-7, and DMI3-2 expression influences MAPK3-1 and MAPK3-2 RTAs under certain circumstances. The results show G. max homologs of the CSP, and defense pathway are linked, apparently involving co-regulated gene expression.

Combining targeted metabolite analyses and transcriptomics to reveal the specific chemical composition and associated genes in the incompatible soybean variety PI437654 infected with soybean cyst nematode HG1.2.3.5.7

BACKGROUND

Soybean cyst nematode, Heterodera glycines, is one of the most devastating pathogens of soybean and causes severe annual yield losses worldwide. Different soybean varieties exhibit different responses to H. glycines infection at various levels, such as the genomic, transcriptional, proteomic and metabolomic levels. However, there have not yet been any reports of the differential responses of incompatible and compatible soybean varieties infected with H. glycines based on combined metabolomic and transcriptomic analyses.

RESULTS

In this study, the incompatible soybean variety PI437654 and three compatible soybean varieties, Williams 82, Zhonghuang 13 and Hefeng 47, were used to clarify the differences in metabolites and transcriptomics before and after the infection with HG1.2.3.5.7. A local metabolite-calibrated database was used to identify potentially differential metabolites, and the differences in metabolites and metabolic pathways were compared between the incompatible and compatible soybean varieties after inoculation with HG1.2.3.5.7. In total, 37 differential metabolites and 20 KEGG metabolic pathways were identified, which were divided into three categories: metabolites/pathways overlapped in the incompatible and compatible soybeans, and metabolites/pathways specific to either the incompatible or compatible soybean varieties. Twelve differential metabolites were found to be involved in predicted KEGG metabolite pathways. Moreover, 14 specific differential metabolites (such as significantly up-regulated nicotine and down-regulated D-aspartic acid) and their associated KEGG pathways (such as the tropane, piperidine and pyridine alkaloid biosynthesis, alanine, aspartate and glutamate metabolism, sphingolipid metabolism and arginine biosynthesis) were significantly altered and abundantly enriched in the incompatible soybean variety PI437654, and likely played pivotal roles in defending against HG1.2.3.5.7 infection. Three key metabolites (N-acetyltranexamic acid, nicotine and D,L-tryptophan) found to be significantly up-regulated in the incompatible soybean variety PI437654 infected by HG1.2.3.5.7 were classified into two types and used for combined analyses with the transcriptomic expression profiling. Associated genes were predicted, along with the likely corresponding biological processes, cellular components, molecular functions and pathways.

CONCLUSIONS

Our results not only identified potential novel metabolites and associated genes involved in the incompatible response of PI437654 to soybean cyst nematode HG1.2.3.5.7, but also provided new insights into the interactions between soybeans and soybean cyst nematodes.

Mitogen activated protein kinase (MAPK)-regulated genes with predicted signal peptides function in the Glycine max defense response to the root pathogenic nematode Heterodera glycines

Glycine max has 32 mitogen activated protein kinases (MAPKs), nine of them exhibiting defense functions (defense MAPKs) to the plant parasitic nematode Heterodera glycines. RNA seq analyses of transgenic G. max lines overexpressing (OE) each defense MAPK has led to the identification of 309 genes that are increased in their relative transcript abundance by all 9 defense MAPKs. Here, 71 of those genes are shown to also have measurable amounts of transcript in H. glycines-induced nurse cells (syncytia) produced in the root that are undergoing a defense response. The 71 genes have been grouped into 7 types, based on their expression profile. Among the 71 genes are 8 putatively-secreted proteins that include a galactose mutarotase-like protein, pollen Ole e 1 allergen and extensin protein, endomembrane protein 70 protein, O-glycosyl hydrolase 17 protein, glycosyl hydrolase 32 protein, FASCICLIN-like arabinogalactan protein 17 precursor, secreted peroxidase and a pathogenesis-related thaumatin protein. Functional transgenic analyses of all 8 of these candidate defense genes that employ their overexpression and RNA interference (RNAi) demonstrate they have a role in defense. Overexpression experiments that increase the relative transcript abundance of the candidate defense gene reduces the ability that the plant parasitic nematode Heterodera glycines has in completing its life cycle while, in contrast, RNAi of these genes leads to an increase in parasitism. The results provide a genomic analysis of the importance of MAPK signaling in relation to the secretion apparatus during the defense process defense in the G. max-H. glycines pathosystem and identify additional targets for future studies.

Exocyst components promote an incompatible interaction between Glycine max (soybean) and Heterodera glycines (the soybean cyst nematode)

Vesicle and target membrane fusion involves tethering, docking and fusion. The GTPase SECRETORY4 (SEC4) positions the exocyst complex during vesicle membrane tethering, facilitating docking and fusion. Glycine max (soybean) Sec4 functions in the root during its defense against the parasitic nematode Heterodera glycines as it attempts to develop a multinucleate nurse cell (syncytium) serving to nourish the nematode over its 30-day life cycle. Results indicate that other tethering proteins are also important for defense. The G. max exocyst is encoded by 61 genes: 5 EXOC1 (Sec3), 2 EXOC2 (Sec5), 5 EXOC3 (Sec6), 2 EXOC4 (Sec8), 2 EXOC5 (Sec10) 6 EXOC6 (Sec15), 31 EXOC7 (Exo70) and 8 EXOC8 (Exo84) genes. At least one member of each gene family is expressed within the syncytium during the defense response. Syncytium-expressed exocyst genes function in defense while some are under transcriptional regulation by mitogen-activated protein kinases (MAPKs). The exocyst component EXOC7-H4-1 is not expressed within the syncytium but functions in defense and is under MAPK regulation. The tethering stage of vesicle transport has been demonstrated to play an important role in defense in the G. max-H. glycines pathosystem, with some of the spatially and temporally regulated exocyst components under transcriptional control by MAPKs.

The mitogen activated protein kinase (MAPK) gene family functions as a cohort during the Glycine max defense response to Heterodera glycines

Mitogen activated protein kinases (MAPKs) play important signal transduction roles. However, little is known regarding how they influence the gene expression of other family members and the relationship to a biological process, including the Glycine max defense response to Heterodera glycines. Transcriptomics have identified MAPK gene expression occurring within root cells undergoing a defense response to a pathogenic event initiated by H. glycines in the allotetraploid Glycine max. Functional analyses are presented for its 32 MAPKs revealing 9 have a defense role, including homologs of Arabidopsis thaliana MAPK (MPK) MPK2, MPK3, MPK4, MPK5, MPK6, MPK13, MPK16 and MPK20. Defense signaling occurring through pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) have been determined in relation to these MAPKs. Five different types of gene expression relate to MAPK expression, influencing PTI and ETI gene expression and proven defense genes including an ABC-G transporter, 20S membrane fusion particle components, glycoside biosynthesis, carbon metabolism, hemicellulose modification, transcription and secretion. The experiments show MAPKs broadly influence defense MAPK gene expression, including the co-regulation of parologous MAPKs and reveal its relationship to proven defense genes. The experiments reveal each defense MAPK induces the expression of a G. max homolog of a PATHOGENESIS RELATED1 (PR1), itself shown to function in defense in the studied pathosystem.

Harpin-inducible defense signaling components impair infection by the ascomycete Macrophomina phaseolina

Soybean (Glycine max) infection by the charcoal rot (CR) ascomycete Macrophomina phaseolina is enhanced by the soybean cyst nematode (SCN) Heterodera glycines. We hypothesized that G. max genetic lines impairing infection by M. phaseolina would also limit H. glycines parasitism, leading to resistance. As a part of this M. phaseolina resistance process, the genetic line would express defense genes already proven to impair nematode parasitism. Using G. max[DT97-4290/PI 642055], exhibiting partial resistance to M. phaseolina, experiments show the genetic line also impairs H. glycines parasitism. Furthermore, comparative studies show G. max[DT97-4290/PI 642055] exhibits induced expression of the effector triggered immunity (ETI) gene NON-RACE SPECIFIC DISEASE RESISTANCE 1/HARPIN INDUCED1 (NDR1/HIN1) that functions in defense to H. glycines as compared to the H. glycines and M. phaseolina susceptible line G. max[Williams 82/PI 518671]. Other defense genes that are induced in G. max[DT97-4290/PI 642055] include the pathogen associated molecular pattern (PAMP) triggered immunity (PTI) genes ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1), NONEXPRESSOR OF PR1 (NPR1) and TGA2. These observations link G. max defense processes that impede H. glycines parasitism to also potentially function toward impairing M. phaseolina pathogenicity. Testing this hypothesis, G. max[Williams 82/PI 518671] genetically engineered to experimentally induce GmNDR1-1, EDS1-2, NPR1-2 and TGA2-1 expression leads to impaired M. phaseolina pathogenicity. In contrast, G. max[DT97-4290/PI 642055] engineered to experimentally suppress the expression of GmNDR1-1, EDS1-2, NPR1-2 and TGA2-1 by RNA interference (RNAi) enhances M. phaseolina pathogenicity. The results show components of PTI and ETI impair both nematode and M. phaseolina pathogenicity.

Identification of Transcriptional and Receptor Networks That Control Root Responses to Ethylene

Transcriptomic analyses with high temporal resolution provide substantial new insight into hormonal response networks. This study identified the kinetics of genome-wide transcript abundance changes in response to elevated levels of the plant hormone ethylene in roots from light-grown Arabidopsis (Arabidopsis thaliana) seedlings, which were overlaid on time-matched developmental changes. Functional annotation of clusters of transcripts with similar temporal patterns revealed rapidly induced clusters with known ethylene function and more slowly regulated clusters with novel predicted functions linked to root development. In contrast to studies with dark-grown seedlings, where the canonical ethylene response transcription factor, EIN3, is central to ethylene-mediated development, the roots of ein3 and eil1 single and double mutants still respond to ethylene in light-grown seedlings. Additionally, a subset of these clusters of ethylene-responsive transcripts were enriched in targets of EIN3 and ERFs. These results are consistent with EIN3-independent developmental and transcriptional changes in light-grown roots. Examination of single and multiple gain-of-function and loss-of-function receptor mutants revealed that, of the five ethylene receptors, ETR1 controls lateral root and root hair initiation and elongation and the synthesis of other receptors. These results provide new insight into the transcriptional and developmental responses to ethylene in light-grown seedlings.

A harpin elicitor induces the expression of a coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene and others functioning during defense to parasitic nematodes

The bacterial effector harpin induces the transcription of the Arabidopsis thaliana NON-RACE SPECIFIC DISEASE RESISTANCE 1/HARPIN INDUCED1 (NDR1/HIN1) coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene. In Glycine max, Gm-NDR1-1 transcripts have been detected within root cells undergoing a natural resistant reaction to parasitism by the syncytium-forming nematode Heterodera glycines, functioning in the defense response. Expressing Gm-NDR1-1 in Gossypium hirsutum leads to resistance to Meloidogyne incognita parasitism. In experiments presented here, the heterologous expression of Gm-NDR1-1 in G. hirsutum impairs Rotylenchulus reniformis parasitism. These results are consistent with the hypothesis that Gm-NDR1-1 expression functions broadly in generating a defense response. To examine a possible relationship with harpin, G. max plants topically treated with harpin result in induction of the transcription of Gm-NDR1-1. The result indicates the topical treatment of plants with harpin, itself, may lead to impaired nematode parasitism. Topical harpin treatments are shown to impair G. max parasitism by H. glycines, M. incognita and R. reniformis and G. hirsutum parasitism by M. incognita and R. reniformis. How harpin could function in defense has been examined in experiments showing it also induces transcription of G. max homologs of the proven defense genes ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1), TGA2, galactinol synthase, reticuline oxidase, xyloglucan endotransglycosylase/hydrolase, alpha soluble N-ethylmaleimide-sensitive fusion protein (α-SNAP) and serine hydroxymethyltransferase (SHMT). In contrast, other defense genes are not directly transcriptionally activated by harpin. The results indicate harpin induces pathogen associated molecular pattern (PAMP) triggered immunity (PTI) and effector-triggered immunity (ETI) defense processes in the root, activating defense to parasitic nematodes.

A Glycine max homolog of NON-RACE SPECIFIC DISEASE RESISTANCE 1 (NDR1) alters defense gene expression while functioning during a resistance response to different root pathogens in different genetic backgrounds

A Glycine max homolog of the Arabidopsis thaliana NON-RACE SPECIFIC DISEASE RESISTANCE 1 (NDR1) coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene (Gm-NDR1-1) is expressed in root cells undergoing a defense response to the root pathogenic nematode, Heterodera glycines. Gm-NDR1-1 overexpression in the H. glycines-susceptible genotype G. max[Williams 82/PI 518671] impairs parasitism. In contrast, Gm-NDR1-1 RNA interference (RNAi) in the H. glycines-resistant genotype G. max[Peking/PI 548402] facilitates parasitism. The broad effectiveness of Gm-NDR1-1 in impairing parasitism has then been examined by engineering its heterologous expression in Gossypium hirsutum which is susceptible to the root pathogenic nematode Meloidogyne incognita. The heterologous expression of Gm-NDR1-1 in G. hirsutum effectively impairs M. incognita parasitism, reducing gall, egg mass, egg and juvenile numbers. In contrast to our prior experiments examining the effectiveness of the heterologous expression of a G. max homolog of the A. thaliana salicyclic acid signaling (SA) gene NONEXPRESSOR OF PR1 (Gm-NPR1-2), no cumulative negative effect on M. incognita parasitism has been observed in G. hirsutum expressing Gm-NDR1-1. The results indicate a common genetic basis exists for plant resistance to parasitic nematodes that involves Gm-NDR1. However, the Gm-NDR1-1 functions in ways that are measurably dissimilar to Gm-NPR1-2. Notably, Gm-NDR1-1 overexpression leads to increased relative transcript levels of its homologs of A. thaliana genes functioning in SA signaling, including NPR1-2, TGA2-1 and LESION SIMULATING DISEASE1 (LSD1-2) that is lost in Gm-NDR1-1 RNAi lines. Similar observations have been made regarding the expression of other defense genes.

Components of the SNARE-containing regulon are co-regulated in root cells undergoing defense

The term regulon has been coined in the genetic model plant Arabidopsis thaliana, denoting a structural and physiological defense apparatus defined genetically through the identification of the penetration (pen) mutants. The regulon is composed partially by the soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) syntaxin PEN1. PEN1 has homology to a Saccharomyces cerevisae gene that regulates a Secretion (Sec) protein, Suppressor of Sec 1 (Sso1p). The regulon is also composed of the β-glucosidase (PEN2) and an ATP binding cassette (ABC) transporter (PEN3). While important in inhibiting pathogen infection, limited observations have been made regarding the transcriptional regulation of regulon genes until now. Experiments made using the model agricultural Glycine max (soybean) have identified co-regulated gene expression of regulon components. The results explain the observation of hundreds of genes expressed specifically in the root cells undergoing the natural process of defense. Data regarding additional G. max genes functioning within the context of the regulon are presented here, including Sec 14, Sec 4 and Sec 23. Other examined G. max homologs of membrane fusion genes include an endosomal bromo domain-containing protein1 (Bro1), syntaxin6 (SYP6), SYP131, SYP71, SYP8, Bet1, coatomer epsilon (ϵ-COP), a coatomer zeta (ζ-COP) paralog and an ER to Golgi component (ERGIC) protein. Furthermore, the effectiveness of biochemical pathways that would function within the context of the regulon ave been examined, including xyloglucan xylosyltransferase (XXT), reticuline oxidase (RO) and galactinol synthase (GS). The experiments have unveiled the importance of the regulon during defense in the root and show how the deposition of callose relates to the process.

Transcriptome Analysis of Cotton (Gossypium hirsutum L.) Genotypes That Are Susceptible, Resistant, and Hypersensitive to Reniform Nematode (Rotylenchulus reniformis)

Reniform nematode is a semi-endoparasitic nematode species causing significant yield loss in numerous crops, including cotton (Gossypium hirsutum L.). An RNA-sequencing analysis was conducted to measure transcript abundance in reniform nematode susceptible (DP90 & SG747), resistant (BARBREN-713), and hypersensitive (LONREN-1) genotypes of cotton (Gossypium hirsutum L.) with and without reniform nematode infestation. Over 90 million trimmed high quality reads were assembled into 84,711 and 80, 353 transcripts using the G. arboreum and the G. raimondii genomes as references. Many transcripts were significantly differentially expressed between the three different genotypes both prior to and during nematode pathogenesis, including transcripts corresponding to the gene ontology categories of cell wall, hormone metabolism and signaling, redox reactions, secondary metabolism, transcriptional regulation, stress responses, and signaling. Further analysis revealed that a number of these differentially expressed transcripts mapped to the G. raimondii and/or the G. arboreum genomes within 1 megabase of quantitative trait loci that had previously been linked to reniform nematode resistance. Several resistance genes encoding proteins known to be strongly linked to pathogen perception and resistance, including LRR-like and NBS-LRR domain-containing proteins, were among the differentially expressed transcripts mapping near these quantitative trait loci. Further investigation is required to confirm a role for these transcripts in reniform nematode susceptibility, hypersensitivity, and/or resistance. This study presents the first systemic investigation of reniform nematode resistance-associated genes using different genotypes of cotton. The candidate reniform nematode resistance-associated genes identified in this study can serve as the basis for further functional analysis and aid in further development of reniform a nematode resistant cotton germplasm.

Whole-genome gene expression profiling revealed genes and pathways potentially involved in regulating interactions of soybean with cyst nematode (Heterodera glycines Ichinohe)

BACKGROUND

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is the most devastating pathogen of soybean. Many gene expression profiling studies have been conducted to investigate the responses of soybean to the infection by this pathogen using primarily the first-generation soybean genome array that covered approximately 37,500 soybean transcripts. However, no study has been reported yet using the second-generation Affymetrix soybean whole-genome transcript array (Soybean WT array) that represents approximately 66,000 predicted soybean transcripts.

RESULTS

In the present work, the gene expression profiles of two soybean plant introductions (PIs) PI 437654 and PI 567516C (both resistant to multiple SCN HG Types) and cultivar Magellan (susceptible to SCN) were compared in the presence or absence of the SCN inoculum at 3 and 8 days post-inoculation using the Soybean WT array. Data analysis revealed that the two resistant soybean lines showed distinctive gene expression profiles from each other and from Magellan not only in response to the SCN inoculation, but also in the absence of SCN. Overall, 1,413 genes and many pathways were revealed to be differentially regulated. Among them, 297 genes were constitutively regulated in the two resistant lines (compared with Magellan) and 1,146 genes were responsive to the SCN inoculation in the three lines, with 30 genes regulated both constitutively and by SCN. In addition to the findings similar to those in the published work, many genes involved in ethylene, protein degradation, and phenylpropanoid pathways were also revealed differentially regulated in the present study. GC-rich elements (e.g., GCATGC) were found over-represented in the promoter regions of certain groups of genes. These have not been observed before, and could be new defense-responsive regulatory elements.

CONCLUSIONS

Different soybean lines showed different gene expression profiles in the presence and absence of the SCN inoculum. Both inducible and constitutive gene expression may contribute to resistance to multiple SCN HG Types in the resistant soybean PI lines. Ethylene, protein degradation, and phenylpropanoid pathways, as well as many other pathways reported previously, may play important roles in mediating the soybean-SCN interactions. The revealed genes, pathways, and promoter elements can be further explored to regulate or engineer soybean for resistance to SCN.

Syntaxin 31 functions in Glycine max resistance to the plant parasitic nematode Heterodera glycines

A Glycine max syntaxin 31 homolog (Gm-SYP38) was identified as being expressed in nematode-induced feeding structures known as syncytia undergoing an incompatible interaction with the plant parasitic nematode Heterodera glycines. The observed Gm-SYP38 expression was consistent with prior gene expression analyses that identified the alpha soluble NSF attachment protein (Gm-α-SNAP) resistance gene because homologs of these genes physically interact and function together in other genetic systems. Syntaxin 31 is a protein that resides on the cis face of the Golgi apparatus and binds α-SNAP-like proteins, but has no known role in resistance. Experiments presented here show Gm-α-SNAP overexpression induces Gm-SYP38 transcription. Overexpression of Gm-SYP38 rescues G. max [Williams 82/PI 518671], genetically rhg1 (-/-), by suppressing H. glycines parasitism. In contrast, Gm-SYP38 RNAi in the rhg1 (+/+) genotype G. max [Peking/PI 548402] increases susceptibility. Gm-α-SNAP and Gm-SYP38 overexpression induce the transcriptional activity of the cytoplasmic receptor-like kinase BOTRYTIS INDUCED KINASE 1 (Gm-BIK1-6) which is a family of defense proteins known to anchor to membranes through a 5' MGXXXS/T(R) N-myristoylation sequence. Gm-BIK1-6 had been identified previously by RNA-seq experiments as expressed in syncytia undergoing an incompatible reaction. Gm-BIK1-6 overexpression rescues the resistant phenotype. In contrast, Gm-BIK1-6 RNAi increases parasitism. The analysis demonstrates a role for syntaxin 31-like genes in resistance that until now was not known.

Overexpression of a soybean salicylic acid methyltransferase gene confers resistance to soybean cyst nematode

Salicylic acid plays a critical role in activating plant defence responses after pathogen attack. Salicylic acid methyltransferase (SAMT) modulates the level of salicylic acid by converting salicylic acid to methyl salicylate. Here, we report that a SAMT gene from soybean (GmSAMT1) plays a role in soybean defence against soybean cyst nematode (Heterodera glycines Ichinohe, SCN). GmSAMT1 was identified as a candidate SCN defence-related gene in our previous analysis of soybean defence against SCN using GeneChip microarray experiments. The current study started with the isolation of the full-length cDNAs of GmSAMT1 from a SCN-resistant soybean line and from a SCN-susceptible soybean line. The two cDNAs encode proteins of identical sequences. The GmSAMT1 cDNA was expressed in Escherichia coli. Using in vitro enzyme assays, E. coli-expressed GmSAMT1 was confirmed to function as salicylic acid methyltransferase. The apparent Km value of GmSAMT1 for salicylic acid was approximately 46 μM. To determine the role of GmSAMT1 in soybean defence against SCN, transgenic hairy roots overexpressing GmSAMT1 were produced and tested for SCN resistance. Overexpression of GmSAMT1 in SCN-susceptible backgrounds significantly reduced the development of SCN, indicating that overexpression of GmSAMT1 in the transgenic hairy root system could confer resistance to SCN. Overexpression of GmSAMT1 in transgenic hairy roots was also found to affect the expression of selected genes involved in salicylic acid biosynthesis and salicylic acid signal transduction.

A kinetic analysis of the auxin transcriptome reveals cell wall remodeling proteins that modulate lateral root development in Arabidopsis

To identify gene products that participate in auxin-dependent lateral root formation, a high temporal resolution, genome-wide transcript abundance analysis was performed with auxin-treated Arabidopsis thaliana roots. Data analysis identified 1246 transcripts that were consistently regulated by indole-3-acetic acid (IAA), partitioning into 60 clusters with distinct response kinetics. We identified rapidly induced clusters containing auxin-response functional annotations and clusters exhibiting delayed induction linked to cell division temporally correlated with lateral root induction. Several clusters were enriched with genes encoding proteins involved in cell wall modification, opening the possibility for understanding mechanistic details of cell structural changes that result in root formation following auxin treatment. Mutants with insertions in 72 genes annotated with a cell wall remodeling function were examined for alterations in IAA-regulated root growth and development. This reverse-genetic screen yielded eight mutants with root phenotypes. Detailed characterization of seedlings with mutations in cellulase3/glycosylhydrolase9b3 and leucine rich extensin2, genes not normally linked to auxin response, revealed defects in the early and late stages of lateral root development, respectively. The genes identified here using kinetic insight into expression changes lay the foundation for mechanistic understanding of auxin-mediated cell wall remodeling as an essential feature of lateral root development.

Distinct and conserved transcriptomic changes during nematode-induced giant cell development in tomato compared with Arabidopsis: a functional role for gene repression

Root-knot nematodes (RKNs) induce giant cells (GCs) from root vascular cells inside the galls. Accompanying molecular changes as a function of infection time and across different species, and their functional impact, are still poorly understood. Thus, the transcriptomes of tomato galls and laser capture microdissected (LCM) GCs over the course of parasitism were compared with those of Arabidopsis, and functional analysis of a repressed gene was performed. Microarray hybridization with RNA from galls and LCM GCs, infection-reproduction tests and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) transcriptional profiles in susceptible and resistant (Mi-1) lines were performed in tomato. Tomato GC-induced genes include some possibly contributing to the epigenetic control of GC identity. GC-repressed genes are conserved between tomato and Arabidopsis, notably those involved in lignin deposition. However, genes related to the regulation of gene expression diverge, suggesting that diverse transcriptional regulators mediate common responses leading to GC formation in different plant species. TPX1, a cell wall peroxidase specifically involved in lignification, was strongly repressed in GCs/galls, but induced in a nearly isogenic Mi-1 resistant line on nematode infection. TPX1 overexpression in susceptible plants hindered nematode reproduction and GC expansion. Time-course and cross-species comparisons of gall and GC transcriptomes provide novel insights pointing to the relevance of gene repression during RKN establishment.

The expression of a naturally occurring, truncated allele of an α-SNAP gene suppresses plant parasitic nematode infection

Transcriptional mapping experiments of the major soybean cyst nematode resistance locus, rhg1, identified expression of the vesicular transport machinery component, α soluble NSF attachment protein (α-SNAP), occurring during defense. Sequencing the α-SNAP coding regions from the resistant genotypes G. max ([Peking/PI 548402]) and G. max ([PI 437654]) revealed they are identical, but differ from the susceptible G. max ([Williams 82/PI 518671]) by the presence of several single nucleotide polymorphisms. Using G. max ([Williams 82/PI 518671]) as a reference, a G → T(2,822) transversion in the genomic DNA sequence at a functional splice site of the α-SNAP([Peking/PI 548402]) allele produced an additional 17 nucleotides of mRNA sequence that contains an in-frame stop codon caused by a downstream G → A(2,832) transition. The G. max ([Peking/PI 548402]) genotype has cell wall appositions (CWAs), structures identified as forming as part of a defense response by the activity of the vesicular transport machinery. In contrast, the 17 nt α-SNAP([Peking/PI 548402]) mRNA motif is not found in G. max ([PI 88788]) that exhibits defense to H. glycines, but lack CWAs. The α-SNAP([PI 88788]) promoter contains sequence elements that are nearly identical to the α-SNAP([Peking/PI 548402]) allele, but differs from the G. max ([Williams 82/PI 518671]) ortholog. Overexpressing the α-SNAP([Peking/PI 548402]) allele in the susceptible G. max ([Williams 82/PI 518671]) genotype suppressed H. glycines infection. The experiments indicate a role for the vesicular transport machinery during infection of soybean by the soybean cyst nematode. However, increased GmEREBP1, PR1, PR2, PR5 gene activity but suppressed PR3 expression accompanied the overexpression of the α-SNAP([Peking/PI 548402]) allele prior to infection.

Biotechnological application of functional genomics towards plant-parasitic nematode control

Plant-parasitic nematodes are primary biotic factors limiting the crop production. Current nematode control strategies include nematicides, crop rotation and resistant cultivars, but each has serious limitations. RNA interference (RNAi) represents a major breakthrough in the application of functional genomics for plant-parasitic nematode control. RNAi-induced suppression of numerous genes essential for nematode development, reproduction or parasitism has been demonstrated, highlighting the considerable potential for using this strategy to control damaging pest populations. In an effort to find more suitable and effective gene targets for silencing, researchers are employing functional genomics methodologies, including genome sequencing and transcriptome profiling. Microarrays have been used for studying the interactions between nematodes and plant roots and to measure both plants and nematodes transcripts. Furthermore, laser capture microdissection has been applied for the precise dissection of nematode feeding sites (syncytia) to allow the study of gene expression specifically in syncytia. In the near future, small RNA sequencing techniques will provide more direct information for elucidating small RNA regulatory mechanisms in plants and specific gene silencing using artificial microRNAs should further improve the potential of targeted gene silencing as a strategy for nematode management.

Mapping cell fate decisions that occur during soybean defense responses

The soybean defense response to the soybean cyst nematode was used as a model to map at cellular resolution its genotype-defined cell fate decisions occurring during its resistant reactions. The defense responses occur at the site of infection, a nurse cell known as the syncytium. Two major genotype-defined defense responses exist, the G. max ([Peking])- and G. max ([PI 88788])-types. Resistance in G. max ([Peking]) is potent and rapid, accompanied by the formation of cell wall appositions (CWAs), structures known to perform important defense roles. In contrast, defense occurs by a potent but more prolonged reaction in G. max ([PI 88788]), lacking CWAs. Comparative transcriptomic analyses with confirmation by Illumina® deep sequencing were organized through a custom-developed application, Pathway Analysis and Integrated Coloring of Experiments (PAICE) that presents gene expression of these cytologically and developmentally distinct defense responses using the Kyoto Encyclopedia of Genes and Genomes (KEGG) framework. The analyses resulted in the generation of 1,643 PAICE pathways, allowing better understanding of gene activity across all chromosomes. Analyses of the rhg1 resistance locus, defined within a 67 kb region of DNA demonstrate expression of an amino acid transporter and an α soluble NSF attachment protein gene specifically in syncytia undergoing their defense responses.

Differences in gene expression amplitude overlie a conserved transcriptomic program occurring between the rapid and potent localized resistant reaction at the syncytium of the Glycine max genotype Peking (PI 548402) as compared to the prolonged and potent resistant reaction of PI 88788

Glycine max L. Merr. (soybean) resistance to Heterodera glycines Ichinohe occurs at the site of infection, a nurse cell known as the syncytium. Resistance is classified into two cytologically-defined responses, the G. max ([Peking])- and G. max ([PI 88788])-types. Each type represents a cohort of G. max genotypes. Resistance in G. max ([Peking]) occurs by a potent and rapid localized response, affecting parasitic second stage juveniles (p-J2). In contrast, resistance occurs by a potent but more prolonged reaction in the genotype G. max ([PI 88788]) that affects nematode development at the J3 and J4 stages. Microarray analyses comparing these cytologically and developmentally distinct resistant reactions reveal differences in gene expression in pericycle and surrounding cells even before infection. The differences include higher relative levels of the differentially expressed in response to arachidonic acid 1 gene (DEA1 [Gm-DEA1]) (+224.19-fold) and a protease inhibitor (+68.28-fold) in G. max ([Peking/PI 548402]) as compared to G. max ([PI 88788]). Gene pathway analyses compare the two genotypes (1) before, (2) at various times during, (3) constitutively throughout the resistant reaction and (4) at all time points prior to and during the resistant reaction. The amplified levels of transcriptional activity of defense genes may explain the rapid and potent reaction in G. max ([Peking/PI 548402]) as compared to G. max ([PI 88788]). In contrast, the shared differential expression levels of genes in G. max ([Peking/PI 548402]) and G. max ([PI 88788]) may indicate a conserved genomic program underlying the G. max resistance on which the genotype-specific gene expression programs are built off.