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

Root-Knot Nematode Early Infection Suppresses Immune Response and Elicits the Antioxidant System in Tomato

The immune response in plants is regulated by several phytohormones and involves the overexpression of defense genes, including the pathogenesis-related (PR-) genes. The data reported in this paper indicate that nematodes can suppress the immune response by inhibiting the expression of defense genes. Transcripts from nine defense genes were detected by qRT-PCR in the roots of tomato plants at three and seven days post-inoculation (dpi) with living juveniles (J2s) of Meloidogyne incognita (root-knot nematodes, RKNs). All the salicylic acid (SA)-responsive genes tested (PR-1, PR-2, PR-4b, PR-5) were down-regulated in response to nematode infection. On the contrary, the expression of jasmonic acid (JA)-responsive genes, including ACO (encoding the enzyme 1-aminocyclopropane-1-carboxylic acid oxidase, which catalyzes the last step of ethylene (ET) biosynthesis) and JERF3 (Jasmonate Ethylene Response Factor 3), was unaffected by the infection. Conversely, the effect of nematode attack on the activities of the defense enzymes endoglucanase and endochitinase, encoded by PR-2 and PR-3, respectively, changed depending on the tested dpi. At 5 dpi, both enzymes were inhibited in inoculated plants compared to healthy controls. The genes encoding glutathione peroxidase (GPX) and catalase (CAT), both part of the antioxidant plant system, were highly overexpressed. Additionally, the activity of the antioxidant enzymes superoxide dismutase (SOD), CAT, and ascorbate peroxidase (APX) was enhanced in infected roots. Isoelectrofocusing of root extracts revealed novel SOD isoforms in samples from inoculated plants. Furthermore, plants were pre-treated with an array of key compounds, including hormone generators, inhibitors of SA or JA-mediated defense pathways, reactive oxygen species (ROS) scavengers and generators, inhibitors of ROS generation, and compounds that interfere with calcium-mediated metabolism. After treatments, plants were inoculated with RKNs, and nematodes were allowed to complete their life cycle. Factors of plant growth and infection level in treated plants were compared with those from untreated inoculated plants. Generally, compounds that decreased SA and/or ROS levels increased infection severity, while those that reduced JA/ET levels did not affect infection rates. ROS generators induced resistance against the pests. Compounds that silence calcium signaling by preventing its intake augmented infection symptoms. The data shown in this paper indicate that SA-mediated plant immune responses are consistently suppressed during the early stages of nematode infection, and this restriction is associated with the activation of the antioxidant ROS-scavenging system.

GmSABP2-1 encodes methyl salicylate esterase and functions in soybean defense against soybean cyst nematode

KEY MESSAGE

The soybean gene GmSABP2-1 encodes methyl salicylate esterase and its overexpression led to significant reduction in development of pathogenic soybean cyst nematode. Soybean cyst nematode (SCN, Heterodera glycines) is one of the most devastating pests of soybean (Glycine max L. Merr.). In searching for SCN-defense genes, a soybean gene of the methylesterase (MES) family was found to be upregulated in an SCN-resistant soybean line and downregulated in an SCN-susceptible line upon SCN infection. This gene was designated as GmSABP2-1. Here, we report on biochemical and overexpression studies of GmSABP2-1 to examine its possible function in SCN resistance. The protein encoded by GmSABP2-1 is closely related to known methyl salicylate esterases. To determine the biochemical function of GmSABP2-1, a full-length cDNA of GmSABP2-1 was cloned into a protein expression vector and expressed in Escherichia coli. The resulting recombinant GmSABP2-1 was demonstrated to catalyze the demethylation of methyl salicylate. The biochemical properties of GmSABP2-1 were determined. Its apparent Km value was 46.2 ± 2.2 μM for methyl salicylate, comparable to those of the known methyl salicylate esterases. To explore the biological significance of GmSABP2-1 in soybean defense against SCN, we first overexpressed GmSABP2-1 in transgenic hairy roots of an SCN-susceptible soybean line. When infected with SCN, GmSABP2-1-overexpressing hairy roots showed 84.5% reduction in the development of SCN beyond J2 stage. To provide further genetic evidence for the role of GmSABP2-1 in SCN resistance, stable transgenic soybean plants overexpressing GmSABP2-1 were produced. Analysis of the GmSABP2-1-overexpressing lines showed a significant reduction in SCN development compared to non-transgenic plants. In conclusion, we demonstrated that GmSABP2-1 encodes methyl salicylate esterase and functions as a resistance-related gene against SCN.

Transcriptome Analysis of GmPUB20A Overexpressing and RNA-Interferencing Transgenic Hairy Roots Reveals Underlying Negative Role in Soybean Resistance to Cyst Nematode

Ubiquitination genes are key components of plant responses to biotic stress. GmPUB20A, a ubiquitination gene, plays a negative role in soybean resistance to soybean cyst nematode (SCN). In this study, we employed high-throughput sequencing to investigate transcriptional changes in GmPUB20A overexpressing and RNA-interfering transgenic hairy roots. Totally, 7661 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that DEGs were significantly enriched in disease resistance and signal transduction pathways. In addition, silencing Glyma.15G021600 and Glyma.09G284700 by siRNA, the total number of nematodes was decreased by 33.48% and 27.47% than control plants, respectively. Further, GUS activity and reactive oxygen species (ROS) assays revealed that GmPUB20A, Glyma.15G021600, and Glyma.09G284700 respond to SCN parasitism and interfere with the accumulation of ROS in plant roots, respectively. Collectively, our study provides insights into the molecular mechanism of GmPUB20A in soybean resistance to SCN.

Overexpression of soybean trypsin inhibitor genes decreases defoliation by corn earworm (Helicoverpa zea) in soybean (Glycine max) and Arabidopsis thaliana

Trypsin inhibitors (TIs) are widely distributed in plants and are known to play a protective role against herbivores. TIs reduce the biological activity of trypsin, an enzyme involved in the breakdown of many different proteins, by inhibiting the activation and catalytic reactions of proteins. Soybean (Glycine max) contains two major TI classes: Kunitz trypsin inhibitor (KTI) and Bowman-Birk inhibitor (BBI). Both genes encoding TI inactivate trypsin and chymotrypsin enzymes, which are the main digestive enzymes in the gut fluids of Lepidopteran larvae feeding on soybean. In this study, the possible role of soybean TIs in plant defense against insects and nematodes was investigated. A total of six TIs were tested, including three known soybean trypsin inhibitors (KTI1, KTI2 and KTI3) and three genes encoding novel inhibitors identified in soybean (KTI5, KTI7, and BBI5). Their functional role was further examined by overexpression of the individual TI genes in soybean and Arabidopsis. The endogenous expression patterns of these TI genes varied among soybean tissues, including leaf, stem, seed, and root. In vitro enzyme inhibitory assays showed significant increase in trypsin and chymotrypsin inhibitory activities in both transgenic soybean and Arabidopsis. Detached leaf-punch feeding bioassays detected significant reduction in corn earworm (Helicoverpa zea) larval weight when larvae fed on transgenic soybean and Arabidopsis lines, with the greatest reduction observed in KTI7 and BBI5 overexpressing lines. Whole soybean plant greenhouse feeding bioassays with H. zea on KTI7 and BBI5 overexpressing lines resulted in significantly reduced leaf defoliation compared to non-transgenic plants. However, bioassays of KTI7 and BBI5 overexpressing lines with soybean cyst nematode (SCN, Heterodera glycines) showed no differences in SCN female index between transgenic and non-transgenic control plants. There were no significant differences in growth and productivity between transgenic and non-transgenic plants grown in the absence of herbivores to full maturity under greenhouse conditions. The present study provides further insight into the potential applications of TI genes for insect resistance improvement in plants.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Weighted gene co-expression network analysis identifies genes related to HG Type 0 resistance and verification of hub gene GmHg1

INTRODUCTION

The soybean cyst nematode (SCN) is a major disease in soybean production thatseriously affects soybean yield. At present, there are no studies on weighted geneco-expression network analysis (WGCNA) related to SCN resistance.

METHODS

Here, transcriptome data from 36 soybean roots under SCN HG Type 0 (race 3) stresswere used in WGCNA to identify significant modules.

RESULTS AND DISCUSSION

A total of 10,000 differentially expressed genes and 21 modules were identified, of which the module most related to SCN was turquoise. In addition, the hub gene GmHg1 with high connectivity was selected, and its function was verified. GmHg1 encodes serine/threonine protein kinase (PK), and the expression of GmHg1 in SCN-resistant cultivars ('Dongnong L-204') and SCN-susceptible cultivars ('Heinong 37') increased significantly after HG Type 0 stress. Soybean plants transformed with GmHg1-OX had significantly increased SCN resistance. In contrast, the GmHg1-RNAi transgenic soybean plants significantly reduced SCN resistance. In transgenic materials, the expression patterns of 11 genes with the same expression trend as the GmHg1 gene in the 'turquoise module' were analyzed. Analysis showed that 11genes were co-expressed with GmHg1, which may be involved in the process of soybean resistance to SCN. Our work provides a new direction for studying the Molecular mechanism of soybean resistance to SCN.

Identification of QTL, QTL-by-environment interactions, and their candidate genes for resistance HG Type 0 and HG Type 1.2.3.5.7 in soybean using 3VmrMLM

Introduction

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is an important disease affecting soybean yield in the world. Potential SCN-related QTLs and QTL-by-environment interactions (QEIs) have been used in SCN-resistant breeding.

Methods

In this study, a compressed variance component mixed model, 3VmrMLM, in genome-wide association studies was used to detect QTLs and QEIs for resistance to SCN HG Type 0 and HG Type 1.2.3.5.7 in 156 different soybean cultivars materials.

Results and discussion

The results showed that 53 QTLs were detected in single environment analysis; 36 QTLs and 9 QEIs were detected in multi-environment analysis. Based on the statistical screening of the obtained QTLs, we obtained 10 novel QTLs and one QEI which were different from the previous studies. Based on previous studies, we identified 101 known genes around the significant/suggested QTLs and QEIs. Furthermore, used the transcriptome data of SCN-resistant (Dongnong L-10) and SCN-susceptible (Suinong 14) cultivars, 10 candidate genes related to SCN resistance were identified and verified by Quantitative real time polymerase chain reaction (qRT-PCR) analysis. Haplotype difference analysis showed that Glyma.03G005600 was associated with SCN HG Type 0 and HG Type 1.2.3.5.7 resistance and had a haplotype beneficial to multi-SCN-race resistance. These results provide a new idea for accelerating SCN disease resistance breeding.

Resistance induction and nematicidal activity of certain monoterpenes against tomato root-knot caused by Meloidogyne incognita

This research was performed to evaluate the potential of carvone, cuminaldehyde, cineole, and linalool for the control of root-knot of tomato. The tested control agents were evaluated for their ability to stimulate systemic resistance to Meloidogyne incognita in tomato by monitoring the transcription levels of defense-related genes. Moreover, the ability of the tested agents to induce nematicidal activity concerning second-stage juveniles (J2) hatching and mortality was evaluated. Furthermore, the effect of the tested agents on certain tomato growth and yield parameters was assessed. The tested monoterpenes showed high nematicidal activity against M. incognita concerning J2 hatching inhibition and mortality. Carvone, cuminaldhyde, linalool, and cineole had LC50 values of 123.5, 172.2, 354.9, 466.4, and 952.3 μg/mL, respectively. Carvone was found to be the most efficient hatching inhibitor. The tested monoterpenes showed a high potential against root-knot under greenhouse and field conditions with respect to root-galling, egg masses, and the number of J2. Carvone was the most effective treatment. The growth and yield characters of treated tomato were significantly increased in monoterpenes treatments compared to untreated control. Treated tomato plants showed expression of defense-related genes (PR1 and PAL) 5-8 folds higher than the control. The results also showed that cuminaldhyde, followed by carvone, linalool, and cineole, had the greatest levels of expression in tomato plants. Taken together, the selected monoterpenes could be used as alternatives to control the root-knot of tomato.

Glutathione contributes to plant defence against parasitic cyst nematodes

Cyst nematodes (CNs) are an important group of root-infecting sedentary endoparasites that severely damage many crop plants worldwide. An infective CN juvenile enters the host's roots and migrates towards the vascular cylinder, where it induces the formation of syncytial feeding cells, which nourish the CN throughout its parasitic stages. Here, we examined the role of glutathione (l-γ-glutamyl-l-cysteinyl-glycine) in Arabidopsis thaliana on infection with the CN Heterodera schachtii. Arabidopsis lines with mutations pad2, cad2, or zir1 in the glutamate-cysteine ligase (GSH1) gene, which encodes the first enzyme in the glutathione biosynthetic pathway, displayed enhanced CN susceptibility, but susceptibility was reduced for rax1, another GSH1 allele. Biochemical analysis revealed differentially altered thiol levels in these mutants that was independent of nematode infection. All glutathione-deficient mutants exhibited impaired activation of defence marker genes as well as genes for biosynthesis of the antimicrobial compound camalexin early in infection. Further analysis revealed a link between glutathione-mediated plant resistance to CN infection and the production of camalexin on nematode infection. These results suggest that glutathione levels affect plant resistance to CN by fine-tuning the balance between the cellular redox environment and the production of compounds related to defence against infection.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

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.

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.

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.

Soybean root transcriptome profiling reveals a nonhost resistant response during Heterodera glycines infection

Heterodera glycines (soybean cyst nematode, SCN) is one of the most devastating pathogens of soybean worldwide. The compatible and in compatible interactions between soybean and SCN have well documented. Nevertheless, the molecular mechanism of a nonhost resistant response in soybean against SCN infection remains obscure. Toward this end, a global transcriptional comparison was conducted between susceptible and resistant reactions of soybean roots infected by taking advantage of finding a new pathotype of SCN (SCNT). The soybean cultivar Lee, which exhibits resistant to SCNT and susceptible to HG 1.2.3.4.7 (SCNs) was utilized in the expriments. The results highlighted a nonhost resistant response of soybean. Transcriptome analysis indicated that the number of differentially expressed genes (DEGs) in the resistant interaction (3746) was much larger than that in the susceptible interaction (602). A great number of genes acting as intrinsic component of membrane, integral component of membrane, cell periphery and plasma membrance were remarkably enriched only in the resistant interaction, while the taurine and hypotaurine, phenylpropanoid pathway, plant-pathogen interaction and transcript factors were modulated in both interactions. This is the first study to examine genes expression patterns in a soybean genotype in response to invasion by a virulent and avirulent SCN population at the transcriptional level, which will provide insights into the complicate molecular mechanism of the nonhost resistant interaction.

Genome-wide association study of Gossypium arboreum resistance to reniform nematode

BACKGROUND

Reniform nematode (Rotylenchulus reniformis) has emerged as one of the most destructive root pathogens of upland cotton (Gossypium hirsutum) in the United States. Management of R. reniformis has been hindered by the lack of resistant G. hirsutum cultivars; however, resistance has been frequently identified in germplasm accessions from the G. arboreum collection. To determine the genetic basis of reniform nematode resistance, a genome-wide association study (GWAS) was performed using 246 G. arboreum germplasm accessions that were genotyped with 7220 single nucleotide polymorphic (SNP) sequence markers generated from genotyping-by-sequencing.

RESULTS

Fifteen SNPs representing 12 genomic loci distributed over eight chromosomes showed association with reniform nematode resistance. For 14 SNPs, major alleles were shown to be associated with resistance. From the 15 significantly associated SNPs, 146 genes containing or physically close to these loci were identified as putative reniform nematode resistance candidate genes. These genes are involved in a broad range of biological pathways, including plant innate immunity, transcriptional regulation, and redox reaction that may have a role in the expression of resistance. Eighteen of these genes corresponded to differentially expressed genes identified from G. hirsutum in response to reniform nematode infection.

CONCLUSIONS

The identification of multiple genomic loci associated with reniform nematode resistance would indicate that the G. arboreum collection is a significant resource of novel resistance genes. The significantly associated markers identified from this GWAS can be used for the development of molecular tools for breeding improved reniform nematode resistant upland cotton with resistance introgressed from G. arboreum. Additionally, a greater understanding of the molecular mechanisms of reniform nematode resistance can be determined through genetic structure and functional analyses of candidate genes, which will aid in the pyramiding of multiple resistance genes.

Transcriptomic and metabolomic analyses reveal that bacteria promote plant defense during infection of soybean cyst nematode in soybean

BACKGROUND

Soybean cyst nematode (SCN) is the most devastating pathogen of soybean. Our previous study showed that the plant growth-promoting rhizobacterium Bacillus simplex strain Sneb545 promotes soybean resistance to SCN. Here, we conducted a combined metabolomic and transcriptomic analysis to gain information regarding the biological mechanism of defence enhancement against SCN in Sneb545-treated soybean. To this end, we compared the transcriptome and metabolome of Sneb545-treated and non-treated soybeans under SCN infection.

RESULTS

Transcriptomic analysis showed that 6792 gene transcripts were common in Sneb545-treated and non-treated soybeans. However, Sneb545-treated soybeans showed a higher concentration of various nematicidal metabolites, including 4-vinylphenol, methionine, piperine, and palmitic acid, than non-treated soybeans under SCN infection.

CONCLUSIONS

Overall, our results validated and expanded the existing models regarding the co-regulation of gene expression and metabolites in plants, indicating the advantage of integrated system-oriented analysis.

RNA-seq data comparisons of wild soybean genotypes in response to soybean cyst nematode (Heterodera glycines)

Soybean [Glycine max (L.) Merr.] is an important crop rich in vegetable protein and oil, and is a staple food for human and animals worldwide. However, soybean plants have been challenged by soybean cyst nematode (SCN, Heterodera glycines), one of the most damaging pests found in soybean fields. Applying SCN-resistant cultivars is the most efficient and environmentally friendly strategy to manage SCN. Currently, soybean breeding and further improvement in soybean agriculture are hindered by severely limited genetic diversity in cultivated soybeans. G. soja is a soybean wild progenitor with much higher levels of genetic diversity compared to cultivated soybeans. In this study, transcriptomes of the resistant and susceptible genotypes of the wild soybean, Glycine soja Sieb & Zucc, were sequenced to examine the genetic basis of SCN resistance. Seedling roots were treated with infective second-stage juveniles (J2s) of the soybean cyst nematode (HG type 2.5.7) for 3, 5, 8 days and pooled for library construction and RNA sequencing. The transcriptome sequencing generated approximately 245 million (M) high quality (Q > 30) raw sequence reads (125 bp in length) for twelve libraries. The raw sequence reads were deposited in NCBI sequence read archive (SRA) database, with the accession numbers SRR5227314-25. Further analysis of this data would be helpful to improve our understanding of the molecular mechanisms of soybean-SCN interaction and facilitate the development of diverse SCN resistance cultivars.

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.

Comparative RNA-Seq Analysis Uncovers a Complex Regulatory Network for Soybean Cyst Nematode Resistance in Wild Soybean (Glycine soja)

Soybean cyst nematode (SCN) is the most damaging pest of soybean worldwide. The molecular mechanism of SCN resistance remains largely unknown. We conducted a global RNA-seq comparison between a resistant genotype (S54) and a susceptible genotype (S67) of Glycine soja, the wild progenitor of soybean, to understand its regulatory network in SCN defense. The number of differentially expressed genes (DEGs) in S54 (2,290) was much larger than that in S67 (555). A number of defense-related genes/pathways were significantly induced only in S54, while photosynthesis and several metabolic pathways were affected in both genotypes with SCN infection. These defense-associated DEGs were involved in pathogen recognition, calcium/calmodulin-mediated defense signaling, jasmonic acid (JA)/ethylene (ET) and sialic acid (SA)-involved signaling, the MAPK signaling cascade, and WRKY-involved transcriptional regulation. Our results revealed a comprehensive regulatory network involved in SCN resistance and provided insights into the complex molecular mechanisms of SCN resistance in wild soybean.

Genome-wide identification of soybean microRNA responsive to soybean cyst nematodes infection by deep sequencing

BACKGROUND

The soybean cyst nematode (SCN), Heterodera glycines, is one of the most devastating diseases limiting soybean production worldwide. It is known that small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), play important roles in regulating plant growth and development, defense against pathogens, and responses to environmental changes.

RESULTS

In order to understand the role of soybean miRNAs during SCN infection, we analyzed 24 small RNA libraries including three biological replicates from two soybean cultivars (SCN susceptible KS4607, and SCN HG Type 7 resistant KS4313N) that were grown under SCN-infested and -noninfested soil at two different time points (SCN feeding establishment and egg production). In total, 537 known and 70 putative novel miRNAs in soybean were identified from a total of 0.3 billion reads (average about 13.5 million reads for each sample) with the programs of Bowtie and miRDeep2 mapper. Differential expression analyses were carried out using edgeR to identify miRNAs involved in the soybean-SCN interaction. Comparative analysis of miRNA profiling indicated a total of 60 miRNAs belonging to 25 families that might be specifically related to cultivar responses to SCN. Quantitative RT-PCR validated similar miRNA interaction patterns as sequencing results.

CONCLUSION

These findings suggest that miRNAs are likely to play key roles in soybean response to SCN. The present work could provide a framework for miRNA functional identification and the development of novel approaches for improving soybean SCN resistance in future studies.

An (E,E)-α-farnesene synthase gene of soybean has a role in defence against nematodes and is involved in synthesizing insect-induced volatiles

Plant terpene synthase genes (TPSs) have roles in diverse biological processes. Here, we report the functional characterization of one member of the soybean TPS gene family, which was designated GmAFS. Recombinant GmAFS produced in Escherichia coli catalysed the formation of a sesquiterpene (E,E)-α-farnesene. GmAFS is closely related to (E,E)-α-farnesene synthase gene from apple, both phylogenetically and structurally. GmAFS was further investigated for its biological role in defence against nematodes and insects. Soybean cyst nematode (SCN) is the most important pathogen of soybean. The expression of GmAFS in a SCN-resistant soybean was significantly induced by SCN infection compared with the control, whereas its expression in a SCN-susceptible soybean was not changed by SCN infection. Transgenic hairy roots overexpressing GmAFS under the control of the CaMV 35S promoter were generated in an SCN-susceptible soybean line. The transgenic lines showed significantly higher resistance to SCN, which indicates that GmAFS contributes to the resistance of soybean to SCN. In soybean leaves, the expression of GmAFS was found to be induced by Tetranychus urticae (two-spotted spider mites). Exogenous application of methyl jasmonate to soybean plants also induced the expression of GmAFS in leaves. Using headspace collection combined with gas chromatography-mass spectrometry analysis, soybean plants that were infested with T. urticae were shown to emit a mixture of volatiles with (E,E)-α-farnesene as one of the most abundant constituents. In summary, this study showed that GmAFS has defence roles in both below-ground and above-ground organs of soybean against nematodes and insects, respectively.

Soybean Resistance to the Soybean Cyst Nematode Heterodera glycines: An Update

The soybean cyst nematode (SCN), Heterodera glycines, remains a serious threat to soybean production throughout the world. A lack of genetic diversity in resistant soybean cultivars has led to a widespread shift toward virulence in SCN populations, leaving farmers with few proven options other than nonhost rotation to manage this nematode. Recent advances in our understanding of the genes controlling resistance to the nematode have led to improved molecular markers, which are, in turn, increasing the efficiency and precision of the breeding pipeline. A better understanding of the molecular and biochemical basis of SCN resistance and nematode virulence will provide information useful for the development of a long-term strategic plan for diversification and the deployment of cultivars that protect current sources of natural resistance while identifying new targets for engineering novel resistance.

Proteomic Profiling and the Predicted Interactome of Host Proteins in Compatible and Incompatible Interactions Between Soybean and Fusarium virguliforme

Sudden death syndrome (SDS) is a complex of two diseases of soybean (Glycine max), caused by the soil borne pathogenic fungus Fusarium virguliforme. The root rot and leaf scorch diseases both result in significant yield losses worldwide. Partial SDS resistance has been demonstrated in multiple soybean cultivars. This study aimed to highlight proteomic changes in soybean roots by identifying proteins which are differentially expressed in near isogenic lines (NILs) contrasting at the Rhg1/Rfs2 locus for partial resistance or susceptibility to SDS. Two-dimensional gel electrophoresis resolved approximately 1000 spots on each gel; 12 spots with a significant (P < 0.05) difference in abundance of 1.5-fold or more were picked, trypsin-digested, and analyzed using quadruple time-of-flight tandem mass spectrometry. Several spots contained more than one protein, so that 18 distinct proteins were identified overall. A functional analysis performed to categorize the proteins depicted that the major pathways altered by fungal infection include disease resistance, stress tolerance, and metabolism. This is the first report which identifies proteins whose abundances are altered in response to fungal infection leading to SDS. The results provide valuable information about SDS resistance in soybean plants, and plant partial resistance responses in general. More importantly, several of the identified proteins could be good candidates for the development of SDS-resistant soybean plants.

Transgenic soybean overexpressing GmSAMT1 exhibits resistance to multiple-HG types of soybean cyst nematode Heterodera glycines

Soybean (Glycine max (L.) Merr.) salicylic acid methyl transferase (GmSAMT1) catalyses the conversion of salicylic acid to methyl salicylate. Prior results showed that when GmSAMT1 was overexpressed in transgenic soybean hairy roots, resistance is conferred against soybean cyst nematode (SCN), Heterodera glycines Ichinohe. In this study, we produced transgenic soybean overexpressing GmSAMT1 and characterized their response to various SCN races. Transgenic plants conferred a significant reduction in the development of SCN HG type 1.2.5.7 (race 2), HG type 0 (race 3) and HG type 2.5.7 (race 5). Among transgenic lines, GmSAMT1 expression in roots was positively associated with SCN resistance. In some transgenic lines, there was a significant decrease in salicylic acid titer relative to control plants. No significant seed yield differences were observed between transgenics and control soybean plants grown in one greenhouse with 22 °C day/night temperature, whereas transgenic soybean had higher yield than controls grown a warmer greenhouse (27 °C day/23 °C night) temperature. In a 1-year field experiment in Knoxville, TN, there was no significant difference in seed yield between the transgenic and nontransgenic soybean under conditions with negligible SCN infection. We hypothesize that GmSAMT1 expression affects salicylic acid biosynthesis, which, in turn, attenuates SCN development, without negative consequences to soybean yield or other morphological traits. Thus, we conclude that GmSAMT1 overexpression confers broad resistance to multiple SCN races, which would be potentially applicable to commercial production.

Systemic above- and belowground cross talk: hormone-based responses triggered by Heterodera schachtii and shoot herbivores in Arabidopsis thaliana

Above- and belowground plant parts are simultaneously attacked by different pests and pathogens. The host mediates these interactions and physiologically reacts, e.g. with local and systemic alterations of endogenous hormone levels coupled with coordinated transcriptional changes. This in turn affects attractiveness and susceptibility of the plant to subsequent attackers. Here, the model plant Arabidopsis thaliana is used to study stress hormone-based systemic responses triggered by simultaneous root parasitism by the cyst nematode Heterodera schachtii and shoot herbivory by the thrips Frankliniella occidentalis and the spider mite Tetranychus urticae. First, HPLC/MS and quantitative reverse transcriptase PCR are used to show that nematode parasitism strongly affects stress hormone levels and expression of hormone marker genes in shoots. Previous nematode infection is then demonstrated to affect the behavioural and life history performance of both arthropods. While thrips explicitly avoid nematode-infected plants, spider mites prefer them. In addition, the life history performance of T. urticae is significantly enhanced by nematode infection. Finally, systemic changes triggered by shoot-feeding F. occidentalis but not T. urticae are shown to make the roots more attractive for H. schachtii. This work emphasises the importance of above- and belowground signalling and contributes to a better understanding of plant systemic defence mechanisms against plant-parasitic nematodes.

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.

RNA-Seq Based Identification of Candidate Parasitism Genes of Cereal Cyst Nematode (Heterodera avenae) during Incompatible Infection to Aegilops variabilis

One of the reasons for the progressive yield decline observed in cereals production is the rapid build-up of populations of the cereal cyst nematode (CCN, Heterodera avenae). These nematodes secrete so-call effectors into their host plant to suppress the plant defense responses, alter plant signaling pathways and then induce the formation of syncytium after infection. However, little is known about its molecular mechanism and parasitism during incompatible infection. To gain insight into its repertoire of parasitism genes, we investigated the transcriptome of the early parasitic second-stage (30 hours, 3 days and 9 days post infection) juveniles of the CCN as well as the CCN infected tissue of the host Aegilops variabilis by Illumina sequencing. Among all assembled unigenes, 681 putative genes of parasitic nematode were found, in which 56 putative effectors were identified, including novel pioneer genes and genes corresponding to previously reported effectors. All the 681 CCN unigenes were mapped to 229 GO terms and 200 KEGG pathways, including growth, development and several stimulus-related signaling pathways. Sixteen clusters were involved in the CCN unigene expression atlas at the early stages during infection process, and three of which were significantly gene-enriched. Besides, the protein-protein interaction network analysis revealed 35 node unigenes which may play an important role in the plant-CCN interaction. Moreover, in a comparison of differentially expressed genes between the pre-parasitic juveniles and the early parasitic juveniles, we found that hydrolase activity was up-regulated in pre J2s whereas binding activity was upregulated in infective J2s. RT-qPCR analysis on some selected genes showed detectable expression, indicating possible secretion of the proteins and putative role in infection. This study provided better insights into the incompatible interaction between H. avenae and the host plant Ae. varabilis. Moreover, RNAi targets with potential lethality were screened out and primarily validated, which provide candidates for engineering-based control of cereal cyst nematode in crops breeding.

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.

Genetic characteristics of soybean resistance to HG type 0 and HG type 1.2.3.5.7 of the cyst nematode analyzed by genome-wide association mapping

BACKGROUND

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most fatal pests of soybean (Glycine max (L.) Merr.) worldwide and causes huge loss of soybean yield each year. Multiple sources of resistance are urgently needed for effective management of SCN via the development of resistant cultivars. The aim of the present study was to investigate the genetic architecture of resistance to SCN HG Type 0 (race 3) and HG Type 1.2.3.5.7 (race 4) in landraces and released elite soybean cultivars mostly from China.

RESULTS

A total of 440 diverse soybean landraces and elite cultivars were screened for resistance to SCN HG Type 0 and HG Type 1.2.3.5.7. Exactly 131 new sources of SCN resistance were identified. Lines were genotyped by SNP markers detected by the Specific Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 36,976 SNPs were identified with minor allele frequencies (MAF) > 4% that were present in 97% of all the genotypes. Genome-wide association mapping showed that a total of 19 association signals were significantly related to the resistance for the two HG Types. Of the 19 association signals, eight signals overlapped with reported QTL including Rhg1 and Rhg4 genes. Another eight were located in the linked regions encompassing known QTL. Three QTL were found that were not previously reported. The average value of female index (FI) of soybean accessions with resistant alleles was significantly lower than those with susceptible alleles for each peak SNP. Disease resistance proteins with leucine rich regions, cytochrome P450s, protein kinases, zinc finger domain proteins, RING domain proteins, MYB and WRKY transcription activation families were identified. Such proteins may participate in the resistant reaction to SCN and were frequently found in the tightly linked genomic regions of the peak SNPs.

CONCLUSIONS

GWAS extended understanding of the genetic architecture of SCN resistance in multiple genetic backgrounds. Nineteen association signals were obtained for the resistance to the two Hg Types of SCN. The multiple beneficial alleles from resistant germplasm sources will be useful for the breeding of cultivars with improved resistance to SCN. Analysis of genes near association signals may facilitate the recognition of the causal gene(s) underlying SCN resistances.

Role of stress-related hormones in plant defence during early infection of the cyst nematode Heterodera schachtii in Arabidopsis

Heterodera schachtii, a plant-parasitic cyst nematode, invades host roots and induces a specific syncytial feeding structure, from which it withdraws all required nutrients, causing severe yield losses. The system H. schachtii-Arabidopsis is an excellent research model for investigating plant defence mechanisms. Such responses are suppressed in well-established syncytia, whereas they are induced during early parasitism. However, the mechanisms by which the defence responses are modulated and the role of phytohormones are largely unknown. The aim of this study was to elucidate the role of hormone-based defence responses at the onset of nematode infection. First, concentrations of main phytohormones were quantified and the expression of several hormone-related genes was analysed using quantitative real-time (qRT)-PCR or GeneChip. Further, the effects of individual hormones were evaluated via nematode attraction and infection assays using plants with altered endogenous hormone concentrations. Our results suggest a pivotal and positive role for ethylene during nematode attraction, whereas jasmonic acid triggers early defence responses against H. schachtii. Salicylic acid seems to be a negative regulator during later syncytium and female development. We conclude that nematodes are able to impose specific changes in hormone pools, thus modulating hormone-based defence and signal transduction in strict dependence on their parasitism stage.

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.

Integrated signaling networks in plant responses to sedentary endoparasitic nematodes: a perspective

Sedentary plant endoparasitic nematodes can cause detrimental yield losses in crop plants making the study of detailed cellular, molecular, and whole plant responses to them a subject of importance. In response to invading nematodes and nematode-secreted effectors, plant susceptibility/resistance is mainly determined by the coordination of different signaling pathways including specific plant resistance genes or proteins, plant hormone synthesis and signaling pathways, as well as reactive oxygen signals that are generated in response to nematode attack. Crosstalk between various nematode resistance-related elements can be seen as an integrated signaling network regulated by transcription factors and small RNAs at the transcriptional, posttranscriptional, and/or translational levels. Ultimately, the outcome of this highly controlled signaling network determines the host plant susceptibility/resistance to nematodes.

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.

Computational discovery of soybean promoter cis-regulatory elements for the construction of soybean cyst nematode-inducible synthetic promoters

Computational methods offer great hope but limited accuracy in the prediction of functional cis-regulatory elements; improvements are needed to enable synthetic promoter design. We applied an ensemble strategy for de novo soybean cyst nematode (SCN)-inducible motif discovery among promoters of 18 co-expressed soybean genes that were selected from six reported microarray studies involving a compatible soybean-SCN interaction. A total of 116 overlapping motif regions (OMRs) were discovered bioinformatically that were identified by at least four out of seven bioinformatic tools. Using synthetic promoters, the inducibility of each OMR or motif itself was evaluated by co-localization of gain of function of an orange fluorescent protein reporter and the presence of SCN in transgenic soybean hairy roots. Among 16 OMRs detected from two experimentally confirmed SCN-inducible promoters, 11 OMRs (i.e. 68.75%) were experimentally confirmed to be SCN-inducible, leading to the discovery of 23 core motifs of 5- to 7-bp length, of which 14 are novel in plants. We found that a combination of the three best tools (i.e. SCOPE, W-AlignACE and Weeder) could detect all 23 core motifs. Thus, this strategy is a high-throughput approach for de novo motif discovery in soybean and offers great potential for novel motif discovery and synthetic promoter engineering for any plant and trait in crop biotechnology.

Expression of tomato salicylic acid (SA)-responsive pathogenesis-related genes in Mi-1-mediated and SA-induced resistance to root-knot nematodes

The expression pattern of pathogenesis-related genes PR-1, PR-2 and PR-5, considered as markers for salicylic acid (SA)-dependent systemic acquired resistance (SAR), was examined in the roots and shoots of tomato plants pre-treated with SA and subsequently infected with root-knot nematodes (RKNs) (Meloidogyne incognita). PR-1 was up-regulated in both roots and shoots of SA-treated plants, whereas the expression of PR-5 was enhanced only in roots. The over-expression of PR-1 in the whole plant occurred as soon as 1 day after SA treatment. Up-regulation of the PR-1 gene was considered to be the main marker of SAR elicitation. One day after treatment, plants were inoculated with active juveniles (J2s) of M. incognita. The number of J2s that entered the roots and started to develop was significantly lower in SA-treated than in untreated plants at 5 and 15 days after inoculation. The expression pattern of PR-1, PR-2 and PR-5 was also examined in the roots and shoots of susceptible and Mi-1-carrying resistant tomato plants infected by RKNs. Nematode infection produced a down-regulation of PR genes in both roots and shoots of SA-treated and untreated plants, and in roots of Mi-carrying resistant plants. Moreover, in resistant infected plants, PR gene expression, in particular PR-1 gene expression, was highly induced in shoots. Thus, nematode infection was demonstrated to elicit SAR in shoots of resistant plants. The data presented in this study show that the repression of host defence SA signalling is associated with the successful development of RKNs, and that SA exogenously added as a soil drench is able to trigger a SAR-like response to RKNs in tomato.

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.

Partial resistance to clubroot in Arabidopsis is based on changes in the host primary metabolism and targeted cell division and expansion capacity

To date, studies of the molecular basis of disease resistance mainly focused on qualitative resistance. However, deciphering mechanisms underlying quantitative resistance could lead to insights into the relationship between qualitative and quantitative resistance and guide the utilization of these two types of resistance to produce durably resistant cultivars. A functional genomics approach, using the CATMA whole-genome microarray, was used to detect changes in gene expression associated with partial quantitative resistance in the Arabidopsis thaliana-Plasmodiophora brassicae pathosystem. The time course of transcript abundance during partial clubroot resistance response was monitored at the whole plant level, and direct comparisons between partial resistance and susceptibility responses were made using the same host genotype. An increasingly complex host response was revealed, as was the differential influence of P. brassicae infection on the transcription of Arabidopsis genes according to the isolate used. We observed, at the transcriptomic level, that metabolic diversion by the pathogen was reduced or delayed, classical plant defense responses were induced earlier and/or more strongly, and cell enlargement and proliferation were actively inhibited in the partial quantitative resistance response compared to the susceptible one.

Recent developments of genomic research in soybean

Soybean is an important cash crop with unique and important traits such as the high seed protein and oil contents, and the ability to perform symbiotic nitrogen fixation. A reference genome of cultivated soybeans was established in 2010, followed by whole-genome re-sequencing of wild and cultivated soybean accessions. These efforts revealed unique features of the soybean genome and helped to understand its evolution. Mapping of variations between wild and cultivated soybean genomes were performed. These genomic variations may be related to the process of domestication and human selection. Wild soybean germplasms exhibited high genomic diversity and hence may be an important source of novel genes/alleles. Accumulation of genomic data will help to refine genetic maps and expedite the identification of functional genes. In this review, we summarize the major findings from the whole-genome sequencing projects and discuss the possible impacts on soybean researches and breeding programs. Some emerging areas such as transcriptomic and epigenomic studies will be introduced. In addition, we also tabulated some useful bioinformatics tools that will help the mining of the soybean genomic data.

Transcriptome analysis of resistant soybean roots infected by Meloidogyne javanica

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

De novo assembly and characterization of the root transcriptome of Aegilops variabilis during an interaction with the cereal cyst nematode

BACKGROUND

Aegilops variabilis No.1 is highly resistant to cereal cyst nematode (CCN). However, a lack of genomic information has restricted studies on CCN resistance genes in Ae. variabilis and has limited genetic applications in wheat breeding.

RESULTS

Using RNA-Seq technology, we generated a root transcriptome at a sequencing depth of 4.69 gigabases of Ae. variabilis No. 1 from a pooled RNA sample. The sample contained equal amounts of RNA extracted from CCN-infected and untreated control plants at three time-points. Using the Trinity method, nearly 52,081,238 high-quality trimmed reads were assembled into a non-redundant set of 118,064 unigenes with an average length of 500 bp and an N50 of 599 bp. The total assembly was 59.09 Mb of unique transcriptome sequences with average read-depth coverage of 33.25×. In BLAST searches of our database against public databases, 66.46% (78,467) of the unigenes were annotated with gene descriptions, conserved protein domains, or gene ontology terms. Functional categorization further revealed 7,408 individual unigenes and three pathways related to plant stress resistance.

CONCLUSIONS

We conducted high-resolution transcriptome profiling related to root development and the response to CCN infection in Ae. variabilis No.1. This research facilitates further studies on gene discovery and on the molecular mechanisms related to CCN resistance.