key: cord-0011778-7ygdbk5a authors: Li, Meng; Zhao, Jing; Su, Yun-Lin title: Transcriptome Analysis of Gene Expression Profiles of Tomato Yellow Leaf Curl Virus-Infected Whiteflies over Different Viral Acquisition Access Periods date: 2020-05-11 journal: Insects DOI: 10.3390/insects11050297 sha: 60967d96c232cdc1fdfb12392f09811eafbc7117 doc_id: 11778 cord_uid: 7ygdbk5a Tomato yellow leaf curl virus (TYLCV), which is transmitted by Bemisia tabaci in a persistent-circulative manner, threatens tomato production worldwide. Little is known about the complicated interaction during this process at the molecular level. In this study, viral AAPs at 0 h, 2 h, 6 h, 12 h and 48 h were investigated using a comparative transcriptome analysis to uncover the transcriptional responses of whiteflies to virus infection. Our results have shown that 755, 587, 1140 and 1347 differentially expressed genes (DEGs) were identified in the comparisons of the data of 0 h vs. 2 h, 0 h vs. 6 h, 0 h vs. 12 h and 0 h vs. 48 h, respectively. KEGG analysis showed that DEGs associated with metabolisms and signal transduction were down-regulated in virus-infected whiteflies. Additionally, 16 up-regulated putative transporter genes and 10 down-regulated genes associated with IL-17 signaling pathway were identified by time-associated gene cluster analysis. These data boost our comprehensions on whitefly-TYLCV interactions associated with different viral AAPs. Plant virus diseases threaten crop yield and cause tremendous losses of agricultural crops worldwide [1] [2] [3] . In general, plant viruses are transmitted by insect vectors in four manners as follows: non-persistent, semi-persistent, persistent-circulative and persistent-propagative [4] [5] [6] . Begomoviruses (genus Begomovirus, family Geminiviridae) are transmitted in a persistent-circulative manner by the whitefly Bemisia tabaci [7] [8] [9] [10] . They are a group of destructive plant viruses that affect many important crops such as tomato [2] and cassava [11] . Tomato yellow leaf curl virus (TYLCV) is one of the most devastating begomoviruses in the 'top 10' plant viruses list [2] . Among the known 39 cryptic species of B. tabaci [12] [13] [14] , as far as we know, there are only two invasive species, Middle East Asia Minor 1 (MEAM1) and Mediterranean (MED), that acquire and transmit TYLCV with high efficiency [15] . After being ingested by whiteflies, TYLCV virions can travel from the midgut lumen into the haemolymph, then move to the salivary gland with the flow of haemolymph and finally be inoculated into healthy plants with saliva excretion [5] . In this process, numerous genes are regulated in the whiteflies. Li et al. [16] and Luan et al. [17] suggested that acquisition of Tomato yellow leaf curl China virus (TYLCCNV) influenced the expression of many genes related to metabolism, cell cycle and immune reactions in whiteflies. Hasegawa et al. [18] and Geng et al. [19] also demonstrated that the infection of TYLCV induced transcriptional changes of the genes involved in cell cycle and antiviral immunity. Furthermore, expression of the genes of receptors, trafficking proteins, cuticle proteins and defense related proteins could also be interfered with by begomoviruses in vectors [18] [19] [20] . Previous study indicated that TYLCV could accomplish the circulation from midgut to salivary gland within 8 h [21] . Furthermore, viruses could be detected in midgut, haemolymph and salivary glands in whiteflies after an acquisition access period (AAP) of 1 h, 2 h and 7 h, respectively [21, 22] . After a viral AAP of 12-48 h, the amount of viruses in whiteflies reached a very high level of abundance [23] . Gene expression of whiteflies were extensively studied to unravel the mechanisms of whitefly-virus interaction. However, as far as we know, those studies only conducted the comparisons of gene expression profiles between the whiteflies feeding on healthy plants and those feeding on virus-infected plants over 24 h [16-20] . Therefore, expression patterns of whitefly genes in the early stage of virus-vector interaction and the temporal regulation in whiteflies with the gradual accumulation of virus in vectors remains unclear. In this study, viral AAPs at 0 h, 2 h, 6 h, 12 h and 48 h were used for comparative transcriptome analysis to further explore the mechanisms in virus-vector interaction. Our work sought to answer the following two questions: (1) regardless of viral amounts, does gene regulation change in TYLCV-infected whiteflies; (2) does any temporally up-regulated or down-regulated genes or pathways exist in TYLCV-infected whiteflies? These answers may help to reveal the begomeovirus-whitefly interactions, especially in the early stage, and to identify novel vector factors in the transmission of TYLCV. MEAM1 whitefly colony (mtCOI GneBank Accession code: GQ332577) was maintained on cotton plants (Gossypium hirsutum L. cv. Zhemian 1793) in insect proof cages. TYLCV-infected tomato (Solanum lycopersicum Mill. cv. Hezuo903) plants were obtained using the procedures as described before [15] . Five-to seven-day-old MEAM1 adults were collected and transferred to TYLCV-infected tomato plants for different periods. Insect rearing, plant cultivation and sample collection were conducted at 26 ± 1 • C, a 14:10 photo: scoto period and a relative humidity of 60 ± 10%. To determine the acquisition efficiency and virus titers in whiteflies after different AAPs on TYLCV-infected plants, approximately 500 whitefly adults were released onto TYLCV-infected tomato plants for 1 h, 2 h, 6 h, 8 h, 12 h, 24 h and 48 h. Following each AAP, 23-24 female adults were collected and assayed individually for the detection of TYLCV DNA by PCR as described in [15] . Whiteflies with an AAP of 2 h, 6 h, 12 h and 48 h were subjected to DNA extraction for quantitative RT-PCR (qPCR) analysis, as described previously [24] . The primers used in PCR/qPCR analyses were listed in Table 1 . As to RNA-seq analysis, each AAP (2 h, 6 h, 12 h and 48 h) consisted of three biological replicates and each biological replicate included 200 whitefly adults. Whiteflies with no access to virus-infected plants were used as control (0 h). In total, 15 samples (3 replicates × 5 AAPs) were collected for transcriptome sequencing. Total RNA of each sample was extracted using Trizol (Invitrogen, Carlsbad, CA, USA) according to the manual instructions. Genomic DNA contamination was removed using RNase-free DNaseI (Takara, Kusatsu, Shiga, Japan). The quantity and quality of total RNA were evaluated by Nano Drop (Thermo Fisher Scientific, Waltham, MA, USA) and Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), respectively. When the RNA integrity number (RIN) of a sample was greater than six, 1000 ng of total RNA was used as input. The total RNA was firstly treated by Oligo (dT)-attached magnetic beads for mRNA purification [25] . Then, cDNA library was constructed as follows: (1) fragmenting mRNA into small piece; (2) generating first-strand cDNA with random hexamer-primed reverse transcription, followed by a second-strand cDNA synthesis; (3) adding A-Tailing Mix AND RNA Index adapters to end repair by incubation; (4) amplifying the cDNA fragments obtained from last steps by PCR, purifying PCR products by Ampure XP Beads, and dissolving the PCR products in EB solution; (5) validating products on the Agilent Technologies 2100 bioanalyzer for quality control; (6) denaturing the double strand PCR products from step (5) by heating and circularizing by the splint oligo sequence to produce the final library; (7) amplifying the final library with phi29 to make DNA nanoball (DNB), which had more than 300 copies of one molecule. DNBs were loaded into nanoarray and then sequenced for 100 bp paired-end reads on BGISEQ-500 platform (BGI-Shenzhen, China). The raw reads generated from sequencing were cleaned by removing adaptor-polluted reads, reads with unknown sequences "N" accounting for more than 5%, and low-quality reads (more than 20% bases with Phred threshold score < 10). The clean reads were mapped to the MEAM1 whitefly genome (http://www.whiteflygenomics.org) using the software of Hierarchical indexing for spliced alignment of transcripts (HISAT) v2.1.0 [26] . Reads per kilobase million mapped reads (RPKM) was used to estimate gene expression levels. It could eliminate the bias caused by effects of sequencing depth and gene length. Fold changes of gene expressions were analyzed by the software of DEGseq [27] and differentially expressed genes (DEGs) were identified with the threshold of Q-value (adjusted p-value) ≤0.001 and log2 (fold change) ≥1. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) functional classifications were performed using Blast2GO and KOBAS (http://kobas.cbi.pku.edu.cn). KEGG enrichment analysis was conducted by using R phyer. p-value was corrected by false discovery rate (FDR). The pathways with FDR <0.01 was considered to be significantly enriched. Data mining and figure presentation process were analyzed by BGI in-house customized data mining system called Dr.Tom (http://biosys.bgi.com). Ten DEGs were selected for qPCR analysis to validate the data obtained by RNA-seq. Thirty whitefly adults were collected as one sample at each AAP. Three replicates were performed. Total RNA was extracted using Trizol (Invitrogen, Carlsbad, CA, USA) following the manual protocol. Reverse transcription was conducted by using PrimeScript RT reagent Kit (TaKaRa, Kusatsu, Shiga, Japan). The qPCR was performed on the CFX96 TM Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) with SYBR Premix Ex ExTaq II (Takara, Kusatsu, Shiga, Japan). We used β-actin as a reference gene. Relative expression level of DEGs was calculated by 2 −∆Ct . Primers used in this experiment were listed in Table 1 . As shown in Figure 1 , the percentage of the whiteflies infected with TYLCV increased with the AAP up to 48 h. Additionally, the raised viruses quantity in whiteflies paralleled the increasing viral feeding time. To investigate the gene expression patterns, the whiteflies after an AAP of 0 h, 2 h, 6 h, 12 h and 48 h were subjected to RNA-seq analysis. The relevant data were deposited in the SRA as BioProject PRJNA562795. Fifteen cDNA libraries were constructed and~21 M clean reads per library were obtained from the transcriptome sequencing (Table S1 ). On average, 91.81% and 77.95% of these reads were mapped and uniquely mapped to the MEAM1 whitefly genome, (http://www.whiteflygenomics.org), respectively. Furthermore, the reads were mapped to 15,662 genes of the whitefly genome. Pearson's correlation coefficients indicated the high reproduce of the data in different replications ( Figure S1 ). To reveal the regulation of gene expression upon acquisition of TYLCV, we compared the data in four ways, 0 h vs. (Figure 2 ). Among the DEGs, the percentage of down-regulated genes increased, while the up-regulated genes decreased throughout the entire AAPs ( Figure 2 ). This suggested that TYLCV had significant influence on the gene expression of its vector, and, moreover, the effect became stronger as the virus amount increased in the whitefly population. GO assignments were performed to reveal the putative functions of DEGs. DEGs were categorized into 38 secondary GO categories in three biological concepts including biological process, cellular component and molecular function (Figure 3 ). Cellular process, membrane and catalytic activity were the most presented categories in biological process, cellular component and molecular function, respectively ( Figure 3 ). The analysis of GO categories suggested that the DEGs To reveal the regulation of gene expression upon acquisition of TYLCV, we compared the data in four ways, 0 h vs. Figure 2 ). Among the DEGs, the percentage of down-regulated genes increased, while the up-regulated genes decreased throughout the entire AAPs ( Figure 2 ). This suggested that TYLCV had significant influence on the gene expression of its vector, and, moreover, the effect became stronger as the virus amount increased in the whitefly population. To reveal the regulation of gene expression upon acquisition of TYLCV, we compared the data in four ways, 0 h vs. (Figure 2 ). Among the DEGs, the percentage of down-regulated genes increased, while the up-regulated genes decreased throughout the entire AAPs ( Figure 2 ). This suggested that TYLCV had significant influence on the gene expression of its vector, and, moreover, the effect became stronger as the virus amount increased in the whitefly population. GO assignments were performed to reveal the putative functions of DEGs. DEGs were categorized into 38 secondary GO categories in three biological concepts including biological process, cellular component and molecular function (Figure 3 ). Cellular process, membrane and catalytic activity were the most presented categories in biological process, cellular component and molecular function, respectively (Figure 3 ). The analysis of GO categories suggested that the DEGs GO assignments were performed to reveal the putative functions of DEGs. DEGs were categorized into 38 secondary GO categories in three biological concepts including biological process, cellular component and molecular function (Figure 3 ). Cellular process, membrane and catalytic activity were the most presented categories in biological process, cellular component and molecular function, respectively ( Figure 3 ). The analysis of GO categories suggested that the DEGs in TYLCV-infected whiteflies at different viral AAPs had similar functional classifications in the biological process, cellular component and molecular function (Figure 3 ). A total of 235 common DEGs were found among different comparisons ( Figure S2 ). Hierarchical cluster analysis generated by pheatmap R indicated a consistent up-or down-regulation of these common DEGs among different comparisons ( Figure S2 ). The result of the KEGG analysis suggested Insects 2020, 11, 297 7 of 14 that 118 common down-regulated DEGs were related to metabolisms, signal transduction, organismal systems, genetic or environmental information processing, cellular process and human diseases ( Figure 4 ). Among the down-regulated DGEs, a large proportion of them were classified in metabolism groups, especially for lipid metabolism, carbohydrate metabolism and amino acid ( Figure 4) . As for the category of signal transduction, 27 common down-regulated DEGs were found to be implicated in 13 signaling pathways such as AMPK signaling pathway, ras signaling pathwaysig, insulin signaling pathway and others ( Figure 4 , Table S2 ). Insects 2020, 11, x 7 of 14 A total of 235 common DEGs were found among different comparisons ( Figure S2 ). Hierarchical cluster analysis generated by pheatmap R indicated a consistent up-or downregulation of these common DEGs among different comparisons ( Figure S2 ). The result of the KEGG analysis suggested that 118 common down-regulated DEGs were related to metabolisms, signal transduction, organismal systems, genetic or environmental information processing, cellular process and human diseases (Figure 4 ). Among the down-regulated DGEs, a large proportion of them were classified in metabolism groups, especially for lipid metabolism, carbohydrate metabolism and amino acid ( Figure 4) . As for the category of signal transduction, 27 common down-regulated DEGs were found to be implicated in 13 signaling pathways such as AMPK signaling pathway, ras signaling pathwaysig, insulin signaling pathway and others ( Figure 4 , Table S2 ). To examine the temporal shift in gene expression during viral acquisition, time-associated gene cluster of all 15,662 genes was conducted by Mfuzz v2.34.0 [25, 28] . It was notable that 1168 genes were found to be continually up-regulated ( Figure 5A ) and 16 genes encoding transporters were identified ( Table 2 ). In addition, 1312 genes were continually down-regulated ( Figure 5B ) and 10 genes associated with interleukin (IL)-17 signaling pathway were identified (Table 3) . Hierarchical cluster analysis of 127 up-regulated DEGs and 355 down-regulated DEGs in the comparison of 0 h vs. 48 h was shown in Figure S3 . Insects 2020, 11, x 8 of 14 To examine the temporal shift in gene expression during viral acquisition, time-associated gene cluster of all 15,662 genes was conducted by Mfuzz v2.34.0 [25, 28] . It was notable that 1168 genes were found to be continually up-regulated ( Figure 5A ) and 16 genes encoding transporters were identified ( Table 2 ). In addition, 1312 genes were continually down-regulated ( Figure 5B ) and 10 genes associated with interleukin (IL)-17 signaling pathway were identified (Table 3) . Hierarchical cluster analysis of 127 up-regulated DEGs and 355 down-regulated DEGs in the comparison of 0 h vs. 48 h was shown in Figure S3 . Expression levels of 10 DEGs, i.e., 4 up-regulated DEGs ( Figure 6A ) and 6 down-regulated DEGs ( Figure 6B ), were verified by qPCR. The results showed that no significant difference in the trends of gene regulation compared to the results generated from RNA-seq analysis, suggesting that the transcriptome data were reliable in this study. Several studies had analyzed the gene expression profiles of TYLCV-infected whiteflies (with the AAPs longer than 24 h) using high-throughput RNA-sequencing and comparative transcriptome analysis [17] [18] [19] [20] . To reveal more details in the early-stage interactions between whiteflies and TYLCV, gene expression profiles of TYLCV-infected whiteflies with the AAPs less Several studies had analyzed the gene expression profiles of TYLCV-infected whiteflies (with the AAPs longer than 24 h) using high-throughput RNA-sequencing and comparative transcriptome analysis [17] [18] [19] [20] . To reveal more details in the early-stage interactions between whiteflies and TYLCV, gene expression profiles of TYLCV-infected whiteflies with the AAPs less than 12 h were investigated in this study. Our results indicated that hundreds of genes were up-or down-regulated as a consequence of viral infection. Numbers of DEGs tend to increase with the duration of AAPs, which may be due to the accumulation of TYLCV in whitefly populations. Interestingly, as the percentage of down-regulated genes increased, while the up-regulated genes decreased during the extending AAPs, we deduced that the repression of vector genes might facilitate the infection of TYLCV in B. tabaci. It should be pointed out that, in addition to the response to virus infection, the changes in gene expressions of whiteflies may also be caused by the adaptation to new hosts [29] . In particular, we noticed that the number of DEGs identified in the comparison of 0 h vs. 2 h (755) is higher than that of 0 h vs. 6 h (587). The number of DEGs identified by comparing expression profiles of whiteflies with different AAPs with that of control (0 h) may be higher than the actual. As was previously observed in [18] , our results showed that the genes clustered in the groups of the metabolisms involving lipid, carbohydrate and amino acid metabolism were down-regulated in virus-infected whiteflies. Similarly, disturbances of amino acid metabolisms in whitefly were found to be associated with TYLCCNV infection [17] and manipulations of metabolisms by plant viruses were reported in other insect vectors [30] . TYLCV infection of MEAM1 reduced the whitefly's fecundity and longevity [31] . Longevity of MED was also found to be shortened by TYLCV infection [32] . These deleterious effects of TYLCV could be caused by the dramatic repression of the genes involved in basic metabolisms including lipid, carbohydrate and amino acid metabolisms. On the other hand, interference of signaling pathway is a universal phenomenon in vector-virus interactions [18, 33] . Herein, 13 signaling pathways were predicted to be down-regulated in TYLCVinfected whiteflies. Genes involved in the 20 signaling pathway were also found to be down-regulated in Tomato chlorosis virus-infected whiteflies [33] . Among the down-regulated pathways, AMPK and PI3K-Akt signaling pathways are associated with cellular apoptosis and autophagy [34] [35] [36] [37] [38] [39] . The regulation of PI3K-Akt signaling pathway by animal viruses have been extensively studied in mammals [40] . Viruses may stimulate [41] [42] [43] [44] or inhibit [45, 46] PI3K-Akt signaling pathway to reduce or activate apoptosis in host cells to facilitate their infection. However, roles of pi3k-akt signaling pathway in the interactions between the viruses and whiteflies remain unknown. Time-associated gene cluster analysis revealed the continuous down-regulation of the genes involved in IL-17 signaling pathway over the AAPs. IL-17 signaling has been shown to play critical role in host defense during bacterial infection [47] . However, as for virus, two studies reported that the IL-17 acted as a pathogenic factor in human [48, 49] . The function of IL-17 signaling in viral infection to insect vector was unclear to date due to insufficient study. The suppression of IL-17 signaling pathway by TYLCV might assist the persistence of virus through the suppression of the down-stream immune reaction in the vector. After the long-term coadaption and coevolution of begomoviruses and whiteflies, whiteflies have developed defenses systems to prevent virus invasion, and meanwhile, viruses have evolved the capacity of utilizing cellular proteins in whitefly to assist their transmission. Viral transportation plays an important role in virus transmission, and whitefly cellular proteins acts complicatedly in the interplays between whitefly and begomoviruses. Up to now, several whitefly proteins had been uncovered their favors in viral transport [50] [51] [52] [53] [54] [55] . A midgut protein, cyclophilin B and collagen, has been characterized to act positively in viral transport [50] [51] [52] . On the contrary, heat shock protein 70, knottin-1 and vesicle associated membrane protein-associated protein B have been proved to cause the inhibition in viral transport [53] [54] [55] . Additionally, the increasing omics-data uncover more and more whitefly proteins involved in viral circulative transport. For example, a large number of receptors and transporters, including sugar transporter ERD6-like 6, proton-coupled amino acid transporter, facilitated trehalose transporter, genes with cargo receptor activity, ATP-binding cassette (ABC) transporters, cuticle proteins, laminin subunit alpha, dystroglycan and integrin alpha-PS2, were found to be associated with viral transport in whiteflies [18] [19] [20] . In this study, in addition to facilitated trehalose transporter, some uncharacterized proteins predicted to be with transporter activity were identified in virus-infected whiteflies. These genes might also be crucial in viral transport. The functions of these identified genes/proteins in viral transmission could be further characterized by RNAi. In summary, we uncovered the transcriptional responses of B. tabaci after feeding on TYLCV-infected plants for 2 h, 6 h, 12 h and 48 h. Our results indicated that hundreds of genes were up-or down-regulated as a consequence of viral infection. Genes involved in metabolisms and signaling transduction were down-regulated in virus-infected whiteflies. Furthermore, 16 up-regulated putative transporter genes and 10 down-regulated genes associated with IL-17 signaling pathway were identified by time-associated gene cluster analysis. This study provided new clues to unravel the mechanisms of whitefly-begomovirus interactions. Table S1 : Summary of RNA-Seq datasets generated from whiteflies fed on TYLCV-infected tomato plants for 0, 2, 6, 12 and 48 h, Table S2 : Common DEGs involved in signal transduction in virus-infected whiteflies. Emerging infectious diseases of plants: Pathogen pollution, climate change and agrotechnology drivers Top 10 plant viruses in molecular plant pathology García-Arenal, F. Evolution and emergence of plant viruses Virus-vector interactions mediating nonpersistent and semipersistent transmission of plant viruses Insect vector interactions with persistently transmitted viruses Insect vector-mediated transmission of plant viruses In Distribution and Dissemination of Begomoviruses in Latin America and the Caribbean Emerging virus diseases transmitted by whiteflies Natural genomic and antigenic variation in whitefly-transmitted geminiviruses (begomoviruses) Persistent, circulative transmission of begomoviruses by whitefly vectors Biology and management of Bemisia whitefly vectors of cassava virus pandemics in Africa Bemisia tabaci: A statement of species status Species concepts as applied to the whitefly Bemisia tabaci systematics: How many species are there? Transmission of Tomato yellow leaf curl virus by two invasive biotypes and a Chinese indigenous biotype of the whitefly Bemisia tabaci 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 Global analysis of the transcriptional response of whitefly to Tomato yellow leaf curl China virus reveals the relationship of coevolved adaptations Comparative transcriptome analysis reveals networks of genes activated in the whitefly, Bemisia tabaci, when fed on tomato plants infected with Tomato yellow leaf curl virus Transcriptome profiling of whitefly guts in response to Tomato yellow leaf curl virus infection Comparative proteomic analysis provides new insight into differential transmission of two begomoviruses by a whitefly Rate of Tomato yellow leaf curl virus translocation in the circulative transmission pathway of its vector, the whitefly Bemisia tabaci The circulative pathway of begomoviruses in the whitefly vector Bemisia tabaci-insights from studies with Tomato yellow leaf curl virus The incredible journey of begomoviruses in their whitefly vector Clathrin-mediated endocytosis is involved in Tomato yellow leaf curl virus transport across the midgut barrier of its whitefly vector Transcription profiles reveal the regulatory synthesis of phenols during the development of lotus rhizome (Nelumbo nucifera Gaertn) HISAT: A fast spliced aligner with low memory requirements DEGseq: An R package for identifying differentially expressed genes from RNA-seq data A software package for soft clustering of microarray data Transcriptional responses of invasive and indigenous whiteflies to different host plants reveal their disparate capacity of adaptation Harmful and beneficial effects of plant viruses in insects Long-term association of Tomato yellow leaf curl virus with its whitefly vector Bemisia tabaci: Effect on the insect transmission capacity, longevity and fecundity Upregulation of temperature susceptibility in Bemisia tabaci upon acquisition of Tomato yellow leaf curl virus (TYLCV) Transcriptome analysis of the whitefly, Bemisia tabaci MEAM1 during feeding on tomato infected with the crinivirus, Tomato chlorosis virus, identifies a temporal shift in gene expression and differential regulation of novel orphan genes AMPK and autophagy get connected The AMPK signaling pathway coordinates cell growth, autophagy and metabolism Akt phosphorylation of BAD couples survival signals to the cell-intrinsic death machinery Akt regulates cell survival and apoptosis at a postmitochondrial level Programmed autophagy in the Drosophila fat body is induced by ecdysone through regulation of the PI3K pathway Growth arrest and autophagy are required for salivary gland cell degradation in Drosophila PI3K-Akt signaling and viral infection Regulation of hepatitis B virus replication by the phosphatidylinositol 3-kinase-akt signal transduction pathway Activation of the N-Ras-PI3K-Akt-mTOR pathway by hepatitis C virus, control of cell survival and viral replication Papillomavirus virus-like particles activate the PI3-kinase pathway via alpha-6 beta-4 integrin upon binding Human papillomavirus type 16 E7 up-regulates Akt activity through the etinoblastoma protein VP1 of foot-and-mouth disease virus induces apoptosis via the Akt signaling pathway The SARS coronavirus nucleocapsid protein induces actin reorganization and apoptosis in COS-1 cells in the absence of growth factors IL-17 signaling in host defense and inflammatory diseases Interleukin-17A modulates human airway epithelial responses to human rhinovirus infection Critical role of IL-17RA in immunopathology of influenza infection Bemisia tabaci midgut protein interacts with begomoviruses and plays a role in virus transmission Implication of the whitefly Bemisia tabaci cyclophilin B protein in the transmission of Tomato yellow leaf curl virus Implication of the whitefly, Bemisia tabaci, collagen protein in begomoviruses acquisition and transmission Implication of Bemisia tabaci heat shock protein 70 in begomovirus-whitefly interactions The whitefly Bemisia tabaci knottin-1 gene is implicated in regulating the quantity of Tomato yellow leaf curl virus ingested and transmitted by the insect Implication of whitefly vesicle associated membrane protein-associated protein B in the transmission of Tomato yellow leaf curl virus This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license We thank Yu Dong from Oregon State University, USA for comments on the manuscript. The authors declare no conflict of interest.