Skip to main content

Comparative transcriptome analysis reveals differences in gene expression in whitefly following individual or combined applications of Akanthomyces attenuatus (Zare & Gams) and matrine

Abstract

Background

Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) is a serious pest of crops in different regions of the world. Our recent studies on the joint application of Akanthomyces attenuatus (a pathogenic insect fungus) and matrine (a botanical insecticide) against B. tabaci have shown promising results. Using RNA sequencing (RNA-Seq), we identified differentially expressed genes involved in whitefly responses to single or mixed applications of A. attenuatus and matrine.

Methods

In this study, we compared the transcriptome profiles of B. tabaci treated with individual and combined treatments of A. attenuatus and matrine to determine variations in gene expression among whiteflies in response to different treatments.

Results

Transcriptomic data analysis showed differential expression of 71, 1194, and 51 genes in response to A. attenuatus (BtA), matrine (BtM), and A. attenuatus + matrine (BtAM) treatment, respectively. A total of 65 common differentially expressed genes (DEGs) were identified between whiteflies treated with A. attenuatus (BtA) and matrine (BtM). A comparison of DEGs across the three treatments (BtA, BtM, and BtAM) revealed two common DEGs. The results also revealed that AMPK signaling, apoptosis, and drug metabolism pathways are likely involved in whitefly defense responses against A. attenuatus and matrine infection. Furthermore, a notable suppression of general metabolism and immune response genes was observed in whiteflies treated with A. attenuatus + matrine (BtAM) compared to whiteflies treated with individual A. attenuatus (BtA) or matrine (BtM) treatments.

Conclusion

Dynamic changes in the number of differentially expressed genes were observed in B. tabaci subjected to different treatments (BtA, BtM, and BtAM). To the best of our knowledge, this is the first report on the molecular interactions between whitefly and individual or combined treatments of A. attenuatus and matrine. These results will further improve our knowledge of the infection mechanism and complex biochemical processes involved in the synergistic action of A. attenuatus and matrine against B. tabaci.

Peer Review reports

Background

The whitefly, Bemisia tabaci, is an invasive pest and disease vector with a rapidly increasing global population of over 50 cryptic species known to date [1]. Over the course of the past two decades, the cryptic species B. tabaci Middle East-Asia Minor 1 (MEAM1) (formerly known as 'Biotype B') replaced the indigenous B. tabaci populations in China [2]. Bemisia tabaci MEAM1 can cause direct (via feeding on phloem sap) or indirect damage (by secreting honeydew) to its host plants [2, 3]. B. tabaci MEAM1 adults have also been found to transmit over 150 plant viruses to commercial crops [4]. The use of synthetic chemicals for the B. tabaci management has resulted in development of insecticide resistance [5]. Furthermore, the heavy application of chemical insecticides has been affecting the natural ecosystem (including air, non-target organisms, and humans as well) which has caused a great awareness about finding sustainable alternates to the synthetic chemical [6]. Recently, the joint use of insect pathogenic fungi and botanical insecticides has emerged as a possible alternative to chemical insecticides [7]. Over the last century, entomopathogenic fungi have been used for the biological control of insect pests [7]. Entomopathogenic fungi are more appealing than other pest control strategies owing to their diversity, metabolite production, and biosafety levels [3]. Among these fungi, Akanthomyces attenuatus Zare & Gams is a well-known pathogen of insect pests, including whiteflies, aphids, and thrips [8,9,10,11,12]. However, the field application of A. attenuatus and other entomopathogenic fungi has received little attention owing to their slow action [13]. Recently, many studies have investigated the joint application of entomopathogenic fungi with natural enemies, synthetic oils, bacterial metabolite avermectins, and other biopesticides [14,15,16,17]. Our recent studies on the combined application of matrine (a botanical insecticide) and entomopathogenic fungi (A. attenuatus and Beauveria brongniartii) against different insect pests have shown promising results [7, 17, 18].

Innate immunity is the first line of protection against numerous pathogenic illnesses. Conserved signal cascades within the innate immune system against insecticides and natural pathogens have been well documented in whiteflies [19,20,21]. However, to the best of our knowledge, no detailed investigations have been conducted on the mechanisms responsible for controlling the whitefly immune system when defending against the individual or combined applications of A. attenuatus and matrine. The recently available whitefly transcriptome sequences [22, 23], in combination with RNA sequencing (RNA-Seq) technology, a revolutionary tool for measuring levels of gene expression [24, 25], have provided unprecedented opportunities to investigate the transcriptional response of B. tabaci to single or mixed applications of A. attenuatus and matrine. Using RNA-Seq, we identified differentially expressed genes in B. tabaci nymphs treated with single or mixed applications of A. attenuatus and matrine and analyzed how the whitefly orchestrates its defense responses to these applications. Altogether, our work provides the first report to reveal the molecular mechanisms of B. tabaci responses to single or mixed applications of A. attenuatus and matrine and new insights into the whitefly infection process during the synergistic action of the interaction between A. attenuatus and matrine.

Results

Insect bioassays

The mortality of the 2nd instar B. tabaci nymphs treated with A. attenuatus and matrine alone or in combination was significantly higher than that of the control. After 3 d of application, B. tabaci mortality in response to the combined treatment with A. attenuatus and matrine increased substantially compared to the individual treatments. After 5 d of treatment, the mortalities of 2nd instar B. tabaci nymphs treated with A. attenuatus (1×106 conidia/mL), matrine (1.0 mg/L), and A. attenuatus + matrine (1×106 conidia/mL + 1.0 mg/L) were 65.67%, 43.67%, and 100%, respectively (Fig. 1).

Fig. 1
figure 1

Concentration mortality response of Bemisia tabaci to individual or joint treatments of Akanthomyces attenuatus and matrine. Error bars indicate standard error of means based on three replicates

Summary of RNA-seq datasets

A total of 659,350,532 raw reads were generated from 12 samples (Additional File 1A). The N reads, low-quality reads, and reads with adapter sequences were removed to ensure the quality and accuracy of sequencing data, and 650,606,120 clean reads were generated that were successfully mapped to 73.3-79.19% of the reference genome; subsequent analysis was based on clean reads. Contigs were obtained first, followed by transcripts, resulting in 35,895 unigenes, of which 23,694 (66.06%) were longer than 1,000 bp. According to the length distribution of B. tabaci unigenes, the largest gene length was >1,800 bp (17,742 genes; 49.47 %), the second largest gene length was 201–400 bp (4,835 genes, 13.48 %), and the smallest gene length was 1,601–1,800 bp (1,309 genes, 3.64%) (Additional file 1B). The Pearson’s correlation coefficients showed highly reproducible data across different replicates (Additional file 1C).

The unigenes were compared using six databases (NR, Swiss-Prot, Pfam, Eggnog, gene ontology [GO], and Kyoto Encyclopedia of Genes and Genomes [KEGG]). A total of 16,016 (0.73%) unigenes were annotated in this study through GO database matching (7,199; 33%), followed by KEGG (7,510; 34%), COG (13,165; 60%), NR (15,910; 73%), Swiss-Prot (10,846; 50%), and Pfam (11,880; 54%) databases (Fig. 2).

Fig. 2
figure 2

Functional annotation of genes to different databases

Overview of B. tabaci differentially expressed genes (DEGs) in response to different treatments

Our results showed differential expression of 71 (69 upregulated and two downregulated) DEGs in whiteflies treated with A. attenuatus compared to whiteflies from the control treatment (BtA-Vs-Bt) (Fig. 3A), among which the majority of genes had Log2FC values within the range of 1 to 2 (69 genes; 97.18%). The number of genes with Log2FC values <-1 was 2 (2.81%) (Additional file 2A). A total of 1,194 (936 upregulated and 258 downregulated) DEGs were observed in whiteflies treated with matrine compared with those from the control treatment group (BtM-Vs-Bt) (Fig. 3B). Among these genes, the majority had Log2FC values in the range of 1 to 2 (790 genes; 66.16%). There were 146 genes (12.22%) with Log2FC values >2, and 258 genes (21.60%) with Log2FC values <-1 were observed (Additional file 2B). Our findings revealed differential expression of 51 (35 upregulated and 16 downregulated) DEGs in whiteflies treated with A. attenuatus + matrine compared to whiteflies from the control treatment (BtAM-Vs-Bt) (Fig. 3C), among which the majority of genes had Log2FC values in the range of 1 to 2 (35 genes; 68.62%), whereas 16 genes (31.37%) with Log2Ratio values <-1 were observed (Additional file 2C).

Fig. 3
figure 3

Differently expressed genes (DEGs) in B. tabaci treated with different treatments. A A. attenuatus versus control; B matrine versus control; C A. attenuatus + matrine versus control and D Venn diagram of the differentially expressed genes

The cluster analysis of DEGs showed that the Log2FC values of most genes in the treatment groups were positive, which indicated that these genes were upregulated compared with those in the control groups (Additional file 2D).

DEGs following treatment with A. attenuatus

Among the 69 genes that were upregulated in whiteflies treated with A. attenuatus compared to whiteflies from the control treatment, 19 were classified as orphan genes because they did not show any homology to known proteins. Another major category of upregulated DEGs in whiteflies treated with A. attenuatus consisted of five cuticle proteins (Bta07904, cuticle protein; Bta07415, cuticle protein; Bta00253, cuticular protein 73D; Bta07899, cuticular protein BJ; Bta09581, cuticle protein). A group of upregulated DEGs in whiteflies treated with A. attenuatus was fatty acyl CoA reductase I consisting of three genes (Bta07543, Bta01237, Bta14907) (Additional file 2A). The main upregulated DEGs involved in B. tabaci response to A. attenuatus were vitellogenin, cuticle protein, fatty acyl-CoA reductase 1, and catalase (Table 1).

Table 1 Top differentially expressed genes (DEGs) in whitefly treated with A. attenuates (A), matrine (B), and A. attenuates + matrine (C)

Only two downregulated DEGs were found in whiteflies treated with A. attenuatus compared to whiteflies from the control treatment. The two downregulated DEGs included von Willebrand factor C domain-containing protein 2-like (Bta01761) and CTLMA2 (Bta05757) (Table 1).

DEGs following treatment with matrine

Numerous upregulated genes (936) were observed in whiteflies treated with matrine compared with those from the control treatment. Among the 936 upregulated DEGs, five major categories were identified, namely 245 orphan genes, 21 cathepsins, 19 cuticle proteins, 17 cytochrome P450, and five fatty acyl-CoA reductase 1 genes. Two categories, orphan genes and cuticle proteins, were also differentially regulated in whiteflies treated with A. attenuatus compared to those from the control treatment. Only 18 out of 245 orphan genes that were differentially expressed between whiteflies treated with matrine were similar to those differentially expressed in whiteflies treated with A. attenuatus. Furthermore, four cuticle protein genes (Bta07904, Bta00253, Bta07899, and Bta09581) were similar to those differentially expressed in whiteflies treated with A. attenuatus (Additional file 2 B). The major upregulated DEGs involved in B. tabaci response to matrine were alpha-glucosidase I, cytochrome P450, vitellogenin, fatty acyl-CoA reductase 1, protein LEA isoform M, and cathepsin F-like protease, which were also upregulated with Log2FC values in the range of 2–3 (Table 1).

Among the 258 genes downregulated in whiteflies treated with matrine compared to whiteflies from the control treatment, 60 were classified as orphan genes. Another category of downregulated DEGs in whiteflies treated with matrine consisted of three cuticle proteins (Bta08284, cuticle protein 6; Bta14107, cuticle protein 1; and Bta01013, cuticular protein, putative). One group of downregulated DEGs in whiteflies treated with matrine was cytochrome P450, which consisted of three genes (Bta13899, Bta11266, and Bta06109). The two downregulated DEGs involved in the response of B. tabaci to matrine were von Willebrand factor C domain-containing protein 2-like (Bta01761) and CTLMA2 (Bta05757) (Table 2).

Table 2 Primers used in the study

DEGs following treatment with A. attenuatus + matrine

In total, 35 DEGs were upregulated in whiteflies treated with A. attenuatus + matrine compared to those from the control treatment. Among the 35 upregulated DEGs, three major categories were identified: 12 orphan genes, three cuticle proteins, and two histone H1 genes. As described above, cuticle proteins were also differentially regulated in whiteflies treated individually with A. attenuatus and matrine compared to whiteflies from the control treatment. Only one of three cuticle protein genes (Bta09545) was similar to those differentially expressed in whiteflies treated with matrine. Furthermore, two upregulated DEGs (Bta09545: Protein LEA-1and Bta05911: Cathepsin F-like protease) were also differentially regulated in whiteflies treated with individual treatments of A. attenuatus and matrine compared to whiteflies from the control treatment (Additional file 2C). The major upregulated DEGs involved in B. tabaci response to A. attenuatus + matrine were transmembrane protease serine 3, cuticle protein, tubulin alpha-3 chain, and histone H1.1 (Table 2).

Among the 16 genes that were downregulated in whiteflies treated with A. attenuatus + matrine compared to whiteflies from the control treatment, four were classified as orphan genes. The twodownregulated DEGs involved in the response of B. tabaci to matrine were von Willebrand factor C domain-containing protein 2-like (Bta01761) and CTLMA2 (Bta05757) (Table 2).

Common DEGs between whiteflies treated with individual treatments of A. attenuatus and matrine

Of the 72 and 1,194 DEGs found in whiteflies treated with individual treatment of A. attenuatus and matrine, respectively, there were 65 common genes that were found to be differentially expressed in both treatments. These common DEGs included several major categories: 1) 18 orphan genes; 2) four cuticle proteins (Bta07904, Bta00253, Bta09581, and Bta07899); 3) three cathepsins (Bta02553 and Bta01771, cathepsin B; Bta05911, cathepsin F-like proteases); 4) three fatty acyl-CoA reductase 1 (Bta01237, Bta07543,and Bta14907); and 5) two glucose dehydrogenase genes (Bta11629; Bta10911) (Additional file 2D).

Common DEGs between whiteflies treated with individual and combined treatments of A. attenuatus and matrine

Of the 72, 1,194, and 51 DEGs found in whiteflies treated with individual or combined treatments of A. attenuatus and matrine, respectively, two genes (one protein LEA-1 isoform gene (Bta00818) and 2) cathepsin F-like protease (Bta05911)) were differentially expressed between all three treatments (Additional file 2E).

Correlation between DEGs and mortality of whiteflies treated with individual and combined treatments of A. attenuatus and matrine

The correlation analysis between DEGs and mortality of whiteflies treated with individual and combined treatments of A. attenuatus and matrine revelaed a negative correlation (-0.768) between both variables.

KEGG pathway analysis

KEGG pathway enrichment analysis of DEGs was performed to identify the potential pathways that were upregulated/ downregulated in whiteflies treated with different treatments (A. attenuatus; matrine; A. attenuatus + matrine). In whiteflies treated with A. attenuatus, the most enriched pathways were “Longevity regulating pathway - worm,” “Renin secretion,” and “Antigen processing and presentation” of Organismal Systems; “Peroxisome,” “Lysosome,” and “Apoptosis” of Cellular Processes, and “Glycosphingolipid biosynthesis - ganglio series, Glycosphingolipid biosynthesis - globo and isoglobo series of Metabolism (Fig. 4A, Additional file 3). Furthermore, 80% (57 of 71) of the total number of upregulated DEGs in whiteflies treated with A. attenuatus were annotated to KEGG pathways. The three major categories of pathways most represented as upregulated in whiteflies treated with A. attenuatus were 1) metabolism, 2) organismal systems, and 3) cellular processes (Additional file 4A). KEGG analysis was not performed for downregulated DEGs in whiteflies treated with A. attenuatus because two DEGs were downregulated.

Fig. 4
figure 4

Most enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in B. tabaci treated with different treatments. A A. attenuatus versus control; B matrine versus control; and C A. attenuatus + matrine versus control

In whiteflies treated with matrine, 47 KEGG pathways were significantly enriched. Within the top 10 enriched pathways, most of them were associated with the Metabolism category, including the “Retinol metabolism,” “Antigen processing and presentation,” “Starch and sucrose metabolism,” “Pentose and glucuronate interconversions,” “Ascorbate and aldarate metabolism,” “Porphyrin and chlorophyll metabolism,” and “Amino sugar and nucleotide sugar metabolism” (Fig. 4B, Additional file 3). Furthermore, 90% (1,077 of 1,194) DEGs were annotated to KEGG pathways among the upregulated DEGs, and accounted for four main categories: 1) metabolism, 2) human diseases, 3) organismal systems, and 4) cellular processes (Additional file 4B). The downregulated DEGs in whiteflies treated with matrine were associated with the following main categories: 1) organismal system, 2) metabolism, 3) environmental information processing, and 4) human diseases (Additional file 4C).

In whiteflies treated with A. attenuatus + matrine, the most enriched pathways were “Systemic lupus erythematosus” of the Immune disease subcategory and “Apoptosis” of the Cell growth and death subcategory (Figure 4C, Additional file 3). In total, 62% (32 out of 51) of the total number of upregulated DEGs in whiteflies treated with A. attenuatus + matrine were annotated to KEGG pathways. The three major categories of pathways most represented as upregulated in whiteflies treated with A. attenuatus were metabolism, human diseases, and cellular processes (Additional file 4D).

Quantitative real-time polymerase chain reaction (qRT- PCR) validation of DEGs

To validate the digital gene expression map (DGE) data, the transcript levels of eight randomly selected genes (six upregulated and two downregulated) were quantified by qRT-PCR, using RNA extracted from three independent biological replicates. In all three biological replicates, ten genes showed concordant changes between the DGE and qRT-PCR data (Figure 5 and Additional file 5), suggesting that the DGE data were reliable.

Fig. 5
figure 5

qRT-PCR validation of differently expressed genes (DEGs) expressed by B. tabaci in response to different treatments. (A) A. attenuatus; (B) matrine; and (C) A. attenuatus + matrine. Error bars represent standard deviation of the mean of three biological replicates

Discussion

Akanthomyces attenuatus Zare & Gams is a well-known pathogen of insect pests, including whiteflies, aphids, and thrips. Recently, many studies have investigated the joint application of pathogenic insect fungi and natural enemies, synthetic oils, bacterial metabolite avermectins, and other biopesticides. Our recent studies on the combined application of A. attenuatus and matrine (a botanical insecticide) against B. tabaci have shown promising results. To date, no study has been conducted to elucidate the immune response of B. tabaci to the combined application of A. attenuatus and matrine at the molecular level. In the present study, we obtained a transcriptome dataset containing 35,895 unigenes of B. tabaci using RNA-Seq. The analysis of DEGs involved in the immune response to single or joint treatments of A. attenuatus and matrine revealed the expression of 5, 54, and 0 immunity-related genes for BtA, BtM, and BtAM, respectively. Most immunity-related genes are also involved in NOD-like receptor signaling, Toll and Imd signaling, Fc gamma R-mediated phagocytosis, and antigen processing and presentation pathways. Our findings differ from those of Xia et al. [21], who examined the expression of B. tabaci genes in response to infection with Beauveria bassiana. Their findings revealed that 654 and 1,681 genes were differentially expressed at 48 and 72 h post-infection, respectively, compared to the control. Chen et al. [26] also observed a lower number of DEGs in Megalurothrips usitatus post-A. attenuatus infection when compared to infection with Beauveria brongniartii. The lower number of DEGs and zero expression of immunity-related genes by B. tabaci treated with combined treatments of A. attenuatus and matrine can be related to the disruption of the defense/immune system of the host [7]. The DEGs were grouped according to their functional roles.

Differential expression of metabolism and biogenesis-related genes

Insects respond to external stress by increasing carbohydrate metabolism and energy production [27]. Major differences were observed within the metabolism category between whiteflies subjected to different treatments. In total, 19 metabolic pathway genes (all upregulated) were found to be differentially expressed in whiteflies treated with A. attenuatus compared with 359 (324 upregulated and 35 downregulated) metabolic pathway DEGs in whiteflies treated with matrine and 12 metabolic pathway genes (all upregulated) in whiteflies treated with A. attenuatus + matrine. We observed an overall lower differential expression of metabolism-related genes (after treatment with A. attenuatus or matrine + A. attenuatus), which may play an important role in the complex interactions between B. tabaci and A. attenuatus. B. tabaci may reduce general metabolism and conserve energy to fight fungal infections with alternative immune responses or reproductive strategies [28,29,30,31]; however, this global downregulation could be regulated by fungal effectors that suppress host genetic information processing and metabolism.

Membrane-bound acyl-lipid desaturases and membrane-bound acyl-CoA desaturases are two types of membrane-bound fatty acid desaturases. Insect acyl-CoA desaturases play a critical role in semiochemical biosynthesis and sensing, cold tolerance, defensive fatty acid production, lipid metabolism, feeding behavior, and larval development [32,33,34,35,36]. Our results showed the upregulation of the acyl-CoA desaturase gene (Bta11879) in response to individual or combined treatments with pathogens and botanicals. The higher expression levels of the acyl-CoA desaturase gene in B. tabaci following different treatments suggest the acceleration of lipid metabolism to produce the necessary energy required for the induction of a strong immune response [37].

Cytochrome P450 monooxygenases are enzymes found throughout the animal kingdom and are mainly involved in metabolism [38]. Toxins can trigger cytochrome P450 production in the diet, which makes this adaptation noteworthy [39]. In addition, cytochrome P450 activities are known to be involved in insecticide metabolism, resulting in molecular activation or, more broadly, detoxification [40, 41]. Cytochrome P450, which is responsible for the specificity of the reaction, is the primary source of proteins in each of these enzymatic systems. In this study, cytochrome P450 was upregulated in response to individual treatments with matrine and A. attenuatus + matrine, confirming the increased detoxification response of B. tabaci to the different treatments.

Differential expression of signaling-related genes

Among the genes annotated to signal transduction pathways, four upregulated DEGs were observed in whiteflies treated with A. attenuatus. In whiteflies treated with matrine, 123 DEGs (89 upregulated and 34 downregulated) were associated with signal transduction pathways, whereas only one upregulated DEG was annotated to signal transduction pathways.

In insects, FAR genes involved in sex pheromone production (pheromone-gland-specific fatty acyl reductase) form a distinct clade with significant changes in gene structure, fatty acyl substrate specificity, and selectivity [42]. Fatty acyl-CoA reductase (FAR) is the key enzyme responsible for the biosynthesis of fatty alcohols and can catalyze the transformation of fatty acids to fatty alcohols and NADPH [43, 44]. Our results showed upregulation of the fatty acyl-CoA reductase 1 gene (Bta11879) in response to individual treatments with pathogens and botanicals. These findings suggest that fatty acid metabolism may be involved in the detoxification of pathogens and botanicals in whiteflies. Upregulated expression of fatty acyl-CoA reductase 1 gene may produce fatty alcohol and wax-ester in the insect epicuticle [45].

Histone H1 is one of the important components of eukaryotic chromatin. This protein is known to play a key role in the chromatin folding into higher-order structure (Allan et al. 1980). It was shown to be a nonspecific repressor of gene expression [46]. Besides, histone HI appears to compete with specific gene activators for regulatory sites of tissue-specific genes, thus determining their ability for further transcription [46]. Moreover, HI is a member of the regulatory cascade involved in the action of nuclear receptors which leads to selective activation of chromatin. Several kinds of indirect evidence suggest that molecular variations in these proteins are likely to have pleiotropic effects, including effects on quantitative traits (growth, fecundity etc.) [47]. The upregulation of Histone H1 gene in whitefly nymphs treated with matrine+A. attenuatus might be related to their effort to cope with external stress through higher transcription or growth activities.

Differential Expression of immunity-related genes

The analysis of DEGs that may participate in the immune response to single or combined treatments of A. attenuatus and matrine (BtA, BtM, and BtAM) revealed some interesting findings. Among whiteflies treated with A. attenuatus, 5 immunity-related (five upregulated, zero downregulated) DEGs were observed. These DEGs were mainly involved in the NOD-like receptor signaling pathway (map04621), the Toll and Imd signaling pathway (map04624), and the antigen processing and presentation pathway (04612). In whiteflies treated with matrine, 54 immunity-related DEGs (51 upregulated and three downregulated) were observed. The majority of immunity DEGs were annotated to the NOD-like receptor signaling pathway (map04621), Toll and Imd signaling pathway (map04624), Fc gamma R-mediated phagocytosis (map 04666), and antigen processing and presentation pathway (04612). In whiteflies treated with A. attenuatus + matrine, no immunity-related gene was expressed, which confirmed the complete disruption of the B. tabaci immune response by the combined application of A. attenuatus and matrine Vitellogenins are insect proteins that provide nutrients for yolk formation [48, 49]. Many studies have revealed that vitellogenins play an essential role in the immune system of different insect species [48, 50, 51]. In this study, it was also found that the expression of some vitellogenins was upregulated after treatment with pathogens (A. attenuatus), botanicals (matrine), and their combination (A. attenuatus + matrine), thereby confirming the role of vitellogenins in B. tabaci defense against different treatments.

In cysteine proteases, cysteine residues modulate catalytic processes [27]. These endopeptidases can regulate digestion, embryonic vitellin degradation, and metamorphosis in insects [28, 29]. Cysteine cathepsins (B, F, L, and O) are the most important cysteine proteases observed in insects, and are well known for their endopeptidase activity [30]. Recent studies have shown that cathepsin F-like enzymes are, quantitatively, the most important cysteine proteases in insects [31]. Miyaji et al. [32] showed that CathF was involved in the maintenance of endogenous cysteine protease activity in the larval hemolymph of Manduca sexta. Our results showed the upregulation of cathepsin F-like protease gene (Bta05911) in response to individual or combined treatments of pathogen and botanical. The upregulated expression of cathepsin F-like protease genes in response to different treatments shows the importance of CathF in the immune defense of B. tabaci through antigen processing and presentation [32].

The insect cuticle is the first barrier against pathogen infection and maintains the form and movement of the insects [52]. The cuticle is a protective envelope that can completely enclose and protect insect larvae from the external environment [32]. Cuticles mainly comprise chitin, lipids, catecholamines, minerals, and proteins [33]. Cuticle proteins are induced to improve the stability of the cuticle required for insect survival under external stress [52]. One of the routes of insecticide resistance is cuticle thickening [34]. Environmental stress enhances the expression of cuticle proteins in Leptinotarsa decemlineata in response to stress, allowing the insects to adapt to new or changing settings [35]. Silva et al. [36] found that cuticle protein genes of aphids were upregulated following insecticide application. In this study, the expression of cuticle protein was upregulated in response to individual treatments of pathogen (A. attenuatus) and botanical (matrine), whereas in the case of A. attenuatus + matrine, there was no or minimal expression of cuticle protein genes, which indicated disruption of the B. tabaci immune process.

Late embryogenesis abundant (LEA) proteins are involved in insect response to environmental stress and thus act as water deficiency stress protectants [53]. The LEA proteins act as a chemical chaperone for trehalose (the most abundant hemolymph sugar in insects) and may preserve cell components from desiccation by replacing water or forming glass [54]. Kikawada et al. [38] reported that LEA protein 1 was accumulated in response to desiccation or NaCl treatment in Polypedilum vanderplanki larvae. Our results showed the upregulation of the LEA protein-1 gene (Bta000818) in response to individual or combined treatments of pathogen and botanical. After different treatments, the upregulation of the LEA protein-1 gene by B. tabaci confirmed a strong immunological/detoxification response through trehalose accumulation to survive drying in response to infection [39].

C-type lectins (CTLs) are involved in the binding site of glycans in direct coordination with lectin-bound Ca2+ [55]. CTLs are divided into two classes based on their binding preference. CTLs with the EPN sequence at the calcium-binding site prefer mannose-type sugars (fucose, glucose), whereas CTLs with the QPD sequence prefer galactose-type sugars [56]. The CTL proteins CTL4 and CTLMA2 have been reported to influence the insect immune response during microbial infection by stabilizing phenol oxidase activity and reducing the downstream melanization response [57, 58]. Our results showed that the CLTMA2 gene (Bta05757) is downregulated in response to individual or combined treatments with pathogens and botanicals. Lower expression levels of CLTMA2 in B. tabaci following different treatments could have resulted in a weaker immune response and increased melanization, implying susceptible whitefly individuals [59].

Single von Willebrand factor C-domain proteins (SVWCs) are single-domain von Willebrand factor type C proteins known primarily from arthropods [60, 61]. Kleinjung et al. [62] discovered the first SVWC in Drosophila melanogaster. SVWC proteins are involved in response to environmental stressors, changes in nutritional status, and microbial infections [63]. Our results showed downregulation of the von Willebrand factor type C protein-2 gene (Bta01761) in response to individual or combined treatments of pathogens and botanicals. The lower expression levels of the von Willebrand factor type C protein-2 gene in B. tabaci following different treatments resulted in a reduced immune response during pathogen invasion through reduced phagocytosis and antimicrobial responses [55, 56].

A large amount of energy is normally required for a rapid and effective immune response [64]. In insects, energy is derived from carbohydrate sources such as glycogen and trehalose. These two stored sugars are enzymatically converted into glucose, which is then used in the production of high-energy ATP through glycolysis and the tricarboxylic acid (TCA) cycle [65]. This metabolic switch is clearly essential for antimicrobial responses, because starvation or blocking the signal transduction for metabolic reprogramming has been found to decrease the immune response against pathogens [66]. Our results further support the possible significance metabolic switching in antimicrobial responses as the lower differential expression of metabolic pathways related genes in whiteflies treated with A. attenuatus and matrine+A. attenuatus treatments ultimately led to the lower expression of immunity related genes. In short, reduced/disrupted metabolic response might have led to the lesser number of diferentially expressed genes in whiteflies exposed to matrine+A. attenuatus when compared to the remain treatments.

Conclusion

In summary, we investigated the complex interactions between the whitefly and individual or combined treatments of A. attenuatus and matrine using RNA sequencing. Our data indicated that genes such as cytochrome P450, cuticle protein, vitellogenin, cysteine proteases, late embryogenesis abundant (LEA) proteins, desaturases, CTLMA2, and von Willebrand factor type C protein-2 were mainly involved in whitefly defense responses against A. attenuatus and matrine. The results also revealed that AMPK signaling, apoptosis, and drug metabolism pathways are likely involved in whitefly defense responses against A. attenuatus and matrine infection. To the best of our knowledge, this is the first report on the molecular interactions occurring between whiteflies and individual or combined treatments of A. attenuatus and matrine. These results will further improve our knowledge of the infection mechanism and complex biochemical processes involved in the synergistic action of A. attenuatus and matrine against B. tabaci.

Methods

Insect cultures

Whitefly cultures were reared on cotton plants according to the method described by Wang et al. [3] at the Engineering Research Centre of Biological Control, Ministry of Education, South China Agricultural University, under standard laboratory conditions (26 °C, 70% relative humidity (RH), and 12:12 h light/dark photoperiod).

Preparation of fungal inoculum and whitefly treatment

Akanthomyces attenuatus isolate SCAUDCL53 (deposited at the Key Laboratory of Biopesticides Innovation and Application of Guangdong Province, South China Agricultural University, Guangzhou, P.R. China) was cultured according to the method described by Ali et al. [67] to calibrate the basal fungal suspension of 1×106 conidia/ml.

Matrine (98% purity) was purchased from Guangdong New Scene Bioengineering Company Ltd., Yangjiang, P.R. China, and the basal concentration of matrine (1.0 mg/mL) was prepared following the methodology of Wu et al. [17].

Cotton leaves having 2nd instar B. tabaci nymphs were individually dipped in A. attenuatus fungal suspension (1×06 conidia/ml), matrine (1.0 mg/L), and A. attenuatus and matrine (1×06 conidia/ml + 1.0 mg/L) for 10 seconds. The leaves were air dried before being placed in a 9 cm agar medium with 20% agar to maintain moisture. Bemisia tabaci dipped in 0.05% Tween 80 served as control. Each treatment consisted of four cotton leaves with 100 B. tabaci nymphs per leaf. The experimental setup was maintained at 25 °C, 85% RH, and a 16:8 h (light/dark photoperiod). The entire experiment was repeated three times (with a fresh batch of insects and new conidial suspension). The mortality of the insects was tracked daily for 5 d after treatment.

RNA extraction for sequencing

The surviving B. tabaci nymphs (200 nymphs pooled per sample and three samples for each treatment considered as three replicates) from different treatments were collected in 1.5 mL Eppendorf tubes after approximately 3 d of treatment, frozen in liquid nitrogen, and stored at -80 °C. Total RNA was extracted using TRIzol reagent (Thermo Fisher Scientific, Germany), according to the manufacturer’s instructions (Invitrogen). The quality of the extracted RNA was determined using 1% agar-gel electrophoresis and a NanoDrop One ultraviolet spectrophotometer (Thermo Fisher Scientific, Germany), and the RNA integrity number (IN) value was calculated using an Agilent 2100 bioanalyzer (Agilenet, USA).

Construction and sequencing of cDNA library

The mRNAs (post RNA extraction) were enriched with Oligo (dT) having magnetic beads treated with DNase I. The mixture was then supplemented with fragmentation buffer to prepare the mRNA templates for the sysnthesis of first cDNA chain whereas RNA enzyme, DNA polymerase I, and dNTPs were added to synthesize the second cDNA chain. AMPure XP beads were used for cDNA (double-stranded) purification followed by the selection of cDNA library through PCR enrichment.. The library was sequenced with Illumina sequencing platform.

Assembly and annotation of the transcriptome

Raw reads were filtered and assembled using Cufflinks (http://cole-trapnell-lab.github.io/cufflinks/). The assembled sequences were removed from redundancy and spliced using TGICL to obtain the longest nonredundant unigene set. Further statistical analysis and quality control were performed using the unigene set. The unigene set obtained by de Novo and clustering was compared with the database using BLAST for functional annotation. The NR, Swiss-Prot, Pfam, EggNOG, GO, and KEGG databases were used for annotation. The annotation information of genes or transcripts could be comprehensively obtained, and the annotation of each database was determined. To understand the transcriptome database of the response of B. tabaci to A. attenuatus and matrine alone or in combination, 12 DGE marker libraries were constructed, including the transcription information of B. tabaci treated with A. attenuatus and matrine alone or in combination with a blank control group. The treatment groups were named as BtA-1, BtA-2, BtA-3, BtM-1, BtM-2, BtM-3, BtAM-1, BtAM-2, BtAM-3, Bt-1, Bt-2, and Bt-3. Table 2 summarizes the basic statistics of the 12 libraries.

Digital gene expression profile analysis

Using the constructed cDNA library data, a DGE was used to identify the DEGs. We regarded genes with a false discovery rate ≤ 0.001 and |fold change| ≥ 2 as DEGs. Genes with similar expression patterns often exhibit similar functions. The R toolkit pheatMap software was used for cluster analysis of DEGs and experimental conditions. GO-seq was used for the GO enrichment analysis. The main difference between this method and the ordinary hypergeometric distribution is that the influence of gene length preference can be eliminated so that the GO term of true enrichment can be calculated more accurately. Combined pathways (KEGG) with the above pathways [68] and GO functional categories with significant enrichment of DEGs were obtained for further analysis.

qRT-PCR verification of DGEs

Primer Premier 5 software (Premier Biosoft, San Francisco, USA) was used to create primers for the DEGs and a reference gene (Table 2). To validate the DEGs, 18 DEGs were selected for qRT-PCR analysis of the RNA samples used for sequencing. The expression analysis was performed using three biological replicates. SYBR Premix Ex Taq (Perfect Real Time, Takara Bio Inc., Japan) was used for the final reaction. The qRT-PCR system had a total volume of 20 μL, which included 0.4 μL of forward/reverse primers, 1 μL of 10-fold cDNA, 8.2 μL of nuclease-free Water, and 10 μL of 2×Perfect Start TM Green qPCR Super Mix. The qRT-PCR was conducted on a Bio-RAD iQ5 thermocycler with the following temperature settings: initial denaturation at 94 °C for 30 minutes, followed by 39 cycles of 94 °C for 5 seconds, 55 °C for 30 minutes, and 72 °C for 10 seconds, followed by a melt curve of 65-95 °C, incrementing 0.5 °C for 0.05 seconds. As described earlier, the 2-ΔΔCt technique was used to calculate relative expression levels using β-actin as an internal control [21, 25, 26, 69, 70].

Availability of data and materials

The sequencing data generated in this study were deposited in the NCBI SRA database. The raw datasets generated in this study are available at http://www.ncbi.nlm.nih.gov/bioproject/PRJNA797870 (NCBI accession number: PRJNA797870).

References

  1. Sun TF, Wu JH, Ali S. Morphological and molecular identification of four Purpureocillium isolates and evaluating their efficacy against the sweet potato whitefly, Bemisia tabaci (Genn.) (Hemiptera: Aleyrodidae). Egypt J Biol Pest Cont. 2021;31:27.

    Article  Google Scholar 

  2. Can Z, Ali S, Musa PD, Wang XM, Qiu BL. Evaluation of the pathogenecity of Aschersonia aleyrodis on Bemisia tabaci in the laboratory and greenhouse. Biocont Sci Technol. 2017;27:210–21.

    Article  Google Scholar 

  3. Wang XS, Xu J, Wang XM, Qiu BL, Cuthbertson AGS, Du CL, et al. Isaria fumosorosea-based-zero-valent iron nanoparticles affect the growth and survival of sweet potato whitefly, Bemisia tabaci (Gennadius). Pest Manag Sci. 2019;75:2174–81.

    Article  CAS  Google Scholar 

  4. Stansly PA, Naranjo SE. Bemisia: Bionomics and Management of a Global Pest. Dordrecht, Heidelberg, London and New York: Springer; 2010. p. 558.

    Book  Google Scholar 

  5. He S, Wu XY, Zheng JZ, Han ZQ, Zhao CX, Wang YY, et al. A study of control effect of Bemisia tabaci (Gennadius) on greenhouse tomatoes using parasitoid Encarsia formosa Gahan. J Anhui Agri Sci. 2013;41:6244–8 (in Chinese).

    Google Scholar 

  6. Revathi K, Chandrasekaran R, Tanigaivel A, Kirubakaran SA, Nathan SS. Biocontrol efficacy of protoplast fusants between Bacillus thuringiensis and Bacillus subtilis against Spodoptera litura Fabr. Arch Phytopathol Plant Prot. 2014;47:1365–75.

    Article  CAS  Google Scholar 

  7. Ali S, Zhang C, Wang ZQ, Wang XM, Wu JH, Cuthbertson AGS, et al. Toxicological and biochemical basis of synergism between the entomopathogenic fungus Lecanicillium muscarium and the insecticide matrine against Bemisia tabaci (Gennadius). Sci Rep. 2017;7:46558.

    Article  CAS  Google Scholar 

  8. Kovac M, Gorczak M, Wrzosek M, Tkaczuk C, Pernek M. Identification of entomopathogenic fungi as naturally occurring enemies of the invasive oak lace bug, Corythucha arcuata (Say) (Hemiptera: Tingidae). Insects. 2020;11:679.

    Article  Google Scholar 

  9. Yang B, Du CL, Ali S, Wu JH. Molecular characterization and virulence of fungal isolates against the bean flower thrips, Megalurothrips usitatus Bagnall (Thysanopera: Thripidea). Egypt J Biol Pest Cont. 2020;30:50.

    Article  Google Scholar 

  10. Du C, Yang B, Wu JH, Ali S. Identification and virulence characterization of two Akanthomyces attenuatus isolates against Megalurothrips usitatus (Thysanopera: Thripidea). Insects. 2019;10:168.

    Article  Google Scholar 

  11. Mohammad AA, Khadim JH, Hassan AMH. Laboratory evaluation of entomopathogenic fungi for the control of khapra beetle (Coleoptera: Dermestidae) and their effects on the beetles' fecundity and longevity. J Agri Urban Entomol. 2019;1:1–11.

    Google Scholar 

  12. Bitencourt RDOB, dos Santos Mallet JR, Mesquita E, Golo PS, Fiorotti J, Bittencourt VREP, et al. Larvicidal activity, route of interaction and ultrastructural changes in Aedes aegypti exposed to entomopathogenic fungi. Acta Tropica. 2021;213:105732.

    Article  Google Scholar 

  13. Koodalingam A, Dayanidhi MK. Studies on biochemical and synergistic effects of immunosuppressive concentration of imidacloprid with Beauveria bassiana and Metarhizium anisopliae for enhancement of virulence against vector mosquito Culex quinquefasciatus. Pesti. Biochem Physiol. 2021;176:104882.

    Article  CAS  Google Scholar 

  14. Bukhari T, Takken W, Koenraadt CJ. Development of Metarhizium anisopliae and Beauveria bassiana formulations for control of malaria mosquito larvae. Parasite Vectors. 2011;4:23.

    Article  Google Scholar 

  15. Deng SQ, Zou WH, Li DL, Chen JT, Huang Q, Zhou LJ, et al. Expression of Bacillus thuringiensis toxin Cyt2Ba in the entomopathogenic fungus Beauveria bassiana increases its virulence. PLoS Negl Trop Dis. 2019;13:e0007590.

    Article  CAS  Google Scholar 

  16. Noskov YA, Polenogova OV, Yaroslavtseva ON, Belevich OE, Yurchenko YA, Chertkova EA, et al. Combined effect of the entomopathogenic fungus Metarhizium robertsii and avermectins on the survival and immune response of Aedes aegypti larvae. Peer J. 2019;7:e7931.

    Article  Google Scholar 

  17. Wu JH, Yu X, Wang XS, Tang LD, Ali S. Matrine enhances the pathogenicity of Beauveria brongniartii against Spodoptera litura (Lepidoptera: Noctuidae). Front Microbiol. 2019;10:1812.

    Article  Google Scholar 

  18. Wu JH, Yang B, Zhang XC, Cuthbertson AGS, Ali S. Synergistic interaction between the entomopathogenic fungus Akanthomyces attenuatus (Zare & Gams) and the botanical insecticide matrine against Megalurothrips usitatus (Bagrall). J Fungi. 2021;7:536.

    Article  CAS  Google Scholar 

  19. Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS. De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010;11:400.

    Article  Google Scholar 

  20. Wang XW, Luan JB, Li JM, Su YL, Xia J, Liu SS. Transcriptome analysis and comparison reveal divergence between two invasive whitefly cryptic species. BMC Genomics. 2011;12:458.

    Article  CAS  Google Scholar 

  21. Xia J, Zhang CR, Zhang S, Li FF, Feng MG, Wang XW, et al. Analysis of whitefly transcriptional responses to Beauveria bassiana infection reveals new insights into insect-fungus interactions. PLoS ONE. 2013;8(7):e68185.

    Article  CAS  Google Scholar 

  22. Yang N, Xie W, Yang X, Wang S, Wu Q, Li R, et al. Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam. PLoS ONE. 2013;8(5):e61820.

    Article  CAS  Google Scholar 

  23. Upadhyay SK, Sharma S, Singh H, Dixit S, Kumar J, Verma PC. Whitefly genome expression reveals host-symbiont interaction in amino acid biosynthesis. PLoS ONE. 2015;10(5):e0126751.

    Article  Google Scholar 

  24. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18:1509–17.

    Article  CAS  Google Scholar 

  25. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.

    Article  CAS  Google Scholar 

  26. Chen Y, Yang B, Li Z, Yue Y, Tian Q, Chen W, Ali S. Wu J. Immune-Related Genes of Megalurothrips usitatus (Bagrall) against Beauveria brongniartii and Akanthomyces attenuatus Identified using RNA Sequencing. Front. Physiol. 2021;12: 671599.

  27. Merritt TJS, Bewick AJ. Genetic diversity in insect metal tolerance. Front Genetics. 2017;8:172.

    Article  Google Scholar 

  28. Alleyne M, Chappell MA, Gelman DB, Beckage NE. Effects of parasitism by the braconid wasp Cotesia congregata on metabolic rate in host larvae of the tobacco hornworm, Manduca sexta. J Insect Physiol. 1997;43:143–54.

    Article  CAS  Google Scholar 

  29. Altincicek B, Gross J, Vilcinskas A. Wounding-mediated gene expression and accelerated viviparous reproduction of the pea aphid Acyrthosiphon pisum. Insect Mol Biol. 2008;17:711–6.

    Article  CAS  Google Scholar 

  30. Gwynn D, Callaghan A, Gorham J, Walters K, Fellowes M. Resistance is costly: trade-offs between immunity, fecundity and survival in the pea aphid. Proc R Soc B. 1803;2005:272.

    Google Scholar 

  31. Li JM, Ruan YM, Li FF, Liu SS, Wang XW. Gene expression profiling of the whitefly (Bemisia tabaci) Middle East-Asia Minor 1 feeding on healthy and Tomato yellow leaf curl China virus-infected tobacco. Insect Sci. 2011;18:11–22.

    Article  Google Scholar 

  32. Miyaji T, Murayama S, Kouzuma Y, Kimura N, Kanost MR, Kramer KJ, et al. Molecular cloning of a multidomain cysteine protease and protease inhibitor precursor gene from the tobacco hornworm (Manduca sexta) and functional expression of the cathepsin F-like cysteine protease domain. Insect Biochem. Mol. Biol. 2010;40:835e846.

    Article  Google Scholar 

  33. Yu S, Ding L, Luo R, Li XJ, Yang J, Liu H. Identification of immunity-related genes in Dialeurodes citri against entomopathogenic fungus Lecanicillium attenuatum by RNA-Seq analysis. PLoS ONE. 2016;11:e0162659.

    Article  Google Scholar 

  34. Delon I, Payre F. Evolution of larval morphology in flies: get in shape with shaven baby. Trends Genet. 2004;20:305–13.

    Article  CAS  Google Scholar 

  35. Shang F, Ding BY, Ye C, Yang L, Chen TY, Xie JQ. Evaluation of a cuticle protein gene as a potential RNAi target in aphids. Pest Manag. Sci. 2020;76:134–40.

    Article  CAS  Google Scholar 

  36. Silva AX, Bacigalupe LD, Luna-Rudloffff M, Figueroa CC. Insecticide resistance mechanisms in the green peach aphid Myzus persicae (Hemiptera: Aphididae) II: costs and benefits. PLoS ONE. 2012;7:e36810.

    Article  CAS  Google Scholar 

  37. Barrett AJ, Rawlings ND, Wasner JF. Handbook of Proteolytic Enzymes. Academic press, London, 1998; 3932pp.

  38. Heidel-Fischer HM, Vogel H. Molecular mechanisms of insect adaptation to plant secondary compounds. Curr. Opin. Insect Sci. 2015;8:8–14.

    Article  Google Scholar 

  39. Hung CF, Prapaipong H, Berenbaum MR, Schuler MA. Differential induction of cytochrome P450 transcripts in Papilio polyxenes by linear and angular furanocoumarins. Insect Biochem. Molec. Biol. 1995;25:89–99.

    Article  CAS  Google Scholar 

  40. Li X, Baudry J, Berenbaum MR, Schuler MA. Structural and functional divergence of insect CYP6B proteins: from specialist to generalist cytochrome P450. Proc. Natl. Acad. Sci. U.S.A. 2004;101:2939–44.

    Article  CAS  Google Scholar 

  41. Taylor M, Feyereisen R. Molecular biology and evolution of resistance to toxicants. Mol. Biol. Evol. 1996;13:719–34.

    Article  CAS  Google Scholar 

  42. Antony B, Ding BJ, Moto K, Aldosari SA, Aldawood AS. Two fatty acyl reductases involved in moth pheromone biosynthesis. Sci. Rep. 2016;6:29927.

    Article  CAS  Google Scholar 

  43. Yang P, Zhu JY, Gong ZJ, Xu DL, Chen XM, Liu WW, et al. Transcriptome analysis of the Chinese white wax scale Ericerus pela with focus on genes involved in wax biosynthesis. PLoS One. 2012;7:e35719.

    Article  CAS  Google Scholar 

  44. Nguyen SH, Webb HK, Mahon PJ, Crawford RJ, Ivanova EP. Natural insect and plant micro-/nanostructsured surfaces: an excellent selection of valuable templates with superhydrophobic and selfcleaning properties. Molecules. 2014;19:13614–2630.

    Article  Google Scholar 

  45. Li XL, Zheng TX, Zheng XW, Han N, Chen XX, Zhang D. Molecular characterization of two fatty Acyl-CoA reductase genes from Phenacoccus solenopsis (Hemiptera: Pseudococcidae). J. Insect Sci. 2016;16:49.

    Article  Google Scholar 

  46. Zlatanova J. Histone HI and the regulation of transcription of eukaryotic genes. Trends in Biochem. Sci. 1990;15:273–6.

    Article  CAS  Google Scholar 

  47. Berdnikov VA, Rozov SM, Temnykh SV, Gorel FL, Kosterin OE. Adaptive nature of interspecies variation of histone H1 in insects. J. Mol. Evol. 1993;497-507.

  48. Kikawada T, Nakahara Y, Kanamori Y, Iwata K, Watanabe M, McGee B, et al. Dehydration-induced expression of LEA proteins in an anhydrobiotic chironomid. Biochem. Biophys. Res. Commun. 2006;348:56e61.

    Article  Google Scholar 

  49. Sakurai M, Furuki T, Akao K, Tanaka D, Nakahara Y, Kikawada T, et al. Vitrification is essential for anhydrobiosis in an African chironomid, Polypedilum vanderplanki. Proc. Natl. Acad. Sci. USA. 2008;105:5093e5098.

    Article  Google Scholar 

  50. Antúnez K, Mendoza Y, Santos E, Invernizzi C. Differential expression of vitellogenin in honey bees (Apis mellifera) with different degrees of Nosema ceranae infection. J. Apic. Res. 2013;52:227–34.

    Article  Google Scholar 

  51. Kodrík D, Ibrahim E, Gautam UK, Capková Frydrychová R, Bednárová A, Krištufek V. Changes in vitellogenin expression caused by nematodal and fungal infections in insects. J. Exp. Biol. 2019;222:jeb202853.

    Article  Google Scholar 

  52. Terra WR, Ferreira C. Biochemistry of digestion. In: Gilbert LI, Iatrov K, Gill SS, editors. Comprehensive Molecular Insect Science, vol. 4. Oxford: Biochem. Mol. Biology. Elsevier; 2005. p. 171–224.

    Chapter  Google Scholar 

  53. Zhang J, Goyer C, Pelletier Y. Environmental stresses induce the expression of putative glycine-rich insect cuticular protein genes in adult Leptinotarsa decemlineata (Say). Insect Mol. Biol. 2008;17:209–16.

    Article  Google Scholar 

  54. Silva AX, Jander G, Samaniego H, Ramsey JS, Figueroa CC. Insecticide resistance mechanisms in the green peach aphid Myzus persicae (Hemiptera: Aphididae) I: a transcriptomic survey. PLoS ONE. 2021;7:e36366.

    Article  Google Scholar 

  55. Bishnoi R, Sousa GL, Contet A, Day CJ, Hou CFD, Proftt LA, Singh D, Jennings MP. Valentine AM, Povelones M, Baxter RHG. Solution structure, glycan specifcity and of phenol oxidase inhibitory activity of Anopheles C-type lectins CTL4 and CTLMA2. Sci. Rep. 2019; 9: 15191.

  56. Weis WI, Taylor ME, Drickamer K. Te C-type lectin superfamily in the immune system. Immunol. Rev. 1998;163:19–34.

    Article  CAS  Google Scholar 

  57. Osta MA, Christophides GK, Kafatos FC. Effects of mosquito genes on Plasmodium development. Science. 2004;303:2030–2.

    Article  CAS  Google Scholar 

  58. Volz J, Müller H-M, Zdanowicz A, Kafatos FC, Osta MAA. Genetic module regulates the melanization response of Anopheles to Plasmodium. Cell Microbiol. 2006;8:1392–405.

    Article  CAS  Google Scholar 

  59. Shiao S-H, Whitten MMA, Zachary D, Hofmann JA, Levashina EA. Fz2 and cdc42 mediate melanization and actin polymerization but are dispensable for Plasmodium killing in the mosquito midgut. PLoS Pathog. 2006;2:e13.

    Article  Google Scholar 

  60. Qin N, Sun HH, Lu MK, Wang JH, Tang T, Liu FS. A single von Willebrand factor C-domain protein acts as an extracellular pattern-recognition receptor in the river prawn Macrobrachium nipponense. J. Biol. Chem. 2020;295(30):10468–77.

    Article  Google Scholar 

  61. Sheldon TJ, Miguel-Aliaga I, Gould AP, Taylor WR, Conklin D. A novel family of single VWC-domain proteins in invertebrates. FEBS Lett. 2007;581:5268–74.

    Article  CAS  Google Scholar 

  62. Kleinjung J, Romein J, Lin K, Heringa J. Contact-based sequence alignment. Nucleic Acids Res. 2004;32:2464–73.

    Article  CAS  Google Scholar 

  63. Han F, Lu A, Yuan Y, Huang W, Beerntsen BT, Huang J, et al. Characterization of an entomopathogenic fungi target integument protein, Bombyx mori single domain von Willebrand factor type C, in the silkworm. Bombyx mori. Insect Mol. Biol. 2017;26:308–16.

    Article  CAS  Google Scholar 

  64. Dolezal T, Krejcova G, Bajgar A, Nedbalova P, Strasser P. Molecular regulations of metabolism during immune response in insects. Insect Biochem. Mol. Biol. 2019;109:31–42.

    Article  CAS  Google Scholar 

  65. Tsai CH, Chuang YC, Lu YH, Lin CY, Tang CK, Wei SC, et al. Carbohydrate metabolism is a determinant for the host specificity of baculovirus infections. iScience. 2022;25:103648.

    Article  CAS  Google Scholar 

  66. Bajgar A, Kucerova K, Jonatova L, Tomcala A, Schneedorferova I, Okrouhlik J, et al. Extracellular adenosine mediates a systemic metabolic switch during immune response. Plos Biol. 2015;13:e1002135.

    Article  Google Scholar 

  67. Ali S, Huang Z, Ren SX. Media composition influences on growth, enzyme activity, and virulence of the entomopathogen hyphomycete Isaria fumosoroseus. Entomol. Exp. Appl. 2009;131:30–8.

    Article  CAS  Google Scholar 

  68. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  Google Scholar 

  69. Upadhyay SK, Dixit S, Sharma S, Singh H, Kumar J, Verma PC. siRNA Machinery in Whitefly (Bemisia tabaci). PLoS ONE. 2013;8:e83692.

    Article  Google Scholar 

  70. Thakur N, Upadhyay SK, Verma PC, Chandrashekar K, Tuli R, Singh PK. Enhanced whitefly resistance in transgenic tobacco plants expressing double stranded RNA of v-ATPase a gene. PloS ONE. 2014;9:e87235.

    Article  Google Scholar 

Download references

Acknowledgments

The authors would like to thank the editor and anonymous reviewers for their constructive comments and suggestions.

Funding

This research was funded by grants from the Key Science and Technology Project of CNTC (110202101053: LS-13), Science and Technology Plan Project of YNTC (2021530000242001), and the Science and Technology Program of Guangzhou, P.R. China (201804020070).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, J.H.W. and S.A.; methodology, X.M.W. and S.A; software, S.T.F.; validation, S.T.F.; investigation, S. T. F.; resources, S.A.; data curation, S.T.F.; writing—original draft preparation, S.A; writing—review and editing, S.A. and M.H.B; supervision, S.A.; project administration, S.A.; funding acquisition, X.M.W. and S.A. The authors read and approved the final manuscript.

Author information

Not applicable.

Corresponding authors

Correspondence to Xingmin Wang or Shaukat Ali.

Ethics declarations

Ethics approval and consent to participate

No specific permits were required for the samples used in this study. We complied with the relevant institutional, national, and international guidelines and legislation for this study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

a. The number of Illumina raw and processed reads produced per RNA-Seq library from whiteflies (Bemisia tabaci) treated with individual or combined treatments of Akanthomyces attenuatus and matrine. b. Length distribution of B. tabaci unigenes. c. Correlation matrix analysis for multiple biological replicates obtained from RNA-Seq libraries prepared from whiteflies (Bemisia tabaci) treated with individual or combined treatments of Akanthomyces attenuatus and matrine.

Additional file 2.

a. List of differentially expressed genes in whiteflies treated with Akanthomyces attenuatus. b. List of differentially expressed genes in whiteflies treated with matrine. c. List of differentially expressed genes in whiteflies treated with Akanthomyces attenuatus + matrine. d. Cluster analysis of differentially expressed genes in whiteflies treated with individual or combined treatments of Akanthomyces attenuatus and matrine. e. List of 66 common differentially expressed genes between whiteflies treated with individual treatments of Akanthomyces attenuatus and matrine (BtA and BtM). f. List of 66 common differentially expressed genes between whiteflies treated with individual and combined treatments of Akanthomyces attenuatus and matrine (BtA, BtM, and BtAM).

Additional file 3.

The main Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway in B. tabaci treated with Akanthomyces attenuatus matrine, Akanthomyces attenuatus + matrine, and controls.

Additional file 4.

Summary of Kyoto Encyclopedia of Genes and Genomes (KEGG) reference pathways of whiteflies associated with different treatments. (A) Akanthomyces attenuatus upregulated pathways; (B) matrine upregulated pathways; (C) matrine downregulated pathways; and (D) Akanthomyces attenuatus + matrine upregulated pathways.

Additional file 5.

 Results of quantitative real time PCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, J., Sun, T., Bashir, M.H. et al. Comparative transcriptome analysis reveals differences in gene expression in whitefly following individual or combined applications of Akanthomyces attenuatus (Zare & Gams) and matrine. BMC Genomics 23, 808 (2022). https://doi.org/10.1186/s12864-022-09048-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-09048-9

Keywords