Identifying Virulence-Associated Genes Using Transcriptomic and Proteomic Association Analyses of the Plant Parasitic Nematode Bursaphelenchus mucronatus

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.

madd-4 plays a critical role in light against Bursaphelenchus xylophilus

Bursaphelenchus xylophilus is a notorious invasive species, causing extensive losses to pine ecosystems globally. Previous studies had shown that the development of B. xylophilus was seriously suppressed by light. However, the mechanism involved in the inhibition is unknown. Here, it is the first report that Bxy-madd-4 is a light-regulated gene, plays a potential role in B. xylophilus in responding to the blue light. Transcriptome sequencing revealed that the expression level of Bxy-madd-4 declined by 86.39% under blue light. The reverse transcription quantitative real-time PCR results were in accord with the transcriptome sequencing, confirming the expression level of Bxy-madd-4 was suppressed by blue light. Bxy-madd-4 promoter::mCherry reporter constructed in Caenorhabditis elegans were utilized to mimic the spatiotemporal expression patterns of Bxy-madd-4. Bxy-madd-4A promoter activity had a strong continuity throughout all development stages in C. elegans. Further RNA interference indicated that only 36.8% of the Bxy-madd-4 dsRNA treated embryos were hatched. Moreover, 71.6% of the hatched nematodes were abnormal, such as particles on the body surface and concave tissues. Our findings contribute towards a better understanding of the mechanism of light against the destructive invasive nematode, providing a promising hint for control of the destructive invasive nematode.

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.

Function of lipid binding proteins of parasitic helminths: still a long road

Infections with parasitic helminths cause severe debilitating and sometimes lethal diseases in humans and domestic animals on a global scale. Unable to synthesize de novo their own fatty acids and sterols, helminth parasites (nematodes, trematodes, cestodes) rely on their hosts for their supply. These organisms produce and secrete a wide range of lipid binding proteins that are, in most cases, structurally different from the ones found in their hosts, placing them as possible novel therapeutic targets. In this sense, a lot of effort has been made towards the structure determination of these proteins, but their precise function is still unknown. In this review, we aim to present the current knowledge on the functions of LBPs present in parasitic helminths as well as novel members of this highly heterogeneous group of proteins.

The role of mab-3 in spermatogenesis and ontogenesis of pinewood nematode, Bursaphelenchus xylophilus

BACKGROUND

Bursaphelenchus xylophilus is one of the most destructive invasive species, causing extensive economic losses worldwide. The sex ratio of female to male of B. xylophilus plays an important role in the nematode infestation. However, little is known about the processes of its sex determination. The double sex/mab-3-related family of transcription factors are highly conserved in animals, playing crucial roles in sex determination, spermatogenesis and ontogenesis. We therefore investigated its orthologue, Bxy-mab-3, in B. xylophilus.

RESULTS

Bxy-mab-3 has two typical conserved DNA-binding domains. It was observed in J2 (the second-stage of juveniles), J3, J4 and male adults (specifically on the spicules), but not in eggs or female adults via mRNA in situ hybridization. RNA-Seq indicated significantly higher expression in males. RNAi showed that the body size and sperm size of male adults were markedly smaller than those of the controls. Meanwhile, almost all the RNAi-treated males failed to mate with the normal females, even 26.34% of interfered males did not produce sperm. However, RNAi of Bxy-mab-3 had no effect on the sex ratio of B. xylophilus.

CONCLUSION

Bxy-mab-3 is indispensable for spermatogenesis, ontogenesis and mating behavior. It is a typical sex-determination gene with differential expression in males and females. However, knocking down Bxy-mab-3 expression could not alter the sex ratio as seen in other species. Our findings contribute towards a better understanding of the molecular events of Bxy-mab-3 in B. xylophilus, which provides promising hints for control of pine wilt disease by blocking ontogenesis and decreasing nematode fecundity.

Comparative Transcriptome Analysis of Pinus densiflora Following Inoculation with Pathogenic (Bursaphelenchus xylophilus) or Non-pathogenic Nematodes (B. thailandae)

Pinus densiflora (Korean red pine) is a species of evergreen conifer that is distributed in Korea, Japan, and China, and of economic, scientific, and ecological importance. Korean red pines suffer from pine wilt disease (PWD) caused by Bursaphelenchus xylophilus, the pinewood nematode (PWN). To facilitate diagnosis and prevention of PWD, studies have been conducted on the PWN and its beetle vectors. However, transcriptional responses of P. densiflora to PWN have received less attention. Here, we inoculated Korean red pines with pathogenic B. xylophilus, or non-pathogenic B. thailandae, and collected cambium layers 4 weeks after inoculation for RNA sequencing analysis. We obtained 72,864 unigenes with an average length of 869 bp (N50 = 1,403) from a Trinity assembly, and identified 991 differentially expressed genes (DEGs). Biological processes related to phenylpropanoid biosynthesis, flavonoid biosynthesis, oxidation–reduction, and plant-type hypersensitive response were significantly enriched in DEGs found in trees inoculated with B. xylophilus. Several transcription factor families were found to be involved in the response to B. xylophilus inoculation. Our study provides the first evidence of transcriptomic differences in Korean red pines inoculated with B. xylophilus and B. thailandae, and might facilitate early diagnosis of PWD and selection of PWD-tolerant Korean red pines.

Proteomics (SWATH-MS) informed by transcriptomics approach of tropical herb Persicaria minor leaves upon methyl jasmonate elicitation

BACKGROUND

Jasmonic acid (JA) and its derivative, methyl JA (MeJA) are hormonal cues released by plants that signal defense response to curb damages from biotic and abiotic stresses. To study such response, a tropical herbal plant, Persicaria minor, which possesses pungent smell and various bioactivities including antimicrobial and anticancer, was treated with MeJA. Such elicitation has been performed in hairy root cultures and plants such as Arabidopsis and rice, yet how MeJA influenced the proteome of an herbal species like P. minor is unknown.

METHOD

In this study, P. minor plants were exogenously elicited with MeJA and leaf samples were subjected to SWATH-MS proteomics analysis. A previously published translated transcriptome database was used as a reference proteome database for a comprehensive protein sequence catalogue and to compare their differential expression.

RESULTS

From this proteomics informed by transcriptomics approach, we have successfully profiled 751 proteins of which 40 proteins were significantly different between control and MeJA-treated samples. Furthermore, a correlation analysis between both proteome and the transcriptome data sets suggests that significantly upregulated proteins were positively correlated with their cognate transcripts (Pearson's r = 0.677) while a weak correlation was observed for downregulated proteins (r = 0.147).

DISCUSSION

MeJA treatment induced the upregulation of proteins involved in various biochemical pathways including stress response mechanism, lipid metabolism, secondary metabolite production, DNA degradation and cell wall degradation. Conversely, proteins involved in energy expensive reactions such as photosynthesis, protein synthesis and structure were significantly downregulated upon MeJA elicitation. Overall protein-transcript correlation was also weak (r = 0.341) suggesting the existence of post-transcriptional regulation during such stress. In conclusion, proteomics analysis using SWATH-MS analysis supplemented by the transcriptome database allows comprehensive protein profiling of this non-model herbal species upon MeJA treatment.

Selection of Reliable Reference Genes for RT-qPCR Analysis of Bursaphelenchus mucronatus Gene Expression From Different Habitats and Developmental Stages

Quantitative reverse transcription polymerase chain reaction (RT-qPCR), a sensitive technique for gene expression analysis, depends on the stability of the reference genes used for data normalization under different experimental conditions. Bursaphelenchus mucronatus, a pine-parasitic nematode varying in virulence, is widely distributed in natural pine forests throughout the northern hemisphere, but has not been investigated with respect to the identification of reference genes suitable for the normalization of RT-qPCR data. In the present study, eight candidate reference genes were analyzed in B. mucronatus under different habitat conditions and at different developmental stages. The expression stability of these genes was assessed by geNorm, NormFinder, BestKeeper, delta Cq, and RefFinder algorithms. In general, our results identified encoding beta-tubulin as the most stable gene. Moreover, pairwise analysis showed that three reference genes were sufficient to normalize the gene expression data under each set of conditions, with genes encoding beta-tubulin, 18S ribosomal RNA and ubiquitin-conjugating enzyme being the most suitable reference genes for different habitat conditions, whereas genes encoding beta-tubulin, histone, and 18S ribosomal RNA exhibited the most stable expression at different developmental stages. Validation of the selected reference genes was performed by profiling the expression of the fatty acid- and retinol-binding protein gene in different habitats, and by profiling the expression of the arginine kinase gene at different developmental stages. This first systematic analysis for the selection of suitable reference genes for RT-qPCR in B. mucronatus will facilitate future functional analyses and deep mining of genetic resources in this nematode.