key: cord-0004186-cr0zskcw authors: Lu, Chien-Yi; Lin, Chen-Sheng; Lai, Hsueh-Chou; Yu, Ya-Wen; Liao, Chih-Yi; Su, Wen-Chi; Ko, Bo-Han; Chang, Young-Sheng; Huang, Su-Hua; Lin, Cheng-Wen title: The Rescue and Characterization of Recombinant, Microcephaly-Associated Zika Viruses as Single-Round Infectious Particles date: 2019-10-31 journal: Viruses DOI: 10.3390/v11111005 sha: 05f1c02ecdb1798e85e750ddfa81ede861f94cf6 doc_id: 4186 cord_uid: cr0zskcw Zika virus (ZIKV) is transmitted by Aedes mosquitoes and exhibits genetic variation with African and Asian lineages. ZIKV Natal RGN strain, an Asian-lineage virus, has been identified in brain tissues from fetal autopsy cases with microcephaly and is suggested to be a neurotropic variant. However, ZIKV Natal RGN strain has not been isolated; its biological features are not yet illustrated. This study rescued and characterized recombinant, single-round infectious particles (SRIPs) of the ZIKV Natal RGN strain using reverse genetic and synthetic biology techniques. The DNA-launched replicon of ZIKV Natal RGN was constructed and contains the EGFP reporter, lacks prM-E genes, and replicates under CMV promoter control. The peak in the ZIKV Natal RGN SRIP titer reached 6.25 × 10(6) TCID50/mL in the supernatant of prM-E-expressing packaging cells 72 h post-transfection with a ZIKV Natal RGN replicon. The infectivity of ZIKV Natal RGN SRIPs has been demonstrated to correlate with the green florescence intensity of the EGFP reporter, the SRIP-induced cytopathic effect, and ZIKV’s non-structural protein expression. Moreover, ZIKV Natal RGN SRIPs effectively self-replicated in rhabdomyosarcoma/muscle, glioblastoma/astrocytoma, and retinal pigmented epithelial cells, displaying unique cell susceptibility with differential attachment activity. Therefore, the recombinant ZIKV Natal RGN strain was rescued as SRIPs that could be used to elucidate the biological features of a neurotropic strain regarding cell tropism and pathogenic components, apply for antiviral agent screening, and develop vaccine candidates. Zika virus (ZIKV) was first isolated in monkeys in Uganda in 1947, and then detected in Aedes africanus mosquitoes in 1948. ZIKV belongs to a mosquito-borne flavivirus of the family Flavivirus, which spread from Africa to south-eastern Asia through transmission by Aedes mosquitoes, such as A. to produce single-round infectious particles (SRIPs). The packaging cells carrying the recombinant plasmid encoding the viral structural proteins are transfected with DNA-launched replicons and then self-replicate viral subgenomes and express all viral proteins, allowing viral assembly and release as SRIPs. Several flavivirus, replicon-based SRIPs, including JEV, dengue virus, West Nile virus, and tick-borne encephalitis virus have been generated and demonstrated as safer vaccine candidates [15] [16] [17] . Subsequently, replicon-based SRIPs are suitable systems for investigating the mechanism of the ZIKV life cycle, ZIKV-host interaction, and virulence. The ZIKV Natal RGN strain that is associated with microcephaly was recognized as the representative, epidemic, ZIKV Asian strain in 2015-2016 [18] . The ZIKV Natal RGN strain is detected in brain tissues from fetal autopsy cases with microcephaly in the Natal region of Brazil; its genome has been sequenced using next-generation sequencing (GenBank accession number KU527068) [15] . The ZIKV Natal RGN strain contains a unique genotype and phenotype; therefore, this study aimed to generate SRIPs of the ZIKV Natal RGN strain using synthetic and reverse genetic technologies. This paper characterizes the infectivity of ZIKV Natal RGN SRIPs in different cell types. Moreover, the paper analyzes the expression pattern of type I interferons and apoptosis-related genes induced by ZIKV Natal RGN SRIPs. TE671 (human, rhabdomyosarcoma/muscle), SF268 (human brain, glioblastoma/astrocytoma), and ARPE-19 (human, retinal pigmented epithelium) cells were cultured in minimum essential medium containing 10% fetal bovine serum, 2 mM glutamine, 1 mM pyruvate, and 1× penicillin-streptomycin at 37 • C with 5% CO 2 . Components of two DNA segments we synthesized, I and II, consisted of the CMV immediate-early promoter (CMVp), cDNA fragments of the ZIKV Natal RGN genome (GenBank accession number KU527068), enhanced green fluorescent protein (EGFP), FMDV-2A (F-2A), hepatitis delta virus ribozyme (HDVr), and bovine growth hormone polyA signal (BGH-pA), which were purchased from Bio Basic Inc. (Ontario, Canada). The construction of synthesized DNA fragments, plasmid information, and restriction enzyme digestion analysis of the DNA fragments we purchased are shown in Supplementary Figures S1 and S2 and Figure 1A . To construct the ZIKV Natal RGN replicon, the plasmid pBR322-Linker, described in our prior report [15] , was chosen as the vector to assemble fragments A, B, C, and D of the CMVp-driven ZIKV Natal RGN replicon with the EGFP reporter via the cloning sites of EcoRI, NotI, ClaI, RsrII, and XhoI ( Figure 1B ). Fragments A, C, and D were generated using the Platinum ® PCR Supermix High Fidelity reaction with the DNA segments we synthesized, I and II, as the templates (Life Technology, Carlsbad, CA, USA). Fragment B was produced through two-round PCR with the Gibson Assembly reaction of two first-round PCR products (B1 and B2) as the template (New England Biolabs, Ipswich, MA, USA). The pairs containing the restriction enzyme site(s) indicated, are listed in Supplementary Table S1 . The resultant plasmid carrying the CMVp-driven ZIKV Natal RGN replicon with the EGFP reporter was propagated in E. coli DH5α and then sequenced by Sanger sequencing assays with sequencing primers (Supplementary Table S2 ). The nucleotide and deduced amino acid sequence alignment analysis of the CMVp-driven ZIKV Natal RGN replicon and its parent strain (GenBank accession number KU527068) was performed using the Lasergene DNASTAR Megalign software. Figure 1 . Construction of the pBR322-based ZIKV Natal RGN replicon and pcDNA3.1-ZIKV prME. Two synthetic DNA segments in the pUC18 plasmid contained the entire ZIKV Natal RGN strain genome: CMVp, EGFP, FMDV-2A, and BGH-pA sequences (A). Four PCR products (Fragments A-D) were cloned into the indicated restriction sites (EcoRI, NotI, ClaI, RsrII, and XhoI) of the pBR322 plasmid and then assembled as the in-frame fusion of the ZIKV Natal RGN replicon with the EGFP reporter under CMV-promoter control (B). Those four PCR products were analyzed using agarose gel electrophoresis (C). The PCR product of ZKIV prM and E genes was cloned into restriction sites (EcoRI and XhoI) of the pcDNA3.1-His-C plasmid (D). Initially, TE671 cells were transfected with pcDNA3.1-ZIKV prM-E, selected with G418 for 2 weeks, analyzed to examine ZIKV prM and E protein expression, and then recognized as the packaging cell line ( Figures 1D,2A middle, 3A middle). Subsequently, the packaging cells were transfected with the ZIKV Natal RGN replicon, and then CPE and the EGFP reporter expressions were examined ( Figure 2 ). In the packaging cells, EGFP reporter expression and CPE showed slight levels 24 h post-transfection, but exhibited significant changes 72 h after transfection with the ZIKV Natal RGN replicon (Figure 2 bottom) . Moreover, in vitro translation and transcription of the ZIKV Natal RGN replicon in packaging cells were further characterized ( Figure 3 ). Immunofluorescent staining indicated that the ZIKV E C-terminus, NS1, and NS5 proteins were massively expressed in replicon-transfected packaging cells ( Figure 3A bottom), but not in mock and packaging cells ( Figure 3A top and middle). The real-time, reverse transcription PCR assay showed large amounts of positive and negative-sense ZIKV subgenomes in replicon-transfected packaging cells 72 h post-transfection, but not in mock cells ( Figure 3B top and bottom) . The results demonstrated efficient translation and replication of the ZIKV subgenomes in the packaging cells using the DNA-launched ZIKV Natal RGN replicon. Figure 1 . Construction of the pBR322-based ZIKV Natal RGN replicon and pcDNA3.1-ZIKV prME. Two synthetic DNA segments in the pUC18 plasmid contained the entire ZIKV Natal RGN strain genome: CMVp, EGFP, FMDV-2A, and BGH-pA sequences (A). Four PCR products (Fragments A-D) were cloned into the indicated restriction sites (EcoRI, NotI, ClaI, RsrII, and XhoI) of the pBR322 plasmid and then assembled as the in-frame fusion of the ZIKV Natal RGN replicon with the EGFP reporter under CMV-promoter control (B). Those four PCR products were analyzed using agarose gel electrophoresis (C). The PCR product of ZKIV prM and E genes was cloned into restriction sites (EcoRI and XhoI) of the pcDNA3.1-His-C plasmid (D). To detect self-amplifying RNA genomes of the CMVp-driven ZIKV Natal RGN replicon, the synthesis of positive and negative-sense RNA subgenomes in vitro was examined using SYBR Green-based real-time PCR. Total RNAs of TE-671 cells transfected with the ZIKV Natal RGN replicon were extracted using the PureLink Mini Total RNA Purification Kit (Thermo Fisher Scientific, Waltham, MA, USA), reverse transcribed into cDNA with specific-capture primers, and measured using real-time PCR with ZIKV NS5-specific primer pairs (Supplementary Table S3 ), according to our prior report [19] . Relative levels of self-amplifying ZIKV sense and antisense genomes were normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The absolute copy number of the ZIKV genome was determined according to the standard curve of real-time PCR for serial dilutions of the plasmid containing the ZIKV NS5 gene at a known concentration. To explore the expression of replicon-based EGFP and ZIKV reporter proteins, replicon-transfected cells were initially photographed using immunofluorescence microscopy; then, an immunofluorescence assay (IFA) was performed with primary antibodies against ZIKV NS1 and NS5 (GeneTex, Inc., Taiwan) and secondary AF546 goat anti-rabbit IgG (Thermo Fisher Scientific). The immunofluorescent staining assay was carried out as described in our prior report [15] . The fluorescence intensity of replicon-based EGFP and ZIKV proteins in ZIKV Natal RGN replicon-transfected cells was counted by Image J. To establish the packaging cells expressing ZIKV structural proteins, the prM/M-E gene was amplified using PCR with specific primers (Supplementary Table S1 ) and synthesized DNA segment I as the template. The PCR product of the ZIKV prM-E gene was digested with EcoRI and XhoI, and then cloned into the EcoRI and XhoI sites of the expression plasmid pcDNA3.1-HisC. The resultant plasmid pcDNA-ZIKV prME was transfected with TE-671 cells at 90% confluence in a 6-well plate with Lipofectamine LTX (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's guidelines. The transfected cells were selected in the culture media with 500 µg/mL G418 for 2 weeks; the expression of ZIKV prM and E proteins in a stably transfected cell line (packaging cell line) was validated by real-time RT-PCR with ZIKV E-specific primers (Supplementary Table S3 ) and immunofluorescence staining with primary antibodies against ZIKV E protein (GeneTex, Inc.), as described above. The packaging cells grown at 90% confluence in 6-well plates were transfected with or without the ZIKV Natal RGN replicon using Lipofectamine LTX. The cytopathic effect (CPE), replicon-based EGFP and ZIKV proteins, and viral sense and antisense genomes in mock and transfected cells, were measured 24 and 72 h post-transfection using immunofluorescent staining assays, as described above. Moreover, real-time RT-PCR assays with NS5-specific primer pairs, listed in Supplementary Table S3 , were performed to examine the positive and negative-sense viral genomes in replicon-transfected packaging cells, as well as the viral genome in SRIPs. ZIKV Natal RGN SRIPs in the cultured media of transfected packaging cells were concentrated using the PEG Virus Precipitation Kit (Abcam, Cambridge, MA, USA). The cultured media were centrifuged at 3200× g for 15 min at 4 • C; then, the supernatant (10 mL) was collected from each and incubated with 2.5 mL PEG solution overnight at 4 • C. The pellets of ZIKV Natal RGN SRIPs were harvested after centrifugation at 3200× g for 30 min at 4 • C, re-suspended in 100 µL Virus Re-suspension Solution, and then stored at −80 • C. To analyze the antigenicity of ZIKV Natal RGN SRIPs, 10 µL each of one-fold and 10-fold SRIP stock dilution was spotted onto a nitrocellulose membrane for the dot-blot assay with anti-ZIKV E antibodies. The membrane was subsequently blocked with 5% skim milk in TBST (Tris-Buffered Saline, 0.1% Tween-20) buffer at 4 • C for 2 h, incubated overnight with anti-ZKIV E antibodies (GeneTex, Inc., Taiwan), and then reacted with HRP-conjugated anti-mouse IgG antibodies (Invitrogen, Carlsbad, CA, USA) after washing with TBST. Immunoreactive signals for the E proteins of ZIKV Natal RGN SRIPs were amplified with ECLTM Western Blotting Detection Reagents (GE Healthcare, Chicago, IL, USA), and then imaged by the Multi-function Gel Image System (MultiGel-21) (Gentaur, San Jose, CA, USA). The infectious titer of the ZIKV Natal RGN SRIP stock was determined by a median tissue culture infectious dose (TCID50) assay. Serial dilutions of the ZIKV Natal RGN SRIP stock were added and incubated on the 90% confluent monolayer of packaging cells in 96-well plates. After incubating for 72 h at 37 • C, the CPE in each well was observed and recorded to determine the TCID50 titer of the ZIKV Natal RGN SRIP stock. In the infectivity assay, packaging cells were cultured in 6-well plates overnight and infected with different doses of ZIKV Natal RGN SRIPs at multiplicity of infections (MOIs) of 3, 0.3, 0.03, and 0.003. CPE, replicon-driven EGFP, and ZIKV NS1 expression in the SRIP-infected cells were examined 72 h post-infection, and fluorescence microscopy and immunofluorescent staining were performed as described above. In cell line susceptibility assays, TE-671, SF-268, and ARPE-19 cells were infected with ZIKV Natal RGN SRIPs at a MOI of 0.15 TCID50/cell. After incubating for 72 h, the cells were photographed to examine CPE and explore replicon-driven EGFP reporter expression using immunofluorescence microscopy. Moreover, the cells were fixed with paraformaldehyde at 4 • C for 1 h, and then immunofluorescence staining was performed with anti-ZIKV NS1 antibodies, as mentioned above. In the attachment assay, ZIKV Natal RGN SRIPs at a dose of 10 5 TCID50/well (MOI = 0.5 TCID50/cell) were added onto the monolayers of ARPE-19, TE671, and SF-268 cells and incubated at 4 • C for 1 h. After washing with cold PBS, ZIKV Natal RGN SRIPs attached on the cell surface were harvested, and the relative levels of extracted viral genomes were quantitated using two-step, SYBR, Green I, real-time RT-PCR with ZIKV NS5-specific primer pairs, as described above. All data from 3 independent experiments were analyzed using Student's t-tests or χ 2 tests. Statistical significance for each assay was judged at p < 0.05. The ZIKV Asian Natal RGN prototype strain, a microcephaly-associated ZIKV, was not isolated from the brain tissues from fetal autopsy cases; therefore, this study intended to recover the recombinant ZIKV Natal RGN strain using reverse genetics and synthetic biology techniques based on the published genome sequence (GenBank accession number KU527068). Two synthetic DNA segments in the pUC18 plasmids comprised the entire ZIKV Natal RGN strain genome, CMVp, EGFP, FMDV-2A, and BGH-pA sequences (Supplementary Figures S1 and S2 and Figure 1A ,B). The synthetic DNA segments were used as templates via PCR to amplify the gene fragments (A-D) that were cloned into the low-copy-number plasmid pBR322 to construct a DNA-launched ZIKV Natal RGN replicon under the control of the CMV promoter ( Figure 1B,C) . Fragment A consisted of CMVp, ZIKV 5'-UTR, ZIKV C protein, the reporter EGFP, and FMDV-2A. Fragment B containing the 73-residue C-terminal region of ZIKV E protein, NS1, NS2A, NS2B, and the N-terminal region of NS3, was produced using Gibson Assembly cloning ( Figure 1C ). In order to join fragments B and C, the reverse primer for fragment B and the forward primer for fragment C contained the restriction enzyme site ClaI (Supplementary Table S1 ), which is recognized as the marker nt5156 ATCGAT 5161 of the replicon in which thymidine at nucleotide 5156 of the ZIKV Natal RGN genome was exchanged with adenosine, but it caused a silent mutation ( Figure 1B) . Interestingly, the RsrII restriction site within the NS4B gene ( nt7473 CGGACCG 7479 ) of the ZIKV Natal RGN genome was a unique restriction site for ligation between fragments C and D. In addition, HDVr was fused with the 3'-UTR of ZIKV in the replicon to verify the accuracy of the 3'-end of the transcribed ZIKV Natal RGN RNA genomes ( Figure 1B,C) . Furthermore, sequencing analysis of the ZIKV Natal RGN replicon indicated that only two amino acid substitutions appeared within the NS5 protein (L3277R and E3280N) (Supplementary Figure S3) . Therefore, the DNA-launched ZIKV Natal RGN replicon was effectively constructed to generate the in-frame fusion of ZIKV C protein, reporter EGFP, FMDV-2A, the 73-residue C-terminal region of ZIKV E protein, and all non-structural ZIKV proteins under the control of the CMV promoter. Initially, TE671 cells were transfected with pcDNA3.1-ZIKV prM-E, selected with G418 for 2 weeks, analyzed to examine ZIKV prM and E protein expression, and then recognized as the packaging cell line ( Figure 1D To examine the production of ZIKV Natal RGN SRIPs, the supernatant of replicon-transfected packaging cells was harvested 72 h post-transfection and analyzed using dot-blotting, real-time RT-PCR, and TCID50 assays (Figure 4) . The dot-blotting assay with anti-ZIKV E protein demonstrated the antigenic structure of the E protein ( Figure 4A) , and real-time RT-PCR assay elucidated the viral subgenomes of ZIKV Natal RGN SRIPs from the supernatant of the replicon-transfected packaging cells ( Figure 4B) . Moreover, the titer of ZIKV Natal RGN SRIPs in the supernatant was determined by TCID50 assay with the packaging cells. The titer of the ZIKV Natal RGN SRIP stock was 6.25 × 10 6 TCID50/mL ( Figure 4C ). To examine the production of ZIKV Natal RGN SRIPs, the supernatant of replicon-transfected packaging cells was harvested 72 h post-transfection and analyzed using dot-blotting, real-time RT-PCR, and TCID50 assays (Figure 4) . The dot-blotting assay with anti-ZIKV E protein demonstrated the antigenic structure of the E protein ( Figure 4A) , and real-time RT-PCR assay elucidated the viral subgenomes of ZIKV Natal RGN SRIPs from the supernatant of the replicon-transfected packaging cells ( Figure 4B) . Moreover, the titer of ZIKV Natal RGN SRIPs in the supernatant was determined by TCID50 assay with the packaging cells. The titer of the ZIKV Natal RGN SRIP stock was 6.25 × 10 6 TCID50/mL ( Figure 4C ). To examine the infectivity of ZIKV Natal RGN SRIPs in vitro, the packaging cells were infected with 10-fold dilutions (MOIs of 3, 0.3, 0.03, and 0.003) of ZIKV Natal RGN SRIPs to determine SRIPinduced CPE and SRIP-driven expression of EGFP reporter and ZIKV NS1 ( To examine the cell line susceptibility to ZIKV Natal RGN SRIPs in vitro, the three cell lines TE-671, SF-268, and ARPE-19 were used ( Figure 6 ). ZIKV Natal RGN SRIPs at a MOI of 0.15 TCID50/cell induced more observable and higher CPE levels in TE-671 cells than those in ARPE-19 and SF268 cells (Figure 6A top) . The CPE level induced by ZIKV Natal RGN SRIPs correlated with the fluorescence intensity of the EGFP reporter encoded within the SRIP genome in the following order: TE-671 > ARPE-19 > SF268 cells infected with ZIKV Natal RGN SRIPs ( Figure 6A middle) . In addition, immunofluorescent staining with anti-ZIKV NS1 antibodies indicated that ZIKV NS1 protein was expressed in these three ZIKV Natal RGN SRIP-infected cell lines, in which the fluorescence intensity of the NS1 protein was consistent with the CPE level induced by ZIKV Natal RGN SRIPs ( Figure 6A bottom). Meanwhile, 10 5 TCID50 of the SRIPs per well was added onto the monolayer of these three cell lines to examine the attachment activity of ZIKV Natal RGN SRIPs ( Figure 6B) . Interestingly, the SRIPs showed the highest level of attachment activity to TE671 cells, in which the highest viral genome number was discovered 1 h after attachment at 4 °C. However, the SRIPs showed the lowest level of attachment activity to SF268 cells. Compared to ARPE-19 and SF268 cells, it was found that TE671 cells exhibited the highest susceptibility to ZIKV Natal RGN SRIPs and allowed the highest number of SRIPs to attach to the surface. To examine the cell line susceptibility to ZIKV Natal RGN SRIPs in vitro, the three cell lines TE-671, SF-268, and ARPE-19 were used ( Figure 6 ). ZIKV Natal RGN SRIPs at a MOI of 0.15 TCID50/cell induced more observable and higher CPE levels in TE-671 cells than those in ARPE-19 and SF268 cells (Figure 6A top) . The CPE level induced by ZIKV Natal RGN SRIPs correlated with the fluorescence intensity of the EGFP reporter encoded within the SRIP genome in the following order: TE-671 > ARPE-19 > SF268 cells infected with ZIKV Natal RGN SRIPs ( Figure 6A middle) . In addition, immunofluorescent staining with anti-ZIKV NS1 antibodies indicated that ZIKV NS1 protein was expressed in these three ZIKV Natal RGN SRIP-infected cell lines, in which the fluorescence intensity of the NS1 protein was consistent with the CPE level induced by ZIKV Natal RGN SRIPs ( Figure 6A bottom). Meanwhile, 10 5 TCID50 of the SRIPs per well was added onto the monolayer of these three cell lines to examine the attachment activity of ZIKV Natal RGN SRIPs ( Figure 6B) . Interestingly, the SRIPs showed the highest level of attachment activity to TE671 cells, in which the highest viral genome number was discovered 1 h after attachment at 4 • C. However, the SRIPs showed the lowest level of attachment activity to SF268 cells. Compared to ARPE-19 and SF268 cells, it was found that TE671 cells exhibited the highest susceptibility to ZIKV Natal RGN SRIPs and allowed the highest number of SRIPs to attach to the surface. The ZIKV Natal RGN strain found in brain tissues from fetal autopsy cases with microcephaly in the Natal region of Brazil has not been isolated [18] . This study, firstly, generated SRIPs of the ZIKV Natal RGN strain using synthetic and reverse genetic technologies (Figures 1-4 and Supplementary Figures S1-S3) . The viral yield of ZIKV Natal RGN SRIPs in prM-E-expressing packaging cells transfected with the ZIKV replicon was 6.25 × 10 6 TCID50/mL (Figure 4) , which indicated the replication capacity of the CMV promoter-launched ZIKV Natal RGN replicon. In addition, a green fluorescence reporter, EGFP, encoded within the ZIKV Natal RGN replicon, represented an observable marker to examine replication of the ZIKV Natal RGN replicon and SRIPs (Figures 2,3, and 5 ). In our laboratory, the pBR322 plasmid has been used to construct a JEV replicon that is genetically stable in E. coli [15] . The JEV genome sequence encoding for structural proteins (C, prM/M, and E) was toxic in E. coli, which resulted in genetic instability of the JEV infectious clones and replicons [15] . Similarly, construction of ZIKV Natal RGN infectious clones was not accomplished due to the genetic instability of ZIKV Natal RGN infectious clones, particularly within the NS1 gene region (Supplementary Figures S1 and S2) . However, the recombinant pcDNA3.1 plasmid containing prM and E genes and the pBR322 plasmid containing 5' and 3'-UTRs, C, and all NS genes, showed genetic stability. Two ZIKV infectious clones obtained by cloning the entire fulllength native cDNA have been reported [20, 21] , and other ZIKV infectious clones were reported using strategies to avoid the toxicity of the ZIKV genome sequence, including in vitro assembly of viral cDNA fragments and insertion of eukaryotic introns [22] [23] [24] [25] [26] [27] . Particularly, the ZIKV Natal RGN infectious clone was produced using bacterial artificial chromosomes (BACs) [27] , in which a singlecopy pBeloBAC11 plasmid containing the full-length genome of the ZIKV Natal RGN strain under the CMV promoter control was used to maintain ZIKV genome stability and reduce toxicity in bacteria. In addition, our prior report demonstrated that JEV SRIPs carrying the EGFP reporter were a more sensitive, effective, and efficient platform for quantitating antiviral efficacy using flow cytometry compared to the traditional plaque assay [15] . Like JEV SRIPs, the ZIKV Natal RGN SRIPs The ZIKV Natal RGN strain found in brain tissues from fetal autopsy cases with microcephaly in the Natal region of Brazil has not been isolated [18] . This study, firstly, generated SRIPs of the ZIKV Natal RGN strain using synthetic and reverse genetic technologies (Figures 1-4 and Supplementary Figures S1-S3) . The viral yield of ZIKV Natal RGN SRIPs in prM-E-expressing packaging cells transfected with the ZIKV replicon was 6.25 × 10 6 TCID50/mL (Figure 4) , which indicated the replication capacity of the CMV promoter-launched ZIKV Natal RGN replicon. In addition, a green fluorescence reporter, EGFP, encoded within the ZIKV Natal RGN replicon, represented an observable marker to examine replication of the ZIKV Natal RGN replicon and SRIPs ( Figure 2, Figure 3 , and Figure 5 ). In our laboratory, the pBR322 plasmid has been used to construct a JEV replicon that is genetically stable in E. coli [15] . The JEV genome sequence encoding for structural proteins (C, prM/M, and E) was toxic in E. coli, which resulted in genetic instability of the JEV infectious clones and replicons [15] . Similarly, construction of ZIKV Natal RGN infectious clones was not accomplished due to the genetic instability of ZIKV Natal RGN infectious clones, particularly within the NS1 gene region ( Supplementary Figures S1 and S2) . However, the recombinant pcDNA3.1 plasmid containing prM and E genes and the pBR322 plasmid containing 5' and 3'-UTRs, C, and all NS genes, showed genetic stability. Two ZIKV infectious clones obtained by cloning the entire full-length native cDNA have been reported [20, 21] , and other ZIKV infectious clones were reported using strategies to avoid the toxicity of the ZIKV genome sequence, including in vitro assembly of viral cDNA fragments and insertion of eukaryotic introns [22] [23] [24] [25] [26] [27] . Particularly, the ZIKV Natal RGN infectious clone was produced using bacterial artificial chromosomes (BACs) [27] , in which a single-copy pBeloBAC11 plasmid containing the full-length genome of the ZIKV Natal RGN strain under the CMV promoter control was used to maintain ZIKV genome stability and reduce toxicity in bacteria. In addition, our prior report demonstrated that JEV SRIPs carrying the EGFP reporter were a more sensitive, effective, and efficient platform for quantitating antiviral efficacy using flow cytometry compared to the traditional plaque assay [15] . Like JEV SRIPs, the ZIKV Natal RGN SRIPs carrying the EGFP reporter could be applicable as a drug screening platform. Therefore, the reverse genetics system for the production of ZIKV Natal RGN SRIPs provided an alternative, applicable, and rapid approach to studying the biological features of the unique ZIKV strains. The human rhabdomyosarcoma/muscle TE671 cells were used as packaging cells expressing ZIKV prM and E proteins (Figure 2 ), in which transfection with the ZIKV Natal RGN replicon significantly generated recombinant SRIPs with the high yield of 6.25 × 10 6 TCID50/mL (Figure 4) . Interestingly, TE671 cells expressing JEV C, prM, and E proteins were also used as a packaging cell line that was efficient at producing JEV SRIPs post-transfection with the JEV replicon [15] . The recombinant ZIKV Natal RGN SRIPs showed infectivity and self-replication in the packaging cells ( Figure 5 ). Two amino acid substitutions within NS5 (L3277R and E3280N) (Supplementary Figure S3) seemed to have no influence on the activities of NS5 methyltransferase and RNA-dependent RNA polymerase. A recent study indicated that human embryonic kidney 293T cells had been used to establish different flavivirus prM-E-expressing packaging cells, including DENV1-4, ZIKV, JEV, West Nile virus, yellow fever virus, and tick-borne encephalitis virus, and the cells were transfected with the DENV1 reporter replicon to generate a chimeric flavivirus SRIP-based neutralization assay [28] . Therefore, recombinant ZIKV Natal RGN SRIPs generated in this study could be used to analyze the ZIKV life cycle, elucidate viral pathogenesis, screen antiviral agents, design faster and safer diagnostic tools, and even develop attenuated vaccines against ZIKV infection. The full-length ZIKV Natal RGN strain genome was directly sequenced using next-generation sequencing with the total RNA extracted from brain tissues from fetal autopsy cases with microcephaly (GenBank accession number KU527068) [18] . In this study, recombinant ZIKV Natal RGN SRIPs were rescued using reverse genetic and synthetic biology techniques (Figures 1-5 and Supplementary Figures S1-S3) . Three cell lines, TE671, SF268, and ARPE-19, were further tested for cell susceptibility to ZIKV Natal RGN SRIPs ( Figure 6 ). The order of cell susceptibility to ZIKV Natal RGN SRIPs was TE671 > ARPE-19 > SF268 cells. The CPE level, EGFP reporter intensity, and ZIKV NS1 protein expression were most obvious (>50% involvement) in TE671 cells, but were mild in ARPE-19 cells and limited in SF268 cells 72 h post-infection with SRIPs at an MOI of 0.15 TCID50/cell. The results indicated that ZIKV Natal RGN SRIPs had a broad range of human tissue tropism with differential cell susceptibilities ( Figure 6A ). In addition, ZIKV Natal RGN SRIPs had the highest attachment activity to TE671 cells among those three cell lines ( Figure 6B ). SRIP attachment to TE671 cells was suggested to be involved in the infectivity of ZIKV Natal RGN in TE671 cells. This finding implied that ZIKV Natal RGN SRIPs had tropism for rhabdomyosarcoma/muscle and retinal pigmented epithelial cells, which was consistent with myalgia, arthralgia, retinitis, chorioretinal atrophy, and pigmentary maculopathy in symptomatic ZIKV infection [29, 30] . Therefore, the cell line susceptibility results revealed the unique features of ZIKV Natal RGN SRIPs and suggested the possible pathogenesis of ZIKV infection. In conclusion, the approach for generating recombinant ZIKV SRIPs with an EGFP reporter showed the genetic stability of ZIKV prM-E genes in the mammalian expression vector and ZIKV replicon sequence in the low-copy plasmid pBR322, providing a faster and safer platform to study the biological aspects of ZIKV infection. ZIKV Natal RGN SRIPs were rescued by the combination of prM-E-expressing packaging cells and the ZIKV Natal RGN replicon under CMV promoter control, which mimicked self-replicating, positive-sense, single strand RNA viruses. The biological features of ZIKV Natal RGN SRIPs were revealed, including CPE, infectivity, cell susceptibility, and attachment activity. Therefore, recombinant ZIKV Natal RGN SRIPs could be further used to investigate cell tropism, persistent infection, and vaccine candidates, and to establish reliable, effective, and efficient assays for screening antiviral agents and diagnosing viral infections. Table S3 . Primer pairs for real-time quantitative RT-PCR of ZIKV E and NS5 RNA copies; Supplemental Figure S1 . The construction of in-frame gene fusion (A), plasmid information and restriction enzyme digestion analysis (B) of the synthesized DNA segment I; Supplemental Figure S2 . The construction of in-frame gene fusion (A), plasmid information and restriction enzyme digestion analysis (B) of the synthesized DNA segment II; Supplemental Figure S3 . Sequencing analysis of ZIKV Asian Natal RGN replicon with the primers listed in Supplemental Table S2A . The amino acid substitutions were also shown in the Figure. The Diversification of Zika Virus: Are There Two Distinct Lineages? Molecular evolution of Zika virus during its emergence in the 20(th) century Evolution of Two Major Zika Virus Lineages: Implications for Pathology, Immune Response, and Vaccine Development Virus Pathogenesis and Tissue Tropism The origin and spread of a mosquito-borne virus Zika Virus-Associated Neurological Disease in the Adult: Guillain-Barré Syndrome, Encephalitis, and Myelitis An update on Zika virus infection Comparative genomic analysis of pre-epidemic and epidemic Zika virus strains for virological factors potentially associated with the rapidly expanding epidemic An Evolutionary Insight into Zika Virus Strains Isolated in the Latin American Region A single mutation in the prM protein of Zika virus contributes to fetal microcephaly An evolutionary NS1 mutation enhances Zika virus evasion of host interferon induction Development of Dengue type-2 virus replicons expressing GFP reporter gene in study of viral RNA replication Development and characterization of the replicon system of Japanese encephalitis live vaccine virus SA14-14-2 Noncytopathic flavivirus replicon RNA-based system for expression and delivery of heterologous genes Single-Round Infectious Particle Antiviral Screening Assays for the Japanese Encephalitis Virus Production of single-round infectious chimeric flaviviruses with DNA-based Japanese encephalitis virus replicon Increased expression of capsid protein in trans enhances production of single-round infectious particles by West Nile virus DNA vaccine candidate SARS coronavirus papain-like protease induces Egr-1-dependent up-regulation of TGF-β1 via ROS/p38 MAPK/STAT3 pathway An infectious cDNA clone of zika virus to study viral virulence, mosquito transmission, and antiviral inhibitors Novel genetically stable infectious clone for a Zika virus clinical isolate and identification of RNA elements essential for virus production A robust method for the rapid generation of recombinant Zika virus expressing the GFP reporter gene Development and characterization of recombinant virus generated from a new world Zika virus infectious clone A reverse genetics platform that spans the Zika virus family tree Characterization of cis-acting RNA elements of Zika virus by using a self-splicing ribozyme-dependent infectious clone A full-length infectious cDNA clone of Zika virus from the 2015 epidemic in Brazil as a genetic platform for studies of virus-host interactions and vaccine development Rescue of Recombinant Zika Virus from a Bacterial Artificial Chromosome cDNA Clone High-throughput neutralization assay for multiple flaviviruses based on single-round infectious particles using dengue virus type 1 reporter replicon Rescue of the 1947 Zika virus prototype strain with a cytomegalovirus promoter-driven cDNA clone Risk factors associated with the ophthalmoscopic findings identified in infants with presumed Zika virus congenital infection All authors declare no potential conflict of interest.