Background The fungal pathogen, Beauveria bassiana, is an efficient biocontrol agent against a variety of agricultural pests. A thorough understanding of the basic principles of insect-fungus interactions may enable the genetic modification of Beauveria bassiana to enhance its virulence. However, the molecular mechanism of insect response to Beauveria bassiana infection is poorly understood, let alone the identification of fungal virulent factors involved in pathogenesis. Methodology/Principal Findings Here, next generation sequencing technology was applied to examine the expression of whitefly (Bemisia tabaci) genes in response to the infection of Beauveria bassiana. Results showed that, compared to control, 654 and 1,681genes were differentially expressed at 48 hours and 72 hours post-infected whiteflies, respectively. Functional and enrichment analyses indicated that the DNA damage stimulus response and drug metabolism were important anti-fungi strategies of the whitefly. Mitogen-activated protein kinase (MAPK) pathway was also likely involved in the whitefly defense responses. Furthermore, the notable suppression of general metabolism and ion transport genes observed in 72 hours post-infected B. tabaci might be manipulated by fungal secreted effectors. By mapping the sequencing tags to B. bassiana genome, we also identified a number of differentially expressed fungal genes between the early and late infection stages. These genes are generally associated with fungal cell wall synthesis and energy metabolism. The expression of fungal cell wall protein genes might play an important role in fungal pathogenesis and the dramatically up-regulated enzymes of carbon metabolism indicate the increasing usage of energy during the fungal infection. Conclusions/Significance To our knowledge, this is the first report on the molecular mechanism of fungus-whitefly interactions. Our results provide a road map for future investigations on insect-pathogen interactions and genetically modifying the fungus to enhance its efficiency in whitefly control.
References
[1]
Brown JK (2007) The Bemisia tabaci complex: genetic and phenotypic variability drives begomovirus spread and virus diversification. Plant Dis 1: 25–56.
[2]
De Barro PJ, Liu SS, Boykin LM, Dinsdale AB (2011) Bemisia tabaci: a statement of species status. Ann Rev Entomol 56: 1–19.
[3]
Sun D, Xu J, Luan J, Liu SS (2011) Reproductive incompatibility between the B and Q biotypes of the whitefly Bemisia tabaci in China: genetic and behavioural evidence. Bull Entomol Res 1: 1–10.
[4]
Wang P, Ruan YM, Liu SS (2010) Crossing experiments and behavioral observations reveal reproductive incompatibility among three putative species of the whitefly Bemisia tabaci. Insect Sci 17: 508–516.
[5]
Brown JK, Czosnek H (2002) Whitefly transmission of plant viruses. Adv Bot Res 36: 65–100.
[6]
Dalton R (2006) Whitefly infestations: the Christmas invasion. Nature 443: 898–900.
[7]
Jiu M, Zhou XP, Tong L, Xu J, Yang X, et al. (2007) Vector-virus mutualism accelerates population increase of an invasive whitefly. PLoS ONE 2: e182.
[8]
Seal SE, Bosch FVD, Jeger MJ (2006) Factors influencing begomovirus evolution and their increasing global significance: implications for sustainable control. Crit Rev Plant Sci 25: 23–46.
[9]
Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, et al. (2010) De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics 11: 400.
[10]
Horowitz AR, Antignus Y, Gerling D (2011) Management of Bemisia tabaci Whiteflies. Springer, Netherlands. 293–322.
[11]
Nauen R, Konanz S (2005) Spiromesifen as a new chemical option for resistance management in whiteflies and spider mites. Pflanzenschutz-Nachr 58: 485–502.
[12]
Nauen R, Reckmann U, Thomzik J, Thielert W (2008) Biological profile of spirotetramat (Movento?)–a new two-way systemic (ambimobile) insecticide against sucking pest species. Bayer CropSci J 61: 245–278.
[13]
Sattelle DB, Cordova D, Cheek TR (2008) Insect ryanodine receptors: molecular targets for novel pest control chemicals. Invertebr Neurosci 8: 107–119.
[14]
Vassiliou V, Emmanouilidou M, Perrakis A, Morou E, Vontas J, et al. (2011) Insecticide resistance in Bemisia tabaci from Cyprus. Insect Sci 18: 30–39.
[15]
Farenhorst M, Mouatcho JC, Kikankie CK, Brooke BD, Hunt RH, et al. (2009) Fungal infection counters insecticide resistance in African malaria mosquitoes. Proc Natl Acad Sci U S A 106: 17443–17447.
[16]
Faria M, Wraight SP (2001) Biological control of Bemisia tabaci with fungi. Crop Prot 20: 767–778.
[17]
Hajek A, St. Leger R (1994) Interactions between fungal pathogens and insect hosts. Annu Rev Entomol 39: 293–322.
[18]
Saito T, Sugiyama K (2005) Pathogenicity of three Japanese strains of entomopathogenic fungi against the silverleaf whitefly, Bemisia argentifolii. Appl Entomol Zool 40: 169–172.
[19]
Daniel C, Wyss E (2010) Field applications of Beauveria bassiana to control the European cherry fruit fly Rhagoletis cerasi. J Appl Entomol 134: 675–681.
[20]
Feng M, Poprawski T, Khachatourians G (1994) Production, formulation and application of the entomopathogenic fungus Beauveria bassiana for insect control: current status. Biocontrol Sci Technol 4: 3–34.
[21]
Torrado-Leon E, Montoya-Lerma J, Valencia-Pizo E (2006) Sublethal effects of Beauveria bassiana (Balsamo) Vuillemin (Deuteromycotina: Hyphomycetes) on the whitefly Bemisia tabaci (Gennadius)(Hemiptera: Aleyrodidae) under laboratory conditions. Mycopathologia 162: 411–419.
[22]
Wraight S, Carruthers R, Jaronski S, Bradley C, Garza C, et al. (2000) Evaluation of the entomopathogenic fungi Beauveria bassiana and Paecilomyces fumosoroseus for microbial control of the silverleaf whitefly, Bemisia argentifolii. Biol Control 17: 203–217.
[23]
Toledo AV, de Remes Lenicov AMM, Lastra CCL (2010) Histopathology caused by the entomopathogenic fungi, Beauveria bassiana and Metarhizium anisopliae, in the adult planthopper, Peregrinus maidis, a maize virus vector. J Insect Sci 10: 35.
[24]
Hou C, Qin G, Liu T, Mei X, Zhang R, et al. (2011) Differential gene expression in silkworm in response to Beauveria bassiana infection. Gene 484: 35–41.
[25]
Qin Y, Ying SH, Chen Y, Shen ZC, Feng MG (2010) Integration of insecticidal protein Vip3Aa1 into Beauveria bassiana enhances fungal virulence to Spodoptera litura larvae by cuticle and per os infection. Appl Environ Microbiol 76: 4611–4618.
[26]
Fan Y, Fang W, Guo S, Pei X, Zhang Y, et al. (2007) Increased insect virulence in Beauveria bassiana strains overexpressing an engineered chitinase. Appl Environ Microbiol 73: 295–302.
[27]
Lu D, Pava-Ripoll M, Li Z, Wang C (2008) Insecticidal evaluation of Beauveria bassiana engineered to express a scorpion neurotoxin and a cuticle degrading prot ease. Appl Microbiol and Biotech 81: 515–522.
[28]
Wang XW, Luan JB, Li JM, Su YL, Xia J, et al. (2011) Transcriptome analysis and comparison reveal divergence between two invasive whitefly cryptic species. BMC Genomics 12: 458.
[29]
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y (2008) RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res 18: 1509–1517.
[30]
Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63.
[31]
Shi WB, Feng MG (2009) Effect of fungal infection on reproductive potential and survival time of Tetranychus urticae (Acari: Tetranychidae). Exp Appl Acarol 48: 229–237.
[32]
Chen J, Xie C, Tian L, Hong L, Wu X, et al. (2010) Participation of the p38 pathway in Drosophila host defense against pathogenic bacteria and fungi. Proc Natl Acad Sci U S A 107: 20774–20779.
[33]
Irving P, Troxler L, Heuer TS, Belvin M, Kopczynski C, et al. (2001) A genome-wide analysis of immune responses in Drosophila. Proc Natl Acad Sci U S A 98: 15119–15124.
[34]
Luan JB, Li JM, Varela N, Wang YL, Li FF, et al. (2011) Global analysis of the transcriptional response of whitefly to Tomato yellow leaf curl China virus reveals the relationship of coevolved adaptations. J Virol 85: 3330–3340.
[35]
Xiao G, Ying SH, Zheng P, Wang ZL, Zhang S, et al. (2012) Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci. Rep. 2: 483.
[36]
Veitch NJ, Johnson PCD, Trivedi U, Terry S, Wildridge D, et al. (2010) Digital gene expression analysis of two life cycle stages of the human-infective parasite, Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes. BMC Genomics 11: 124.
[37]
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, et al. (2011) Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature 474: 380–384.
[38]
Bos JIB, Prince D, Pitino M, Maffei ME, Win J, et al. (2010) A functional genomics approach identifies candidate effectors from the aphid species Myzus persicae (green peach aphid). PLoS Genet 6: e1001216.
[39]
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔC(T) method. Methods 25: 402–408.
[40]
Hoen PAC, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RHAM, et al. (2008) Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res 36: e141.
[41]
Hegedus Z, Zakrzewska A, ágoston VC, Ordas A, Rácz P, et al. (2009) Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol 46: 2918–2930.
[42]
VanGuilder HD, Vrana KE, Freeman WM (2008) Twenty-five years of quantitative PCR for gene expression analysis. BioTechniques 44: 619–626.
[43]
Mao X, Cai T, Olyarchuk JG, Wei L (2005) Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21: 3787–3793.
[44]
Yin Y, Martin J, Abubucker S, Scott AL, McCarter JP, et al. (2008) Intestinal transcriptomes of nematodes: comparison of the parasites Ascaris suum and Haemonchus contortus with the free-living Caenorhabditis elegans. PLoS Negl Trop Dis 2: e269.
[45]
Ciccia A, Elledge SJ (2010) The DNA damage response: making it safe to play with knives. Mol Cell 40: 179–204.
[46]
Toller IM, Neelsen KJ, Steger M, Hartung ML, Hottiger MO, et al. (2011) Carcinogenic bacterial pathogen Helicobacter pylori triggers DNA double-strand breaks and a DNA damage response in its host cells. Proc Natl Acad Sci U S A 108: 14944–14949.
[47]
Song J, Durrant WE, Wang S, Yan S, Tan EH, et al. (2011) DNA repair proteins are directly involved in regulation of gene expression during plant immune response. Cell Host Microbe 9: 115–124.
[48]
Wu M, Huang H, Zhang W, Kannan S, Weaver A, et al. (2011) Host DNA Repair Proteins in Response to Pseudomonas aeruginosa in Lung Epithelial Cells and in Mice. Infect Immun 79: 75–87.
[49]
Li X, Schuler MA, Berenbaum MR (2007) Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol 52: 231–253.
[50]
Aguilar R, Jedlicka AE, Mintz M, Mahairaki V, Scott AL, et al. (2005) Global gene expression analysis of Anopheles gambiae responses to microbial challenge. Insect Biochem Mol Biol 35: 709–719.
[51]
Huang L, Cheng T, Xu P, Cheng D, Fang T, et al. (2009) A genome-wide survey for host response of silkworm, Bombyx mori during pathogen Bacillus bombyseptieus infection. PLoS ONE 4: e8098.
[52]
Zibaee A, Bandani AR, Tork M (2009) Effect of the entomopathogenic fungus, Beauveria bassiana, and its secondary metabolite on detoxifying enzyme activities and acetylcholinesterase (AChE) of the Sunn pest, Eurygaster integriceps (Heteroptera: Scutellaridae). Biocontrol Science and Technology 19: 485–498.
[53]
Clarkson JM, Charnley AK (1996) New insights into the mechanisms of fungal pathogenesis in insects. Trends in Microbiol 4: 197–203.
[54]
Serebrov V, Gerber O, Malyarchuk A, Martemyanov V, Alekseev A, et al. (2006) Effect of entomopathogenic fungi on detoxification enzyme activity in greater wax moth Galleria mellonella L.(Lepidoptera, Pyralidae) and role of detoxification enzymes in development of insect resistance to entomopathogenic fungi. Biology Bulletin 33: 581–586.
[55]
Nelson N (1999) Metal ion transporters and homeostasis. EMBO J 18: 4361–4371.
[56]
Perraud AL, Knowles H, Schmitz C (2004) Novel aspects of signaling and ion-homeostasis regulation in immunocytes: The TRPM ion channels and their potential role in modulating the immune response. Mol Immunol 41: 657–673.
[57]
Krattinger SG, Lagudah ES, Spielmeyer W, Singh RP, Huerta-Espino J, et al. (2009) A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat. Science 323: 1360.
[58]
Tanabe M, Atkins HS, Harland DN, Elvin SJ, Stagg AJ, et al. (2006) The ABC transporter protein OppA provides protection against experimental Yersinia pestis infection. Infect Immun 74: 3687.
[59]
Kobae Y, Sekino T, Yoshioka H, Nakagawa T, Martinoia E, et al. (2006) Loss of AtPDR8, a plasma membrane ABC transporter of Arabidopsis thaliana, causes hypersensitive cell death upon pathogen infection. Plant Cell Physiol 47: 309–318.
[60]
Dong C, Davis RJ, Flavell RA (2002) MAP kinases in the immune response. Annu Rev Immunol 20: 55–72.
[61]
Scanga CA, Bafica A, Feng CG, Cheever AW, Hieny S, et al. (2004) MyD88-deficient mice display a profound loss in resistance to Mycobacterium tuberculosis associated with partially impaired Th1 cytokine and nitric oxide synthase 2 expression. Infect Immun 72: 2400–2404.
[62]
Tauszig-Delamasure S, Bilak H, Capovilla M, Hoffmann JA, Imler JL (2001) Drosophila MyD88 is required for the response to fungal and Gram-positive bacterial infections. Nat Immunol 3: 91–97.
[63]
Dong Y, Morton JC, Ramirez JL, Souza-Neto JA, Dimopoulos G (2012) The entomopathogenic fungus Beauveria bassiana activate toll and JAK-STAT pathway-controlled effector genes and anti-dengue activity in Aedes aegypti. Insect Biochem Mol Biol 42: 126–132.
[64]
Pendland J, Hung S, Boucias D (1993) Evasion of host defense by in vivo-produced protoplast-like cells of the insect mycopathogen Beauveria bassiana. J Bacteriol 175: 5962–5969.
[65]
Boyce JD, Harper M, Michael FS, John M, Aubry A, et al. (2009) Identification of novel glycosyltransferases required for assembly of the Pasteurella multocida A: 1 lipopolysaccharide and their involvement in virulence. Infect Immun 77: 1532–1542.
[66]
Monge RA, Roman E, Nombela C, Pla J (2006) The MAP kinase signal transduction network in Candida albicans. Microbiology 152: 905–912.
[67]
Dunn MF, Ramírez-Trujillo J, Hernández-Lucas I (2009) Major roles of isocitrate lyase and malate synthase in bacterial and fungal pathogenesis. Microbiology 155: 3166–3175.
[68]
Liu P, Wood D, Nester EW (2005) Phosphoenolpyruvate carboxykinase is an acid-induced, chromosomally encoded virulence factor in Agrobacterium tumefaciens. J Bacteriol 187: 6039–6045.
[69]
Mu?oz-Elías EJ, McKinney JD (2005) Mycobacterium tuberculosis isocitrate lyases 1 and 2 are jointly required for in vivo growth and virulence. Nat Med 11: 638–644.
[70]
Sch?r J, Stoll R, Schauer K, Loeffler DIM, Eylert E, et al. (2010) Pyruvate carboxylase plays a crucial role in carbon metabolism of extra-and intracellularly replicating Listeria monocytogenes. J Bacteriol 192: 1774–1784.
[71]
Padilla-Guerrero IE, Barelli L, González-Hernández GA, Torres-Guzmán JC, Bidochka MJ (2011) Flexible metabolism in Metarhizium anisopliae and Beauveria bassiana: role of the glyoxylate cycle during insect pathogenesis. Microbiology 157: 199–208.
[72]
Wall DM, Duffy PS, DuPont C, Prescott JF, Meijer WG (2005) Isocitrate lyase activity is required for virulence of the intracellular pathogen Rhodococcus equi. Infect Immun 73: 6736–6741.
[73]
Piper PW, Ortiz-Calderon C, Holyoak C, Coote P, Cole M (1997) Hsp30, the integral plasma membrane heat shock protein of Saccharmyces cerevisiae, is a stress-inducible regulator of plasma membrane H+-ATPase. Cell Stress Chaperones 2: 12–24.
[74]
Collette JR, Lorenz MC (2011) Mechanisms of immune evasion in fungal pathogens. Curr Opin Microbiol14: 668–675.
[75]
Stergiopoulos I, de Wit PJGM (2009) Fungal effector proteins. Annu Rev Phytopathol 47: 233–263.
[76]
Mayer AM, Staples RC (2002) Laccase: new functions for an old enzyme. Phytochem 60: 551–565.
[77]
Fang W, Fernandes éK, Roberts DW, Bidochka MJ, St Leger RJ (2010) A laccase exclusively expressed by Metarhizium anisopliae during isotropic growth is involved in pigmentation, tolerance to abiotic stresses and virulence. Fungal Genet and Biol 47: 602–607.
[78]
Gerardo NM, Altincicek B, Anselme C, Atamian H, Barribeau SM, et al. (2010) Immunity and other defenses in pea aphids, Acyrthosiphon pisum. Genome Biol 11: R21.
[79]
Alleyne M, Chappell MA, Gelman DB, Beckage NE (1997) Effects of parasitism by the braconid wasp Cotesia congregata on metabolic rate in host larvae of the tobacco hornworm, Manduca sexta. J Insect Physiol 43: 143–154.
[80]
Altincicek B, Gross J, Vilcinskas A (2008) Wounding-mediated gene expression and accelerated viviparous reproduction of the pea aphid Acyrthosiphon pisum. Insect Mol Biol 17: 711–716.
[81]
Gwynn D, Callaghan A, Gorham J, Walters K, Fellowes M (2005) Resistance is costly: trade-offs between immunity, fecundity and survival in the pea aphid. Proc R Soc B 272: 1803.
[82]
Li JM, Ruan YM, Li FF, Liu SS, Wang XW (2011) 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 18: 11–22.
[83]
Karatolos N, Pauchet Y, Wilkinson P, Chauhan R, Denholm I, et al. (2011) Pyrosequencing the transcriptome of the greenhouse whitefly, Trialeurodes vaporariorum reveals multiple transcripts encoding insecticide targets and detoxifying enzymes. BMC Genomics 12: 56.
[84]
Leshkowitz D, Gazit S, Reuveni E, Ghanim M, Czosnek H, et al. (2006) Whitefly (Bemisia tabaci) genome project: analysis of sequenced clones from egg, instar, and adult (viruliferous and non-viruliferous) cDNA libraries. BMC Genomics 7: 79.
[85]
Xie W, Meng QS, Wu QJ, Wang SL, Yang X, et al. (2012) Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance. PLoS ONE 7: e35181.
[86]
Gao Q, Jin K, Ying SH, Zhang Y, Xiao G, et al. (2011) Genome sequencing and comparative transcriptomics of the model entomopathogenic fungi Metarhizium anisopliae and M. acridum. PLoS Genet 7: e1001264.
[87]
Wang C, Leger RJS (2005) Developmental and transcriptional responses to host and nonhost cuticles by the specific locust pathogen Metarhizium anisopliae var. acridum. Eukaryot Cell 4: 937–947.