Comparative analysis of pre- and post-parasitic transcriptomes and mining pioneer effectors of Heterodera avenae

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.

Distribution of Bacterial Endosymbionts of the Cardinium Clade in Plant-Parasitic Nematodes

Bacteria of the genus "Candidatus Cardinium" and related organisms composing the Cardinium clade are intracellular endosymbionts frequently occurring in several arthropod groups, freshwater mussels and plant-parasitic nematodes. Phylogenetic analyses based on two gene sequences (16S rRNA and gyrB) showed that the Cardinium clade comprised at least five groups: A, B, C, D and E. In this study, a screening of 142 samples of plant-parasitic nematodes belonging to 93 species from 12 families and two orders using PCR with specific primers and sequencing, revealed bacteria of Cardinium clade in 14 nematode samples belonging to 12 species of cyst nematodes of the family Heteroderidae. Furthermore, in this study, the genome of the Cardinium cHhum from the hop cyst nematode, Heterodera humuli, was also amplified, sequenced and analyzed. The comparisons of the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values for the strain Cardinium cHhum with regard to related organisms with available genomes, combined with the data on 16S rRNA and gyrB gene sequence identities, showed that this strain represents a new candidate species within the genus "Candidatus Paenicardinium". The phylogenetic position of endosymbionts of the Cardinium clade detected in nematode hosts was also compared to known representatives of this clade from other metazoans. Phylogenetic reconstructions based on analysis of 16S rRNA, gyrB, sufB, gloEL, fusA, infB genes and genomes and estimates of genetic distances both indicate that the endosymbiont of the root-lesion nematode Pratylenchus penetrans represented a separate lineage and is designated herein as a new group F. The phylogenetic analysis also confirmed that endosymbionts of ostracods represent the novel group G. Evolutionary relationships of bacterial endosymbionts of the Cardinium clade within invertebrates are presented and discussed.

Description of Rotylenchus zhongshanensis sp. nov. (Tylenchomorpha: Hoplolaimidae) and discovery of its endosymbiont Cardinium

A new bisexual species of Rotylenchus is described and illustrated based on morphological, morphometric and molecular characterizations. Rotylenchus zhongshanensis sp. nov. is characterized by having a conoid lip region complying with the basic pattern for Hoplolaimidae, but with pharyngeal glands slightly overlapping intestine dorsally and cuticle thickened abnormally in female tail terminus. Females have robust stylet (30.1-33.8 μm). The pharyngeal gland has short dorsal (11.2-16.8 μm) overlap on the intestine. The vulva is located at 48.0-56.5% of body length, and phasmids are pore-like, 4-6 annuli posterior to the anus. For males, phasmids are pore-like, 11-17 annuli posterior to cloaca. The spicules are ventrally arcuate (21.0-28.5 μm) with gubernaculum in 5-8 μm length. The rRNA and mitochondrial COI genes were successfully sequenced from the assembled whole-genome sequences of the new species, and were used for reconstructing the phylogenetic relationships of the new species. A new strain of cyto-endosymbiont Cardinium was also discovered from the genome sequences of R. zhongshanensis sp. nov. The 16S rRNA phylogeny analyses revealed that this new bacterial strain is closed to that from cyst and root-lesion nematodes.

A Bursaphelenchus xylophilus Effector, BxSCD3, Suppresses Plant Defense and Contributes to Virulence

Bursaphelenchus xylophilus is the most economically important species of migratory plant-parasitic nematodes (PPNs) and causes severe damage to forestry in China. The successful infection of B. xylophilus relies on the secretion of a repertoire of effector proteins. The effectors, which suppress the host pine immune response, are key to the facilitation of B. xylophilus parasitism. An exhaustive list of candidate effectors of B. xylophilus was predicted, but not all have been identified and characterized. Here, an effector, named BxSCD3, has been implicated in the suppression of host immunity. BxSCD3 could suppress pathogen-associated molecular patterns (PAMPs) PsXEG1- and INF1-triggered cell death when it was secreted into the intracellular space in Nicotiana benthamiana. BxSCD3 was highly up-regulated in the early infection stages of B. xylophilus. BxSCD3 does not affect B. xylophilus reproduction, either at the mycophagous stage or the phytophagous stage, but it contributes to the virulence of B. xylophilus. Moreover, BxSCD3 significantly influenced the relative expression levels of defense-related (PR) genes PtPR-3 and PtPR-6 in Pinus thunbergii in the early infection stage. These results suggest that BxSCD3 is an important toxic factor and plays a key role in the interaction between B. xylophilus and host pine.

On the role of dauer in the adaptation of nematodes to a parasitic lifestyle

Nematodes are presumably the most abundant Metazoa on Earth, and can even be found in some of the most hostile environments of our planet. Various types of hypobiosis evolved to adapt their life cycles to such harsh environmental conditions. The five most distal major clades of the phylum Nematoda (Clades 8-12), formerly referred to as the Secernentea, contain many economically relevant parasitic nematodes. In this group, a special type of hypobiosis, dauer, has evolved. The dauer signalling pathway, which culminates in the biosynthesis of dafachronic acid (DA), is intensively studied in the free-living nematode Caenorhabditis elegans, and it has been hypothesized that the dauer stage may have been a prerequisite for the evolution of a wide range of parasitic lifestyles among other nematode species. Biosynthesis of DA is not specific for hypobiosis, but if it results in exit of the hypobiotic state, it is one of the main criteria to define certain behaviour as dauer. Within Clades 9 and 10, the involvement of DA has been validated experimentally, and dauer is therefore generally accepted to occur in those clades. However, for other clades, such as Clade 12, this has hardly been explored. In this review, we provide clarity on the nomenclature associated with hypobiosis and dauer across different nematological subfields. We discuss evidence for dauer-like stages in Clades 8 to 12 and support this with a meta-analysis of available genomic data. Furthermore, we discuss indications for a simplified dauer signalling pathway in parasitic nematodes. Finally, we zoom in on the host cues that induce exit from the hypobiotic stage and introduce two hypotheses on how these signals might feed into the dauer signalling pathway for plant-parasitic nematodes. With this work, we contribute to the deeper understanding of the molecular mechanisms underlying hypobiosis in parasitic nematodes. Based on this, novel strategies for the control of parasitic nematodes can be developed.

Comparative transcriptomic analysis of candidate effectors to explore the infection and survival strategy of Bursaphelenchus xylophilus during different interaction stages with pine trees

BACKGROUND

The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages.

RESULTS

In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana.

CONCLUSIONS

The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors.

Expression interference of a number of Heterodera avenae conserved genes perturbs nematode parasitic success in Triticum aestivum

The cereal cyst nematode, Heterodera avenae is distributed worldwide and causes substantial damage in bread wheat, Triticum aestivum. This nematode is extremely difficult to manage because of its prolonged persistence as unhatched eggs encased in cysts. Due to its sustainable and target-specific nature, RNA interference (RNAi)-based strategy has gained unprecedented importance for pest control. To date, RNAi strategy has not been exploited to manage H. avenae in wheat. In the present study, 40 H. avenae target genes with different molecular function were rationally selected for in vitro soaking analysis in order to assess their susceptibility to RNAi. In contrast to target-specific downregulation of 18 genes, 7 genes were upregulated and 15 genes showed unaltered expression (although combinatorial soaking showed some of these genes are RNAi susceptible), suggesting that a few of the target genes were refractory or recalcitrant to RNAi. However, RNAi of 37 of these genes negatively altered nematode behavior in terms of reduced penetration, development and reproduction in wheat. Subsequently, wheat plants were transformed with seven H. avenae target genes (that showed greatest abrogation of nematode parasitic success) for host-induced gene silencing (HIGS) analysis. Transformed plants were molecularly characterized by PCR, RT-qPCR and Southern hybridization. Production of target gene-specific double- and single-stranded RNA (dsRNA/siRNA) was detected in transformed plants. Transgenic expression of galectin, cathepsin L, vap1, serpin, flp12, RanBPM and chitinase genes conferred 33.24-72.4 % reduction in H. avenae multiplication in T₁ events with single copy ones exhibiting greatest reduction. A similar degree of resistance observed in T₂ plants indicated the consistent HIGS effect in the subsequent generations. Intriguingly, cysts isolated from RNAi plants were of smaller size with translucent cuticle compared to normal size, dark brown control cysts, suggesting H. avenae developmental retardation due to HIGS. Our study reinforces the potential of HIGS to manage nematode problems in crop plant.

Plant Proteins and Processes Targeted by Parasitic Nematode Effectors

Sedentary endoparasitic nematodes, such as root-knot nematodes (RKN; Meloidogyne spp.) and cyst nematodes (CN; Heterodera spp. and Globodera spp.) cause considerable damage to agricultural crops. RKN and CN spend most of their life cycle in plant roots, in which they induce the formation of multinucleate hypertrophied feeding cells, called "giant cells" and "syncytia," respectively. The giant cells result from nuclear divisions of vascular cells without cytokinesis. They are surrounded by small dividing cells and they form a new organ within the root known as a root knot or gall. CN infection leads to the fusion of several root cells into a unique syncytium. These dramatically modified host cells act as metabolic sinks from which the nematode withdraws nutrients throughout its life, and they are thus essential for nematode development. Both RKN and CN secrete effector proteins that are synthesized in the oesophageal glands and delivered to the appropriate cell in the host plant via a syringe-like stylet, triggering the ontogenesis of the feeding structures. Within the plant cell or in the apoplast, effectors associate with specific host proteins, enabling them to hijack important processes for cell morphogenesis and physiology or immunity. Here, we review recent findings on the identification and functional characterization of plant targets of RKN and CN effectors. A better understanding of the molecular determinants of these biotrophic relationships would enable us to improve the yields of crops infected with parasitic nematodes and to expand our comprehension of root development.

Heterodera avenae GLAND5 Effector Interacts With Pyruvate Dehydrogenase Subunit of Plant to Promote Nematode Parasitism

Heterodera avenae mainly infects cereal crops and causes severe economic losses. Many studies have shown that parasitic nematodes can secrete effector proteins to suppress plant immune responses and then promote parasitism. In this study, we showed that HaGland5, a novel effector of H. avenae, was exclusively expressed in dorsal esophageal gland cell of nematode, and up-regulated in the early parasitic stage. Transgenic Arabidopsis thaliana lines expressing HaGland5 were significantly more susceptible to H. schachtii than wild-type control plants. Conversely, silencing of HaGland5 through barley stripe mosaic virus-medicated host-induced gene silencing technique substantially reduced the infection of H. avenae in wheat. Moreover, HaGland5 could suppress the plant defense responses, including the repression of plant defense-related genes, reducing deposition of cell wall callose and the burst of reactive oxygen species. Mass spectrometry, co-immunoprecipitation, and firefly luciferase complementation imaging assays confirmed that HaGland5 interacted specifically with Arabidopsis pyruvate dehydrogenase subunit (AtEMB3003).

RNAi-Based Biocontrol of Wheat Nematodes Using Natural Poly-Component Biostimulants

With the growing global demands on sustainable food production, one of the biggest challenges to agriculture is associated with crop losses due to parasitic nematodes. While chemical pesticides have been quite successful in crop protection and mitigation of damage from parasites, their potential harm to humans and environment, as well as the emergence of nematode resistance, have necessitated the development of viable alternatives to chemical pesticides. One of the most promising and targeted approaches to biocontrol of parasitic nematodes in crops is that of RNA interference (RNAi). In this study we explore the possibility of using biostimulants obtained from metabolites of soil streptomycetes to protect wheat (Triticum aestivum L.) against the cereal cyst nematode Heterodera avenae by means of inducing RNAi in wheat plants. Theoretical models of uptake of organic compounds by plants, and within-plant RNAi dynamics, have provided us with useful insights regarding the choice of routes for delivery of RNAi-inducing biostimulants into plants. We then conducted in planta experiments with several streptomycete-derived biostimulants, which have demonstrated the efficiency of these biostimulants at improving plant growth and development, as well as in providing resistance against the cereal cyst nematode. Using dot blot hybridization we demonstrate that biostimulants trigger a significant increase of the production in plant cells of si/miRNA complementary with plant and nematode mRNA. Wheat germ cell-free experiments show that these si/miRNAs are indeed very effective at silencing the translation of nematode mRNA having complementary sequences, thus reducing the level of nematode infestation and improving plant resistance to nematodes. Thus, we conclude that natural biostimulants produced from metabolites of soil streptomycetes provide an effective tool for biocontrol of wheat nematode.

Two venom allergen-like proteins, HaVAP1 and HaVAP2, are involved in the parasitism of Heterodera avenae

Despite the fact that venom allergen-like proteins (VAPs) have been identified in many animal- and plant-parasitic nematodes, studies on VAPs in Heterodera avenae, which is an important phytonematode, are still in their infancy. Here, we isolated, cloned and characterized two VAPs, named HaVAP1 and HaVAP2, from H. avenae. The two encoded proteins, HaVAP1 and HaVAP2, harbour an SCP-like domain each, but share only 38% identity with each other. HaVAP1 and HaVAP2 are expressed in subventral and dorsal oesophageal glands, respectively. HaVAP1 is expressed mainly at the early stages, whereas HaVAP2 accumulates principally at the late stages. Both HaVAP1 and HaVAP2 are secreted when expressed in Nicotiana benthamiana leaves, but HaVAP1 is delivered into chloroplasts, whereas HaVAP2 is translocated to the nucleus without signal peptides. Knocking down HaVAP1 increased the virulence of H. avenae. In contrast, silencing of HaVAP2 hampered the parasitism of H. avenae. Both HaVAP1 and HaVAP2 suppressed the cell death induced by BAX in N. benthamiana leaves. Moreover, HaVAP2 physically interacted with a CYPRO4-like protein (HvCLP) of Hordeum vulgare in the nucleus of the plant. It is reasonable to speculate that the changes in the transcript of HvCLP are associated with HaVAP2 during the parasitism of H. avenae. All results obtained in this study show that both HaVAP1 and HaVAP2 are involved in the parasitism of H. avenae, but they possess different functions, broadening our understanding of the parasitic mechanism of H. avenae.

A Novel G16B09-Like Effector From Heterodera avenae Suppresses Plant Defenses and Promotes Parasitism

Plant parasitic nematodes secrete effectors into host plant tissues to facilitate parasitism. In this study, we identified a G16B09-like effector protein family from the transcriptome of Heterodera avenae, and then verified that most of the members could suppress programmed cell death triggered by BAX in Nicotiana benthamiana. Ha18764, the most homologous to G16B09, was further characterized for its function. Our experimental evidence suggested that Ha18764 was specifically expressed in the dorsal gland and was dramatically upregulated in the J4 stage of nematode development. A Magnaporthe oryzae secretion system in barley showed that the signal peptide of Ha18764 had secretion activity to deliver mCherry into plant cells. Arabidopsis thaliana overexpressing Ha18764 or Hs18764 was more susceptible to Heterodera schachtii. In contrast, BSMV-based host-induced gene silencing (HIGS) targeting Ha18764 attenuated H. avenae parasitism and its reproduction in wheat plants. Transient expression of Ha18764 suppressed PsojNIP, Avr3a/R3a, RBP-1/Gpa2, and MAPK kinases (MKK1 and NPK1Nt)-related cell death in Nicotiana benthamiana. Co-expression assays indicated that Ha18764 also suppressed cell death triggered by four H. avenae putative cell-death-inducing effectors. Moreover, Ha18764 was also shown strong PTI suppression such as reducing the expression of plant defense-related genes, the burst of reactive oxygen species, and the deposition of cell wall callose. Together, our results indicate that Ha18764 promotes parasitism, probably by suppressing plant PTI and ETI signaling in the parasitic stages of H. avenae.

Comparative Genomics of Wolbachia-Cardinium Dual Endosymbiosis in a Plant-Parasitic Nematode

Wolbachia and Cardinium are among the most important and widespread of all endosymbionts, occurring in nematodes and more than half of insect and arachnid species, sometimes as coinfections. These symbionts are of significant interest as potential biocontrol agents due to their abilities to cause major effects on host biology and reproduction through cytoplasmic incompatibility, sex ratio distortion, or obligate mutualism. The ecological and metabolic effects of coinfections are not well understood. This study examined a Wolbachia-Cardinium coinfection in the plant-parasitic nematode (PPN), Pratylenchus penetrans, producing the first detailed study of such a coinfection using fluorescence in situ hybridization (FISH), polymerase chain reaction (PCR), and comparative genomic analysis. Results from FISH and single-nematode PCR showed 123/127 individuals in a focal population carried Cardinium (denoted strain cPpe), and 48% were coinfected with Wolbachia strain wPpe. Both endosymbionts showed dispersed tissue distribution with highest densities in the anterior intestinal walls and gonads. Phylogenomic analyses confirmed an early place of cPpe and long distance from a sister strain in another PPN, Heterodera glycines, supporting a long history of both Cardinium and Wolbachia in PPNs. The genome of cPpe was 1.36 Mbp with 35.8% GC content, 1,131 predicted genes, 41% having no known function, and missing biotin and lipoate synthetic capacity and a plasmid present in other strains, despite having a slightly larger genome compared to other sequenced Cardinium. The larger genome revealed expansions of gene families likely involved in host-cellular interactions. More than 2% of the genes of cPpe and wPpe were identified as candidate horizontally transferred genes, with some of these from eukaryotes, including nematodes. A model of the possible Wolbachia-Cardinium interaction is proposed with possible complementation in function for pathways such as methionine and fatty acid biosynthesis and biotin transport.

Endosymbionts of Plant-Parasitic Nematodes

Some of the most agriculturally important plant-parasitic nematodes (PPNs) harbor endosymbionts. Extensive work in other systems has shown that endosymbionts can have major effects on host virulence and biology. This review highlights the discovery, development, and diversity of PPN endosymbionts, incorporating inferences from genomic data. Cardinium, reported from five PPN hosts to date, is characterized by its presence in the esophageal glands and other tissues, with a discontinuous distribution across populations, and genomic data suggestive of horizontal gene exchange. Xiphinematobacter occurs in at least 27 species of dagger nematode in the ovaries and gut epithelial cells, where genomic data suggest it may serve in nutritional supplementation. Wolbachia, reported in just three PPNs, appears to have an ancient history in the Pratylenchidae and displays broad tissue distribution and genomic features intermediate between parasitic and reproductive groups. Finally, a model is described that integrates these insights to explain patterns of endosymbiont replacement.

Characterization of Putative Effectors from the Cereal Cyst Nematode Heterodera avenae

Few molecular details of effectors of Heterodera avenae parasitism are known. We performed a high-throughput sequencing analysis of the H. avenae transcriptome at five developmental stages. A total of 82,549 unigenes were ultimately obtained, and 747 transcripts showed best hits to genes putatively encoding carbohydrate-active enzymes in plant-parasitic nematodes that play an important role in the invasion process. A total of 1,480 unigenes were homologous to known phytonematode effectors, and 63 putative novel effectors were identified in the H. avenae transcriptomes. Twenty-three unigenes were analyzed by qRT-PCR and confirmed to be highly expressed during at least one developmental stage. For in situ hybridization, 17 of the 22 tested putative effectors were specifically expressed and located in the subventral gland cells, and five putative novel effectors were specifically expressed in the dorsal gland. Furthermore, 115 transcripts were found to have putative lethal RNA interference (RNAi) phenotypes. Three target genes with lethal RNAi phenotypes and two of the four tested putative effectors were associated with a decrease in the number of cysts through in vitro RNAi technology. These transcriptomic data lay a foundation for further studies of interactions of H. avenae with cereal and H. avenae parasitic control.

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.

Cereal Cyst Nematodes: A Complex and Destructive Group of Heterodera Species

Small grain cereals have served as the basis for staple foods, beverages, and animal feed for thousands of years. Wheat, barley, oats, rye, triticale, rice, and others are rich in calories, proteins, carbohydrates, vitamins, and minerals. These cereals supply 20% of the calories consumed by people worldwide and are therefore a primary source of energy for humans and play a vital role in global food and nutrition security. Global production of small grains increased linearly from 1960 to 2005, and then began to decline. Further decline in production is projected to continue through 2050 while global demand for these grains is projected to increase by 1% per annum. Currently, wheat, barley, and oat production exceeds consumption in developed countries, while in developing countries the consumption rate is higher than production. An increasing demand for meat and livestock products is likely to compound the demand for cereals in developing countries. Current production levels and trends will not be sufficient to fulfill the projected global demand generated by increased populations. For wheat, global production will need to be increased by 60% to fulfill the estimated demand in 2050. Until recently, global wheat production increased mostly in response to development of improved cultivars and farming practices and technologies. Production is now limited by biotic and abiotic constraints, including diseases, nematodes, insect pests, weeds, and climate. Among these constraints, plant-parasitic nematodes alone are estimated to reduce production of all world crops by 10%. Cereal cyst nematodes (CCNs) are among the most important nematode pests that limit production of small grain cereals. Heavily invaded young plants are stunted and their lower leaves are often chlorotic, forming pale green patches in the field. Mature plants are also stunted, have a reduced number of tillers, and the roots are shallow and have a "bushy-knotted" appearance. CCNs comprise a number of closely-related species and are found in most regions where cereals are produced.

Large-Scale Identification and Characterization of Heterodera avenae Putative Effectors Suppressing or Inducing Cell Death in Nicotiana benthamiana

Heterodera avenae is one of the most important plant pathogens and causes vast losses in cereal crops. As a sedentary endoparasitic nematode, H. avenae secretes effectors that modify plant defenses and promote its biotrophic infection of its hosts. However, the number of effectors involved in the interaction between H. avenae and host defenses remains unclear. Here, we report the identification of putative effectors in H. avenae that regulate plant defenses on a large scale. Our results showed that 78 of the 95 putative effectors suppressed programmed cell death (PCD) triggered by BAX and that 7 of the putative effectors themselves caused cell death in Nicotiana benthamiana. Among the cell-death-inducing effectors, three were found to be dependent on their specific domains to trigger cell death and to be expressed in esophageal gland cells by in situ hybridization. Ten candidate effectors that suppressed BAX-triggered PCD also suppressed PCD triggered by the elicitor PsojNIP and at least one R-protein/cognate effector pair, suggesting that they are active in suppressing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). Notably, with the exception of isotig16060, these putative effectors could also suppress PCD triggered by cell-death-inducing effectors from H. avenae, indicating that those effectors may cooperate to promote nematode parasitism. Collectively, our results indicate that the majority of the tested effectors of H. avenae may play important roles in suppressing cell death induced by different elicitors in N. benthamiana.