# Transgenic refractory Aedes aegypti lines are resistant to multiple serotypes of dengue virus

Wei-Liang Liu, Chia-Wei Hsu, Shih-Peng Chan, Pei-Shi Yen, Matthew P. Su, Jian-Chiuan Li, Hsing-Han Li, Lie Cheng, Cheng-Kang Tang, Shih-Hsun Ko, et al.

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

To cite this version:

Wei-Liang Liu, Chia-Wei Hsu, Shih-Peng Chan, Pei-Shi Yen, Matthew P. Su, et al. Transgenic refractory Aedes aegypti lines are resistant to multiple serotypes of dengue virus. Scientific Reports, 2021, 11 (1), pp.23865. 10.1038/s41598-021-03229-4 . pasteur-03554699

HAL Id: pasteur-03554699

https://pasteur.hal.science/pasteur-03554699

Submitted on 3 Feb 2022

**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>

Check for updates

OPEN

# Transgenic refractory *Aedes aegypti* lines are resistant to multiple serotypes of dengue virus

Wei-Liang Liu<fnref n="1" />,¹¹, Chia-Wei Hsu²,¹¹, Shih-Peng Chan³, Pei-Shi Yen⁶, Matthew P. Su⁵,⁶, Jian-Chiuan Li⁷, Hsing-Han Li⁷, Lie Cheng<fnref n="1" />, Cheng-Kang Tang⁷, Shih-Hsun Ko<fnref n="1" />, Huai-Kuang Tsai⁸, Zing Tsung-Yeh Tsai⁸, Omar S. Akbari⁹, Anna-Bella Failloux⁴ & Chun-Hong Chen<fnref n="1" />,²,⁷,¹⁰

The areas where dengue virus (DENV) is endemic have expanded rapidly, driven in part by the global spread of *Aedes* species, which act as disease vectors. DENV replicates in the mosquito midgut and is disseminated to the mosquito's salivary glands for amplification. Thus, blocking virus infection or replication in the tissues of the mosquito may be a viable strategy for reducing the incidence of DENV transmission to humans. Here we used the mariner *Mos1* transposase to create an *Aedes aegypti* line that expresses virus-specific miRNA hairpins capable of blocking DENV replication. These microRNA are driven by the blood-meal-inducible carboxypeptidase A promoter or by the polyubiquitin promoter. The transgenic mosquitoes exhibited significantly lower infection rates and viral titers for most DENV serotypes 7 days after receiving an infectious blood meal. The treatment was also effective at day 14 post infection after a second blood meal had been administered. In viral transmission assay, we found there was significantly reduced transmission in these lines. These transgenic mosquitoes were effective in silencing most of the DENV genome; such an approach may be employed to control a dengue fever epidemic.

Dengue fever is a major mosquito-borne viral disease affecting some parts of the world. The World Health Organization has estimated that 50 million cases of dengue virus (DENV) infection occur worldwide per year, and that over 2.5 billion people are now at risk of infection<fnref n="1" />⁻³. Dengue fever is caused by DENV, with four closely related yet antigenically distinct serotypes (DENV-1, DENV-2, DENV-3, and DENV-4)⁴,⁵. This serotype diversity hampers the prevention of dengue fever and creates difficulties for producing a single effective vaccine. The four serotypes differ from each other by 25-40% at the amino acid level. Moreover, the DENV-1 to -4 serotypes comprise multiple genotypes. This diversity makes it particularly difficult to develop effective anti-DENV treatment methods.

DENV is transmitted by female mosquitoes mainly of the *Aedes aegypti* and *Ae. albopictus* species. Once mosquito acquired the DENV, virus would infect midgut in 3-5 days, disseminate into circulatory system to reach salivary glands in 7-10 days, and transmit virus to next host in 10-14 days⁶⁻⁸. Therefore, blocking virus acquisition or replication in the tissues of the mosquito may be a strategy for reducing the incidence of DENV transmission to humans. Several approaches can be employed to control mosquito-borne diseases, for example, using chemical pesticides, releasing *Wolbachia*-infected mosquitoes⁹, and using genetic engineering for mosquito control. Nevertheless, blocking virus acquisition or replication in the tissues of the mosquito may be more efficient, and could be a crucial strategy for reducing the incidence of DENV transmission to humans. Studies

<footnote n="1">¹National Mosquito-Borne Diseases Control Research Center, National Health Research Institutes, Miaoli, Taiwan. ²Institutes of Molecular and Cellular Biology, National Taiwan University, Taipei, Taiwan. ³Graduate Institute of Microbiology, College of Medicine, National Taiwan University, Taipei, Taiwan. ⁴Unit of Arboviruses and Insect Vectors, Department of Virology, Institute Pasteur, Paris, France. ⁵Department of Biological Science, Nagoya University, Nagoya, Japan. ⁶Institute for Advanced Research, Nagoya University, Nagoya, Japan. ⁷National Institute of Infectious Diseases and Vaccinology, National Health Research Institutes, Miaoli, Taiwan. ⁸Institute of Information Science, Academia Sinica, Taipei, Taiwan. ⁹Section of Cell and Developmental Biology, Division of Biological Sciences, University of California, San Diego, La Jolla, CA 92093, USA. ¹⁰National Institute of Infectious Diseases and Vaccinology, National Mosquito-Borne Diseases Control Research Center, National Health Research Institutes, Zhunan, Taiwan. ¹¹These authors contributed equally: Wei-Liang Liu and Chia-Wei Hsu. ²email: chunhong@gmail.com</footnote>
have revealed that a gene drive to carry a refractory gene into a population is now feasible in some cases¹⁰⁻¹². Driving the anti-DENV gene into the mosquito population is thus also an option for DENV control. Such effector genes can help prevent the relevant disease, and gene-drive technology can be used to drive the introgression of these genes into a population of *Ae. aegypti*, replacing the disease-competent population with a refractory one.

In recent years, more promising genetic-based strategies have been developed to control mosquito populations¹³⁻¹⁶. Among them, RNA interference (RNAi) has been shown to be effective against mosquito-borne diseases, including DENV, chikungunya virus (CHIKV), and Zika virus (ZIKV)¹⁷⁻¹⁹. RNAi can be triggered by conventional inverted repeats or synthetic microRNAs (miRNAs)²⁰,²¹. miRNAs are 20 to 24 nucleotides (nt) small RNAs that regulates many functions, such as gene expression, post-transcriptionally through translational repression and mRNA cleavage. In the miRNA processing pathway, the miRNA precursors (pri-miRNA) are processed into ~70 nt hairpins by Drosha and folded into stem-loop structures in the nucleus. They are subsequently exported into the cytoplasm by Exportin5. In the cytoplasm, the stem-loop structures are cleaved into ~22 nt miRNA duplexes by Dicer-1, which are then loaded into either Ago-1 or Ago-2 proteins in miRNA-induced silencing complexes (RISCs) according to their different structural properties²²⁻²⁴. By recognizing a complementary sequence contained within a target RNA, RISCs execute silencing through RNA degradation, translational inhibition, or both methods²⁵,²⁶. To construct the artificial antiviral synthetic miRNA, the consensus sequence region of the virus was selected as the target site for the antiviral miRNA²⁷,²⁸. The corresponding miRNAs were designed and cloned in tandem to make a miRNA-based RNAi cassette. Using a miRNA-based genetically engineered model to produce a polycistronic artificial miRNA permits the simultaneous expression of multiple synthetic miRNAs that mediate the degradation of various regions of one target RNA sequence or multiple RNA sequences from different DENV genes.

Transgenic mosquitoes also require functional promoters to drive expression of the synthetic miRNAs following a blood meal (BM). Active promoters have been identified in the cells of the *Ae. aegypti* mosquito. The midgut is reportedly the primary tissue for DENV replication. Both the polyubiquitin (*PUB*) promoter²⁷ and the blood-meal-inducible carboxypeptidase A (*CPA*) promoter²⁸ are active in the midgut epithelial cells of *Ae. aegypti*. Moreover, multiple blood feedings in female *Ae. aegypti* during a single gonotrophic cycle are associated with variation in the pattern of DENV transmission²⁹. Therefore, whether multiple blood feeds may lead to increases in density or associated changes in DENV blocking ability of refractory transgenic mosquitoes should be investigated.

Here, we demonstrate anti-DENV miRNAs under the control of *AeCPA* and *AePUB* promoters that drive the expression of a polycistronic cluster of eight synthetic miRNAs; these miRNAs were engineered to target conserved sequences in the four DENV genomes. Our experiments revealed that our transgenic mosquitoes had an effective antiviral capacity. The antiviral capability of transgenic lines may present new opportunities for DENV control.

## Results

**Designing the anti-DENV miRNA and targeting sites on the DENV genome.** To identify the most highly conserved regions of each serotype of the DENV genome for constructing the silencing cassette, we analyzed the viral genome collection of the Viral Bioinformatic Resources Center (VBRC; http://athena.bioc.uvic.ca/), and the serotype sequences were aligned and blasted separately using CLC bio-Genomics software (QIAGEN Digital Insights). The miRNA targeting sequences were designed based on DENV consensus sequences. Eight regions of DENV-1 to DENV-4 flanked by the 5' untranslated region (5'UTR), capsid, and 3'UTR were selected as the miRNA targeting sites (Fig. 1a). The sequence coverage of the anti-DENV miRNAs relative to the viruses used to generate the consensus sequence was 66.9–100% (Fig. 1b).

All antiviral miRNA clusters were constructed in an effort to place them under the *Ae. aegypti* polyubiquitin (*AePUB*) or *Ae. aegypti* carboxypeptidase A (*AeCPA*) gene promoters to elicit either constitutive or blood-meal-inducible midgut-specific expression of the effector molecules. The polycistronic artificial miRNAs would be expressed simultaneously under the control of the same transgenic expression cassettes; therefore, multiple copies of miRNA specific to one DENV serotype or to all four serotypes could be produced in one transgenic mosquito to induce antiviral RNAi machinery to reduce DENV replication as illustrated in Fig. 1c. We used the mariner *Mos1* transposase³⁰ to generate transgenic *Aedes* mosquitoes that could express anti-DENV microRNAs, and we used a genomic integration system to produce a series of stable hereditary transgenic mosquitoes containing the enhanced green fluorescent protein (EGFP) reporter gene (Supplementary Fig. 1). The EGFP reporter was driven by the 3×P3 promoter expressed in the eyes or anal papillae of transgenic larvae lines (Fig. 1d). Several transgenic lines were selected according to miRNA expression level. One heterozygous line for each construct (*AeCPA-8miR* and *AePUB-8miR*) was retained for further experiments (Supplementary Table 1). The transgenic cassettes in mosquito chromosomes were confirmed by Southern blot analysis using restriction enzymes that have a EcoRI/HindIII cutting site within the transgene line and [P³²]-labeled probes complementary to the EGFP reporter gene. The results indicated that the reporter gene had been integrated into the genomes of the transgenic mosquitoes (Supplementary Fig. 2) and also confirmed that the mosquitoes carried the artificial miRNA gene.

**MicroRNA expression levels in transgenic mosquitoes.** Clustered miRNAs typically exhibit widely varying expression levels³¹⁻³³. To investigate the expression levels of each miRNA in the cluster in each transgenic mosquito line, we evaluated the fold change in miRNA expression in the transgenic mosquitoes using quantitative polymerase chain reaction (qPCR) analyses. The miRNA expression in the midgut and carcass tissue 16 h post blood meal (PBM) was compared with miRNA expression level in the absence of a blood meal (no-PBM) in *AeCPA-8miR* transgenic lines. We found that miR-8-2 (targets DENV-2 and -4) expression had the highest fold increase, with a fold change near 70 relative to no-PBM, and the second highest expression was
<figure><img src="image_4.png" /><figcaption>(b)</figcaption></figure>

<table><thead><tr><th rowspan="2">Generation</th><th rowspan="2">miR</th><th rowspan="2">Target Site</th><th rowspan="2">Position</th><th colspan="4">Targeting Frequency</th></tr><tr><th>Type 1 (651)</th><th>Type 2 (615)</th><th>Type 3 (356)</th><th>Type 4 (44)</th></tr></thead><tbody><tr><td rowspan="7">DENV-8miR</td><td>miR-1</td><td>TAGAGAGCAGATCTCTGATGAA</td><td>5'UTR</td><td>578/651 (88.8%)</td><td>612/615 (99.5%)</td><td>238/356 (66.9%)</td><td>0/44 (0%)</td></tr><tr><td>miR-2</td><td>TATGCTGAAACGCGAGAGAAAC</td><td>Capsid</td><td>0/651 (0%)</td><td>615/615 (100%)</td><td>0/356 (0%)</td><td>44/44 (100%)</td></tr><tr><td>miR-3</td><td>CGTCTTTCAATATGCTGAAACG</td><td>Capsid</td><td>651/651 (100%)</td><td>0/615 (0%)</td><td>0/356 (0%)</td><td>0/44 (0%)</td></tr><tr><td>miR-4</td><td>CGCGAGAAACCGCGTGTCAACT</td><td>Capsid</td><td>642/651 (98.6%)</td><td>0/615 (0%)</td><td>0/356 (0%)</td><td>0/44 (0%)</td></tr><tr><td>miR-5</td><td>CGTCTATCAATATGCTGAAACG</td><td>Capsid</td><td>0/651 (0%)</td><td>0/615 (0%)</td><td>355/356 (99.7%)</td><td>0/44 (0%)</td></tr><tr><td>miR-6</td><td>CGTGAGAAACCGTGTGTCAACT</td><td>Capsid</td><td>0/651 (0%)</td><td>0/615 (0%)</td><td>354/356 (99.4%)</td><td>0/44 (0%)</td></tr><tr><td>miR-7</td><td>AGCAAACCGTGCTGCCTGTAGC</td><td>3'UTR</td><td>491/651 (75.4%)</td><td>0/615 (0%)</td><td>330/356 (92.7%)</td><td>44/44 (100%)</td></tr><tr><td></td><td>miR-8</td><td>CACAGTTGGCGAAGAGATTCTC</td><td>Capsid</td><td>649/651 (99.7%)</td><td>0/615 (0%)</td><td>352/356 (98.9%)</td><td>0/44 (0%)</td></tr></tbody></table>

**Figure 1.** Sequence of the designed antiviral miRNAs was based on DENV consensus genomic sequence. (a) DENV target genes. Through the mariner (Mos1) transposase, ubiquitous and midgut-specific induction promoters drive the expression of a polycistronic cluster of eight synthetic small RNAs that were engineered to target conserved genes in the DENV genome. A genomic integration system was used to produce a series of transgenic mosquitoes that can express antiviral miRNAs. The DENV-8miR plasmid contains eight miRNAs: miR_8-1 targeting the 5'UTR, miR_8-7 targeting the 3'UTR, and miR_8-2 to miR_8-6 and miR_8-8 targeting the capsid, as indicated with colored circles. (b) The binding site coverage and similarity are presented in the table. Five types of miRNAs for DENV-1 (Myanmar strain), two types of miRNAs for DENV-2 (NGC strain), five types of miRNAs for DENV-3 (98TW503 strain), and two types of miRNAs for DENV-4 (H241 strain) were complementary targets to DENV. The viral genome sequences of each serotype were downloaded from the Viral Bioinformatic Resources Center (http://athena.bioc.uvic.ca/) and then aligned using CLC bio-Genomics software version 9.5.3 (CLC Bio, https://www.qiagenbioinformatics.com/). (c) Schematic of anti-DENV miRNA production. Schematic of the anti-DENV transgene, consisting of an *AeCPA/AePUB* promoter that drives expression of the polycistronic cluster of eight synthetic small RNAs that were engineered to target conserved genes in the DENV genome. (d) Screening of transgenic mosquitoes under fluorescence microscopy. The EGFP reporter driven by the 3×P3 promoter was expressed in the eyes or anal papillae of transgenic larvae lines. Examples of EGFP reporter expression in larval and adult stages under natural (left) and blue light (right) conditions. Arrows indicate eyes and asterisks indicate anal papillae regions with EGFP expression.
<figure><img src="image_5.png" /><figcaption>(d)</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 1.** (continued)</figcaption></figure>

noted in miR_8-8 (targets DENV-1 and -3), which increased approximately 50-fold. Expression levels of miR_8-3, -4, -5, and -6 remained modest, whereas those of miR_8-1 and 8-7 were low (Fig. 2a, left panel). The antiviral miRNAs of AeCPA-8miR mosquitoes were detectable with similar results in the carcass, and an increased expression level was observed 16 h PBM (Fig. 2a, right panel). As the AePUB-8miR transgenic line miR_8-2 was by far the most highly expressed miRNA, with a fold change about 40 times relative to no-PBM. The expression levels of miR_8-6, -7 and -8 were modest, whereas those of miR_8-1, -3, -4, and -5 were relatively low (Fig. 2b left panel). In the carcass tissue, the antiviral miRNAs were also detected (Fig. 2b, right panel). We interpret these data to indicate that the expression of the antiviral miRNAs in the midgut and carcass remained inducible after receiving a blood meal.
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Detection of microRNA expression levels in transgenic mosquitoes. Artificial miRNAs expressed in the midgut and carcass of female transgenic mosquitoes from (a) *AeCPA-8miR* and (b) *AePUB-8miR* lines detected through qPCR. Total RNA was isolated from mosquito midguts dissected at 0 and 16 h post blood meal (PBM). RT and qPCR were conducted as described in the “Methods” section. Data are presented as relative expression levels compared with aae-miR-100. Statistical comparisons of expression levels before and after taking a blood meal for each miRNA were made using t-tests (with adjustments for multiple comparisons). Each bar represents as mean ± SEM. Significant *p* values are indicated by asterisks: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.</figcaption></figure>

**Transgenic mosquitoes suppressed DENV titer and infection rate.** To test the antiviral capacity of the transgenic mosquitoes, we challenged the aforementioned two lines, *AeCPA-8miR* and *AePUB-8miR*, with all four DENV serotypes (DENV-1 of the Myanmar strain, DENV-2 of the NGC strain, DENV-3 of the 98TW503 strain, and DENV-4 of the H241 strain). We conducted the test after an incubation period of 7 or 14 days after oral feeding with a virus-infected BM. We assessed the infection rate (IR) and viral titer, and we defined IR as the number of DENV-infected mosquitoes divided by the total number of mosquitoes and measured the viral titer per mosquito. Seven days after a BM containing 10⁷ plaque-forming units (PFU)/mL of virus, the *AeCPA-8miR* line exhibited diverse antiviral capacity against the four strains (Fig. 3a–d). In the control sample (Orlando strain mosquitoes), the viral titer from the midgut of DENV-2 was approximately 1.6 × 10⁶ PFU/mL, but in *AeCPA-8miR* the average titer was approximately 10² PFU/mL (Fig. 3b). For DENV-3, the average viral titer in the control sample was approximately 3.3 × 10⁴ PFU/mL, compared with approximately 2.8 × 10³ PFU/mL for the *AeCPA-8miR* line (Fig. 3c). Similarly, *AeCPA-8miR* was still effective against DENV-4 virus (10² PFU/mL vs. 3.8 × 10³ PFU/mL for the control, *p* < 0.01; Fig. 3d). However, for DENV-1, no difference was noted in the viral titer between the transgenic and control mosquitoes (Fig. 3a). With respect to IR (Fig. 3e–h), the *AeCPA-8miR* line exhibited a substantially reduced IR relative to the control line for DENV-2 (20% vs. 55%), DENV-3 (22% vs. 47%), and DENV-4 (33% vs. 55%), but not for DENV-1 (33% vs. 43%). This suggests that the *AeCPA* promoter driving a cassette containing the eight anti-DENV miRNAs can suppress three of the four serotypes but not DENV-1.
<figure><img src="image_8.png" /></figure>

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

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

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

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

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

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

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

**Figure 3.** DENV challenge experiments on *AeCPA-8miR* mosquitoes at day 7. We conducted DENV challenge experiments on the transgenic *Aedes* mosquito line *AeCPA-8miR*. Seven days after a blood meal containing DENV, we collected the DENV-infected mosquitoes' midgut and performed plaque assays using BHK21 cells. We performed the virus challenges using (a) DENV-1 of the Myanmar strain, (b) DENV-2 of the NGC strain, (c) DENV-3 of the 98TW503 strain, and (d) DENV-4 of the H241 strain on the *AeCPA-8miR* and control mosquitoes (Orlando strain). For each mosquito line, the mean viral titer is plotted, and the standard error of the mean is indicated. Uninfected mosquito samples were not plotted or included when determining the mean or standard error of the mean (n = 32–40). The statistical test is the Mann–Whitney rank sum test of GraphPad Prism (version 5.0, http://www.graphpad.com) to calculate and analyze the difference in the virus titers. (e–h) The infection rates (IRs) of the mosquito lines are indicated as percentages. IR was defined as the number of virus positive midgut samples out of the total number tested. A *t*-test was used to analyze the difference in the mosquito IRs. The error bars correspond to the confidence intervals (95%). All data are represented as mean ± SEM. Significant *p* values are indicated by asterisks: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.
In contrast to AeCPA, which is induced by a BM and specifically expressed in the midgut, PUB is expressed
in all tissues. The AePUB-8miR strain had a lower virus titer for DENV-2 after infection (2.1 × 10⁴ PFU/mL vs.
9.7 × 10⁵ PFU/mL for the control, *p* < 0.05; Fig. 4b). However, its viral titers for DENV-1, DENV-3, and DENV-4
were similar to those of the control line (Fig. 4a,c,d). The IR for the AePUB-8miR line was lower than that of
the control for DENV-2 and DENV-3 (13% vs. 53% and 15% vs. 47%, respectively, *p* < 0.001; Fig. 4f,g). It had a
slight suppressive effect against and had an even higher IR for DENV-1 (78% vs. 43%, *p* < 0.01; Fig. 4e), whereas
no effect was noted for DENV-4 (35% vs. 55%; Fig. 4h).

For the longer incubation time, 14 days post infection (dpi), our analysis revealed that the viral titers of DENV-2 and DENV-4 in the AeCPA-8miR line were significantly lower than in the control groups (DENV-2: 3.2 × 10⁴ PFU/mL vs. 4.6 × 10⁵ PFU/mL, *p* < 0.01; DENV-4: 1.3 × 10⁵ PFU/mL vs. 1.7 × 10⁴ PFU/mL; Supplementary Fig. 3b,d). However, for DENV-1 and DENV-3, the viral titers of the AeCPA-8miR line showed no differences from those of the control mosquitoes (Supplementary Fig. 3a,c). In terms of IR, AeCPA-8miR mosquitoes have lower IR after DENV-2 (32% vs. 65%) or DENV-4 (45% vs. 72%) infection, and higher IR after DENV-3 infection (45% vs. 97%) but no significant inhibition of DENV-1 (Supplementary Fig. 3e-h). In AePUB-8miR transgenic mosquitoes, the viral titer of DENV-2 was also significantly lower than in the control group (7.5 × 10³ vs. 9.3 × 10⁵ PFU/mL, *p* < 0.05, Supplementary Fig. 4b), but significant differences were not observed for the other three serotypes of DENV. IR analysis indicated that DENV-3 infection in AePUB-8miR mosquitoes was slightly lower than in the control line (25% vs. 45%, Supplementary Fig. 4g), but no effect was observed for the other three serotypes. Overall, our results indicate that the transgenic mosquitoes using either promoter had the ability to defend against DENV replication and infection.

**Repeat BM affected infection rate and viral titers.** In the field, female mosquitoes may obtain blood meals multiple times, and this behavior increases the risk of DENV transmission and outbreaks. Since the AeCPA promoter is activated by blood meal uptake, we next investigated the defending efficiency of the eight miRNAs after a 2nd blood meal. We first infected mosquitoes with 10⁷ PFU/mL DENV-4 at the first BM and then incubated the successfully blood-fed mosquitoes for seven days. After incubation, mosquitoes consumed another BM which was virus-free. These two BMs activated the AeCPA promoter to express the 8 miRNAs to silence DENV genomes (Fig. 5a). qPCR results show that miRNAs were expressed after both BMs in midgut and carcass tissue (Supplementary Fig. 5). The antiviral capacity test revealed that, compared with the transgenic mosquitoes that did not have a second BM, the mosquitoes that received a second BM exhibited a slight but significant reduction in viral titer (5.8 × 10⁴ PFU/mL vs. 7.1 × 10³ PFU/mL, *p* < 0.01; Fig. 5b) and a 22% lower IR (from 62.5 to 40.0%, *p* < 0.05; Fig. 5c). These results indicate that the second blood meal potentially increases the blocking efficiency of DENV in refractory transgenic mosquitoes.

**Virus transmission efficiency in transgenic mosquitoes.** To evaluate the transmission ability of these transgenic lines, the saliva of orally infected mosquitoes was examined. Mosquitoes were infected orally with DENV-1 to -4 and then incubated for 7 or 14 days, after which the salivary glands and saliva of these mosquitoes were collected to test the virus titers and transmission efficiency (Fig. 6). In the salivary glands of the AeCPA-8miR line, DENV-2 and DENV-4 viral titers were significantly lower both 7 and 14 dpi when compared to controls (Fig. 6d,j). However, for DENV-1 and DENV-3, no difference was noted in the viral titer between the transgenic and control mosquitoes (Fig. 6a,g). For the AePUB-8miR line, we found reductions in only DENV-2 at 14 dpi (Fig. 6d), but showed no significant changes for any of the other three serotypes of DENV. After testing for potential suppressive effects in the salivary glands, we next tested viral titer in saliva. In the saliva, we found that AeCPA-8miR showed significantly reduced viral titers for DENV-2 and DENV-4 at 7 and 14 dpi (Fig. 6e,k). On the other hand, AePUB-8miR was shown to reduce the virus titer of DENV-2 at 14 dpi, but there was no significant difference in virus titer at 7 dpi (Fig. 6e). For DENV-1 and DENV-3, neither AeCPA-8miR and AePUB-8miR line showed a significant suppressive effect (Fig. 6b,h). At 7 dpi, the DENV-2 transmission efficiency was, on average, 19.4% ± 14.1% and 38.9% ± 22.3% in the AeCPA-8miR and AePUB-8miR lines, respectively, and the control (Orlando) mosquitoes had an average transmission efficiency of 53.3% ± 14.4% (Fig. 6f). At 14 dpi, we observed that the transgenic lines had a significantly lower transmission efficiency than the control mosquitoes (Orlando: 54.2% ± 8.3%, AeCPA-8miR: 8.3% ± 11.8%, AePUB-8miR: 34.9% ± 5.9%) (Fig. 6f). We also found similar results in the transmission efficiency for the AeCPA-8miR line was lower than that of the control for DENV-4 at 7 and 14 dpi (22.2% ± 15.4% vs. 52.2% ± 11.9% and 8.3% ± 11.8% vs. 51.9% ± 11.5%, respectively, *p* < 0.01; Fig. 6l). However, for DENV-1 and DENV-3, the transmission efficiency of the AeCPA-8miR and AePUB-8miR line showed no differences from those of the control mosquitoes (Fig. 6c,i). In these transgenic lines, the transmission efficiency of DENV-2 and DENV-4 is significantly affected by AeCPA-8 miRNA, while that of DENV-1 and DENV-3 is not affected.

**Discussion**

Dengue fever is one of the most common and widespread mosquito-borne diseases in the world. It contains four antigenically distinct serotypes (DENV-1 to -4), any one of which could be the major serotype in an epidemic. Different approaches for preventing mosquitoes from being infected with DENV had been demonstrate³⁻³⁶. However, each approach only targeted a specific serotype, which makes these previous methods insufficient and impractical. In this study, we generated two transgenic mosquito lines that expressed antiviral microRNAs, which were designed to target the genomes of the four DENV serotypes under the control of the Ae. aegypti PUB and CPA promoters. After viral challenge, the viral titer and IR suggested a significant suppression of three serotypes of DENV infection in AeCPA-8miRNA transgenic lines and one serotype of DENV infection in AePUB-8miRNA transgenic lines at day 7 PBM compared to wild-type (Orlando strain) mosquitoes. After a second blood meal,
<figure><img src="image_16.png" /><figcaption>(b)</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>(c)</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>(d)</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>(e)</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>(f)</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>(g)</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>(h)</figcaption></figure>

<figure><img src="image_23.png" /></figure>
►**Figure 4.** DENV challenge experiments on *AePUB*-8miR mosquitoes at day 7. We conducted DENV challenge experiments on the transgenic *Aedes* mosquito line *AePUB*-8miR. Seven days after a blood meal containing DENV, we collected the DENV-infected mosquitoes and performed plaque assays using BHK21 cells. We performed virus challenges using (a) DENV-1 of the Myanmar strain, (b) DENV-2 of the NGC strain, (c) DENV-3 of the 98TW503 strain, and (d) DENV-4 of the H241 strain on the *AePUB*-8miR and control mosquitoes (Orlando strain). For each mosquito line, the mean viral titer is plotted, and the standard error of the mean is indicated. Uninfected mosquito samples were not plotted or included when determining the mean or standard error of the mean (n = 32–40). The statistical test is the Mann–Whitney rank sum test was of GraphPad Prism (version 5.0, http://www.graphpad.com) to calculate and analyze the difference in the virus titers. (e–h) The infection rates (IRs) of the mosquito lines are indicated as percentages. IR was defined as the number of virus positive midgut samples divided by the total number tested. A t-test was used to analyze the difference in the mosquito IRs. The error bars correspond to the confidence intervals (95%). All data are represented as mean ± SEM. Significant *p* values are indicated by asterisks: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.

the *AeCPA*-8miRNA line retained its effective antiviral activity at day 14 post infection. The virus transmission efficiencies in the *AeCPA*-8miRNA line were significantly reduced compared with both the *AePUB*-8miRNA and the control mosquitoes at 14 days PBM. In summary, our study reveals that miRNA-based genetically engineered mosquitoes had reduced viral titers, IRs, and transmission efficiency for the most of the four DENV serotypes.

Both of the transgenic mosquito lines that were generated here had one copy of the DENV miRNA cassette (heterogeneous, data not shown). Our results demonstrate that the heterogeneous DENV-8 miRNA transgenic mosquito lines were able to obtain refractoriness to DENV infection and transmission, with the result that the infection rates of some of the four serotypes (types 2, 3, and 4 in *AeCPA*-8miR; 2 and 3 in *AePUB*-8miR; Figs. 3 and 4) of DENV were significantly lower than the infection rates of the corresponding serotype in the Orlando strain. Furthermore, we found that the DENV-4 infection in the salivary gland of the *AeCPA*-8miR line was significantly lower than in the Orlando strain, resulting in the viral load in the saliva being reduced by 93% (from 36.1 to 2.7%) when compared to Orlando mosquitoes. Similar results have also been found in previous research³², which demonstrated that a heterogeneous line of miRNA transgenic mosquitoes were able to reduce ZIKV infection and transmission. Lower infection and transmission rates in mosquitoes represent a lower risk of virus dissemination, and this can be achieved in a heterogeneous background. We believe modifying the miRNA sequences or promoter strength may provide a higher suppression ability against DENV infection and transmission.

After mating, female mosquitoes are able to ingest multiple blood meals³⁷,³⁸. Once infected by DENV, and after the virus has entered the salivary glands, females are able to transmit DENV during subsequent biting events for the rest of their lives³⁹. The CPA promoter is not a continuous expression promoter; its activity is triggered by a blood meal. From the single oral DENV infection/blood meal results (Fig. 3), we found that viral load and infection rate of DENV in the CPA promoter transgenic line were significantly lower than those of Orlando mosquitoes. In the *AeCPA*-8miR line, after the blood meal is digested, antiviral protection is likely to be lost because the miRNA expression will cease. Virus infection in mosquito body will resume a natural state. We then stimulated the expression of the miRNAs again (Supplementary Fig. 5) by feeding another blood meal to these mosquitoes seven days after first blood meal. In these mosquitoes, we found that the viral load and infection rate were significantly lower than in those that were only fed a blood meal once (Fig. 5). These results demonstrated that whether these miRNAs were expressed during or after DENV infection, they were capable of suppressing DENV replication.

In our result, the expression level of each miRNA of DENV-8miR cluster were different (Fig. 2). In a cell, several miRNAs existed as a cluster was not hard to see, each cluster was under controlled by a promoter, however, the expression levels of these miRNAs are different³¹,⁴⁰⁻⁴², even though they are all produced from the same transcript and driven by one promoter. There are some possible reasons to make this variation of expression, the architecture of miRNAs, including the dimensional structure of the hairpin, the spacing sequence and the length of spacing are all crucial importance for the knockdown efficiency. In the miRNA processing, the guide strand is the one that is used in the RISC for RNA interference⁴³,⁴⁴. However, how the RISC chooses the guide strand (the functional strand for silencing) from the individual miRNA precursor was dependent on their individual sequence. Based on these results, although the knockdown efficacy in this study did not achieve to the same level for all four different genomes of DENV. In the future, we may optimize biophysical properties and/or sequence of the spacer, spacer, and select the more efficient target sequence to optimize knockdown efficiency for all DENV genome.

The results of virus challenge and plaque assay experiments revealed that DENV infection was reduced to around 50% in the transgenic mosquitoes. Both of the transgenic lines showed effective antiviral activity against DENV-2 and DENV-4 infection either at 7- or 14-days post infection, but not to DENV-1 infection (Figs. 3 and 4 and Supplementary Figs. S3 and S4). This may be due to the high expression level of miR_8-2 in both transgenic lines (Fig. 2). The sequence of miR_8-2 is designed to target genomic sequences present in both DENV-2 and DENV-4, resulting in the effective suppression of DENV replication. From the viral load results and DENV-1 infection rates, we can conclude that the function of miR_8-3, miR_8-4, and miR_8-8 may not be sufficient to provide antiviral activity, even though the expression level of miR_8-8 was the second highest among the eight miRNAs. miR_8-8 targets both DENV-1 and DENV-3, however DENV-3 is also targeted by two additional miRNAs, miR_8-5 and miR_8-6, which do not have targets within the DENV-1 sequence. Antiviral activity against DENV-3 in the transgenic lines may have contributions from these two miRNAs. miR_8-1 and miR_8-7 were designed to target multiple serotypes of DENV, including DENV-1, however, their expression was the lowest of these eight miRNAs, which may have restricted their efficiency in anti-DENV activity. The actual conditions of
<figure><img src="image_24.png" /></figure>

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

<figure><img src="image_26.png" /><figcaption>**Figure 5.** DENV challenge experiments using *AeCPA-8miR* mosquitoes administered with a second blood meal. (a) Schematic of the study design. A mixture containing DENV-4 supernatant and mouse blood was used to feed the *CPA* > 8miR transgenic mosquitoes at the first BM. Seven days after the blood feeding, these transgenic mosquitoes received another feeding with mouse blood not containing the virus; the DENV-infected mosquitoes were collected at day 14, and a plaque assay using BHK21 cells was performed to compare these mosquitoes with the *AeCPA-8miR* mosquitoes that received only one BM. (b) Viral titers of the two groups of mosquitoes. Viral titer of repeat BM versus day 14 in *AeCPA-8miR* mosquitoes. For each mosquito line, the mean number of viral titers is plotted, and the standard error of the mean is indicated. Uninfected mosquito samples were not plotted or included when determining the mean or standard error of the mean (n = 40). The statistical test is the Mann-Whitney rank sum test was of GraphPad Prism (version 5.0, http://www graphpad com) to calculate and analyze the difference in the virus titers. (c) The infection rates of the two groups of mosquitoes. Infection rate was defined as the number of positive midgut samples divided by the total number tested. The infection ratio of mosquito lines is indicated as percentages. A t-test was used to analyze the difference in the mosquito IRs. The error bars correspond to the confidence intervals (95%). All data are represented as mean ± SEM. Significant *p* values are indicated by asterisks: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.</figcaption></figure>
<figure><img src="image_27.png" /></figure>

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

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

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

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

<figure><img src="image_32.png" /><figcaption>**Figure 6.** Transgenic mosquitoes have reduced virus transmission potential. At 7- and 14-dpi with DENV-1 to -4, the salivary glands and the saliva were collected in female transgenic and control (Orlando strain) mosquitoes. Viral suppression was determined by plaque assay. (a,d,g,j) In salivary glands. (b,e,h,k) In saliva. Saliva samples were collected after 30 min in a P20 tip containing 5 µL of FBS and then expelled into 45 µL of L-15 media for analysis. For each mosquito line, the mean viral titer is plotted, and the standard error of the mean is indicated. Uninfected mosquito samples were not plotted or included when determining the mean or standard error of the mean. Each sample corresponds to three replicates (3 × 12 mosquitoes). The statistical test is the Mann-Whitney rank sum test was of GraphPad Prism (version 5.0, http://www.graphpad.com) to calculate and analyze the difference in the virus titers. (c,f,i,l) The transmission efficiencies of infection of the mosquito lines are indicated as percentages, which are defined as the number of positive saliva samples (i.e., successful transmission events) divided by the number of tested samples. A t-test was used to analyze potential differences in the mosquito transmission efficiency assay. The error bars correspond to the 95% confidence intervals. All data are represented as mean ± SEM. Significant *p* values are indicated by asterisks: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.</figcaption></figure>
<figure><img src="image_33.png" /></figure>

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

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

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

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

<figure><img src="image_38.png" /><figcaption>**Figure 6.** (continued)</figcaption></figure>

how these eight miRNAs function in the suppression of viral load and infection rate in DENV infection are far more complicated and may include unexpected interactions, for example, mutations generated in the DENV genome sequence, the targeting positions of these eight miRNAs within the DENV genome, interactions between exogenous and endogenous miRNAs, secondary structures or protein-RNA interactions, or even transcription or translation regulation⁴⁵ by these eight miRNAs. Therefore, to improve the efficiency of suppression of DENV infection in transgenic mosquitoes, future experiments should focus on optimizing the sequences of anti-DENV miRNAs, as well as their expression and processing.

Another factor that affects DENV replication is temperature, which may also affect the efficacy of RNA silencing⁴⁶,⁴⁷. Previous studies have revealed that the extrinsic DENV-2 incubation period in *Ae. aegypti*, *Ae. albopictus*, and *Culex pipiens* was reduced at higher temperatures⁴⁸⁻⁵⁰. In other words, a short extrinsic incubation period may increase the capacity of a vector to transmit the virus. In the present context of global warming, the
DENV may have a shorter extrinsic incubation period in mosquitoes. In future studies, virus challenge experiments should therefore be performed at higher temperatures to elucidate the antiviral capacity of transgenic mosquitoes more comprehensively.

This study demonstrated that the expression of DENV sequence specific miRNAs in transgenic *Ae. aegypti* decreased DENV infectivity. These eight miRNAs were designed to target genome sequences in all four serotypes of DENV separately. Their expression improved the RNAi pathway's effectiveness in the fight against DENV and resulted in an effectively reduced viral load and infection rate in transgenic mosquitoes. The transmission efficiency of DENV in transgenic mosquitoes was decreased as well. This strategy may be useful for combating the spread of DENV or other arboviruses.

## Methods

**Mosquito rearing.** *Aedes aegypti* (Orlando strain, which was from the University of Maryland Insect Transformation Facility) were reared in an insectary under the following conditions: 28 °C, 70%-80% relative humidity, and a 12-h dark-light regime. Approximately 100 adults were pooled and maintained in a 30 × 30 × 30 cm acrylic cage with 10% sucrose as a source of dietary sugar. Adult mosquitoes shifted to feeding on blood from a mouse donor (BALB/c; BioLASCO Taiwan Co., Ltd., Taipei, Taiwan). Three days after the BM, the mosquitoes were allowed to lay eggs on wet filter paper, which was then dried and stored.

**miRNA-based RNAi design.** In order to design artificial miRNAs targeting the viral genomic RNA of the four serotypes (1, 2, 3, and 4) of DENV, the viral genome sequences of each serotype were downloaded from the Viral Bioinformatic Resources Center (http://athena.bioc.uvic.ca/) and then aligned using CLC bio-Genomics software version 9.5.3 (CLC Bio, QIAGEN, https://www.qiagenbioinformatics.com/). The downloaded viral genome sequences included 651 sequences of serotype 1, 615 sequences of serotype 2, 356 sequences of serotype 3, and 44 sequences of serotype 4. Each serotype was aligned separately. The alignment results revealed several highly conserved regions located within the 5'UTR, capsid, or 3'UTR in the viral genomes. We selected miRNA sequences from these highly conserved regions based on their GC content and the length of the conserved region. These miRNA sequences were then verified in their sequence identity to each serotype via alignment with genome sequences using UltraEdit (IDM Computer Solutions). The chosen miRNA sequences were then cloned into the miR6.1 backbone, which was derived from *Drosophila melanogaster*<sup>51</sup>, to generate the miR-based RNAi cassette for silencing DENV in mosquitoes. The sequences of these eight miRNAs are listed in Supplementary Table 2.

**Plasmid construction.** We first generated tandem miRNA coding DNA fragments. Double-stranded DNA fragments of eight DENV miRNAs were composed using paired oligonucleotides (Supplementary Table 3) through annealing and extension, followed by another PCR amplification with a primer set: miR6.1_5' EcoRI/BglII and miR6.1_3' XhoI/BamHI, to produce miRNA units that contained restriction enzyme cutting sites. These eight miRNA units were assembled to become a tandem miRNA coding DNA fragment using BglII and BamHI sites and then subcloned into the EcoRI/XhoI sites of pMOS1_AePUB-CLEC18A-2xHA_3xP3-eGFP<sup>40</sup> to create pMOS1_nanos-DENV-8miR_3xP3-EGFP constructs. To finish the construction, promoter sequences of AeCPA and AePUB were amplified from the genomic DNA of *Ae. aegypti* using a CPA primer set (AeCPA-F and AeCPA-R) and a PUB primer set (AePUB-F and AePUB-R). The amplified DNA products were then cloned into the FseI and EcoRI sites of pMOS1_nanos-DENV-8miR_3xP3-EGFP constructs with In-Fusion HD Cloning technology (Clontech), resulting in pMOS1-AeCPA-DENV-8miR-3xP3-EGFP and pMOS1-AePUB-DENV-8miR-3xP3-EGFP constructs. Both plasmids were then used as donor plasmids for embryo microinjection.

**Generation of transgenic mosquitoes.** To establish the *AeCPA*-8miR and *AePub*-8miR transgenic mosquito lines, donor and helper plasmids were mixed with injection buffer (5 mM KCl and 0.1 mM NaH₂PO₄, pH 6.8) at concentrations of 1000 ng/µL and 600 ng/µL, respectively, and then microinjected into preblastoderm-stage embryos, as described in a related study<sup>52</sup>. Injected embryos were then hatched and reared under standard conditions. All the generation 0 (G0) eggs were collected. We used a fluorescent microscope (SZX12; Olympus) to screen G0 mosquito larvae to identify those with the EGFP fluorescent marker in their eyes and anal papillae, a pattern driven by the 3xP3 promoter. Positive G0 males were outcrossed to wild type females (Orlando strain) in a ratio of 1:3. G1 mosquitoes positive for the reporter were self-crossed by placing a single male with a single female. Eggs from G1 were collected, and the genotypes of the G1 parents were ascertained to distinguish the progeny genotypes. The progeny larvae of the following two generations were screened for eye-specific EGFP reporter gene expression to confirm their genotype. Further, we confirmed the integration of the transgenic cassettes into mosquito chromosomes by sequencing and Southern blot and used the heterozygous transgenic lines (expressing the eight miRNAs) as the experimental groups in the subsequent studies.

**Southern blot analysis.** The genomic DNA extractions were conducted in accordance with the manufacturer's protocol (Invitrogen). Approximately 2.5 µg of extracted total genomic DNA was digested by EcoRI/HindIII, followed by DNA separation on 0.8% agarose gel. The separated DNA was then transferred onto a nylon membrane and hybridized overnight at 42 °C with the oligonucleotide's probes, whose cover miR-1 to miR-8 coding sequences, respectively (Supplementary Table 4). We mixed these eight DNA oligonucleotide with equal molar and performed 5' <sup>32</sup>P-labeling by T4 polynucleotide kinase for hybridization. The expected size of hybridization pattern was 1970 bp digested in *AeCPA*-8miR and *AePUB*-8miR mosquitoes, respectively.
**DENV.** DENV serotypes 1–4 were propagated in C6/36 (from *Ae. albopictus*) cells. Briefly, C6/36 cells, grown in a 175 mm² flask at a confluency of approximately 80% at 28 °C, were infected at a multiplicity of infection of 0.5 for 3 h at 28 °C. Later, fresh medium was added and the infection proceeded for 5 days at 28 °C. Finally, the cell culture supernatant was filtered through a 0.22-µm filter, aliquoted, and stored at – 80 °C until use.

**Oral DENV infection of mosquitoes and 2nd blood feeding.** For DENV oral infection, virus supernatant from each serotype of DENV-infected C6/36 cells was mixed 1:1 with mouse blood to yield 10⁷ PFU/mL virus blood. This virus-containing blood was then provided to 150 female mosquitoes (Orlando, AeCPA-8miR, or AePUB-8miR lines at 7 to 10 days post-emergence) by artificial blood feeding for 40 min. Successfully blood-fed mosquitoes, which can be identified easily by visual cues, were separated from non-blood-fed mosquitoes and then incubated for 7 or 14 days at 28 °C and 80% relative humidity.

For the 2nd blood feeding experiment, Orlando or AeCPA-8miR mosquitoes at seven days post DENV infection were provided with another blood meal lacking the virus supernatant through artificial blood feeding for 40 min. Mosquitoes that were successfully blood-fed were transferred to a new cage and then incubated for the designated days.

**Plaque assay.** Viral plaque assays for both virus stocks and infected mosquitoes were performed using previously published protocols⁵³. BHK-21 cells were seeded into six-well plates at 2.5 × 10⁵ cells per well, with the addition of 2 mL of 5% fetal bovine serum (FBS) in RPMI medium. Plates were incubated at 37 °C for 16–24 h before the medium was removed, and 100 µL of a tenfold serial dilution of the virus was layered on top of the cells. Plates were shaken gently and incubated at 37 °C for 2 h. After absorption, we removed the infectious medium, and 3 mL of 1% methylcellulose (MC) medium (with 2% FBS, 100 mM sodium pyruvate, 1 M HEPES in RPMI medium, pH 7.4) was added to each well. Samples were incubated at 37 °C with 5% CO₂ supplemented for 5 days. Plates were harvested when plaques were visible to the naked eye. When the incubation was completed, the MC medium was removed and the plates were stained through the addition of 1% crystal violet solution to each well. Staining was performed for at least 3 h. Plates were subsequently rinsed and plaques were counted.

**DENV transmission assay.** Female mosquitoes of each experimental groups were infected orally with 10⁷ PFU/mL of DENV and then incubated for 7 or 14 days. Saliva of these infected mosquitoes were collected accordingly³¹,⁵⁴. Briefly, legs and wings of the infected mosquito was removed, followed by inserted its mouth-piece into a P200 tip which was filled with 5 µL of FBS. 30 min later, FBS in the tip was added into 45 µL of serum-free medium in RPMI and then analyzed with plaque assay. Tissues of this mosquitoes were subjected to the plaque assay as well. Transmission efficiency was calculated as follow:

$$ \text{Transmission efficiency} = \frac{\text{Number of mosquitoes with infectious saliva}}{\text{Number of examined mosquito}} $$

**qPCR.** Mosquito small RNAs were extracted using Tri-reagent (Merck) following the manufacturer's protocol. The reverse transcription (RT) of extracted RNAs (20 µL) proceeded as follows: the extracted RNA was mixed with the stem-loop RT primer (Supplementary Table 5), and the dNTPs were denatured at 65 °C for 5 min and then immediately cooled on ice. The treated RNA mixtures were incubated at 16 °C for 30 min followed by pulsed RT (30 °C for 30 s, 42 °C for 30 s, and 50 °C for 1 s for 60 cycles; terminated by 5 min at 85 °C) to increase detection sensitivity. The miRNA analysis used 2 µL of cDNA as a sample. Analysis was performed using a Universal Probe Library probe assay (Universal Probe Library probe #21) with qPCR (ABI PRISM 7900HT Sequence Detection System). Amplification proceeded with an initial step at 50 °C for 2 min; a denaturing step at 95 °C for 10 min; and 40 cycles of amplification at 95 °C for 15 s and 60 °C for 1 min. qPCR data was normalized against aae-U6 (as the reference gene) and used the endogenous miRNA aae-miR-100 as an endogenous control to obtain the expression profiles of anti-DENV miRNAs. Primers that were used in the miRNA analysis are listed in Supplementary Table 6.

**Statistical analysis.** Statistical Package for GraphPad Prism version 5.0 (http://www.graphpad.com) and Microsoft Excel were used for data entry, processing, and statistical analysis. Analysis of variance with a post-hoc Tukey's test was used to compare differences in means between different treatments. The Mann–Whitney rank sum test was used to analyze the difference in the titer of each virus between the wild-type (Orlando strain) and transgenic (AeCPA-8miR and AePUB-8miR) mosquitoes. The t test was used to analyze the difference in mosquito IR. A p value < 0.05 was considered significant.

**Data availability**

All data needed to evaluate the conclusions in the paper are present of the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Received: 22 March 2021; Accepted: 25 November 2021

Published online: 13 December 2021

**References**

1. Gubler, D. J. Dengue and dengue hemorrhagic fever. *Clin. Microbiol. Rev.* **11**, 480–496 (1998).
2. Brady, O. J. et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. *PLoS Negl. Trop. Dis.* **6**, e1760. https://doi.org/10.1371/journal.pntd.0001760 (2012).

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

4. Mairuhu, A. T., Wagenaar, J., Brandjes, D. P. & Gorp, E. C. Dengue: An arthropod-borne disease of global importance. *Eur. J. Clin. Microbiol. Infect. Dis.* **23**, 425–433 (2004).

5. Rigau-Pérez, J. et al. Dengue and dengue haemorrhagic fever. *Lancet* **352**, 971–977 (1998).

6. Cheng, L. et al. Releasing intracellular NS1 from mosquito cells for the detection of dengue virus-infected mosquitoes. *Viruses* **12**, 1105 (2020).

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

8. Salazar, M. I. et al. Dengue virus type 2: Replication and tropisms in orally infected *Aedes aegypti* mosquitoes. *BMC Microbiol.* **7**, 9 (2007).

9. Iturbe-Ormaetxe, I., Walker, T. & O'Neill, S. L. Wolbachia and the biological control of mosquito-borne disease. *EMBO Rep.* **12**, 508–518 (2011).

10. Franz, A. W. et al. Stability and loss of a virus resistance phenotype over time in transgenic mosquitoes harbouring an antiviral effector gene. *Insect Mol. Biol.* **18**, 661–672 (2009).

11. Mathur, G. et al. Transgene-mediated suppression of dengue viruses in the salivary glands of the yellow fever mosquito, *Aedes aegypti*. *Insect Mol. Biol.* **19**, 753–763 (2010).

12. Champer, J., Buchman, A. & Akbari, O. S. Cheating evolution: Engineering gene drives to manipulate the fate of wild populations. *Nat. Rev. Genet.* **17**, 146–159 (2016).

13. Travanty, E. A. et al. Using RNA interference to develop dengue virus resistance in genetically modified *Aedes aegypti*. *Insect Biochem. Mol. Biol.* **34**, 607–613 (2004).

14. Carvalho, D. O. et al. Mass production of genetically modified *Aedes aegypti* for field releases in Brazil. *J. Vis. Exp.* **83**, e3579. https://doi.org/10.3791/3579 (2014).

15. Wilke, A. B. B., Beier, J. C. & Benelli, G. Transgenic mosquitoes—Fact or fiction?. *Trends Parasitol.* **34**, 456–465 (2018).

16. Catteruccia, F., Crisanti, A. & Wimmer, E. A. Transgenic technologies to induce sterility. *Malar. J.* **8**, S7 (2009).

17. Saldana, M. A. et al. Zika virus alters the microRNA expression profile and elicits an RNAi response in *Aedes aegypti* mosquitoes. *PLoS Negl. Trop. Dis.* **11**, e0005760 (2017).

18. Campbell, C. L., Harrison, T., Hess, A. M. & Ebel, G. D. MicroRNA levels are modulated in *Aedes aegypti* after exposure to Dengue-2. *Insect Mol. Biol.* **23**, 132–139 (2014).

19. Maharaj, P. D., Widen, S. G., Huang, J., Wood, T. G. & Thangamani, S. Discovery of mosquito saliva microRNAs during CHIKV infection. *PLoS Negl. Trop. Dis.* **9**, e0003386 (2015).

20. Fire, A. et al. Potent and specific genetic interference by double-stranded RNA in *Caenorhabditis elegans*. *Nature* **391**, 806–811 (1998).

21. Ding, S. W. & Li, W. X. Viral suppressors of RNA silencing. *Curr. Opin. Biotechnol.* **12**, 150–154 (2001).

22. Ghildiyal, M., Xu, J., Seitz, H., Weng, Z. & Zamore, P. D. Sorting of Drosophila small silencing RNAs partitions microRNA strands into the RNA interference pathway. *RNA* **16**, 43–56 (2010).

23. Forstemann, K., Horwich, M. D., Wee, L., Tomari, Y. & Zamore, P. D. Drosophila microRNAs are sorted into functionally distinct argonaute complexes after production by dicer-1. *Cell* **130**, 287–297 (2007).

24. Yang, J. S. et al. Intertwined pathways for Argonaute-mediated microRNA biogenesis in *Drosophila*. *Nucleic Acids Res.* **42**, 1987–2002 (2014).

25. Li, Z. & Rana, T. M. Therapeutic targeting of microRNAs: Current status and future challenges. *Nat. Rev. Drug Discov.* **13**, 622–638 (2014).

26. Wilczynska, A. & Bushell, M. The complexity of miRNA-mediated repression. *Cell Death. Differ.* **22**, 22–33 (2015).

27. Anderson, M. A., Gross, T. L., Myles, K. M. & Adelman, Z. N. Validation of novel promoter sequences derived from two endogenous ubiquitin genes in transgenic *Aedes aegypti*. *Insect Mol. Biol.* **19**, 441–449 (2010).

28. Moreira, L. A. et al. Robust gut-specific gene expression in transgenic *Aedes aegypti* mosquitoes. *Proc. Natl. Acad. Sci.* **97**, 10895–10898 (2000).

29. Scott, T. W. et al. Longitudinal studies of *Aedes aegypti* (Diptera: Culicidae) in Thailand and Puerto Rico: Blood feeding frequency. *J. Med. Entomol.* **37**, 89–101 (2000).

30. Lobo, N. F., Clayton, J. R., Fraser, M. J., Kafatos, F. C. & Collins, F. H. High efficiency germ-line transformation of mosquitoes. *Nat. Protoc.* **1**, 1312–1317 (2006).

31. Yen, P. S., James, A., Li, J. C., Chen, C. H. & Failoux, A. B. Synthetic miRNAs induce dual arboviral-resistance phenotypes in the vector mosquito *Aedes aegypti*. *Commun. Biol.* **1**, 11–19 (2018).

32. Buchman, A. et al. Engineered resistance to Zika virus in transgenic *Aedes aegypti* expressing a polycistronic cluster of synthetic small RNAs. *Proc. Natl. Acad. Sci.* **116**, 3656–3661 (2016).

33. Hussain, M. & Asgari, S. MicroRNA-like viral small RNA from Dengue virus 2 autoregulates its replication in mosquito cells. *Proc. Natl. Acad. Sci.* **111**, 2746–2751 (2014).

34. Franz, A. W. et al. Engineering RNA interference-based resistance to dengue virus type 2 in genetically modified *Aedes aegypti*. *Proc. Natl. Acad. Sci.* **103**, 4198–4203 (2006).

35. Olson, K. et al. Genetically engineered resistance to dengue-2 virus transmission in mosquitoes. *Science* **272**, 884–886 (1996).

36. Gaines, P. et al. Pathogen-derived resistance to dengue type 2 virus in mosquito cells by expression of the premembrane coding region of the viral genome. *J. Virol.* **70**, 2132–2137 (1996).

37. Ruktanonchai, N. W., Lounibos, L. P., Smith, D. L. & Allan, S. A. Experience-and age-mediated oviposition behavior in the yellow fever mosquito *Stegomyia aegypti* (*Aedes aegypti*). *Med. Vet. Entomol.* **29**, 255–262 (2015).

38. Zeller, M. & Koella, J. C. Effects of food variability on growth and reproduction of *Aedes aegypti*. *Ecol. Evol.* **6**, 552–559 (2016).

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

40. Cheng, L. et al. Transgenic expression of human C-type lectin protein CLEC18A reduces dengue virus type 2 infectivity in *Aedes aegypti*. *Front. Immunol.* **12**, 640367 (2021).

41. Fabian, M. & Sonenberg, N. The mechanics of miRNA-mediated gene silencing: A look under the hood of miRISC. *Nat. Struct. Mol. Biol.* **19**, 586 (2012).

42. Park, S. et al. MicroRNA clustering on the biogenesis of suboptimal microRNAs. *Appl. Biol. Chem.* **64**, 51 (2021).

43. Zhang, X. et al. The expression of microRNA and microRNA clusters in the aging heart. *PLoS ONE* **7**, e34688 (2012).

44. Seyhan, A. A. A multiplexed miRNA and transgene expression platform for simultaneous repression and expression of protein coding sequences. *Mol. Biosyst.* **12**, 295 (2016).

45. O'Brien, J., Hayder, H., Zayed, Y. & Peng, C. Overview of microRNA biogenesis, mechanisms of actions, and circulation. *Front. Endocrinol.* **9**, 402 (2018).

46. Kumiko, U. T. Is the efficiency of RNA silencing evolutionarily regulated? *Int. J. Mol. Sci.* **17**, e719 (2016).

47. Patil, B. L. & Fauquet, C. M. Light intensity and temperature affect systemic spread of silencing signal in transient agroinfiltration studies. *Mol. Plant Pathol.* **16**, 484–494 (2015).
48. Xiao, F. Z. et al. The effect of temperature on the extrinsic incubation period and infection rate of dengue virus serotype 2 infection in *Aedes albopictus*. *Arch. Virol.* **159**, 3053–3057 (2014).

49. Rohani, A., Wong, Y. C., Zamre, I., Lee, H. L. & Zurainee, M. N. The effect of extrinsic incubation temperature on development of dengue serotype 2 and 4 viruses in *Aedes aegypti*. *Southeast Asian J. Trop. Med. Public Health* **40**, 942–950 (2009).

50. Richards, S. L., Mores, C. N., Lord, C. C. & Tabachnick, W. J. Impact of extrinsic incubation temperature and virus exposure on vector competence of *Culex pipiens* quinquefasciatus Say (Diptera: Culicidae) for West Nile virus. *Vector Borne Zoonotic Dis.* **7**, 629–636 (2007).

51. Chen, C. H. et al. A synthetic maternal-effect selfish genetic element drives population replacement in *Drosophila*. *Science* **316**, 597–600 (2007).

52. Li, H. H. et al. C-type lectins link immunological and reproductive processes in *Aedes aegypti*. *iScience*. **23**, 101486 (2020).

53. Roehrig, J. T., Hombac, J. & Barrett, A. D. Guidelines for plaque-reduction neutralization testing of human antibodies to dengue viruses. *Viral Immunol.* **21**, 123–132 (2008).

54. Dubrulle, M., Mousson, L., Moutailler, S., Vazeille, M. & Failloux, A. B. Chikungunya virus and Aedes mosquitoes: Saliva is infectious as soon as two days after oral infection. *PLoS ONE* **4**, e5895 (2009).

**Author contributions**

W.L.L., C.W.H., S.P.C. performed data curation, formal analysis and assisted in writing the original and revised manuscript. P.S.Y., M.P.S., J.C.L., H.H.L. L. C., C.K.T., S.H.K., S.P.C. supervised data generation, and assisted the experiments. H.K.T. and Z.T.YT. provided subject matter expertise and supervised data generation. C.H.C., A.B.F., O.S.A. conceptualized the project, C.H.C. performed data analysis, secured funding, and wrote the paper. All authors reviewed the manuscript.

**Funding**

This study also was supported by the National Health Research Institutes (NHRI), Taiwan (Grant No. NHRI-MR-110-GP-12) and Ministry of Health and Welfare (Grant No. MOHW103-TDU-B-212-000005) awarded to Chun-Hong Chen. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**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-021-03229-4.

**Correspondence** and requests for materials should be addressed to C.-H.C.

**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_39.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) 2021, corrected publication 2022