# Comparison of chikungunya viruses generated using infectious clone or the Infectious Subgenomic Amplicons (ISA) method in Aedes mosquitoes

Souand Mohamed Ali, Anubis Vega-Rúa, Jean-Sélim Driouich, Xavier de Lamballerie, Anna-Bella Failloux, Antoine Nougairède

► To cite this version:

Souand Mohamed Ali, Anubis Vega-Rúa, Jean-Sélim Driouich, Xavier de Lamballerie, Anna-Bella Failloux, et al.. Comparison of chikungunya viruses generated using infectious clone or the Infectious Subgenomic Amplicons (ISA) method in Aedes mosquitoes. PLoS ONE, 2018, 13 (6), pp.e0199494. 10.1371/journal.pone.0199494. pasteur-01843209

HAL Id: pasteur-01843209

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

Submitted on 18 Jul 2018

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
RESEARCH ARTICLE

# Comparison of chikungunya viruses generated using infectious clone or the Infectious Subgenomic Amplicons (ISA) method in *Aedes* mosquitoes

Souand Mohamed Ali¹*, Anubis Vega-Rúa²*, Jean-Sélim Driouich¹, Xavier de Lamballerie¹, Anna-Bella Failloux², Antoine Nougairède¹*

1 Unité des Virus Émergents (UVE: Aix-Marseille Univ-IRD 190-Inserm 1207-IHU Méditerranée Infection),

Marseille, France, 2 Department of Virology, Arboviruses and Insect Vectors, Institut Pasteur, Paris, France

* These authors contributed equally to this work.

² Current address: Laboratoire d'Entomologie Médicale, Institut Pasteur de la Guadeloupe, Abymes, Guadeloupe

* antoine.nougairede@univ-amu.fr

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

## Abstract

Reverse genetics systems provide the opportunity to manipulate viral genomes and have been widely used to study RNA viruses and to develop new antiviral compounds and vaccine strategies. The recently described method called ISA (Infectious Subgenomic Amplicons) gives the possibility to rescue RNA viruses in days. We demonstrated in cell culture that the use of the ISA method led to a higher genetic diversity of viral populations than that observed using infectious clone technology. However, no replicative fitness difference was observed. In the present study, we used the chikungunya virus as a model to compare in *Aedes aegypti* and *Aedes albopictus* mosquitoes the genotypic and phenotypic characteristics of viruses produced either from an infectious clone or using the ISA method. We confirmed the results found in *cellulo* corroborating that the use of the ISA method was associated with higher genetic diversity of viral populations in mosquitoes but did not affect the vector competence validating its use for *in vivo* experiments.

## Introduction

Reverse genetics systems that give the possibility to rescue infectious viruses from DNA copies of their genomes, are important tools to explore viral life cycle and to contribute to the development of new antiviral compounds and vaccine candidates [1-4].

Numerous reverse genetics systems that allow producing wild type and genetically modified viruses have been previously developed, each associated to specific benefits and limitations [1]. Infectious clones are widely used but commonly difficult to create, in particular because of the instability and toxicity of some viral sequences into bacteria [5-7]. To overcome these difficulties, bacterium-free approaches alternative methods have been developed, such as the ISA (infectious subgenomic amplicons) method that was recently applied to a large panel of single-stranded positive-sense RNA viruses [6, 8-10].

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

### OPEN ACCESS

**Citation:** Mohamed Ali S, Vega-Rúa A, Driouich J-S, de Lamballerie X, Failloux A-B, Nougairède A (2018) Comparison of chikungunya viruses generated using infectious clone or the Infectious Subgenomic Amplicons (ISA) method in *Aedes* mosquitoes. PLoS ONE 13(6): e0199494. https://doi.org/10.1371/journal.pone.0199494

**Editor:** Lisa F. P. Ng, Singapore Immunology Network, SINGAPORE

**Received:** January 31, 2018

**Accepted:** June 10, 2018

**Published:** June 28, 2018

**Copyright:** © 2018 Mohamed Ali et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the French "Agence Nationale de la Recherche" under grant agreement no. ANR-14-CE14-0001 (RNA Cacci-Code) and by the European Union's Horizon 2020 Research and Innovation Programme under grants agreements no. 653316 (European Virus Archive goes global project: https://www.european-virus-
archive.com/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

Using the Chikungunya virus (CHIKV; family *Togaviridae*; genus *Alphavirus*) as a model, we previously demonstrated *in cellulo* that the use of the ISA method generated a higher genetic diversity of viral populations than that observed with the use of an infectious clone [11]. We concluded that PCR-based reverse genetics methods -such as the ISA method- are associated with artificial viral heterogeneity because the polymerases employed during PCR amplification generate a number of undesired nucleotide changes. Since the structure of viral populations can shape virus phenotype, we previously investigated the impact of using the ISA method on viral phenotype *in cellulo*, in comparison with using an infectious clone [12–14]. We did not identify any difference in the replicative fitness of the viruses produced by both methods [11].

Accordingly, the main objective of the current study was to compare *in vivo* the impact on genotype and phenotype of using either the ISA method or an infectious clone. As in our previous work, we used CHIKV as a model virus and investigated the genotypic and phenotypic characteristics of viruses produced by both methods in the two main mosquito vectors of CHIKV, *Aedes aegypti* and *Aedes albopictus*.

# Results and discussion

## Production of viruses

A previously described infectious clone of the LR2006 strain of CHIKV was used (i) to transfect cells and generate infectious particles (named IC_virus) and (ii) to produce subgenomic amplicons used during the ISA procedure to rescue infectious virus (named ISA_virus) [15]. Infectious clone and subgenomic amplicons were transfected into HEK-293 cells. Viruses were passaged twice in Vero cells. Virus stocks were used to infect mosquitoes.

## Mosquito oral infections

*Ae. aegypti* and *Ae. albopictus* females were orally infected with an infectious blood meal (final infectious titre: 10⁶.⁸ TCID₅₀/mL). After removing non-engorged females, mosquitoes were incubated 6/9 days. For each experimental condition, a group of 30 mosquitoes was used to assess vector competence.

## Evaluation of genetic diversity in virus stocks and *Aedes* mosquitoes

The genetic diversity was evaluated using complete viral genome sequences from virus stocks used to infect mosquitoes and from positive mosquito heads (presence of virus was assessed using a TCID₅₀ assay; see below). Data were generated by next generation sequencing and we considered in our analyses all the substitutions with a frequency ≥2%.

We first evaluated the genetic diversity in virus stocks used to infect mosquitoes (the sequence of the infectious clone was used as a reference). Not surprisingly, transfection of the infectious clone led to production of clonal viral population and no mutation was detected in the IC_virus. In contrast, transfection of amplicons led to a higher genetic diversity and 13 mutations were detected in ISA_virus (Table 1). This result is in accordance with our previous work [11]. Of note, only minority mutations (substitution frequency <50%) were detected in ISA_virus and therefore both viruses had identical consensus sequence.

We then evaluated the genetic diversity from grounded heads of *Aedes* mosquito. Based on our experiment, it is difficult to generate complete genome sequences from individual mosquito heads. Thereby, we initiated analyses with eight heads per condition. Finally, data from three to eight heads per condition were used. To evaluate the impact of the reverse genetics method used, we compared the number of minority and majority mutations (substitution
<table><caption>Table 1. Substitutions generated using the ISA method and detected in *Aedes* mosquitoes infected by the ISA_virus.</caption><thead><tr><th rowspan="2">Nucleotide position</th><th rowspan="2">Original substitution frequency*</th><th rowspan="2">Region</th><th rowspan="2">nucleotide change</th><th rowspan="2">aa change</th><th colspan="3">Number of mosquitoes with the substitution detected on viral genome</th></tr><tr><th>A. aegypti 6 dpi (n = 8)</th><th>A. aegypti 9 dpi (n = 6)</th><th>A. albopictus 6 dpi (n = 5)</th></tr></thead><tbody><tr><td>502</td><td>0,21</td><td>nsP1</td><td>A→G</td><td>-</td><td>1</td><td>0</td><td>1</td></tr><tr><td>710</td><td>0,03</td><td>nsP1</td><td>T→C</td><td>-</td><td>0</td><td>0</td><td>0</td></tr><tr><td>881</td><td>0,03</td><td>nsP1</td><td>A→C</td><td>K→Q</td><td>0</td><td>0</td><td>0</td></tr><tr><td>1105</td><td>0,03</td><td>nsP1</td><td>C→A</td><td>-</td><td>0</td><td>0</td><td>0</td></tr><tr><td>1590</td><td>0,22</td><td>nsP1</td><td>T→G</td><td>L→R</td><td>1</td><td>0</td><td>1</td></tr><tr><td>2458</td><td>0,28</td><td>nsP2</td><td>T→C</td><td>-</td><td>4</td><td>0</td><td>2</td></tr><tr><td>6923</td><td>0,23</td><td>nsP4</td><td>A→G</td><td>T→A</td><td>3</td><td>2</td><td>4</td></tr><tr><td>7628</td><td>0,28</td><td>C</td><td>C→T</td><td>P→L</td><td>4</td><td>0</td><td>2</td></tr><tr><td>8612</td><td>0,02</td><td>E2</td><td>A→G</td><td>E→G</td><td>0</td><td>0</td><td>0</td></tr><tr><td>8766</td><td>0,15</td><td>E2</td><td>A→G</td><td>-</td><td>4</td><td>0</td><td>2</td></tr><tr><td>9100</td><td>0,1</td><td>E2</td><td>A→G</td><td>N→D</td><td>0</td><td>3</td><td>1</td></tr><tr><td>9690</td><td>0,02</td><td>E2</td><td>T→C</td><td>-</td><td>0</td><td>0</td><td>0</td></tr><tr><td>11644</td><td>0,1</td><td>3'UTR</td><td>A→G</td><td>-</td><td>0</td><td>3</td><td>1</td></tr></tbody></table>

List of the substitutions detected in virus stock (ISA_virus) used to infect *Aedes* mosquitoes. None of these mutations were detected in mosquitoes infected by the IC_virus. “dpi” means days post-infection.

https://doi.org/10.1371/journal.pone.0199494.t001

frequency < and ≥ 50% respectively) found for each condition. No clear trend can be found with minority mutations (Fig 1A): no impact of the reverse genetic method was observed in *Ae. aegypti* and *Ae. albopictus* heads at day 6 post-infection and the ISA method was associated with a slightly (not significant) lower number of mutations at day 9 post-infection in *Ae. aegypti*. Rather, the ISA method was associated with a higher number of majority mutations in all conditions (Fig 1B). This difference was significant in *Ae. aegypti* mosquitoes (Wilcoxon test; *p* value = 0.012 and 0.047 at 6 and 9 days post-infection, respectively). These findings suggested that in mosquito infected by ISA_virus, some of minority variants present in the blood meals were probably selected during the passage of the midgut barrier. To confirm this hypothesis, we specifically followed the mutations detected on ISA_virus stock (Table 1). We observed that none of the mutation initially present with a frequency <10% were found in mosquitoes. In contrast, all the mutations initially present with a frequency ≥10% were found in at least 2 mosquitoes (ranging between 2 and 9). Most of them were majority mutations in mosquitoes (21/39 of which 21 were totally fixed). Altogether, these findings demonstrated that genetic diversity in mosquitoes was higher when the ISA method was used to produce the virus.

In all conditions, the proportion of transition ranged between 56% and 100% (Fig 1C). The ISA method was associated with a significant higher number of transition at day 9 post-infection in *Ae. aegypti* (Fischer test; *p* value = 0.005). The reasons for this discrepancy remained unclear. No significant difference in the proportion of synonymous and non-synonymous mutations was found (Fig 1D). In all conditions, mutations detected in mosquitoes were distributed throughout the whole genome (Fig A in S1 File).

## Evaluation of vector competence in *Aedes* mosquitoes

To assess the ability of *Ae. aegypti* and *Ae. albopictus* mosquitoes to allow IC_virus and ISA_virus to pass through the midgut barrier, dissemination efficiency (DE; proportion of

<footnote>* frequency found in virus stock.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Genetic diversity in *Aedes* mosquitoes. Panels A and B represent the number of minority and majority mutations respectively detected in mosquito heads (substitution frequency < and ≥ 50% respectively). Mutation characteristics are represented in panels C (Non synonymous; Synonymous mutations) and D (Transition; Transversion). In Panels A and B, the bottom and top of the box represent first and third quartiles, the band inside the box represents median value and the ends of the error bars represent minimum and maximum values. IC and ISA mean group of mosquitoes infected by IC_virus and ISA_virus respectively. “dpi” means days post-infection. * means that a significant difference (p value <0.05) was observed between mosquitoes infected by IC_virus and ISA_virus.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0199494.g001

mosquitoes with infectious virus in the head) was determined by testing individual grounded heads with a TCID₅₀ assay (Fig 2A). Similar DE values ranging from 83 to 93% were observed with *Ae. aegypti* at days 6 and 9 post-infection for both viruses. Similarly, close DE values ranging from 53 to 57% were found with *Ae. albopictus* at day 6 post-infection. The intensity of viral dissemination was also assessed by comparing amounts of infectious particles in collected heads (only positive heads were taken into account). In all conditions, mean values were similar ranging from 4.0 to 4.2 log₁₀ TCID₅₀/head (Fig 2B). All these findings indicated that the virus dissemination in *Aedes* mosquitoes was not significantly impacted by the reverse genetic method used to generate the virus.

To determine the ability of *Ae. aegypti* and *Ae. albopictus* mosquitoes to sustain transmission of IC_virus and ISA_virus, transmission efficiency (TE; proportion of mosquito with infectious virus in the saliva) was assessed by testing individual collected saliva with a TCID₅₀ assay (Fig 2C). TE values ranging from 37 to 47% were found in all cases except with *Ae.*
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Vector competence in *Aedes* mosquitoes. Dissemination efficiency (DE; proportion of mosquito with infectious virus in the head) was determined by testing mosquito heads with a TCID₅₀ assay (Panel A). The intensity of viral dissemination was assessed by comparing amounts of infectious particles in collected heads (Panel B; only positive heads were taken into account). Transmission efficiency (TE; proportion of mosquitoes with infectious virus in the saliva) was determined by testing mosquito saliva with a TCID₅₀ assay (Panel C). The intensity of viral transmission was assessed by comparing amounts of infectious particles in collected saliva (Panel D; only positive saliva were taken into account). In panels B and D, the average values are shown and error bars represent standard deviation. IC and ISA mean group of mosquitoes infected by IC_virus and ISA_virus respectively. "dpi" means days post-infection.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0199494.g002

*aegypti* infected with the ISA_virus at day 6 post-infection (60%). However, this difference was not significant. The intensity of viral transmission was also assessed by comparing amounts of infectious particles in collected saliva (only positive heads were taken into account). In all conditions, mean values were similar ranging from 0.8 to 1.3 log₁₀ TCID50/saliva (Fig 2D). All these findings indicated that the virus transmission was not affected by the reverse genetic method used to generate the virus.

## Conclusion

The ISA method represents a technological breakthrough for the production of recombinant RNA viruses [10]. We previously demonstrated that the ISA method conducted to higher genetic diversity of viral populations in cellulo but this genotypic difference was associated with no significant difference in replicative fitness [11]. In this study, we compared in *Aedes* mosquitoes phenotypic and genotypic characteristics of CHIKVs produced using ISA- and IC-based reverse genetics methods. We observed a higher genetic diversity of viral population in heads of mosquitoes infected by viruses generated using the ISA method. However, we found
that vector competence was not affected by this genotypic difference since it was identical regardless the reverse genetic procedure used validating the use of both methods for *in vivo* experiments. These findings suggest that a certain level of artificial heterogeneity in virus population is easily supported by CHIKV during its replication in *Aedes* mosquitoes highlighting the high level of mutational robustness of this RNA virus. However, it is also possible that deleterious mutants were already negatively selected during the first round of viral replication in cell culture. This hypothesis remains to be tested by analysing the early evolutionary events that occurred *in cellulo*.

# Methods

## Cells

Vero cells (ATCC number CCL-81) were maintained at 37°C/5% CO₂ in minimal essential medium (Life Technologies) with 7% heat-inactivated foetal bovine serum (FBS; Life Technologies), 1% penicillin/streptomycin (PS; 5000 U.ml⁻¹ and 5000 µg.ml⁻¹; Life Technologies) and 1% Glutamine (Gln; 200 mmol.l⁻¹; Life Technologies). HEK-293 cells (ATCC number CCL-1573) were maintained at 37°C/5% CO₂ in the same medium supplemented with 1% non-essential amino acids (Life Technologies).

## Mosquitoes

We used in this study (i) the PAEA colony of *Ae. Aegypti* established from a population collected as larvae in urban settings of French Polynesia in 1994 and (ii) the BP strain of *Ae. Albopictus* collected as larvae in urban settings of Binh Phuoc (Vietnam) in 2011 and maintained in laboratory for 9 generation before experimental infections. After egg hatching, larvae were split in pans filled with one liter tap water supplemented with yeast tablets. Emerging adults were maintained at 28°C, 80% relative humidity in cages with a 14h:10h light:dark cycle and supplied ad libitum with a 10% sucrose solution.

## CHIK infectious clone

We used a previously described CHIKV infectious clone in this study (LR2006 strain of CHIKV [15]; GenBank accession EU224268).

## Preparation of DNA fragments for the ISA procedure

DNA fragments for the ISA procedure were prepared as previously described [11]. Briefly, the complete viral genome was amplified by PCR, using the infectious clone as template, in three overlapping DNA fragments. The 5'-extremity of the first fragment and the 3'-extremity of the last fragment were flanked respectively by the human cytomegalovirus promoter (pCMV) and the hepatitis delta virus ribozyme followed by the simian virus 40 polyadenylation signal (HDR/SV40pA). As previously described, PCR products were purified and treated by a digestion step with the restriction enzyme *Dpn*I (New England Biolabs) to ensure complete removal of the infectious clone used as template [10].

## Cell transfection

Cell transfection was performed as previously described [11]. Briefly, a total amount of 1 µg of DNA (infectious clone or equimolar mix of the three DNA fragments) were transfected into a 12.5 cm² culture flask of subconfluent HEK-293 cells. Each virus was then passaged twice using Vero cells. Clarified cell supernatants (virus stocks) were used to perform vector competence study and whole-genome sequencing.
## Vector competence study

Five to seven day-old females were fed during 30 min with an infectious blood meal containing 2 ml of washed rabbit erythrocytes and 1 ml of viral suspension supplemented with a phagostimulant (ATP) at a final concentration of 5 mM. The final infectious titer in blood meals was $10^{6.8}$ TCID50/mL. Non-engorged females were then discarded. Engorged females were transferred in cardboard containers and maintained at 28°C with 10% sucrose.

Batches of 30 mosquitoes were analyzed at days 6 and 9 post-infection for each condition described in Fig 2. Days post-infection were chosen based on kinetics of CHIKV dissemination and transmission efficiencies in Ae. aegypti and Ae. albopictus mosquitoes (Fig B in S1 File). To assess viral dissemination, mosquito heads were removed and ground in 250 µL of Leibovitz L15 medium (Life Technologies) supplemented with 2% FBS. Clarified supernatant was used to perform a TCID50 assay. To assess viral transmission, saliva was collected from individual females. For saliva collection, wings and legs were removed and the proboscis was inserted into a 20 µL tip containing 5 µL of FBS. After 45 min of salivation, FBS containing saliva was transferred into 45 µL of Leibovitz for titration (TCID50 assay).

## Tissue-culture infectious dose 50 (TCID₅₀) assay

96-well plates culture of confluent Vero cells were used. Serial dilutions of samples were tested as previously described [15]. The procedure varied depending of the nature of the sample tested.

**1/ Mosquito samples.** Mosquito heads: each serial 10-fold dilution was tested four times (first dilution: 1/5). Mosquito saliva: because sample size is very small, each serial 5-fold dilution was tested once (first dilution: 1/2.5). Amphotericine B was added to cell medium (final concentration: 2.50 µg/mL; Life Technologies). Serial dilutions of samples were removed 2 hours after infection and 150 µL of fresh medium (3% FBS) was added in each well.

**2/ Cell supernatants.** Each serial 10-fold dilution was tested six times (first dilution: 1/10) as previously described [15].

In all cases, plates were incubated for 7 days and scored for absence or presence of cytopathic effect (CPE) in each well. Determination of infectious titres TCID₅₀/ml was performed using the method of Reed and Muench [16].

## Complete genome sequencing

Viral RNA was extracted from 80 µL of clarified supernatants from grounded heads using the NucleoSpin 8 / 96 Virus kit (Macherey-Nagel). As previously described, set of specific primer pairs was used to generate overlapping amplicons covering the full length genome [10].

Sequencing was performed using the Ion PGM Sequencer [17] (Life Technologies). Read analyses were performed as previously described [18]. To assess the genetic diversity of viral populations, mutation frequency for each position was calculated as the number of mutated reads divided by the total number of reads at that site. Only substitution with a frequency ≥2% were taken into account for analysis (Table A in S1 File).

## Statistical analysis

All the tests were carried out with the R software [19]. This includes Shapiro-Wilk test, Wilcoxon test, Fischer test, and Student t-test.

## Supporting information

**S1 File.** Figure A: Mutation distribution.

Figure B: kinetics of IC_virus dissemination and transmission efficiencies in *Aedes*
mosquitoes.

Table A: List of the mutations detected in this study.

(PDF)

## Author Contributions

**Conceptualization:** Anubis Vega-Rúa, Xavier de Lamballerie, Anna-Bella Failloux, Antoine Nougairède.

**Data curation:** Souand Mohamed Ali.

**Formal analysis:** Souand Mohamed Ali, Anubis Vega-Rúa, Jean-Sélim Driouich, Antoine Nougairède.

**Funding acquisition:** Xavier de Lamballerie, Anna-Bella Failloux.

**Investigation:** Souand Mohamed Ali, Anubis Vega-Rúa, Antoine Nougairède.

**Methodology:** Souand Mohamed Ali, Anubis Vega-Rúa, Jean-Sélim Driouich.

**Supervision:** Xavier de Lamballerie, Anna-Bella Failloux, Antoine Nougairède.

**Writing – original draft:** Souand Mohamed Ali, Antoine Nougairède.

**Writing – review & editing:** Anubis Vega-Rúa, Jean-Sélim Driouich, Xavier de Lamballerie, Anna-Bella Failloux.

## References

1. Aubry F, Nougairede A, Gould EA, de Lamballerie X. Flavivirus reverse genetic systems, construction techniques and applications: a historical perspective. Antiviral research. 2015; 114:67–85. Epub 2014/12/17. https://doi.org/10.1016/j.antiviral.2014.12.007 PMID: 25512228.

2. Neumann G, Kawaoka Y. Reverse genetics systems for the generation of segmented negative-sense RNA viruses entirely from cloned cDNA. Curr Top Microbiol Immunol. 2004; 283:43–60. Epub 2004/08/10. PMID: 15298167.

3. Yang CC, Hu HS, Wu RH, Wu SH, Lee SJ, Jiaang WT, et al. A novel dengue virus inhibitor, BP13944, discovered by high-throughput screening with dengue virus replicon cells selects for resistance in the viral NS2B/NS3 protease. Antimicrob Agents Chemother. 2014; 58(1):110–9. https://doi.org/10.1128/AAC.01281-13 PMID: 24145533; PubMed Central PMCID: PMCPMC3910792.

4. Chambers TJ, Nestorowicz A, Mason PW, Rice CM. Yellow fever/Japanese encephalitis chimeric viruses: construction and biological properties. J Virol. 1999; 73(4):3095–101. PubMed Central PMCID: PMCPMC104070. PMID: 10074160

5. Rice CM, Grakoui A, Galler R, Chambers TJ. Transcription of infectious yellow fever RNA from full-length cDNA templates produced by in vitro ligation. New Biol. 1989; 1(3):285–96. Epub 1989/12/01. PMID: 2487295.

6. Gritsun TS, Gould EA. Infectious transcripts of tick-borne encephalitis virus, generated in days by RT-PCR. Virology. 1995; 214(2):611–8. Epub 1995/12/20. doi: S0042-6822(85)70072-4 [pii] https://doi.org/10.1006/viro.1995.0072 PMID: 8553563.

7. Ruggli N, Rice CM. Functional cDNA clones of the Flaviviridae: strategies and applications. Adv Virus Res. 1999; 53:183–207. Epub 1999/12/03. PMID: 10582099.

8. Edmonds J, van Grinsven E, Prow N, Bosco-Lauth A, Brault AC, Bowen RA, et al. A novel bacterium-free method for generation of flavivirus infectious DNA by circular polymerase extension reaction allows accurate recapitulation of viral heterogeneity. J Virol. 2013; 87(4):2367–72. Epub 2012/12/14. doi: JVI.03162-12 [pii] https://doi.org/10.1128/JVI.03162-12 PMID: 23236063; PubMed Central PMCID: PMC3571472.

9. Siridechadilok B, Gomutsukhavadee M, Sawaengpol T, Sangiambut S, Puttikhunt C, Chin-inmanu K, et al. A simplified positive-sense-RNA virus construction approach that enhances analysis throughput. J Virol. 2013; 87(23):12667–74. Epub 2013/09/21. doi: JVI.02261-13 [pii] https://doi.org/10.1128/JVI.02261-13 PMID: 24049164; PubMed Central PMCID: PMC3838137.

10. Aubry F, Nougairede A, de Fabritus L, Querat G, Gould EA, de Lamballerie X. Single-stranded positive-sense RNA viruses generated in days using infectious subgenomic amplicons. The Journal of general
virology. 2014; 95(Pt 11):2462-7. Epub 2014/07/24. https://doi.org/10.1099/vir.0.068023-0 PMID: 25053561; PubMed Central PMCID: PMCPmc4202267.

11. Driouich JS, Ali SM, Amroun A, Aubry F, de Lamballerie X, Nougairede A. SuPReMe: a rapid reverse genetics method to generate clonal populations of recombinant RNA viruses. Emerging microbes & infections. 2018; 7(1):40. https://doi.org/10.1038/s41426-018-0040-2 PMID: 29559627; PubMed Central PMCID: PMCPMC5861084.

12. Holmes EC, Moya A. Is the quasispecies concept relevant to RNA viruses? J Virol. 2002; 76(1):460-5. Epub 2001/12/12. https://doi.org/10.1128/JVI.76.1.460-462.2002 PMID: 11739715; PubMed Central PMCID: PMC135735.

13. Ciota AT, Ngo KA, Lovelace AO, Payne AF, Zhou Y, Shi PY, et al. Role of the mutant spectrum in adaptation and replication of West Nile virus. J Gen Virol. 2007; 88(Pt 3):865-74. Epub 2007/02/28. doi: 88/3/865 [pii] https://doi.org/10.1099/vir.0.82606-0 PMID: 17325359; PubMed Central PMCID: PMC3249657.

14. Domingo E, Gonzalez-Lopez C, Pariente N, Airaksinen A, Escarmis C. Population dynamics of RNA viruses: the essential contribution of mutant spectra. Arch Virol Suppl. 2005;(19):59-71. PMID: 16355868.

15. Nougairede A, De Fabritus L, Aubry F, Gould EA, Holmes EC, de Lamballerie X. Random codon re-encoding induces stable reduction of replicative fitness of Chikungunya virus in primate and mosquito cells. PLoS pathogens. 2013; 9(2):e1003172. Epub 2013/02/26. https://doi.org/10.1371/journal.ppat.1003172 PMID: 23436995; PubMed Central PMCID: PMCPmc3578757.

16. Reed LJAMH. A simple method of estimating fifty per cent endpoint. American Journal of Epidemiology. 1938; 27 (3):493-7.

17. Rothberg JM, Hinz W, Rearick TM, Schultz J, Mileski W, Davey M, et al. An integrated semiconductor device enabling non-optical genome sequencing. Nature. 2011; 475(7356):348-52. https://doi.org/10.1038/nature10242 PMID: 21776081.

18. Aubry F, Nougairede A, de Fabritus L, Piorkowski G, Gould EA, de Lamballerie X. "ISA-Lation" of Single-Stranded Positive-Sense RNA Viruses from Non-Infectious Clinical/Animal Samples. PloS one. 2015; 10(9):e0138703. Epub 2015/09/26. https://doi.org/10.1371/journal.pone.0138703 PMID: 26407018; PubMed Central PMCID: PMCPmc4583506.

19. Team RC. R: A Language and Environment for Statistical Computing. 2013.