Comparative Transcriptome Analysis of the Pinewood Nematode Bursaphelenchus xylophilus Reveals the Molecular Mechanism Underlying Its Defense Response to Host-Derived α-pinene

Engineered pine endophytic Bacillus toyonensis with nematocidal and colonization abilities for pine wilt disease control

Introduction

The pinewood nematode (PWN) is responsible for causing pine wilt disease (PWD), which has led to the significant decline of conifer species in Eurasian forests and has become a globally invasive quarantine pest. Manipulating plant-associated microbes to control nematodes is an important strategy for sustainable pest management. However, it has proven difficult to find pine-associated bacteria that possess both nematocidal activity and the ability to colonize pine tissues.

Methods

The stress experiments with turpentine and pine tissue extract were carried out to screen for the desired target strain that could adapt to the internal environment of pine trees. This strain was used to construct an engineered nematocidal strain. Additionally, a fluorescent strain was constructed to determine its dispersal ability in Pinus massoniana seedlings through plate separation, PCR detection, and fluorescence microscopy observations. The engineered nematocidal strain was tested in the greenhouse experiment to assess its ability to effectively protect P. massoniana seedlings from nematode infection.

Results

This study isolated a Bacillus toyonensis strain Bxy19 from the healthy pine stem, which showed exceptional tolerance in stress experiments. An engineered nematocidal strain Bxy19P3C6 was constructed, which expressed the Cry6Aa crystal protein and exhibited nematocidal activity. The fluorescent strain Bxy19GFP was also constructed and used to test its dispersal ability. It was observed to enter the needles of the seedlings through the stomata and colonize the vascular bundle after being sprayed on the seedlings. The strain was observed to colonize and spread in the tracheid after being injected into the stems. The strain could colonize the seedlings and persist for at least 50 days. Furthermore, the greenhouse experiments indicated that both spraying and injecting the engineered strain Bxy19P3C6 had considerable efficacy against nematode infection.

Discussion

The evidence of the colonization ability and persistence of the strain in pine advances our understanding of the control and prediction of the colonization of exogenously delivered bacteria in pines. This study provides a promising approach for manipulating plant-associated bacteria and using Bt protein to control nematodes.

Transcriptomic response of Pinus massoniana to infection stress from the pine wood nematode Bursaphelenchus xylophilus

The pinewood nematode (PWN) Bursaphelenchus xylophilus is a forestry quarantine pest and causes an extremely dangerous forest disease that is spreading worldwide. Due to the complex pathogenic factors of pine wood nematode disease, the pathogenesis is still unknown. B. xylophilus ultimately invades a host and causes death. However, little is known about the defence-regulating process of host pine after infection by B. xylophilus at the molecular level. Therefore, we wanted to understand how Pinus massoniana regulates its response to invasion by B. xylophilus. P. massoniana were artificially inoculated with B. xylophilus solution, while those without B. xylophilus solution were used as controls. P. massoniana inoculated with B. xylophilus solution for 0 h, 6 h, 24 h, and 120 h was subjected to high-throughput sequencing to obtain transcriptome data. At various time points (0 h, 6 h, 24 h, 120 h), gene transcription was measured in P. massoniana inoculated with PWN. At different time points, P. massoniana gene transcription differed significantly, with a response to early invasion by PWN. According to Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, P. massoniana response to PWN invasion involves a wide range of genes, including plant hormone signal transformation, flavonoid biosynthesis, amino sugar and nucleoside sugar metabolism, and MAPK signalling pathways. Among them, inoculation for 120 hours had the greatest impact on differential genes. Subsequently, weighted gene coexpression network analysis (WGCNA) was used to analyse transcriptional regulation of P. massoniana after PWN infection. The results showed that the core gene module of P. massoniana responding to PWN was "MEmagenta", enriched in oxidative phosphorylation, amino sugar and nucleotide sugar metabolism, and the MAPK signalling pathway. MYB family transcription factors with the highest number of changes between infected and healthy pine trees accounted for 20.4% of the total differentially expressed transcription factors. To conclude, this study contributes to our understanding of the molecular mechanism of initial PWN infection of P. massoniana. Moreover, it provides some important background information on PWN pathogenic mechanisms.

Effects of Colletotrichum gloeosporioides and Poplar Secondary Metabolites on the Composition of Poplar Phyllosphere Microbial Communities

Poplar anthracnose caused by Colletotrichum gloeosporioides is a common disease affecting poplars globally that causes the destruction and alteration of poplar phyllosphere microbial communities; however, few studies have investigated these communities. Therefore, in this study, three species of poplar with different resistances were investigated to explore the effects of Colletotrichum gloeosporioides and poplar secondary metabolites on the composition of poplar phyllosphere microbial communities. Evaluation of the phyllosphere microbial communities before and after inoculation of the poplars with C. gloeosporioides revealed that both bacterial and fungal OTUs decreased after inoculation. Among bacteria, the most abundant genera were Bacillus, Plesiomonas, Pseudomonas, Rhizobium, Cetobacterium, Streptococcus, Massilia, and Shigella for all poplar species. Among fungi, the most abundant genera before inoculation were Cladosporium, Aspergillus, Fusarium, Mortierella, and Colletotrichum, while Colletotrichum was the main genus after inoculation. The inoculation of pathogens may regulate the phyllosphere microorganisms by affecting the secondary metabolites of plants. We investigated metabolite contents in the phyllosphere before and after the inoculation of the three poplar species, as well as the effects of flavonoids, organic acids, coumarins, and indoles on poplar phyllosphere microbial communities. We speculated that coumarin had the greatest recruitment effect on phyllosphere microorganisms, followed by organic acids through regression analysis. Overall, our results provide a foundation for subsequent screening of antagonistic bacteria and fungi against poplar anthracnose and investigations of the mechanism by which poplar phyllosphere microorganisms are recruited. IMPORTANCE Our findings revealed that the inoculation of Colletotrichum gloeosporioides has a greater effect on the fungal community than the bacterial community. In addition, coumarins, organic acids, and flavonoids may have recruitment effects on phyllosphere microorganisms, while indoles may have inhibitory effects on these organisms. These findings may provide the theoretical basis for the prevention and control of poplar anthracnose.

Transcriptome analysis reveals the high ribosomal inhibitory action of 1,4-naphthoquinone on Meloidogyne luci infective second-stage juveniles

The root-knot nematode (RKN) Meloidogyne luci presents a threat to the production of several important crops. This nematode species was added to the European Plant Protection Organization Alert list in 2017. The scarce availability of efficient nematicides to control RKN and the phasing out of nematicides from the market have intensified the search for alternatives, such as phytochemicals with bionematicidal properties. The nematicidal activity of 1,4-naphthoquinone (1,4-NTQ) against M. luci has been demonstrated; however, knowledge of the potential mode(s) of action of this compound is still scarce. In this study, the transcriptome profile of M. luci second-stage juveniles (J2), the infective stage, in response to 1,4-NTQ exposure was determined by RNA-seq to identify genes and pathways that might be involved in 1,4-NTQ's mode(s) of action. Control treatments, consisting of nematodes exposed to Tween^® 80 (1,4-NTQ solvent) and to water, were included in the analysis. A large set of differentially expressed genes (DEGs) was found among the three tested conditions, and a high number of downregulated genes were found between 1,4-NTQ treatment and water control, reflecting the inhibitory effect of this compound on M. luci, with a great impact on processes related to translation (ribosome pathway). Several other nematode gene networks and metabolic pathways affected by 1,4-NTQ were also identified, clarifying the possible mode of action of this promising bionematicide.

Importance of ABC Transporters in the Survival of Parasitic Nematodes and the Prospect for the Development of Novel Control Strategies

ABC transporters, a family of ATP-dependent transmembrane proteins, are responsible for the active transport of a wide range of molecules across cell membranes, including drugs, toxins, and nutrients. Nematodes possess a great diversity of ABC transporters; however, only P-glycoproteins have been well-characterized compared to other classes. The ABC transport proteins have been implicated in developing resistance to various classes of anthelmintic drugs in parasitic nematodes; their role in plant and human parasitic nematodes still needs further investigation. Therefore, ABC transport proteins offer a potential opportunity to develop nematode control strategies. Multidrug resistance inhibitors are becoming more attractive for controlling nematodes due to their potential to increase drug efficacy in two ways: (i) by limiting drug efflux from nematodes, thereby increasing the amount of drug that reaches its target site, and (ii) by reducing drug excretion by host animals, thereby enhancing drug bioavailability. This article reviews the role of ABC transporters in the survival of parasitic nematodes, including the genes involved, their regulation and physiological roles, as well as recent developments in their characterization. It also discusses the association of ABC transporters with anthelmintic resistance and the possibility of targeting them with next-generation inhibitors or nutraceuticals (e.g., polyphenols) to control parasitic infections.

Effect of Associated Bacteria GD1 on the Low-Temperature Adaptability of Bursaphelenchus xylophilus Based on RNA-Seq and RNAi

To explore the effect of associated bacteria on the low-temperature adaptability of pinewood nematodes (PWNs), transcriptome sequencing (RNA-seq) of PWN AH23 treated with the associated bacterial strain Bacillus cereus GD1 was carried out with reference to the whole PWN genome. Bioinformatic software was utilized to analyze the differentially expressed genes (DEGs). This study was based on the analysis of DEGs to verify the function of daf-11 by RNAi. The results showed that there were 439 DEGs between AH23 treated with GD1 and those treated with ddH₂O at 10 °C. There were 207 pathways annotated in the KEGG database and 48 terms annotated in the GO database. It was found that after RNAi of daf-11, the survival rate of PWNs decreased significantly at 10 °C, and fecundity decreased significantly at 15 °C. It can be concluded that the associated bacteria GD1 can enhance the expression of genes related to PWN low-temperature adaptation and improve their adaptability to low temperatures.

Comprehensive Analysis of Copy Number Variations on Glycoside Hydrolase 45 Genes among Different Bursaphelenchus xylophilus Strains

Bursaphelenchus xylophilus is considered the most dangerous quarantine pest in China. It causes enormous economic and ecological losses in many countries from Asia and Europe. The glycoside hydrolase 45 gene family has been demonstrated in early studies to contribute to the cell wall degradation ability of B. xylophilus during its infection. However, the copy number variation (CNV) of the GH45 gene and its association with B. xylophilus pathogenicity were not fully elucidated. In this study, we found that the GH45 gene with two copies is the most predominant type among 259 B. xylophilus strains collected from China and Japan. Additionally, 18 strains are identified as GH45 genes with a single copy, and only two strains are verified to have three copies. Subsequent expression analysis and inoculation test suggest that the copy numbers of the GH45 gene are correlated with gene expression as well as the B. xylophilus pathogenicity. B. xylophilus strains with more copies of the GH45 gene usually exhibit more abundant expression and cause more severe wilt symptoms on pine trees. The aforementioned results indicated the potential regulatory effects of CNV in B. xylophilus and provided novel information to better understand the molecular pathogenesis of this devastating pest.

Effects on community composition and function Pinus massoniana infected by Bursaphelenchus xylophilus

Pine wilt disease (PWD) is a worldwide forest disease caused by pine wood nematode (PWN). In this article, we investigated the composition, organization, correlation, and function of the endophytic microbial community in Pinus massoniana field with and without PWN. Samples were taken from branches, upper, middle, and lower trunks, as well as soil, from both healthy and infected trees. The results showed that the fungal diversity of healthy pines is around 1.1 times that of infected pines, while the bacterial diversity is about 0.75 times that of infected pines at the OTUs level. An increase of the abundance of pathogenic fungus such as Saitozyma, Graphilbum, Diplodia, Candida, Pseudoxanthomonas, Dyella and Pantoea was witnessed in infected pines according to the result of LEfSe. Furthermore, Ophiostoma and saprophytic fungus such as Entomocorticium, ganoderma, tomentella, entomocorticium were exclusively prominent in infected pines, which were substantially and highly connected with other species (p < 0.05), indicating the trees' vulnerability and making the wood blue. In healthy pines, the top three functional guilds are parasites, plant pathogens, and saprotrophs. Parasites (36.52%) are primarily found in the branches, plant pathogens (29.12%) are primarily found in the lower trunk, and saprotrophs (67.88%) are primarily found in the upper trunk of disease trees. Pines' immunity is being eroded due to an increase in the quantity and types of diseases. PICRUSt2 research revealed that NADH or NADPH, as well as carbon-nitrogen bonds, were more abundant in healthy pines, but acid anhydrides and transferring phosphorus-containing groups were more abundant in infected pines. The shift in resin secretion lowers the tree's potential and encourages pine wilt and mortality. In total, PWN may have disrupted the microbiological ecology and worked with the community to hasten the demise of pines.

Selection and Validation of Reliable Reference Genes for qRT-PCR Normalization of Bursaphelenchus xylophilus from Different Temperature Conditions and Developmental Stages

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a powerful technique for studying gene expression. The key to quantitative accuracy depends on the stability of the reference genes used for data normalization under different experimental conditions. Pine wood nematode (PWN; Bursaphelenchus xylophilus) is the causal agent of the devastating pine wilt disease (PWD). Extensive and prompt research is needed to understand the molecular mechanism of PWD, but identification of the reference PWN genes for standardized qRT-PCR has not been reported yet. We have analyzed eight candidate reference genes of PWN across different temperature conditions and developmental stages. Delta Ct method, GeNorm, NormFinder, BestKeeper, and RefFinder algorithms were used to evaluate the stability of expression of these genes. Finally, we use arginine kinase gene (AK) in different temperatures and heat shock protein 90 (HSP90) in different developmental stages to confirm the stability of these genes. UBCE and EF1γ were most stable across different temperature treatments, whereas EF1γ and Actin were most stable across different developmental stages. In general, these results indicate that EF1γ is the most stable gene for qRT-PCR under different conditions. The systematic analysis of qRT-PCR reference gene selection will be helpful for future functional analysis and exploration of B. xylophilus genetic resources.

Enhanced accumulation of pinosylvin stilbenes and related gene expression in Pinus strobus after infection of pine wood nematode

Pine wood nematodes (PWNs; Bursaphelenchus xylophilus) infect pine trees and cause serious pine wilt disease. Eastern white pine (Pinus strobus) has resistance to PWN. However, the detailed defense mechanisms of P. strobus against PWN are not well known. When P. strobus plants were infected with PWNs, the accumulation of stilbenoids, dihydropinosylvin monomethyl ether (DPME) and pinosylvin monomethyl ether (PME) was increased remarkably. Both DPME and PME had high nematicidal activity. The nematicidal activity of the two compounds was resulted in a developmental stage-dependent manner. Pinosylvin monomethyl ether was more toxic to adult PWNs than juveniles, whereas DPME was found more toxic to juvenile PWNs than the adults. The genes involved in PME and DPME biosynthesis such as phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), pinosylvin synthase (STS) and pinosylvin O-methyltransferase (PMT) were isolated using de novo sequencing of the transcriptome in P. strobus. In addition, transcription factors (TFs; bHLH, MYB and WRKY) related to stilbene biosynthesis were isolated. qPCR analyses of the selected genes (PAL, 4CL, STS and PMT) including TFs (bHLH, MYB and WRKY) revealed that the expression level of the selected genes highly enhanced after PWN infection. Our results suggest that pinosylvin-type stilbenoid biosynthesis is highly responsive to PWN infection and plays an important role in PWN resistance of P. strobus trees.

Copy Number Variations of Glycoside Hydrolase 45 Genes in Bursaphelenchus xylophilus and Their Impact on the Pathogenesis of Pine Wilt Disease

The pine wood nematode Bursaphelenchus xylophilus parasitizes millions of pine trees worldwide each year, causing severe wilt and the death of host trees. Glycoside hydrolase 45 genes of B. xylophilus are reported to have been acquired by horizontal gene transfer from fungi and are responsible for cell wall degradation during nematode infection. Previous studies ignored the possibility of copy number variations of such genes. In this study, we determined that two of the glycoside hydrolase 45 genes evolved to maintain multiple copies with distinct expression levels, enabling the nematode to infect a variety of pine hosts. Additionally, tandem repeat variations within coding regions were also detected between different copies of glycoside hydrolase 45 genes that could result in changes in protein sequences and serve as an effective biological marker to detect copy number variations among different B. xylophilus populations. Consequently, we were able to further identify the copy number variations of glycoside hydrolase 45 genes among B. xylophilus strains with different virulence. Our results provide new insights into the pathogenicity of B. xylophilus, provide a practical marker to genotype copy number variations and may aid in population classification.

Comparative Analysis of Bursaphelenchus xylophilus Secretome Under Pinus pinaster and P. pinea Stimuli

The pinewood nematode (PWN), Bursaphelenchus xylophilus, the pine wilt disease's (PWD) causal agent, is a migratory endoparasitic nematode skilled to feed on pine tissues and on fungi that colonize the trees. In order to study B. xylophilus secretomes under the stimulus of pine species with different susceptibilities to disease, nematodes were exposed to aqueous pine extracts from Pinus pinaster (high-susceptible host) and P. pinea (low-susceptible host). Sequential windowed acquisition of all theoretical mass spectra (SWATH-MS) was used to determine relative changes in protein amounts between B. xylophilus secretions, and a total of 776 secreted proteins were quantified in both secretomes. From these, 22 proteins were found increased in the B. xylophilus secretome under the P. pinaster stimulus and 501 proteins increased under the P. pinea stimulus. Functional analyses of the 22 proteins found increased in the P. pinaster stimulus showed that proteins with peptidase, hydrolase, and antioxidant activities were the most represented. On the other hand, gene ontology (GO) enrichment analysis of the 501 proteins increased under the P. pinea stimulus revealed an enrichment of proteins with binding activity. The differences detected in the secretomes highlighted the diverse responses from the nematode to overcome host defenses with different susceptibilities and provide new clues on the mechanism behind the pathogenicity of this plant-parasitic nematode. Proteomic data are available via ProteomeXchange with identifier PXD024011.

Physiological effect of colloidal carbon quantum dots on Bursaphelenchus xylophilus

Bursaphelenchus xylophilus (B. xylophilus) is a dangerous plant pest which could result in Pine Wild Disease (PWD).

Molecular variation among virulent and avirulent strains of the quarantine nematode Bursaphelenchus xylophilus

Bursaphelenchus xylophilus is an emerging pathogenic nematode that is responsible for a devastating epidemic of pine wilt disease worldwide, causing severe ecological damage and economic losses to forestry. Two forms of this nematode have been reported, i.e., with strong and weak virulence, commonly referred as virulent and avirulent strains. However, the pathogenicity-related genes of B. xylophilus are not sufficiently characterized. In this study, to find pathogenesis related genes we re-sequenced and compared genomes of two virulent and two avirulent populations. We identified genes affected by genomic variation, and functional annotation of those genes indicated that some of them might play potential roles in pathogenesis. The performed analysis showed that both avirulent populations differed from the virulent ones by 1576 genes with high impact variants. Demonstration of genetic differences between virulent and avirulent strains will provide effective methods to distinguish these two nematode virulence forms at the molecular level. The reported results provide basic information that can facilitate development of a better diagnosis for B. xylophilus isolates/strains which present different levels of virulence and better understanding of the molecular mechanism involved in the development of the PWD.

Adaptation of pine wood nematode Bursaphelenchus xylophilus to β-pinene stress

BACKGROUND

The pine wood nematode (PWN; Bursaphelenchus xylophilus) is the most damaging biological pest in pine forest ecosystems in China. However, the pathogenic mechanism remains unclear. Tracheid cavitation induced by excess metabolism of volatile terpenes is a typical characteristic of pine trees infected by B. xylophilus. β-pinene, one of the main volatile terpenes, influences PWN colonization and reproduction, stimulating pathogenicity during the early stages of infection. To elucidate the response mechanism of PWN to β-pinene, pathogenesis, mortality, and reproduction rate were investigated under different concentrations of β-pinene using a transcriptomics approach.

RESULTS

A low concentration of β-pinene (BL, C < 25.74 mg/ml) inhibited PWN reproduction, whereas a high concentration (BH, C > 128.7 mg/ml) promoted reproduction. Comparison of PWN expression profiles under low (BL, 21.66 mg/ml) and high (BH, 214.5 mg/ml) β-pinene concentrations at 48 h identified 659 and 418 differentially expressed genes (DEGs), respectively, compared with controls. Some key DEGs are potential regulators of β-pinene via detoxification metabolism (cytochrome P450, UDP-glucuronosyltransferases and short-chain dehydrogenases), ion channel/transporter activity (unc and ATP-binding cassette families), and nuclear receptor -related genes. Gene Ontology enrichment analysis of DEGs revealed metabolic processes as the most significant biological processes, and catalytic activity as the most significant molecular function for both BL and BH samples. Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO) analysis showed that xenobiotics biodegradation and metabolism, carbohydrate metabolism, lipid metabolism, amino acid metabolism, metabolism of cofactors and vitamins, and transport and catabolism were the dominant terms in metabolism categories.

CONCLUSION

In addition to detoxification via reduction/oxidation (redox) activity, PWN responds to β-pinene through amino acid metabolism, carbohydrate metabolism, and other pathways including growth regulation and epidermal protein changes to overcome β-pinene stress. This study lays a foundation for further exploring the pathogenic mechanism of PWN.

Potential Molecular Mimicry Proteins Responsive to α-pinene in Bursaphelenchus xylophilus

: Bursaphelenchus xylophilus is a nematode species that has damaged pine trees worldwide, but its pathogenesis has not been fully characterized. α-pinene helps protect host species during the early B. xylophilus infection and colonization stages. In this study, we identified potential molecular mimicry proteins based on a comparative transcriptomic analysis of B. xylophilus. The expression levels of three genes encoding secreted B. xylophilus proteins were influenced by α-pinene. We cloned one gene encoding a thaumatin-like protein, Bx-tlp-2 (accession number MK000287), and another gene encoding a cysteine proteinase inhibitor, Bx-cpi (accession number MK000288). Additionally, α-pinene appeared to induce Bx-tlp-1 expression, but had the opposite effect on Bx-cpi expression. An analysis of the expression of the potential molecular mimicry proteins in B. xylophilus infecting pine trees revealed that the α-pinene content was consistent with the expression levels of Bx-tlp-1 (Bx-cpi) and Pm-tlp (Pm-cpi) over time. Thus, these genes likely have important roles contributing to the infection of pine species by B. xylophilus. The results of this study may be relevant for future investigations of the functions of Bx-tlp-1, Bx-tlp-2 and Bx-cpi, which may provide a point to explore the relationship between B. xylophilus and host pines.

The anti-phytoalexin gene Bx-cathepsin W supports the survival of Bursaphelenchus xylophilus under Pinus massoniana phytoalexin stress

Background

Pine trees challenged by Bursaphelenchus xylophilus invasion produce phytoalexins to combat this nematode. Nevertheless, the phytoalexins of Asian pine trees are ineffective against B. xylophilus. The anti-phytoalexin genes of B. xylophilus disable almost all Asian pine phytoalexins, which has allowed B. xylophilus to devastate pine forests in eastern Asia over the last four decades. However, to date, the factors that stimulate anti-phytoalexin gene expression and the mechanisms by which these genes act are not well understood.

Results

Here, we described anti-phytoalexin genes in B. xylophilus using transcriptomic and bioinformatics analyses. The genes that were induced by both Pinus massoniana and carvone and had similarly elevated expression trends were considered anti-phytoalexin genes. Altogether, 187 anti-phytoalexin genes were identified, including 4 cathepsin genes. KEGG pathway enrichment indicated that those cathepsins were related to the Lysosome pathway. Since cathepsins help to maintain metabolic homeostasis by participating in the degradation of heterophagic and autophagic material, the lysosomal cathepsin gene Bx-cathepsin W was cloned and characterized. The results of the RNAi assessment indicated that the knockdown of Bx-cathepsin W reduced the survival rates of B. xylophilus under carvone or P. massoniana stress. The correlation between Bx-cathepsin W and the susceptibility of pines showed that Bx-cathepsin W might help improve the anti-phytotoxin ability of B. xylophilus.

Conclusions

The results indicated that the anti-phytoalexin gene Bx-cathepsin W supported the survival of B. xylophilus under P. massoniana phytoalexin stress. The cDNA library sequencing, differentially expressed gene identification, and WGCNA algorithm analysis provided insight at a systemic level into the gene regulation of B. xylophilus in response to the immune reaction of P. massoniana. These results will lead to a better understanding of the function of nematode defenses in host innate immunity.

Molecular Characterization and Functional Analysis of Three Autophagy Genes, BxATG5, BxATG9, and BxATG16, in Bursaphelenchus xylophilus

The pine wood nematode (PWN), Bursaphelenchus xylophilus, is the pathogen responsible for pine wilt disease (PWD), a devastating forest disease with a pathogenic mechanism that remains unclear. Autophagy plays a crucial role in physiological and pathological processes in eukaryotes, but its regulatory mechanism and significance in PWN are unknown. Therefore, we cloned and characterized three autophagy genes, BxATG5, BxATG9, and BxATG16, in PWN. BxATG9 and BxATG16 were efficiently silenced through RNA interference, and we found that BxATG16 positively regulated the expression of BxATG5. Silencing BxATG9 and BxATG16 severely inhibited feeding and reproduction in PWN, indicating that autophagy is essential for these processes. We then examined the expression patterns of these three autophagy genes in PWN under the stresses of α-pinene and H₂O₂, the main defense substances of pine trees, and during the development of PWD using quantitative reverse transcription polymerase chain reaction. The expression levels of BxATG5, BxATG9, and BxATG16 all significantly increased after nematodes were stressed with α-pinene and H₂O₂ and inoculated into pine trees, suggesting that autophagy plays an important role in the defense and pathogenesis of PWN. In this study, the molecular characteristics and functions of the autophagy genes BxATG5, BxATG9, and BxATG16 in PWN were elucidated.