Next Article in Journal
An Atypical Etiology of Fungal Keratitis Caused by Roussoella neopustulans
Next Article in Special Issue
Regulatory Roles of Histone Modifications in Filamentous Fungal Pathogens
Previous Article in Journal
Diagnosis from Tissue: Histology and Identification
Previous Article in Special Issue
The Arf-GAP Proteins AoGcs1 and AoGts1 Regulate Mycelial Development, Endocytosis, and Pathogenicity in Arthrobotrys oligospora
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 8
Issue 5
10.3390/jof8050506
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Entomopathogenic Fungus Metarhizium anisopliae Affects Feeding Preference of Sogatella furcifera and Its Potential Targets’ Identification

by
Yirong Wang
1,2,
Lijuan Han
1,2,
Yuxian Xia
1,2,* and
Jiaqin Xie
1,2,*
1
School of Life Sciences, Genetic Engineering Research Center, Chongqing University, Chongqing 401331, China
2
Chongqing Engineering Research Center for Fungal Insecticides, and Key Laboratory of Gene Function and Regulation Technology under Chongqing Municipal Education Commission, Chongqing 401331, China
*
Authors to whom correspondence should be addressed.
J. Fungi 2022, 8(5), 506; https://doi.org/10.3390/jof8050506
Submission received: 25 April 2022 / Revised: 11 May 2022 / Accepted: 12 May 2022 / Published: 15 May 2022
(This article belongs to the Special Issue Interactions between Filamentous Fungal Pathogens and Hosts)
Browse Figures
Versions Notes

Abstract

:
The rice planthopper Sogatella furcifera is a unique vector of the southern rice black-streaked dwarf virus (SRBSDV). The feeding behavior of S. furcifera should directly affect the diffusion of this virus. In this study, we noted that the infection of Metarhizium anisopliae CQMa421 on S. furcifera disturbed the feeding behavior of this pest to SRBSDV-infected rice, from preference to non-preference. Then, we further investigated the potential targets of M. anisopliae CQMa421 on the feeding behavior of S. furcifera after 0 h, 24 h and 48 h of infection by transcriptomic analysis via Illumina deep sequencing. A total of 93.27 GB of data was collected after sequencing, from which 91,125 unigenes were annotated, including 75 newly annotated genes. There were 1380 vs. 2187 and 137 vs. 106 upregulated and downregulated differentially expressed genes (DEGs) detected at 24 h and 48 h, respectively. The biological functions and associated metabolic processes of these genes were determined with the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The results suggested that major of DEGs are involved in energy metabolism, biosynthesis, immune response, the FoxO signaling pathway, the MAPK signaling pathway and apoptosis in response to the fungal infection. Noteworthily, several olfactory-related genes, including odorant receptors and odorant binding proteins, were screened from these differentially expressed genes, which played critical roles in regulating the olfactory behavior of insects. Taken together, these results provide new insights for understanding the molecular mechanisms underlying fungus and host insect interaction, especially for olfactory behavior regulated by fungus.

1. Introduction

Rice is a staple crop, feeding more than half of the global population worldwide [1]. During the rice-growing period, it should face a large number of pest and disease threats, which seriously damage the quality and quantity of rice food [2,3]. The rice planthoppers, including Sogatella furcifera, Nilaparvata lugens and Laodelphax striatellus, are one of the most destructive pests in rice plants [2,4,5]. Those pests not only suck the sap out of rice straw to result in the death of rice plants, but also transmit several plant viruses, such as the rice ragged stunt virus, the rice grassy stunt virus and the southern rice black-streaked dwarf virus (SRBSDV) [3,6]. Moreover, these pests have strong reproductive and dispersal abilities, allowing them to damage the rice plants for the duration of the growing period. A previous study has reported that the occurrence of the rice virus SRBSDV keeps in line with the dispersal and outbreak of the S. furcifera population [7]. The damages caused by the plant virus is more severe than that caused by pests in some cases [8]. Thus, an effective method to control the proliferation of rice viruses is to block out the transmission of vectors (e.g., the rice planthoppers).
Entomopathogens commonly exist under natural conditions and play key roles in regulating the population of insects [9,10]. Compared to chemical insecticides, the fungal pesticides are environmentally friendly biocontrol agents with less risks for nontarget organisms [10]. Some entomopathogenic fungi, such as Metarhizium anisopliae and Beauveria bassiana, have been used for the control of agricultural and forestry pests, including Chironomus riparius [11], Locusta migratoria [12], Helicoverpa armigera [13], Alphitobius diaperinus [14] and Nilaparvata lugens [15,16]. Combined use of M. anisopliae and dsRNA or chemicals has shown good potential for the control of the rice pest N. lugens [16,17]. A specific M. anisopliae strain has also been applied for the control of resistant insect pests [18]. Moreover, industrial products of M. anisopliae have been reported that could suppress the population of rice planthoppers to low levels in large-scale applications [15].
Entomopathogenic fungi could directly penetrate the cuticles of insect pests by secreting a few proteases and chitinases under the function of turgor pressure [19]. After penetrating into the host body, rapid proliferation of fungal hyphal bodies in haemocoel would deprive the host of nutrients and result in its death [20]. Due to the immune response of the insect host, fungi would secrete several toxins to interfere with the host nervous system and suppress such immune responses, so as to benefit the reproduction of the hyphal body of the fungus [21,22]. Correspondingly, the number and type of host insect hemocytes would change to better interact with the invaded pathogens during the processes [23]. Meanwhile, the energy metabolism of the host is intensive, and many antimicrobial peptides are generated to function on the pathogens [24]. Although a few studies have examined the interaction between insect hosts and pathogens under different conditions by transcriptomic profiling [25,26,27], the means by which the rice planthopper S. furcifera responds to fungal infection is less known, especially in terms of their behaviors.
The olfactory behavior of insects is critical for feeding, mating, oviposition and avoiding natural enemies [28,29]. Those insects have evolved sophisticated, sensitive and specific olfactory systems to detect and discriminate amongst an enormous variety of odorants under complex environments [30]. They could identify special chemical odor molecules and locate the odor source by olfactory sensation and exhibit a behavioral response correspondingly [28,30]. For instance, the tobacco hawkmoths Manduca sexta preferentially oviposit on the jimson weed Datura wrightii due to the α-copaene induced via potato beetle Lema daturaphila infection [31]; the 4-vinylanisole is an important pheromone that could cause an aggregation behavior of the migratory locust Locusta migratoria [32]. The rice planthopper S. furcifera also shows a feeding preference for rice plants infected by the SRBSDV virus [33]. Generally, the olfactory behavior of insects is related to the expression level of certain genes, including the odorant binding proteins (OBPs), the odorant receptors (ORs) and the chemosensory proteins (CSPs) [28,34,35]. The knockdown or RNAi of those olfactory-related genes (e.g., OBPs, ORs or CSPs) would directly affect the feeding or oviposition behavior in insects [36].
The white-backed S. furcifera is a major vector of SRBSDV, and outbreaks of this pest may cause the prevalence of this virus [37]. Interestingly, this pest shows a preference to those rice plants that have been infected by SRBSDV, compared to non-SRBSDV-infected plants [33], which would benefit the transmission of this virus. Our previous study has reported that the fungal M. anisopliae CQMa421 could infect adults and nymphs of rice planthoppers [16,38]. However, whether infection of M. anisopliae would affect the behavior of S. furcifera is unknown. Thus, we first examined the feeding behavior of the rice pest S. furcifera after M. anisopliae infection. Transcriptomic profiling is a useful tool to investigate and screen for potential targets under different conditions [39,40]. To understand the potential mechanisms, and to seek the potential olfactory-related targets, we further studied the response of S. furcifera to M. anisopliae infection under different periods by transcriptomic analysis. The potential olfactory-related genes that might affect the feeding behavior of S. furcifera also have been screened from the DEGs. This study provides new insights for understanding the underlying mechanisms of an insect host interacting with a fungal infection, especially for the olfactory behavior regulated by the fungus. These results would improve the potential control strategies for the vector S. furcifera and its transmitted virus by integrated pest management.

2. Methods and Materials

2.1. Insect and Plant Culture

The rice planthopper S. furcifera used in this study was maintained and obtained from the Genetic Engineering Research Center, Chongqing University, Chongqing, China. S. The furcifera individuals were reared on rice seedings in cages (35 cm × 35 cm × 40 cm) at an insectary (T = 27 ± 1 °C; light:dark = 14 h:10 h). The SRBSDV rice plant was originally collected from a rice paddy in Nanxiong, Guangdong, China. To cultivate the SRBSDV rice plant, we first collected S. furcifera nymphs from SRBSDV rice plants to acquire this virus and then transferred them to a new normal (non-SRBSDV) rice plant in a tube for 2 or 3 days. In this way, the new normal rice plant would be infected by SRBSDV via S. furcifera nymphs. Prior to the next experiments, S. furcifera females and males were transferred to a cage (35 cm × 35 cm × 40 cm) containing new rice seedings for 48 h. We then collected the newly hatched nymphs for the feeding choice experiment and transcriptomic analysis.

2.2. Effects of M. anisopliae on S. furcifera Feeding Behavior

We collected S. furcifera adults after 5-day emergence for two feeding preference choices and random and Y-tube olfactory choice tests. In the random choice test, two differently treated rice plants (SRBSDV rice and non-SRBSDV rice plant) were transferred in a catercorner of flowerpots sealed with a net. Then, S. furcifera individuals were put in these flowerpots, and the number of insects on each rice plant was recorded at 3 h, 6 h, 18 h and 24 h. In the Y-tube olfactory choice experiment, we referenced the commonly used methods of a similar, previous study [41]. The number of S. furcifera in each arm of the Y-tube was recorded after staying more than 15 s. We then counted the total number of S. furcifera in each arm within 2 h.
The M. anisopliae CQMa421 strain (China General Microbiological Culture Collection Center, CGMCC No. 460) used in these experiments was obtained from the Genetic Engineering Research Center, Chongqing University, Chongqing, China. The fungal strain M. anisopliae CQMa421 was originally isolated from Cnaphalocrocis medinalis and stored in a −80 °C freezer in our laboratory. To evaluate the effects of M. anisopliae CQMa421 on the feeding behavior of S. furcifera, the S. furcifera adults were initially treated with 1 × 107 conidia/mL suspension of M. anisopliae, which was prepared with the methods mentioned in our previous study [16]. Briefly, the M. anisopliae conidia were harvested after 15 days of growth in ¼ SDAY medium, which included 10 g glucose, 5 g yeast extract, 2.5 g peptone and 18 g agar/liter of sterilized water. Then, 1 × 107 conidia/mL of M. anisopliae conidia was prepared using 0.1% Tween 80 by hemocytometer. A total of 1 mL of such conidial suspension was applied as a spray to the rice seedings by a Potter Precision Spray Tower (Burkard Manufacturing, UK). Sequentially, S. furcifera adults were transferred into the cage maintaining M. anisopliae-treated rice seedings and inoculated on them. Then, we used random choice and Y-tube olfactory choice to examine the feeding choice of S. furcifera for SRBSDV and non-SRBSDV rice plants at 48 h post-inoculation. The consequences of feeding choice in two tests were recorded according to the above methods.

2.3. Total RNA Isolation, Quantification and Sequencing

To further identify the potential olfactory-related targets, the S. furcifera individuals from the M. anisopliae-treated or control groups were collected at 0, 24 h and 48 h for transcriptomic analysis. Five sampling groups of S. furcifera were flash-frozen in liquid nitrogen and stored at −80 °C in a refrigerator prior to RNA extraction. Total RNA was extracted from whole S. furcifera individuals using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). The purity, concentration and integrity of the RNA samples were tested using a NanoPhotometer® spectrophotometer (IMPLEN, Westlake Village, CA, USA) and agarose gels to ensure the use of high-quality samples for transcriptome sequencing. Then, a total amount of 1 μg of RNA per sample was used for RNA sample preparation. Sequencing libraries were generated using a NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA) following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample.
Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under an elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5×). First-strand cDNA was synthesized using random hexamer primers and M-MuLV reverse transcriptase. Second-strand cDNA synthesis was subsequently obtained using DNA polymerase I and RNase H. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of the 3′ ends of the DNA fragments, NEBNext adaptors with hairpin loop structures were ligated to prepare the samples for hybridization. To preferentially select cDNA fragments that were 240 bp in length, the library fragments were purified with the AMPure XP system (Beckman Coulter, Brea, CA, USA). Then, 3 μL of USER enzyme (NEB, Ipswich, MA, USA) was used with size-selected, adaptor-ligated cDNA at 37 °C for 15 min, followed by 5 min at 95 °C before PCR. Then, PCR was performed with Phusion High-Fidelity DNA polymerase, universal PCR primers and the index (X) primer. Finally, PCR products were purified by the AMPure XP system, and library quality was assessed on an Agilent Bioanalyzer 2100 system (Agilent, Santa Clara, CA, USA).
Clustering of the index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v4-cBot-HS (Illumina, San Diego, CA, USA), according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina platform and paired-end reads were generated.

2.4. Data Analysis

Raw data/raw reads in FASTQ format were first processed through in-house Perl scripts. In this step, clean data/clean reads were obtained by removing reads containing adapters, reads containing poly-N sequences and low-quality reads from the raw data. All downstream analyses were based on clean data with high quality. The adaptor sequences and low-quality sequence reads were removed from the data sets. At the same time, the Q20, Q30, GC content and sequence duplication level of the clean data were calculated. Raw sequences were transformed into clean reads after data processing. These clean reads were then mapped to the reference genome sequence [42]. Only reads with a perfect match or one mismatch were further analyzed and annotated based on the reference genome. Hisat2 software was used to map reads to the reference genome. The Illumina sequence reads have been submitted to the NCBI SRA database (accession No. PRJNA786731).
Gene functional annotation was based on the following databases: Nr (NCBI nonredundant protein sequences); Pfam (protein family); Nt (NCBI nonredundant nucleotide sequences); KO (KEGG orthologue database); SWISS-PROT (a manually annotated and reviewed protein sequence database); and GO (gene ontology). Gene expression levels were estimated by the fragments per kilobase of transcript per million fragments mapped (FPKM) values, with the following formula:
FPKM = cDNA Fragments/Mapped Fragments (Millions) * Transcript Length (kb)
Differential expression analysis was performed using DESeq2, which provided statistical analyses for determining differential expression in digital gene expression (DGE) data using a model based on the negative binomial distribution. The resulting p values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted p value < 0.01 found by DESeq2 were considered to be differentially expressed. Gene ontology (GO) enrichment analysis of the differentially expressed genes (DEGs) was implemented by the GOseq R package based on the Wallenius noncentral hypergeometric distribution, which adjusts for gene length bias in DEGs. KEGG is a database resource for understanding the high-level functions and utilities of biological systems, such as cells, organisms, and ecosystems, from molecular information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technology (http://www.genome.jp/kegg/, 17 November 2021). KOBAS software was then used to test the statistical enrichment of DEGs in KEGG pathways.

2.5. Validation of DEGs Library and Olfactory-Related Genes Identification

To validate the DEGs in the libraries, 20 DEGs (i.e., control vs. treatment) were randomly selected for a comparison using real-time quantitative PCR (qPCR). Moreover, the olfactory-related genes (e.g., ORs, OBPs and CSPs) in DEGs were also screened from the DEGs and analyzed by qPCR. The qPCR was performed on an iCycler iQ real-time PCR system (Bio-Rad, Hercules, CA, USA) with a QuantiNova SYBR Green PCR Kit (QIAGEN, Dusseldorf, Germany) according to the manufacturer’s instructions. The cycling parameters were as follows: initial denaturation at 94 °C for 30 s, followed by 40 cycles of 94 °C for 10 s, 55 °C for 30 s and 95°C for 1 s. The β-actin was selected as the reference gene for the normalization of the expression of the DEGs, according to the 2−ΔΔCt method. The primers designed for qPCR in this experiment are listed in Table S1.

3. Results

3.1. Effects of M. anisopliae on Feeding Choice of S. furcifera

We first checked the feeding behavior of S. furcifera on rice plants infected with SRBSDV or non-SRBSDV by random choice and Y-tube olfactory choice tests. In the random feeding choice test, the pest S. furcifera showed a feeding preference for the rice plant infected with SRBSDV during 24 h, except the time at 6 h (Figure 1A). Similarly, we also noted that this pest exhibited a feeding preference for the arm of the Y-tube connected with the SRBSDV-treated rice plant (Figure 1B, p < 0.001; t-test). However, when the pests were challenged by fungal M. anisopliae CQMa421, they showed no apparent feeding preference for any of the two rice plants (i.e., SRBSDV or non-SRBSDV rice plants), both in the random (except 3 h; see Figure 1C) and the Y-tube olfactory choice tests (Figure 1D).

3.2. Summary of Digital Gene Expression by Transcriptomic Profiling

To investigate potential targets affecting the feeding preference of the rice planthopper S. furcifera induced by the fungal M. anisopliae challenge, we further conducted the transcriptomic responses of S. furcifera to M. anisopliae inoculation after 0 h, 24 h and 48 h. A total of 15 DGE tag libraries were constructed, including the T-24 h vs. blank (0 h), T-24 h vs. CK-24 h, T-48 h vs. blank (0 h) and T-48 h vs. CK-48 h comparison for different treatments. By Illumina deep sequencing, 38.39~46.64 million raw reads were obtained from each treatment sample (Table 1). Prior to mapping, low-quality and adapter reads were filtered and 35.83~44.91 million clean-sequence reads per library were retained (Table 1). All samples had Q30 values greater than 92%, and the GC content ranged from 39.42% to 42.44% (Table 1), suggesting the reliability of the sequencing results.

3.3. Transcriptomic Comparison and Analysis of M. anisopliae Inoculation

To compare the DEGs among 15 different libraries, the gene expression levels were first determined from the FPKM values. Global analysis of the transcriptomic changes in S. furcifera at different time points after M. anisopliae inoculation demonstrated up- or downregulated genes in the control and treatment groups. DESeq.2 was selected to test DEGs with p < 0.01. These results showed that 1932 genes were upregulated at 24 h after fungal inoculation, and 1599 genes were downregulated by over two folds (|log2 (FoldChange)| > 2) in T-24 h vs. 0 h (Figure 2A). These results showed that 1380 genes were upregulated at 24 h after fungal inoculation, and 2187 genes were downregulated by over two folds (|log2 (FoldChange)| > 2) in T-24 h vs. CK-24 h (Figure 2B). At 48 h, there were 565 upregulated and 1345 downregulated genes (Figure 2C) in T-48 h vs. 0 h. By contrast, there were 137 upregulated and 106 downregulated genes after 48 h inoculation in T-48 h vs. CK-48 h (Figure 2D). These results showed that many different genes would be involved in the S. furcifera response to M. anisopliae CQMa421 infection over time (Figure 2E). We also found 1051, 44, 929 and 54 DEGs commonly expressed for T-24 h vs. CK-24 h and T-24 h vs. 0 h, for T-24 h vs. CK-24 h and T-48 h vs. CK-48 h, for T-24 h vs. 0 h and T-48 h vs. 0 h and for T-48 h vs. CK-48 h and T-48 h vs. CK-0 h, respectively (Figure S1). Specifically, the common 44 DEGs of the 24 h and 48 h treatment vs. control were majorly involved in RNA-directed DNA polymerase (g17653, g22334, g36185 and g48810), fatty acyl-CoA reductase (g82369), sphingosine kinase (g36165), glycerol-3-phosphate dehydrogenase (g14910) and malate dehydrogenase (g86484) (Table S2).

3.4. Functional Classification and Pathway Analysis

To examine the functions of the DEGs after M. anisopliae inoculation, we used the GO databases to map them. Specifically, after 24 h inoculation, the biological function category was assigned to 2008 DEGs (Figure 3A) compared with 0 h, and the biological function category was assigned to 2084 DEGs with the control after 24 h (Figure 3B). By contrast, the biological function category was assigned to 1107 DEGs (Figure 3C) after 48 h, compared with 0 h. However, the biological function category included 101 DEGs after 48 h with the control treatment, indicating different physiological responses of host S. furcifera to fungal infection over time (Figure 3D). Totally, the annotated GO terms included 5300 DEGs (2008, 2084, 1107 and 101 DEGs) in the BLAST database within the categories of biological process, cell component and molecular function (Figure 3A–D). For different periods, it also showed different abundance among the categories, and most of the DEGs were enriched in energy metabolism and antimicrobial activity during this period (Figure 3B,D).
In this study, we also selected KEGG to identify the metabolic and signal transduction pathways associated with these DEGs. The DEGs were mapped to 133 pathways in the KEGG database between the blank control (0 h) and treatment groups at 24 h inoculation (Figure 4A), and 110 pathways were mapped in the KEGG database between the control and treatment groups at 24 h post-inoculation (Figure 4B). By contrast, the DEGs between the blank control (0 h) and treatment groups at 48 h post-inoculation were mapped to 107 pathways in the KEGG database (Figure 4C), and the DEGs between the control and treatment groups at 48 h post-inoculation were mapped to 25 pathways in the KEGG database (Figure 4D). The top 20 pathways of each group in the richness analysis were exhibited in Figure 4, majorly including the pathways of carbon metabolism at 24 h vs. 0 h and the pathways of ribosome and its biogenesis at 24 h vs. control. Meanwhile, we noted that many pathways were involved in the peroxisome, amino sugar and nucleotide sugar metabolism at 48 h vs. 0 h, and in apoptosis and the MAPK signaling pathway at 48 h vs. control. These pathways are important for maintaining the physiological functions of the host insect and for enhancing its response to pathogenic invasion.
From the top 20 highly or lowly DEGs in host S. furcifera after M. anisopliae infection, the identified DEGs by the Nr annotation were mostly included in the cuticular proteins (g53545 and g51851), the nose resistant to fluoxetine protein 6 (g39720) and the elongation of the very long chain fatty acids protein (g80315). The upregulated targets were involved in retrovirus-related Pol polyprotein (g39683) and RNA-directed DNA polymerase from mobile element jockey (g36211 and g23323). After 24 h inoculation, the vitellogenin (g89906, g89903 and g89900) and serine protease (g89852) were downregulated, but the RNA-directed DNA polymerase (g31405 and g87480) was upregulated. We also noted that some of the targets, including the cardioaccelerator peptide receptor (g56305) and the nucleic-acid-binding protein (g48675) were downregulated after 48 h inoculation. In contrast, the putative RNA-directed DNA polymerase (g57847), adult-specific cuticular protein (g71543) and dopamine receptor (g88013) were upregulated during this period (Table S3).

3.5. Validation of DEGs Using RT-qPCR and Olfactory-Related Genes Identification

To further validate the RNAs identified through sequencing, we used qPCR to evaluate the results of RNA sequencing. We selected 20 of the highly or lowly expressed DEGs based on the Illumina sequencing results. The results showed that the expression trends of the selected RNAs showed a slight discrepancy from the findings of the sequencing analysis (Figure 5A), which might be due to the differences in the sensitivity, specificity and applied algorithms between the two techniques.
The expression level of olfactory-related genes plays key roles in regulating the feeding behavior of insects. Thus, we further examined the potential genes that might affect the feeding preference of S. furcifera after M. anisopliae infection. We screened a few olfactory-related genes, including OPBs, ORs, IRs and CSPs, from the DEGs after M. anisopliae inoculation (Figure 5B–E). Interestingly, we noted that the OR_23 and OR_43 were downregulated at 48 h post-inoculation, which may be potential target genes affecting the feeding behavior of S. furcifera.

4. Discussion

Damage by insect pests, including rice planthoppers, Cnaphalocrocis medinalis and Chilo suppressalis, are one of major factors resulting in losses of rice yields [2,43]. The population of rice planthoppers has evolved a high resistance to many commonly used insecticides [44]. On the other hand, those pests are important vectors of rice plant viruses (e.g., the rice ragged stunt virus, the rice grassy stunt virus and the southern rice black-streaked dwarf virus) for indirectly affecting the growth of rice plants further [45]. Previous studies have reported that the olfactory behavior of insect vectors would affect the efficiency of virus transmission and dispersal [46,47]. Moreover, carrying plant viruses could regulate the feeding behavior of insect vectors [33,48]. Thus, avoiding or disturbing the feeding behavior of vectors is a potential strategy to suppress the plant virus. In this study, we noted that the feeding preference of S. furcifera for SRBSDV rice plants has changed after fungal M. anisopliae infection. A further transcriptomic analysis suggests that the infection of M. anisopliae on S. furcifera would affect the expression of many potential targets, including several olfactory-related genes.
The rice virus SRBSDV is a serious plant virus, which could result in considerable yield losses of rice [6]. After its first report in Goungdong, China in 2008, the virus has spread into Japan and Vietnam, and has become one of the most important plant viruses for rice [37]. The rice planthopper S. furcifera is the unique vector of SRBSDV with a long-life transmitting ability. The occurrence of SRBSDV is related to the breakout of the S. furcifera population under field conditions [7]. Interestingly, this pest showed a feeding preference for those rice plants infected with SRBSDV, due to some specific odors [39]. Another insect pest, the whitefly, also showed a steady preference for Tomato yellow leaf curl virus-infected plants [49]. The preference of vectors would enhance the probability of the virus to enter uninfected hosts, and eventually benefit the virus’ spread among the plant community [49,50]. However, we found that the feeding preference of S. furcifera for certain target plants disappeared at 48 h post-inoculation.
Transcriptome profiling is a beneficial method to identify the potential genes that regulate physiology, growth and behavior. In this study, to understand how the M. anisopliae regulated the feeding behavior of S. furcifera, we identified the DEGs after fungal inoculation at 0 h, 24 h and 48 h. During fungal infection, a series of metabolic pathways are screened, including amide and peptide biosynthesis, carbohydrate catabolic processes and immune responses (GO0002429, GO0002757 and GO0002764, respectively). In fact, the recognition of hosts by pathogens occurred in the initial inoculation stage within 4 h. We also noted the cell surface receptor signaling pathway (GO0002768) in S. furcirera after 24 h inoculation. The recognition molecules, such as peptidoglycan recognition proteins (PGRPs), β-1,3- glucan recognition proteins (βGRPs), galectins, C-type lectins (CTLs) and scavenger receptors (SCRs) [24,27,51], play vital roles during this period. However, we did not investigate the initial stages after this fungal inoculation in this study, investigating the stages after 24 h instead. The high-expression genes are majorly involved in RNA-directed DNA polymerase (g36211, g31405 and g87480) during this period. In contrast, the suppressed metabolisms are involved in serine protease (g89852) and farnesol dehydrogenase (g28615). Moreover, the MAPK signaling pathway (dnx04013), apoptosis (dnx04214), Toll and IMD signaling pathway (dnx04624) and FoxO signaling pathway (dnx04068) were activated at 48 h post-inoculation.
Importantly, we identified several DEGs of olfactory-related genes after the fungal M. anisopliae inoculation under different periods. The olfactory behaviors of insects are important clues for targeting foods, mating and avoiding natural enemies [28,29]. The oviposition behavior was disturbed after the Or35 deficiency in the tobacco hawkmoths Manduca sexta [31] and the Or31 deficiency in Helicoverpa assulta [52]. In this study, we noted that the feeding preferences have been disturbed, from preference to non-preference for SRBSDV rice plants, after M. anisopliae infection. This phenomenon indicated that M. anisopliae might affect the olfactory capability by interacting with these targets. Indeed, we found several differentially expressed OBPs and ORs after M. anisopliae infection. An altered expression of CSPs and OBPs in response to fungal infection has also been reported in the red fire ant Solenopsis invicta [53]. Those genes play important roles in regulating the host behavior, especially feeding behaviors. Moreover, feeding behaviors would affect the transmission of rice viruses due to their dependency on the insect vector [33,54]. The variation of feeding preferences to targeting rice plants would have a further impact on the dispersal of SRBSDV from plant to plant. Thus, we may further develop strategies (e.g., RNAi to target OR OBP) to regulate the efficiency of vector S. furcifera transmitting this virus.
In conclusion, we noted that the infection of M. anisopliae could disturb the feeding behavior of host S. furcifera and further identified several olfactory-related genes (i.e., ORs and OBPs) by transcriptomic profiling. These results not only indicate that the fungal M. anisopliae can regulate host behavior, but also suggest such regulations may be involved in the expression level of olfactory-related genes. Future work will warrant the specific functions of several ORs and OBPs, and develop effective approaches for the control of this pest and its transmitted SRBSDV.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jof8050506/s1, Figure S1: The Venn diagram of DEGs after the fungal infection. Table S1: The primers used for RT-qPCR in this study. Table S2: The descriptions of commonly expressed DEGs in Venn diagram after M. anisopliae infection. Table S3: Top 20 highly or lowly DEGs in different period of M. anisopliae infection.

Author Contributions

Conceptualization, J.X. and Y.X.; methodology, Y.W. and L.H.; project administration, J.X. and Y.X.; investigation, Y.W. and L.H.; writing—review and editing, J.X., Y.X., Y.W. and L.H.; funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the grants from the Natural Science Foundation of Chongqing (Project No. cstc2021jcyj-msxmX0561), the Program for Chongqing Talents of China (Project No. cstc2021ycjh-bgzxm0313) and the Fundamental Research Funds for the Central Universities (Chongqing University) with Project No. 2021CDJQY-049.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequence data are available in NCBI GenBank following the accession number No. PRJNA786731 in the manuscript.

Acknowledgments

We thank Deyu Zeng and Chunfang Huang for assistance with rearing the insects in this experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Muthayya, S.; Sugimoto, J.D.; Montgomery, S.; Maberly, G.F. An overview of global rice production, supply, trade, and consumption. Ann. N. Y. Acad. Sci. 2014, 1324, 7–14. [Google Scholar] [CrossRef] [PubMed]
  2. Ma, K.C.; Lee, S.C. Occurrence of major rice insect pests at different transplanting times and fertilizer levels in paddy field. Korean J. Appl. Entomol. 1996, 35, 132–136. [Google Scholar]
  3. Wei, T.Y.; Li, Y. Rice Reoviruses in insect vectors. Annu. Rev. Phytopathol. 2016, 54, 99–120. [Google Scholar] [CrossRef] [PubMed]
  4. Lou, Y.G.; Zhang, G.R.; Zhang, W.Q.; Hu, Y.; Zhang, J. Biological control of rice insect pests in China. Biol. Control 2013, 67, 8–20. [Google Scholar] [CrossRef]
  5. Cheng, J. Rice planthopper problems and relevant causes in China. In Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia; International Rice Research Institute: Los Baños, Philippines, 2009; Chapter 1; pp. 157–178. [Google Scholar]
  6. Zhou, G.H.; Xu, D.L.; Xu, D.G.; Zhang, M.X. Southern rice black-streaked dwarf virus: A white-backed planthopper-transmitted fijivirus threatening rice production in Asia. Front. Microbiol. 2013, 4, 270. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, Y.; Zhang, G.; Chen, X.; Li, X.J.; Xiong, K.; Cao, S.P.; Hu, Y.Y.; Lu, M.H.; Liu, W.C.; Tuan, H.A.; et al. The influence of Sogatella furcifera (Hemiptera: Delphacidae) migratory events on the southern rice black-streaked dwarf virus epidemics. J. Econ. Entomol. 2017, 110, 854–864. [Google Scholar] [CrossRef]
  8. Savary, S.; Willocquet, L.; Pethybridge, S.J.; Esker, P.; McRoberts, N.; Nelson, A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 2019, 3, 430–439. [Google Scholar] [CrossRef]
  9. Hamiduzzaman, M.M.; Sinia, A.; Guzman-Novoa, E.; Goodwin, P.H. Entomopathogenic fungi as potential biocontrol agents of the ecto-parasitic mite, Varroa destructor, and their effect on the immune response of honey bees (Apis mellifera L.). J. Invertebr. Pathol. 2012, 111, 237–243. [Google Scholar] [CrossRef]
  10. Sánchezpeña, S.R.; Lara, J.S.; Medina, R.F. Occurrence of entomopathogenic fungi from agricultural and natural ecosystems in Saltillo, México, and their virulence towards thrips and whiteflies. J. Insect Sci. 2011, 11, 1. [Google Scholar] [CrossRef]
  11. Bordalo, M.D.; Gravato, C.; Beleza, S.; Campos, D.; Lopes, I.; Pestana, J.L.T. Lethal and sublethal toxicity assessment of Bacillus thuringiensis var. israelensis and Beauveria bassiana based bioinsecticides to the aquatic insect Chironomus riparius. Sci. Total Environ. 2020, 698, 134155. [Google Scholar] [CrossRef]
  12. Jiang, W.; Peng, Y.; Ye, J.; Wen, Y.; Liu, G.; Xie, J. Effects of the entomopathogenic fungus Metarhizium anisopliae on the mortality and immune response of Locusta migratoria. Insects 2019, 11, 36. [Google Scholar] [CrossRef] [Green Version]
  13. Fite, T.; Tefera, T.; Negeri, M.; Damte, T.; Sori, W. Evaluation of Beauveria bassiana, Metarhizium anisopliae, and Bacillus thuringiensis for the management of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) under laboratory and field conditions. Biocontrol Sci. Technol. 2019, 30, 278–295. [Google Scholar] [CrossRef]
  14. Rice, S.J.; Baker, D.K.; Leemon, D.M. Development of mycoinsecticide formulations with Beauveria bassiana and Metarhizium anisopliae for the control of lesser mealworm, Alphitobius diaperinus, in chicken broiler houses. Biocontrol 2019, 64, 489–500. [Google Scholar] [CrossRef]
  15. Peng, G.X.; Xie, J.Q.; Guo, R.; Keyhani, N.O.; Zeng, D.Y.; Yang, P.Y.; Xia, Y.X. Long-term field evaluation and large-scale application of a Metarhizium anisopliae strain for controlling major rice pests. J. Pest Sci. 2021, 94, 969–980. [Google Scholar] [CrossRef]
  16. Tang, J.; Liu, X.; Ding, Y.; Jiang, W.; Xie, J. Evaluation of Metarhizium anisopliae for rice planthopper control and its synergy with selected insecticides. Crop Prot. 2019, 121, 132–138. [Google Scholar] [CrossRef]
  17. Pan, C.Y.; Cai, Y.J.; Li, T.C.; Zhang, W.Q. Control effect of Metarhizium anisopliae and its mixture with dsRNA on the brown planthopper. J. Environ. Entomol. 2016, 38, 1071–1077. [Google Scholar]
  18. Fang, W.; Vega-Rodriguez, J.; Ghosh, A.K.; Jacobs-Lorena, M.; Kang, A.; St. Leger, R.J. Development of transgenic fungi that kill human malaria parasites in mosquitoes. Science 2011, 331, 1074–1077. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, C.; Wang, S. Insect pathogenic fungi: Genomics, molecular interactions, and genetic improvements. Annu. Rev. Entomol. 2017, 62, 73. [Google Scholar] [CrossRef]
  20. Duressa, T.F.; Vanlaer, R.; Huybrechts, R. Locust cellular defense against inoculations: Sites of pathogen clearance and hemocyte proliferation. Dev. Comp. Immunol. 2015, 48, 244–253. [Google Scholar] [CrossRef]
  21. Olombrada, M.; Lazaro-Gorines, R.; Lopez-Rodriguez, J.C.; Martinez-Del-Pozo, A.; Onaderra, M.; Maestro-Lopez, M.; Lacadena, J.; Gavilanes, J.G.; Garcia-Ortega, L. Fungal Ribotoxins: A review of potential biotechnological applications. Toxins 2017, 9, 71. [Google Scholar] [CrossRef] [Green Version]
  22. Shakeel, M.; Xu, X.; Xu, J.; Li, S.; Yu, J.; Zhou, X.; Xu, X.; Hu, Q.; Yu, X.; Jin, F. Genome-Wide identification of destruxin A-responsive immunity-related microRNAs in diamondback moth, Plutella xylostella. Front. Immunol. 2018, 9, 185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yu, Y.; Cao, Y.; Xia, Y.; Liu, F. Wright-Giemsa staining to observe phagocytes in Locusta migratoria infected with Metarhizium acridum. J. Invertebr. Pathol. 2016, 139, 19–24. [Google Scholar] [CrossRef] [PubMed]
  24. Lu, H.L.; St. Leger, R.J. Insect immunity to entomopathogenic fungi. Adv. Genet. 2016, 94, 251–285. [Google Scholar] [PubMed]
  25. Xie, J.; Peng, Y.; Xia, Y. Genome-Wide identification and analysis of Nilaparvata lugens microRNAs during challenge with the entomopathogenic fungus Metarhizium anisopliae. J. Fungi 2021, 7, 295. [Google Scholar] [CrossRef] [PubMed]
  26. Junges, A.; Boldo, J.T.; Souza, B.K.; Guedes, R.L.; Sbaraini, N.; Kmetzsch, L.; Thompson, C.E.; Staats, C.C.; de Almeida, L.G.; de Vasconcelos, A.T.; et al. Genomic analyses and transcriptional profiles of the glycoside hydrolase family 18 genes of the entomopathogenic fungus Metarhizium anisopliae. PLoS ONE 2014, 9, e107864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Xu, J.; Xu, X.X.; Shakeel, M.; Li, S.Z.; Wang, S.; Zhou, X.Q.; Yu, J.L.; Xu, X.J.; Yu, X.Q.; Jin, F.L. The Entomopathogenic fungi Isaria fumosorosea plays a vital role in suppressing the immune system of Plutella xylostella: RNA-Seq and DGE analysis of immunity-related genes. Front. Microbiol. 2017, 8, 1421. [Google Scholar] [CrossRef] [Green Version]
  28. Leal, W.S. Odorant reception in insects: Roles of receptors, binding proteins, and degrading enzymes. Annu. Rev. Entomol. 2013, 58, 373–391. [Google Scholar] [CrossRef]
  29. Pelosi, P.; Iovinella, I.; Zhu, J.; Wang, G.; Dani, F.R. Beyond chemoreception: Diverse tasks of soluble olfactory proteins in insects. Biol. Rev. 2018, 93, 184–200. [Google Scholar] [CrossRef] [Green Version]
  30. Del Marmol, J.; Yedlin, M.A.; Ruta, V. The structural basis of odorant recognition in insect olfactory receptors. Nature 2021, 597, 126–131. [Google Scholar] [CrossRef]
  31. Zhang, J.; Komail Raza, S.A.; Wei, Z.; Keesey, I.W.; Parker, A.L.; Feistel, F.; Chen, J.; Cassau, S.; Fandino, R.A.; Grosse-Wilde, E.; et al. Competing beetles attract egg laying in a hawkmoth. Curr. Biol. 2022, 32, 861–869. [Google Scholar] [CrossRef]
  32. Guo, X.J.; Yu, Q.Q.; Chen, D.F.; Wei, J.N.; Yang, P.C.; Yu, J.; Wang, X.H.; Kang, L. 4-Vinylanisole is an aggregation pheromone in locusts. Nature 2020, 584, 584–588. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, H.; He, X.; Zheng, X.; Yang, Y.; Tian, J.; Lu, Z. Southern rice black-streaked dwarf virus (SRBSDV) directly affects the feeding and reproduction behavior of its vector, Sogatella furcifera (Horvath) (Hemiptera: Delphacidae). Virol. J. 2014, 11, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wei, H.S.; Qin, J.H.; Cao, Y.Z.; Li, K.B.; Yin, J. Two classic OBPs modulate the responses of female Holotrichia oblita to three major ester host plant volatiles. Insect Mol. Biol. 2021, 30, 390–399. [Google Scholar] [CrossRef] [PubMed]
  35. He, P.; Chen, G.L.; Li, S.; Wang, J.; Ma, Y.F.; Pan, Y.F.; He, M. Evolution and functional analysis of odorant-binding proteins in three rice planthoppers: Nilaparvata lugens, Sogatella furcifera, and Laodelphax striatellus. Pest Manag. Sci. 2019, 75, 1606–1620. [Google Scholar] [CrossRef]
  36. Zhan, H.; Dewer, Y.; Zhang, J.; Tian, J.; Li, D.; Qu, C.; Yang, Z.; Li, F.; Luo, C. Odorant-binding protein 1 plays a crucial role in the olfactory response of Bemisia tabaci to R-curcumene. J. Agric. Food Chem. 2021, 69, 12785–12793. [Google Scholar] [CrossRef]
  37. Zhou, G.H.; Wen, J.J.; Cai, D.J.; Li, P.; Xu, D.L.; Zhang, S.G. Southern rice black-streaked dwarf virus: A new proposed Fijivirus species in the family Reoviridae. Chin. Sci. Bull. 2008, 53, 3677–3685. [Google Scholar] [CrossRef]
  38. Peng, Y.; Tang, J.; Hong, M.; Xie, J. Suppression of rice planthopper populations without affecting the rice microbiota by the entomopathogenic fungus Metarhizium anisopliae. Appl. Environ. Microbiol. 2020, 86, e01337-20. [Google Scholar] [CrossRef]
  39. Xu, Y.; Zhou, W.; Zhou, Y.; Wu, J.; Zhou, X. Transcriptome and comparative gene expression analysis of Sogatella furcifera (Horvath) in response to southern rice black-streaked dwarf virus. PLoS ONE 2012, 7, e36238. [Google Scholar] [CrossRef]
  40. Lin, J.; He, J.; Liang, A.; Wang, F. Transcriptome profiling and dimorphic expression of sex-related genes in fifth-instar nymphs of Sogatella furcifera, an important rice pest. Genomics 2020, 112, 1105–1111. [Google Scholar] [CrossRef]
  41. Yang, S.; Xu, R.; Yang, S.-Y.; Kuang, R.-P. Olfactory responses of Aphidius gifuensis to odors of host plants and aphid-plant complexes. Insect Sci. 2009, 16, 503–510. [Google Scholar] [CrossRef]
  42. Wang, L.; Tang, N.; Gao, X.; Chang, Z.; Zhang, L.; Zhou, G.; Guo, D.; Zeng, Z.; Li, W.; Akinyemi, I.A.; et al. Genome sequence of a rice pest, the white-backed planthopper (Sogatella furcifera). Gigascience 2017, 6, giw004. [Google Scholar]
  43. Cheng, J. Rice Planthoppers in the Past Half Century in China. In Rice Planthoppers; Springer: Dordrecht, The Netherlands, 2015; pp. 1–32. [Google Scholar]
  44. Liao, X.; Xu, P.F.; Gong, P.P.; Wan, H.; Li, J.H. Current susceptibilities of brown planthopper Nilaparvata lugens to triflumezopyrim and other frequently used insecticides in China. Insect Sci. 2021, 28, 115–126. [Google Scholar] [CrossRef] [PubMed]
  45. Hibino, H. Biology and epidemiology of rice viruses. Annu. Rev. Phytopathol. 1996, 34, 249–274. [Google Scholar] [CrossRef] [PubMed]
  46. Hu, Z.; Chai, R.; Liu, X.; Dong, Y.; Su, D.; Desneux, N.; Tan, X.; Luo, C. Barley yellow dwarf virus-infected wheat plant modulated selection behavior of vector aphids. J. Pest Sci. 2021, 95, 1273–1285. [Google Scholar] [CrossRef]
  47. Mauck, K.E.; De Moraes, C.M.; Mescher, M.C. Deceptive chemical signals induced by a plant virus attract insect vectors to inferior hosts. Proc. Natl. Acad. Sci. USA 2010, 107, 3600–3605. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, S.; Guo, H.; Zhu-Salzman, K.; Ge, F.; Sun, Y. PEBP balances apoptosis and autophagy in whitefly upon arbovirus infection. Nat. Commun. 2022, 13, 846. [Google Scholar] [CrossRef]
  49. Wang, S.; Guo, H.; Ge, F.; Sun, Y. Apoptotic neurodegeneration in whitefly promotes the spread of TYLCV. Elife 2020, 9, e56168. [Google Scholar] [CrossRef]
  50. Ingwell, L.L.; Eigenbrode, S.D.; Bosque-Perez, N.A. Plant viruses alter insect behavior to enhance their spread. Sci. Rep. 2012, 2, 578. [Google Scholar] [CrossRef] [Green Version]
  51. Dan, H. Drosophila immunity: Paths and patterns. Curr. Opin. Immunol. 2003, 15, 12–19. [Google Scholar]
  52. Li, R.T.; Huang, L.Q.; Dong, J.F.; Wang, C.Z. A moth odorant receptor highly expressed in the ovipositor is involved in detecting host-plant volatiles. Elife 2020, 9, e53706. [Google Scholar] [CrossRef]
  53. Wei, Z.; Ortiz-Urquiza, A.; Keyhani, N.O. Altered expression of chemosensory and odorant binding proteins in response to fungal infection in the red imported rire ant, Solenopsis invicta. Front. Physiol. 2021, 12, 596571. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, J.; Zheng, X.; Chen, Y.; Hu, J.; Dong, J.; Su, X.; Zhang, Z. Southern rice black-streaked dwarf virus infection improves host suitability for its insect vector, Sogatella furcifera (Hemiptera: Delphacidae). J. Econ. Entomol. 2014, 107, 92–97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jof 08 00506 g001 550
Figure 1. The effects of M. anisopliae on host S. furcifera feeding behavior. (A) Feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under random feeding test; (B) feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under Y-tube olfactory test; (C) after 48 h infection, the feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under random feeding test; (D) after 48 h infection, the feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under the Y-tube olfactory test. The asterisks indicate the significant difference. “*” indicates p < 0.05; “**” indicates p < 0.01 and “***” indicates p < 0.001.
Figure 1. The effects of M. anisopliae on host S. furcifera feeding behavior. (A) Feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under random feeding test; (B) feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under Y-tube olfactory test; (C) after 48 h infection, the feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under random feeding test; (D) after 48 h infection, the feeding choice of S. furcifera for SRBSDV or non-SRBSDV rice plants under the Y-tube olfactory test. The asterisks indicate the significant difference. “*” indicates p < 0.05; “**” indicates p < 0.01 and “***” indicates p < 0.001.
Jof 08 00506 g001
Jof 08 00506 g002 550
Figure 2. The volcano plot and numbers of differentially expressed genes (DEGs) of S. furcifera identified after M. anisopliae inoculation. (A) The volcano plot of DEGs at 24 h post inoculation vs. 0 h; (B) the volcano plot of DEGs at 24 h post inoculation vs. 24 h control; (C) the volcano plot of DEGs at 48 h post inoculation vs. 0 h; (D) the volcano plot of DEGs at 48 h post inoculation vs. 48 h control; and (E) the number of S. furcifera DEGs during M. anisopliae infection from 0 to 48 h.
Figure 2. The volcano plot and numbers of differentially expressed genes (DEGs) of S. furcifera identified after M. anisopliae inoculation. (A) The volcano plot of DEGs at 24 h post inoculation vs. 0 h; (B) the volcano plot of DEGs at 24 h post inoculation vs. 24 h control; (C) the volcano plot of DEGs at 48 h post inoculation vs. 0 h; (D) the volcano plot of DEGs at 48 h post inoculation vs. 48 h control; and (E) the number of S. furcifera DEGs during M. anisopliae infection from 0 to 48 h.
Jof 08 00506 g002
Jof 08 00506 g003 550
Figure 3. Gene ontology (GO) categories of S. furcifera DEGs after challenge with M. anisopliae at different periods. (A) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 24 h vs. 0 h; (B) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 24 h; (C) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 48 h vs. 0 h; and (D) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 48 h. The asterisks (*) indicate the significant enrichment of GO terms.
Figure 3. Gene ontology (GO) categories of S. furcifera DEGs after challenge with M. anisopliae at different periods. (A) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 24 h vs. 0 h; (B) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 24 h; (C) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 48 h vs. 0 h; and (D) GO categories of S. furcifera DEGs after challenge with M. anisopliae at 48 h. The asterisks (*) indicate the significant enrichment of GO terms.
Jof 08 00506 g003
Jof 08 00506 g004 550
Figure 4. Enrichment and dispersion point map of differentially expressed genes (DEGs) in KEGG pathways. (A) The DEGs involved in the pathways after challenge for 24 h vs. 0 h; (B) the DEGs involved in the pathways after challenge for 24 h; (C) the DEGs involved in the pathways after challenge for 48 h vs. 0 h; (D) the DEGs involved in the pathways after challenge for 48 h. The circles in the graph indicate that the KEGG pathway with the gene number and enrichment factor (q-value, different colors) are displayed on the y and x axes, respectively.
Figure 4. Enrichment and dispersion point map of differentially expressed genes (DEGs) in KEGG pathways. (A) The DEGs involved in the pathways after challenge for 24 h vs. 0 h; (B) the DEGs involved in the pathways after challenge for 24 h; (C) the DEGs involved in the pathways after challenge for 48 h vs. 0 h; (D) the DEGs involved in the pathways after challenge for 48 h. The circles in the graph indicate that the KEGG pathway with the gene number and enrichment factor (q-value, different colors) are displayed on the y and x axes, respectively.
Jof 08 00506 g004
Jof 08 00506 g005 550
Figure 5. Transcriptomic validation by qPCR and the olfactory-related genes identification. (A) A comparison of qPCR and RNA-seq of 20 genes; (B) the DEGs involved in the regulation of host olfactory behavior at 24 h vs. 0 h; (C) the DEGs involved in the regulation of host olfactory behavior at 24 h; (D) the DEGs involved in the regulation of host olfactory behavior at 48 h vs. 0 h; and (E), the DEGs involved in the regulation of host olfactory behavior at 48 h.
Figure 5. Transcriptomic validation by qPCR and the olfactory-related genes identification. (A) A comparison of qPCR and RNA-seq of 20 genes; (B) the DEGs involved in the regulation of host olfactory behavior at 24 h vs. 0 h; (C) the DEGs involved in the regulation of host olfactory behavior at 24 h; (D) the DEGs involved in the regulation of host olfactory behavior at 48 h vs. 0 h; and (E), the DEGs involved in the regulation of host olfactory behavior at 48 h.
Jof 08 00506 g005
Table
Table 1. Statistical summary of S. furcifera groups, after having inoculated M. anisopliae or not.
Table 1. Statistical summary of S. furcifera groups, after having inoculated M. anisopliae or not.
Sample NameRaw ReadsRaw BasesClean ReadsClean BasesQ30
T_24 h_138,699,9405.80G35,830,0305.37G92.91%
T_24 h_242,519,4766.37G40,606,9486.09G93.49%
T_24 h_338,627,1485.79G36,601,0085.49G93.25%
CK_24 h_141,683,1066.25G39,847,9285.98G93.70%
CK_24 h_244,749,8026.71G42,760,6666.41G93.92%
CK_24 h_346,649,4306.99G44,910,8186.74G93.08%
T_48 h_140,096,2626.01G38,405,3765.76G93.02%
T_48 h_242,917,0826.43G41,055,7746.16G93.73%
T_48 h_340,289,4446.04G38,627,4785.79G93.26%
CK_48 h_143,887,0286.58G42,403,9146.36G92.82%
CK_48 h_238,657,4065.79G36,866,6225.53G93.27%
CK_48 h_340,229,5306.03G38,382,2425.76G93.54%
Blank_141,866,6146.27G40,000,4446.00G93.20%
Blank_238,396,5325.75G36,692,4185.50G94.56%
Blank_343,095,9046.46G41,381,5686.21G94.07%
T, M. anisopliae treatment group; CK, control group.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Han, L.; Xia, Y.; Xie, J. The Entomopathogenic Fungus Metarhizium anisopliae Affects Feeding Preference of Sogatella furcifera and Its Potential Targets’ Identification. J. Fungi 2022, 8, 506. https://doi.org/10.3390/jof8050506

AMA Style

Wang Y, Han L, Xia Y, Xie J. The Entomopathogenic Fungus Metarhizium anisopliae Affects Feeding Preference of Sogatella furcifera and Its Potential Targets’ Identification. Journal of Fungi. 2022; 8(5):506. https://doi.org/10.3390/jof8050506

Chicago/Turabian Style

Wang, Yirong, Lijuan Han, Yuxian Xia, and Jiaqin Xie. 2022. "The Entomopathogenic Fungus Metarhizium anisopliae Affects Feeding Preference of Sogatella furcifera and Its Potential Targets’ Identification" Journal of Fungi 8, no. 5: 506. https://doi.org/10.3390/jof8050506

APA Style

Wang, Y., Han, L., Xia, Y., & Xie, J. (2022). The Entomopathogenic Fungus Metarhizium anisopliae Affects Feeding Preference of Sogatella furcifera and Its Potential Targets’ Identification. Journal of Fungi, 8(5), 506. https://doi.org/10.3390/jof8050506

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop