# Whole Transcriptome Analysis of Aedes albopictus Mosquito Head and Thorax Post-Chikungunya Virus Infection

Ravi Kiran Vedururu, Matthew J. Neave, Vinod Sundaramoorthy, Diane Green, Jennifer A. Harper, Paul R. Gorry, Jean-Bernard Duchemin, Prasad Paradkar

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

To cite this version:

Ravi Kiran Vedururu, Matthew J. Neave, Vinod Sundaramoorthy, Diane Green, Jennifer A. Harper, et al.. Whole Transcriptome Analysis of Aedes albopictus Mosquito Head and Thorax Post-Chikungunya Virus Infection. Pathogens, 2019, 8 (3), pp.132. 10.3390/pathogens8030132. hal-03235171

HAL Id: hal-03235171

https://hal.science/hal-03235171

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_2.png" /></figure>

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

Article

# Whole Transcriptome Analysis of *Aedes albopictus*
Mosquito Head and Thorax Post-Chikungunya
Virus Infection

Ravi kiran Vedururu <fnref n="1" />,<fnref n="2" /> ID, Matthew J. Neave <fnref n="3" />, Vinod Sundaramoorthy <fnref n="3" /> ID, Diane Green <fnref n="3" />, Jennifer A. Harper <fnref n="3" />, Paul R. Gorry <fnref n="4" />, Jean-Bernard Duchemin <fnref n="1" />,† ID and Prasad N. Paradkar <fnref n="1" />,∗ ID

Received: 24 July 2019; Accepted: 23 August 2019; Published: 27 August 2019

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

check for
updates

**Abstract:** Chikungunya virus (CHIKV) is transmitted by *Aedes* mosquitoes and causes prolonged arthralgia in patients. After crossing the mosquito midgut barrier, the virus disseminates to tissues including the head and salivary glands. To better understand the interaction between *Aedes albopictus* and CHIKV, we performed RNAseq analysis on pools of mosquito heads and parts of the thorax 8 days post infection, which identified 159 differentially expressed transcripts in infected mosquitos compared to uninfected controls. After validation using RT-qPCR (reverse transcriptase-quantitative polymerase chain reaction), inhibitor of Bruton's tyrosine kinase (*BTKi*), which has previously been shown to be anti-inflammatory in mammals after viral infection, was further evaluated for its functional significance. Knockdown of *BTKi* using double-stranded RNA in a mosquito cell line showed no significant difference in viral RNA or infectivity titer. However, *BTKi* gene knocked-down cells showed increased apoptosis 24 hours post-infection compared with control cells, suggesting involvement of *BTKi* in the mosquito response to viral infection. Since *BTK* in mammals promotes an inflammatory response and has been shown to be involved in osteoclastogenesis, a hallmark of CHIKV pathogenesis, our results suggest a possible conserved mechanism at play between mosquitoes and mammals. Taken together, these results will add to our understanding of *Aedes Albopictus* interactions with CHIKV.

**Keywords:** Chikungunya; *Aedes albopictus*; RNAseq; host-pathogen interactions; Bruton's tyrosine kinase; *BTKi*

## 1. Introduction

Chikungunya virus (CHIKV) is an enveloped, positive-sense RNA virus belonging to the Alphavirus genus in the Togaviridae family [1,2]. The virus causes a disease characterised by a high febrile illness and debilitating and protracted arthralgia [3–5].

Since the 2004–2005 outbreak in Reunion island, east of Africa, *Aedes albopictus* has played a significant role in expansion of CHIKV outbreaks around the world. A single mutation in the E1 viral gene has been shown to confer increased vector competence to this mosquito [6]. Phylogenetic analysis of CHIKV genome sequences has identified three distinct lineages, namely West African, East/Central/south African (ECSA) and Asian. The Indian ocean sub-lineage (IOL), which caused the 2004 Reunion Island outbreak, is considered to have evolved from within the ECSA clade. An American sub-lineage was identified to have emerged from within the Asian clade

<footnote n="1">1. CSIRO Health & Biosecurity, Australian Animal Health Laboratory, Geelong 3220, Australia</footnote>

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

<footnote n="3">3. CSIRO, Australian Animal Health Laboratory, Geelong 3220, Australia</footnote>

<footnote n="4">4. School of Health and Biomedical Science, RMIT University, Bundoora 3083, Australia</footnote>

<footnote>* Correspondence: Prasad.Paradkar@csiro.au</footnote>

<footnote>† Current Address: Mosquitoes and Emerging Arboviruses, Institut Pasteur de la Guyane, Cayenne 97306, French Guiana.</footnote>
from an outbreak in the Caribbean and American states from around 2013. Virulence and pathogenicity of CHIKV varies among different lineages and clades with the West African lineage virus, which is considered to be more virulent than the Asian and American groups [7,8]. While most Alphaviruses are genetically similar, differences in 3' untranslated regions and in proteins including E2 and nSP3 are thought to contribute significantly to their individual virulence and pathogenesis characteristics [8–11]. Antisera raised against one lineage of CHIKV has also been found to be protective against other lineages of CHIKV, suggesting a high level of serological conservation [8,12].

In the mosquito after blood-feeding, the virus needs to cross two critical barrier tissues, the midgut and salivary glands. The virus must infect and cross the epithelial cells before digestion of the blood meal and then overcome the midgut escape barriers to disseminate into the haemocoel, from where it reaches and infects other mosquito tissues, including the salivary glands [13]. The mosquito is considered to be infective once the virus is detected in the saliva [14,15].

Previous research has identified a number of genes involved in the early stages of *Aedes albopictus* infection with CHIKV [16]. However, to understand the interactions between the virus and mosquito at the transcriptomic level after dissemination of the virus, we performed RNASeq analysis of the head and anterior one-third part of the thorax (containing the salivary glands), collected 8 days post CHIKV infection (dpi). The results identified multiple differentially expressed genes involved in various biological and molecular cellular functions. Our results identified *BTKi* to be functionally significant in the mosquito response to CHIKV infection.

## 2. Results

### 2.1. RNASeq and Differential Gene Expression (DGE) Analysis

To determine differentially expressed genes, two pools of head/thorax (six per pool) collected at 8 dpi from infected mosquitoes were created. For controls, one pool of head/thorax from uninfected mosquitoes was used. The three libraries yielded between 40.6 million and 47.15 million reads. After quality trimming and removing the reads mapping to the chicken genome (to remove read contributions by undigested chicken blood contamination), the remaining reads were aligned to the *Aedes albopictus* reference genome (Foshan strain genome sequence (AaloF1) from Vectorbase). The reads that did not map to the genome were then used to build a custom de novo transcriptome. The results confirmed the infection status of the two pools, and the control library contained no CHIKV reads as expected (reference sequence: MH229986) (Table 1).

<table><caption>**Table 1.** RNASeq data summary.</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 HT1</td><td>47.15</td><td>63.28%</td><td>1.12</td></tr><tr><td>Infected HT2</td><td>40.60</td><td>60.35%</td><td>2.09</td></tr><tr><td>Control HT</td><td>42.10</td><td>62.91%</td><td>0.00</td></tr></tbody></table>

Total reads obtained from HT (Head and Thorax), their alignment percentage to the reference *Aedes albopictus* genome and the percentage of reads aligned to the chikungunya virus genome.

Several differentially expressed genes were identified using both DESeq2 (for reads aligned to the RefSeq genome) and edgeR (for the de novo transcriptome built with unaligned reads) (Table 2), and the data was visualised as Volcano plots (Figure 1). The full list of genes and transcripts differentially expressed with *p*-values (and FDR (False Discovery Rate)) of less than 0.05 is provided in Supplementary File S2 (Total: 159, Up: 74 & Down: 85). The fasta file obtained as output for the custom transcriptome assembly is provided in Supplementary Information S3.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Volcano plots from DESeq2 and edgeR.</figcaption></figure>

Volcano plots from DESeq2 and edgeR of differentially expressed genes in CHIKV-infected *Aedes albopictus* head/thorax compared with uninfected controls. DESeq2 was performed by aligning reads to the reference *Aedes albopictus* genome; while edgeR analysis was done on reads that did not align to the reference genome and were aligned to the custom transcriptome.

<table><caption>**Table 2.** Differential Gene Expression analysis.</caption><thead><tr><th rowspan="2"></th><th colspan="3">No of Differentially Expressed Genes</th></tr><tr><th>Total</th><th>Up</th><th>Down</th></tr></thead><tbody><tr><td>DESeq2</td><td>96</td><td>51</td><td>45</td></tr><tr><td>edgeR</td><td>63</td><td>23</td><td>40</td></tr></tbody></table>

Number of genes found to be differentially expressed in mosquito head/thorax 8 dpi using DESeq2 (*p*-value: <0.05) and edgeR (FDR: <0.05) analysis.

## 2.2. Ontology Analysis

Gene set enrichment analysis and ontology performed using topGO revealed that several biological and molecular processes were significantly modified during CHIKV infection (Figure 2). The complete output of the topGO analysis is provided in Supplementary Information S4.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** topGO enrichment comparison in differentially expressed genes. Enrichment analysis of up- and down-regulated genes in the heads and thorax of *Aedes albopictus* mosquitoes in response to CHIKV at 8 dpi. Enrichment percentage is calculated as the ratio of "significant" (number of times the gene ontology number was observed as differentially expressed) to "expected" (number of times the gene ontology number was expected based on observation in control samples) gene numbers.</figcaption></figure>

In the head/thorax of *Aedes albopictus*, up-regulated molecular functions were RNA (*p*-value: 0.046) and mRNA binding (*p*-value: 0.049), and molecular functions that were down-regulated were lysozyme activity (*p*-value: 0.0038), serine-type endopeptidase activity (*p*-value: 0.0058), protein heterodimerisation (*p*-value: 0.00181) and chitin binding (*p*-value: 0.00532). Biological processes that were down-regulated were defence response (*p*-value: 9.70 × 10⁻⁵), proteolysis (*p*-value: 0.0075) and chitin metabolic process (*p*-value: 0.0106). Up-regulated biological processes were homophilic cell adhesion via the plasma membrane (*p*-value: 0.0056), UMP (uridine monophosphate) biosynthetic
process (*p*-value: 0.0245), regulation of transport (*p*-value: 0.0245) and spliceosomal complex assembly (*p*-value: 0.0342).

## 2.3. RT-qPCR (Reverse Transcriptase-Quantitative Polymerase Chain Reaction) Based Validation of RNASeq Data

Fourteen differentially expressed genes were selected, based on previously described involvement in the immune response, for validation using RT-qPCR (11 aligned to the reference genome and 3 from custom transcriptome) (Table 3). For this, heads/thoraxes were dissected from five individual CHIKV-infected adult *Aedes albopictus* female mosquitoes (8 dpi). Heads/thoraxes from uninfected mosquitoes were used as controls.

The RT-qPCR results, using gene-specific primers, compared with the RNAseq data (Table 3) showed that 7 out of 14 target genes were concordant.

<table><caption>**Table 3.** List of genes validated using qRT-PCR and a comparison of the expression fold change between RT-qPCR and RNAseq.</caption><thead><tr><th></th><th>Gene Annotation</th><th>Expression Fold Change (RT-qPCR)</th><th>LogFC (RNASeq)</th></tr></thead><tbody><tr><td>1</td><td>Quaking protein A</td><td>0.8(↓)</td><td>2.54</td></tr><tr><td>2</td><td>PDZ and LIM domain protein 7-like isoform X1</td><td>0.43(↓)</td><td>1.53</td></tr><tr><td>3</td><td>Putative ecdysteroid-regulated 16 kDa protein</td><td>0.96(↓)</td><td>2.07</td></tr><tr><td>4</td><td>Fat-like cadherin-related tumor suppressor homolog</td><td>0.01(↓)</td><td>3.51</td></tr><tr><td>5</td><td>Peptidylprolyl isomerase</td><td>0.75(↓)</td><td>2.24</td></tr><tr><td>6</td><td>Fasciclin-2-like isoform X1</td><td>2.97(↑)</td><td>3.55</td></tr><tr><td>7</td><td>Putative glycine-rich RNA binding protein</td><td>0</td><td>2.96</td></tr><tr><td>8</td><td>Leucine-rich immune protein (Long)</td><td>3.65(↑)</td><td>-1.65</td></tr><tr><td>9</td><td>Inhibitor of Bruton tyrosine kinase</td><td>17.54(↑)</td><td>2.49</td></tr><tr><td>10</td><td>Phosrestin i (arrestin b) (arrestin 2)</td><td>1149.59(↑)</td><td>1.51</td></tr><tr><td>11</td><td>Putative uncharacterised protein containing CysCysHisCys (CCHC) zinc finger domain</td><td>0.09(↓)</td><td>-3.12</td></tr><tr><td>12</td><td>Protein no-on-transient A isoform X2</td><td>0.01(↓)</td><td>-5.08</td></tr><tr><td>13</td><td>Ficolin-3-like</td><td>0.07(↓)</td><td>-6.42</td></tr><tr><td>14</td><td>PDZ and LIM domain protein Zasp-like isoform X5</td><td>2.9(↑)</td><td>10</td></tr></tbody></table>

List of genes selected for validation by qRT-PCR and their annotation based on BlastX and BlastN and results of RT-qPCR on the 14 gene targets selected from RNAseq DGE analysis. Expression fold change for RT-qPCR was calculated using 2⁻ΔΔCt and the up- or down-regulation compared to controls is indicated by arrows, while LogFC (from RNA-Seq data) was calculated using DESeq2 and edgeR.

## 2.4. Functional Significance of BTKi in CHIKV-Infected RML12 Cells

Inhibitor of Bruton's tyrosine kinase (BTKi) was identified to be one of the significantly up-regulated transcripts, both in RT-qPCR and RNAseq. To assess the functional significance of BTKi, an *Aedes albopictus* larval cell line (RML12) was transfected with anti-BTKi dsRNA. Twenty-four hours post transfection, the cells were infected with CHIKV (MOI: 1). Twenty-four hours post infection, the successful knockdown of BTKi gene was confirmed with RT-qPCR returning a ΔΔCt value of 2.39 (expression fold change: 0.1913, decrease in expression of BTKi: 80.87%). There were no significant changes in the viral RNA (ΔΔCt: -0.022, expression fold change: 1.015), or TCID₅₀
(Figure 3b,d). The results indicate that inhibition of *BTKi*, while detrimental to host cells, may not have a direct impact on viral replication.

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

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

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

<figure><img src="image_10.png" /><figcaption>**Figure 3.** (a) Comparison of BTKi expression between BTKi-KD and GFP (green fluorescent protein)-KD RML12 cells; (b) Comparison of intracellular CHIKV RNA via RT-qPCR in BTKi and GFP knockdown RML12 cells; (c) Comparison of % of apoptosis detected in BTKi-KD and GFP-KD RML12 cells; (d) Comparison of TCID<sub>50</sub> of extracellular CHIKV between BTKi and GFP knockdown RML12 cells.</figcaption></figure>

BTK has previously been shown to be involved in inducing apoptosis. To determine whether similar pathways are involved in CHIKV-infected mosquito cells, anti-*BTKi* dsRNA-transfected and CHIKV-infected RML12 cells were grown on glass cover slips and 24 hpi TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) staining was performed to quantify apoptosis. RML12 cells transfected with anti-GFP dsRNA and infected with CHIKV were used as controls. Cell nuclei were stained using DAPI and apoptosis was measured as the number of cells showing terminal deoxynucleotidyl transferase dUTP nick end labelling. The percentage of apoptosis was measured as the ratio of apoptotic cells to DAPI (4',6-diamidino-2-phenylindole) -stained cells (i.e., nuclei) (Figure 4). Significant apoptosis was detected in BTKi-KD cells compared to control cells, 3.07% in BTKi-KD RML12 cells versus 0.186% in control RML12 cells (t: 2.0965, df: 16 and *p*: 0.0261) (Figure 3c).
<figure><img src="image_11.png" /><figcaption>**Figure 4.** TUNEL staining of *BTKi* knocked down RML-12 cells 24 hours post infection with CHIKV showing increased apoptosis compared to control RML12 cells (Scale bar: 50 µM).</figcaption></figure>

**3. Discussion**

Multiple previous studies have been performed to understand the interactions between chikungunya virus and *Aedes aegypti*, the traditional vector [17–19]. Previously, RNASeq was performed on *Aedes albopictus* midguts 2 days post CHIKV infection, which identified a number of differentially expressed genes [16]. These results helped toward understanding the host-pathogen interaction at the first critical barrier site, the midgut. Yet, our understanding of transcriptome level interactions in *Aedes albopictus* after viral dissemination is preliminary. We attempted to study the virus-vector interaction in the head/thorax region 8 dpi as a proxy for the dissemination site. This site is important as successful infection of salivary glands (included in the thorax) ensures the presence of the virus in the saliva of the mosquito and in turn enables successful transmission.

Our results show differential regulation of biological processes including RNA and mRNA binding, lysosomal and serpin pathways and down-regulation of defensin genes. These could be due to either mosquito immune responses to CHIKV or viral modulation of immunity. Processes, such as the regulation of transport and homophilic cell adhesion via the plasma membrane, could be involved in viral assembly and export [20–23]. Functional studies to determine whether these genes are pro- or anti-viral are needed and could explain the role of these genes in the infection process.

Interestingly, two odorant binding proteins (OBPs) were found to be differentially expressed. While Obp25 (AALF018602) was up-regulated, D7 protein (AALF024478) was down-regulated. In *Aedes aegypti*, salivary glands infected with dengue virus showed differential expression of odorant binding proteins [24]. It was shown that this down-regulation of OBPs reduced the chemosensory abilities of the mosquitoes, which may affect host-seeking and hence there is a reduced exposure to the virus via feeding, which may hinder transmission capabilities. Our results suggest that similar mechanisms may be in play here as well, although this remains to be verified with functional studies.

The RNAi pathway is one of the major pathways involved in antiviral responses in insects [25,26]. Consistent with our previous study on *Aedes albopictus* midguts post CHIKV infection, our data does
not show statistically significant changes in the expression of genes involved in RNAi pathways in
the head and thorax. This is consistent with a previous study involving the chikungunya virus and
*Aedes aegypti* [17]. While it is possible that regulation of the RNAi pathway does not occur at the
transcriptional level or that RNAi proteins are not differentially regulated but ubiquitously expressed,
more focused validation may provide a clearer conclusion to this observation. Expression levels at
different time points post infection may also need to be assessed before reaching this conclusion.
The results are also consistent with other known antiviral genes like *JAK/STAT*, *IMD* and *Toll* pathways,
which were also not found to be differentially regulated [17].

Our results showed that the concordance between our RNASeq and RT-qPCR data was low (50%). This was possibly because there were not enough biological replicates for RNA-Seq analysis. Higher concordance can be achieved by increasing the number of infected and control samples. It is possible that due to the low number of replicates and reduced statistical power, some differentially expressed transcripts could be missed. Also, to be noted is the fact that tissue from six mosquitos were pooled for testing. There could be biological variation, both in terms of infection status and transcriptional response of the mosquitos. Hence, the RT-qPCR-based validation was performed on individual mosquitoes. This variation in pools could also explain the reduced concordance between RNASeq and RT-qPCR. For functional analysis, further research was focused on *BTKi*, which was validated through RT-qPCR.

BTKi is an inhibitor of Bruton's tyrosine kinase (Btk) [27]. Btk is heavily involved in innate immune responses [28]. In higher level organisms, inactivating mutations in Btk results in a condition known as X-linked agammaglobulinemia (XLA), characterised by the inability to produce mature B-cells and gamma globulins (including antibodies). Increased susceptibility to apoptotic death on exposure to pathogen-associated inflammatory signals is observed in BTK-deficient macrophages. BTK is involved in NK cell activation via the up-regulated NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) pathway [29]. Btk, along with Tec kinases, are heavily involved in osteoclast differentiation and regulation. Severe osteopetrosis is observed in mice with a non-functional *Btk* gene. Reduced osteoclastic bone resorption was observed in osteoporosis and inflammation-induced bone destruction on inhibition of Tec kinases. Btk and especially Tec kinases, are selectively expressed in osteoclasts and not in osteoblasts. Increased bone mass in mouse without Btk is due to defective osteoclastic bone resorption because of defective osteoclast differentiation [30]. CHIKV is known to infect osteoblasts and up-regulate IL-6 (interleukin-6) and RANKL (Receptor activator of nuclear factor kappa-B ligand) generation and decrease OPG (osteoprotegerin) production. An altered RANKL/OPG ratio gives rise to arthralgia, a characteristic morbidity of chikungunya fever [31]. Literature also presents evidence of the benefit of Btk inhibition in viral inflammatory disease. Specifically, BTK inhibition with ibrutinib has shown a major protective effect in the lung tissue of mice during influenza viral infection by ameliorating excessive inflammatory response [32].

Our data suggests that BTKi is significantly upregulated in the head and thorax tissue of *Aedes albopictus* mosquitoes after chikungunya virus infection. While the knock-down of *BTKi* in RML12 cells did not result in significant differences in viral titers or RNA levels, there was increased apoptosis in cells. These observations are consistent with a recently published study that showed that inhibition of BTK is protective in mammalian host tissue during influenza viral infection [32]. In the context of CHIKV-induced arthritis in patients, it is interesting to note that BTK is involved in osteoclastogenesis and the inhibition of BTK on bone resorption is protective. Our data suggests a possible conserved mechanism at play for BTK and associated Tec kinases in the context of chikungunya infections in both mosquitoes and mammals.

Overall, our results showed significant changes in the transcriptome of *Aedes albopictus* mosquitoes after CHIKV dissemination. This study, for the first time, examined differential gene expression in the head and thorax (a dissemination site containing salivary glands) of infected mosquitoes. Our results can be utilized for determining potential proviral and antiviral host factors, and in turn, will be helpful in reducing the high impact of CHIKV infections by targeting the vector, *Aedes albopictus*.
# 4. Materials and Methods

## 4.1. Chikungunya Virus and Aedes Mosquito Rearing, Infection, RNA Extraction and cDNA Preparation

*Aedes albopictus* mosquitoes were grown in the BSL-3 (Bio-Safety level-3) insectary at the Australian Animal Health Labs under the following conditions: 27.5 °C and 70% in relative humidity with a 12-h light and dark cycle from eggs obtained from Hammond Island in the Torres Straits. The mosquitoes used for experiments were from multiple generations (between 5–15 generations after collection). Chikungunya virus isolated in 2006 from a viremic patient, who returned from a trip to Mauritius, at the Victorian infectious disease reference laboratory, Melbourne, was used for this study. Our previous studies showed that this was a viral isolate from an ECSA clade (IOL sub-lineage), containing the E1-226 variation with improved fitness to *albopictus* mosquitoes and was closest in sequence similarity to the LR2006 isolate [16]. Five-to eight-day-old female mosquitoes were infected with CHIKV (1:100 dilution of stock; final TCID<sub>50</sub>: 1.5 × 10<sup>7</sup>/mL) spiked in fresh (not more than 2 days post bleeding) whole heparinized chicken blood collected from specific pathogen-free chicken collected on-site, through the Hemotek membrane feeding systems for one hour. Post feeding, blood-fed mosquitoes were separated and maintained in 200 mL cardboard containers at 27.5 °C, 70% humidity and in a 12:12 day:night photoperiod for 8 days with a 10% sugar solution ad libitum. For controls, chicken blood mixed with media supernatant from an uninfected Vero cell culture was fed to female mosquitoes. Eight days post exposure to CHIKV, the heads and anterior one-third of the thorax of the mosquitoes were dissected using 18-gauge needles, and tissues from six mosquitoes each were pooled into tubes, each with 50 µL of Qiagen RLTplus buffer (Qiagen, Chadstone, VIC, Australia) along with 5–10 silica beads (1 mm) and stored at –80 °C.

Tissue homogenisation was performed on a MP Biomedicals FastPrep -24™ homogeniser (MP Biomedicals Australia Pty Ltd, Seven Hills, NSW, Australia) (3 cycles, speed: 6.5 m/s, 45 s/cycle) via bead beating. RNA was then extracted using the RNeasy™ kit (Qiagen, Chadstone, VIC, Australia). Complementary DNA was generated by using random hexamers and Superscript-III reverse transcriptase (Thermo Fisher Scientific Inc., Scoresby, VIC, Australia) following the manufacturer's protocols. cDNA generated from the RNA was tested for CHIKV RNA using an in-house-designed qRT-PCR, using E1 gene specific primers (Supplementary Table S1). For RNASeq data validation, adult *Aedes albopictus* female mosquitoes were infected with CHIKV as described above. RNA was extracted from the HTs of five infected mosquitoes 8 dpi, and cDNA was generated individually using previously described protocols. cDNA from the head/thorax of five uninfected mosquitoes was used as controls.

## 4.2. RNASeq and RT-qPCR

RNASeq libraries were prepared using Nugen's Ovation Universal RNASeq kit (Tecan Genomics, Inc, Redwood City, CA, USA) following the manufacturer's specification with a minor modification in the HL-dsDNAse treatment. HL-dsDNase from Thermo Fisher Scientific was used in our library preparation, along with the 10× 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).

Gene-specific primers were designed for the 14 targets to be validated, and using 18s rRNA as an internal control, qPCR-based validation 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, VIC, Australia) (cycling conditions: 30 s at 95 °C, 40 cycles of 5 s at 95 °C and 30 s at 60 °C and a melt curve). Threshold cycle values were automatically calculated on the QuantStudio™ Software (Thermo Fisher Scientific Inc., Scoresby, VIC, Australia). ΔΔCt values were calculated using the average ΔCt value of controls and 18s rRNA as a reference. All primer sequences are provided in Supplementary Table S1.
### 4.3. Differential Gene Expression and Gene Ontology Analysis

Trimmomatic v0.36 was used for quality trimming of raw sequences [33]. Using Hisat2 v2.0.5, the reads were aligned to the CHIKV reference sequence (GenBank ID: MH229986) to assess the infection status [34]. The reads were aligned to the *Gallus gallus* genome and the aligned reads were removed bioinformatically before proceeding to the next step. This was performed to remove any RNA reads that could be from undigested chicken blood meal given to the mosquitos during infection. The reads were then aligned to the *Aedes albopictus* Foshan strain genome sequence (AaloF1) from Vectorbase (www vectorbase.org) using Hisat2 [35]. The sorted Bam files were uploaded into Galaxy virtual lab v1.4.6.p5 and aligned transcripts were quantified using featureCounts [36]. Using default parameters, DESeq2 v2.11.38 was used to obtain differentially expressed genes [37].

A custom de novo transcriptome was generated using Trinity v2.3.2, by combining the unaligned reads [38]. edgeR was used to align the reads, measure the transcript counts and quantify the differentially expressed genes using the custom transcriptome [39]. BlastX and BlastN were used to annotate the genes [40,41].

Using topGO and based on *p*-values (Fisher's exact method) differentially expressed genes were grouped based on their ontologies [42]. Enrichment percentage was calculated as the ratio of the number of times genes in the pathway were differentially expressed compared to the expected number by chance.

### 4.4. BTKi Knockdown and Functional Studies

Anti-BTKi dsRNA was generated by using a 422-bp-long segment (including the T7 promoter tag) of the coding sequence as a template via the Invitrogen™ MEGAscript™ RNAi Kit (Thermo Fisher Scientific Inc., Scoresby, VIC, Australia) by following the manufacturer's protocol (reaction incubation time: 4 hours) (primer details in Supplementary Table S1). RML12 cells were grown in a modified Leibovitz's L-15 medium (for 400 mL of L15 medium, 65 mL of FCS (fetal calf serum), 40 mL of TPB (Tryptose Phosphate Broth), 1.25 mL of Pen-strep and 5 mL of glutamine were added). Cells were maintained at 27 °C with 1% CO2.

RML12 cells of albopictus origin [43–45], grown in 24-well plates, were transfected with anti-BTKI dsRNA using Cellfectin®II (Thermo Fisher Scientific Inc., Scoresby, VIC, Australia) by following the manufacturer's protocol. Briefly, cells were incubated for 5 hours with a mixture of 800 ng of purified dsRNA, 8 µL of Cellfectin®II reagent and 240 µL of transfection media (L15 medium with TPB and glutamine but no FCS and antibiotics) per well. The mixture was replaced with normal L15 growth medium at the end of incubation. For controls, RML12 cells were transfected with anti-GFP dsRNA.

Twenty-four hours post transfection with dsRNA, the cells were infected with CHIKV (MOI: 1, stock TCID₅₀: 1.5 × 10⁷/mL). Briefly, each well (in a 24-well plate) of RML12 cells were incubated for one hour with 250 µL of infection medium (L15 medium with TPB and glutamine but no antibiotics and 2% FCS) and CHIKV. At the end of incubation, the infection media was replaced with regular L15 growth media.

Twenty-four hours post infection, the cell culture supernatant was collected and TCID₅₀ was performed on Vero cells. RNA was extracted by lysing the monolayer of RML12 cells with 350 µL of RLT Plus buffer. RNA extraction and cDNA generation were performed as described before. qPCR was performed to assess the knockdown of BTKi and change in viral RNA.

### 4.5. Confocal Microscopy and TUNEL Staining

For detection of apoptosis, RML12 cells were grown on glass cover slips (13 mm diameter) compatible with 24-well plates. The cells were transfected with anti-BTKI dsRNA or anti-GFP dsRNA (for controls) followed by CHIKV infection, as described before. BD Cytofix/Cytoperm™ (BD Biosciences, North Ryde, NSW, Australia) was used to fix the cells for one hour and stored in PBS at 4 °C.
TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labelling) staining was performed using an *in situ* cell death detection kit, TMR red (Sigma-Aldrich, Castle Hill, NSW, Australia) and the nuclei was stained with 4',6-Diamidine-2'-phenylindole dihydrochloride (DAPI) (Sigma-Aldrich, Castle Hill, NSW, Australia). The cells were then imaged in a Zeiss LSM 800 confocal system (Carl Zeiss Pty Ltd, North Ryde, NSW, Australia) using 10× objective covering at least 500 cells per image. Percentage of TUNEL positive apoptotic cells against the total number of cells (DAPI stained nuclei) was quantified using ImageJ (particle analyser plugin, https://imagej.nih.gov/ij/index.html) from each image.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2076-0817/8/3/132/s1, S1: List of primers, S2: List of differentially expressed genes, S3: Custom transcriptome, S4 topGO output.

**Author Contributions:** R.k.V., P.R.G., J.-B.D. and P.N.P. conceived and designed the experiments. R.k.V. 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. performed the cell culture and functional studies. V.S., J.H. and D.G. performed staining and microscopy. 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 Chris Freebairn for the collection of *Aedes albopictus* eggs from Torres Strait Islands. Thanks to Adam Foord from CSIRO (AAHL) for discussions about the optimisation of qRT-PCR and Maria Doyle from the Peter MacCallum Cancer Centre for guidance with the RNAseq data analysis. We thank the Australian Microscopy & Microanalysis Research Facility (AMMRF) for supporting the confocal imaging capability at AAHL that was utilised in this study.

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

## References

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

2. 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]

3. Kamath, S.; Das, A.K.; Parikh, F.S. Chikungunya. *J. Assoc. Phys. India* **2006**, *54*, 725–726.

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

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

6. Tsetsarkin, K.A.; Vanlandingham, D.L.; McGee, C.E.; Higgs, S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. *PLoS Pathog.* **2007**, *3*, e201. [CrossRef] [PubMed]

7. Lanciotti, R.S.; Lambert, A.J. Phylogenetic Analysis of Chikungunya Virus Strains Circulating in the Western Hemisphere. *Am. J. Trop. Med. Hyg.* **2016**, *94*, 800–803. [CrossRef]

8. Langsjoen, R.M.; Haller, S.L.; Roy, C.J.; Vinet-Oliphant, H.; Bergren, N.A.; Erasmus, J.H.; Livengood, J.A.; Powell, T.D.; Weaver, S.C.; Rossi, S.L. Chikungunya Virus Strains Show Lineage-Specific Variations in Virulence and Cross-Protective Ability in Murine and Nonhuman Primate Models. *mBio* **2018**, *9*, e02449-17. [CrossRef]

9. Teo, T.-H.; Her, Z.; Tan, J.J.L.; Lum, F.-M.; Lee, W.W.L.; Chan, Y.-H.; Ong, R.-Y.; Kam, Y.-W.; Leparc-Goffart, I.; Gallian, P.; et al. Caribbean and La Réunion Chikungunya Virus Isolates Differ in Their Capacity To Induce Proinflammatory Th1 and NK Cell Responses and Acute Joint Pathology. *J. Virol.* **2015**, *89*, 7955–7969. [CrossRef]

10. Chen, R.; Wang, E.; Tsetsarkin, K.A.; Weaver, S.C. Chikungunya virus 3' untranslated region: Adaptation to mosquitoes and a population bottleneck as major evolutionary forces. *PLoS Pathog.* **2013**, *9*, e1003591. [CrossRef]

11. Arankalle, V.A.; Shrivastava, S.; Cherian, S.; Gunjikar, R.S.; Walimbe, A.M.; Jadhav, S.M.; Sudeep, A.B.; Mishra, A.C. Genetic divergence of Chikungunya viruses in India (1963–2006) with special reference to the 2005–2006 explosive epidemic. *J. Gen. Virol.* **2007**, *88*, 1967–1976. [CrossRef]
12. Roy, C.J.; Adams, A.P.; Wang, E.; Plante, K.; Gorchakov, R.; Seymour, R.L.; Vinet-Oliphant, H.; Weaver, S.C. Chikungunya vaccine candidate is highly attenuated and protects nonhuman primates against telemetrically monitored disease following a single dose. *J. Infect. Dis.* **2014**, *209*, 1891–1899. [CrossRef]

13. 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]

14. 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]

15. 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]

16. Vedururu, R.K.; Neave, M.J.; Tachedjian, M.; Klein, M.J.; Gorry, P.R.; Duchemin, J.-B.; Paradkar, P.N. RNAseq Analysis of Aedes albopictus Mosquito Midguts after Chikungunya Virus Infection. *Viruses* **2019**, *11*, 513. [CrossRef]

17. 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]

18. 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]

19. 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]

20. Mateo, M.; Generous, A.; Sinn, P.L.; Cattaneo, R. Connections matter—How viruses use cell–cell adhesion components. *J. Cell Sci.* **2015**, *128*, 431–439. [CrossRef]

21. Bhella, D. The role of cellular adhesion molecules in virus attachment and entry. *Philos. Trans. R. Soc. B Biol. Sci.* **2015**, *370*, 20140035. [CrossRef]

22. Chisenhall, D.M.; Christofferson, R.C.; McCracken, M.K.; Johnson, A.-M.F.; Londono-Renteria, B.; Mores, C.N. Infection with dengue-2 virus alters proteins in naturally expectorated saliva of Aedes aegypti mosquitoes. *Parasites Vectors* **2014**, *7*, 252. [CrossRef]

23. Gulley, M.M.; Zhang, X.; Michel, K. The roles of serpins in mosquito immunology and physiology. *J. Insect Physiol.* **2013**, *59*, 138–147. [CrossRef]

24. Sim, S.; Ramirez, J.L.; Dimopoulos, G. Dengue Virus Infection of the Aedes aegypti Salivary Gland and Chemosensory Apparatus Induces Genes that Modulate Infection and Blood-Feeding Behavior. *PLoS Pathog.* **2012**, *8*, e1002631. [CrossRef]

25. Campbell, C.L.; Keene, K.M.; Brackney, D.E.; Olson, K.E.; Blair, C.D.; Wilusz, J.; Foy, B.D. Aedes aegyptiuses RNA interference in defense against Sindbis virus infection. *BMC Microbiol.* **2008**, *8*, 47. [CrossRef]

26. 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]

27. Liu, W.; Quinto, I.; Chen, X.; Palmieri, C.; Rabin, R.L.; Schwartz, O.M.; Nelson, D.L.; Scala, G. Direct inhibition of Brutons tyrosine kinase by IBtk, a Btk-binding protein. *Nat. Immunol.* **2001**, *2*, 939. [CrossRef]

28. Lopez-Herrera, G.; Vargas-Hernandez, A.; Gonzalez-Serrano, M.E.; Berron-Ruiz, L.; Rodriguez-Alba, J.C.; Espinosa-Rosales, F.; Santos-Argumedo, L. Bruton's tyrosine kinase–an integral protein of B cell development that also has an essential role in the innate immune system. *J. Leukoc. Biol.* **2014**, *95*, 243–250. [CrossRef]

29. Bao, Y.; Zheng, J.; Han, C.; Jin, J.; Han, H.; Liu, Y.; Lau, Y.L.; Tu, W.; Cao, X. Tyrosine kinase Btk is required for NK cell activation. *J. Biol. Chem.* **2012**, *287*, 23769–23778. [CrossRef]

30. Shinohara, M.; Koga, T.; Okamoto, K.; Sakaguchi, S.; Arai, K.; Yasuda, H.; Takai, T.; Kodama, T.; Morio, T.; Geha, R.S.; et al. Tyrosine kinases Btk and Tec regulate osteoclast differentiation by linking RANK and ITAM signals. *Cell* **2008**, *132*, 794–806. [CrossRef]

31. Amdekar, S.; Parashar, D.; Alagarasu, K. Chikungunya Virus-Induced Arthritis: Role of Host and Viral Factors in the Pathogenesis. *Viral Immunol.* **2017**, *30*, 691–702. [CrossRef] [PubMed]
32. Florence, J.M.; Krupa, A.; Booshehri, L.M.; Davis, S.A.; Matthay, M.A.; Kurdowska, A.K. Inhibiting Bruton's tyrosine kinase rescues mice from lethal influenza-induced acute lung injury. *Am. J. Physiol. Lung Cell. Mol. Physiol.* **2018**, *315*, L52–L58. [CrossRef]

33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics* **2014**, *30*, 2114–2120. [CrossRef] [PubMed]

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

35. Giraldo-Calderon, 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.* **2015**, *43*, D707–D713. [CrossRef] [PubMed]

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

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

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

39. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* (Oxford, England) **2009**, *26*, 139–140. [CrossRef]

40. 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]

41. 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]

42. Adrian Alexa, J.R. topGO: Enrichment Analysis for Gene Ontology. *R package version* 2.32.0. 2016. Available online: https://bioconductor.org/packages/release/bioc/html/topGO.html (accessed on 6 May 2018).

43. Voronin, D.; Tran-Van, V.; Potier, P.; Mavingui, P. Transinfection and growth discrepancy of Drosophila Wolbachia strain wMel in cell lines of the mosquito Aedes albopictus. *J. Appl. Microbiol.* **2010**, *108*, 2133–2141. [CrossRef] [PubMed]

44. Frentiu, F.D.; Robinson, J.; Young, P.R.; McGraw, E.A.; O'Neill, S.L. Wolbachia-Mediated Resistance to Dengue Virus Infection and Death at the Cellular Level. *PLoS ONE* **2010**, *5*, e13398. [CrossRef] [PubMed]

45. McLean, B.J.; Hobson-Peters, J.; Webb, C.E.; Watterson, D.; Prow, N.A.; Nguyen, H.D.; Hall-Mendelin, S.; Warrilow, D.; Johansen, C.A.; Jansen, C.C.; et al. A novel insect-specific flavivirus replicates only in Aedes-derived cells and persists at high prevalence in wild Aedes vigilax populations in Sydney, Australia. *Virology* **2015**, *486*, 272–283. [CrossRef] [PubMed]

<figure><img src="image_12.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/).