Next Generation Sequencing Reveals Regulation of Distinct Aedes microRNAs during Chikungunya Virus Development

Background Application of genomics and Next Generation sequencing has led to the identification of new class of cellular functional molecules, namely, small RNAs. Of the several classes of ncRNAs (non-coding RNA), microRNAs have been demonstrated to exert determinative influence on various cellular processes. It is becoming abundantly clear that host/vector/pathogen encoded microRNAs impact eventual pathogenesis. In this context, the participation of vector based microRNAs in disease transmission and pathogen development is being investigated intensively. A few studies have highlighted the role of vector encoded microRNAs in pathogen infection. We conducted this study to evaluate the role of host miRNAs upon CHIKV (Chikungunya Virus) infection in an important vector, Aedes albopictus. Findings We identified 88 and 79 known miRNAs in uninfected and CHIKV infected Ae. albopictus Singh's cell line respectively. We further identified nine novel miRNAs in Ae. albopictus. Comparison of the two libraries revealed differential expression of 77 common miRNAs between them. CHIKV infection specifically altered the miRNA profile of a specific set of eight miRNAs. Putative targets of these regulated miRNAs were identified and classified into their pathways. Conclusions In our study we have identified and described the profiles of various miRNAs upon CHIKV infection in Ae. albopictus. This investigation provides an insight about cellular modification by miRNAs during CHIKV infection and the results provide leads for identifying potential candidates for vector based antiviral strategies.


Introduction
Vector borne diseases like malaria, dengue, chikungunya are a major burden to the economy of countries in the tropical and subtropical regions. The situation in Asian countries like India are especially alarming with recent studies predicting India to have the maximum global burden of dengue [1] and encompassing nearly 46% of global population at risk for vivax malaria [2]. Currently, no vaccines are available for any of the above mentioned infections and vector control with intense public awareness is the only effective measure for disease control. Transmission blocking approaches provide hope for disease control in laboratory conditions [3] which clearly emphasizes the importance of understanding vector-pathogen dynamics for effective control strategies.
Insects play an important role in disease transmission in plants, animals and humans. Mosquitoes are especially significant in disease transmission in humans. Fundamental biological process of blood feeding in female mosquitoes for egg development is utilized by pathogens either for completing their life cycle or amplifying within the mosquito before infecting the mammalian hosts. Until now, mechanisms involved in disease transmission were obscure. With the recent onslaught of several arboviral epidemics and increase of vector transmitted diseases, several studies, both genetic and non-genetic, have established the molecular mechanisms that may be involved in the host seeking behavior of mosquitoes [4,5], their response to pathogens [6,7] and those factors of pathogen development within the vector affecting transmission [8,9].
RNA interference is the most important phenomena used by insects as the first line of defense against microbial invaders and has been extensively studied in the recent times [10][11][12][13]. Infection of mosquitoes by virus results in generation of different kind of small RNA populations of viral origin, namely 22 nt siRNA population [14], longer nt piRNA population [15,16] as well as modulation of host miRNAs [17][18][19]. MicroRNAs are a class of small (19-24 nt in length), non-coding, single stranded RNAs known to post-transcriptionally regulate gene expression by binding to 39 untranslated regions (39UTRs) of the target mRNAs, resulting in translational repression/cleavage of the cognate mRNAs. First identified in Caenorhabditis elegans, miRNAs have since been identified in a variety of organisms including plants, animals, insects and viruses [20]. Their function primarily is regulation of gene expression and play crucial roles in cell development, proliferation, differentiation, host defense, fat metabolism and apoptosis [21].
We undertook the present study to identify mosquito host miRNAs that may be impacted during chikungunya infection. For this purpose, we performed deep sequencing of small RNA population in Aedes albopictus Singh's cell line upon chikungunya virus infection. We identified several Aedes known and novel miRNAs and profiled their expression during CHIKV infection. In a separate study we had identified microRNAs in Anopheles stephensi mosquitoes and profiled these miRNAs upon Plasmodium infection [22]. During the course of analysis, we recognized there were miRNAs common between Aedes and Anopheles; however, there was discrete regulation of some of the common miRNAs upon infection by CHIKV and Plasmodium. In this report, we describe Ae. Albopictus miRNAs and their regulation upon CHIKV infection. Furthermore, we highlight the repertoire of common and distinct host miRNAs regulated during Plasmodium and CHIKV infections in their respective hosts. We also have predicted the targets of some of these miRNAs in a bid to understand the role these miRNAs may be playing during pathogen infection in mosquitoes.

Cell culture and virus infection
Ae. albopictus cell line Singh's line (ATCC-CCL-126) was maintained at 28uC in DMEM media supplemented with 10% FBS and antibiotics. The cells were infected with MOI = 10 of the wild type chikungunya virus (Accession no. JF950631.1) isolated from India during an outbreak in 2010 [23]. Cells were infected in triplicates and harvested 24 hours post infection and processed further for small RNA sequencing.

Small RNA sequencing
Illumina Truseq small RNA libraries for two samples, namely, Ae. Albopictus Singh's cell line uninfected (SL) and Ae. albopictus Singh's cell line infected with CHIKV (iSL) were made using manufacturer's instructions (Illumina Inc). Briefly, 1 ug of total RNA was ligated with 39and 59adaptors in sequential steps and reverse transcribed using RT primers. Following amplification of the adaptor enriched fragments, population of small RNA within length 140-160 bps were eluted from 6% TBE PAGE gel. Eluted products were precipitated using sodium acetate and ethanol and dissolved in RNAse free water. The samples were quantified using the Total RNA Nano chip in Bioanalyzer. Deep sequencing was performed using Illumina Genome Analyzer II.

Data analysis pipeline
Identification of known microRNAs. Data was analyzed using an in-house pipeline developed for this purpose ( Figure S1). Mature miRNA sequences and pre-miRNA sequences of all insect species available including Aedes aegypti (Aae), Anopheles gambiae (Aga), Acyrthosiphon pisum (Api), Apis mellifera (Ame), Bombyx mori (Bmo), Culex quinquefasciatus (Cqu), Drosophila melanogaster (Dme) and Tribolium castaneum(Tca) were downloaded from miRBase database V.19 [24]. All the sequences were combined for the analysis and 100% similar sequences were pooled together using CD-HIT web tool [25] to reduce data redundancy. Other non-coding RNAs (ncRNA) were downloaded from ncRNA database [26]. Sequences representing coding region of Aedes aegypti, Anopheles stephensi and An. gambiae genomes were also downloaded from Vectorbase [27] using BioMart tool. The tools Bowtie, RNAfold and RNAplot used in this study were downloaded and incorporated in our inhouse PERL based pipeline developed for small RNA analysis. All the downloaded sequences were indexed using Bowtie [28] and used for the analysis. Analysis for all libraries was performed separately. Reads derived from deep sequencing were trimmed and filtered to fetch sequences having length. = 18 bases. Hundred percent matched sequences were pooled together to generate fasta files (expression file) having the unique number along with read counts for each of the unique sequence as an ID, for further analysis. The filtered reads (length. = 18) were aligned against the mature miRNA sequences using Bowtie allowing one mismatches and taking other parameters as default.
Identification of novel miRNA. For identification of novel miRNAs, unmatched sequences after the known miRNA prediction were mapped to pre-miRNA sequences keeping the same parameters as used for known miRNA prediction. For prediction of novel miRNAs, the pre-miRNA unmatched reads were mapped to ncRNA sequences to filter out any other non-coding RNAs from the libraries. The unmatched sequences after this step were then matched to coding region of Ae. aegypti, An. gambie and An. stephensi to filter out the sequences falling in coding regions. The final set of unmatched sequences were then matched to respective genomes and matched sequences were subjected to novel miRNA prediction pipeline which uses RNAfold and RNAplot [29] for secondary structure prediction of pre-miRNA and mapping of small RNA into the pre-miRNA. Those sequences matching to the genome, 75 nt flanking region from both sides of the small RNA sequences were fetched as precursor sequences. These presequences were then folded and energies were calculated by RNAfold and RNAplot. Novel miRNA sequences were then filtered on the basis three main criteria: the hairpin loop structure of the fetched precursor sequence, the presence of the small RNA sequence in any of the arm of precursor sequence, energy , = 220 KJ/mol.
Relative abundance and expression profiling of microRNAs from small RNA sequencing libraries. miRNA abundance in both libraries were calculated after normalizing the reads. The samples were normalized by calculating tag per million of total RNA reads (TPM) and was used for comparing relative abundance of specific microRNAs within each data set. For the purpose of identification of known miRNAs from different

Author Summary
Small, non-coding, RNAs such as microRNAs (miRNA) of length 18-24 nucleotide are reported in wide range of organisms and known to play a determinative role in regulation of various cellular processes. They have also reported to have regulatory roles during pathogen infection, maturation and transmission. Chikungunya virus is an important alphavirus transmitted through two important vectors, Aedes aegypti and Ae. albopictus. While Ae. aegypti has been the preferred vector for transmission of this virus, recently, Ae. albopictus has gained notoriety for disease transmission. Using Illumina platform, we sequenced the small RNA population an Aedes albopictus cell line infected with chikungunya virus. We studied the expression modulation of miRNAs upon CHIKV infection and found distinct set of miRNAs regulated in Aedes upon CHIKV infection. We further predicted the targets of these miRNAs and performed pathway analysis of these targets. The cellular functional targets of some of these miRNAs overlapped suggesting aggressive participation of the targeted pathways in establishing CHIKV infection. organisms, our pipeline was designed in a manner to treat reads with even single nucleotide variation and also reads with different lengths as a distinct small RNA. However, for profiling the relative abundance of miRNAs, we combined all read counts with same miRNA names so as to reduce redundancy in data and to increase accuracy of the reads count of the miRNAs. The miRNAs in iSL (Singh's cell line infected with CHIKV) were further profiled against the SL (Singh's cell line uninfected) to identify the differentially expressed miRNAs upon CHIKV infection.
Expression profiling of miRNA by real time PCR. From the list of miRNAs, two miRNAs namely miR-281-5p and miR-75p were selected at random. Real time PCR analysis was performed using ncode real time PCR kit (cat. MIRCQ 100). Real time PCR was set up in triplicates for miRNAs of uninfected and 24 hours CHIKV infected cells as per the manufacturer's instructions. 5.8 s rRNA was used as an endogenous control for miRNA expression profiling. Expression levels were then calculated against SL as a calibrator using 2 2DDC T method.
Comparison of miRNA expression profiles in different vectors upon pathogen infection. miRNAs identified in the present study were compared with miRNAs identified in a recent study [22] to understand the behavior of these miRNAs in different mosquitoes infected by pathogens specific to them. miRNAs that were common were analyzed to understand mosquito/pathogen specific miRNA regulation. Statistical analysis. Statistical tests for identifying significant differentially expressed miRNAs were performed using edgeR module with few modifications in the script. The pvalue cutoff was performed on the data with the significance threshold selected as 0.05.

Target prediction of microRNAs and pathway
analysis. For identifying putative targets of differentially up regulated miRNAs, 39 UTR sequences of Ae. aegypti and An. gambiae were downloaded from VectorBase and were subjected to target prediction using RNA hybrid tool [30]. The targets were filtered on the basis of complementarity of the miRNA with the targets and energy of the miRNA:target duplex, = 220 Kcal/mol. The targets of selected down-regulated and up-regulated miRNAs were further subjected to KOBAS analysis for the identification of significant pathways, selecting the species as Ae. aegypti and An. gambiae and KEGG as pathway database keeping the threshold p value #0.05 [31]. The vectorbase ID for the targets of significant pathways were then used to generate the miRNA:mRNA interaction network. Cytoscape was used for visualizing the networks [32].

Characterization of Aedes miRNAs in Ae. albopictus cell line
To evaluate the impact of CHIKV infection on Ae. albopictus miRNAs, we carried out high throughput small RNA sequencing using Illumina platform on small RNA libraries obtained from CHIKV infected Ae. albopictus cell line Singh's line and compared it with uninfected cell line. Details of small RNA sequencing data generation of the two libraries and subsequent analyses are detailed in Table 1. In the libraries, there was an overall reduced number of reads in iSL as opposed to SL which could be due to the impact of infection. There was a bimodal pattern of reads, one distinct peak at 22 nt and another peak at 27-30 nt in both libraries. A total of 79.19% reads of SL, 81.94% reads of iSL were mappable and these reads were mapped against Ae. aegypti genome. Analysis pipeline for identification of known microRNAs included zero and one mismatch to accommodate the sequence differences between Ae. albopictus and Ae. aegypti. The remaining raw reads were then mapped against chikungunya genome to discover CHIKV encoded small RNA population. We found a total of 18589 reads (. = 16 nt) mapping against chikungunya genome which was too low. We suspect that the overall lesser number of reads in the infected library to be a reason for this reduced number of virus derived small RNA population.
Mapping of reads from both SL and iSL against known mature miRNA library resulted in identification of 88 and 79 miRNAs of SL and iSL respectively ( Table 2). Among the identified miRNAs, 36 were conserved in three species of mosquitoes (Aedes, Anopheles and Culex) and 12 were conserved in all the eight species we used for our analysis.
Furthermore, a total of nine novel miRNAs were identified in both the libraries.Out of these miRNAs, novel-miR-2 was found to be having miR*.The structures of these novel miRNAs are shown in Figure S2.
Using the read counts, we further calculated the abundance of the miRNAs in the libraries. Using TPM values, the miRNAs were categorized into (i) significantly abundant (SA) with TPM greater than equal to 1000, (ii) moderately Abundant (MA) with TPM 100-999, (iii) abundant (A) with values ranging from 10-99 and (iv) rare (R) with values less than 10. In case of SL, 53 were found to be abundant, while in iSL, 51 were found to be abundant. Several miRNAs in iSL was rare (n = 29) while in SL, the number of miRNAs with counts less than 10 were 36 ( Table 2). miR-184-3p was found to be the most abundant in SL as well as in iSL.

Regulation of Aedes miRNAs upon CHIKV infection
Log fold change of miRNAs used in our study is depicted in Fig. 1. Statistical analysis performed using the edgeR package revealed 41 miRNAs to be significant differentially expressed in Aedes (p value,0.05). It is noteworthy that most of the miRNAs were down-regulated upon CHIKV infection. Based on p-value and log fold change, we identified four miRNAs (miR-100, miR-283, miR-305-3p and miR-927) to be significantly over-expressed and four (miR-1000, miR-2b, miR-2c-3p and miR-190-5p) to be under-expressed.
To further validate the expression of these miRNAs, real-time PCR was performed on SL and iSL with miR-281-5p and miR-7 (Fig. 2). The miRNA expression profiles were found to be comparable to that derived from small RNA sequencing. The amplified products were cloned and sequenced.

Target prediction and network analysis
Further to miRNA profiling, efforts were taken to identify putative targets of all the miRNAs that were regulated in Aedes upon CHIKV infection. Four significantly over-expressed miR-NAs (miR-100, miR-283, miR-305-3p and miR-927) and 4 underexpressed miRNAs (miR-1000,miR-2b,miR-2c-3p and miR-190-5p) were analyzed for their putative targets. For the over-expressed miRNAs, 375 targets were identified for miR-100, 469 for miR-305-5p, 108 for miR-927, 186 for miR-283 and 395 targets for aal-miR-2944b-5p. For down-regulated miRNAs, 460 targets were identified for miR-1000, 226 for miR-2b, 263 for miR-2c-3p and 193 putative targets for miR-190-5p. These targets were further analyzed and clustered using KOBAS and the targets of the significantly represented pathways (p value,0.05) studied (Figure 3a). KOBAS analysis for up-regulated miRNAs revealed that the four up-regulated miRNAs targets the mRNAs playing role in seventeen different pathways. Out of all these pathways, Natural killer cell mediated cytotoxicity and Protein processing in endoplasmic reticulum pathways were common among miR-100, miR-283 and miR-305-3p. Citrate cycle (TCA cycle),dorso-ventral axis formation and valine, leucine and isoleucine degradation were found to be common between miR-100 and miR-305-3p and SNARE interactions in vesicular transport was commonly targeted by miR-927 and miR-305-3p. miR-305-3p was found to be mainly targeting the pathways essential for viral entry such as,ECMreceptor interaction,endocytosis,SNARE interactions in vesicular transport. Likewise in the case of down-regulated miRNAs, ribosome pathway was found to be common among miR-1000, miR-2b and miR-2c.
The interaction of these targets was studied and the miR-NA:mRNA interaction network is shown in Fig. 3b. Several transcripts were targeted by two or more up-regulated and downregulated miRNAs. However, three transcripts, AAEL008528, AAEL013939 and AAEL009042, were found to be common targets among the up-regulated miRNAs and transcript, AAEL003427 was found to be a common target among the down-regulated miRNAs. Detailed investigations are important to understand the role of these mRNAs during CHIKV development in Aedes.

miRNA regulation is specific to vector and pathogen
An ongoing study in our laboratory identified miRNAs from An. stephensi and studied their temporal modulation upon Plasmodium infection. While analyzing the data, we recognized several common miRNAs between the two vectors namely, Ae. albopictus and An. stephensi. We now asked this question whether the modulation of these common miRNAs were species/pathogen specific. For this purpose, we studied the differential expression of these miRNAs. Upon analyzing the common miRNAs between the two vectors, 26 miRNA were found to be significant differentially expressed (p value,0.05). We also observed reper- Table 1. Details of small RNA sequencing information and subsequent data analysis.  Amongst these miRNAs, one miRNA, miR-2944-5p showed upregulation upon both CHIKV and Plasmodium infections. Similarly, miR-2b showed significant down-regulation in both insects. Targets for these two miRNAs were also predicted, subjected to KOBAS analysis and pathway analysis performed. It was seen that target prediction of aal-miR-2944b-5p revealed 395 targets in Aedes and ast-miR-2944b-5p showed 392 targets in Anopheles. Four pathways namely, protein processing in endoplasmic reticulum, citrate cycle (TCA cycle),ribosome and ubiquitin mediated proteolysis were found to be common among them (Fig. 4a). Network analysis of the targets showed 13 transcripts are shared between the two insects (Fig. 4b).
In addition to the above mentioned pathways, targets in Immunity cluster were studied in detail in order to understand role of insect immunity during different infections. It was observed that four target genes belonging to Immunity were common to both the miRNAs and only 2 genes interacting with each other in the network ( Figure 5).
In case of miR-2b, target prediction revealed that aal-miR-2b targeted 226 transcripts, whereas ast-miR-2b had 321 targets. KOBAS analysis of these targets showed two significant pathways for aal-miR-2b and five pathways for ast-miR-2b (Figure 6a), of which, ribosome pathway was found to be common. Interaction network of these targets with miRNA showed three common transcripts (AAEL010168, AAEL006511 and AAEL003427) were targeted by both the miRNA homologs (Figure 6b). No significant targets of immunity pathways were found for miR-2b.

Discussion
In this study we have identified miRNAs of a major vector, Aedes, using an Ae. albopictus cell line and profiled their expression upon CHIKV infection using NextGen sequencing strategies. We identified 90 miRNAs in Ae. albopictus of which eight miRNAs common to both uninfected and infected Singh's line were differentially regulated upon CHIKV infection. Through computational biology, we have further predicted targets of the modulated miRNAs and studied the pathways of these targets. Using information of miRNAs that are modulated in Anopheles upon Plasmodium infection [22], we compared the expression profiles of the common miRNAs in these two distinct vectors upon specific pathogen infection. Even though the data sets derived are from two diverse conditions, i.e., whole insect in case of Anopheles and cell line in the case of Aedes which is an over-simplified derivation of the whole insect, the findings provide insight to understanding the role of miRNAs as seen in previous reports [33]. Our analysis show different miRNAs that are pathogen/ insect specific are up-regulated in the two insects and several of the targets of these miRNAs interact with one another to perform their functions.
Deep sequencing has revealed that distinct repertoire of miRNAs are expressed abundantly in Aedes. It was observed that miR-184 was found to be most abundant in both libraries. Studies have shown that this miRNA multiple roles in female germinal development in Drosophila [34]. We also found miRNA-2940-3p to  be expressed in abundance in Aedes as reported by Skalsky et al. [33]. A recent study has shown the involvement of this miRNA in Wolbachia colonization in Aedes [17]. It is interesting to note that the expression of miR-2940 was down-regulated upon CHIKV infection. Detailed studies are needed to understand the interplay between Wolbachia and other viruses in Aedes and the role this miRNA might play in their growth and development.
Upon CHIKV infection, we showed that distinct miRNAs, namely, miR-100, miR-305-3p, miR-283 and miR-927-5p were up-regulated in Aedes. Likewise, we also showed that miR-1000, miR-2b, miR-2c and miR-190-5p were down-regulated upon CHIKV infection.We used computational biology to predict the targets of these miRNAs. Using KOBAS, we identified that several metabolic and signaling pathways were altered upon CHIKV infection, with protein processing in endoplasmic reticulum being the common pathway in three of the four up-regulated miRNAs analyzed. Viral replication is active in the first 24 hours of infection and most of the host transcription machinery is hijacked by the virus which could be the reason for these pathways being affected upon infection. Another pathway that was significantly altered was part of the immunity pathway showcasing the involvement of insect immunity upon pathogen infection and the role microRNAs may be playing in regulation of insect immune system.
We further proceeded to evaluate the interactions of the predicted transcripts in these pathways with the miRNAs. We observed that several of the transcripts were common among the miRNAs. Of interest were three transcripts that were common among three of the four up-regulated miRNAs. It was seen that these three transcripts encoded for a protein tyrosine phosphatase   Similarly, one transcript, a 40S ribosomal protein S16, was common is three of the four miRNAs down-regulated. Each of these targets is important in different aspects of host response to virus pathogenesis as is evidenced by several reports. It is known that viruses activate extracellular signal regulated kinases(ERKs) during infection [35,36]. A recent study has highlighted the importance of this kinase in antiviral defense in insects during arbovirus infection [37].Similarly, it is known virus entry is regulated by SHP2 and it upstream modulators [38].Similarly, ribosomal proteins (AAE003427, AAEL010168) were common targetsin miRNAs down-regulated during CHIKV infection.
Further, we hypothesized that comparison of different vectors upon different pathogen infections will yield pathogen-vector specific modulators. For this purpose, we compared the miRNA population of Aedes infected with CHIKV and Anopheles infected with Plasmodium. We identified 76 common miRNAs in the two insect datasets, of which one miRNA, miR-2944-5p, was significantly up-regulated and one miRNA, miR-2b, that was significantly down-regulated. We argued that predicting the targets of these miRNA homologs in the two insects will result in distinct targets and provide information as to the host cell regulators to specific pathogens. As seen in the other up-regulated miRNAs, signaling and metabolic pathways were most altered in miR-2944-5p upon Plasmodium or CHIKV infection. In addition, immunity pathway was also studied in detail and revealed interesting information. Pathway analysis revealed four genes, AAEL007494, AAEL010470 (both coding for a calcium binding protein), AAEL008528 (a tyrosine protein phosphatase nonreceptor type 11 (SHP2), AAEL013786 (growth factor receptor binding protein 2). Literature search for these proteins revealed that a growth factor receptor binding protein 2 (GRB 2) binds to NS5A of HCV (Hepatitis C virus) in a homology 3 domain/liganddependent manner and perturbs mitogenic signaling [39]. Another study has implicated GRB2 (Growth factor receptor-bound protein 2) in murine leukemia virus entry through its interaction with ectotropic MLV receptor [40]. A more recent study has shown that Grb2 controls phosphorylation of FGFR2 (Fibroblast growth factor receptor 2) by inhibiting receptor kinase and SHP2 phosphatase activity [38]. It is significant to note that SHP2 phosphatase is another target of miR-2944-5p as well as the three miRNAs upregulated that was discussed earlier in this study.
Our present work has identified known and novel miRNAs in Ae. albopictus. We have analyzed the distribution and abundance of these miRNAs and evaluated their modulation upon CHIKV infection. Target prediction and pathway analysis of these miRNAs were also Figure 4. miRNA-targets KOBAS analysis and interaction network of miR-2944b-5p of Ae. aegypti and An. stephensi. (A) Graph represents significant pathways of aal-miR-2944b-5p and ast-miR-2944b-5p after KOBAS analysis. (B) miRNA-target interaction network of miR-2944b-5p of both the insects showing the shared common targets. doi:10.1371/journal.pntd.0002616.g004 Figure 5. miRNA-targets interaction network of immunity genes for aal-miR-2944b-5p and ast-miR-2944b-5p. Common targets between both the miRNAs are marked as green star. doi:10.1371/journal.pntd.0002616.g005 performed throwing light to the possible role of these miRNAs. Further comparison of some of these Aedes miRNAs with data generated from An. stephensi miRNAs regulate upon Plasmodium infection has revealed interesting information on microRNAs of these two important vectors, upon pathogen infection. While An. stephensi is an important urban vector of South East Asia, Ae.albopictus is gaining notoriety for the spread of chikungunya and more recently dengue. Analysis of the expression of miRNAs and their putative targets in our study have provided important action points for studying the basis of pathogenesis in these vectors.