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. 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. 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.
Vector borne diseases like malaria, dengue, chikungunya are a major burden to the economy of countries in the tropical and sub-tropical 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]–[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]–[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 3′ untranslated regions (3′UTRs) 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.
Ae. albopictus cell line Singh' s line (ATCC-CCL-126) was maintained at 28°C 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. 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 3′and 5′adaptors 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.
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. Mapping of known miRNAs was performed using the genome of Ae. aegypti as the genome of Aedes albopictus is not available. To accommodate the difference in the genome sequences, we extended the analysis to include miRNA reads having upto 2 mismatches with Ae. Aegypti genome. Our analysis with one mismatch revealed additionally29 miRNAs in the libraries, namely, miR-11-5p, miR-12-3p, miR-184-5p, miR-1889-3p, miR-1889-5p, miR-275-5p, miR-277-5p, miR-278-5p, miR-286-3p, miR-2940-5p, miR-2941, miR-2945-5p, miR-296, miR-2951-3p, miR-2a-5p, miR-306-5p, miR-308-5p, miR-309b-3p, miR-316, miR-317-3p, miR-317-5p, miR-34-3p, miR-745-3p, miR-87, miR-92a-5p, miR-92b-5p, miR-980-5p, miR-998-5p and miR-9b. The identification of 29 new miRNA after allowing 1 mismatch also indicates the differences in Aedes aegypti and Aedes albopictus genomes (Table 3). Further to identification of Aedes miRNAs, we analyzed the distribution and abundance of these miRNAs within the libraries. Among 90 miRNAs identified in this study, majority of miRNAs were found in both the libraries (n = 77). Eleven miRNAs were found only in uninfected SL (miR-1175-3p, miR-988-5p, miR-137, miR-219, miR-1, miR-193, miR-210-5p, miR-2796, miR-2944a-3p, miR-932-5p, miR-iab-8-5p), while two miRNAs (miR-286b and miR-1891) were found to be unique to iSL. When we performed the analysis with 1 mismatch, we were able to identify miR-286b reads in SL; however, miR-1891 was absent in uninfected library in 1 mismatch data (Table 3). With respect to novel miRNAs, all nine novel miRNA were found to be present in both SL and iSL libraries. Out of these miRNAs, two miRNAs, novel-miR-4 and novel-miR-9 showed significant up-regulation in iSL and 5 miRNAs, novel-miR1-5p, novel-miR-2-3p, novel-miR-5, novel-miR-7 and novel-miR-8 showed significant down-regulation in iSL library (Table 4). 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. 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. 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 miRNAs (miR-100, miR-283, miR-305-3p and miR-927) and 4 under-expressed 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, ECM-receptor 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 miRNA: mRNA interaction network is shown in Fig. 3b. Several transcripts were targeted by two or more up-regulated and down-regulated 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. 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 repertoires of miRNAs showing distinct expression pattern upon CHIKV infection in Aedes and Plasmodium infection in Anopheles. Amongst these miRNAs, one miRNA, miR-2944-5p showed up-regulation 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.
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 (SHP2), ERK1/2 and an ubiquitin fusion degradation protein. 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 non-receptor 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/ligand-dependent 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 up-regulated 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 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.
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.
biology
Next Generation Sequencing Reveals Regulation of Distinct Aedes microRNAs during Chikungunya Virus Development