The genome and lifestage-specific transcriptomes of a plant-parasitic nematode and its host reveal susceptibility genes involved in trans-kingdom synthesis of vitamin B5.

Plant-parasitic nematodes are a major threat to crop production in all agricultural systems. The scarcity of classical resistance genes highlights a pressing need to find new ways to develop nematode-resistant germplasm. Here, we sequence and assemble a high-quality phased genome of the model cyst nematode Heterodera schachtii to provide a platform for the first system-wide dual analysis of host and parasite gene expression over time, covering all major parasitism stages. Analysis of the hologenome of the plant-nematode infection site identified metabolic pathways that were incomplete in the parasite but complemented by the host. Using a combination of bioinformatic, genetic, and biochemical approaches, we show that a highly atypical completion of vitamin B5 biosynthesis by the parasitic animal, putatively enabled by a horizontal gene transfer from a bacterium, is required for full pathogenicity. Knockout of either plant-encoded or now nematode-encoded steps in the pathway significantly reduces parasitic success. Our experiments establish a reference for cyst nematodes, further our understanding of the evolution of plant-parasitism by nematodes, and show that congruent differential expression of metabolic pathways in the infection hologenome represents a new way to find nematode susceptibility genes. The approach identifies genome-editing-amenable targets for future development of nematode-resistant crops.

A chromosomal assembly of the soybean cyst nematode genome.

The soybean cyst nematode (Heterodera glycines) is a sedentary plant parasite that exceeds billion USD annually in yield losses. This problem is exacerbated by H. glycines populations overcoming the limited sources of natural resistance in soybean and by the lack of effective and safe alternative treatments. Although there are genetic determinants that render soybeans resistant to nematode genotypes, resistant soybeans are increasingly ineffective because their multiyear usage has selected for virulent H. glycines populations. Successful H. glycines infection relies on the comprehensive re-engineering of soybean root cells into a syncytium, as well as the long-term suppression of host defences to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms that control genomic effector acquisition, diversification, and selection are important insights needed for the development of essential novel control strategies. As a foundation to obtain this understanding, we created a nine-scaffold, 158 Mb pseudomolecule assembly of the H. glycines genome using PacBio, Chicago, and Hi-C sequencing. A Mikado consensus gene prediction produced an annotation of 22,465 genes using short- and long-read expression data. To evaluate assembly and annotation quality, we cross-examined synteny among H. glycines assemblies, and compared BUSCO across related species. To describe the predicted proteins involved in H. glycines' secretory pathway, we contrasted expression between preparasitic and parasitic stages with functional gene information. Here, we present the results from our assembly and annotation of the H. glycines genome and contribute this resource to the scientific community.

Esophageal Gland RNA-Seq Resource of a Virulent and Avirulent Population of the Soybean Cyst Nematode Heterodera glycines.

The soybean cyst nematode Heterodera glycines is the most economically devastating pathogen of soybean in the United States and threatens to become even more damaging through the selection of virulent nematode populations in the field that can overcome natural resistance mechanisms in soybean cultivars. This pathogen, therefore, demands intense transcriptomic/genomic research inquiries into the biology of its parasitic mechanisms. H. glycines delivers effector proteins that are produced in specialized gland cells into the soybean root to enable infection. The study of effector proteins, thus, is particularly promising when exploring novel management options against this pathogen. Here, we announce the availability of a gland cell-specific RNA-seq resource. These data represent an expression snapshot of gland cell activity during early soybean infection of a virulent and an avirulent H. glycines population, providing a unique and highly valuable resource for scientists examining effector biology and nematode virulence.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Large tandem duplications affect gene expression, 3D organization, and plant-pathogen response.

Rapid plant genome evolution is crucial to adapt to environmental changes. Chromosomal rearrangements and gene copy number variation (CNV) are two important tools for genome evolution and sources for the creation of new genes. However, their emergence takes many generations. In this study, we show that in Arabidopsis thaliana, a significant loss of ribosomal RNA (rRNA) genes with a past history of a mutation for the chromatin assembly factor 1 (CAF1) complex causes rapid changes in the genome structure. Using long-read sequencing and microscopic approaches, we have identified up to 15 independent large tandem duplications in direct orientation (TDDOs) ranging from 60 kb to 1.44 Mb. Our data suggest that these TDDOs appeared within a few generations, leading to the duplication of hundreds of genes. By subsequently focusing on a line only containing 20% of rRNA gene copies (20rDNA line), we investigated the impact of TDDOs on 3D genome organization, gene expression, and cytosine methylation. We found that duplicated genes often accumulate more transcripts. Among them, several are involved in plant-pathogen response, which could explain why the 20rDNA line is hyper-resistant to both bacterial and nematode infections. Finally, we show that the TDDOs create gene fusions and/or truncations and discuss their potential implications for the evolution of plant genomes.

SCNBase: a genomics portal for the soybean cyst nematode (Heterodera glycines).

Soybean is an important worldwide crop, and farmers continue to experience significant yield loss due to the soybean cyst nematode (SCN), Heterodera glycines. This soil-borne roundworm parasite is rated the most important pathogen problem in soybean production. The infective nematodes enter into complex interactions with their host plant by inducing the development of specialized plant feeding cells that provide the parasites with nourishment. Addressing the SCN problem will require the development of genomic resources and a global collaboration of scientists to analyze and use these resources. SCNBase.org was designed as a collaborative hub for the SCN genome. All data and analyses are downloadable and can be analyzed with three integrated genomic tools: JBrowse, Feature Search and BLAST. At the time of this writing, a number of genomic and transcriptomic data sets are already available, with 43 JBrowse tracks and 21 category pages describing SCN genomic analyses on gene predictions, transcriptome and read alignments, effector-like genes, expansion and contraction of genomic repeats, orthology and synteny with related nematode species, Single Nucleotide Polymorphism (SNPs) from 15 SCN populations and novel splice sites. Standard functional gene annotations were supplemented with orthologous gene annotations using a comparison to nine related plant-parasitic nematodes, thereby enabling functional annotations for 85% of genes. These annotations led to a greater grasp on the SCN effectorome, which include over 3324 putative effector genes. By designing SCNBase as a hub, future research findings and genomic resources can easily be uploaded and made available for use by others with minimal needs for further curation. By providing these resources to nematode research community, scientists will be empowered to develop novel, more effective SCN management tools.

The genome of the soybean cyst nematode (Heterodera glycines) reveals complex patterns of duplications involved in the evolution of parasitism genes.

BACKGROUND: Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock. RESULTS: Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons. CONCLUSIONS: This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

The cyst nematode effector protein 10A07 targets and recruits host posttranslational machinery to mediate its nuclear trafficking and to promote parasitism in Arabidopsis.

Plant-parasitic cyst nematodes synthesize and secrete effector proteins that are essential for parasitism. One such protein is the 10A07 effector from the sugar beet cyst nematode, Heterodera schachtii, which is exclusively expressed in the nematode dorsal gland cell during all nematode parasitic stages. Overexpression of H. schachtii 10A07 in Arabidopsis thaliana produced a hypersusceptible phenotype in response to H. schachtii infection along with developmental changes reminiscent of auxin effects. The 10A07 protein physically associates with a plant kinase and the IAA16 transcription factor in the cytoplasm and nucleus, respectively. The interacting plant kinase (IPK) phosphorylates 10A07 at Ser-144 and Ser-231 and mediates its trafficking from the cytoplasm to the nucleus. Translocation to the nucleus is phosphorylation dependent since substitution of Ser-144 and Ser-231 by alanine resulted in exclusive cytoplasmic accumulation of 10A07. IPK and IAA16 are highly upregulated in the nematode-induced syncytium (feeding cells), and deliberate manipulations of their expression significantly alter plant susceptibility to H. schachtii in an additive fashion. An inactive variant of IPK functioned antagonistically to the wild-type IPK and caused a dominant-negative phenotype of reduced plant susceptibility. Thus, exploitation of host processes to the advantage of the parasites is one mechanism by which cyst nematodes promote parasitism of host plants.

Members of the Meloidogyne avirulence protein family contain multiple plant ligand-like motifs.

Sedentary plant-parasitic nematodes engage in complex interactions with their host plants by secreting effector proteins. Some effectors of both root-knot nematodes (Meloidogyne spp.) and cyst nematodes (Heterodera and Globodera spp.) mimic plant ligand proteins. Most prominently, cyst nematodes secrete effectors that mimic plant CLAVATA3/ESR-related (CLE) ligand proteins. However, only cyst nematodes have been shown to secrete such effectors and to utilize CLE ligand mimicry in their interactions with host plants. Here, we document the presence of ligand-like motifs in bona fide root-knot nematode effectors that are most similar to CLE peptides from plants and cyst nematodes. We have identified multiple tandem CLE-like motifs conserved within the previously identified Meloidogyne avirulence protein (MAP) family that are secreted from root-knot nematodes and have been shown to function in planta. By searching all 12 MAP family members from multiple Meloidogyne spp., we identified 43 repetitive CLE-like motifs composing 14 unique variants. At least one CLE-like motif was conserved in each MAP family member. Furthermore, we documented the presence of other conserved sequences that resemble the variable domains described in Heterodera and Globodera CLE effectors. These findings document that root-knot nematodes appear to use CLE ligand mimicry and point toward a common host node targeted by two evolutionarily diverse groups of nematodes. As a consequence, it is likely that CLE signaling pathways are important in other phytonematode pathosystems as well.

Mining novel effector proteins from the esophageal gland cells of Meloidogyne incognita.

Meloidogyne incognita is one of the most economically damaging plant pathogens in agriculture and horticulture. Identifying and characterizing the effector proteins which M. incognita secretes into its host plants during infection is an important step toward finding new ways to manage this pest. In this study, we have identified the cDNAs for 18 putative effectors (i.e., proteins that have the potential to facilitate M. incognita parasitism of host plants). These putative effectors are secretory proteins that do not contain transmembrane domains and whose genes are specifically expressed in the secretory gland cells of the nematode, indicating that they are likely secreted from the nematode through its stylet. We have determined that, in the plant cells, these putative effectors are likely to localize to the cytoplasm. Furthermore, the transcripts of many of these novel effectors are specifically upregulated during different stages of the nematode's life cycle, indicating that they function at specific stages during M. incognita parasitism. The predicted proteins showed little to no homology to known proteins from free-living nematode species, suggesting that they evolved recently to support the parasitic lifestyle. On the other hand, several of the effectors are part of gene families within the M. incognita genome as well as that of M. hapla, which points to an important role that these putative effectors are playing in both parasites. With the discovery of these putative effectors, we have increased our knowledge of the effector repertoire utilized by root-knot nematodes to infect, feed on, and reproduce on their host plants. Future studies investigating the roles that these proteins play in planta will help mitigate the effects of this damaging pest.

Isolation of whole esophageal gland cells from plant-parasitic nematodes for transcriptome analyses and effector identification.

Esophageal glands of plant-parasitic nematodes are highly specialized cells whose gene expression products include secreted effector proteins, which govern nematode parasitism of host plants. Therefore, elucidating the transcriptomes of esophageal glands with the goal of identifying nematode effectors is a promising avenue to understanding nematode parasitism and its evolutionary origins as well as to devising nematode control strategies. We have developed a method to separate and isolate individual esophageal gland cells from multiple species of plant-parasitic nematodes while preserving RNA quality. We have used such isolated gland cells for transcriptome analysis via high-throughput DNA sequencing. This method relies on the differential histochemical staining of the gland cells after homogenization of phytonematode tissues. Total RNA was extracted from whole gland cells isolated from eight different plant-parasitic nematode species. To validate this approach, the isolated RNA from three plant-parasitic nematode species-Globodera rostochiensis, Pratylenchus penetrans, and Radopholus similis-was amplified, gel purified, and used for 454 sequencing. We obtained 456,801 total reads with an average read length of 409 bp. Sequence analyses revealed the presence of homologs of previously known nematode effectors in these libraries, thus validating our approach. These data provide compelling evidence that this technical advance can be used to relatively easily and expediently discover effector repertoires of plant-parasitic nematodes.

The Arabidopsis microRNA396-GRF1/GRF3 regulatory module acts as a developmental regulator in the reprogramming of root cells during cyst nematode infection.

The syncytium is a unique plant root organ whose differentiation is induced by plant-parasitic cyst nematodes to create a source of nourishment. Syncytium formation involves the redifferentiation and fusion of hundreds of root cells. The underlying regulatory networks that control this unique change of plant cell fate are not understood. Here, we report that a strong down-regulation of Arabidopsis (Arabidopsis thaliana) microRNA396 (miR396) in cells giving rise to the syncytium coincides with the initiation of the syncytial induction/formation phase and that specific miR396 up-regulation in the developed syncytium marks the beginning of the maintenance phase, when no new cells are incorporated into the syncytium. In addition, our results show that miR396 in fact has a role in the transition from one phase to the other. Expression modulations of miR396 and its Growth-Regulating Factor (GRF) target genes resulted in reduced syncytium size and arrested nematode development. Furthermore, genome-wide expression profiling revealed that the miR396-GRF regulatory system can alter the expression of 44% of the more than 7,000 genes reported to change expression in the Arabidopsis syncytium. Thus, miR396 represents a key regulator for the reprogramming of root cells. As such, this regulatory unit represents a powerful molecular target for the parasitic animal to modulate plant cells and force them into novel developmental pathways.

Arabidopsis spermidine synthase is targeted by an effector protein of the cyst nematode Heterodera schachtii.

Cyst nematodes are sedentary plant parasites that cause dramatic cellular changes in the plant root to form feeding cells, so-called syncytia. 10A06 is a cyst nematode secretory protein that is most likely secreted as an effector into the developing syncytia during early plant parasitism. A homolog of the uncharacterized soybean cyst nematode (Heterodera glycines), 10A06 gene was cloned from the sugar beet cyst nematode (Heterodera schachtii), which is able to infect Arabidopsis (Arabidopsis thaliana). Constitutive expression of 10A06 in Arabidopsis affected plant morphology and increased susceptibility to H. schachtii as well as to other plant pathogens. Using yeast two-hybrid assays, we identified Spermidine Synthase2 (SPDS2), a key enzyme involved in polyamine biosynthesis, as a specific 10A06 interactor. In support of this protein-protein interaction, transgenic plants expressing 10A06 exhibited elevated SPDS2 mRNA abundance, significantly higher spermidine content, and increased polyamine oxidase (PAO) activity. Furthermore, the SPDS2 promoter was strongly activated in the nematode-induced syncytia, and transgenic plants overexpressing SPDS2 showed enhanced plant susceptibility to H. schachtii. In addition, in planta expression of 10A06 or SPDS2 increased mRNA abundance of a set of antioxidant genes upon nematode infection. These data lend strong support to a model in which the cyst nematode effector 10A06 exerts its function through the interaction with SPDS2, thereby increasing spermidine content and subsequently PAO activity. Increasing PAO activity results in stimulating the induction of the cellular antioxidant machinery in syncytia. Furthermore, we observed an apparent disruption of salicylic acid defense signaling as a function of 10A06. Most likely, increased antioxidant protection and interruption of salicylic acid signaling are key aspects of 10A06 function in addition to other physiological and morphological changes caused by altered polyamines, which are potent plant signaling molecules.

Effective and specific in planta RNAi in cyst nematodes: expression interference of four parasitism genes reduces parasitic success.

Cyst nematodes are highly evolved sedentary plant endoparasites that use parasitism proteins injected through the stylet into host tissues to successfully parasitize plants. These secretory proteins likely are essential for parasitism as they are involved in a variety of parasitic events leading to the establishment of specialized feeding cells required by the nematode to obtain nourishment. With the advent of RNA interference (RNAi) technology and the demonstration of host-induced gene silencing in parasites, a new strategy to control pests and pathogens has become available, particularly in root-knot nematodes. Plant host-induced silencing of cyst nematode genes so far has had only limited success but similarly should disrupt the parasitic cycle and render the host plant resistant. Additional in planta RNAi data for cyst nematodes are being provided by targeting four parasitism genes through host-induced RNAi gene silencing in transgenic Arabidopsis thaliana, which is a host for the sugar beet cyst nematode Heterodera schachtii. Here it is reported that mRNA abundances of targeted nematode genes were specifically reduced in nematodes feeding on plants expressing corresponding RNAi constructs. Furthermore, this host-induced RNAi of all four nematode parasitism genes led to a reduction in the number of mature nematode females. Although no complete resistance was observed, the reduction of developing females ranged from 23% to 64% in different RNAi lines. These observations demonstrate the relevance of the targeted parasitism genes during the nematode life cycle and, potentially more importantly, suggest that a viable level of resistance in crop plants may be accomplished in the future using this technology against cyst nematodes.

Cellulose binding protein from the parasitic nematode Heterodera schachtii interacts with Arabidopsis pectin methylesterase: cooperative cell wall modification during parasitism.

Plant-parasitic cyst nematodes secrete a complex of cell wall-digesting enzymes, which aid in root penetration and migration. The soybean cyst nematode Heterodera glycines also produces a cellulose binding protein (Hg CBP) secretory protein. To determine the function of CBP, an orthologous cDNA clone (Hs CBP) was isolated from the sugar beet cyst nematode Heterodera schachtii, which is able to infect Arabidopsis thaliana. CBP is expressed only in the early phases of feeding cell formation and not during the migratory phase. Transgenic Arabidopsis expressing Hs CBP developed longer roots and exhibited enhanced susceptibility to H. schachtii. A yeast two-hybrid screen identified Arabidopsis pectin methylesterase protein 3 (PME3) as strongly and specifically interacting with Hs CBP. Transgenic plants overexpressing PME3 also produced longer roots and exhibited increased susceptibility to H. schachtii, while a pme3 knockout mutant showed opposite phenotypes. Moreover, CBP overexpression increases PME3 activity in planta. Localization studies support the mode of action of PME3 as a cell wall-modifying enzyme. Expression of CBP in the pme3 knockout mutant revealed that PME3 is required but not the sole mechanism for CBP overexpression phenotype. These data indicate that CBP directly interacts with PME3 thereby activating and potentially targeting this enzyme to aid cyst nematode parasitism.

Arabidopsis small RNAs and their targets during cyst nematode parasitism.

Plant-parasitic cyst nematodes induce the formation of specialized feeding cells in infected roots, which involves plant developmental processes that have been shown to be influenced by microRNAs (miRNAs) and other small RNAs. This observation provided the foundation to investigate the potential involvement of small RNAs in plant-cyst nematode interactions. First, we examined the susceptibilities of Arabidopsis DICER-like (dcl) and RNA-dependent RNA polymerase (rdr) mutants to the sugar beet cyst nematode Heterodera schachtii. The examined mutants exhibited a trend of decreased susceptibility, suggesting a role of small RNAs mediating gene regulation processes during the plant-nematode interaction. Second, we generated two small RNA libraries from aseptic Arabidopsis roots harvested at 4 and 7 days after infection with surface-sterilized H. schachtii. Sequences of known miRNAs as well as novel small interfering (si)RNAs were identified. Following this discovery, we used real-time reverse-transcriptase polymerase chain reaction to quantify a total of 15 Arabidopsis transcripts that are known targets of six of the different miRNA families found in our study (miR160, miR164, miR167, miR171, miR396, and miR398) in inoculated and noninoculated Arabidopsis roots. Our analyses showed mostly negative correlations between miRNA accumulation and target gene mRNA abundance, suggesting regulatory roles of these miRNAs during parasitism. Also, we identified a total of 125 non-miRNA siRNAs. Some of these siRNAs perfectly complement protein-coding mRNAs or match transposon or retrotransposon sequences in sense or antisense orientations. We further quantified a group of siRNAs in H. schachtii-inoculated roots. The examined siRNAs exhibited distinct expression patterns in infected and noninfected roots, providing additional evidence for the implication of small RNAs in cyst nematode parasitism. These data lay the foundation for detailed analyses of the functions of small RNAs during phytonematode parasitism.

Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita.

Plant-parasitic nematodes are major agricultural pests worldwide and novel approaches to control them are sorely needed. We report the draft genome sequence of the root-knot nematode Meloidogyne incognita, a biotrophic parasite of many crops, including tomato, cotton and coffee. Most of the assembled sequence of this asexually reproducing nematode, totaling 86 Mb, exists in pairs of homologous but divergent segments. This suggests that ancient allelic regions in M. incognita are evolving toward effective haploidy, permitting new mechanisms of adaptation. The number and diversity of plant cell wall-degrading enzymes in M. incognita is unprecedented in any animal for which a genome sequence is available, and may derive from multiple horizontal gene transfers from bacterial sources. Our results provide insights into the adaptations required by metazoans to successfully parasitize immunocompetent plants, and open the way for discovering new antiparasitic strategies.

GmEREBP1 is a transcription factor activating defense genes in soybean and Arabidopsis.

Ethylene-responsive element-binding proteins (EREBPs) are plant-specific transcription factors, many of which have been linked to plant defense responses. Conserved EREBP domains bind to the GCC box, a promoter element found in pathogenesis-related (PR) genes. We previously identified an EREBP gene from soybean (GmEREBP1) whose transcript abundance decreased in soybean cyst-nematode-infected roots of a susceptible cultivar, whereas it increased in abundance in infected roots of a resistant cultivar. Here, we report further characterization of this gene. Transient expression analyses showed that GmEREBP1 is localized to the plant nucleus and functions as a transcriptional activator in soybean leaves. Transgenic soybean plants expressing GmEREBP1 activated the expression of the ethylene (ET)-responsive gene PR2 and the ET- and jasmonic acid (JA)-responsive gene PR3, and the salicylic acid (SA)-responsive gene PR1 but not the SA-responsive PR5. Similarly, transgenic Arabidopsis plants expressing GmEREBP1 showed elevated mRNA abundance of the ET-regulated gene PR3 and the ET- and JA-regulated defense-related gene PDF1.2 but not the ET-regulated GST2, and the SA-regulated gene PR1 but not the SA-regulated PR2 and PR5. Transgenic soybean and Arabidopsis plants inoculated with cyst nematodes did not display a significantly altered susceptibility to nematode infection. These results collectively show that GmEREBP1 functions as a transacting inducer of defense gene expression in both soybean and Arabidopsis and mediates the expression of both ET- and JA- and SA-regulated defense-related genes in these plant species.