# Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis

Cassandra Modahl, Avisha Chowdhury, Dolyce Low, Menchie Manuel, Dorothée Missé, R Manjunatha Kini, Ian Mendenhall, Julien Pompon

► To cite this version:

Cassandra Modahl, Avisha Chowdhury, Dolyce Low, Menchie Manuel, Dorothée Missé, et al., Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis. Scientific Reports, 2023, 13 (1), pp.11271. 10.1038/s41598-023-38354-9. hal-04911291

HAL Id: hal-04911291

https://hal.science/hal-04911291v1

Submitted on 24 Jan 2025

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>

Check for updates

OPEN

# Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors *Aedes albopictus* and *Aedes malayensis*

Cassandra M. Modahl<sup>1,4,7</sup>, Avisha Chowdhury<sup>1,5,7</sup>, Dolyce H. W. Low<sup>2</sup>, Menchie C. Manuel<sup>2</sup>, Dorothée Missé<sup>6</sup>, R. Manjunatha Kini<sup>1,3</sup>, Ian H. Mendenhall<sup>2</sup> & Julien Pompon<sup>2,6</sup>

Dengue (DENV) and chikungunya (CHIKV) viruses are among the most preponderant arboviruses. Although primarily transmitted through the bite of *Aedes aegypti* mosquitoes, *Aedes albopictus* and *Aedes malayensis* are competent vectors and have an impact on arbovirus epidemiology. Here, to fill the gap in our understanding of the molecular interactions between secondary vectors and arboviruses, we used transcriptomics to profile the whole-genome responses of *A. albopictus* to CHIKV and of *A. malayensis* to CHIKV and DENV at 1 and 4 days post-infection (dpi) in midguts. In *A. albopictus*, 1793 and 339 genes were significantly regulated by CHIKV at 1 and 4 dpi, respectively. In *A. malayensis*, 943 and 222 genes upon CHIKV infection, and 74 and 69 genes upon DENV infection were significantly regulated at 1 and 4 dpi, respectively. We reported 81 genes that were consistently differentially regulated in all the CHIKV-infected conditions, identifying a CHIKV-induced signature. We identified expressed immune genes in both mosquito species, using a de novo assembled midgut transcriptome for *A. malayensis*, and described the immune architectures. We found the JNK pathway activated in all conditions, generalizing its antiviral function to Aedines. Our comprehensive study provides insight into arbovirus transmission by multiple *Aedes* vectors.

Arboviruses (arthropod-borne viruses) are a global public health burden<sup>1,2</sup>. Dengue viruses (DENV) from the *Flaviviridae* family and chikungunya virus (CHIKV) from the *Togaviridae* family are among the most prevalent arboviruses. Both viruses are now endemic to subtropical and tropical areas of Asia, the Americas, Oceania and Africa, and cause sporadic epidemics in Europe and North America<sup>3,4</sup>. Climate change modeling forecasts an expansion of arbovirus geographic distribution into temperate regions, which have previously been relatively spared<sup>5,6</sup>. Currently, DENV infects an estimated 400 million people yearly, causing symptoms ranging from flu-like to severe hemorrhage, leading to shock and death for approximately 25,000 people<sup>7</sup>. CHIKV emerged as a global pathogen in 2004 and has ever since infected more than 6 million people, causing flu-like symptoms associated with severe polyarthralgia, in some cases lasting for years and drastically impairing life quality<sup>8</sup>. Unfortunately, there is no effective curative treatment for these arbovirus infections, only a sub-effective vaccine with associated safety risks for DENV<sup>9,10</sup>, and no licensed vaccine against CHIKV. Vector control remains the primary intervention to limit arbovirus spread and is widely deployed. However, current vector control strategies have modest efficacy in preventing epidemics and limited success against container-ovipositing species<sup>11-13</sup>. A better understanding of the molecular mechanisms that enable viral transmission by mosquito vectors will help to design novel interventions.

Both DENV and CHIKV are transmitted by *Aedes* mosquitoes, mainly by the species *Aedes aegypti*<sup>7,8</sup>. However, multiple *Aedes* species have long been known to be susceptible to DENV and CHIKV infections<sup>14</sup> and can act as vectors of importance, as exemplified by the mild dengue and chikungunya epidemics that occurred in areas without *A. aegypti*<sup>8,15-19</sup>. For instance, *Aedes albopictus* is sufficient to maintain endemic dengue in La

<sup>1</sup>Department of Biological Sciences, National University of Singapore, Singapore, Singapore. <sup>2</sup>Programme in Emerging Infectious Diseases, Duke-NUS Medical School, Singapore, Singapore. <sup>3</sup>Present address: Department of Pharmacology, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore. <sup>4</sup>Present address: Liverpool School of Tropical Medicine, Liverpool, U.K.. <sup>5</sup>Present address: Toronto Centre for Liver Disease, Toronto General Hospital, University Health Network, University of Toronto, Toronto, Canada. <sup>6</sup>Present address: MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France. <sup>7</sup>These authors contributed equally: Cassandra M. Modahl and Avisha Chowdhury. <sup>8</sup>email: Julien.pompon@ird.fr
Réunion island<sup>20,21</sup>, is responsible for sporadic outbreaks of increasing intensity in Southern Europe<sup>22–24</sup>, and caused an outbreak in an urban park in Japan<sup>18</sup>. For CHIKV, its envelope protein evolved to increase transmission by *A. albopictus* without compromising *A. aegypti* vector competence<sup>25</sup>, causing an unprecedented epidemic spanning the Indian Ocean region, East Africa and India<sup>3</sup>. Beside *A. albopictus*, our group identified *Aedes malayensis*, a species closely related to *A. albopictus* and distributed throughout Southeast Asia<sup>26</sup>, as a competent vector for both DENV and CHIKV<sup>27</sup>. Vector capacity of *A. albopictus* and *A. malayensis* is further supported by their anthropophilic feeding habits<sup>28,29</sup>. While evidence for the vector status of these two peridomestic *Aedes* species is mounting, there is a dearth of knowledge on the molecular mechanisms that regulate DENV and CHIKV infection in these secondary vectors.

Arbovirus transmission occurs when viruses imbibed during a mosquito bite on an infectious host propagate from the midgut to the salivary glands, from where they are salivated into the skin of other hosts during subsequent bites³⁰. The initial midgut infection (the first barrier) that determines transmission success is strongly controlled by immunity³⁰,³¹. Mosquitoes mount a potent immune response to infection through the canonical immune pathways, namely, Toll (TOLLPATH), Immune-Deficiency (IMDPATH) and Janus Kinase/Signal Transduction and Activators of Transcription (JAKSTAT)³²,³³. Recently, we added the c-Jun *N*-terminal Kinase (JNK) pathway as another immune arm with broad antiviral activity against DENV, Zika virus and CHIKV³⁴. The multiple immune pathways share a common architecture for activation³⁵. Pathogens are detected by pathogen recognition receptors (PRR), mainly gram-negative binding (GNBP) and peptidoglycan receptor-like proteins (PGRP). The activation signal is transferred through Spaetzle proteins (SPZ) or directly to transmembrane receptors such as Toll and IMD receptors, which trigger a cytoplasmic biochemical cascade resulting in the activation and nuclear translocation of transcription factors like Rel (REL) for TOLLPATH and IMDPATH, STAT for JAK-STAT or Kayak for JNK. The signaling by PRR and through the biochemical cascade are regulated by proteases such as CLIP-domain serine proteases (CLIP) and serine protease inhibitors (SRPN). Translocated transcription factors then upregulate antiviral effectors including anti-microbial peptides (AMP), Thioester-containing proteins (TEP), C-type lectins (CTL), fibrinogen-related proteins (FREP), galectins (GALE), lysozymes (LYS), MD2-like proteins (ML), scavenger receptors (SCR), catalases (CAT), superoxide dismutases (SOD) and peroxidases (HPX), the last three acting by regulating reactive oxygen species (ROS). Moreover, RNAi limits arbovirus infections by degrading viral RNA genomes through siRNA regulator proteins (SRRP)³⁶,³⁷. Viral infection can also be restricted by autophagy (APHAG), by apoptosis mediated by caspases (CASP), caspase activators (CASPA) and inhibitors of apoptosis (IAP), and by melanization triggered by prophenoloxidases (PPO). The resulting multimodal immunity restricts DENV and CHIKV as early as 3 days post-infection (dpi)³⁸⁻⁴⁴.

Most studies on mosquito antiviral immunity have been conducted with *A. aegypti*. Our knowledge regarding the immune response in other *Aedes* vectors is strikingly deficient. Only a few transcriptomic studies have been conducted with *A. albopictus*, and, to our knowledge, none has been performed with *A. malayensis* or any other *Aedes* vector. In *A. albopictus*, a previous study showed that DENV infection induced a low number of differentially expressed genes (DEGs) at 1 and 5 dpi in midguts, but only one DEG at 1 dpi was related to immunity⁴⁵. Another study with DENV and *A. albopictus* identified 18 upregulated DEGs with a majority of AMPs, and 13 downregulated DEGs including immune pathway components such as *Myd88* and *STAT* at 1 dpi⁴⁶. For CHIKV, only one study looked at transcriptomic regulation in *A. albopictus*. The results show 25 DEGs primarily related to metabolism and none to immunity in midguts at 2 dpi, although these regulations could not be confirmed by RT-qPCR⁴⁷. Despite the growing interest in *A. albopictus* and *A. malayensis*, there is a clear lack of understanding of their basic molecular responses to DENV and CHIKV.

Here, to bridge the knowledge gap in our understanding of the molecular response and immunity in secondary vectors, we provide a comprehensive account of the transcriptomic responses at 1 and 4 dpi with CHIKV in both *A. albopictus* and *A. malayensis*, and with DENV in *A. malayensis*. We deployed high-throughput sequencing and validate gene expressions by RT-qPCR in different mosquito batches. We describe infection-induced gene regulations using the chromosome-level *A. albopictus* genome assembly and a de novo assembled *A. malayensis* midgut transcriptome, as this species' genome is not available. By providing insight regarding the transcriptomic responses to two arboviruses in two relevant *Aedes* vectors, our study expands our understanding of arbovirus transmission and identifies common targets to block mosquito transmission of arboviruses.

## Results

**Transcriptome regulations by CHIKV and DENV in *A. albopictus* and *A. malayensis* midguts at 1 and 4 dpi.** To determine the transcriptomic response to arboviral infections in midguts of secondary vectors, *A. albopictus* mosquitoes were orally infected with 10⁷ pfu/ml of CHIKV and *A. malayensis* mosquitoes were orally infected with either the same concentration of CHIKV or 2 × 10⁷ pfu/ml of DENV. Both blood inoculum are within the high-end of viremia measured in patients⁴⁸,⁴⁹ and resulted in 100% infected mosquitoes³⁴. Controls were fed non-infectious blood. Viruses and mosquitoes were sympatric as they were all collected in Singapore, and the CHIKV strain collected before 2010 did not possess the envelope mutation that enhances *A. albopictus* infection⁵⁰. As *Wolbachia* can influence mosquito susceptibility to arboviruses⁵¹, we conducted a PCR detection assay and did not find the bacteria in our colonies (Fig. S1).

We then performed high-throughput RNA-sequencing on three replicates of 20 pooled midguts collected at 1 and 4 dpi for each virus and mosquito combination. We selected these two time points to cover the initial response to infection at 1 dpi and the established infection at 4 dpi. We obtained between 38 and 52 million reads per sample (Table S1). After quality filtering, we confirmed mosquito species and virus infection by mapping reads to species-specific cytochrome c oxidase (COI) genes and each viral genome, respectively (Table S1). We also quantified the percentage of viral reads and observed that CHIKV reads were more abundant than DENV and had plateaued at 1 dpi in both mosquito species, whereas DENV reads increased between 1 and 4
dpi (Fig. 1a). Since CHIKV and DENV have different patterns of infection kinetics⁵,⁵³, this should be taken into account when interpreting gene regulations.

Using the most recent genome assembles for *A. albopictus*⁵⁴ (AalbFP1.0) and *A. aegypti*⁵⁵ (AaegL5.0) (*A. malayensis* genome has not been assembled), we mapped reads from *A. albopictus* and *A. malayensis* and calculated DEGs between infected samples and time-matched non-infected controls. To select the best genome approach and validate DEGs, we quantified expression levels of several genes at 4 dpi by RT-qPCR in midguts of separate mosquito batches for both species. For both *A. albopictus* and *A. malayensis*, DEGs obtained by mapping on the *A. albopictus* genome correlated best with qPCR data (Fig. S2 and S3) and were further analyzed. Of note, *A. malayensis* genes do not have IDs and hereafter, are identified from *A. albopictus* orthologs.

To identify similarities in the transcriptomes between our different conditions, we performed a principal component analysis (PCA) on normalized gene expression values (Fig. 1b). Interestingly, virus infection was not the main source of transcriptome variation, mosquito species and time post blood-feeding were responsible for the greatest variance. The PC1 explained 66% of the total variance and segregated the mosquito species, and the PC2 explained 20% of the total variance and separated the time of collection.

Upon CHIKV infection in *A. albopictus*, there were 1,793 DEGs at 1 dpi and 339 at 4 dpi (Table S2). Upon CHIKV infection in *A. malayensis*, there were 943 DEGs at 1 dpi and 222 at 4 dpi. Upon DENV infection in *A. malayensis*, there were 74 DEGs at 1 dpi and 69 at 4 dpi.

Among all the CHIKV infection conditions (i.e., both mosquito species at 1 and 4 dpi), we identified shared DEGs (Fig. 1c). First, we found 81 DEGs common to all conditions and that 79 were upregulated in all conditions (Table S2). The top eight most regulated were highly expressed (> 48 fold) and included two uncharacterized genes (*AALFPA_051205*, *AALFPA_066142*), one gene putatively involved in antiviral response (*AALFPA_05206*), one serine-rich protein kinase (*AALFPA_068585*), one deubiquitinase (*AALFPA_074126*), one cytochrome P450 (*AALFPA_062825*), one homologue of *hunchback* (*AALFPA_070528*) and one homolog of *canoe* (*AALFPA_068948*). Second, we looked at the effect of time post infection within mosquito species. We

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

<figure><img src="image_4.png" /><figcaption>c Between CHIKV
conditions</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>d Between *A. malayensis*
conditions</figcaption></figure>

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

**Figure 1.** The midgut transcriptomic responses of *A. albopictus* to CHIKV and *A. malayensis* to CHIKV and DENV at 1 and 4 dpi. (a) Quantification of viral reads. ***, *p* < 0.001 as determined by unpaired t-test. (b) Principal Component Analysis (PCA) using DESeq2 normalized gene expression values. (c, d) Venn diagrams of DEGs with conditions with CHIKV infection (c) and with conditions with *A. malayensis*.
found that 223 DEGs were conserved in *A. albopictus* between 1 and 4 dpi (Table S2), representing 12% and 66% of all DEGs for each dpi, respectively. We observed that 153 DEGs were common in *A. malayensis* between 1 and 4 dpi, representing 16% and 69% of within-dpi DEGs, respectively. Third, we reported the shared DEGs among mosquito species at each time point. At 1 dpi, 530 DEGs were common between *A. albopictus* and *A. malayensis*, representing 29% and 56% of all DEGs for each species, respectively. At 4 dpi, 130 DEGs were common, representing 7.3% and 13.8% of all DEGs for *A. albopictus* and *A. malayensis*, respectively. Overall, the transcriptomic impact of CHIKV infection is relatively conserved across the two mosquito species and a majority of DEGs at 4 dpi were already altered at 1 dpi.

We subsequently looked at the shared DEGs in *A. malayensis* infected with DENV or CHIKV at 1 and 4 dpi. Firstly, there was only one conserved DEG between all conditions (Fig. 1d) and it was the uncharacterized *AAFLPA_064235*, which was slightly upregulated in all conditions (Table S2). Secondly, we looked at the effect of time post infection. There were only two DEGs common between 1 and 4 dpi with DENV, which was less than between 1 and 4 dpi with CHIKV, as detailed above. Finally, we observed how the virus influenced the transcriptome within each time point. There were 28 DEGs common between DENV and CHIKV at 1 dpi, representing 38% and 3% of all DEGs for DENV and CHIKV, respectively. There were 19 shared DEGs between the two viruses at 4 dpi, representing 28% and 9% of DENV- and CHIKV-induced DEGs. Overall, in *A. malayensis*, DENV induced a different transcriptomic response at 1 and 4 dpi and the transcriptomic signature of DENV and CHIKV differed.

**Functional annotations of genes responsive to CHIKV and DENV infections in *A. albopictus* and *A. malayensis* midguts.** We annotated the DEGs into different functional categories based on characterized functions, homology, or protein domains. For both comparative purposes and because of the absence of an annotated *A. malayensis* genome, we relied on the orthologs in *A. albopictus* to describe virus-induced gene alterations in *A. malayensis*. Apoptotic-related genes were mostly regulated by CHIKV in *A. albopictus* at 1 dpi (Fig. 2; Table S2) with the up- and downregulation of apoptosis-inducers such as CASP. Interestingly, the expression of an inhibitor of apoptosis (IAP; AALFPA_076435) was increased in all CHIKV conditions. Only one apoptotic gene was induced by DENV at 1 dpi in *A. malayensis*.

Several genes related to blood-feeding were regulated by CHIKV infection at the different time points and in both mosquito species but only one blood-feeding gene was regulated by DENV at 1 dpi in *A. malayensis* (Fig. 2; Table S2). Particularly, one putative odorant binding (*AALFPA_066142*) was strongly (> 48 fold) induced in all CHIKV conditions and one gustatory receptor (*AALFPA_047124*) was similarly highly upregulated by CHIKV at both time points in *A. malayensis*. Genes related to cytoskeleton were mostly upregulated in all the CHIKV conditions. However, contrasting regulation of cytoskeleton genes were observed for DENV in *A. malayensis*, where they were downregulated 1 dpi, but upregulated at 4 dpi (Fig. 2; Table S2). Genes related to digestion were mostly regulated by CHIKV at 1 dpi in both mosquito species and in lower quantity by CHIKV at 4 dpi and DENV at 1 dpi (Fig. 2; Table S2). None of the digestion-related genes were regulated by DENV at 4 dpi. Of interest, two collagen genes and one chitinase gene were upregulated in all CHIKV conditions. A large proportion of genes with diverse functions were regulated in all conditions.

Among the genes related to immunity, in CHIKV-infected *A. albopictus*, 36 DEGs were upregulated and 21 downregulated at 1 dpi, and 12 were upregulated at 4 dpi (Fig. 2; Table S2). In CHIKV-infected *A. malayensis*, 30 DEGs were induced and five were reduced at 1 dpi, while expression of nine were increased at 4 dpi. In DENV-infected *A. malayensis*, there was one upregulated immune DEGs and five downregulated at 1 dpi, and

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Distribution of the DEGs according to their functional group. APO, apoptosis; BF, bite and feeding, CS, cytoskeleton; DIV, diverse; DIG, digestion; IMM, immunity; MET, metabolism; PROT, proteolysis; RSM, redox, stress and mitochondria; RTT, replication, transcription and translation. TRP, transport; UNK, unknown.</figcaption></figure>
four increased DEGs at 4 dpi. Among these DEGs, there were several PRRs. One GNBP and three PGRPs were downregulated by CHIKV at 1 dpi in *A. albopictus*. A variety of serine protease or serine protease inhibitor DEGs were regulated in either direction by CHIKV and by DENV in both mosquito species. For IMD pathway, one activator (*sickie-like*; AALFPA_052752) was inhibited and another one (*hyperplastic discs*; AALFPA_079063) was activated at 1 dpi in CHIKV-infected *A. albopictus*. Another activator (*sickie-like*; AALFPA_058168) was upregulated by CHIKV at 4 dpi in both mosquito species. For Toll pathway, two Toll-like genes (AALFPA_053631 and AALFPA_057415) were induced by CHIKV at 1 dpi in *A. albopictus*. For Jak/STAT pathway, three repressors (three Suppressor of cytokine signaling (SOCS)-like; AALFPA_051177; AALFPA_054706; AALFPA_056404) were activated at 1 dpi with CHIKV and at 4 dpi with DENV in *A. malayensis*, whereas one Bromodomain and WD repeat domain containing 3 gene (*BRWD3*, AALFPA_050507) was induced by CHIKV at 1 dpi in *A. albopictus*. For JNK pathway, two *Kayak-like* genes (AALFPA_079058 and AALFPA_079354) were upregulated by CHIKV at 1 dpi in both mosquito species and by DENV at 4 dpi in *A. malayensis*. At 1 dpi with CHIKV, two other JNK activators were induced in either mosquito species (*A. albopictus*: AALFPA_071724 and AALFPA_078630; *A. malayensis*: AALFPA_071724 and AALFPA_040889). In terms of immune effectors, two complement genes (AALFPA_045642 and C1q, AALFPA_055182) were increased by CHIKV at 1 dpi in *A. albopictus*, but no antimicrobial peptide (AMP) was directly regulated. Moreover, one helicase (AALFPA_080443), involved in RNAi, was induced at 1 dpi with CHIKV in *A. albopictus*. The same helicase plus another (AALFPA_066096) were upregulated at 1 dpi with CHIKV in *A. malayensis*. Overall, although we observed the regulations of components from several immune pathways, there was a consistent induction of JNK activators.

**Immune regulations in *A. albopictus* midgut infected by CHIKV.** Because the immune response chiefly determines vector competence³²,³³, we focused on immunity. We identified the immune genes transcribed in *A. albopictus* midguts in both infected and control conditions, and classified them based on the categories defined by Palatini et al.⁵⁴ (Table 1; Table S3). A majority of the immune genes identified in the *A. albopictus* genome were detected in our RNA-seq dataset (463 out of 664; not including JNK-related genes). Interestingly, most of the PRRs (GNBPs and PGRPs), all the components of the cytoplasmic signaling of the immune path-

<table><caption>**Table 1.** List of immune gene orthologs identified in *A. albopictus* and *A. malayensis* midgut. Gene identification was based on de novo assembled transcripts for *A. malayensis*. <sup>a</sup>Data extracted from⁵⁴.</caption><thead><tr><th rowspan="2">Immune category</th><th rowspan="2">Abbreviations</th><th rowspan="2">A. albopictus genome<sup>a</sup></th><th colspan="2">Number expressed in midguts</th></tr><tr><th>A. albopictus</th><th>A. malayensis</th></tr></thead><tbody><tr><td>Antimicrobial peptides</td><td>AMP</td><td>6</td><td>5</td><td>6</td></tr><tr><td>Autophagy pathway members</td><td>APHAG</td><td>30</td><td>29</td><td>14</td></tr><tr><td>Caspase activators</td><td>CASPA</td><td>3</td><td>3</td><td>3</td></tr><tr><td>Caspases</td><td>CASP</td><td>24</td><td>21</td><td>7</td></tr><tr><td>Catalases</td><td>CAT</td><td>2</td><td>1</td><td>3</td></tr><tr><td>CLIP-domain serine proteases</td><td>CLIP</td><td>118</td><td>63</td><td>9</td></tr><tr><td>C-type lectins</td><td>CTL</td><td>66</td><td>24</td><td>3</td></tr><tr><td>Fibrinogen-related proteins</td><td>FREP</td><td>52</td><td>35</td><td>6</td></tr><tr><td>Galectins</td><td>GALE</td><td>12</td><td>11</td><td>8</td></tr><tr><td>Gram-negative binding proteins</td><td>GNBP</td><td>13</td><td>12</td><td>2</td></tr><tr><td>Heme peroxidases</td><td>HPX</td><td>36</td><td>33</td><td>7</td></tr><tr><td>IMD pathway members</td><td>IMDPATH</td><td>16</td><td>16</td><td>9</td></tr><tr><td>Inhibitors of apoptosis</td><td>IAP</td><td>5</td><td>5</td><td>3</td></tr><tr><td>Jak/STAT pathway members</td><td>JAKSTAT</td><td>3</td><td>3</td><td>1</td></tr><tr><td>JNK pathway members</td><td>JNK</td><td>Not in Paper</td><td>13</td><td>11</td></tr><tr><td>Lysozymes</td><td>LYS</td><td>7</td><td>5</td><td>2</td></tr><tr><td>MD2-like proteins</td><td>ML</td><td>36</td><td>22</td><td>6</td></tr><tr><td>Peptidoglycan receptor-like proteins</td><td>PGRP</td><td>18</td><td>17</td><td>10</td></tr><tr><td>Prophenoloxidases</td><td>PPO</td><td>23</td><td>13</td><td>1</td></tr><tr><td>Rel-like NFkappa-B proteins</td><td>REL</td><td>4</td><td>3</td><td>2</td></tr><tr><td>Scavenger receptors</td><td>SCR</td><td>28</td><td>22</td><td>7</td></tr><tr><td>Serine protease inhibitors</td><td>SRPN</td><td>46</td><td>30</td><td>8</td></tr><tr><td>SiRNA pathway members</td><td>SRRP</td><td>51</td><td>45</td><td>20</td></tr><tr><td>Superoxide dismutases</td><td>SOD</td><td>9</td><td>5</td><td>2</td></tr><tr><td>Spaetzle-like proteins</td><td>SPZ</td><td>13</td><td>9</td><td>5</td></tr><tr><td>Thioester-containing proteins</td><td>TEP</td><td>7</td><td>5</td><td>1</td></tr><tr><td>Toll pathway members</td><td>TOLLPATH</td><td>10</td><td>10</td><td>4</td></tr><tr><td>Toll receptors</td><td>TOLL</td><td>26</td><td>16</td><td>2</td></tr><tr><td>Total</td><td></td><td>664</td><td>476</td><td>162</td></tr></tbody></table>
ways (IMDPATH, TOLLPATH and JAKSTAT), 13 components of the JNK pathway, and three out of four REL downstream transcription factors were expressed in the midguts. Inversely, we detected a lower proportion of transcripts among enzymes involved in the regulation of the immune signaling, i.e., 63 of 118 CLIPs and 30 of 46 SRPNs. SPZ (nine of 13) and TOLL receptors (16 of 26) were also expressed in lower proportions. Effectors such as AMP, TEP and GALE had most of their orthologs expressed in midguts, whereas CTL, FREP, ML and PPO had lower proportions of orthologs detected. Members of ROS system were expressed in midguts with 1 out of 2 CATs, 33 out of 36 HPXs, 22 out of 28 SCRs and 5 out of 9 SODs. Moreover, a large proportion of SRRPs and most of APHAGs and apoptosis factors (CASPAs, CASPs and IAPs) were transcribed in midguts, completing the picture of the immune machinery activated in midguts.

To identify patterns in the immune response to CHIKV infection at 1 and 4 dpi in *A. albopictus*, we clustered and plotted the expressions of the 50 most variable immune genes in a heatmap (Fig. 2; Table S3). We detected 5 clusters. Cluster 1 contained only one C-type lysozyme that also belonged to the significant DEGs identified above and was highly upregulated (> 77 fold) at both time points. Cluster 2 was made of five genes only down-regulated at 4 dpi and included two PRRs (one PGRP and one GNBP), two proteases involved in the regulation of the immune signaling (one CLIP and one SRPN), and one LYS. Cluster 3 consisted of 12 genes induced at 1 dpi and then reduced at 4 dpi. These included several genes implicated in regulating the immune activation, i.e., three SRPNs (two of which were significant DEGs), one SPZ and one TOLL. Their clustering with effectors such as one AMP (Attacin-B), one CTL and two FREPs, suggested their regulations through the above enzymes. Cluster 3 also included two SCRs involved in ROS regulation and one PPO that induces melanization. Inversely to cluster 3, cluster 4 was made of 14 genes downregulated at 1 dpi and upregulated at 4 dpi. Among the proteins that regulate immune signaling, there were one PGRP, one SPZ, one SRPN and one CLIP. Six immune effectors were also grouped in the same cluster and included three FREPs and three MLs. Three CASPs that trigger apoptosis were present in cluster 4, together with one SRRP. Finally, cluster 5 contained 18 genes only downregulated at 1 dpi. Genes involved in immune signaling regulation were two PGRPs, one TOLL, three SRPNs and two CLIPs. One ML effector was also present. Regulations of ROS was apparent due to the clustering of five HPXs, regulation of melanization due to two PPOs, and regulation of RNAi because of two SRRP aubergine homologs.

**Immune regulations in *A. malayensis* midgut infected by CHIKV and DENV.** In *A. malayensis*, we undertook a similar immune-centered approach. However, to provide the most accurate description for this mosquito species that does not have an assembled genome, we de novo assembled transcripts before searching for protein homologs in *A. albopictus* transcriptome. We identified 162 immune proteins (Data S1 details the protein sequences) expressed in the midgut of *A. malayensis* and functionally categorized them as above⁵⁴ (Table 1; Table S4). There were a lower number of expressed homologs for all immune categories than in *A. albopictus*, except for AMP. These lower numbers of transcripts may be due to the conservative parameters (see "Methods") we used for homolog identification. Nonetheless, for immune signaling, we found 12 PRRs (2 GNBPs and 10 PGRPs), 5 SPZs, 2 TOLLS, 9 components of the IMDPATH, 1 of the JAKSTAT, 11 of the JNK, 4 of the TOLLPATH and 2 REL homologs. Among enzymes that regulate the signaling, we detected 9 CLIPs and 8 SRPNs. For effectors, we noted 6 AMPs, 3 CTLs, 6 FREPs, 8 GALEs, 2 LYSs, 6 MLs and 1 TEP. Moreover, our results indicated a functional RNAi with 20 SRRPs. Autophagy, apoptosis, melanization and ROS regulation were also activated with 14 APHAGs, 13 apoptosis-related (3 CASPAs, 7 CASPs, 3 IAPs), 1 PPO, and 19 oxidative stress-related transcripts (3 CATs, 7 HPXs, 7 SCRs and 2 SODs), respectively.

We next clustered the expression of the 50 most variable immune genes upon CHIKV (Fig. 4a; Table S5) and DENV (Fig. 4b; Table S5) infection, separately. With CHIKV-infected samples, we identified four clusters. Cluster 1 contained ten genes only upregulated at 1 dpi. Among the components of the immune pathway signaling, there were two SRPNs, one JNK component, and one REL transcription factor. For effectors, one LYS, one FREP and one GALE were clustered. The detection of two CASPs and Ago2, a major SRRP, suggested apoptosis and RNAi activation, respectively. Cluster 2 was made of 16 genes induced at 1 dpi and then reduced at 4 dpi. For signaling regulation, there were one GNBP, three SPZs, one CLIP and one homolog of Kayak, the JNK transcription factor. Multiple effectors were in cluster 2 and included two AMPs, one TEP, one FREP, one LYS, one GALE and one ML. There were also one APHAG and two SCRs. Cluster 3 was composed of eight genes upregulated only at 4 dpi. For immune signaling, there were one GNBP, one TOLL and one CLIP. For effectors, there was only one FREP. Additionally, there were two ROS-related genes (one SOD and one SCR), one IAP and one APHAG. Cluster 4 consisted of 16 genes only downregulated at 1 dpi. Those related to immune signaling included one CTL, two SRPNs, one IMDPATH and Cactus, the negative regulator of TOLLPATH transcription. Among effectors, there were two MLs. One APHAG and three HPXs were also clustered. Interestingly, five SRRPs were grouped in cluster 4 and included Ago3, loquacious and two aubergine homologs.

With DENV-infected samples, we identified five clusters (Fig. 4b; Table S5). Cluster 1 consisted of six genes moderately upregulated at 1 dpi and strongly upregulated at 4 dpi. Among signaling genes, there was only one JNK component. Effectors included one ML and one GALE. There were also one CASPA, one APHAG and one SCR. Cluster 2 included 11 genes only upregulated at 4 dpi. No immune signaling-related genes were in cluster 2 and there were one ML and one GALE effectors. There were three apoptosis-related genes (two IAPs and one CASP), two APHAGs, one SOD and three RNAi-related genes (SRRPs). Cluster 3 was made of three genes induced at 1 dpi and then reduced at 4 dpi. These included one AMP, one SPZ and one LYS. Cluster 4 consisted of seven genes only upregulated at 1 dpi. Immune signaling genes included one GNBP and one SPZ, while immune effectors were one TEP, three FREPs and one AMP. Eventually, cluster 5 was composed of 23 genes moderately upregulated at 1 dpi and moderately downregulated at 4 dpi. Immune signaling genes included two PGRPs, three SPZs, three CLIPs, two SRPNs and two IMDPATHs such as an IKK-β homolog. Effectors consisted
of two AMPs, three MLs and one CTL. There were also one APHAG, two CASPs and two ROS-related genes (one CAT and one SCR).

Discussion

Public health threats from dengue and chikungunya are expanding geographically and in intensity³,⁵⁶. While *A. aegypti* remains the main vector, other *Aedes* species are competent enough to trigger moderate-scale epidemics and sustain DENV and CHIKV transmissions⁸,⁵⁷. Furthermore, the ongoing deployment of new interventions aimed at *A. aegypti* like the *Wolbachia*-mediated transmission reduction⁵⁸,⁵⁹ will most probably promote the roles of these secondary vectors. Here, we set to prepare for the next steps of vector control by comprehensively describing the transcriptomic responses to CHIKV and DENV in two capable vectors, *A. albopictus* and *A. malayensis*²⁷. Because infection onset in the mosquito midgut dictates vector competence, we analyzed midguts at early times post infection, i.e., 1 and 4 dpi. Our results profile DEGs in two mosquito species infected by two viruses at two time points. Next, we focused on immune-related genes and utilized a de novo transcriptome assembly in *A. malayensis* to identify immune genes for the first time in this species. Overall, our study reveals patterns of transcriptomic responses in two *Aedes* vectors and details the immune responses, shedding new light on shared and divergent mechanisms of arboviral transmission.

The overall transcriptomic response was shaped by viral infections as well as time post infection and varied between the two mosquito species. Although we orally infected mosquitoes with a similar inoculum for both viruses, the transcriptomic response was greater upon CHIKV than DENV infection. Such difference may relate to the faster CHIKV replication kinetic⁵², which reached a plateau as early as 1 dpi in both mosquito species. However, in spite of the similar viral genomic RNA (gRNA) levels between the two time points, the number of CHIKV-induced DEGs did not correlate with gRNA copies across the time points and sharply decreased from 1 to 4 dpi. This suggests a return to homeostasis after quick establishment of viral replication factories. Interestingly, CHIKV-induced DEGs at 4 dpi were mostly already regulated at 1 dpi in both mosquito species. These conserved regulated genes may represent the cellular machinery required for CHIKV multiplication. In DENV-infected *A. malayensis* midguts, the transcriptomic response did not correlate with gRNA quantities, which grew significantly between 1 and 4 dpi whereas DEG numbers were similar between the time points. In contrary to CHIKV infection, there was a very limited number of shared DEGs between 1 and 4 dpi in DENV-infected *A. malayensis*. There was also little overlap between DENV- and CHIKV-induced DEGs in this mosquito species. These transcriptomic differences document divergences in how DENV and CHIKV establish infection in mosquito midguts.

Strikingly, we identified a transcriptomic signature of CHIKV infection. There were 81 DEGs regulated in the same direction upon CHIKV infection in the two mosquito species (comparison was made using homologs from both species) and at both time points (Fig. 1; Table S2). While these common DEGs were related to multiple functional groups, we discuss the top eight most regulated ones as they were highly induced in all conditions. AALFPA_052026 and its *A. malayensis* homolog have a *D. melanogaster* homolog (CG8492) that reduces viral infection⁶⁰. When we detailed the immune response, this gene did not cluster with other immune genes in *A. albopictus* (Fig. 3). However, in *A. malayensis*, it responded similarly to a Relish homolog (Fig. 4a), suggesting a potential regulation through one of the immune pathways that induces Relish. A deubiquitinase was part of the CHIKV transcriptomic signature and may have influenced the immune response and cellular homeostasis, given the function of ubiquitination in both processes⁶¹,⁶². CHIKV infection consistently induced a serine-rich protein kinase (SRPK), which may influence gene expression through chromatin remodeling⁶³. A cytochrome P450 monooxygenase was upregulated and could have metabolized xenobiotic substances⁶⁴,⁶⁵. Interestingly, two highly expressed genes of the CHIKV signature are related to cellular development, i.e., a homolog of hunch-back, a transcription factor that regulates cellular development⁶⁶,⁶⁷, and a homolog of canoe, which maintains adherens junction between cells. Both genes may be required for the cell multiplication that occurs in infected mosquito midguts to compensate for damaged cells⁶⁸. Together with two other uncharacterized homologs, these top induced DEGs shared across all our CHIKV conditions deserve further functional studies.

Given the importance of midgut immunity in determining vector competence³⁰,³²,³³, we detailed the immune response in both mosquito species. In *A. albopictus*, most of the components of the IMD, TOLL, JAK/STAT and JNK pathways were expressed, indicating that the corresponding pathways are functional in midguts. In contrary, a lower proportion of enzymes involved in signaling regulation were detected and this suggests a midgut-specific regulation of immunity as different CLIP or SRPN interact with different immune components⁶⁹,⁷⁰. Apoptosis, autophagy, RNAi, ROS and melanization were also functional in midguts. In *A. malayensis*, smaller proportions of the immune pathway components were detected in midguts. While this may be related to the technical limitations associated with de novo assembly such as high error rates and short assembled contigs⁷¹, each canonical pathway was still represented by several key members. For TOLL and IMD pathways, the two downstream transcription factors, homologs of *D. melanogaster* Dorsal and Relish, respectively, were expressed in *A. malayensis* midguts (Table S4). The negative regulators for Dorsal and Relish, namely Caspar and Cactus respectively, were also found. Similarly, the downstream transcription activators for JNK and JAK/STAT, i.e., Kayak and STAT respectively, were transcribed, confirming the conservation of these pathways across insects⁷² and indicating that the canonical immune pathways function in *A. malayensis* midguts.

Activation of the TOLL, IMD, JAK/STAT and JNK pathways occurs through transcriptional regulation of their components in response to infection³⁴,⁴³,⁷³,⁷⁴. Interestingly, the JNK pathway appeared to be induced in both mosquito species and upon both viral infections. Two *Kayak*-like homologs were significantly upregulated (i.e., DEG) at 1 dpi with CHIKV in both mosquito species and at 4 dpi with DENV in *A. malayensis*. In *A. aegypti* salivary glands, Kayak was previously shown to be induced by both CHIKV and DENV and to be co-regulated with Basket, one of the upstream components of the JNK pathway³⁴. Supporting the activation of JNK, we also
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Immune gene regulation upon CHIKV infection in *A. albopictus* midgut at 1 and 4 dpi. Heatmap of the top 50 most regulated DEGs. Underlined genes were significantly regulated. Hierarchical clustering of gene expression is shown on the left.</figcaption></figure>

observed an induction of the *Basket* homolog (among the 50 most regulated immune genes) in *A. malayensis* at 1 dpi with CHIKV and 4 dpi with DENV. JNK activation clears viral infection by activating apoptosis and the complement system³⁴. Together with the JNK components, we observed co-regulations of multiple CASPs and
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Immune gene regulation upon CHIKV and DENV infection in *A. malayensis* midgut at 1 and 4 dpi. Heatmap of the top 50 most regulated DEGs. None of these gene were significantly regulated. Hierarchical clustering of gene expression is shown on the left.</figcaption></figure>

TEPs. Although our transcriptomic data indicate an activation of the JNK pathway, other immune pathways appeared induced as well. For instance, the induction of a *Relish* homolog together with the reduction of its inhibitor *Cactus* among the most regulated immune genes in *A. malayensis* at 1 dpi with CHIKV suggest the activation of the Toll pathway. The induction of *Ikk-β* in DENV-infected *A. malayensis* at 1 dpi show activation of the IMD pathway<sup>75</sup>. However, while components of several immune pathways were regulated, the activation of the JNK pathway was consistent in both *A. albopictus* and *A. malayensis* upon infection by CHIKV and DENV.

In conclusion, we conducted a high-throughput transcriptomic analysis in two *Aedes* vectors with DENV and CHIKV at two early time points in midguts. Our dataset provides a multidimensional picture of the transcriptomic and more specifically immune regulations in these conditions. In *A. albopictus*, we identified different gene regulation clusters, which may represent different immune pathway responses to CHIKV infection in midguts. In *A. malayensis*, we provided the first description of the immune response to arboviral infections. This knowledge is required to broaden our understanding of arboviral transmission by *Aedes* vectors.

## Material and methods

**Mosquito rearing.** *Aedes albopictus* and *A. malayensis* mosquito colonies were established in 2010 and 2014, respectively, from eggs collected using oviposition traps in the parks of Singapore<sup>27</sup>. Since then, these colonies were reared in the insectary where eggs were hatched in MilliQ water, larvae fed on a mixture of fish food (TetraMin fish flakes), yeast and liver powder (MP Biomedicals) and adults maintained on 10% sucrose and fed pig's blood twice weekly. Mosquito colonies were maintained at 28 °C and 50% relative humidity with a 12 h:12 h light: dark cycle.
**Viruses.** Dengue virus serotype 2 PVP110 was isolated from a patient enrolled in the Early DENGue infection and outcome study (EDEN) conducted in Singapore in 2008⁷⁶. Chikungunya virus SGP011 was isolated from a patient at the National University Hospital in Singapore in 2008⁵⁰. DENV isolates were propagated in C6/36 (CRL-1660) and CHIKV in Vero (CCL-81) cell lines. Virus stocks were titered with BHK-21 cell plaque assay as previously described⁷⁷, aliquoted and stored at -80 °C. Virus stocks were thawed only once for infection.

**Mosquito infection.** Three-to-five day-old female mosquitoes were starved for 24 h before they were fed on an infectious blood meal containing 40% volume of washed erythrocytes from specific pathogen free (SPF) pig's blood (PWG Genetics), 5% 10 mM ATP (Thermo Scientific), 5% human serum (Sigma) and 50% virus solution in RPMI media (Gibco), using Hemotek membrane feeder system (Discovery Workshops) covered with porcine intestine membrane (sausage casing). The virus titers in blood meals were 2 × 10⁷ pfu/ml for DENV and 1 × 10⁷ pfu/ml for CHIKV. Blood titers were validated by plaque assay using BHK-21 cells. Control mosquitoes were fed with the same blood meal composition except for virus solution, which was replaced by RPMI. Following oral feeding, the fully engorged females were selected and kept in a cage with ad libitum access to a 10% sucrose solution in an incubation chamber with conditions similar to insect rearing.

**Detection of Wolbachia.** Groups of 10 mosquitoes were homogenized using a bead mill homogenizer (Biospec Products) and total DNA was extracted with QIAamp DNA extraction kit (Qiagen). *Wolbachia* presence was detected by PCR amplification with GoTaq Master Mix (Promega) and the following primers: Forward (CAT ACC TAT TCG AAG GGA TAG; Rev-AGC TTC GAG TGA AAC CAA TTC⁷⁸). Two replicates were conducted for *A. albopictus* and *A. malayensis* colonies.

**Midgut collection, library preparation and RNA-sequencing.** Midguts from control, DENV- and CHIKV-orally infected *A. albopictus* and *A. malayensis* were dissected at one- and four-days post infection (dpi). The blood was removed from the midguts. Three repeats per condition with 20 midguts per replicate were homogenized using a bead mil homogenizer (Biospec Products, USA). Total RNA was extracted using E.Z.N.A Total RNA kit I (OMEGA Bio-Tek). RNA-sequencing libraries were prepared using a True-Seq Stranded Total RNA with Ribo-Zero Gold kit (Illumina), according to manufacturer's instructions. Following quantification by RT-qPCR with a KAPA Library Quantification Kit (KAPA Biosystems), libraries were pooled in equimolar concentrations for cluster generation on cBOT system (Illumina) and sequenced (150 bp pair-end) on a HiSeq 3000 instrument (Illumina) at the Duke-NUS Genome Biology Facility.

**Identification of differentially expressed genes.** Sequenced read quality was assessed using the Java program FastQC (Babraham Institute Bioinformatics, UK). Contaminating adaptors and low-quality reads (Phred + 33 score < 20) were removed with Trimmomatic v0.39⁷⁹ using a sliding window of 4 bps. To confirm virus infection and mosquito species, reads were aligned with Bowtie2 v2.3.4.3⁸⁰ to the genomes of the viruses used for infection and *A. albopictus* and *A. malayensis* mitochondrial cytochrome c oxidase 1 (CO1) genes, respectively. Reads from both mosquito species were aligned on *A. albopictus* genome (AaloF1.2) and the *A. aegypti* genome (AaegL5.0) with STAR v2.5.4a in two-pass mode⁸¹ and HTseq v0.6.0⁸² used to produce count files, that were then used as input into DESeq²⁸³ to identify DEGs with at least a 1.4-fold (log2 = 0.48) change between control and infected conditions and an adjusted False Discovery Rate (FDR) of 0.05. DEGs were annotated by importing gene names, descriptions and computed GO functions and processes from VectorBase⁸⁴, and functional groups were manually assigned. Predicted protein sequences were also searched against the NCBI Diptera protein database (Blastp, e-value 1.0E-5 threshold) and the FlyBase⁸⁵ *D. melanogaster* peptide database for identification of orthologues.

**Validation of RNA-seq by real-time quantitative PCR.** *Aedes albopictus* and *A. malayensis* mosquitoes from different batches than the ones used for RNA-seq were orally infected with DENV and CHIKV. At four dpi, midguts from eight mosquitoes were dissected and pooled in triplicates. Total RNA was extracted from the midgut samples using a E.Z.N.A. Total RNA kit I (OMEGA), DNase treated using a Turbo DNA-free kit (Thermo Fisher Scientific), and reverse transcribed using iScript cDNA synthesis kit (Biorad). Expression for 11 genes was quantified with qPCR using the SensiFast Sybr no-rox kit (Bioline) and primers designed based on *A. albopictus* genome, as detailed in Table S6. *Actin* expression was used for normalization. Reactions were performed with the following cycle conditions: an initial 95 °C for 10 min, followed by 40 cycles of 95 °C for 5 s, 60 °C for 20 s and ending with a melting curve analysis. The delta delta method was used to calculate relative fold changes. Pearson correlations were calculated between qPCR-based and DEG log2 fold changes with excel.

**De novo identification of immune genes in *A. malayensis* and quantification of their expression.** For the de novo assembled *A. malayensis* transcriptome, reads were assembled with Trinity v2.4.0⁸⁶ and ABySS v2.1.5⁸⁷. Trinity parameters specified a contig size of at least 150 bps with a minimum k-mer coverage of 2. ABySS assemblies were completed for k-mer sizes 21, 31, 41, 51, 61, 71, 81 and 91 bps. Individual assemblies were then combined and redundancy removed with CD-HIT v4.6.8⁸⁸,⁸⁹. Blastn v2.12⁹⁰ was used to identify contigs belonging to the DENV2 and CHIKV genomes used for infection, and these were subsequently removed. Transdecoder from the Trinity pipeline was then used to identify all contigs with open reading frames (ORFs), and only those with complete and nonredundant ORFs were retained. Blastp was then used to identify translated contigs matching (e-value 1.0E-5) to NCBI Diptera proteins (taxid:7147, accessed June 2021), and all others were filtered out. Using RSEM v1.3.0⁹¹, reads were mapped to this final set and used as input into DESeq²⁸³ to deter-
mine log2 Fold-changes between infected and control conditions and DEGs with the same criteria as detailed above.

To confidently annotate all *A. malayensis* immune transcripts, the protein products of *A. albopictus* immune genes identified by Palatini et al.⁵⁴ and an additional 13 JNK pathway genes were used as a database to identify orthologues from the translated *A. malayensis* de novo assembled transcriptome. This was done with Diamond Blastp⁹² and a strict search criteria of an e-value 1.0E-¹⁰ or less, followed by a second search against all annotated *A. albopictus* proteins to confirm identifications. From this identified set, only transcripts that could be fully translated with no ambiguities, covered at least 90% of database subject sequences, and shared less than 98% identity with each other were considered confident *A. malayensis* immune genes. These were also searched against the FlyBase⁸⁵ *D. melanogaster* peptide database for identification of orthologues.

## Data availability

Sequenced reads generated during this study are available under NCBI accessions: SRR14621613-SRR14621644, BioProject PRJNA731987.

Received: 15 December 2022; Accepted: 6 July 2023

Published online: 12 July 2023

## References

1. Mayer, S. V., Tesh, R. B. & Vasilakis, N. The emergence of arthropod-borne viral diseases: A global prospective on dengue, chikungunya and Zika fevers. *Acta Trop.* **166**, 155-163 (2017).

2. Paixão, E. S., Teixeira, M. G. & Rodrigues, L. C. Zika, chikungunya and dengue: The causes and threats of new and re-emerging arboviral diseases. *BMJ Glob. Health* **3**, e000530 (2018).

3. Bettis, A. A. et al. The global epidemiology of chikungunya from 1999 to 2020: A systematic literature review to inform the development and introduction of vaccines. *PLoS Negl. Trop. Dis.* **16**, e0010069 (2022).

4. Messina, J. P. et al. The current and future global distribution and population at risk of dengue. *Nat. Microbiol.* **4**, 1508-1515 (2019).

5. Semenza, J. C. & Suk, J. E. Vector-borne diseases and climate change: a European perspective. *FEMS Microbiology Letters* **365**, fnx244 (2018).

6. Tabachnick, W. J. Climate change and the arboviruses: Lessons from the evolution of the dengue and yellow fever viruses. *Annu. Rev. Virol.* **3**, 125-145 (2016).

7. Bhatt, S. et al. The global distribution and burden of dengue. *Nature* **496**, 504-507 (2013).

8. Silva, L. A. & Dermody, T. S. Chikungunya virus: Epidemiology, replication, disease mechanisms, and prospective intervention strategies. *J. Clin. Investig.* **127**, 737-749 (2017).

9. Halstead, S. B. Dengvaxia sensitizes seronegatives to vaccine enhanced disease regardless of age. *Vaccine* **35**, 6355-6358 (2017).

10. Thomas, S. J. & Yoon, I.-K. A review of Dengvaxia*: Development to deployment. *Hum. Vaccin. Immunother.* **15**, 2295-2314 (2019).

11. Bowman, L. R., Donegan, S. & McCall, P. J. Is dengue vector control deficient in effectiveness or evidence? Systematic review and meta-analysis. *PLoS Negl. Trop. Dis.* **10**, e0004551 (2016).

12. Ooi, E.-E., Goh, K.-T. & Gubler, D. J. Dengue prevention and 35 years of vector control in Singapore. *Emerg. Infect. Dis.* **12**, 887-893 (2006).

13. Reiner, R. C. Jr. et al. Quantifying the epidemiological impact of vector control on dengue. *PLoS Negl. Trop. Dis.* **10**, e0004588 (2016).

14. Rosen, L., Roseboom, L. E., Gubler, D. J., Lien, J. C. & Chaniotis, B. N. Comparative susceptibility of mosquito species and strains to oral and parenteral infection with dengue and Japanese encephalitis viruses. *Am. J. Trop. Med. Hyg.* **34**, 603-615 (1985).

15. Jourdain, F. et al. From importation to autochthonous transmission: Drivers of chikungunya and dengue emergence in a temperate area. *PLoS Negl. Trop. Dis.* **14**, e0008320 (2020).

16. Metselaar, D. et al. An outbreak of type 2 dengue fever in the Seychelles, probably transmitted by *Aedes albopictus* (Skuse). *Bull. World Health Organ.* **58**, 937-943 (1980).

17. Peng, H.-J. et al. A local outbreak of dengue caused by an imported case in Dongguan China. *BMC Public Health* **12**, 83 (2012).

18. Tsuda, Y. et al. Biting density and distribution of *Aedes albopictus* during the September 2014 outbreak of dengue fever in Yoyogi park and the vicinity of Tokyo metropolis, Japan. *Jap. J. Infect. Dis.* **69**, 1-5 (2016).

19. Yue, Y., Liu, Q., Liu, X., Zhao, N. & Yin, W. Dengue fever in Mainland China, 2005-2020: A descriptive analysis of dengue cases and *Aedes* data. *Int. J. Environ. Res. Public Health* **19**, 3910 (2022).

20. Hafsia, S. et al. Overview of dengue outbreaks in the southwestern Indian Ocean and analysis of factors involved in the shift toward endemicity in Reunion Island: A systematic review. *PLoS Negl. Trop. Dis.* **16**, e0010547 (2022).

21. Paupy, C., Girod, R., Salvan, M., Rodhain, F. & Failloux, A.-B. Population structure of *Aedes albopictus* from La Réunion Island (Indian Ocean) with respect to susceptibility to a dengue virus. *Heredity* **87**, 273-283 (2001).

22. Lazzarini, L. et al. First autochthonous dengue outbreak in Italy, August 2020. *Eurosurveillance* **25**, 2001606 (2020).

23. Cochet, A. et al. Autochthonous dengue in mainland France, 2022: Geographical extension and incidence increase. *Eurosurveillance* **27**, 2200818 (2022).

24. Riccò, M., Peruzzi, S., Balzarini, F., Zaniboni, A. & Ranzieri, S. Dengue fever in Italy: The 'Eternal Return' of an emerging arboviral disease. *Trop. Med. Infect. Dis.* **7**, 10 (2022).

25. Tsetsarkin, K. A., Chen, R. & Weaver, S. C. Interspecies transmission and chikungunya virus emergence. *Curr. Opin. Virol.* **16**, 143-150 (2016).

26. Huang, Y.-M. Contributions to the Mosquito Fauna of Southeast Asia. XIV. The Subgenus Stegomyia of Aedes in Southeast Asia I - The Scutellaris Group of Species. https://apps.dtic.mil/sti/citations/ADA510169 (1972).

27. Mendenhall, I. H. et al. Peridomestic *Aedes malayensis* and *Aedes albopictus* are capable vectors of arboviruses in cities. *PLoS Negl. Trop. Dis.* **11**, e0005667 (2017).

28. Kek, R. et al. Feeding host range of *Aedes albopictus* (Diptera: Culicidae) demonstrates its opportunistic host-seeking behavior in rural Singapore. *J. Med. Entomol.* **51**, 880-884 (2014).

29. Miot, E. F. et al. A peridomestic *Aedes malayensis* population in Singapore can transmit yellow fever virus. *PLoS Negl. Trop. Dis.* **13**, e0007783 (2019).

30. Franz, A. W., Kantor, A. M., Passarelli, A. L. & Clem, R. J. Tissue barriers to arbovirus infection in mosquitoes. *Viruses* **7**, 3741-3767 (2015).

31. Ramesh, K. et al. Increased mosquito midgut infection by dengue virus recruitment of plasmin is blocked by an endogenous Kazal-type inhibitor. *iScience* **21**, 564-576 (2019).

32. Cheng, G., Liu, Y., Wang, P. & Xiao, X. Mosquito defense strategies against viral infection. *Trends Parasitol.* **32**, 177-186 (2016).

33. Sim, S., Jupatanakul, N. & Dimopoulos, G. Mosquito immunity against arboviruses. *Viruses* **6**, 4479-4504 (2014).
34. Chowdhury, A. et al. JNK pathway restricts DENV2, ZIKV and CHIKV infection by activating complement and apoptosis in mosquito salivary glands. *PLoS Pathog.* **16**, e1008754 (2020).

35. Lemaitre, B. & Hoffmann, J. The host defense of *Drosophila melanogaster*. *Ann. Rev. Immunol.* **25**, 697–743 (2007).

36. Dong, Y. et al. The *Aedes aegypti* siRNA pathway mediates broad-spectrum defense against human pathogenic viruses and modulates antibacterial and antifungal defenses. *PLoS Biol.* **20**, e3001668 (2022).

37. McFarlane, M. et al. Characterization of *Aedes aegypti* innate-immune pathways that limit chikungunya virus replication. *PLoS Negl. Trop. Dis.* **8**, e2994 (2014).

38. Luplertlop, N. 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.* **7**, e1001252 (2011).

39. Weng, S.-C. et al. A thioester-containing protein controls dengue virus infection in *Aedes aegypti* through modulating immune response. *Front. Immunol.* **12**, 670122 (2021).

40. Bonizzoni, M. et al. Complex modulation of the *Aedes aegypti* transcriptome in response to dengue virus infection. *PLoS ONE* **7**, 1–14 (2012).

41. 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 Genomics* **18**, 382 (2017).

42. 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. *Inter. J. Mol. Sci.* **20**, 3133 (2019).

43. Ramirez, J. L. & Dimopoulos, G. The toll immune signaling pathway control conserved anti-dengue defenses across diverse *Ae. aegypti* strains and against multiple dengue virus serotypes. *Dev. Comp. Immunol.* **34**, 625–629 (2010).

44. Sim, S. et al. Transcriptomic profiling of diverse *Aedes aegypti* strains reveals increased basal-level immune activation in dengue virus-refractory populations and identifies novel virus-vector molecular interactions. *PLoS Neg. Trop. Dis.* **7**, e2295 (2013).

45. Tsujimoto, H. et al. Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito. *Aedes albopictus*. *PLoS ONE* **12**, e0171345 (2017).

46. Liu, Y.-X. et al. Integrated analysis of miRNAs and transcriptomes in *Aedes albopictus* midgut reveals the differential expression profiles of immune-related genes during dengue virus serotype-2 infection. *Insect Sci.* **23**, 377–385 (2016).

47. Vedururu, R. K. et al. RNASeq analysis of *Aedes albopictus* mosquito midguts after chikungunya virus infection. *Viruses* **11**, 513 (2019).

48. Tun, Y. M. et al. Virological, serological and clinical analysis of chikungunya virus infection in Thai patients. *Viruses* **14**, 1805 (2022).

49. Vaughn, D. W. et al. Dengue viremia titer, antibody response pattern, and virus serotype correlate with disease severity. *J. Infect. Dis.* **181**, 2–9 (2000).

50. Her, Z. et al. Active infection of human blood monocytes by chikungunya virus triggers an innate immune response. *J. Immunol.* **184**, 5903–5913 (2010).

51. Aliota, M. T. et al. The wMel strain of Wolbachia reduces transmission of chikungunya virus in *Aedes aegypti*. *PLoS Negl. Trop. Dis.* **10**, e0004677 (2016).

52. Merwaiss, F. et al. Chikungunya virus replication rate determines the capacity of crossing tissue barriers in mosquitoes. *J. Virol.* **95**, e01956-e2020 (2021).

53. Salazar, M. I., Richardson, J., Sanchez-Vargas, I., Olson, K. & Beaty, B. Dengue virus type 2: Replication and tropisms in orally infected *Aedes aegypti* mosquitoes. *BMC Microbiol.* **7**, 9 (2007).

54. Palatini, U. et al. Improved reference genome of the arbovirus vector *Aedes albopictus*. *Genome Biol.* **21**, 215 (2020).

55. Matthews, B. J. et al. Improved reference genome of *Aedes aegypti* informs arbovirus vector control. *Nature* **563**, 501–507 (2018).

56. Pierson, T. C. & Diamond, M. S. The continued threat of emerging flaviviruses. *Nat. Microbiol.* **5**, 796–812 (2020).

57. Fontenille, D. & Powell, J. R. From anonymous to public enemy: How does a mosquito become a feared arbovirus vector?. *Pathogens* **9**, 265 (2020).

58. Lambrechts, L. et al. Assessing the epidemiological effect of wolbachia for dengue control. *Lancet Infect. Dis.* **15**, 862–866 (2015).

59. Indriani, C. et al. Reduced dengue incidence following deployments of Wolbachia-infected *Aedes aegypti* in Yogyakarta, Indonesia: A quasi-experimental trial using controlled interrupted time series analysis. *Gates Open Res.* **4**, 50 (2020).

60. Martins, N. E. et al. Host adaptation to viruses relies on few genes with different cross-resistance properties. *Proceed. Natl. Acad. Sci.* **111**, 5938–5943 (2014).

61. Hu, H. & Sun, S.-C. Ubiquitin signaling in immune responses. *Cell Res.* **26**, 457–483 (2016).

62. Cruz Walma, D. A., Chen, Z., Bullock, A. N. & Yamada, K. M. Ubiquitin ligases: Guardians of mammalian development. *Nat. Rev. Mol. Cell Biol.* **23**, 350–367 (2022).

63. Loh, B. J., Cullen, C. F., Vogt, N. & Ohkura, H. The conserved kinase SRPK regulates karyosome formation and spindle microtubule assembly in Drosophila oocytes. *J. Cell Sci.* **125**, 4457–4462 (2012).

64. Stavropoulou, E., Pircalabioru, G. G. & Bezirtzoglou, E. The role of cytochromes P450 in infection. *Front. Immunol.* **9**, 89 (2018).

65. Scott, J. G. & Wen, Z. Cytochromes P450 of insects: The tip of the iceberg. *Pest Manag. Sci.* **57**, 958–967 (2001).

66. Lehmann, R. & Nüsslein-Volhard, C. hunchback, a gene required for segmentation of an anterior and posterior region of the Drosophila embryo. *Dev. Biol.* **119**, 402–417 (1987).

67. Torres-Oliva, M., Schneider, J., Wiegleb, G., Kaufholz, F. & Posnien, N. Dynamic genome wide expression profiling of Drosophila head development reveals a novel role of Hunchback in retinal glia cell development and blood-brain barrier integrity. *PLoS Genet.* **14**, e1007180 (2018).

68. Janeh, M., Osman, D. & Kambris, Z. Comparative analysis of midgut regeneration capacity and resistance to oral infection in three disease-vector mosquitoes. *Sci. Rep.* **9**, 14556 (2019).

69. Veillard, F., Troxler, L. & Reichhart, J.-M. *Drosophila melanogaster* clip-domain serine proteases: Structure, function and regulation. *Biochimie* **122**, 255–269 (2016).

70. Kanost, M. R. & Jiang, H. B. Clip-domain serine proteases as immune factors in insect hemolymph. *Curr. Opin. Insect Sci.* **11**, 47–55 (2015).

71. Salzberg, S. L. et al. GAGE: A critical evaluation of genome assemblies and assembly algorithms. *Genome Res.* **22**, 557–567 (2012).

72. Waterhouse, R. M. et al. Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. *Science* **316**, 1738–1743 (2007).

73. Xi, Z., Ramirez, J. L. & Dimopoulos, G. The *Aedes aegypti* Toll pathway controls dengue virus infection. *PLoS Pathog.* **4**, e1000098 (2008).

74. Jupatanakul, N. et al. Engineered *Aedes aegypti* JAK/STAT pathway-mediated immunity to dengue virus. *PLoS Negl. Trop. Dis.* **11**, e0005187 (2017).

75. Ertürk-Hasdemir, D. et al. Two roles for the Drosophila IKK complex in the activation of Relish and the induction of antimicrobial peptide genes. *PNAS* **106**, 9779–9784 (2009).

76. Christenbury, J. G. et al. A method for full genome sequencing of all four serotypes of the dengue virus. *J. Virol. Methods* **169**, 202–206 (2010).

77. Manokaran, G. et al. Dengue subgenomic RNA binds TRIM25 to inhibit interferon expression for epidemiological fitness. *Science* **350**, 217–221 (2015).
78. Werren, J. H. & Windsor, D. M. Wolbachia infection frequencies in insects: evidence of a global equilibrium?. *Proc. R. Soc. London Ser. B Biol. Sci.* **267**, 1277-1285 (2000).

79. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).

80. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **9**, 357–359 (2012).

81. Dobin, A. et al. STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics* **29**, 15–21 (2013).

82. Anders, S., Pyl, P. T. & Huber, W. HTSeq—a Python framework to work with high-throughput sequencing data. *Bioinformatics* **31**, 166–169 (2015).

83. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **15**, 550 (2014).

84. Giraldo-Calderón, G. I. et al. VectorBase: An updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. *Nucleic Acids Res.* **43**, D707–713 (2015).

85. Gramates, L. S. et al. FlyBase at 25: Looking to the future. *Nucleic Acids Res.* **45**, D663–D671 (2017).

86. Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. *Nat. Biotechnol* **29**, 644–652 (2011).

87. Birol, I. et al. De novo transcriptome assembly with ABySS. *Bioinformatics* **25**, 2872–2877 (2009).

88. Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. *Bioinformatics* **28**, 3150–3152 (2012).

89. Li, W. & Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. *Bioinformatics* **22**, 1658–1659 (2006).

90. Camacho, C. et al. BLAST+: Architecture and applications. *BMC Bioinform.* **10**, 421 (2009).

91. Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. *BMC Bioinform.* **12**, 323 (2011).

92. Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods* **12**, 59–60 (2015).

## Acknowledgements

We thank the Genome Biology Facility at Duke-NUS Medical School for the RNA sequencing. Support for this research came from a Ministry of Education (Singapore) Tier3 grant (MOE2015-T3-1-003) to RMK and JP, a National Medical Research Council (Singapore) ZRRF grant (ZRRF/007/2017) to JP, an Agence Nationale de la Recherche (France) grant (ANR-20-CE15-0006) to JP, a Singapore Infectious Disease Initiative grant (SIDI/2017/008) to IHM, and the Duke-NUS Medical School Signature Research Programme funded by the Agency for Science Technology and Research (A\*Star Singapore).

## Author contributions

Conceptualization: I.H.M. and J.P. Data curation: C.M.M. Formal analysis: C.M.M., A.C., D.M., R.M.K., I.H.M. and J.P. Funding acquisition: R.M.K., I.H.M. and J.P. Investigation: C.M.M., A.C., D.H.W.L. and M.C. Project administration: I.H.M. and J.P. Software: C.M.M. Supervision: I.H.M. and J.P. Visualization: C.M.M. and J.P. Writing—original draft: C.M.M. and J.P. Writing—review and editing: C.M.M., A.C., D.M., R.M.K., I.H.M. and J.P.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-38354-9.

**Correspondence** and requests for materials should be addressed to J.P.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023