Syncytium gene expression in Glycine max([PI 88788]) roots undergoing a resistant reaction to the parasitic 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.

Biochar application alters soil metabolites and nitrogen cycle-related microorganisms in a soybean continuous cropping system

Biochar application is a promising practice to enhance soil fertility. However, it is unclear how field-aged biochar affects the soil metabolites and microbial communities in soybean fields. Here, the rhizosphere soil performance after amending with biochar addition rates at 0 (CK), 20 (B20), 40 (B40), and 60 t ha-1 (B60) was examined via a five-year in-situ field experiment based on a soybean continuous cropping system. Untargeted metabolomics and metagenomics analysis techniques were applied to study the regulatory mechanism of biochar on soybean growth from metabolomics and N cycle microbiology perspectives. We found that the contents of soil total N (TN), available N (Ava N), NH4+-N, and NO3--N were significantly increased with biochar addition amounts by 20.0-65.7 %, 3.6-10.7 %, 29.5-57.1 %, and 24.4-46.7 %, respectively. The B20, B40, and B60 triggered 259 (236 were up-regulated and 23 were down-regulated), 236 (220 were up-regulated and 16 were down-regulated), and 299 (264 were up-regulated and 35 were down-regulated) differential metabolites, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and topology analysis demonstrated that differential metabolites were highly enriched in seven metabolic pathways such as Oxidative phosphorylation and Benzoxazinoid biosynthesis. Moreover, ten differential metabolites were up-regulated in all three treatments with biochar. Biochar treatments decreased the Nitrospira abundance in soybean rhizosphere soil while increasing Bradyrhizobium abundance significantly in B60. Mantel test revealed that as the biochar addition rate grows, the correlation between Nitrospira and soil properties other than NO3--N became stronger. In conclusion, the co-application of biochar with fertilizers is a feasible and effective way to improve soil N supply, even though biochar has undergone field aging. This work offers new insights into the variations in soil metabolites and microbial communities associated with N metabolism processes under biochar addition in soybean continuous cropping soils.

A genome-wide study of the lipoxygenase gene families in Medicago truncatula and Medicago sativa reveals that MtLOX24 participates in the methyl jasmonate response

BACKGROUND

Lipoxygenase (LOX) is a multifunctional enzyme that is primarily related to plant organ growth and development, biotic and abiotic stress responses, and production of flavor-associated metabolites. In higher plants, the LOX family encompasses several isozymes with varying expression patterns between tissues and developmental stages. These affect processes including seed germination, seed storage, seedling growth, fruit ripening, and leaf senescence. LOX family genes have multiple functions in response to hormones such as methyl jasmonate (MeJA) and salicylic acid.

RESULTS

In this study, we identified 30 and 95 LOX homologs in Medicago truncatula and Medicago sativa, respectively. These genes were characterized with analyses of their basic physical and chemical properties, structures, chromosomal distributions, and phylogenetic relationships to understand structural variations and their physical locations. Phylogenetic analysis was conducted for members of the three LOX subfamilies (9-LOX, type I 13-LOX, and type II 13-LOX) in Arabidopsis thaliana, Glycine max, M. truncatula, and M. sativa. Analysis of predicted promoter elements revealed several relevant cis-acting elements in MtLOX and MsLOX genes, including abscisic acid (ABA) response elements (ABREs), MeJA response elements (CGTCA-motifs), and antioxidant response elements (AREs). Cis-element data combined with transcriptomic data demonstrated that LOX gene family members in these species were most likely related to abiotic stress responses, hormone responses, and plant development. Gene expression patterns were confirmed via quantitative reverse transcription PCR. Several MtLOX genes (namely MtLOX15, MtLOX16, MtLOX20, and MtLOX24) belonging to the type I 13-LOX subfamily and other LOX genes (MtLOX7, MtLOX11, MsLOX23, MsLOX87, MsLOX90, and MsLOX94) showed significantly different expression levels in the flower tissue, suggesting roles in reproductive growth. Type I 13-LOXs (MtLOX16, MtLOX20, MtLOX21, MtLOX24, MsLOX57, MsLOX84, MsLOX85, and MsLOX94) and type II 13-LOXs (MtLOX5, MtLOX6, MtLOX9, MtLOX10, MsLOX18, MsLOX23, and MsLOX30) were MeJA-inducible and were predicted to function in the jasmonic acid signaling pathway. Furthermore, exogenous MtLOX24 expression in Arabidopsis verified that MtLOX24 was involved in MeJA responses, which may be related to insect-induced abiotic stress.

CONCLUSIONS

We identified six and four LOX genes specifically expressed in the flowers of M. truncatula and M. sativa, respectively. Eight and seven LOX genes were induced by MeJA in M. truncatula and M. sativa, and the LOX genes identified were mainly distributed in the type I and type II 13-LOX subfamilies. MtLOX24 was up-regulated at 8 h after MeJA induction, and exogenous expression in Arabidopsis demonstrated that MtLOX24 promoted resistance to MeJA-induced stress. This study provides valuable new information regarding the evolutionary history and functions of LOX genes in the genus Medicago.

Establishment and maintenance of DNA methylation in nematode feeding sites

A growing body of evidence indicates that epigenetic mechanisms, particularly DNA methylation, play key regulatory roles in plant-nematode interactions. Nevertheless, the transcriptional activity of key genes mediating DNA methylation and active demethylation in the nematode feeding sites remains largely unknown. Here, we profiled the promoter activity of 12 genes involved in maintenance and de novo establishment of DNA methylation and active demethylation in the syncytia and galls induced respectively by the cyst nematode Heterodera schachtii and the root-knot nematode Meloidogyne incognita in Arabidopsis roots. The promoter activity assays revealed that expression of the CG-context methyltransferases is restricted to feeding site formation and development stages. Chromomethylase1 (CMT1), CMT2, and CMT3 and Domains Rearranged Methyltransferase2 (DRM2) and DRM3, which mediate non-CG methylation, showed similar and distinct expression patterns in the syncytia and galls at various time points. Notably, the promoters of various DNA demethylases were more active in galls as compared with the syncytia, particularly during the early stage of infection. Mutants impaired in CG or CHH methylation similarly enhanced plant susceptibility to H. schachtii and M. incognita, whereas mutants impaired in CHG methylation reduced plant susceptibility only to M. incognita. Interestingly, hypermethylated mutants defective in active DNA demethylation exhibited contrasting responses to infection by H. schachtii and M. incognita, a finding most likely associated with differential regulation of defense-related genes in these mutants upon nematode infection. Our results point to methylation-dependent mechanisms regulating plant responses to infection by cyst and root-knot nematodes.

Effector-Dependent and -Independent Molecular Mechanisms of Soybean-Microbe Interaction

Soybean is a pivotal staple crop worldwide, supplying the main food and feed plant proteins in some countries. In addition to interacting with mutualistic microbes, soybean also needs to protect itself against pathogens. However, to grow inside plant tissues, plant defense mechanisms ranging from passive barriers to induced defense reactions have to be overcome. Pathogenic but also symbiotic micro-organisms effectors can be delivered into the host cell by secretion systems and can interfere with the immunity system and disrupt cellular processes. This review summarizes the latest advances in our understanding of the interaction between secreted effectors and soybean feedback mechanism and uncovers the conserved and special signaling pathway induced by pathogenic soybean cyst nematode, Pseudomonas, Xanthomonas as well as by symbiotic rhizobium.

Comparisons of constitutive resistances to soybean cyst nematode between PI 88788- and Peking-type sources of resistance in soybean by transcriptomic and metabolomic profilings

Soybean cyst nematode (SCN) is a serious damaging disease in soybean worldwide. Peking- and PI 88788-type sources of resistance are two most important germplasm used in breeding resistant soybean cultivars against this disease. However, until now, no comparisons of constitutive resistances to soybean cyst nematode between these two types of sources had been conducted, probably due to the influences of different backgrounds. In this study, we used pooled-sample analysis strategy to minimize the influence of different backgrounds and directly compared the molecular mechanisms underlying constitutive resistance to soybean cyst nematode between these two types of sources via transcriptomic and metabolomic profilings. Six resistant soybean accessions that have identical haplotypes as Peking at Rgh1 and Rhg4 loci were pooled to represent Peking-type sources. The PI88788-type and control pools were also constructed in a same way. Through transcriptomic and metabolomics anaylses, differentially expressed genes and metabolites were identified. The molecular pathways involved in the metabolism of toxic metabolites were predicted to play important roles in conferring soybean cyst nematode resistance to soybean. Functions of two resistant candidate genes were confirmed by hairy roots transformation methods in soybean. Our studies can be helpful for soybean scientists to further learn about the molecular mechanism of resistance to soybean cyst nematode in soybean.

Integrated Omic Approaches Reveal Molecular Mechanisms of Tolerance during Soybean and Meloidogyne incognita Interactions

The root-knot nematode (RKN), Meloidogyne incognita, is a devastating soybean pathogen worldwide. The use of resistant cultivars is the most effective method to prevent economic losses caused by RKNs. To elucidate the mechanisms involved in resistance to RKN, we determined the proteome and transcriptome profiles from roots of susceptible (BRS133) and highly tolerant (PI 595099) Glycine max genotypes 4, 12, and 30 days after RKN infestation. After in silico analysis, we described major defense molecules and mechanisms considered constitutive responses to nematode infestation, such as mTOR, PI3K-Akt, relaxin, and thermogenesis. The integrated data allowed us to identify protein families and metabolic pathways exclusively regulated in tolerant soybean genotypes. Among them, we highlighted the phenylpropanoid pathway as an early, robust, and systemic defense process capable of controlling M. incognita reproduction. Associated with this metabolic pathway, 29 differentially expressed genes encoding 11 different enzymes were identified, mainly from the flavonoid and derivative pathways. Based on differential expression in transcriptomic and proteomic data, as well as in the expression profile by RT-qPCR, and previous studies, we selected and overexpressed the GmPR10 gene in transgenic tobacco to assess its protective effect against M. incognita. Transgenic plants of the T₂ generation showed up to 58% reduction in the M. incognita reproduction factor. Finally, data suggest that GmPR10 overexpression can be effective against the plant parasitic nematode M. incognita, but its mechanism of action remains unclear. These findings will help develop new engineered soybean genotypes with higher performance in response to RKN infections.

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.

Exploring Soybean Resistance to Soybean Cyst Nematode

Resistance to the soybean cyst nematode (SCN) is a topic incorporating multiple mechanisms and multiple types of science. It is also a topic of substantial agricultural importance, as SCN is estimated to cause more yield damage than any other pathogen of soybean, one of the world's main food crops. Both soybean and SCN have experienced jumps in experimental tractability in the past decade, and significant advances have been made. The rhg1-b locus, deployed on millions of farm acres, has been durable and will remain important, but local SCN populations are gradually evolving to overcome rhg1-b. Multiple other SCN resistance quantitative trait loci (QTL) of proven value are now in play with soybean breeders. QTL causal gene discovery and mechanistic insights into SCN resistance are contributing to both basic and applied disciplines. Additional understanding of SCN and other cyst nematodes will also grow in importance and lead to novel disease control strategies.

Essential Roles of Cupredoxin Family Proteins in Soybean Cyst Nematode Resistance

Soybean cyst nematode (SCN, Heterodera glycines), one of the most devastating soybean pathogens, causes a significant yield loss in soybean production. One of the most effective ways to manage SCN is to grow resistant cultivars. Therefore, comparative study using resistant and susceptible soybean cultivars provides a powerful tool to identify new genes involved in soybean SCN resistance. In the present study, a transcriptome analysis was carried out using both the resistant (PI88788) and susceptible (Williams 82) soybean cultivars to characterize the responses to nematode infection. Various defense-related genes and different pathways involved in nematode resistance were recognized as being highly expressed in resistant cultivar. Promoter-GUS analysis was conducted to monitor the spatial expression pattern of the genes highly induced by nematode infection. Two nematode-inducible promoters for Glyma.05g147000 (encoding caffeoyl-CoA O-methyltransferase) and Glyma.06g036700 (encoding cupredoxin superfamily protein) were characterized, and the promoters could efficiently drive the expression of known nematode resistance genes (α-SNAPRhg1HC or GmSHMT) to affect soybean SCN resistance. Interestingly, expression of the cupredoxin family genes was upregulated not only by SCN, but also by jasmonic acid treatment. DNA sequence analysis identified that a conserved motif (GGTGCATG) with high similarity to SCNbox1 and GC-rich element is enriched in their promoter regions, suggesting its potential to serve as a nematode-responsive regulatory element. Overexpression of Glyma.06g036700 significantly enhanced soybean resistance to cyst nematode. Overall, our findings not only highlight the essential role of cupredoxin family genes in SCN resistance, but also offer potential functional tools to develop nematode resistance in crops.

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.

Kinase-dead mutation: A novel strategy for improving soybean resistance to soybean cyst nematode Heterodera glycines

Protein kinases phosphorylate proteins for functional changes and are involved in nearly all cellular processes, thereby regulating almost all aspects of plant growth and development, and responses to biotic and abiotic stresses. We generated two independent co-expression networks of soybean genes using control and stress response gene expression data and identified 392 differentially highly interconnected kinase hub genes among the two networks. Of these 392 kinases, 90 genes were identified as "syncytium highly connected hubs", potentially essential for activating kinase signalling pathways in the nematode feeding site. Overexpression of wild-type coding sequences of five syncytium highly connected kinase hub genes using transgenic soybean hairy roots enhanced plant susceptibility to soybean cyst nematode (SCN; Heterodera glycines) Hg Type 0 (race 3). In contrast, overexpression of kinase-dead variants of these five syncytium kinase hub genes significantly enhanced soybean resistance to SCN. Additionally, three of the five tested kinase hub genes enhanced soybean resistance to SCN Hg Type 1.2.5.7 (race 2), highlighting the potential of the kinase-dead approach to generate effective and durable resistance against a wide range of SCN Hg types. Subcellular localization analysis revealed that kinase-dead mutations do not alter protein cellular localization, confirming the structure-function of the kinase-inactive variants in producing loss-of-function phenotypes causing significant decrease in nematode susceptibility. Because many protein kinases are highly conserved and are involved in plant responses to various biotic and abiotic stresses, our approach of identifying kinase hub genes and their inactivation using kinase-dead mutation could be translated for biotic and abiotic stress tolerance.

Evaluation of Scopoletin from Penicillium janthinellum Snef1650 for the Control of Heterodera glycines in Soybean

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is responsible for causing a major soybean disease globally. The fungal strain Penicillium janthinellum Snef1650 was evaluated against H. glycines. However, the effective determinants of the P. janthinellum strain are unknown. By performing pot experiments, a functioning compound was isolated from P. janthinellum Snef1650 through organic solvent extraction, semi-preparative HPLC, Sephadex LH-20 column chromatography, and silica gel column chromatography, and the isolated compound was identified to be scopoletin through 1H NMR, 13C NMR, and HPLC–MS. The pot experiments indicated that the treatment of soybean seeds with scopoletin drastically reduced the SCN population. The field experiments performed in 2017 and 2018 revealed that scopoletin decreased over 43.7% juveniles in the roots and over 61.55% cysts in the soil. Scopoletin treatment also promoted soybean growth and improved its yield, with an increase in plot yield by >5.33%. Scopoletin obtained from P. janthinellum Snef1650 could be used as an anti-H. glycines biocontrol agent.

Novel resistance strategies to soybean cyst nematode (SCN) in wild soybean

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is the most damaging soybean pest worldwide and management of SCN remains challenging. The current SCN resistant soybean cultivars, mainly developed from the cultivated soybean gene pool, are losing resistance due to SCN race shifts. The domestication process and modern breeding practices of soybean cultivars often involve strong selection for desired agronomic traits, and thus, decreased genetic variation in modern cultivars, which consequently resulted in limited sources of SCN resistance. Wild soybean (Glycine soja) is the wild ancestor of cultivated soybean (Glycine max) and it's gene pool is indisputably more diverse than G. max. Our aim is to identify novel resistant genetic resources from wild soybean for the development of new SCN resistant cultivars. In this study, resistance response to HG type 2.5.7 (race 5) of SCN was investigated in a newly identified SCN resistant ecotype, NRS100. To understand the resistance mechanism in this ecotype, we compared RNA seq-based transcriptomes of NRS100 with two SCN-susceptible accessions of G. soja and G. max, as well as an extensively studied SCN resistant cultivar, Peking, under both control and nematode J2-treated conditions. The proposed mechanisms of resistance in NRS100 includes the suppression of the jasmonic acid (JA) signaling pathway in order to allow for salicylic acid (SA) signaling-activated resistance response and polyamine synthesis to promote structural integrity of root cell walls. Our study identifies a set of novel candidate genes and associated pathways involved in SCN resistance and the finding provides insight into the mechanism of SCN resistance in wild soybean, advancing the understanding of resistance and the use of wild soybean-sourced resistance for soybean improvement.

Modulation of (Homo)Glutathione Metabolism and H2O2 Accumulation during Soybean Cyst Nematode Infections in Susceptible and Resistant Soybean Cultivars

In plant immune responses, reactive oxygen species (ROS) act as signaling molecules that activate defense pathways against pathogens, especially following resistance (R) gene-mediated pathogen recognition. Glutathione (GSH), an antioxidant and redox regulator, participates in the removal of hydrogen peroxide (H2O2). However, the mechanism of GSH-mediated H2O2 generation in soybeans (Glycine max (L.) Merr.) that are resistant to the soybean cyst nematode (SCN; Heterodera glycines Ichinohe) remains unclear. To elucidate this underlying relationship, the feeding of race 3 of H. glycines with resistant cultivars, Peking and PI88788, was compared with that on a susceptible soybean cultivar, Williams 82. After 5, 10, and 15 days of SCN infection, we quantified γ-glutamylcysteine (γ-EC) and (homo)glutathione ((h)GSH), and a gene expression analysis showed that GSH metabolism in resistant cultivars differed from that in susceptible soybean roots. ROS accumulation was examined both in resistant and susceptible roots upon SCN infection. The time of intense ROS generation was related to the differences of resistance mechanisms in Peking and PI88788. ROS accumulation that was caused by the (h)GSH depletion-arrested nematode development in susceptible Williams 82. These results suggest that (h)GSH metabolism in resistant soybeans plays a key role in the regulation of ROS-generated signals, leading to resistance against nematodes.

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.

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.

The soybean Rhg1 amino acid transporter gene alters glutamate homeostasis and jasmonic acid-induced resistance to soybean cyst nematode

Rhg1 (resistance to Heterodera glycines 1) is an important locus that contributes to resistance against soybean cyst nematode (SCN; Heterodera glycines Ichinohe), which is the most economically damaging disease of soybean worldwide. Simultaneous overexpression of three genes encoding a predicted amino acid transporter, an α-soluble N-ethylmaleimide-sensitive factor attachment protein (α-SNAP) and a predicted wound-induced protein resulted in resistance to SCN provided by this locus. However, the roles of two of these genes (excluding α-SNAP) remain unknown. Here, we report the functional characterization of Glyma.18G022400, a gene at the Rhg1 locus that encodes the predicted amino acid transporter Rhg1-GmAAT. Although the direct role of Rhg1-GmAAT in glutamate transport was not demonstrated, multiple lines of evidence showed that Rhg1-GmAAT impacts glutamic acid tolerance and glutamate transportation in soybean. Transcriptomic and metabolite profiling indicated that overexpression of Rhg1-GmAAT activated the jasmonic acid (JA) pathway. Treatment with a JA biosynthesis inhibitor reduced the resistance provided by the Rhg1-containing PI88788 to SCN, which suggested that the JA pathway might play a role in Rhg1-mediated resistance to SCN. Our results could be helpful for the clarification of the mechanism of resistance to SCN provided by Rhg1 in soybean.

Functions of Flavonoids in Plant⁻Nematode Interactions

Most land plants can become infected by plant parasitic nematodes in the field. Plant parasitic nematodes can be free-living or endoparasitic, and they usually infect plant roots. Most damaging are endoparasites, which form feeding sites inside plant roots that damage the root system and redirect nutrients towards the parasite. This process involves developmental changes to the root in parallel with the induction of defense responses. Plant flavonoids are secondary metabolites that have roles in both root development and plant defense responses against a range of microorganisms. Here, we review our current knowledge of the roles of flavonoids in the interactions between plants and plant parasitic nematodes. Flavonoids are induced during nematode infection in plant roots, and more highly so in resistant compared with susceptible plant cultivars, but many of their functions remain unclear. Flavonoids have been shown to alter feeding site development to some extent, but so far have not been found to be essential for root–parasite interactions. However, they likely contribute to chemotactic attraction or repulsion of nematodes towards or away from roots and might help in the general plant defense against nematodes. Certain flavonoids have also been associated with functions in nematode reproduction, although the mechanism remains unknown. Much remains to be examined in this area, especially under field conditions.

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.

Genome-wide association study uncovers a novel QTL allele of AtS40-3 that affects the sex ratio of cyst nematodes in Arabidopsis

Plant-parasitic cyst nematodes are obligate sedentary parasites that infect the roots of a broad range of host plants. Cyst nematodes are sexually dimorphic, but differentiation into male or female is strongly influenced by interactions with the host environment. Female populations typically predominate under favorable conditions, whereas male populations predominate under adverse conditions. Here, we performed a genome-wide association study (GWAS) in an Arabidopsis diversity panel to identify host loci underlying variation in susceptibility to cyst nematode infection. Three different susceptibility parameters were examined, with the aim of providing insights into the infection process, the number of females and males present in the infected plant, and the female-to-male sex ratio. GWAS results suggested that variation in sex ratio is associated with a novel quantitative trait locus allele on chromosome 4. Subsequent candidate genes and functional analyses revealed that a senescence-associated transcription factor, AtS40-3, and PPR may act in combination to influence nematode sex ratio. A detailed molecular characterization revealed that variation in nematode sex ratio was due to the disturbed common promoter of AtS40-3 and PPR genes. Additionally, single nucleotide polymorphisms in the coding sequence of AtS40-3 might contribute to the natural variation in nematode sex ratio.

The transcriptomic changes of Huipizhi Heidou (Glycine max), a nematode-resistant black soybean during Heterodera glycines race 3 infection

Glycine max (soybean) is an extremely important crop, representing a major source of oil and protein for human beings. Heterodera glycines (soybean cyst nematode, SCN) infection severely reduces soybean production; therefore, protecting soybean from SCN has become an issue for breeders. Black soybean has exhibited a different grade of resistance to SCN. However, the underlying mechanism of Huipizhi Heidou resistance against SCN remains elusive. The Huipizhi Heidou (ZDD2315) and race 3 of Heterodera glycines were chosen to study the mechanism of resistance via examination of transcriptomic changes. After 5, 10, and 15days of SCN infection, whole roots were sampled for RNA extraction, and uninfected samples were simultaneously collected as a control. 740, 1413, and 4925 genes were isolated by padj (p-value adjusted)<0.05 after 5, 10, and 15days of the infection, respectively, and 225 differentially expressed genes were overlapped at all the time points. We found that the differentially expressed genes (DEGs) at 5, 10, and 15days after infection were involved in various biological function categories; in particular, induced genes were enriched in defense response, hormone mediated signaling process, and response to stress. To verify the pathways observed in the GO and KEGG enrichment results, effects of hormonal signaling in cyst-nematode infection were further examined via treatment with IAA (indo-3-acetic acid), salicylic acid (SA), gibberellic acid (GA), jasmonic acid (JA), and ethephon, a precursor of ethylene. The results indicate that five hormones led to a significant reduction of J2 number in the roots of Huipizhi Heidou and Liaodou15, representing SCN-resistant and susceptible lines, respectively. Taken together, our analyses are aimed at understanding the resistance mechanism of Huipizhi Heidou against the SCN race 3 via the dissection of transcriptomic changes upon J2 infection. The data presented here will help further research on the basis of soybean and cyst-nematode interaction.

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.

Anatomical Alterations in Plant Tissues Induced by Plant-Parasitic Nematodes

Plant-parasitic nematodes (PPNs) interact with plants in different ways, for example, through subtle feeding behavior, migrating destructively through infected tissues, or acting as virus-vectors for nepoviruses. They are all obligate biotrophic parasites as they derive their nutrients from living cells which they modify using pharyngeal gland secretions prior to food ingestion. Some of them can also shield themselves against plant defenses to sustain a relatively long lasting interaction while feeding. This paper is centered on cell types or organs that are newly induced in plants during PPN parasitism, including recent approaches to their study based on molecular biology combined with cell biology-histopathology. This issue has already been reviewed extensively for major PPNs (i.e., root-knot or cyst nematodes), but not for other genera (viz. Nacobbus aberrans, Rotylenchulus spp.). PPNs have evolved with plants and this co-evolution process has allowed the induction of new types of plant cells necessary for their parasitism. There are four basic types of feeding cells: (i) non-hypertrophied nurse cells; (ii) single giant cells; (iii) syncytia; and (iv) coenocytes. Variations in the structure of these cells within each group are also present between some genera depending on the nematode species viz. Meloidogyne or Rotylenchulus. This variability of feeding sites may be related in some way to PPN life style (migratory ectoparasites, sedentary ectoparasites, migratory ecto-endoparasites, migratory endoparasites, or sedentary endoparasites). Apart from their co-evolution with plants, the response of plant cells and roots are closely related to feeding behavior, the anatomy of the nematode (mainly stylet size, which could reach different types of cells in the plant), and the secretory fluids produced in the pharyngeal glands. These secretory fluids are injected through the stylet into perforated cells where they modify plant cytoplasm prior to food removal. Some species do not produce specialized feeding sites (viz. Ditylenchus, Subanguina), but may develop a specialized modification of the root system (e.g., unspecialized root galls or a profusion of roots). This review introduces new data on cell types and plant organs stimulated by PPNs using sources varying from traditional histopathology to new holistic methodologies.

Transcriptional responses of wheat and the cereal cyst nematode Heterodera avenae during their early contact stage

Cereal cyst nematode (Heterodera avenae) is attracted to and aggregated around wheat roots to initiate infection, but this interaction between wheat and the nematode is not fully understood. The transcriptional responses of both wheat and H. avenae were examined during their early contact stage by mRNA sequencing analysis; certain numbers of the differentially expressed genes (DEGs) were validated using quantitative real-time PCR. The immobile host wheat root only had 93 DEGs (27 up-regulated and 66 down-regulated), while the mobile plant parasitic nematode reacted much more actively with 879 DEGs (867 up-regulated and 12 down-regulated). Among them, a number of wheat DEGs (mostly down-regulated) were involved in biotic stress pathways, while several putative effector genes were up-regulated in the nematode DEGs. One putative chitinase-like effector gene of H. avenae was able to suppress BAX-triggered programmed cell death in Nicotiana benthamiana. Results of these experiments demonstrated that nematode responded more actively than wheat during the contact stage of parasitism. The parasite's responses mainly involved up-regulation of genes including at least one anti-plant-defence effector gene, whereas the host responses mainly involved down-regulation of certain defence-related genes.

Histopathology combined with transcriptome analyses reveals the mechanism of resistance to Meloidogyne incognita in Cucumis metuliferus

Root-knot nematodes (Meloidogyne spp.) cause serious threat to cucumber production. Cucumis metuliferus, a relative of cucumber, is reported to be resistant to Meloidogyne incognita, yet the underlying resistance mechanism remains unclear. In this study, the response of resistant C. metuliferus accession PI482443 following nematode infection was studied in comparison with susceptible C. sativus cv. Jinlv No.3. Roots of selected Cucumis seedings were analysed using histological and biochemical techniques. Transcriptome changes of the resistance reaction were investigated by RNA-seq. The results showed that penetration and development of the nematode in resistant plants were reduced when compared to susceptible plants. Infection of a resistant genotype with M. incognita resulted in a hypersensitive reaction. The induction of phenylalanine ammonia lyase and peroxidase activities after infection was greater in resistant than susceptible roots. Several of the most relevant genes for phenylpropanoid biosynthesis, plant hormone signal transduction, and the plant-pathogen interaction pathway that are involved in resistance to the nematode were significantly altered. The resistance in C. metuliferus PI482443 to M. incognita was associated with reduced nematode penetration, retardation of nematode development, and hypersensitive necrosis. The expression of genes resulting in the deposition of lignin, toxic compounds synthesis, cell wall reinforcement, suppression of nematode feeding and resistance protein accumulation, and activation of several transcription factors might all contribute to the resistance response to the pest. These results may lead to a better understanding of the resistance mechanism and aid in the identification of potential targets resistant to pests for cucumber improvement.

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.

Comparative Transcriptome Analysis of Resistant and Susceptible Common Bean Genotypes in Response to Soybean Cyst Nematode Infection

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) reproduces on the roots of common bean (Phaseolus vulgaris L.) and can cause reductions in plant growth and seed yield. The molecular changes in common bean roots caused by SCN infection are unknown. Identification of genetic factors associated with SCN resistance could help in development of improved bean varieties with high SCN resistance. Gene expression profiling was conducted on common bean roots infected by SCN HG type 0 using next generation RNA sequencing technology. Two pinto bean genotypes, PI533561 and GTS-900, resistant and susceptible to SCN infection, respectively, were used as RNA sources eight days post inoculation. Total reads generated ranged between ~ 3.2 and 5.7 million per library and were mapped to the common bean reference genome. Approximately 70-90% of filtered RNA-seq reads uniquely mapped to the reference genome. In the inoculated roots of resistant genotype PI533561, a total of 353 genes were differentially expressed with 154 up-regulated genes and 199 down-regulated genes when compared to the transcriptome of non- inoculated roots. On the other hand, 990 genes were differentially expressed in SCN-inoculated roots of susceptible genotype GTS-900 with 406 up-regulated and 584 down-regulated genes when compared to non-inoculated roots. Genes encoding nucleotide-binding site leucine-rich repeat resistance (NLR) proteins, WRKY transcription factors, pathogenesis-related (PR) proteins and heat shock proteins involved in diverse biological processes were differentially expressed in both resistant and susceptible genotypes. Overall, suppression of the photosystem was observed in both the responses. Furthermore, RNA-seq results were validated through quantitative real time PCR. This is the first report describing genes/transcripts involved in SCN-common bean interaction and the results will have important implications for further characterization of SCN resistance genes in common bean.

A review of methods for nematode identification

Nematodes are non-segmented roundworms found in soil, aquatic environment, plants, or animals. Either useful or pathogenic, they greatly influence environmental equilibrium, human and animal health, as well as plant production. Knowledge on their taxonomy and biology are key issues to answer the different challenges associated to these organisms. Nowadays, most of the nematode taxonomy remains unknown or unclear. Several approaches are available for parasite identification, from the traditional morphology-based techniques to the sophisticated high-throughput sequencing technologies. All these techniques have advantages or drawbacks depending on the sample origin and the number of nematodes to be processed. This review proposes an overview of all newly available methods available to identify known and/or unknown nematodes with a specific focus on emerging high-throughput molecular techniques.

Differentially expressed small RNAs in Arabidopsis galls formed by Meloidogyne javanica: a functional role for miR390 and its TAS3-derived tasiRNAs

Root-knot nematodes (RKNs) induce inside the vascular cylinder the giant cells (GCs) embedded in the galls. The distinctive gene repression in early-developing GCs could be facilitated by small RNAs (sRNA) such as miRNAs, and/or epigenetic mechanisms mediated by 24nt-sRNAs, rasiRNAs and 21-22nt-sRNAs. Therefore, the sRNA-population together with the role of the miR390/TAS3/ARFs module were studied during early gall/GC formation. Three sRNA libraries from 3-d-post-inoculation (dpi) galls induced by Meloidogyne javanica in Arabidopsis and three from uninfected root segments were sequenced following Illumina-Solexa technology. pMIR390a::GUS and pTAS3::GUS lines were assayed for nematode-dependent promoter activation. A sensor line indicative of TAS3-derived tasiRNAs binding to the ARF3 sequence (pARF3:ARF3-GUS) together with a tasiRNA-resistant ARF3 line (pARF3:ARF3m-GUS) were used for functional analysis. The sRNA population showed significant differences between galls and controls, with high validation rate and correspondence with their target expression: 21-nt sRNAs corresponding mainly to miRNAs were downregulated, whilst 24-nt-sRNAs from the rasiRNA family were mostly upregulated in galls. The promoters of MIR390a and TAS3, active in galls, and the pARF3:ARF3-GUS line, indicated a role of TAS3-derived-tasiRNAs in galls. The regulatory module miR390/TAS3 is necessary for proper gall formation possibly through auxin-responsive factors, and the abundance of 24-nt sRNAs (mostly rasiRNAs) constitutes a gall hallmark.

Large-scale identification of wheat genes resistant to cereal cyst nematode Heterodera avenae using comparative transcriptomic analysis

Background

Cereal cyst nematode Heterodera avenae, an important soil-borne pathogen in wheat, causes numerous annual yield losses worldwide, and use of resistant cultivars is the best strategy for control. However, target genes are not readily available for breeding resistant cultivars. Therefore, comparative transcriptomic analyses were performed to identify more applicable resistance genes for cultivar breeding.

Methods

The developing nematodes within roots were stained with acid fuchsin solution. Transcriptome assemblies and redundancy filteration were obtained by Trinity, TGI Clustering Tool and BLASTN, respectively. Gene Ontology annotation was yielded by Blast2GO program, and metabolic pathways of transcripts were analyzed by Path_finder. The ROS levels were determined by luminol-chemiluminescence assay. The transcriptional gene expression profiles were obtained by quantitative RT-PCR.

Results

The RNA-sequencing was performed using an incompatible wheat cultivar VP1620 and a compatible control cultivar WEN19 infected with H. avenae at 24 h, 3 d and 8 d. Infection assays showed that VP1620 failed to block penetration of H. avenae but disturbed the transition of developmental stages, leading to a significant reduction in cyst formation. Two types of expression profiles were established to predict candidate resistance genes after developing a novel strategy to generate clean RNA-seq data by removing the transcripts of H. avenae within the raw data before assembly. Using the uncoordinated expression profiles with transcript abundance as a standard, 424 candidate resistance genes were identified, including 302 overlapping genes and 122 VP1620-specific genes. Genes with similar expression patterns were further classified according to the scales of changed transcript abundances, and 182 genes were rescued as supplementary candidate resistance genes. Functional characterizations revealed that diverse defense-related pathways were responsible for wheat resistance against H. avenae. Moreover, phospholipase was involved in many defense-related pathways and localized in the connection position. Furthermore, strong bursts of reactive oxygen species (ROS) within VP1620 roots infected with H. avenae were induced at 24 h and 3 d, and eight ROS-producing genes were significantly upregulated, including three class III peroxidase and five lipoxygenase genes.

Conclusions

Large-scale identification of wheat resistance genes were processed by comparative transcriptomic analysis. Functional characterization showed that phospholipases associated with ROS production played vital roles in early defense responses to H. avenae via involvement in diverse defense-related pathways as a hub switch. This study is the first to investigate the early defense responses of wheat against H. avenae, not only provides applicable candidate resistance genes for breeding novel wheat cultivars, but also enables a better understanding of the defense mechanisms of wheat against H. avenae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2037-8) contains supplementary material, which is available to authorized users.

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.

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.

Two closely related members of Arabidopsis 13-lipoxygenases (13-LOXs), LOX3 and LOX4, reveal distinct functions in response to plant-parasitic nematode infection

The responses of two closely related members of Arabidopsis 13-lipoxygenases (13-LOXs), LOX3 and LOX4, to infection by the sedentary nematodes root-knot nematode (Meloidogyne javanica) and cyst nematode (Heterodera schachtii) were analysed in transgenic Arabidopsis seedlings. The tissue localization of LOX3 and LOX4 gene expression using β-glucuronidase (GUS) reporter gene constructs showed local induction of LOX3 expression when second-stage juveniles reached the vascular bundle and during the early stages of plant-nematode interaction through gall and syncytia formation. Thin sections of nematode-infested knots indicated LOX3 expression in mature giant cells, and high expression in neighbouring cells and those surrounding the female body. LOX4 promoter was also activated by nematode infection, although the GUS signal weakened as infection and disease progressed. Homozygous insertion mutants lacking LOX3 were less susceptible than wild-type plants to root-knot nematode infection, as reflected by a decrease in female counts. Conversely, deficiency in LOX4 function led to a marked increase in females and egg mass number and in the female to male ratio of M. javanica and H. schachtii, respectively. The susceptibility of lox4 mutants was accompanied by increased expression of allene oxide synthase, allene oxide cyclase and ethylene-responsive transcription factor 4, and the accumulation of jasmonic acid, measured in the roots of lox4 mutants. This response was not found in lox3 mutants. Taken together, our results reveal that LOX4 and LOX3 interfere differentially with distinct metabolic and signalling pathways, and that LOX4 plays a major role in controlling plant defence against nematode infection.

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.

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.

Gene expression profiling of resistant and susceptible soybean lines infected with soybean cyst nematode

Soybean cyst nematode (SCN) is the most devastating pathogen of soybean. Information about the molecular basis of soybean-SCN interactions is needed to assist future development of effective management tools against this pathogen. Toward this end, soybean transcript abundance was measured using the Affymetrix Soybean Genome Array in a susceptible and a resistant reaction of soybean to SCN infection. Two genetically related soybean sister lines TN02-226 and TN02-275, which are resistant and susceptible, respectively, to the SCN race 2 infection were utilized in these experiments. Pairwise comparisons followed by false discovery rate analysis indicated that the expression levels of 162 transcripts changed significantly in the resistant line, of which 84 increased while 78 decreased. However, in the susceptible line, 1,694 transcripts changed significantly, of which 674 increased while 1,020 decreased. Comparative analyses of these transcripts indicated that a total of 51 transcripts were in common between resistance and susceptible responses. In this set, 42 transcripts increased in the resistant line, but decreased in the susceptible line. Quantitative real-time reverse-transcription polymerase chain reaction confirmed the results of microarray analysis. Of the transcripts to which a function could be assigned, genes were associated with metabolism, cell wall modification, signal transduction, transcription, and defense. Microarray analyses examining two genetically related soybean lines against the same SCN population provided additional insights into the specific changes in gene expression of a susceptible and a resistant reaction beneficial for identification of genes involved in defense.

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.

Transcriptional analysis of soybean root response to Fusarium virguliforme, the causal agent of sudden death syndrome

Sudden death syndrome (SDS) of soybean can be caused by any of four distinct Fusarium species, with F. virguliforme and F. tucumaniae being the main casual agents in North and South America, respectively. Although the fungal tissue is largely confined to the roots, the fungus releases a toxin that is translocated to leaf tissues, in which it causes interveinal chlorosis and necrosis leading to scorching symptoms and possible defoliation. In this study, we report on an Affymetrix analysis measuring transcript abundances in resistant (PI 567.374) and susceptible (Essex) roots upon infection by F. virguliforme, 5 and 7 days postinoculation. Many of the genes with increased expression were common between resistant and susceptible plants (including genes related to programmed cell death, the phenylpropanoid pathway, defense, signal transduction, and transcription factors), but some genotype-specific expression was noted. Changes in small (sm)RNA levels between inoculated and mock-treated samples were also studied and implicate a role for these molecules in this interaction. In total, 2,467 genes were significantly changing in the experiment, with 1,694 changing in response to the pathogen; 93 smRNA and 42 microRNA that have putative soybean gene targets were identified from infected tissue. Comparing genotypes, 247 genes were uniquely modulating in the resistant host, whereas 378 genes were uniquely modulating in the susceptible host. Comparing locations of differentially expressed genes to known resistant quantitative trait loci as well as identifying smRNA that increased while their putative targets decreased (or vice versa) allowed for the narrowing of candidate SDS defense-associated genes.

Analysis of gene expression in soybean (Glycine max) roots in response to the root knot nematode Meloidogyne incognita using microarrays and KEGG pathways

BACKGROUND

Root-knot nematodes are sedentary endoparasites that can infect more than 3000 plant species. Root-knot nematodes cause an estimated $100 billion annual loss worldwide. For successful establishment of the root-knot nematode in its host plant, it causes dramatic morphological and physiological changes in plant cells. The expression of some plant genes is altered by the nematode as it establishes its feeding site.

RESULTS

We examined the expression of soybean (Glycine max) genes in galls formed in roots by the root-knot nematode, Meloidogyne incognita, 12 days and 10 weeks after infection to understand the effects of infection of roots by M. incognita. Gene expression was monitored using the Affymetrix Soybean GeneChip containing 37,500 G. max probe sets. Gene expression patterns were integrated with biochemical pathways from the Kyoto Encyclopedia of Genes and Genomes using PAICE software. Genes encoding enzymes involved in carbohydrate and cell wall metabolism, cell cycle control and plant defense were altered.

CONCLUSIONS

A number of different soybean genes were identified that were differentially expressed which provided insights into the interaction between M. incognita and soybean and into the formation and maintenance of giant cells. Some of these genes may be candidates for broadening plants resistance to root-knot nematode through over-expression or silencing and require further examination.

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.

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.

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.