Zika virus infects renal proximal tubular epithelial cells with prolonged persistency and cytopathic effects

Keywords: cytopathic effects, prolonged persistency, renal proximal tubular epithelial cells, urinary system, Zika virus

INTRODUCTION

A Zika virus (ZIKV) infection is typically accompanied by mild symptoms that are present in 20% infected individuals.¹ However, the 2015 ZIKV outbreak in Brazil was associated with microcephaly, a severe developmental abnormality of the fetal brain.² In vivo studies of mice^{3, 4, 5, 6} and non-human primates⁷ have established a causal link between ZIKV infection and the birth defects of newborns. In vitro studies have identified neural progenitor cells (NPCs) and neural stem cells (NSCs) as the permissive host cells of ZIKV in the central nervous system, providing insights into the mechanism of ZIKV-induced fetal demise.^{8, 9, 10, 11} In addition, ZIKV can infect and cause damage in the human peripheral nervous system, which may lead to Guillain–Barré syndrome (GBS) in adults.^12, 13

Considerable efforts have been devoted to understanding the viral reservoir of ZIKV in human due to its importance in the prevention and treatment of ZIKV infections. It has been reported that sexual transmission of ZIKV can occur from males to females,^14, 15, 16 males to males¹⁷ and, in one suspected case, a female to a male.¹⁸ ZIKV may also be transmitted from sexual partners without ZIKV-related symptoms.¹⁶ ZIKV can exist in human semen^19, 20 and in the genital tract.^21, 22, 23 In support of these clinical findings, studies using mouse models have shown that ZIKV can infect and replicate in the vagina^24, 25 and that infection of testes led to testis damage and male infertility.^26, 27 These studies uncovered a ZIKV reservoir in the reproductive system apart from that in the nervous system. ZIKV has been reported to exist in the urine of adults^{28, 29, 30, 31, 32} and in neonates²⁹ as well as in the amniotic fluids (composed primarily of fetal urine) of ZIKV-infected fetuses.^33, 34 Interestingly, although in most patients ZIKV RNA cannot be detected in the sera after the first week of illness,^35, 36 their urine samples have detectable levels of ZIKV RNA until at least 2 weeks after the onset of symptoms.^{30, 36, 37, 38} The United States Centers for Disease Control and Prevention (CDC) has recommended the use of urine samples for the non-invasive diagnosis of ZIKV.³⁹ Nevertheless, it is still debatable whether urine should be used for the molecular diagnosis of ZIKV.^40, 41, 42 This uncertainty is, in part, due to the unidentified source of ZIKV in the urine samples. Because urinary ZIKV typically exhibits a higher viral load and longer duration than that in the sera,³² the urinary system may present an independent viral reservoir that is different from the bloodstream.

Renal infections caused by ZIKV are largely unknown. However, a previous study showed that ZIKV could infect renal glomerular cells.⁴³ In this study, we showed that ZIKV can infect renal tubular cells in an immunodeficient mouse model. In addition, we found that the in vitro ZIKV infection of human RPTEpiCs (hRPTEpiCs) resulted in the prolonged persistency of viruses. Both in vivo and in vitro ZIKV infections resulted in caspase-3-induced cell apoptosis. RNA-Seq analyses of ZIKV-infected hRPTEpiCs revealed a notable transcriptional change of a large set of genes, including type I interferon signaling and anti-viral response genes.

Go to:

MATERIALS AND METHODS

Cell lines and viruses

Aedes albopictus C6/36 cells were grown in 30% RPMI-1640 (Gibco-Life Technologies, Carlsbad, CA, USA), 60% Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 10% fetal bovine serum (FBS, Gibco) at 28 °C, as previously described.⁴⁴ HK2 cells were cultured in DMEM/F12 (1:1, Gibco) supplemented with 10% FBS, 100 IU/mL of penicillin and 100 μg/mL of streptomycin and were maintained at 37 °C in a fully humidified atmosphere with 5% CO₂. Primary hRPTEpiCs were purchased from ScienCell (Carlsbad, CA, USA) and maintained on polylysine (ScienCell)-coated plates in Epithelial Cell Medium (EpiCM, ScienCell) supplemented with 2% FBS (ScienCell), 1% epithelial cell growth supplement (EpiCGS, ScienCell) and 1% penicillin/streptomycin solution (P/S, ScienCell). The SZ01 ZIKV stock⁴⁵ was kindly provided by Professor Cheng-Feng Qin at the Department of Virology, State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology. The MR766 ZIKV stock was obtained from ATCC (VR-1838; Manassas, VA, USA).

Animal experiments

All animal experiments described in this study were approved by Fudan University and ShanghaiTech University. Animal studies were carried out in strict accordance with the Institutional Animal Care and Use Committee approved protocols. Inoculation of animals with ZIKV was performed under anesthesia using chloral hydrate, and all efforts were made to minimize animal suffering.

C57BL/6 mice deficient of type I and II interferon (IFN) receptors (AG6 mice) were purchased from the B&K Universal Group Limited (Shanghai, China) and housed under specific pathogen-free (SPF) conditions at the animal facilities of the Shanghai Public Health Clinical Center, Fudan University (Shanghai, China). The mice were transferred to the Animal Biosafety Level 2 Laboratory before infection. Groups of 4- to 8-week-old mice of both sexes were used for all experiments. All mice were intraperitoneally (i.p.) inoculated with 10⁵ plaque-forming unit (PFU) of ZIKV in a 100 μL volume. Kidneys were harvested from both mock and infected mice at 7 days post infection and were immediately fixed for 16 h in 10% neutral buffered formalin (NBF). Fixed kidneys were mounted and sectioned into 10-μm slices for RNA in situ hybridization (ISH) and hematoxylin and eosin (H&E) staining.

RNA ISH was performed using a RNAscope 2.5 HD Reagent Kit-BROWN (Advanced Cell Diagnostics, Newark, CA, USA) according to the manufacturer’s instructions. NBF-fixed tissue slides were hydrated and then immersed successively in 200 mL of 50% ethanol, 200 mL of 70% ethanol and 400 mL of 100% ethanol for 5 min each at room temperature. ZIKV genomic RNA was detected using an RNA probe from Advanced Cell Diagnostics (Cat No. 467771). Hybridization signals were visualized by chromogenic reactions using DAB chromogen, followed by counterstain with Gill’s hematoxylin. RNA ISH and H&E stain samples were photographed using a Tissue FAXS 200 flow-type quantitative tissue analyzer (TissueGnostics GmbH, Vienna, Austria).

Quantification of ZIKV

The ZIKV viral RNA in the supernatant of medium was quantified by RT-qPCR. Viral RNA was extracted using a QIAamp viral RNA mini kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. One microgram of RNA was transcribed into cDNA using random primers using Moloney murine leukemia virus (M-MLV) reverse transcriptase (Promega, Charbonnieres, France) and the cDNA product was amplified by RT-qPCR using primers SZ01-1F (CAA GGA GTG GGA AGC GGA G) and SZ01-1R (CCA TGT GAT GTC ACC TGC TCT) and the T7 RiboMAX Express RNAi System (Promega) on an Applied Biosystems 7300 real-time PCR cycler (Applied Biosystems, Carlsbad, CA, USA). qPCR data were analyzed using SDS software from Applied Biosystems and the viral RNA was quantified by comparing each sample’s threshold cycle (CT) value with a ZIKV RNA standard curve. The TCID₅₀ of each collected supernatant was determined by titrating HK2 cells in a 96-well plate.

Immunofluorescence and immunohistochemistry experiments

For immunofluorescence (IF) experiments, frozen kidney sections (10 μm) and hRPTEpiC slices were fixed with ice-cold acetone for 10 min. The sections were washed with PBS for four times and blocked with blocking buffer (1% BSA and 2% donkey serum diluted in PBS) at room temperature for 30 min. Primary antibodies were diluted in blocking buffer and incubated overnight at 4 °C; the primary antibodies used were anti-flavivirus envelope protein antibody (1: 200, clone D1-4G2-4-15, Millipore, Darmstadt, Germany), rabbit anti-Vimentin (1:200, #5741, Cell Signaling Technology, CST, Danvers, MA, USA) and rabbit anti-cleaved caspase-3 (1:400, #9661, CST). The secondary antibodies used were Alexa Fluor 488 donkey anti-rabbit IgG H&L (1:1 000, ab150073, Abcam) and Alexa Fluor 568 donkey anti-mouse IgG (H+L) (1:1 000, ab175472, Abcam, Cambridge, UK).

For immunohistochemical (IHC) staining, frozen kidney sections (10 μm) were treated as described above. Primary antibodies were diluted in blocking buffer and incubated overnight at 4 °C; the primary antibodies used were rat anti-CD4 (1:200, ab25475, Abcam), rat anti-CD8 (1:400, ab22378, Abcam), rat anti-neutrophile (1:200, ab2557, Abcam) and rat anti-monocyte chemoattractant protein-1 (MCP-1; 1:200, ab8101, Abcam). After rinsing with PBS, the primary antibodies were subsequently detected with a rabbit anti-rat IgG H&L (HRP) secondary antibody. Sections were lightly counterstained with Mayer’s hematoxylin. Images were taken using a Tissue FAXS 200 flow-type tissue quantitative analyzer (TissueGnostics GmbH, Vienna, Austria).

Statistics

For RT-qPCR and TCID₅₀ quantification, five biological replicates were performed, and the results are displayed as the mean±sd. Unless otherwise noted, statistical analyses were performed using two-tailed Student’s t-test.

Whole-transcriptome (RNA-Seq) analysis

Mock and infected hRPTEpiCs (three replicates in each group) were harvested at 48 h post infection. Whole-transcriptome sequencing was performed by Genergy Biotechnology Inc. (Shanghai, China). RNA-Seq short reads were aligned to the human genome (GRCh38) using GSNAP⁴⁶ with a maximum of two mismatches. On average, ~23.7 million reads across all samples were aligned to the reference, which accounted for 95.9% of total reads. Gene expression was determined as number of short reads that fully/partially aligned to the annotated gene model using HTseq.⁴⁷ Expressed genes were defined as those genes having more than 10 total mapped reads in all samples with at least two of three replicates having more than two reads. In total, 17 699 genes met the criteria and were defined as expressed in both mock and ZIKV-infected libraries. Differential expressed genes (DEGs) were identified using the R Bioconductor package ‘edgeR’.⁴⁸ Trimmed mean of M-values normalization was applied to account for differences in library size among samples. Read counts were fitted into generalized linear model with negative binomial distribution and the effect of viral infection was evaluated by a likelihood ratio test. P-values were adjusted for multiple tests using false discovery rate.⁴⁹ Significant DEGs were identified with a FDR ≤0.05 and a log₂(fold change) ≥1. Gene ontology enrichment analysis was performed using GOrilla⁵⁰ by comparing the up/downregulated DEGs to a list of all expressed genes. Significant GO terms with FDR ≤0.05 were reported.

Go to:

RESULTS

Viral production and ZIKV infection in a mouse model

Although ZIKV can be detected in patients’ urine sample, no renal abnormalities or biopsies in humans have been reported to date. To identify the host cells of ZIKV in kidneys, we chose to use type I and II IFN receptor-deficient mice (AG6) to model the in vivo renal infection of ZIKV. SZ01, an Asian-lineage ZIKV,⁴⁵ was produced in C6/36 cells in high quantity and PFU. Buffered saline or 1 × 10⁵ PFU of SZ01 ZIKV was inoculated into AG6 mice via i.p. injection. Kidneys from both the mock and ZIKV groups were harvested and sectioned at day 7 post infection. No bleeding or anatomic abnormalities of kidneys were observed.

ISH, IHC and IF studies of infected kidneys

The nephron is the basic unit of the kidney and is composed of a renal corpuscle and a tubule (Figure 1A). The renal corpuscle is the initial filtering component of the kidney, consisting of glomerulus and Bowman’s capsule. According to the relative positions, renal tubules are categorized into the proximal tubule, the loop of Henle and the distal tubule. The proximal tubule has a brush border, or so-called microvilli, and thus has a larger lumen than the PCT lumen. Most of the reabsorption and secretion occurs in proximal tubule.

 

Figure 1

Renal infection of Zika virus (ZIKV) in AG6 mice. (A) Components of the nephron. (B) Detection of ZIKV RNA in mouse kidneys. Left panel, in situ hybridization (ISH) using a ZIKV-specific probe. Blue and red arrows indicate renal glomerular and tubular cells, respectively. Hybridization signals of ZIKV RNA were developed by a chromogenic reaction and are observed as black particles. Right panel, hematoxylin and eosin (H&E) staining. Scale bar, 50 μm. (C) Immunofluorescence staining of caspase-3 in ZIKV-infected renal tissues of AG6 mice. Left panel, representative images. Scale bar, 200 μm. Right panel, statistical analyses of caspase-3-positive cells. **Significant difference between mock and infected cells (P<0.01; Student’s t-test).

In our studies, the ZIKV RNA in infected mouse kidneys were detected by ISH using a ZIKV-specific probe. We found that ZIKV RNA was detected in both the renal glomerulus and tubule (Figure 1B, left panel). It appeared that ZIKV signals clustered in the renal glomerulus, whereas the renal tubule displayed a lesser degree of ZIKV signals. This result suggests that ZIKV may penetrate the glomeruli and then into tubules during renal filtration. H&E staining revealed ZIKV-induced swelling of kidneys (Figure 1B, right panel). IF staining of infected renal cells using an anti-caspase-3 antibody indicated ZIKV-induced apoptosis (Figure 1C). IHC staining indicated that infiltrating CD8⁺ cells were significantly increased in ZIKV-infected kidneys, whereas no difference in infiltrating neutrophils or MCP-1 expression was observed between the control and ZIKV-infected kidneys (Supplementary Figure S1). Taken together, our results suggest that mouse renal tubular cells are susceptible in vivo host cells of ZIKV.

ZIKV infection in the immortalized hRPTEpiC line HK2

We next sought to assess ZIKV infection in human renal tubular cells. We first evaluated ZIKV infection in an immortalized hRPTEpiC line, HK2, which was derived from a normal adult human kidney.⁵¹ To avoid bias of a single virus strain, we used both SZ01, an Asian-lineage ZIKV, and MR766, an African-lineage ZIKV (ATCC VR-1838), for our study. IF staining using an anti-ZIKV envelop protein (ZIKVE) antibody showed that both SZ01 and MR766 could efficiently infect HK2 cells at 24 and 48 h post inoculation (Figure 2A). The use of MR766 ZIKV appeared to result in higher infection rate than SZ01 ZIKV. In addition, ZIKV infection induced notable cytopathogenic effects (CPEs) in HK2 cells at 48 h post infection (Figure 2B). To determine the production of infectious viral particles, we analyzed the TCID₅₀ of SZ01 and MR766 ZIKV from infected HK2 cells over a period of 96 h. The results indicated that HK2 could produce infectious SZ01 and MR766 particles in a time-dependent manner. In addition, with the same initial viral dose, HK2 produced more infectious MR766 particles than SZ01 at various time points.

Figure 2

ZIKV infection in the HK2 cell line. (A) Immunofluorescence staining of HK2 cells infected with SZ01 or MR766 ZIKV using an anti-ZIKV envelope protein antibody. Upper panel, representative images. Scale bar, 20 μm. Lower panel, quantification of infected cells. (B) Light microscopy images showing SZ01 ZIKV-induced CPEs. (C) Production of infectious SZ01 and MR766 ZIKV in HK2 cells. Cytopathogenic effects, CPEs.

ZIKV infection in primary hRPTEpiCs

Next, we explored ZIKV infection in primary hRPTEpiCs. ZIKV-infected primary hRPTEpiCs were stained with a ZIKV antibody and an antibody against vimentin, a marker of hRPTEpiCs. We found that hRPTEpiCs could be efficiently infected by ZIKV (Figure 3A), which is in line with the results observed in the HK2 cell line. We next assessed viral replication in hRPTEpiCs by quantifying ZIKV RNA in the supernatant of medium using RT-qPCR. ZIKV was found to persist in hRPTEpiCs for more than 30 days after infection (Figure 3B), which is consistent with the prolonged persistency of ZIKV in patients’ urine samples. Moreover, the MR766 ZIKV infection induced CPEs in hRPTEpiCs over a course of 12 days, as evidenced by the increased cell apoptosis (Figure 3C). By contrast, the mock treatment led to little or no cell apoptosis. In addition, Asian-lineage ZIKV SZ01 induced caspase-3-mediated apoptosis in hRPTEpiCs (Supplementary Figure S2). It was noted that SZ01 ZIKV only induced a modest degree of apoptosis after an extended period of elongated infection, whereas MR766 ZIKV induced a notably higher level of apoptosis at all examined time points. Similar to the results obtained with the HK2 cell line, primary hRPTEpiCs could produce both SZ01 and MR766 infectious viral particles over a course of 17 days. It was also noted that primary hRPTEpiCs produced more infectious MR766 than SZ01 at various time points (Figure 3D). Collectively, these results suggest that hRPTEpiCs are a permissive cell host of ZIKV and that MR766 has higher fitness in hRPTEpiCs than SZ01.

Figure 3

ZIKV infection in primary hRPTEpiCs. (A) Immunofluorescence staining of primary hRPTEpiCs infected with SZ01 ZIKV at a MOI of 2.5 for 96 h. Scale bar, 30 μm. (B) Long-term persistence of SZ01 ZIKV in primary hRPTEpiCs. Viral copies in the supernatant of medium were determined by RT-qPCR using ZIKV standard samples. Values are shown as the mean±standard deviation (sd; n=5). (C) MR766 ZIKV-induced apoptosis in primary hRPTEpiCs. Upper panel, representative immunofluorescence images. Scale bar, 200 μm. Lower panel, quantification of ZIKV-positive and caspase-3-positive cells. ***Significant difference between mock and infected cells (P<0.001; Student’s t-test). (D) Production of infectious SZ01 and MR766 ZIKV particles in primary hRPTEpiCs.

RNA-Seq analysis of ZIKV-infected hRPTEpiCs

To understand the cellular response of hRPTEpiCs to ZIKV infection, we performed RNA-Seq analysis on mock and ZIKV-infected primary hRPTEpiCs at 48 h post infection. Analysis of the whole transcriptome revealed a large set of dysregulated genes in ZIKV-infected hRPTEpiCs, including 726 upregulated genes and 624 downregulated genes (Figure 4A). The top 10 up- and downregulated genes exhibited up to a 40- or 20-fold change of mRNA expression, respectively (Supplementary Tables S1 and S2). Gene ontology (GO) analysis uncovered a variety of perturbed biological processes, including upregulated immune response and cytokine-mediated signaling pathways (Figure 4B). The downregulated biological processes are composed primarily of developmental GO terms. Analyses of differentially expressed genes (DEGs) showed that ZIKV infection stimulated the upregulation of a variety of type I IFN signaling, anti-viral responses and inflammatory genes (Supplementary Figure S3).

 

Figure 4

RNA-Seq analysis of mock and ZIKV-infected primary hRPTEpiCs at 48 h post infection. (A) Volcano plot of genes differentially expressed in mock and ZIKV-infected hRPTEpiCs. (B) Genes with differential expression between mock and infected hRPTEpiCs were subjected to gene ontology (GO) analyses. The top 10 most significantly upregulated and downregulated biological processes are shown in red and blue, respectively.

Go to:

DISCUSSION

First discovered in 1947 in the Zika forest of Uganda, ZIKV initially gained little notice as a threat to public health due to the benign symptoms of infection. However, the ZIKV outbreak in Brazil in 2015 was found to be associated with a notable rise of neonate microcephaly, a medical condition characterized by developmental abnormalities of fetal brains.² Clinical and laboratory studies have confirmed the causal link between ZIKV infection and fetal microcephaly.⁵² A recent study has suggested that the primary targets of ZIKV in the central nervous system are astrocytes,⁵³ although neural progenitor cells⁸ and neural stem cells¹¹ can also be infected. However, the transmission pathway of ZIKV within the central nervous system, as well as the physiological relevance of the prolonged persistency of ZIKV in human fetal neural progenitors,⁹ remains unclear. In addition, ZIKV can damage the peripheral nervous system, leading to a progressive neurological disease known as GBS that causes severe muscle weakness in adults.¹³ The mechanism of ZIKV-induced GBS is largely unknown. In addition to the nervous system, recent studies have highlighted the presence of a viral reservoir of ZIKV in the reproductive system.^24, 27

Despite progress, a long-standing unexplored question is why ZIKV persists longer and has a higher load in urine than in serum. It is most likely that an independent viral reservoir exists in the urinary system. A recent study has shown that ZIKV can infect glomerular podocytes and renal glomerular endothelial and mesangial cells.⁴³ In this study, we showed that inoculation with ZIKV in mice deficient of type I and II IFN receptors led to viral replication in both renal glomeruli and tubules. In addition, renal glomeruli presented a higher density of ZIKV signals than renal tubules. This result suggested a possible mechanism of the penetration of plasma ZIKV to renal tubular cells through filtration of renal glomeruli. We also found that ZIKV exhibited long-term persistence in primary hRPTEpiCs. Moreover, two strains of ZIKV, SZ01 and MR766, could both replicate in immortalized and primary hRPTEpiCs as infectious viral particles. These results together suggest that hRPTEpiCs in kidneys may serve as a viral reservoir for persistent viral replication. In addition, both in vitro and in vivo ZIKV infections induced notable caspase-3-mediated cell apoptosis. These observations suggest caution for ZIKV-associated renal abnormalities. It would therefore be very interesting to investigate ZIKV-induced renal pathogenesis in humans in future studies.

We also noted that different African and Asian ZIKV strains exhibited notably different infectious profiles. Compared with Asian-lineage ZIKV SZ01, African ZIKV MR766 produced a higher level of infectious viral particles in both immortalized and primary hRPTEpiCs (Figures 2C and ​and3D).3D). In addition, MR766 induced higher degree of cell apoptosis than SZ01 (Figure 3C). Previous studies have shown that African ZIKV could lead to a higher infection rate and virus-induced cell death in vitro⁵⁴ and in vivo⁵⁵ compared with those treated with Asian ZIKV. Our current findings are in line with previous results and indicate that African- and Asian-lineage ZIKV have strain-specific pathogenesis profiles that should be taken into account in future studies.

We noted that our RNA-Seq results were different from those obtained from ZIKV-infected human neural progenitor cells, whereby the majority of the most significantly dysregulated GO terms were cell death or cell cycle-related biological processes.⁸ This discrepancy could be due to the cell type-specific response to ZIKV infection. In addition, RNA-Seq analyses suggested that ZIKV did not inhibit the anti-viral response of hRPTEpiCs as was seen in other cell types.⁵⁶ This result is likely attributable to the different innate immune response mechanisms in various cell types.

Another important implication of our study is that it provides experimental support for the CDC-approved diagnostic method of ZIKV. The lack of laboratory studies regarding the renal infection of ZIKV, especially the unclear origin of prolonged ZIKV persistency in urine, raised concerns of whether urine samples should be used for the accurate diagnosis of ZIKV.^40, 42 Our study showed that ZIKV can infect RPTEpiCs both in vitro and in vivo, thus rationalizing the use of urine samples for non-invasive ZIKV diagnosis.

Go to:

Acknowledgments

This work was supported by the Natural Science Foundation of China (81430030), ShanghaiTech University (Startup Fund to BJ and JL), National Basic Research Program of China (2014CB542502), Shanghai ShenKang Hospital Development Center (SHDC12014104), Shanghai Public Health Clinical Center (2016-28 ksf0730), and the Technology Service Platform for Detecting High level Biological Safety Pathogenic Microorganism (15DZ2290200, Shanghai Science and Technology Commission).

Go to:

Footnotes

Supplementary Information for this article can be found on the Emerging Microbes & Infections website (http://www.nature.com/emi)

Go to:

Supplementary Material

Supplementary Figure S1

Click here for additional data file.^{(473K, pdf)}

Supplementary Figure S2

Click here for additional data file.^{(250K, pdf)}

Supplementary Figure S3

Click here for additional data file.^{(330K, pdf)}

Supplementary Table S1

Click here for additional data file.^{(157K, pdf)}

Supplementary Table S2

Click here for additional data file.^{(157K, pdf)}

Go to:

References

1. Duffy MR, Chen TH, Hancock WT et al. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med 2009; 360: 2536–2543. [PubMed] [Google Scholar]
2. Brasil P, Pereira JP Jr, Moreira ME et al. Zika virus infection in pregnant women in Rio de Janeiro. N Engl J Med 2016; 375: 2321–2334. [PMC free article] [PubMed] [Google Scholar]
3. Wu KY, Zuo GL, Li XF et al. Vertical transmission of Zika virus targeting the radial glial cells affects cortex development of offspring mice. Cell Res 2016; 26: 645–654. [PMC free article] [PubMed] [Google Scholar]
4. Miner JJ, Cao B, Govero J et al. Zika virus infection during pregnancy in mice causes placental damage and fetal demise. Cell 2016; 165: 1081–1091. [PMC free article] [PubMed] [Google Scholar]
5. Cugola FR, Fernandes IR, Russo FB et al. The Brazilian Zika virus strain causes birth defects in experimental models. Nature 2016; 534: 267–271. [PMC free article] [PubMed] [Google Scholar]
6. Li C, Xu D, Ye Q et al. Zika virus disrupts neural progenitor development and leads to microcephaly in mice. Cell Stem Cell 2016; 19: 120–126. [PubMed] [Google Scholar]
7. Adams Waldorf KM, Stencel-Baerenwald JE, Kapur RP et al. Fetal brain lesions after subcutaneous inoculation of Zika virus in a pregnant nonhuman primate. Nat Med 2016; 22: 1256–1259. [PMC free article] [PubMed] [Google Scholar]
8. Tang H, Hammack C, Ogden SC et al. Zika virus infects human cortical neural progenitors and attenuates their growth. Cell Stem Cell 2016; 18: 587–590. [PMC free article] [PubMed] [Google Scholar]
9. Hanners NW, Eitson JL, Usui N et al. Western Zika virus in human fetal neural progenitors persists long term with partial cytopathic and limited immunogenic effects. Cell Rep 2016; 15: 2315–2322. [PMC free article] [PubMed] [Google Scholar]
10. Dang J, Tiwari SK, Lichinchi G et al. Zika virus depletes neural progenitors in human cerebral organoids through activation of the innate immune receptor TLR3. Cell Stem Cell 2016; 19: 258–265. [PMC free article] [PubMed] [Google Scholar]
11. Liang Q, Luo Z, Zeng J et al. Zika virus NS4A and NS4B proteins deregulate Akt-mTOR signaling in human fetal neural stem cells to inhibit neurogenesis and induce autophagy. Cell Stem Cell 2016; 19: 663–671. [PMC free article] [PubMed] [Google Scholar]
12. Cao-Lormeau VM, Blake A, Mons S et al. Guillain-Barre syndrome outbreak associated with Zika virus infection in French Polynesia: a case-control study. Lancet 2016; 387: 1531–1539. [PMC free article] [PubMed] [Google Scholar]
13. Parra B, Lizarazo J, Jimenez-Arango JA et al. Guillain-Barre syndrome associated with Zika virus infection in Colombia. N Engl J Med 2016; 375: 1513–1523. [PubMed] [Google Scholar]
14. D’Ortenzio E, Matheron S, Yazdanpanah Y et al. Evidence of sexual transmission of Zika virus. N Engl J Med 2016; 374: 2195–2198. [PubMed] [Google Scholar]
15. Turmel JM, Abgueguen P, Hubert B et al. Late sexual transmission of Zika virus related to persistence in the semen. Lancet 2016; 387: 2501. [PubMed] [Google Scholar]
16. Brooks RB, Carlos MP, Myers RA et al. Likely sexual transmission of Zika virus from a man with no symptoms of infection - Maryland, 2016. MMWR Morb Mortal Wkly Rep 2016; 65: 915–916. [PubMed] [Google Scholar]
17. Deckard DT, Chung WM, Brooks JT et al. Male-to-male sexual transmission of Zika virus—Texas, January 2016. MMWR Morb Mortal Wkly Rep 2016; 65: 372–374. [PubMed] [Google Scholar]
18. Davidson A, Slavinski S, Komoto K et al. Suspected female-to-male sexual transmission of Zika virus - New York City, 2016. MMWR Morb Mortal Wkly Rep 2016; 65: 716–717. [PubMed] [Google Scholar]
19. Mansuy JM, Dutertre M, Mengelle C et al. Zika virus: high infectious viral load in semen, a new sexually transmitted pathogen? Lancet Infect Dis 2016; 16: 405. [PubMed] [Google Scholar]
20. Atkinson B, Hearn P, Afrough B et al. Detection of Zika virus in semen. Emerg Infect Dis 2016; 22: 940. [PMC free article] [PubMed] [Google Scholar]
21. Visseaux B, Mortier E, Houhou-Fidouh N et al. Zika virus in the female genital tract. Lancet Infect Dis 2016; 16: 1220. [PubMed] [Google Scholar]
22. Prisant N, Bujan L, Benichou H et al. Zika virus in the female genital tract. Lancet Infect Dis 2016; 16: 1000–1001. [PubMed] [Google Scholar]
23. Nicastri E, Castilletti C, Balestra P et al. Zika virus infection in the central nervous system and female genital tract. Emerg Infect Dis 2016; 22: 2228–2230. [PMC free article] [PubMed] [Google Scholar]
24. Yockey LJ, Varela L, Rakib T et al. Vaginal exposure to Zika virus during pregnancy leads to fetal brain infection. Cell 2016; 166: 1247–1256. [PMC free article] [PubMed] [Google Scholar]
25. Tang WW, Young MP, Mamidi A et al. A mouse model of Zika virus sexual transmission and vaginal viral replication. Cell Rep 2016; 17: 3091–3098. [PMC free article] [PubMed] [Google Scholar]
26. Govero J, Esakky P, Scheaffer SM et al. Zika virus infection damages the testes in mice. Nature 2016; 540: 438–442. [PMC free article] [PubMed] [Google Scholar]
27. Ma W, Li S, Ma S et al. Zika virus causes testis damage and leads to male infertility in mice. Cell 2016; 167: 1511–1524. [PubMed] [Google Scholar]
28. Kutsuna S, Kato Y, Takasaki T et al. Two cases of Zika fever imported from French Polynesia to Japan, December 2013 to January 2014. Euro Surveill 2014; 19: 20683. [PubMed] [Google Scholar]
29. Besnard M, Lastere S, Teissier A et al. Evidence of perinatal transmission of Zika virus, French Polynesia, December 2013 and February 2014. Euro Surveill 2014; 19: 20751. [PubMed] [Google Scholar]
30. Gourinat AC, O’Connor O, Calvez E et al. Detection of Zika virus in urine. Emerg Infect Dis 2015; 21: 84–86. [PMC free article] [PubMed] [Google Scholar]
31. Korhonen EM, Huhtamo E, Smura T et al. Zika virus infection in a traveller returning from the Maldives, June 2015. Euro Surveill 2016; 21: 30107. [PubMed] [Google Scholar]
32. Zhang FC, Li XF, Deng YQ et al. Excretion of infectious Zika virus in urine. Lancet Infect Dis 2016; 16: 641–642. [PubMed] [Google Scholar]
33. Calvet G, Aguiar RS, Melo AS et al. Detection and sequencing of Zika virus from amniotic fluid of fetuses with microcephaly in Brazil: a case study. Lancet Infect Dis 2016; 16: 653–660. [PubMed] [Google Scholar]
34. Oliveira Melo AS, Malinger G, Ximenes R et al. Zika virus intrauterine infection causes fetal brain abnormality and microcephaly: tip of the iceberg? Ultrasound Obstet Gynecol 2016; 47: 6–7. [PubMed] [Google Scholar]
35. Lanciotti RS, Kosoy OL, Laven JJ et al. Genetic and serologic properties of Zika virus associated with an epidemic, Yap State, Micronesia, 2007. Emerg Infect Dis 2008; 14: 1232–1239. [PMC free article] [PubMed] [Google Scholar]
36. Bingham AM, Cone M, Mock V et al. Comparison of test results for Zika virus RNA in urine, serum, and saliva specimens from persons with travel-associated Zika virus disease - Florida, 2016. MMWR Morb Mortal Wkly Rep 2016; 65: 475–478. [PubMed] [Google Scholar]
37. Campos Rde M, Cirne-Santos C, Meira GL et al. Prolonged detection of Zika virus RNA in urine samples during the ongoing Zika virus epidemic in Brazil. J Clin Virol 2016; 77: 69–70. [PubMed] [Google Scholar]
38. Roze B, Najioullah F, Ferge JL et al. Zika virus detection in urine from patients with Guillain-Barre syndrome on Martinique, January 2016. Euro Surveill 2016; 21: 30154. [PubMed] [Google Scholar]
39. Interim guidance for Zika virus testing of urine - United States, 2016. MMWR Morb Mortal Wkly Rep 2016; 65: 474. [PubMed] [Google Scholar]
40. Wiwanitkit V. Urine-based molecular diagnosis of Zika virus. Int Urol Nephrol 2016; 48: 2023. [PubMed] [Google Scholar]
41. Lamb LE, Bartolone SN, Kutluay SB et al. Advantage of urine based molecular diagnosis of Zika virus. Int Urol Nephrol 2016; 48: 1961–1966. [PubMed] [Google Scholar]
42. Lamb LE, Bartolone SN, Kutluay SB et al. Reply to the letter: Urine based molecular diagnosis of Zika virus by Viroj Wiwanitkit. Int Urol Nephrol 2016; 48: 2025. [PubMed] [Google Scholar]
43. Alcendor DJ. Zika virus infection of the human glomerular cells: implications for viral reservoirs and renal pathogenesis. J Infect Dis 2017. [PMC free article] [PubMed]
44. Hamel R, Dejarnac O, Wichit S et al. Biology of zika virus infection in human skin cells. J Virol 2015; 89: 8880–8896. [PMC free article] [PubMed] [Google Scholar]
45. Deng YQ, Zhao H, Li XF et al. Isolation, identification and genomic characterization of the Asian lineage Zika virus imported to China. Sci China Life Sci 2016; 59: 428–430. [PubMed] [Google Scholar]
46. Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics 2010; 26: 873–881. [PMC free article] [PubMed] [Google Scholar]
47. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 2015; 31: 166–169. [PMC free article] [PubMed] [Google Scholar]
48. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010; 26: 139–140. [PMC free article] [PubMed] [Google Scholar]
49. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B 1995; 57: 289–300. [Google Scholar]
50. Eden E, Navon R, Steinfeld I et al. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 2009; 10: 48. [PMC free article] [PubMed] [Google Scholar]
51. Ryan MJ, Johnson G, Kirk J et al. HK-2: an immortalized proximal tubule epithelial cell line from normal adult human kidney. Kidney Int 1994; 45: 48–57. [PubMed] [Google Scholar]
52. Ming GL, Tang H, Song H. Advances in Zika virus research: stem cell models, challenges, and opportunities. Cell Stem Cell 2016; 19: 690–702. [PMC free article] [PubMed] [Google Scholar]
53. van den Pol AN, Mao G, Yang Y et al. Zika virus targeting in the developing brain. J Neurosci 2017; 37: 2161–2175. [PMC free article] [PubMed] [Google Scholar]
54. Simonin Y, Loustalot F, Desmetz C et al. Zika virus strains potentially display different infectious profiles in human neural cells. EBioMedicine 2016; 12: 161–169. [PMC free article] [PubMed] [Google Scholar]
55. Tripathi S, Balasubramaniam VR, Brown JA et al. A novel Zika virus mouse model reveals strain specific differences in virus pathogenesis and host inflammatory immune responses. PLoS Pathog 2017; 13: e1006258. [PMC free article] [PubMed] [Google Scholar]
56. Grant A, Ponia SS, Tripathi S et al. Zika virus targets human STAT2 to inhibit type I interferon signaling. Cell Host Microbe 2016; 19: 882–890. [PMC free article] [PubMed] [Google Scholar]

Articles from Emerging Microbes & Infections are provided here courtesy of Taylor & Francis