Next Article in Journal
microRNA-99a Restricts Replication of Hepatitis C Virus by Targeting mTOR and De Novo Lipogenesis
Next Article in Special Issue
Distinct New York City Aedes albopictus Mosquito Populations Display Differences in Salivary Gland Protein D7 Diversity and Chikungunya Virus Replication
Previous Article in Journal
Multi-Laboratory Comparison of Next-Generation to Sanger-Based Sequencing for HIV-1 Drug Resistance Genotyping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes

Australian Infectious Disease Research Centre, School of Biological Sciences, The University of Queensland, Brisbane QLD 4072, Australia
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(7), 695; https://doi.org/10.3390/v12070695
Submission received: 19 May 2020 / Revised: 19 June 2020 / Accepted: 25 June 2020 / Published: 27 June 2020
(This article belongs to the Special Issue Alphaviruses: Interactions between Arboviruses and Mosquitoes)

Abstract

:
Alphaviruses are globally distributed and predominately transmitted by mosquitoes. Aedes species are common vectors for the clinically important alphaviruses—Chikungunya, Sindbis, and Ross River (RRV) viruses—with Aedes aegypti also being a vector for the flaviviruses dengue, Yellow Fever, and Zika viruses. Ae. aegypti was putatively implicated in the large 1979–1980 South Pacific Islands outbreak of RRV—the leading cause of arboviral disease in Australia today. The RNA interference (RNAi) defense response in mosquitoes involves a number of small RNAs, with their kinetics induced by alphaviruses being poorly understood, particularly at the tissue level. We compared the small RNA profiles between RRV-infected and non-infected Ae. aegypti midgut and fat body tissues at 2, 6, and 12 days post-inoculation (dpi). RRV induced an incremental RNAi response, yielding short interfering and P-element-induced-wimpy-testis (PIWI)-interacting RNAs. Fourteen host microRNAs were differentially expressed due to RRV with the majority in the fat body at 2 dpi. The largely congruent pattern of microRNA regulation with previous reports for alphaviruses and divergence from those for flaviviruses suggests a degree of conservation, whereas patterns of microRNA expression unique to this study provide novel insights into the tissue-specific host-virus attributes of Ae. aegypti responses to this previously unexplored old-world alphavirus.

1. Introduction

Members of the Alphavirus genus (Togaviridae) are mosquito-borne enveloped RNA viruses that are globally distributed, and those of medical significance are attributed to causing encephalitic (e.g., Eastern (EEEV), Western (WEEV), and Venezuelan Equine Encephalitis (VEEV) viruses) or arthritogenic (e.g., Ross River virus (RRV), Chikungunya virus (CHIKV), and the archetypal Sindbis virus (SINV)) diseases [1]. Alphaviruses are maintained through propagation in sylvatic transmission cycles (animal–mosquito–animal) between amplifying vertebrate hosts in which the titer of viremia is high and sustained [2]. Occasionally, transmission cycles encroach into urban environments via short-term amplifying hosts such as peridomestic animals (e.g., bats, birds, cats, dogs, or possums), however, during epidemics, transmission cycles that are entirely urban (human–mosquito–human) can occur in densely populated regions [3].
RRV is an arthritogenic virus endemic to Australia, Papua New Guinea, and South Pacific Islands [4,5,6]. The virus causes epidemic polyarthritis (EPA; otherwise known as Ross River disease or fever) and is the leading cause of viral disease in Australia with a rolling annual mean between 2012 and 2017 exceeding 6000 cases [7]. On the basis of previous estimates and factoring in inflation, the national cost of EPA in 2018 exceeded AUD 9.3 million in medical expenses, whereas the true magnitude of associated costs is unclear when also factoring in mosquito surveillance, control, and response activities (e.g., vector control exceeds AUD 13 million per annum), or outbreak impact on tourism and industry [8,9].
RRV was first isolated from Aedes vigilax near Townsville, Queensland, in 1959 [10]. In Australia, Ae. vigilax is an epidemiologically important saltmarsh mosquito in the coastal regions of the north [11]. Other important species of mosquitoes include Aedes camptorhynchus in the south [4], the inland freshwater species Culex annulirostris [11], and Aedes notoscriptus in urbanized regions [12]. Interestingly, Ae. notoscriptus is associated with transmission cycles of RRV involving the Australian brushtail possum [13,14], both of which are prevalent in New Zealand [15,16]. Although there are no records of RRV circulation in New Zealand, the presence of Ae. notoscriptus and brushtail possums indicates a significant potential for local transmission if RRV were to be introduced [17]. There is no definitive record of the species of mosquitoes involved in the Pacific Islands epidemic of 1979–1980, whereas those implicated on the basis of prevalence and capacity to readily transmit RRV included Ae. vigilax, Cx. annulirostris, Aedes polynesiensis, Aedes oceanicus, Culex quinquefasciatus, and Aedes aegypti [18,19]. RRV has been isolated from 42 species of mosquitoes [4,11,20], whereas laboratory competence studies have only indicted 10 as being capable of transmission [4]. RRV has never been isolated from Ae. aegypti in the field, and this species of mosquito is not considered a vector of RRV; however, multiple competency studies have shown that Ae. aegypti can readily transmit the virus (e.g. [21]).
Ae. aegypti is unrivalled in medical importance, being culpable for our most impactful diseases including CHIKV, dengue virus (DENV), Yellow Fever virus (YFV), and Zika virus (ZIKV), wherein, collectively, tens of millions of cases of disease and tens of thousands of deaths occur globally each year [22,23,24,25]. The medical importance of Ae. aegypti is underscored by its capacity to adapt to a range of ecological niches, both spatial and temporal, to transmit multiple disparately related and highly pathogenic viruses, and for its natural drive to preferentially seek out human blood. The capacity of Ae. aegypti to vector a wide range of viruses and its ability to adapt to foreign environments stresses our need to better understand the molecular host–virus interactions in this mosquito.
Genomic studies analyzing differentially regulated genes have enabled significant progress in elucidating mosquito innate immune responses to arbovirus infection. Sim et al. (2005) pioneered this approach, reporting the modulation of Anopheles gambiae gene expression in response to O’nyong-nyong virus (ONNV: Alphavirus) infection, and later in 2007 they demonstrated that this modulation resulted in the production of heat shock protein cognate 70B (HSC70B), which impeded ONNV replication [26,27]. In response to arboviral infection, Ae. aegypti has been reported to initiate Toll [28], c-Jun N-terminal kinase (JNK) [29], immune deficiency (Imd) [30], and Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathways [31]. The JAK/STAT pathway has been reported to be utilized by Ae. aegypti to inhibit DENV [32] but not CHIKV [31,32], whereas there are contradictory reports for ZIKV [32,33].
RNA interference (RNAi) has been recognized as an important mechanism by which eukaryotes regulate transcriptional and post-transcriptional gene expression, mRNA methylation, and chromatin remodeling [34]. The RNAi pathway produces small non-coding RNAs (sncRNAs: ≈20–30 nt) that are involved in a range of biological processes including host–pathogen interactions and anti-viral responses [34,35]. In insects, the major sncRNAs are short interfering RNA (siRNA), P-element-induced-wimpy-testis (PIWI)-interacting RNA (piRNA), and microRNA (miRNA).
Campbell et al. (2008) showed that Ae. aegypti produced complementary siRNAs in response to a SINV challenge and that silencing Argonaut 2 (Ago2) or Dicer-2 (Dcr2) endonucleases caused a transient increase in SINV replication [36]. Concomitantly, Cirimotich et al. (2009) demonstrated that Flock House virus (FHV: Nodaviridae, genus Alphanodavirus) B2 protein inhibited the siRNA response by shielding SINV double stranded RNA (dsRNA) from RNAi detection, resulting in a fatal viraemia in Ae. aegypti [37]. Additionally, Ae. aegypti has been shown to use siRNA to control DENV [38] and CHIKV [31], and in Ae. aegypti-derived Aag2 cells, the Semliki Forest virus (SFV: Togaviridae, genus Alphavirus) [39]—examples that underscore the versatility of the RNAi response.
Several studies have detected alphavirus-derived piRNAs in mosquitoes infected with SFV, CHIKV, ONNV, and SINV, although the consequence of these piRNA responses remain unclear [40]. In Ae. aegypti, only 19% of total piRNAs are directed against transposons [41], whereas a substantial portion of genes involved in germline and embryonic development in An. gambiae encode piRNA transcripts, suggesting that piRNAs in mosquitoes have multiple functions [34].
The ingestion of viremic blood by a mosquito results in a diverse and complex miRNA response [42]. Liu et al. (2015) showed that miR-1174 and miR-1175 were upregulated in the gut of Ae. aegypti following a blood meal [43]. Bryant et al. (2010) showed that miR-8, miR-14, and miR-275 were upregulated in the fat body of Ae. aegypti following a blood meal and that these are involved in digestion and egg development [44]. Campbell et al. (2014) reported 35 differentially abundant miRNAs in Ae. aegypti in response to DENV2 infection, with the majority of predicted targets involved in transport, signal transduction, and metabolism [45]. In the saliva of Ae. aegypti infected with CHIKV, Maharaj et al. (2015) identified five miRNAs that were predicted to be involved in the regulation of CHIKV infection [46].
There are very few studies demonstrating miRNA responses in mosquitoes infected with RNA viruses at the tissue level. In particular, there are no studies on RRV. This may be because none of the important vectors of RRV have been sequenced and therefore a platform for comparative genomics to an epidemiologically relevant reference cannot be made. The genome of Culex quinquefasciatus has been sequenced [47], although vector competence studies have shown this mosquito to be poorly susceptible to RRV infection and, therefore highly unlikely to be involved in RRV maintenance [48]. Ae. aegypti has been shown to competently transmit and potentially facilitate the global spread of a range of viruses, including RRV, which stresses our need to better understand the repertoire of molecular host–virus interactions in this mosquito. Further, the complete genome, transcriptome, and miRNAome of Ae. aegypti are well characterized. In this study, we revealed the differentially abundant small RNA profile of the pertinent organs (midgut and fat body) in Ae. aegypti infected with RRV. These results provide an important and relevant contribution to our understanding of the molecular host–virus interactions within a medically significant vector infected with an old-world alphavirus.

2. Materials and Methods

2.1. Ethics Statement

The T48 strain of RRV used in this study was kindly provided by Jody Peters of the School of Chemistry and Molecular Biosciences, The University of Queensland. The virus was originally isolated from Ae. vigilax mosquitoes (Townsville, QLD, Australia). RRV protocols were approved by the University of Queensland Biological Sciences Biosafety Committee (Reference number: IBC/231B/SBS approved on 22 May 2018).

2.2. Mosquito Infections with Ross River Virus

Ae. aegypti (Innisfail strain) mosquitoes were propagated from eggs and raised in an insectary maintained at 27 °C, 80% relative humidity, and a 12 h photoperiod. RRV (T48 strain) was amplified in C6/36 cells (Ae. albopictus) cultured at 28 °C in a 1:1 v/v mixture of Schneider’s Drosophila Medium (Life Technologies, California, USA) and Mitsuhashi and Maramorosch (Himedia, Mumbai, India) medium supplemented with 2% fetal bovine serum (FBS). The virus was then amplified in African green monkey kidney (Vero) cells cultured at 37 °C in Opti-MEM (Life Technologies, California, USA) media with 2% FBS for 48 h. Control C6/36 and Vero cells that did not include virus were cultured each time and in the same way. The same volume of C6/36 cell supernatant that was used to inoculate Vero cells was also transferred between control cultures. Virus amplified in Vero cells was mixed 1:1 v/v with whole human blood (Australian Red Cross) and 0.1% ATP (phagostimulant). Females aged 4–6 days old were fed blood mixtures with or without virus at 2 × 105 plaque forming units per mL (pfu/mL). At 2, 6, and 12 days post-feeding (dpf), mosquitoes were cold anaesthetized and transferred to −80 °C storage until use.

2.3. Mosquito Dissection and RNA Extraction

The midgut and fat body tissues were dissected from a total of 180 mosquitoes yielding 360 organs. The organs were pooled in lots of 10 yielding 36 samples comprising the tissues fat body (N = 18) or midgut (N = 18) that were infected (N = 9) or uninfected (N = 9) at day 2 (N = 3), day 6 (N = 3), or day 12 (N = 3). Total RNA was extracted from samples (N = 36) using the miRNeasy Mini Kit following the manufacturer’s protocol (Qiagen, Hilden, Germany). For practical reasons, we did not confirm infection or non-infection in all mosquitoes (N = 180). During dissection, the RNA from the carcass remnants of 10 randomly selected mosquitoes sampled from each day and treatment was used to detect virus presence or absence (control) by reverse transcription quantitative PCR (RT-qPCR) on a Rotor-Gene Q (Qiagen, Hilden, Germany) machine using ribosomal protein subunit 17 (RPS17) and in-house-designed RRV-specific primers (Table S1). Mean normalized expression (MNE) of samples relative to RPS17 was calculated using the second derivative maximum take-off and amplification efficiency values produced by the Rotor-Gene Q software.

2.4. Library Preparations and Sequencing

The TruSeq Small RNA protocol (Illumina) was used to generate and sequence complementary DNA (cDNA) libraries from small-input RNA on a short-read sequencing platform with an Illumina Hi-Seq machine (Illumina). Briefly, cDNA libraries were amplified from 1 µg of size-fractionated total RNA ligated with indexed 3’ and 5’ adapters (Table S1). Libraries were size-selected by polyacrylamide gel, purified, and concentrated. An Agilent Technologies 2100 Bioanalyzer (California, USA) was used to assess purity for quality control. Illumina RTA software (1.17.21.3) was used for base-calling for quality control, and Illumina bcl2fastq 2.17 software was used for demultiplexing. Final raw reads were in FASTQ format with the 3’ adapter retained. The raw data for infected and uninfected Ae. aegypti samples generated in this study were deposited in the Sequence Read Archive (SRA) under the accession PRJNA635740.

2.5. Small RNA Analysis

The software tool Fastq_quality_filter (FASTX-toolkit) (v0.0.6) [49] was used to earmark raw reads that had at least 95% of their sequence content with a phred quality score above 20. Using Flexbar [50], we trimmed the 3’ adaptors (Table S1), allowing a minimum overlap of 4 nucleotides (nt), and remaining reads with less than 16 nt were discarded. The length distribution of reads was visualized using scales [51], plyr [52], and ggplot2 [53] in the R package (v3.6). The miRPro pipeline [54] was used to enumerate and categorize the samples as miRNA content and other RNA biotypes. Cleaned reads were mapped to the most recent miRBase repository (v21) [55] using Novoalign (Novocraft, Selangor, Malaysia) and both mature and novel mature miRNAs were enumerated with htseq-count (HTSEQ 1.11) [56] and Samtools algorithm [57]. Remaining reads were mapped to the most current Ae. aegypti genome (AaegL5; VectorBase) [58] for RNA reads categorization. Differential expression of miRNAs (Table S2) was calculated from concatenated mature and novel mature miRNA counts using DESeq2 (v3.11) [59] following various namesake Bioconductor vignettes in the R package [60,61]. The counts from all samples were included in the experimental design for differential expression, which included the three main factors (1) tissue (tsu), (2) time (day), and (3) condition (cnd) and interaction terms between these factors (design = tsu + tsu:day + tsu:day:cnd). The levels for each factor were (1) tsu: fat body (fb) and midgut (mg); (2) day: 2, 6, and 12; and (3) cnd: control (con) and infected (rrv). Fat body, day 2, and control was set as the reference levels by which change was compared. The coefficients derived from the design enabled comparison between tissues mg and fb irrespective of infection (tsu_mg_vs_fb), the effect of age between days with reference to deviance from day 2 (reference) irrespective of infection (e.g., tsu.fb.day6 or tsu.mg.day12), the effect size of condition within tissues within day (e.g., tsufb.day2.cndrrv), within tissue between days (e.g., contrast = tsufb.day2.cndrrv and tsufb.day6.cndrrv), or between tissues within day (e.g., contrast = tsufb.day2.cndrrv and tsumg.day2.cndrrv). The effect size of condition between tissues between days (e.g., tsufb.day2.cndrrv and tsumg.day12.cndrrv) was not explored. To test for significance, we used the Wald test (alpha, 0.05) for the pairwise statistical comparison of probability calculated between the coefficients generated by our experimental design. The threshold for significance was set at a false discovery rate (FDR) of 0.05. The FDR implemented in DESeq2 is an adjusted p-value given by the Benjamini and Hochberg algorithm [62]. We used the approximate posterior estimation for generalized linear model (apeglm) function [63] to exclude results from low or highly variable counts and further excluded results with a baseMean < 100. In DESeq2, the term baseMean refers to the mean normalized count across all samples. Remaining miRNAs were considered significantly differentially expressed (SDE). In summary, our criteria qualifying miRNA significance were for a normalized mean count of samples higher than 100, a FDR less than 0.05, and an effect size fold change (FC) of 1.5.

2.6. RNAi Activity Analysis

For siRNA and piRNA identification, reads that did not align with either the Ae. aegypti genome (AaegL5) or miRBase records were mapped to the genome of the T48 strain of RRV (GenBank: GQ433359.1) [64] using Novoalign (Novocraft, Selangor, Malaysia). BAM alignment files were filtered by strand and lengths of 21 nt or 27–29 nt for siRNA and piRNA reads analyses, respectively. We used kpPlotBAMCoverage in karyoploteR (v1.14.0) within the R package (v3.6) [65] to visualize the genome coverage of siRNA and piRNA reads. We used WebLogo software (v3) [66] to visualize U1 and A10 piRNA signatures, and 10 nt corroborative overlap signatures were assessed and scored using the public Galaxy server at https://mississippi.sorbonne-universite.fr [67].

2.7. The miRNA Target Identification

The algorithms miRanda [68], RNAhybrid [69], TargetScan [70], PITA [71], and RNA22 [72] were used to predict miRNA targets against Ae. aegypti 3’UTR sequences. For this, a bed file of headers containing the phrase three_prime_utr and corresponding coordinates was extracted from an annotation file (Aedes-aegypti-LVP_AGWG_BASEFEATURES_AaegL5.2.gtf) with grep and used as a template for BEDTools [73] to extract those corresponding 3’UTR sequences from file (Aedes-aegypti-LVP_AGWG_CHROMOSOMES_AaegL5.fa). Sequences shorter than 30 nt were discarded, and to remove duplicates for genes with splice variants, we used only the longest 3’UTR. All target prediction programs were operated locally on a Linux machine using default command line parameters. We referred to the results exclusion criteria used by Oliveira et al. (2017) to guide the exclusion criteria that we used in this study. For miRanda, we used a minimum score of 150 and a minimum folding energy (mfe) of < −18 (kcal/mol) [74]. For RNA22, we used a p-value of < 0.1 and mfe values < −12 (kcal/mol). For RNAhybrid, we set the p-value at 0.05 and mfe values < −20 (kcal/mol). For PITA, we used mfe values < −18 and excluded 6mer types from both PITA and TargetScan predictions. The intersect of targets predicted by at least three programs were used for downstream analysis (Table S3).

2.8. Gene Ontology

The Module Functional Analysis in OmicsBox (Biobam, Valencia, Spain) was used for gene ontology (GO) annotation and analysis (Table S4) of the miRNA targets predicted by at least three programs [75]. This module used National Centre for Biotechnology Information (NCBI) Blast against homologous Ae. aegypti sequences to establish the protein identity of the query sequences from this study. The corresponding GO IDs and GO names of the mapped and annotated proteins were retrieved using OmicsBox InterPro scan and these were linked with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Table S5). The statistics and graphical representation of all the results was produced by the OmicsBox program.

2.9. RT-qPCR Analysis

The miScript II RT kit (Qiagen, Hilden, Germany) was used following the manufacturer’s protocol to synthesize cDNA from 500 ng RNA per 20 µL reaction using the 5× HiSpec Buffer and the proprietary-sequence universal primer on a Veriti thermal cycler (Applied Biosystems, California, USA). cDNA samples were diluted in 200 µL nuclease-free water and 2 µL thereafter amplified in 10 µL duplicate reactions with miRNA-specific primers (Integrated DNA Technologies, Iowa, USA) by RT-qPCR using miScript SYBR Green PCR kit reagents (Qiagen, Hilden, Germany) on a 72-well Rotor-Gene Q thermocycler (Qiagen). The thermocycling parameters used for amplification were 95 °C for 15 min, followed by 40 cycles of 94 °C for 15 s, 60 °C for 30 s, and 70 °C for 30 s. The mean normalized expression (MNE) of samples was calculated using the second derivative maximum and corresponding amplification efficiencies generated by the Rotor-Gene Q software [76]. The Ae. aegypti small nuclear RNA U6 was used as the calibrator. The log2 ratio between control and treatment MNE values were used to calculate fold change. Three biological replicates were used per treatment.

3. Results and Discussion

3.1. Illumina Sequencing of Small RNAs

An Illumina sequencing platform was used to produce small RNA profiles of RRV-infected and non-infected Ae. aegypti mosquitoes. To investigate host miRNA and RNAi responses to RRV infection, we extracted RNA samples from midgut and fat body tissues of mosquitoes collected at 2, 6, and 12 days post-infection (dpi). The rates of viral infection were determined from a subset of individual mosquitoes (N = 10) from days 2, 6, and 12 and the relative viral load in those mosquitoes with a detectable degree of RRV infection (N = 6) was confirmed by RT-qPCR, which demonstrated a progressive increase in viral load over time (Figure 1). No virus was detected in uninfected negative control mosquitoes. Considering that carcasses were used for determination of virus infection after dissection of tissues, 100% dissemination of virus at 12 dpi suggested high rate of infection even if it was detected in about 50% of mosquito carcasses at 2 dpi when virus was at early stages of infection in the midgut, which was already removed from mosquito bodies.
The combined yield of raw data obtained from control and infected small RNA libraries was 200 and 232 million reads, respectively (Table S2). Of the reads, we discarded 6–28% in different libraries due to their low-quality score or lack of adapter sequence. A total of 94% of the annotated Ae. aegypti miRNAs referenced by miRBase were found among all the samples in our data, representing 7–22% of clean reads in different libraries. The sequence length distributions in all libraries expressed peaks at 21–23 nt, consistent with the characteristic lengths of miRNA and short interfering RNA (siRNA), and peaks at 27–29 nt representing piRNA—a common feature of insect small RNA libraries (Figure 2).

3.2. Differential Expression of Ae. aegypti miRNAs in Response to RRV Infection

Small RNA library analysis of RRV-infected Ae. aegypti fat body and midgut tissues identified 14 SDE miRNAs at different time points compared with uninfected mosquitoes (Figure 3). Interestingly, miRNA modulation was prominent early in infection at 2 dpi, and by day 12, when the viral load was higher, most miRNAs in both tissues were stably expressed (Figure 4). In contrast to siRNAs, there is very little evidence that miRNAs directly target viral genomes [77], and therefore we would not expect to see the positive correlation of changes in miRNA expression with viral load. It is then perhaps not surprising to see the dramatic shift in miRNA expression during the early adjustment phase to infection when the greatest number of changes in target gene expression occur. The identified miRNAs were predominately downregulated (60%) and from the fat body (93%) tissue (Table 1). An intriguing observation is the preponderance for the downregulation of Ae. aegypti miRNAs reported for both flavivirus (DENV and ZIKV) [45,78] and alphavirus (CHIKV and SINV) [46,79,80] infections, although direct comparison between these studies and our results is difficult considering the different time points or samples used (whole mosquito/tissues/cell line). Myles et al. (2008) reported an overall depletion of miRNA expression in Ae. aegypti infected with SINV, which was the first record in the literature investigating the miRNA expression of Ae. aegypti infected with an old-world alphavirus, when Ae. aegypti miRNAs had not yet been annotated and the authors’ analysis for differential expression used Drosophila melanogaster miRNAs as reference [80].
miR-9b-5p was upregulated in the midgut (FC 1.67) at 2 dpi and was the only SDE miRNA found in midgut samples (Figure 5A). This miRNA was also upregulated in the fat body at 2 dpi (FC 2.29), which was the most significant result from this study (FDR, 9.55e-08). Maharaj et al. (2015) reported downregulation of this miRNA from a read count of 64 to 34 in the saliva of Ae. aegypti infected with CHIKV at 10 dpi, although the change in miR-9b expression was not considered statistically significant [46]. By contrast, in Ae. albopictus, Liu et al. (2015) found miR-9b to be significantly downregulated at 7 dpi with DENV. Considering about 50% virus dissemination was detectable in mosquitoes at 2 dpi (Figure 1A), we speculate that the early changes in fat body miRNAs could be part of an early systemic response to virus infection following the receipt of cytokines from infected midgut cells via hemolymph. Crosstalk between tissues for immune induction following blood feeding has been documented (reviewed in [81]). In contrast to the significant downregulation of miR-2940-3p in the study by Liu et al. (2015), in our study, miR-2940-3p was upregulated. By comparison, miR-184-3p and miR-2940-5p were significantly downregulated in Liu et al. (2015) [82] and these miRNAs were also downregulated in our study. The continued enquiry of miRNA profile changes elicited by disparate pathogens in this mosquito or other dipterans is required to better clarify their evolutionary conservation and functional significance in response to infection.
The most abundant SDE miRNAs found in this study—miR-184-3p, miR-275-5p, miR-2940-5p, miR-317-3p, and miR-989-3p—are commonly reported as highly abundant and widely conserved [83]. Additionally, miR-989-3p (FC 5.85) and miR-275-5p (FC 0.34), which were both from fat body tissue at 2 dpi, had the most extreme fold changes. By comparison, in An. gambiae infected with Plasmodium bergei, miR-989 was also significantly upregulated (FC 4) and deemed exclusive to midgut tissue [84], whereas other studies had contended that, with reference to results comparing the conserved homologues in both Anopheles stephensi and Ae. aegypti, miR-989-3p is predominantly expressed in ovaries [85], and rarely in the midgut tissue [86]. We detected miR-989-3p in the midgut tissue at 2 dpi, whereas its change in expression compared to uninfected samples was not statistically significant. Other studies reporting the significant differential expression of miR-989-3p involved the vectors Ae. albopictus and Culex quinquefasciatus infected with West Nile virus (WNV) [87], Ae. aegypti infected with DENV [88] or ZIKV [78], or Ae. aegypti cells (Aag2) infected with Wolbachia [89]. According to our criteria for significance, there were no SDE miRNAs detected from mosquitoes at 6 dpi in either tissue.
The only SDE miRNAs from day 12 mosquitoes were miR-71-3p and miR-2940-3p in the fat body tissue. An analysis of the change in expression of miR-71-3p revealed that in control mosquitoes this miRNA became progressively more abundant throughout the time course, whereas in the fat body at 12 dpi, miR-71-3p was significantly downregulated. By comparing samples across time, it appeared that the condition effect due to RRV was an impediment to the natural age-related upregulation of miR-71-3p by day 12 (Figure 5B). It would be interesting to investigate the target genes of miR-71-3p, the consequence of their not being silenced at that age, and whether their continued expression was part of a deliberate anti-viral defense response or caused by the virus.
Saldaña et al. (2017) also noted a disproportionately low number of SDE miRNAs in Ae. aegypti mosquitoes infected with ZIKV, despite the pronounced viral load reached by day 14 [78]. Comparing results for miR-2940-3p from this study and Saldaña et al. (2017), the miRNA shows remarkably similar measures of abundance and degree of upregulation (FC ≈ 2). While whole Ae. aegypti mosquitoes were infected with a flavivirus and expression had peaked earlier at day 7 in Saldaña et al. (2017), the measures of high abundance are consistent with published results for this miRNA from both flaviviruses [78,82] and alphaviruses [79]. Wolbachia is an endosymbiotic bacterium that is currently used to control DENV spread as its presence in Ae. aegypti blocks replication of a number of viruses [90]. In Ae. aegypti, Wolbachia induces and utilizes miR-2940 to suppress the expression of pelo protein [91], upregulate host metalloprotease f41 ftsh transcripts [89], and downregulate the expression of Ae. aegypti DNA methyltransferase 2 [92], activities that aid its preservation in the host and contribute to inhibition of replication and subsequent transmission of co-infected DENV.
Mosquito miRNA induction by CHIKV has been studied extensively, however, research into other alphaviruses is needed to supplement the literature for comparing between species to further develop a characterization of their distinct molecular responses and elucidate with greater clarity the functional reasons for miRNA involvement and why conservation or divergence in miRNA form and function has evolved. Because this project is a foray into the Ae. aegypti miRNA responses to an old-world alphavirus, we prioritized our analysis of miRNA function to those also published with reference to CHIKV or Ae. aegypti fat body.

3.3. RT-qPCR Validation of Differentially Expressed miRNAs

The relative expression of 10 miRNAs was ascertained by RT-qPCR from RNA samples extracted from fat body and midgut tissues from control and RRV-infected samples at 2 dpi to validate the RNA-Seq data and subsequent assessment of differential expression of miRNAs (Figure 6). The 10 miRNA candidates for this experiment were selected to represent the scope of fold changes determined by DESeq2, which included those showing the most (miR-989-3p) and least (miR-999-3p) divergence. The direction of change was mostly consistent between methods with discrepancies for miR-1-3p and miR-2940-3p from midgut and miR-999-3p, bantam-5p, and miR-989-3p from fat body, which notably also showed the largest discrepancy between the reported magnitudes of fold change.

3.4. The miRNA Target Analysis

The putative Ae. aegypti 3’UTR mRNA targets of the SDE miRNAs determined in this study (Table 1) were predicted using the default parameters of miRanda, RNAhybrid, TargetScan, RNA22, and PITA software. From the intersect of targets predicted by at least four programs, we found 788 genes as putative binding sites for at least one miRNA. Among those, 65 genes were associated with at least 2 miRNAs, 6 genes with 3 miRNAs, and AAEL003953 (LOC5579800, ec:2.7.1.123) was found to be a putative target for miR-9b-5p, miR-275-5p, miR-2940-3p, and miR-2940-5p (Figure 7). The greatest number of putative binding sites were for miR-989-3p (367), miR-9b-5p (287), and miR-2940-3p (73), which were also determined in this study to have the highest upregulation and by the same rank (Table 1).
An analysis of the KEGG pathways revealed a number of binding sites that were for genes putatively involved in mosquito immunity and signaling (Table S5). There were 328 sequences with the molecular function annotation (GO:0003674) comprising genes ascribed to hydrolase (78), polymerase (82), and iron binding (57) activities (Figure 8 and Figure 9). Both miR-989-3p and miR-9b-5p were associated with genes involved in zinc-finger protein and the mitogen-activated protein kinase (MAPK) signaling pathway as described by the response to stress ontology (GO:0006350). It has been shown that the expression of zinc-finger antiviral protein inhibits the replication of members of the Alphavirus genus including SINV, SFV, VEEV, and RRV [93]. AAEL013596 is a phosphatidylinositol 3-kinase (PI3K) regulatory subunit and is involved in the Toll pathway. In Drosophila, the Toll-7 transmembrane receptor activates antiviral autophagy facilitated by the phosphatidylinositol-3-kinase and the mammalian target of rapamycin (PI3K-Akt-mTOR) signalling pathway [94], whereas for SINV, which forms membrane bound replication complexes, the activation of this pathway in Drosophila becomes pro-viral, facilitating its cap-dependent viral RNA translation [95]. In human cells, the metabolic changes induced by the hyperactivation of PI3K-Akt by SFV and RRV is also pro-viral [96].
The beta subunit kinase, AAEL003245, is an inhibitor of nuclear factor kappa B and involved in the mTOR signaling pathway, which relates to vitellogenesis, protein synthesis, and motility. Genes were associated with the Wingless signaling pathway (Wnt), which works synergistically with the mTOR pathway in regulating vitellogenesis [97]. Following a blood meal, miR-8 is significantly upregulated in Ae. aegypti fat body and was shown to target the Wnt pathway to promote the production of yolk protein precursors (YPP) and lipid accumulation used for vitellogenesis [98]. Conversely, in this study, miR-8 was significantly downregulated in the fat body of RRV-infected Ae. aegypti mosquitoes, which may signify an active response to stall vitellogenesis. Schwenke et al. (2016) showed that pleiotropic signaling mechanisms regulate the allocation of metabolic resources between reproductive and immune processes in female insects, and that egg production is stopped to reallocate energy resources to the defense response during infection [99]. In humans, elements of the mTOR signaling pathway are modulated by the Influenza A virus [100] and ZIKV [101], which promote viral replication during late stage cell stress.
To corroborate our findings with published evidence, we compared our list of predicted miRNA target genes with studies involving changes in the Ae. aegypti transcriptome due to arbovirus infection. We only included transcripts that were reported as significantly modulated and predicted as a target for at least one miRNA by no less than three algorithms. In total, we found 18 unique genes, representing the predicted targets for 12 out of the 14 miRNAs from our study, which were shown to be significantly modulated in different studies involving whole Ae. aeygpti, midgut and carcass with CHIKV, DENV, WNV, or ZIKV at 1 dpi, 2 dpi, 4 dpi, 7 dpi, or 14 dpi [102,103,104,105]. We focused on significantly downregulated transcripts that were predicted as targets for the significantly upregulated miRNAs in this study. Although no genes that were predicted by all five algorithms were found, our most significantly upregulated miRNAs, miR-9b-5p and miR-989-5p, which had most of their targets predicted by four algorithms, were strongly represented (Table S6). Among the upregulated miRNAs, 80% were paired with downregulated transcripts but only 22% of the downregulated miRNAs were paired with upregulated transcripts. The top 10 most downregulated genes included studies of DENV, WNV, and ZIKV and were from abundances measured early during the course of infection predominately from 1 and 2 dpi. The greatest number of pairings was from transcript abundances measured early during the course of infection at 1, 2, and 4 dpi.

3.5. RRV is a Target of the Ae. aegypti RNAi Response

The dsRNA intermediates produced during viral replication become a target of the mosquito RNAi defense response, as reported for a number of alphaviruses including SINV [36,80,106], SFV [39], and CHIKV [31]. The ribonuclease Dicer-2 cleaves these intermediates into 21 nt double-stranded short interfering RNAs (siRNAs), which when bound to the RNA-induced silencing complex (RISC) become a template of the viral RNA sequence and guide for subsequent cleavage. To investigate potential RNAi activity against RRV, we mapped reads to the virus genome that had not aligned with either the Ae. aegypti genome or miRBase records. There were 30,652, 238,857, and 420,187 reads in total from days 2, 6, and 12, respectively, which mapped to the viral genome with a length distribution of 15–51 nt. The ratio between reads from fat body and midgut were 1:1.67, 1:0.56, and 1:1.26 for days 2, 6, and 12, respectively. The production of 21 nt sequences, a hallmark of an siRNA response to infection, was evident from day 2, and these increased throughout the time course for both tissues (Figure 10). The distribution of reads mapped somewhat evenly across the genome with a greater number in each case mapping to the sense strand (Figure 10). These results confirm an active RNAi response of Ae. aegypti to RRV infection.

3.6. Production of RRV-Derived vpiRNAs

In Ae. aegypti, there are eight members of the PIWI protein clade (Piwi 1–7 and Ago3) that are involved in piRNA production from transposons, mRNA or viral RNA. The Tudor protein Veneno recruits Yb, Piwi5, and Ago3 [107] to mediate virus-derived piRNA (vpiRNA) production [108] by ping-pong amplification [109], where sense and anti-sense vpiRNA pairs show a nucleotide bias for uridine at position 1 (U1) and arginine at position 10 (A10), respectively [110,111]. Alphavirus-derived piRNAs have been detected in mosquitoes infected with SFV, CHIKV, ONNV, and SINV [40].
To investigate if vpiRNAs are produced in Ae. aegypti infected with RRV, we mapped 27–29 nt reads from 2, 6, and 12 dpi to the RRV genome. The number of reads that mapped to the viral genome increased throughout the time course as the infection progressed with a bias for reads mapping to the sense strand (Figure 11). At 6 and 12 dpi for both fat body and midgut samples, there was a moderate clustering of reads at the 3’ end of the genome, which encodes the subgenomic RNA for the production of non-structural proteins (Figure 11). The subgenomic promoters of the alphaviruses SINV, CHIKV, and SFV have been shown to drive subgenomic RNA production in mosquitoes, and in those studies, reads mapped in distinct hotspots to the 3’ non-structural end of the viral genome spanning the subgenomic RNA, and almost exclusively to the sense strand [112]. In our study, clusters of reads were also formed in distinct hotspots spanning the 3’ end of the genome and predominantly mapping to the sense strand. Therefore, it appears that the subgenomic RNA of RRV drives vpiRNA production in Ae. aegypti.
To investigate if the vpiRNAs in this study showed the U1/A10 nucleotide bias characteristic of ping-pong amplification, we used Weblogo software (v3) to determine the nucleotide frequency at each position for reads grouped by the same length, strand, and time-point post-infection. There was a modest A10 bias in fat body anti-sense sequences of 28 nt at 6 dpi (Figure 12A), whereas the corresponding U1 bias was more readily apparent (Figure 12B).
We used the Small RNA Signatures tool from the Galaxy instance Mississippi.sorbonne-universite.fr to investigate 10 nt overlap probabilities in sequences in the range of 27–29 nt that mapped to the RRV genome. The software pairs sequences with their reverse complement and tally the number of pairs of sequences by the number of their overlapping nucleotides. In fat body tissue at 6 dpi, there was a high frequency of sequence pairs with a 10 nt overlap (Figure 13). The prominent 10 nt peaks for both the number of pairs counted and the probability of seeing pairs with a 10 nt overlap are signatures of piRNA production by ping-pong amplification.

4. Conclusions

In summary, we found that RRV promoted the significant differential expression of 14 miRNAs in Ae. aegypti, with most occurring in the fat body tissue at 2 dpi. Several genes related to immunity were predicted as targets for the SDE miRNAs from this study, suggesting that RRV may elicit an active defense response by Ae. aegypti. A prominent increase in 21 nt sequences was apparent at 6 and 12 days post-infection, indicating a strong siRNA defense response in Ae. aegypti infected with RRV. There was a detectable vpiRNA response to RRV in Ae. aegypti fat body, and those sequences showed a distinct U1 bias, although the A10 bias in the corresponding sequences was not as distinct. The overlap probabilities of the 27–29 nt sequences from fat body piRNA-like reads at 6 dpi did show a z-score spike at 10 nt, which provided a complement to the evidence from this study for RRV-derived vpiRNA production in Ae. aegypti. This study provides insight into the small RNA responses of Ae. aegypti to RRV, and in particular we provide the first characterization of the miRNA responses in Ae. aegypti to a previously unexplored old-world alphavirus in this context.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/7/695/s1, Table S1: Primers for RT-qPCR, miScript, and the NGS 3’adapter sequence; Table S2: DESeq2 results of miRNA differential expression; Table S3: miRNA Target prediction; Table S4: OmicsBox Gene Ontology analysis; Table S5: Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways; Table S6: Comparison of transcriptomes.

Author Contributions

Experiments, J.B.S.; conceptualization, S.A.; formal analysis, J.B.S.; data curation, J.B.S.; writing—original draft preparation, J.B.S.; writing—review and editing, S.A.; visualization, J.B.S.; supervision, S.A.; project administration, S.A.; funding acquisition, S.A. All authors have read and agreed to the published version of the manuscript.

Funding

Sassan Asgari is supported by the Australian Research Council, grant number DP190102048. James Sinclair is supported by a University of Queensland scholarship.

Acknowledgments

We wish to thank Rhys Parry from the University of Queensland for helpful discussions in regard to data analysis, and Jody Peters from the University of Queensland for providing the virus.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Contigiani, M.S.; Diaz, L.A. In Arthropod Borne Diseases; Springer: Berlin/Heidelberg, Germany, 2017; pp. 115–135. [Google Scholar]
  2. Weaver, S.C.; Barrett, A.D.T. Transmission cycles, host range, evolution, and emergence of arboviral disease. Nat. Rev. Genet. 2004, 2, 789–801. [Google Scholar] [CrossRef] [PubMed]
  3. Koolhof, I.S.; Carver, S. Epidemic host community contribution to mosquito-borne disease transmission: Ross River virus. Epidemiol. Infect. 2016, 145, 656–666. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Russell, R.C. Ross River virus: Ecology and distribution. Annu. Rev. Entomol. 2002, 47, 1–31. [Google Scholar] [CrossRef]
  5. Lau, C.L.; Aubry, M.; Nhan, D.M.T.-X.; Teissier, A.; Paulous, S.; Desprès, P.; de Lamballerie, X.; Pastorino, B.; Cao-Lormeau, V.-M.; Weinstein, P. New evidence for endemic circulation of Ross River virus in the Pacific Islands and the potential for emergence. Int. J. Infect. Dis. 2017, 57, 73–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Togami, E.; Gyawali, N.; Ong, O.; Kama, M.; Cao-Lormeau, V.-M.; Aubry, M.; Ko, A.I.; Nilles, E.J.; Collins-Emerson, J.M.; Devine, G.J.; et al. First evidence of concurrent enzootic and endemic transmission of Ross River virus in the absence of marsupial reservoirs in Fiji. Int. J. Infect. Dis. 2020, 96, 94–96. [Google Scholar] [CrossRef]
  7. NNDSS. National Notifiable Diseases Surveillance Fortnightly Summary Notes—4 August to 17 August 2018. Available online: https://www1.health.gov.au/internet/main/publishing.nsf/Content/cdnareport-fn16-18.htm (accessed on 18 May 2020).
  8. Ratnayake, J.T.B. The valuation of social and economic costs of mosquito-transmitted Ross River virus; Griffith University: Queensland, Australia, 2006. [Google Scholar]
  9. Tomerini, D.M. The impact of local government mosquito control programs on Ross River virus disease in Queensland, Australia; Griffith University: Queensland, Australia, 2008. [Google Scholar]
  10. Doherty, R.; Whitehead, R.; Gorman, B.; O’gower, A. The isolation of a third group A arbovirus in Australia, with preliminary observations on its relationship to epidemic polyarthritis. Aust. J. Sci. 1963, 26, 183–184. [Google Scholar]
  11. Russell, R. Ross River virus: Disease trends and vector ecology in Australia. Bull. Soc. Vector Ecol. 1994, 19, 73–81. [Google Scholar]
  12. Ritchie, S.A.; Phillips, D.A.; Standfast, H.A.; Fanning, I.D.; McGinn, D.; Kay, B.H. Ross River Virus in Mosquitoes (Diptera: Culicidae) During the 1994 Epidemic Around Brisbane, Australia. J. Med Entomol. 1997, 34, 156–159. [Google Scholar] [CrossRef]
  13. Kay, B.H.; Ryan, P.A.; Hall, R.A.; Boyd, A.M. Mosquito feeding patterns and natural infection of vertebrates with Ross River and Barmah Forest viruses in Brisbane, Australia. Am. J. Trop. Med. Hyg. 2007, 76, 417–423. [Google Scholar] [CrossRef]
  14. Boyd, A.M.; Hall, R.A.; Gemmell, R.T.; Kay, B.H. Experimental infection of Australian brushtail possums, Trichosurus vulpecula (Phalangeridae: Marsupialia), with Ross River and Barmah Forest viruses by use of a natural mosquito vector system. Am. J. Trop. Med. Hyg. 2001, 65, 777–782. [Google Scholar] [CrossRef] [Green Version]
  15. Holder, P.; Browne, G.; Bullians, M. The mosquitoes of New Zealand and their animal disease significance. Surveillance 1999, 26, 12–15. [Google Scholar]
  16. Derraik, J.G. A survey of the mosquito (Diptera: Culicidae) fauna of the Auckland Zoological Park. New Zealand Èntomol. 2004, 27, 51–55. [Google Scholar] [CrossRef]
  17. Kelly-Hope, L.A.; Kay, B.H.; Purdie, D.M.; Williams, G. Infectious Disease: The risk of Ross River and Barmah Forest virus disease in Queensland: Implications for New Zealand. Aust. New Zealand J. Public Heal. 2002, 26, 69–77. [Google Scholar] [CrossRef]
  18. Gubler, D.J. Transmission of Ross River virus by Aedes polynesiensis and Aedes aegypti. Am. J. Trop. Med. Hyg. 1981, 30, 1303–1306. [Google Scholar] [CrossRef] [PubMed]
  19. Nasci, R.S.; Mitchell, C.J. Larval Diet, Adult Size, and Susceptibility of Aedes aegypti (Diptera: Culicidae) to Infection with Ross River Virus. J. Med. Èntomol. 1994, 31, 123–126. [Google Scholar] [CrossRef] [PubMed]
  20. Whelan, P.I.; Merianos, A.; Patel, M.; Tai, K.; Currie, B. The epidemiology of arbovirus infection in the Northern Territory 1980–1992. Arbovirus Res. Aust. 1993, 6, 266–270. [Google Scholar]
  21. Taylor, W.P.; Marshall, I.D. Adaptation studies with Ross River virus: Retention of field level virulence. J. Gen. Virol. 1975, 28, 73–83. [Google Scholar] [CrossRef] [Green Version]
  22. Bhatt, S.; Gething, P.W.; Brady, O.J.; Messina, J.P.; Farlow, A.W.; Moyes, C.L.; Drake, J.M.; Brownstein, J.S.; Hoen, A.G.; Sankoh, O.; et al. The global distribution and burden of dengue. Nature 2013, 496, 504–507. [Google Scholar] [CrossRef] [PubMed]
  23. Nhan, D.M.T.-X.; Cao-Lormeau, V.-M.; Gubler, D.J. Zika virus: Following the path of dengue and chikungunya? Lancet 2015, 386, 243–244. [Google Scholar] [CrossRef]
  24. Powell, J.R. Mosquito-borne human viral diseases: Why Aedes aegypti? Am. J. Trop. Med. Hyg. 2018, 98, 1563–1565. [Google Scholar] [CrossRef]
  25. Tabachnick, W.J. Evolutionary genetics and arthropod-borne disease: The Yellow Fever mosquito. Am. Èntomol. 1991, 37, 14–26. [Google Scholar] [CrossRef]
  26. Sim, C.; Hong, Y.S.; van Landingham, D.L.; Harker, B.W.; Christophides, G.K.; Kafatos, F.C.; Higgs, S.; Collins, F.H. Modulation of Anopheles gambiae gene expression in response to o’nyong-nyong virus infection. Insect Mol. Biol. 2005, 14, 475–481. [Google Scholar] [CrossRef] [PubMed]
  27. Sim, C.; Hong, Y.S.; Ttsetsarkin, K.; van Landingham, D.L.; Higgs, S.; Collins, F.H. Anopheles gambiae heat shock protein cognate 70B impedes o’nyong-nyong virus replication. BMC Genom. 2007, 8, 231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Xi, Z.; Ramirez, J.L.; Dimopoulos, G. The Aedes aegypti Toll pathway controls dengue virus infection. PLoS Pathog. 2008, 4, e1000098. [Google Scholar] [CrossRef] [PubMed]
  29. Sanders, H.R.; Foy, B.D.; Evans, A.M.; Ross, L.S.; Beaty, B.J.; Olson, K.E.; Gill, S.S. Sindbis virus induces transport processes and alters expression of innate immunity pathway genes in the midgut of the disease vector, Aedes aegypti. Insect Biochem. Mol. Biol. 2005, 35, 1293–1307. [Google Scholar] [CrossRef] [PubMed]
  30. Luplertlop, N.; Surasombatpattana, P.; Patramool, S.; Dumas, E.; Wasinpiyamongkol, L.; Sauné, L.; Hamel, R.; Bernard, E.; Sereno, D.; Thomas, F.; et al. Induction of a peptide with activity against a broad spectrum of pathogens in the Aedes aegypti salivary gland, following infection with dengue virus. PLoS Pathog. 2011, 7, e1001252. [Google Scholar] [CrossRef]
  31. McFarlane, M.; Arias-Goeta, C.; Martin, E.; O’Hara, Z.; Lulla, A.; Mousson, L.; Rainey, S.M.; Misbah, S.; Schnettler, E.; Donald, C.L.; et al. Characterization of Aedes aegypti Innate-Immune Pathways that Limit Chikungunya Virus Replication. PLoS Neglect. Trop. Dis. 2014, 8, e2994. [Google Scholar] [CrossRef]
  32. Jupatanakul, N.; Sim, S.; Angleró-Rodríguez, Y.I.; Souza-Neto, J.A.; Das, S.; Poti, K.E.; Rossi, S.L.; Bergren, N.; Vasilakis, N.; Dimopoulos, G. Engineered Aedes aegypti JAK/STAT pathway-mediated immunity to dengue virus. PLoS Neglect. Trop. Dis. 2017, 11, e0005187. [Google Scholar] [CrossRef]
  33. Angleró-Rodríguez, Y.I.; MacLeod, H.J.; Kang, S.; Carlson, J.S.; Jupatanakul, N.; Dimopoulos, G. Aedes aegypti molecular responses to Zika virus: Modulation of infection by the Toll and Jak/Stat immune pathways and virus host factors. Front. Microbiol. 2017, 8, 2050. [Google Scholar] [CrossRef] [Green Version]
  34. Hussain, M.; Etebari, K.; Asgari, S. Functions of small RNAs in mosquitoes. In Advances in Insect Physiology; Elsevier BV: Amsterdam, The Netherlands, 2016; Volume 51, pp. 189–222. [Google Scholar]
  35. Blair, C.D.; Olson, K.E. The Role of RNA interference (RNAi) in arbovirus-vector interactions. Viruses 2015, 7, 820–843. [Google Scholar] [CrossRef] [Green Version]
  36. Campbell, C.L.; Keene, K.M.; Brackney, D.E.; Olson, K.E.; Blair, C.D.; Wilusz, J.; Foy, B.D. Aedes aegypti uses RNA interference in defense against Sindbis virus infection. BMC Microbiol. 2008, 8, 47. [Google Scholar] [CrossRef] [Green Version]
  37. Cirimotich, C.; Scott, J.C.; Phillips, A.T.; Geiss, B.J.; Olson, K.E. Suppression of RNA interference increases alphavirus replication and virus-associated mortality in Aedes aegypti mosquitoes. BMC Microbiol. 2009, 9, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Franz, A.W.; Sanchez-Vargas, I.; Adelman, Z.N.; Blair, C.D.; Beaty, B.J.; James, A.A.; Olson, K.E. Engineering RNA interference-based resistance to dengue virus type 2 in genetically modified Aedes aegypti. Proc. Natl. Acad. Sci. USA 2006, 103, 4198–4203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Siu, R.W.C.; Fragkoudis, R.; Simmonds, P.; Donald, C.L.; Chase-Topping, M.E.; Barry, G.; Attarzadeh-Yazdi, G.; Rodriguez-Andres, J.; Nash, A.A.; Merits, A.; et al. Antiviral RNA interference responses induced by Semliki Forest virus infection of mosquito cells: Characterization, origin, and frequency-dependent functions of virus-derived small interfering RNAs. J. Virol. 2010, 85, 2907–2917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Bronkhorst, A.W.; van Rij, R.P. The long and short of antiviral defense: Small RNA-based immunity in insects. Curr. Opin. Virol. 2014, 7, 19–28. [Google Scholar] [CrossRef]
  41. Arensburger, P.; Hice, R.H.; Wright, J.A.; Craig, N.L.; Atkinson, P.W. The mosquito Aedes aegypti has a large genome size and high transposable element load but contains a low proportion of transposon-specific piRNAs. BMC Genom. 2011, 12, 606. [Google Scholar] [CrossRef] [Green Version]
  42. Hussain, M.; Walker, T.; O’Neill, S.L.; Asgari, S. Blood meal induced microRNA regulates development and immune associated genes in the Dengue mosquito vector, Aedes aegypti. Insect Biochem. Mol. Biol. 2013, 43, 146–152. [Google Scholar] [CrossRef]
  43. Liu, S.; Lucas, K.J.; Roy, S.; Ha, J.; Raikhel, A.S. Mosquito-specific microRNA-1174 targets serine hydroxymethyltransferase to control key functions in the gut. Proc. Natl. Acad. Sci. USA 2014, 111, 14460–14465. [Google Scholar] [CrossRef] [Green Version]
  44. Bryant, B.; Macdonald, W.; Raikhel, A.S. microRNA miR-275 is indispensable for blood digestion and egg development in the mosquito Aedes aegypti. Proc. Natl. Acad. Sci. USA 2010, 22391–22398. [Google Scholar] [CrossRef] [Green Version]
  45. Campbell, C.L.; Harrison, T.; Hess, A.M.; Ebel, G.D. MicroRNA levels are modulated in Aedes aegypti after exposure to D engue-2. Insect. Mol. Biol. 2014, 23, 132–139. [Google Scholar] [CrossRef] [Green Version]
  46. Maharaj, P.D.; Widen, S.G.; Huang, J.; Wood, T.G.; Thangamani, S. Discovery of mosquito saliva microRNAs during CHIKV infection. PLoS Neglect. Trop. Dis. 2015, 9, e0003386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Arensburger, P.; Megy, K.; Waterhouse, R.M.; Abrudan, J.; Amedeo, P.; Antelo, B.; Bartholomay, L.; Bidwell, S.; Caler, E.; Camara, F.; et al. Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics. Science 2010, 330, 86–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kay, B.H.; Fanning, I.D.; Carley, J.G. Vector competence of Culex pipiens quinquefasciatus for Murray Valley Encephalitis, Kunjin, and Ross River viruses from Australia. Am. J. Trop. Med. Hyg. 1982, 31, 844–848. [Google Scholar] [CrossRef] [PubMed]
  49. Pearson, W.; Wood, T.; Zhang, Z.; Miller, W. Comparison of DNA sequences with protein sequences. Genomics 1997, 46, 24–36. [Google Scholar] [CrossRef] [Green Version]
  50. Dodt, M.; Roehr, J.T.; Ahmed, R.; Dieterich, C. FLEXBAR—Flexible barcode and adapter processing for next-generation sequencing platforms. Biology 2012, 1, 895–905. [Google Scholar] [CrossRef] [Green Version]
  51. Wickham, H.; Seidel, D. Scales: scale functions for visualization. Available online: https://CRAN.R-project.org/package=scales (accessed on 18 May 2020).
  52. Wickham, H. The split-apply-combine strategy for data analysis. J. Stat. Softw. 2011, 40, 1–29. [Google Scholar] [CrossRef] [Green Version]
  53. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  54. Shi, J.; Dong, M.; Li, L.; Luz-Madrigal, A.; Tsonis, P.A.; del Rio-Tsonis, K.; Liang, C. mirPRo—A novel standalone program for differential expression and variation analysis of miRNAs. Sci. Rep. 2015, 5, 14617. [Google Scholar] [CrossRef] [Green Version]
  55. Griffiths-Jones, S. miRBase: microRNA sequences, targets, and gene nomenclature. Nucleic Acids Res. 2006, 34, D140–D144. [Google Scholar] [CrossRef]
  56. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2014, 31, 166–169. [Google Scholar] [CrossRef]
  57. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.R.; Durbin, R. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  58. Giraldo-Calderón, G.I.; Emrich, S.J.; Maccallum, R.M.; Maslen, G.; Dialynas, E.; Topalis, P.; Ho, N.; Gesing, S.; Madey, G.; Collins, F.H.; et al. VectorBase: An updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res. 2014, 43, D707–D713. [Google Scholar] [CrossRef]
  59. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 002832. [Google Scholar] [CrossRef] [Green Version]
  60. Huber, W.; Carey, V.J.; Gentleman, R.; Anders, S.; Carlson, M.; Carvalho, B.S.; Bravo, H.C.; Davis, S.; Gatto, L.; Girke, T.; et al. Orchestrating high-throughput genomic analysis with Bioconductor. Nat. Methods 2015, 12, 115–121. [Google Scholar] [CrossRef] [PubMed]
  61. Anonymous. The R project for statistical computing. Available online: http://www.r-project.org/ (accessed on 13 February 2012).
  62. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  63. Zhu, A.; Ibrahim, J.G.; Love, M.I. Heavy-tailed prior distributions for sequence count data: Removing the noise and preserving large differences. Bioinformatics 2019, 35, 2084–2092. [Google Scholar] [CrossRef] [PubMed]
  64. Jones, A.; Lowry, K.; Aaskov, J.; Holmes, E.; Kitchen, A. Molecular evolutionary dynamics of Ross River virus and implications for vaccine efficacy. J. Gen. Virol. 2009, 91, 182–188. [Google Scholar] [CrossRef] [Green Version]
  65. Gel, B.; Serra, E. KaryoploteR: An R/Bioconductor package to plot customizable genomes displaying arbitrary data. Bioinformatics 2017, 33, 3088–3090. [Google Scholar] [CrossRef]
  66. Crooks, G.E.; Hon, G.; Chandonia, J.-M.; Brenner, S. WebLogo: A Sequence Logo Generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [Green Version]
  67. Afgan, E.; Baker, D.; Batut, B.; Beek, M.V.D.; Bouvier, D.; Čech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy platform for accessible, reproducible, and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018, 46, W537–W544. [Google Scholar] [CrossRef] [Green Version]
  68. Enright, A.J.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D.S. MicroRNA targets in Drosophila. Genome Biol. 2003, 5, R1. [Google Scholar] [CrossRef] [Green Version]
  69. Krüger, J.; Rehmsmeier, M. RNAhybrid: microRNA target prediction easy, fast, and flexible. Nucleic Acids Res. 2006, 34, W451–W454. [Google Scholar] [CrossRef] [PubMed]
  70. Agarwal, V.; Subtelny, A.O.; Thiru, P.; Ulitsky, I.; Bartel, B. Predicting microRNA targeting efficacy in Drosophila. Genome Biol. 2018, 19, 152. [Google Scholar] [CrossRef] [Green Version]
  71. Kertesz, M.; Iovino, N.; Unnerstall, U.; Gaul, U.; Segal, E. The role of site accessibility in microRNA target recognition. Nat. Genet. 2007, 39, 1278–1284. [Google Scholar] [CrossRef] [PubMed]
  72. Loher, P.; Rigoutsos, I. Interactive exploration of RNA22 microRNA target predictions. Bioinformatics 2012, 28, 3322–3323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Oliveira, A.C.; Bovolenta, L.A.; Nachtigall, P.G.; Herkenhoff, M.; Lemke, N.; Pinhal, D. Combining results from distinct microrna target prediction tools enhances the performance of analyses. Front. Genet. 2017, 8, 59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef]
  76. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29, 45. [Google Scholar] [CrossRef]
  77. Trobaugh, D.W.; Klimstra, W.B. MicroRNA regulation of rna virus replication and pathogenesis. Trends Mol. Med. 2016, 23, 80–93. [Google Scholar] [CrossRef]
  78. Saldana, M.A.; Etebari, K.; Hart, C.E.; Widen, S.G.; Wood, T.G.; Thangamani, S.; Asgari, S.; Hughes, G.L. Zika virus alters the microRNA expression profile and elicits an RNAi response in Aedes aegypti mosquitoes. PLoS Neglect. Trop. Dis. 2017, 11, e0005760. [Google Scholar] [CrossRef]
  79. Shrinet, J.; Jain, S.; Jain, J.; Bhatnagar, R.K.; Sunil, S. Next generation sequencing reveals regulation of distinct Aedes microRNAs during chikungunya virus development. PLoS Neglect. Trop. Dis. 2014, 8, e2616. [Google Scholar] [CrossRef] [PubMed]
  80. Myles, K.M.; Wiley, M.R.; Morazzani, E.M.; Adelman, Z.N. Alphavirus-derived small RNAs modulate pathogenesis in disease vector mosquitoes. Proc. Natl. Acad. Sci. USA 2008, 105, 19938–19943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Pakpour, N.; Riehle, M.A.; Luckhart, S. Effects of ingested vertebrate-derived factors on insect immune responses. Curr. Opin. Insect Sci. 2014, 3, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Liu, Y.; Zhou, Y.; Wu, J.; Zheng, P.; Li, Y.; Zheng, X.; Puthiyakunnon, S.; Tu, Z.; Chen, X.-G. The expression profile of Aedes albopictus miRNAs is altered by dengue virus serotype-2 infection. Cell Biosci. 2015, 5, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Etebari, K.; Asgari, S. Accuracy of microrna discovery pipelines in non-model organisms using closely related species genomes. PLoS ONE 2014, 9, e84747. [Google Scholar] [CrossRef] [Green Version]
  84. Winter, F.; Edaye, S.; Hüttenhofer, A.; Brunel, C. Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucleic Acids Res. 2007, 35, 6953–6962. [Google Scholar] [CrossRef] [Green Version]
  85. Akbari, O.S.; Antoshechkin, I.; Amrhein, H.; Williams, B.; Diloreto, R.; Sandler, J.; Hay, B.A. The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector. G3 Genes/Genomes/Genetics 2013, 3, 1493–1509. [Google Scholar] [CrossRef] [Green Version]
  86. Mead, E.A.; Tu, Z. Cloning, characterization, and expression of microRNAs from the Asian malaria mosquito, Anopheles stephensi. BMC Genom. 2008, 9, 244. [Google Scholar] [CrossRef] [Green Version]
  87. Skalsky, R.L.; van Landingham, D.L.; Scholle, F.; Higgs, S.; Cullen, B.R. Identification of microRNAs expressed in two mosquito vectors, Aedes albopictus and Culex quinquefasciatus. BMC Genom. 2010, 11, 119. [Google Scholar] [CrossRef] [Green Version]
  88. Etebari, K.; Osei-Amo, S.; Blomberg, S.P.; Asgari, S. Dengue virus infection alters post-transcriptional modification of microRNAs in the mosquito vector Aedes aegypti. Sci. Rep. 2015, 5, 15968. [Google Scholar] [CrossRef] [Green Version]
  89. Mayoral, J.G.; Etebari, K.; Hussain, M.; Khromykh, A.A.; Asgari, S. Wolbachia infection modifies the profile, shuttling and structure of microRNAs in a mosquito cell line. PLoS ONE 2014, 9, e96107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Hedges, L.M.; Brownlie, J.C.; O’Neill, S.L.; Johnson, K.; Pörtner, H.-O.; Farrell, A.P. Wolbachia and virus protection in insects. Science 2008, 322, 702. [Google Scholar] [CrossRef] [PubMed]
  91. Asad, S.; Hussain, M.; Hugo, L.E.; Osei-Amo, S.; Zhang, G.; Watterson, D.; Asgari, S. Suppression of the pelo protein by Wolbachia and its effect on dengue virus in Aedes aegypti. PLoS Neglect. Trop. Dis. 2018, 12, e0006405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Zhang, G.; Hussain, M.; O’Neill, S.L.; Asgari, S. Wolbachia uses a host microRNA to regulate transcripts of a methyltransferase, contributing to dengue virus inhibition in Aedes aegypti. Proc. Natl. Acad. Sci. USA 2013, 110, 10276–10281. [Google Scholar] [CrossRef] [Green Version]
  93. Bick, M.J.; Carroll, J.-W.N.; Gao, G.; Goff, S.P.; Rice, C.M.; Macdonald, M.R. expression of the zinc-finger antiviral protein inhibits alphavirus replication. J. Virol. 2003, 77, 11555–11562. [Google Scholar] [CrossRef] [Green Version]
  94. Nakamoto, M.; Moy, R.; Xu, J.; Bambina, S.; Yasunaga, A.; Shelly, S.S.; Gold, B.; Cherry, S. Virus recognition by Toll-7 activates antiviral autophagy in Drosophila. Immunity 2012, 36, 658–667. [Google Scholar] [CrossRef] [Green Version]
  95. Patel, R.K.; Hardy, R.W. Role for the phosphatidylinositol 3-Kinase-Akt-TOR pathway during Sindbis virus Replication in arthropods. J. Virol. 2012, 86, 3595–3604. [Google Scholar] [CrossRef] [Green Version]
  96. Mazzon, M.; Castro, C.; Thaa, B.; Liu, L.; Mutso, M.; Liu, X.; Mahalingam, S.; Griffin, J.L.; Marsh, M.; McInerney, G. Alphavirus-induced hyperactivation of PI3K/AKT directs pro-viral metabolic changes. PLoS Pathog. 2018, 14, e1006835. [Google Scholar] [CrossRef] [Green Version]
  97. Gulia-Nuss, M.; Robertson, A.E.; Brown, M.R.; Strand, M.R. insulin-like peptides and the target of rapamycin pathway coordinately regulate blood digestion and egg maturation in the mosquito Aedes aegypti. PLoS ONE 2011, 6, e20401. [Google Scholar] [CrossRef]
  98. Lucas, K.J.; Roy, S.; Ha, J.; Gervaise, A.L.; Kokoza, V.A.; Raikhel, A.S. MicroRNA-8 targets the Wingless signaling pathway in the female mosquito fat body to regulate reproductive processes. Proc. Natl. Acad. Sci. USA 2015, 112, 1440–1445. [Google Scholar] [CrossRef] [Green Version]
  99. Schwenke, R.A.; Lazzaro, B.P.; Wolfner, M.F. Reproduction-Immunity Trade-Offs in Insects. Annu. Rev. Entomol. 2015, 61, 239–256. [Google Scholar] [CrossRef] [Green Version]
  100. Kuss-Duerkop, S.K.; Wang, J.; Mena, I.; White, K.M.; Metreveli, G.; Sakthivel, R.; Mata, M.A.; Muñoz-Moreno, R.; Chen, X.; Krammer, F.; et al. Influenza virus differentially activates mTORC1 and mTORC2 signaling to maximize late stage replication. PLoS Pathog. 2017, 13, e1006635. [Google Scholar] [CrossRef] [PubMed]
  101. Liang, Q.; Luo, Z.; Zeng, J.; Chen, W.; Foo, S.-S.; Lee, S.-A.; Ge, J.; Wang, S.; Goldman, S.A.; Zlokovic, B.V.; et al. Zika virus NS4A and NS4B proteins deregulate Akt-mTOR signaling in human fetal neural stem cells to inhibit neurogenesis and induce autophagy. Cell Stem Cell 2016, 19, 663–671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Dong, S.; Behura, S.K.; Franz, A.W.E. The midgut transcriptome of Aedes aegypti fed with saline or protein meals containing chikungunya virus reveals genes potentially involved in viral midgut escape. BMC Genom. 2017, 18, 382. [Google Scholar] [CrossRef] [Green Version]
  103. Etebari, K.; Hegde, S.; Saldaña, M.A.; Widen, S.G.; Wood, T.G.; Asgari, S.; Hughes, G.L. global transcriptome analysis of Aedes aegypti mosquitoes in response to zika virus infection. mSphere 2017, 2, e00456-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Raquin, V.; Merkling, S.H.; Gausson, V.; Moltini-Conclois, I.; Frangeul, L.; Varet, H.; Dillies, M.-A.; Saleh, M.C.; Lambrechts, L. Individual co-variation between viral RNA load and gene expression reveals novel host factors during early dengue virus infection of the Aedes aegypti midgut. PLoS Neglect. Trop. Dis. 2017, 11, e0006152. [Google Scholar] [CrossRef] [PubMed]
  105. Zhao, L.; Alto, B.W.; Jiang, Y.; Yu, F.; Zhang, Y. Transcriptomic analysis of Aedes aegypti innate immune system in response to ingestion of Chikungunya virus. Int. J. Mol. Sci. 2019, 20, 3133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Khoo, C.C.H.; Piper, J.; Sanchez-Vargas, I.; Olson, K.E.; Franz, A.W. The RNA interference pathway affects midgut infection- and escape barriers for Sindbis virus in Aedes aegypti. BMC Microbiol. 2010, 10, 130. [Google Scholar] [CrossRef] [Green Version]
  107. Joosten, J.; Miesen, P.; Taşköprü, E.; Pennings, B.; Jansen, P.W.T.C.; Huynen, M.A.; Vermeulen, M.; van Rij, R.P. The Tudor protein Veneno assembles the ping-pong amplification complex that produces viral piRNAs in Aedes mosquitoes. Nucleic Acids Res. 2019, 47, 2546–2559. [Google Scholar] [CrossRef] [Green Version]
  108. Miesen, P.; Girardi, E.; van Rij, R.P. Distinct sets of PIWI proteins produce arbovirus and transposon-derived piRNAs in Aedes aegypti mosquito cells. Nucleic Acids Res. 2015, 43, 6545–6556. [Google Scholar] [CrossRef] [Green Version]
  109. Vodovar, N.; Bronkhorst, A.W.; van Cleef, K.W.R.; Miesen, P.; Blanc, H.; van Rij, R.P.; Saleh, M.C. Arbovirus-derived piRNAs exhibit a ping-pong signature in mosquito cells. PLoS ONE 2012, 7, e30861. [Google Scholar] [CrossRef] [PubMed]
  110. Brennecke, J.; Aravin, A.A.; Stark, A.; Dus, M.; Kellis, M.; Sachidanandam, R.; Hannon, G.J. Discrete Small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell 2007, 128, 1089–1103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Gunawardane, L.S.; Saito, K.; Nishida, K.M.; Miyoshi, K.; Kawamura, Y.; Nagami, T.; Siomi, H.; Siomi, M.C. A slicer-mediated mechanism for repeat-associated siRNA 5’ end formation in Drosophila. Science 2007, 315, 1587–1590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Miesen, P.; Joosten, J.; van Rij, R.P. PIWIs Go Viral: Arbovirus-derived piRNAs in vector mosquitoes. PLoS Pathog. 2016, 12, e1006017. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Measurement of Ross River virus (RRV) infection in Aedes aegypti mosquitoes at 2, 6, and 12 days post-infection (dpi) by RT-qPCR. (A) Rate of infection: percentage of mosquitoes with a detectable degree of RRV infection (N = 10), and (B) their relative viral load (N = 6). Error bars indicate the standard error of the mean of normalized expression values. The Ae. aegypti ribosomal protein subunit 17, RPS17, gene was used as the internal calibrator.
Figure 1. Measurement of Ross River virus (RRV) infection in Aedes aegypti mosquitoes at 2, 6, and 12 days post-infection (dpi) by RT-qPCR. (A) Rate of infection: percentage of mosquitoes with a detectable degree of RRV infection (N = 10), and (B) their relative viral load (N = 6). Error bars indicate the standard error of the mean of normalized expression values. The Ae. aegypti ribosomal protein subunit 17, RPS17, gene was used as the internal calibrator.
Viruses 12 00695 g001
Figure 2. Length distribution of quality-filtered and adapter trimmed reads from Ae. aegypti fat body and midgut tissues at 2, 6, and 12 dpi. Error bars indicate the standard error of the mean of biological replicates (N = 3).
Figure 2. Length distribution of quality-filtered and adapter trimmed reads from Ae. aegypti fat body and midgut tissues at 2, 6, and 12 dpi. Error bars indicate the standard error of the mean of biological replicates (N = 3).
Viruses 12 00695 g002
Figure 3. RRV caused differential expression of Ae. aegypti miRNAs. Log2 fold changes of miRNA expression following infection with RRV at 2 dpi in fat body (green) and midgut (yellow) and in fat body at 12 dpi (purple). Error bars represent the standard error of the mean of biological replicates (N = 3).
Figure 3. RRV caused differential expression of Ae. aegypti miRNAs. Log2 fold changes of miRNA expression following infection with RRV at 2 dpi in fat body (green) and midgut (yellow) and in fat body at 12 dpi (purple). Error bars represent the standard error of the mean of biological replicates (N = 3).
Viruses 12 00695 g003
Figure 4. RRV caused differential expression of Ae. aegypti miRNAs mostly early in infection. Volcano plots show miRNA log2 fold changes due to RRV infection in fat body at (A) 2 dpi, (B) 6 dpi, and (C) 12 dpi, and in midgut at (D) 2 dpi, (E) 6 dpi, and (F) 12 dpi. Significantly differentially expressed miRNAs (orange) were found in fat body samples at 2 dpi and 12 dpi and in midgut samples at 2 dpi.
Figure 4. RRV caused differential expression of Ae. aegypti miRNAs mostly early in infection. Volcano plots show miRNA log2 fold changes due to RRV infection in fat body at (A) 2 dpi, (B) 6 dpi, and (C) 12 dpi, and in midgut at (D) 2 dpi, (E) 6 dpi, and (F) 12 dpi. Significantly differentially expressed miRNAs (orange) were found in fat body samples at 2 dpi and 12 dpi and in midgut samples at 2 dpi.
Viruses 12 00695 g004
Figure 5. Ae. aegypti miRNA differential expression due to RRV infection between tissues and across time. The expression profile of (A) miR-9b-5p and (B) miR-71-3p, in control (blue) and infected (red) fat body and midgut samples at 2, 6, and 12 dpi. Significant changes are labelled with the fold change (FC) and adjusted p-value (padj). The boxplot and whiskers show the mean of the normalized reads counts of the biological replicates (N = 3).
Figure 5. Ae. aegypti miRNA differential expression due to RRV infection between tissues and across time. The expression profile of (A) miR-9b-5p and (B) miR-71-3p, in control (blue) and infected (red) fat body and midgut samples at 2, 6, and 12 dpi. Significant changes are labelled with the fold change (FC) and adjusted p-value (padj). The boxplot and whiskers show the mean of the normalized reads counts of the biological replicates (N = 3).
Viruses 12 00695 g005
Figure 6. RRV induced differential expression of Ae. aegypti miRNAs. Contrasting RNA-Seq and RT-qPCR measurement of RRV-induced differential expression of 10 miRNAs in Ae. aegypti fat body and midgut RNA samples at 2 dpi. Error bars represent the standard error of the mean of the biological replicates (N = 3).
Figure 6. RRV induced differential expression of Ae. aegypti miRNAs. Contrasting RNA-Seq and RT-qPCR measurement of RRV-induced differential expression of 10 miRNAs in Ae. aegypti fat body and midgut RNA samples at 2 dpi. Error bars represent the standard error of the mean of the biological replicates (N = 3).
Viruses 12 00695 g006
Figure 7. Venn diagram showing the number of predicted targets for each miRNA, wherein the intersect shows the number of predicted targets in common.
Figure 7. Venn diagram showing the number of predicted targets for each miRNA, wherein the intersect shows the number of predicted targets in common.
Viruses 12 00695 g007
Figure 8. Molecular function. Tree map of gene ontology descriptions ascribed to genes from this study showing the interrelationship between the parent (blue) and children (shades of orange) ontologies involved in the molecular functions of binding and catalytic activities.
Figure 8. Molecular function. Tree map of gene ontology descriptions ascribed to genes from this study showing the interrelationship between the parent (blue) and children (shades of orange) ontologies involved in the molecular functions of binding and catalytic activities.
Viruses 12 00695 g008
Figure 9. Molecular function. The number and percentage of genes from this study ascribed to various molecular functions.
Figure 9. Molecular function. The number and percentage of genes from this study ascribed to various molecular functions.
Viruses 12 00695 g009
Figure 10. The proliferating 21 nt short interfering RNA (siRNA) count over time was indicative of a progressively increasing viral replication and RNAi response in Ae. aegypti mosquitoes. These time-series plots show the (A) length distribution of reads mapping to the RRV (T48) genome at 2, 6, and 12 dpi, and (B) the 21 nt siRNA coverage of sense (blue lines) and anti-sense (red lines) RRV (T48) genome strands. Error bars represent the SEM of the biological replicates (N = 3).
Figure 10. The proliferating 21 nt short interfering RNA (siRNA) count over time was indicative of a progressively increasing viral replication and RNAi response in Ae. aegypti mosquitoes. These time-series plots show the (A) length distribution of reads mapping to the RRV (T48) genome at 2, 6, and 12 dpi, and (B) the 21 nt siRNA coverage of sense (blue lines) and anti-sense (red lines) RRV (T48) genome strands. Error bars represent the SEM of the biological replicates (N = 3).
Viruses 12 00695 g010
Figure 11. RRV-derived PIWI-interacting RNAs (piRNAs) generated by the fat body and midgut tissues of Ae. aegypti mosquitoes. Distribution of 27–29 nt small RNAs that mapped across the sense (blue) and anti-sense (red) strands of the RRV genome at 2, 6, and 12 dpi. Distinct 3’ hotspots indicate virus-derived piRNA (vpiRNA) production against the subgenomic RNA.
Figure 11. RRV-derived PIWI-interacting RNAs (piRNAs) generated by the fat body and midgut tissues of Ae. aegypti mosquitoes. Distribution of 27–29 nt small RNAs that mapped across the sense (blue) and anti-sense (red) strands of the RRV genome at 2, 6, and 12 dpi. Distinct 3’ hotspots indicate virus-derived piRNA (vpiRNA) production against the subgenomic RNA.
Viruses 12 00695 g011
Figure 12. RRV vpiRNA signatures in Ae. aegypti. (A) There was an A10 nucleotide bias, indicative of ping-pong amplification, present in the 28 nt sense vpiRNA-like strands produced in Ae. aegypti fat body tissue at 6 dpi with RRV, although there were other biases at positions 1 and 26–28. (B) The U1 bias, a signature of piRNA ping-pong amplification, was evident in 28 nt piRNAs produced in Ae. aegypti fat body at 6 dpi.
Figure 12. RRV vpiRNA signatures in Ae. aegypti. (A) There was an A10 nucleotide bias, indicative of ping-pong amplification, present in the 28 nt sense vpiRNA-like strands produced in Ae. aegypti fat body tissue at 6 dpi with RRV, although there were other biases at positions 1 and 26–28. (B) The U1 bias, a signature of piRNA ping-pong amplification, was evident in 28 nt piRNAs produced in Ae. aegypti fat body at 6 dpi.
Viruses 12 00695 g012
Figure 13. Sequences in the range 27–29 nt from Ae. aegypti fat body at 6 days post infection with RRV exhibited a high 10 nt overlap probability. (A) A tally of sequence pairs by the number of their overlapping nucleotides, and (B) the probability that sequence pairs overlap by the number of intersecting nucleotides. The z-scores (axis on the right) show the mean +/− the number of standard deviations away from the mean. The number of sequences that overlap by 10 nucleotides and the probability that sequence pairs overlap by 10 nt were higher than two standard deviations from the mean. These peaks at 10 nt are a signature of piRNA production by ping-pong amplification.
Figure 13. Sequences in the range 27–29 nt from Ae. aegypti fat body at 6 days post infection with RRV exhibited a high 10 nt overlap probability. (A) A tally of sequence pairs by the number of their overlapping nucleotides, and (B) the probability that sequence pairs overlap by the number of intersecting nucleotides. The z-scores (axis on the right) show the mean +/− the number of standard deviations away from the mean. The number of sequences that overlap by 10 nucleotides and the probability that sequence pairs overlap by 10 nt were higher than two standard deviations from the mean. These peaks at 10 nt are a signature of piRNA production by ping-pong amplification.
Viruses 12 00695 g013
Table 1. Differentially expressed Ae. aegypti miRNAs upon RRV infection.
Table 1. Differentially expressed Ae. aegypti miRNAs upon RRV infection.
Tissue 1dpimiRNAControl 2RRV 2FC 3p-ValueFDR
FB2miR-9b-5p475210,9392.309.10 × 10−109.55 × 10−8
FB2miR-317-3p22,08413,362−1.654.81 × 10−62.52 × 10−4
FB2miR-275-3p12,4656290−1.983.74 × 10−51.31 × 10−3
FB2miR-275-5p2208752−2.944.39 × 10−49.22 × 10−3
FB2miR-8-5p99265071−1.962.41 × 10−46.32 × 10−3
FB2miR-2945-3p907014,2921.581.40 × 10−30.02
FB2miR-9a-5p477982761.731.56 × 10−30.02
FB2miR-988-3p13402062−1.531.68 × 10−30.02
FB2miR-34-3p852523−1.631.69 × 10−30.02
FB2miR-184-3p149,60293,549−1.601.83 × 10−30.02
FB2miR-2940-5p29,03416,911−1.721.89 × 10−30.02
FB2miR-989-3p11,46867,1385.852.18 × 10−30.02
FB12miR-71-3p498366−2.013.85 × 10−59.58 × 10−3
FB12miR-2940-3p11,49522,2031.934.19 × 10−40.05
MG2miR-9b-5p913015,2361.661.67 × 10−40.02
1 FB: fat body, MG: midgut; 2 base mean counts; 3 FC: fold change.

Share and Cite

MDPI and ACS Style

Sinclair, J.B.; Asgari, S. Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes. Viruses 2020, 12, 695. https://doi.org/10.3390/v12070695

AMA Style

Sinclair JB, Asgari S. Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes. Viruses. 2020; 12(7):695. https://doi.org/10.3390/v12070695

Chicago/Turabian Style

Sinclair, James B., and Sassan Asgari. 2020. "Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes" Viruses 12, no. 7: 695. https://doi.org/10.3390/v12070695

APA Style

Sinclair, J. B., & Asgari, S. (2020). Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes. Viruses, 12(7), 695. https://doi.org/10.3390/v12070695

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop