# RNASeq Analysis of Aedes albopictus Mosquito Midguts after Chikungunya Virus Infection

Ravi Kiran Vedururu, Matthew J. Neave, Mary Tachedjian, Melissa J. Klein, Paul R. Gorry, Jean-Bernard Duchemin, Prasad Paradkar

► To cite this version:

Ravi Kiran Vedururu, Matthew J. Neave, Mary Tachedjian, Melissa J. Klein, Paul R. Gorry, et al. RNASeq Analysis of Aedes albopictus Mosquito Midguts after Chikungunya Virus Infection. Viruses, 2019, 11 (6), pp.513. 10.3390/v11060513. hal-03235177

HAL Id: hal-03235177

https://hal.science/hal-03235177

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# RNASeq Analysis of *Aedes albopictus* Mosquito Midguts after Chikungunya Virus Infection

Ravi kiran Vedururu <fnref n="1" />,<fnref n="2" /> ID, Matthew J. Neave <fnref n="1" />, Mary Tachedjian <fnref n="1" />, Melissa J. Klein <fnref n="1" />, Paul R. Gorry <fnref n="3" />, Jean-Bernard Duchemin <fnref n="1" />,† and Prasad N. Paradkar <fnref n="1" />,* ID

Received: 4 May 2019; Accepted: 2 June 2019; Published: 4 June 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Chikungunya virus (CHIKV) is an emerging pathogen around the world and causes significant morbidity in patients. A single amino acid mutation in the envelope protein of CHIKV has led to a shift in vector preference towards *Aedes albopictus*. While mosquitoes are known to mount an antiviral immune response post-infection, molecular interactions during the course of infection at the tissue level remain largely uncharacterised. We performed whole transcriptome analysis on dissected midguts of *Aedes albopictus* infected with CHIKV to identify differentially expressed genes. For this, RNA was extracted at two days post-infection (2-dpi) from pooled midguts. We initially identified 25 differentially expressed genes (*p*-value < 0.05) when mapped to a reference transcriptome. Further, multiple differentially expressed genes were identified from a custom *de novo* transcriptome, which was assembled using the reads that did not align with the reference genome. Thirteen of the identified transcripts, possibly involved in immunity, were validated by qRT-PCR. Homologues of seven of these genes were also found to be significantly upregulated in *Aedes aegypti* midguts 2 dpi, indicating a conserved mechanism at play. These results will help us to characterise the molecular interaction between *Aedes albopictus* and CHIKV and can be utilised to reduce the impact of this viral infection.

**Keywords:** Chikungunya; *Aedes albopictus*; RNASeq; Host-pathogen interactions

## 1. Introduction

Arboviruses, such as the dengue, chikungunya, and Zika viruses, are a significant burden on public health systems worldwide. These viruses are transmitted by mosquitoes and can cause high morbidity and mortality, with dengue alone causing more than 300 million infections per year [1]. First identified and described in 1955 in a report on an outbreak at the border of Tanzania and Mozambique in 1952, chikungunya virus (CHIKV) is an enveloped, positive sense RNA virus belonging to the alphavirus genus in the *Togaviridae* family [2,3]. CHIKV infection in humans causes high febrile illness, cutaneous exanthema and debilitating and often prolonged arthralgia [4–6].

While *Aedes aegypti* is the traditional vector for CHIKV, *Aedes albopictus* has been observed to be implicated in viral transmission ever since the Reunion Island outbreak in 2005–2006 [7]. A single amino acid change in codon 226 of the E1 gene, which encodes for the envelope protein of the virus, has improved the fitness of the Indian ocean lineage CHIKV in *Aedes albopictus* [8]. As an invasive

<footnote n="1">¹ CSIRO Health & Biosecurity, Australian Animal Health Laboratory, Geelong 3220, Australia; Ravikiran.Vedururu@csiro.au (R.k.V.); Matthew.Neave@csiro.au (M.J.N.); Mary.Tachedjian@csiro.au (M.T.); Melissa.Klein@csiro.au (M.J.K.); Jean-bernard.Duchemin@csiro.au (J.-B.D.)</footnote>

<footnote n="2">² School of Applied Sciences, RMIT University, Bundoora 3083, Australia</footnote>

<footnote n="3">³ School of Health and Biomedical Science, RMIT University, Bundoora 3083, Australia; Paul.Gorry@rmit.edu.au</footnote>

<footnote>* Correspondence: Prasad.Paradkar@csiro.au; Tel.: +61-3-52275462</footnote>

<footnote>† Current Address: Mosquitoes and Emerging Arboviruses, Institut Pasteur de la Guyane, Cayenne 97306, French Guiana.</footnote>
species, *Aedes albopictus* has been expanding its traditional habitat of tropical and sub-tropical regions to much cooler temperate regions. *Aedes albopictus* also survives in favourable microhabitats, even in winter and freezing temperatures [9]. These factors have further increased the risk of CHIKV to cause outbreaks in areas where mosquito-borne viral diseases are uncommon, such as Northern America and temperate Europe [10,11].

For a mosquito to become infective, the virus needs to cross two critical barrier tissues: the midgut and salivary glands. These infection barriers can be influenced by multiple factors, including viral factors, such as viral glycoproteins, or vector factors, such as the presence of a viral receptor, host replication factors and the microbiome composition of the midgut [12]. When an adult female *Aedes albopictus* mosquito is exposed to CHIKV in the process of blood feeding, the virus infects the midgut usually in a matter of hours [13–15]. Inside mosquitoes, after feeding, the blood meal moves down to the midgut where the virus must contact epithelial cells before digestion of the blood meal and formation of the peritrophic matrix takes place. Following the successful infection of the midgut, the virus must then overcome midgut escape barriers to disseminate to other tissues, such as the haemocoel [16]. From the haemocoel, the virus makes its way to the salivary glands of the mosquito. Once the virus is detected in the saliva, the mosquito is considered to be infective and a competent vector [17,18].

Since vector competence of the mosquito is initially dependent on infection of the midgut (Virus crossing the first critical barrier), understanding the molecular interactions between virus and mosquito midgut becomes essential. Although previous studies have used transcriptome-based approaches to identify mosquito-virus interactions, tissue-specific responses during the course of infection are not well-understood [19–24]. Here, using next-generation sequencing, we characterised the whole transcriptome response at the midgut in *Aedes albopictus* in response to CHIKV infection.

## 2. Materials and Methods

### 2.1. Chikungunya Virus

Chikungunya virus isolate 06113879 (Mauritius strain), isolated from a viraemic traveller who returned to Australia in 2006, was obtained from the Victorian Infectious Diseases Reference Laboratory (VIDRL), Melbourne [25]. The isolate was passaged in Vero cells (ATCC, Manassas, VA, USA) four times, followed by once in C6/36 (*Aedes albopictus* larval cell line) followed again by Vero cells, and was then used for the experiments. A TCID50 assay was performed on Vero cells to determine viral titer.

### 2.2. *Aedes* Mosquito Rearing, Infection and RNA Extraction

All experiments were performed under biosafety level 3 (BSL-3) conditions in the insectary at the Australian Animal Health Laboratory, CSIRO. Insectary conditions were maintained at 27.5 °C and 70% in relative humidity with a 12 h light and dark cycle. Female mosquitoes (5–8 days old) were challenged with a chicken blood meal spiked with CHIKV (1 in 100 dilution of stock virus, TCID₅₀ 1.5 × 10⁹/mL) through chicken skin membrane feeding. After one hour, the mosquitoes were anesthetised with CO₂, and the blood-fed females sorted and kept in 200 mL cardboard cup containers at 27.5 °C, 70% humidity and a 14:10 day:night photoperiod for 2 days with 10% sugar solution ad libitum. For controls, females were fed with blood mixed with media supernatant from an uninfected Vero cell culture. Midguts were dissected at 2 dpi and were stored in 50 µL of Qiagen RLTplus buffer with 5–10 silica beads (1 mm) at –80 °C.

Bead beating was performed on MP Biomedicals FastPrep-24™ homogeniser, 3 cycles, speed: 6.5 m/s, 45 s each cycle. RNA was extracted using the RNeasy™ kit (Qiagen, Chadstone, Victoria, Australia), and cDNA was generated by using random hexamers and Superscript-III reverse transcriptase (Thermo Fisher Scientific Inc. Scoresby, Victoria, Australia) following the manufacturer's protocols.
Complementary DNA (cDNA) generated from the RNA extracted from the midguts pools was tested for CHIKV viral RNA using an in-house-designed qRT-PCR, using primers specific for the E1 gene (Table S1). For RNASeq data validation, adult *Aedes albopictus* female mosquitoes were infected with CHIKV as described above. RNA was extracted from the midguts of 5 infected mosquitoes 2 dpi, and cDNA was generated individually by previously described protocols. cDNA from the midguts of 5 uninfected mosquitoes was used as controls. Mosquitoes from multiple generations but approximately of the same age were used for the experiments.

For comparison of differential gene expression in *Aedes aegypti*, RNA was collected from the midguts of 6 mosquitoes and infected with CHIKV through blood feeding at 2 dpi. RNA that was extracted from the midguts of 5 uninfected mosquitoes was used as a control. cDNA was generated as described above and qPCR was performed with *Aedes aegypti* gene-specific primers.

## 2.3. qPCR

Quantitative PCR (qPCR) was performed using gene-specific primers and 18s rRNA specific primers as internal controls were used to validate the expression changes of 8 targets. Midgut tissue from three infected and control mosquitoes was used for validation of 5 long non-coding RNAs (lncRNAs) by the above-mentioned method.

qPCR was performed on an Applied Biosystems QuantStudio™ 6 using the SYBR Green Master Mix: SYBR Premix Ex Taq II (Tli RNase H Plus) (Takara-Scientifix Pty Ltd., Clayton, Victoria, Australia). The following cycling conditions were used with a melt curve at the end, 30 s at 95 °C, 40 cycles of 5 s at 95 °C and 30 s at 60 °C. The baseline and Ct values were calculated automatically using the supplied QuantStudio™ Software (Thermo Fisher Scientific Inc. Scoresby, Victoria, Australia), and the ΔΔCt values were calculated using the average ΔCt value of controls and 18s rRNA as a reference.

## 2.4. RNASeq and Viral Genome Sequencing

To identify the genes involved in the initial infection stage, midguts from CHIKV-infected *Aedes albopictus* mosquitoes at 2 dpi from 6 mosquitoes were pooled together for RNA extraction to obtain sufficient material. This was also required to avoid using a low-RNA input RNAseq kit, which would likely introduce bias during the PCR amplification stage. The pool size was kept as low as possible to retain information on biological variations.

Libraries for RNASeq were prepared using Nugen's Ovation Universal RNASeq kit, following manufacturer's specification with a minor modification in the HL-dsDNAse treatment. During first strand synthesis with DNase treatment, HL-dsDNase from Thermo Fisher Scientific was used in our library preparation, along with the 10x buffer supplied that their protocol used. The libraries were pooled and sequenced on a single lane of Hiseq-2500 (Macrogen Inc., Seoul, South Korea) to generate 2 × 100 bp reads. The fastq files were deposited in NCBI's (National Center for Biotechnology Information) Sequence Read Archive (SRA Accession ID: SRP140387).

The Qiagen QIAseq FX Single Cell RNA Library kit was used for Illumina library preparation from total RNA extracted from CHIKV-infected Vero cell culture supernatant using the RNeasy™ kit (Qiagen, Chadstone, Victoria, Australia). The library was sequenced on a Miniseq (Illumina, Scoresby, Victoria, Australia), with the mid-output kit (300 cycles) generating 2 × 150 bp paired-end reads. The resultant fastq files were quality-trimmed and assembled to a consensus sequence on the CLC Genomics workbench v9.5.2. The sequence was annotated and submitted to GenBank (MH229986).

## 2.5. Differential Gene Expression and Gene Ontology Analysis

Quality trimming of the raw sequences was performed using Trimmomatic v0.36. The reads were aligned to the CHIKV reference sequence (GenBank ID: MH229986) to assess the infection status using Hisat2 v2.0.5 [26].

Following removal of *Gallus gallus* (GenBank assembly accession: GCA_000002315.3) reads (due to chicken blood feeding) using SAMtools v1.3.1, the remaining reads were aligned to the *Aedes albopictus*
Foshan strain genome sequence (AaloF1) from Vectorbase using Hisat2, and the resultant SAM file was sorted and converted into a BAM file using SAMtools [27,28].

On Galaxy virtual lab v1.4.6.p5, featureCounts v1.4.6-p5 was used to quantify aligned transcripts from the sorted BAM files with default parameters for paired-end reads, and DESeq2 v2.11.38 was used to obtain differentially expressed genes between the controls and infected samples by using default parameters [29,30].

Using Trinity v2.3.2, a custom de novo transcriptome was built by combining the unaligned reads from midguts (D2) [31,32]. This transcriptome was used as a reference genome. edgeR was used to align the reads, measure the transcript counts and quantify the differentially expressed genes. The differentially expressed genes were annotated using BlastX [33,34].

The gene ontology (GO) IDs of the differentially expressed genes were obtained using the Biomart tool, and topGO analysis was performed to identify the Molecular Functions (MF), Biological Processes (BP) and Cellular Components (CC) that were either enriched or depleted in the differentially expressed genes [35,36]. Using topGO's classic algorithm, and based on *p*-values generated using Fisher's exact method, differentially expressed genes were grouped based on their ontologies. Enrichment percentage was calculated as the ratio of the number of times particular genes in the pathway were differentially expressed compared to the expected number by chance.

## 3. Results

### 3.1. Whole Genome Sequencing of CHIKV

The chikungunya virus isolate 06113879 was used to infect *Aedes albopictus* mosquitoes. Although a 559 bp portion of the E1 gene has been published (GenBank ID: EU404186.1), full genome sequencing of this isolate has not been performed, which can inform us about the genotype of this virus as well as be used as a reference for removing the virus reads that are detected during the RNAseq analysis.

Whole genome sequencing performed on the CHIKV isolate, after single Vero passaging, using the Illumina MiniSeq, resulted in about 8.5 million quality-trimmed, paired-end reads. The assembly resulted in an 11,929 bp long consensus sequence (MH229986), which perfectly matched the previously published 559 bp portion of the E1 gene from this isolate.

The viral consensus sequence was most similar (Identity: 11,705/11,985 (97.7%), Similarity: 11,705/11,985 (97.7%), Gaps: 245/11,985 (2.0%)) to the CHIKV strain LR2006_OPY1 (GenBank: KT449801.1) and had the E1-A226V mutation, indicating that this isolate also belongs to the same lineage as the CHIKV that caused the La Reunion outbreak in the Indian ocean in 2006.

### 3.2. RNAseq

To determine differentially expressed genes, five pools of midguts (six mosquitoes per pool) collected at 2 dpi from infected mosquitoes were created. For controls, three pools of midguts were used. Initially, to determine the infection status of these tissues, qRT-PCR was performed using CHIKV-specific primers. Based on these results and qualitative and quantitative requirements of the Nugen Ovation universal RNAseq kit, two control and three infected midgut pools were used to prepare sequencing libraries.

The sequencing of five libraries on a single lane of HiSeq-2500 resulted in 37 million to 170 million reads each. After quality trimming, reads mapping to the chicken genome were discarded to remove those originating from undigested chicken blood. The remaining reads were aligned to the *Aedes albopictus* reference genome with an average alignment of 62.25%. The read alignment to the CHIKV reference sequence (MH229986) also confirmed that all the three infected libraries contained viral reads, while the two control libraries did not, as shown in Table 1.
<table><caption>**Table 1.** RNASeq next-generation sequencing (NGS) data summary. The total number of obtained reads, the reads' alignment percentage to the reference genome and the percentage of reads aligned to the chikungunya virus genome.</caption><thead><tr><th></th><th>Total Reads (millions)</th><th>% Mapped to RefSeq Genome</th><th>% of CHIKV Reads</th></tr></thead><tbody><tr><td>Infected MG 1</td><td>170.65</td><td>60.83%</td><td>0.01%</td></tr><tr><td>Infected MG 2</td><td>37.13</td><td>60.46%</td><td>0.02%</td></tr><tr><td>Infected MG 3</td><td>38.02</td><td>65.28%</td><td>0.10%</td></tr><tr><td>Control MG1</td><td>44.77</td><td>61.53%</td><td>0.00%</td></tr><tr><td>Control MG2</td><td>44.72</td><td>63.15%</td><td>0.00%</td></tr></tbody></table>

### 3.3. Differential Expression and TopGO Analysis

Differentially expressed genes were identified using DESeq2 (for reads aligned to the RefSeq genome) and edgeR (for the de novo transcriptome built with aligned reads) and were plotted as Volcano plots (Figure 1). The results showed 25 genes to be differentially expressed as detected by DESeq2 (Up: 14 and Down: 11). edgeR did not show statistically significant (False Discovery Rate < 0.05) differential expression of genes due to a sub-optimal number of control libraries. However, targets from this dataset were chosen based on raw *p*-values for validation. The complete list of genes and transcripts that were differentially expressed with *p*-values of less than 0.05 is provided in Supplementary Information S2. The fasta files that were obtained as output for the custom transcriptome assembly are provided in Supplementary Information S3.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Volcano plots from DESeq2 and edgeR for 2 dpi differential gene expression; Volcano plots from DESeq2 (top panel) and edgeR (bottom panel) of differentially expressed genes from 2 dpi. DESeq2 was performed by aligning reads to the *Aedes albopictus* reference genome. edgeR analysis was done on reads that did not align to the reference genome and were aligned to the custom transcriptome.</figcaption></figure>
To determine the biological processes and molecular functions of the differentially expressed genes, gene set enrichment analysis and ontology was performed using topGO (Figure 2). As expected, several biological and molecular processes were significantly affected during CHIKV infection.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** topGo enrichment comparison in differentially expressed genes. Enrichment analysis of down- and upregulated genes in the midguts of *Aedes albopictus* mosquitoes in response to Chikungunya virus (CHIKV) at 2 dpi. Enrichment % is calculated as the ratio of ‘significant’ (Number of times the gene ontology number is observed as differentially expressed) to ‘expected’ (Number of times the gene ontology number is expected based on observation in control samples) gene numbers.</figcaption></figure>

The differentially regulated biological processes were metabolic (p-value = 0.012) and the serine family amino acid biosynthetic process (p-value = 0.026). The differentially regulated molecular functions, likely of interest from an innate immune response point of view, were lysozyme activity (downregulated, p-value = 0.028), alkaline phosphatase (downregulated, p-value = 0.033) and carboxypeptidase activity (upregulated, p-value = 0.042).

Apart from protein-coding genes, five long non-coding RNA (lnc RNA) were also detected to be differentially expressed. These were >200 bp long sequences with no apparent reading frames and, when searched using Blast against the published *Aedes albopictus* genomes on NCBI, showed high similarity to sequences annotated as long non-coding RNAs.

### 3.4. RNAseq Data Validation on qRT-PCR

Based on previously known immune involvement, we selected eight genes (one aligned to the reference genome and seven from the custom transcriptome) and five long non-coding RNA (from the custom transcriptome) for validation by qRT-PCR (Table 2). For this, midguts were dissected from five individual CHIKV-infected adult *Aedes albopictus* female mosquitoes (2 dpi). Midguts from uninfected mosquitoes were used as controls.

As the results showed, for *Aedes albopictus* mosquitoes, among the eight protein coding transcripts chosen for validation, the expression pattern for six targets was concordant, and the expression pattern for two targets was discordant, between RNAseq and qRT-PCR. Among the five tested long non-coding RNAs, two were concordant. All seven identified homologue targets, tested in *Aedes aegypti* midguts at 2 dpi, were upregulated as compared to uninfected controls. In *Aedes albopictus*, only five of these targets were upregulated as per qRT-PCR.
<table><caption>**Table 2.** List of genes validated by qRT-PCR and comparison of expression fold changes.</caption><thead><tr><th rowspan="2">Genes</th><th colspan="2">Aedes albopictus</th><th>Aedes aegypti</th></tr><tr><th>LogFC (RNAseq)</th><th>Expression Fold Change (qRT-PCR)</th><th>Expression Fold Change (qRT-PCR)</th></tr></thead><tbody><tr><td>ARF GTPase-activating protein GIT2</td><td>-8.7</td><td>0.76(↓)</td><td>39.74(↑)</td></tr><tr><td>NPC2 homolog</td><td>6.29</td><td>5.35(↑)</td><td>11.63(↑)</td></tr><tr><td>Mucin-22/FLO-11-like</td><td>-8.61</td><td>0.16(↓)</td><td>7.34(↑)</td></tr><tr><td>Translocon-associated protein subunit delta</td><td>8.82</td><td>1.88(↑)</td><td>40.18(↑)</td></tr><tr><td>ATP-dependent RNA helicase dbp2</td><td>-9.07</td><td>6.53(↑)</td><td>26.22(↑)</td></tr><tr><td>Uncharacterized gene coding for Sina and RING_Ubox domains containing protein</td><td>8.84</td><td>1.51(↑)</td><td>41.07(↑)</td></tr><tr><td>E3 ubiquitin-protein ligase MARCH6</td><td>-9.04</td><td>1.03(↑)</td><td>28.65(↑)</td></tr><tr><td>Ankyrin repeat domain-containing protein 44</td><td>-5.56</td><td>0.5(↓)</td><td>No homologue/Not tested</td></tr><tr><td>PREDICTED: Aedes albopictus uncharacterized LOC109424229 (LOC109424229), ncRNA</td><td>-11.22</td><td>2.14(↑)</td><td>No homologue/Not tested</td></tr><tr><td>PREDICTED: Aedes albopictus uncharacterized LOC109622934 (LOC109622934), ncRNA</td><td>-7.13</td><td>0.06(↓)</td><td>No homologue/Not tested</td></tr><tr><td>PREDICTED: Aedes albopictus uncharacterized LOC109423409 (LOC109423409), ncRNA</td><td>-9.07</td><td>0.73(↓)</td><td>No homologue/Not tested</td></tr><tr><td>PREDICTED: Aedes albopictus uncharacterized LOC109414360 (LOC109414360), ncRNA</td><td>-9.53</td><td>1.02(↑)</td><td>No homologue/Not tested</td></tr><tr><td>PREDICTED: Aedes albopictus uncharacterized LOC109424229 (LOC109424229), ncRNA</td><td>-8.98</td><td>0.04(↓)</td><td>No homologue/Not tested</td></tr></tbody></table>

## 4. Discussion

Chikungunya virus is a re-emerging alphavirus causing a high morbidity with long-term arthralgia in infected patients. Previous studies have taken approaches to understand the interaction between chikungunya virus and the *Aedes aegypti* vector, although these focused only on certain known genes and pathways [37]. However, considering the switch in vector preference towards *Aedes albopictus* by the Indian Ocean strains and invasive nature of this mosquito species, it is paramount to characterise the interaction between CHIKV and the new vector.

Previous studies with whole transcriptome analysis in mosquito vectors have used either whole mosquitoes or a cell culture [19,21–23,38]. Our objective here was to study the vector-virus interaction specifically at the midgut, which is the first barrier site to understand the factors that play a critical role in determining mosquito vector competence.

In the current study, unbiased transcriptional analysis was performed on midgut tissues collected from lab-reared adult female *Aedes albopictus* mosquitoes post CHIKV infection. The gene expression patterns, when compared to uninfected samples, revealed the transcriptional changes that are likely to be in response to the viral infection. Our analysis revealed that at 2 dpi in the midgut, most of the transcriptional changes were related to metabolism. An analysis of molecular functions revealed that while lysozyme activity and alkaline phosphatase were downregulated, carboxypeptidase activity was upregulated. Indeed, lysozymal and carboxypeptidase pathways have previously been implicated in innate immune responses [39–43].
One of the genes that was identified to be upregulated, and later validated by qRT-PCR, was a homolog of the Neimann Pick 2 (NPC2) gene, which in humans encodes for an intracellular cholesterol transporter. Loss of function mutations in the NPC2 gene lead to a lysosomal disorder known as Niemann–Pick disease type C, which causes an increased accumulation of lipids in cellular compartments and leads to cell death [44,45]. In many viral infections, replication and assembly of the viral particles occur at the lysosomal and other intracellular membrane-bound organelles, including DENV infection in its vector *Aedes aegypti* [46]. NPC2 has previously been shown to be involved in the replication and transport of several viruses, such as ebolavirus, vesicular stomatitis virus and influenza A virus, into and out of cells [47–49]. Interestingly, we detected a significant upregulation of the NPC2 homolog both in *Aedes albopictus* as well as *Aedes aegypti* midguts 2 days post infection with CHIKV. Recent studies have shown that Imipramine-based inhibition of the NPC2 protein results in severely diminished CHIKV replication in human fibroblast skin cells [50]. The fact that the NPC2 gene is significant and essential for CHIKV replication in humans and is significantly upregulated in the mosquito midgut after infection implies that an evolutionarily conserved mechanism may be at play and presents possible therapeutic opportunities in clinical treatment.

Multiple prior publications have also shown that the RNAi pathway is one of the major pathways involved in antiviral responses in insects [16,37,51,52]. We did not detect differential expression of the traditional innate-immune response pathways, such as JAK/STAT, IMD and Toll, consistent with previous studies with CHIKV infection of *Aedes aegypti* [37]. It is possible that the regulation of this pathway either does not occur at the transcriptional level or the proteins involved are ubiquitously expressed and not differentially regulated. It is also possible that the time point we selected did not coincide with RNAi activation.

Long non-coding RNAs are RNA molecules that are over 200 bp long and do not contain an open reading frame. They are produced by RNA Polymerase II and are processed like mRNAs by the cellular mechanism, including polyadenylation. lncRNAs have been implicated with involvement in multiple viral infections, including dengue and Influenza [39,53,54]. The mechanism of their involvement is complex, with pro or anti-viral activity [55–60]. Our understanding of lncRNA involvement in host viral interactions is still limited. Our results showed an increase in the number of lncRNAs in *Aedes albopictus* midguts after CHIKV challenge, suggesting a role during the infection. The functional significance of this differential regulation remains to be seen.

While the DESeq2 analysis identified statistically significant differential gene expression, edgeR, performed on the custom transcriptome that was assembled using reads that did not align to the published *Aedes albopictus* reference genome, was not able to find significantly altered genes with the set FDR < 0.05. This may be due to the small sample size.

To address this issue, we performed RT-qPCR-based validation of transcripts that appeared to be differentially expressed based on raw *p*-values, on a fresh set of individual mosquito midguts (five controls vs. five infected). RT-qPCR showed significant differential expression among all the transcripts tested, including in the targets from the edgeR dataset. These genes also showed differential expression when tested in CHIKV-infected *Aedes aegypti*, further validating their significance during infection.

The targets from the edgeR dataset, chosen for qPCR-based validation, were selected based on evidence from previously published data regarding their involvement in innate immune pathways. Seven out of eight protein coding targets and all of the five lncRNAs validated on RT-qPCR were from the edgeR dataset, and most showed varying levels of differential expression. Functional characterisation of the identified genes may help us to decipher the results and understand their role in mosquito-virus interactions.

Our results showed differences in gene expression pattern between *Aedes aegypti* and *Aedes albopictus*. While all of the seven tested genes in *Aedes aegypti* were upregulated, only five were upregulated in *Aedes albopictus* with the remaining two being downregulated. This may indicate differences in host-pathogen interactions between the two species of mosquitoes when infected with CHIKV.
Overall, our results showed significant changes in the transcriptome of *Aedes albopictus* mosquitoes after CHIKV infection, with the identified genes being involved in multiple cellular processes. This study examined differential gene expression at the midgut (the first critical barrier site) in infected *Aedes albopictus* mosquitoes. This study can be utilised in determining potential pro-viral and antiviral host factors and, in turn, will be helpful in reducing the high impact of CHIKV infections by targeting the vector *Aedes albopictus*.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/1999-4915/11/6/513/s1, Table S1: List of primers, S2: List of differentially expressed genes, S3: Custom transcriptome output.

**Author Contributions:** R.k.V., P.G., J.B.D., and P.N.P. conceived and designed the experiments. R.k.V., M.K., and J.B.D. performed the mosquito experiments. R.k.V. performed the RNAseq library prep. R.k.V. and M.N. performed the data analysis for RNAseq. R.k.V. and M.T. performed the whole genome sequencing of the virus. R.k.V. performed the qRT-PCR-based validations. R.k.V. and P.N.P performed the overall data analysis. R.k.V. wrote the manuscript and prepared the figures. All authors reviewed and approved the final manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** The chikungunya virus isolate was obtained from the Victorian Infectious disease reference laboratory, Melbourne. Thanks to Kim Blasdell for helping with the CLC Genomics workbench and to Chris Freebairn for collection of *Aedes albopictus* eggs from Torres Strait islands. Thanks to Omar Akbari of the University of California, San Diego for providing the *Aedes aegypti* Liverpool strain. Thanks to Adam Foord from CSIRO (AAHL) for discussions about optimisation of qRT-PCR and Dr Maria Doyle from the Peter MacCallum Cancer centre for guidance with the RNAseq data analysis.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. 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. [CrossRef]

2. Lumsden, W.H. An epidemic of virus disease in Southern Province, Tanganyika Territory, in 1952-53. II. General description and epidemiology. *Trans. R. Soc. Trop. Med. Hyg.* **1955**, *49*, 33-57. [CrossRef]

3. Robinson, M.C. An epidemic of virus disease in Southern Province, Tanganyika Territory, in 1952-53. I. Clinical features. *Trans. R. Soc. Trop. Med. Hyg.* **1955**, *49*, 28-32. [CrossRef]

4. Kamath, S.; Das, A.K.; Parikh, F.S. Chikungunya. *J. Assoc. Physicians India* **2006**, *54*, 725-726. [PubMed]

5. Higgs, S.; Vanlandingham, D. Chikungunya virus and its mosquito vectors. *Vector Borne Zoonotic Dis.* **2015**, *15*, 231-240. [CrossRef] [PubMed]

6. Simon, F.; Javelle, E.; Oliver, M.; Leparc-Goffart, I.; Marimoutou, C. Chikungunya Virus Infection. *Curr. Infect. Dis. Rep.* **2011**, *13*, 218-228. [CrossRef] [PubMed]

7. Reiter, P.; Fontenille, D.; Paupy, C. *Aedes albopictus* as an epidemic vector of chikungunya virus: Another emerging problem? *Lancet Infect. Dis.* **2006**, *6*, 463-464. [CrossRef]

8. Tsetsarkin, K.A.; Chen, R.; Weaver, S.C. Interspecies transmission and chikungunya virus emergence. *Curr. Opin. Virol.* **2016**, *16*, 143-150. [CrossRef]

9. Waldock, J.; Chandra, N.L.; Lelieveld, J.; Proestos, Y.; Michael, E.; Christophides, G.; Parham, P.E. The role of environmental variables on *Aedes albopictus* biology and chikungunya epidemiology. *Pathog. Glob. Health* **2013**, *107*, 224-241. [CrossRef]

10. Rezza, G.; Nicoletti, L.; Angelini, R.; Romi, R.; Finarelli, A.C.; Panning, M.; Cordioli, P.; Fortuna, C.; Boros, S.; Magurano, F.; et al. Infection with chikungunya virus in Italy: An outbreak in a temperate region. *Lancet* **2007**, *370*, 1840-1846. [CrossRef]

11. Enserink, M. Chikungunya: No Longer a Third World Disease. *Science* **2007**, *318*, 1860. [CrossRef] [PubMed]

12. Jupatanakul, N.; Sim, S.; Dimopoulos, G. The insect microbiome modulates vector competence for arboviruses. *Viruses* **2014**, *6*, 4294-4313. [CrossRef] [PubMed]

13. Arias-Goeta, C.; Mousson, L.; Rougeon, F.; Failloux, A.-B. Dissemination and Transmission of the E1-226V Variant of Chikungunya Virus in *Aedes albopictus* Are Controlled at the Midgut Barrier Level. *PLoS ONE* **2013**, *8*, e57548. [CrossRef] [PubMed]
14. Dong, S.; Kantor, A.M.; Lin, J.; Passarelli, A.L.; Clem, R.J.; Franz, A.W.E. Infection pattern and transmission potential of chikungunya virus in two New World laboratory-adapted *Aedes aegypti* strains. *Sci. Rep.* **2016**, *6*, 24729. [CrossRef] [PubMed]

15. Dubrulle, M.; Mousson, L.; Moutailler, S.; Vazeille, M.; Failloux, A.-B. Chikungunya Virus and Aedes Mosquitoes: Saliva Is Infectious as soon as Two Days after Oral Infection. *PLoS ONE* **2009**, *4*, e5895. [CrossRef] [PubMed]

16. Franz, A.W.E.; Kantor, A.M.; Passarelli, A.L.; Clem, R.J. Tissue Barriers to Arbovirus Infection in Mosquitoes. *Viruses* **2015**, *7*, 3741–3767. [CrossRef] [PubMed]

17. Kuno, G.; Chang, G.-J.J. Biological Transmission of Arboviruses: Reexamination of and New Insights into Components, Mechanisms, and Unique Traits as Well as Their Evolutionary Trends. *Clin. Microbiol. Rev.* **2005**, *18*, 608–637. [CrossRef] [PubMed]

18. Chompoosri, J.; Thavara, U.; Tawatsin, A.; Boonserm, R.; Phumee, A.; Sangkitporn, S.; Siriyasatien, P. Vertical transmission of Indian Ocean Lineage of chikungunya virus in *Aedes aegypti* and *Aedes albopictus* mosquitoes. *Parasites Vectors* **2016**, *9*, 227. [CrossRef] [PubMed]

19. 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. [CrossRef] [PubMed]

20. 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-00417. [CrossRef] [PubMed]

21. Shrinet, J.; Srivastava, P.; Sunil, S. Transcriptome analysis of *Aedes aegypti* in response to mono-infections and co-infections of dengue virus-2 and chikungunya virus. *Biochem. Biophys. Res. Commun.* **2017**, *492*, 617–623. [CrossRef] [PubMed]

22. Bonizzoni, M.; Dunn, W.A.; Campbell, C.L.; Olson, K.E.; Marinotti, O.; James, A.A. Complex Modulation of the *Aedes aegypti* Transcriptome in Response to Dengue Virus Infection. *PLoS ONE* **2012**, *7*, e50512. [CrossRef] [PubMed]

23. Fragkoudis, R.; Chi, Y.; Siu, R.W.; Barry, G.; Attarzadeh-Yazdi, G.; Merits, A.; Nash, A.A.; Fazakerley, J.K.; Kohl, A. Semliki Forest virus strongly reduces mosquito host defence signaling. *Insect Mol. Biol.* **2008**, *17*, 647–656. [CrossRef] [PubMed]

24. 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. [CrossRef] [PubMed]

25. Pyke, A.T. Partial E1 gene Sequence of Mauritius 2006 isolate 06113879. Unpublished work. 2008.

26. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. *Nat. Methods* **2015**, *12*, 357. [CrossRef]

27. Lawson, D.; Arensburger, P.; Atkinson, P.; Besansky, N.J.; Bruggner, R.V.; Butler, R.; Campbell, K.S.; Christophides, G.K.; Christley, S.; Dialynas, E.; et al. VectorBase: A data resource for invertebrate vector genomics. *Nucleic Acids Res.* **2009**, *37*, D583–D587. [CrossRef]

28. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. *Bioinformatics* **2009**, *25*, 2078–2079. [CrossRef]

29. Afgan, E.; Sloggett, C.; Goonasekera, N.; Makunin, I.; Benson, D.; Crowe, M.; Gladman, S.; Kowsar, Y.; Pheasant, M.; Horst, R.; et al. Genomics Virtual Laboratory: A Practical Bioinformatics Workbench for the Cloud. *PLoS ONE* **2015**, *10*, e0140829. [CrossRef]

30. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef]

31. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. *Nat. Biotechnol.* **2011**, *29*, 644. [CrossRef]

32. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-Seq: Reference generation and analysis with Trinity. *Nat. Protoc.* **2013**, *8*, 1494–1512. [CrossRef] [PubMed]

33. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. *BMC Bioinform.* **2009**, *10*, 421. [CrossRef] [PubMed]
34. Boratyn, G.M.; Schäffer, A.A.; Agarwala, R.; Altschul, S.F.; Lipman, D.J.; Madden, T.L. Domain enhanced lookup time accelerated BLAST. *Biol. Direct* **2012**, *7*, 12. [CrossRef] [PubMed]

35. Alexa, A.; Rahnenfuhrer, J. topGO: Enrichment Analysis for Gene Ontology, R package Version 2.33.0; 2016. Available online: https://bioconductor.org/packages/release/bioc/html/topGO.html (accessed on 6 May 2018).

36. Smedley, D.; Haider, S.; Durinck, S.; Pandini, L.; Provero, P.; Allen, J.; Arnaiz, O.; Awedh, M.H.; Baldock, R.; Barbiera, G.; et al. The BioMart community portal: An innovative alternative to large, centralized data repositories. *Nucleic Acids Res.* **2015**, *43*, W589–W598. [CrossRef] [PubMed]

37. 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 Negl. Trop. Dis.* **2014**, *8*, e2994. [CrossRef] [PubMed]

38. Colpitts, T.M.; Cox, J.; Vanlandingham, D.L.; Feitosa, F.M.; Cheng, G.; Kurscheid, S.; Wang, P.; Krishnan, M.N.; Higgs, S.; Fikrig, E. Alterations in the *Aedes aegypti* Transcriptome during Infection with West Nile, Dengue and Yellow Fever Viruses. *PLoS Pathog.* **2011**, *7*, e1002189. [CrossRef] [PubMed]

39. Hasan, M.; Koch, J.; Rakheja, D.; Pattnaik, A.K.; Brugarolas, J.; Dozmorov, I.; Levine, B.; Wakeland, E.K.; Lee-Kirsch, M.A.; Yan, N. Trex1 regulates lysosomal biogenesis and interferon-independent activation of antiviral genes. *Nat. Immunol.* **2013**, *14*, 61–71. [CrossRef]

40. Isoe, J.; Zamora, J.; Miesfeld, R.L. Molecular Analysis of the *Aedes aegypti* Carboxypeptidase Gene Family. *Insect Biochem. Mol. Biol.* **2009**, *39*, 68–73. [CrossRef]

41. Johnston, C.; Jiang, W.; Chu, T.; Levine, B. Identification of Genes Involved in the Host Response to Neurovirulent Alphavirus Infection. *J. Virol.* **2001**, *75*, 10431–10445. [CrossRef]

42. Merkling, S.H.; van Rij, R.P. Beyond RNAi: Antiviral defense strategies in Drosophila and mosquito. *J. Insect Physiol.* **2013**, *59*, 159–170. [CrossRef]

43. Swevers, L.; Liu, J.; Smagghe, G. Defense Mechanisms against Viral Infection in Drosophila: RNAi and Non-RNAi. *Viruses* **2018**, *10*, 230. [CrossRef] [PubMed]

44. Vance, J.E.; Karten, B. Niemann-Pick C disease and mobilization of lysosomal cholesterol by cyclodextrin. *J. Lipid Res.* **2014**, *55*, 1609–1621. [CrossRef] [PubMed]

45. Wang, Y.H.; Twu, Y.C.; Wang, C.K.; Lin, F.Z.; Lee, C.Y.; Liao, Y.J. Niemann-Pick Type C2 Protein Regulates Free Cholesterol Accumulation and Influences Hepatic Stellate Cell Proliferation and Mitochondrial Respiration Function. *Int. J. Mol. Sci.* **2018**, *19*, 1678. [CrossRef] [PubMed]

46. Jupatanakul, N.; Sim, S.; Dimopoulos, G. *Aedes aegypti* ML and Niemann-Pick type C family members are agonists of dengue virus infection. *Dev. Comp. Immunol.* **2014**, *43*, 1–9. [CrossRef] [PubMed]

47. Herbert, A.S.; Davidson, C.; Kuehne, A.I.; Bakken, R.; Braigen, S.Z.; Gunn, K.E.; Whelan, S.P.; Brummelkamp, T.R.; Twenhafel, N.A.; Chandran, K.; et al. Niemann-pick C1 is essential for ebolavirus replication and pathogenesis in vivo. *mBio* **2015**, *6*, e00565-00515. [CrossRef] [PubMed]

48. Infante, R.E.; Wang, M.L.; Radhakrishnan, A.; Kwon, H.J.; Brown, M.S.; Goldstein, J.L. NPC2 facilitates bidirectional transfer of cholesterol between NPC1 and lipid bilayers, a step in cholesterol egress from lysosomes. *Proc. Natl. Acad. Sci.* **2008**, *105*, 15287–15292. [CrossRef]

49. Amini-Bavil-Olyaee, S.; Choi, Y.J.; Lee, J.H.; Shi, M.; Huang, I.C.; Farzan, M.; Jung, J.U. The antiviral effector IFITM3 disrupts intracellular cholesterol homeostasis to block viral entry. *Cell Host Microbe* **2013**, *13*, 452–464. [CrossRef] [PubMed]

50. Wichit, S.; Hamel, R.; Bernard, E.; Talignani, L.; Diop, F.; Ferraris, P.; Liegeois, F.; Ekchariyawat, P.; Lupertlop, N.; Surasombatpattana, P.; et al. Imipramine Inhibits Chikungunya Virus Replication in Human Skin Fibroblasts through Interference with Intracellular Cholesterol Trafficking. *Sci. Rep.* **2017**, *7*, 3145. [CrossRef] [PubMed]

51. Cirimotich, C.M.; 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. [CrossRef] [PubMed]

52. Saleh, M.C.; Tassetto, M.; van Rij, R.P.; Goic, B.; Gausson, V.; Berry, B.; Jacquier, C.; Antoniewski, C.; Andino, R. Antiviral immunity in Drosophila requires systemic RNA interference spread. *Nature* **2009**, *458*, 346–350. [CrossRef] [PubMed]

53. Wang, J.; Wang, Y.; Zhou, R.; Zhao, J.; Zhang, Y.; Yi, D.; Li, Q.; Zhou, J.; Guo, F.; Liang, C.; et al. Host Long Noncoding RNA lncRNA-PAAN Regulates the Replication of Influenza A Virus. *Viruses* **2018**, *10*, 330. [CrossRef] [PubMed]
54. Etebari, K.; Asad, S.; Zhang, G.; Asgari, S. Identification of *Aedes aegypti* Long Intergenic Non-coding RNAs and Their Association with Wolbachia and Dengue Virus Infection. *PLoS Negl. Trop. Dis.* **2016**, *10*, e0005069.
[CrossRef] [PubMed]

55. Wang, X.J.; Jiang, S.C.; Wei, H.X.; Deng, S.Q.; He, C.; Peng, H.J. The Differential Expression and Possible Function of Long Noncoding RNAs in Liver Cells Infected by Dengue Virus. *Am. J. Trop. Med. Hyg.* **2017**, *97*, 1904–1912. [CrossRef] [PubMed]

56. Liu, W.; Ding, C. Roles of LncRNAs in Viral Infections. *Front. Cell. Infect. Microbiol.* **2017**, *7*, 205. [CrossRef] [PubMed]

57. Qiu, L.; Wang, T.; Tang, Q.; Li, G.; Wu, P.; Chen, K. Long Non-coding RNAs: Regulators of Viral Infection and the Interferon Antiviral Response. *Front. Microbiol.* **2018**, *9*, 1621. [CrossRef] [PubMed]

58. Du, M.; Yuan, L.; Tan, X.; Huang, D.; Wang, X.; Zheng, Z.; Mao, X.; Li, X.; Yang, L.; Huang, K.; et al. The LPS-inducible lncRNA Mirt2 is a negative regulator of inflammation. *Nat. Commun.* **2017**, *8*, 2049. [CrossRef] [PubMed]

59. Nikhil, S.; Singh, S. Implications of non-coding RNAs in viral infections. *Rev. Med. Virol.* **2016**, *26*, 356–368. [CrossRef]

60. Valadkhan, S.; Gunawardane, L.S. lncRNA-mediated regulation of the interferon response. *Virus Res.* **2016**, *212*, 127–136. [CrossRef]

<figure><img src="image_6.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).