Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2017 Jun 15;546(7658):406-410.
doi: 10.1038/nature22401. Epub 2017 May 24.

Establishment and cryptic transmission of Zika virus in Brazil and the Americas

N R Faria  1   2 J Quick  3 I M Claro  4 J Thézé  1 J G de Jesus  5 M Giovanetti  5   6 M U G Kraemer  1   7   8 S C Hill  1 A Black  9   10 A C da Costa  4 L C Franco  2 S P Silva  2 C-H Wu  11 J Raghwani  1 S Cauchemez  12   13 L du Plessis  1 M P Verotti  14 W K de Oliveira  15   16 E H Carmo  17 G E Coelho  18   19 A C F S Santelli  18   20 L C Vinhal  18 C M Henriques  17 J T Simpson  21 M Loose  22 K G Andersen  23 N D Grubaugh  23 S Somasekar  24 C Y Chiu  24 J E Muñoz-Medina  25 C R Gonzalez-Bonilla  25 C F Arias  26 L L Lewis-Ximenez  27 S A Baylis  28 A O Chieppe  29 S F Aguiar  29 C A Fernandes  29 P S Lemos  2 B L S Nascimento  2 H A O Monteiro  2 I C Siqueira  5 M G de Queiroz  30 T R de Souza  30   31 J F Bezerra  30   32 M R Lemos  33 G F Pereira  33 D Loudal  33 L C Moura  33 R Dhalia  34 R F França  34 T Magalhães  34   35 E T Marques Jr  34 T Jaenisch  36 G L Wallau  34 M C de Lima  37 V Nascimento  37 E M de Cerqueira  37 M M de Lima  38 D L Mascarenhas  39 J P Moura Neto  40 A S Levin  4 T R Tozetto-Mendoza  4 S N Fonseca  41 M C Mendes-Correa  4 F P Milagres  42 A Segurado  4 E C Holmes  43 A Rambaut  44   45 T Bedford  9 M R T Nunes  2   46 E C Sabino  4 L C J Alcantara  5 N J Loman  3 O G Pybus  1   47
Affiliations

Establishment and cryptic transmission of Zika virus in Brazil and the Americas

N R Faria et al. Nature. .

Abstract

Transmission of Zika virus (ZIKV) in the Americas was first confirmed in May 2015 in northeast Brazil. Brazil has had the highest number of reported ZIKV cases worldwide (more than 200,000 by 24 December 2016) and the most cases associated with microcephaly and other birth defects (2,366 confirmed by 31 December 2016). Since the initial detection of ZIKV in Brazil, more than 45 countries in the Americas have reported local ZIKV transmission, with 24 of these reporting severe ZIKV-associated disease. However, the origin and epidemic history of ZIKV in Brazil and the Americas remain poorly understood, despite the value of this information for interpreting observed trends in reported microcephaly. Here we address this issue by generating 54 complete or partial ZIKV genomes, mostly from Brazil, and reporting data generated by a mobile genomics laboratory that travelled across northeast Brazil in 2016. One sequence represents the earliest confirmed ZIKV infection in Brazil. Analyses of viral genomes with ecological and epidemiological data yield an estimate that ZIKV was present in northeast Brazil by February 2014 and is likely to have disseminated from there, nationally and internationally, before the first detection of ZIKV in the Americas. Estimated dates for the international spread of ZIKV from Brazil indicate the duration of pre-detection cryptic transmission in recipient regions. The role of northeast Brazil in the establishment of ZIKV in the Americas is further supported by geographic analysis of ZIKV transmission potential and by estimates of the basic reproduction number of the virus.

PubMed Disclaimer

Figures

Extended Data Figure 1
Extended Data Figure 1
a. The distribution of CT-values for the RT-qPCR+ samples tested during the ZiBRA journey in Brazil (n=181 samples; median CT = 35.96). b. shows the distribution of the temporal lag between the date of onset of clinical symptoms and the date of sample collection of RT-qPCR+ samples (median lag = 2 days). Red dashed lines represent the median of the distributions. (c) Validation of sequencing approaches. A phylogeny of the ZIKV Asian genotype estimated using PhyML is shown. The expanded clade highlighted in blue contains the WHO reference ZIKV sequence (accession number KX369547), which was generated using Illumina MiSeq. Sequences generated using MinION chemistries R9.4 2D, R9.4 1D, R9 1D, R9 2D and R7.3 2D contain no nucleotide differences and hence were also placed in this clade. Scale bars represent expected nucleotide substitutions per site (s/s). Am-ZIKV=American Zika virus lineage.
Extended Data Figure 2
Extended Data Figure 2
Temporal signal of the ZIKV Asian genotype. The correlation between sampling dates and genetic distances from the tips to the root of a maximum likelihood (ML) tree, estimated using PhyML, was explored using TempEst. a. Estimates for the dataset used in the phylogenetic analysis presented in Fig. 3c, and b. estimates for the same dataset with the addition of the P6-740 strain sampled in 1966 (accession number HQ234499).
Extended Data Figure 3
Extended Data Figure 3
A non-clock maximum likelihood phylogeny of our ZIKV data set. Bootstrap branch support values are shown at each node. The phylogeny was estimated using PhyML. Sequences generated in this study are highlighted in red. Scale bar represents expected nucleotide substitutions per site.
Extended Data Figure 4
Extended Data Figure 4
Ancestral node location posterior probabilities (ANLPP), for nodes A, B and C, estimated using the complete dataset (top row) and ten replicate subsampled data sets (other rows). See Methods for details. ANLPPs were calculated using two approaches: DTA=discrete trait analysis method (left side columns) and BASTA=Bayesian structured coalescent approximation method (right side columns). For each method, we employed an asymmetric model of location exchange to estimate ancestral node locations and to infer patterns of virus spread among regions.
Extended Data Figure 5
Extended Data Figure 5
Epidemic growth rates estimated from weekly ZIKV notified cases in Brazil. Time series show the number of ZIKV notified cases in each region of Brazil. Periods from which exponential growth were estimated are highlighted in grey.
Extended Data Figure 6
Extended Data Figure 6
Seasonal suitability for ZIKV transmission in the Americas. These maps were estimated by collating data on Aedes mosquitoes, temperature, relative humidity and precipitation, and are the basis of the trends in suitability for different regions shown in main text Figs. 1 and 4. For method details, see ,.
Extended Data Figure 7
Extended Data Figure 7
Partial autocorrelation functions for the linear model associating climatic suitability and ZIKV notified cases in each geographic region in Brazil. The residuals for the North, Northeast, Centre-West and Southeast regions show no autocorrelation, while a small amount of autocorrelation cannot be excluded for the South region.
Fig. 1
Fig. 1. Geographic and temporal distribution of ZIKV in Brazil
a. Sampling location of genome sequences from Brazil and the Americas. Federal states in Brazil are coloured according to 5 geographic regions (lower inset). A red line surrounds the states surveyed by the ZiBRA mobile lab in 2016. State codes are PA=Pará, MA=Maranhão, CE=Ceará, TO=Tocantins, RN=Rio Grande do Norte, PB=Paraíba, PE=Pernambuco, AL=Alagoas, BA=Bahia, RJ=Rio de Janeiro, SP=São Paulo. Underlined states represent those from which sequences in this study were generated (upper inset). Publicly available sequences were also collated from non-underlined states. b. Confirmed and notified ZIKV cases in NE Brazil. Upper panel shows the temporal distribution of RT-qPCR+ cases detected during ZiBRA fieldwork. Only samples with known collection dates are included (n=138 out of 181 confirmed cases). Lower panel shows notified ZIKV cases in NE Brazil between 01 Jan 2015 and 19 Nov 2016 (n=122,779). The dashed line represents the average climatic vector suitability score for NE Brazil (Methods). The vertical arrow indicates date of ZIKV confirmation in NE Brazil/Americas. c. Notified ZIKV cases in the Centre-West, Southeast, North, and South regions of Brazil (clockwise from top left). The dashed lines represent the average climatic vector suitability score for each region.
Fig. 2
Fig. 2. Zika virus genetic diversity and sequencing statistics
a. The percentage of ZIKV genome sequenced plotted against RT-qPCR Ct-value, for each sample. Each circle represents a sequence recovered from an infected individual in Brazil and is coloured by sampling location. b. Illustration of sequencing coverage across the ZIKV genome for the ZiBRA sequences, including data generated by both mobile and static laboratories. c. Regression of sequence sampling dates against root-to-tip genetic distances in a maximum likelihood phylogeny of the Asian-ZIKV lineage. Extended Data Fig. 2b contains a comparable analysis that also includes P6-740 (the oldest Asian-ZIKV strain collected in 1966). d. Average pairwise genetic diversity of the PreAm-ZIKV strains (grey line) and of the Am-ZIKV lineage (black line), calculated using a sliding window of 300 nucleotides with a step size of 50 nucleotides.
Fig. 3
Fig. 3. Phylogeography of ZIKV in the Americas
Maximum clade credibility phylogeny, estimated from complete and partial Am-ZIKV genomes using a molecular clock phylogeographic approach (Methods). Terminal branches with yellow circles indicate sequences reported in this study. Terminal branches with no circles and reduced opacity are those reported in a companion paper. Thin vertical grey boxes indicate statistical uncertainty of estimated dates of nodes A, B and C (Extended Data Table 3c). Branch colours indicate the most probable ancestral lineage locations. Diamonds at internal nodes are sized in proportion to clade posterior probabilities. For selected nodes, coloured numbers show the posterior probabilities of ancestral locations and numbers in grey are clade posterior probabilities. Asterisks indicate the three available genomes from microcephaly cases. A black arrow indicates the oldest Brazilian ZIKV sequence. The grey arrow and dotted line denotes when ZIKV was first confirmed in the Americas. Nodes A and B are equivalent to the nodes named identically in. Text labels along the bottom of the figure denote clades of sequences from regions outside of NE Brazil. RJ1 to RJ4 are clades from Rio de Janeiro state, TO from Tocantins, and SP1 from São Paulo state. Clades from outside Brazil are denoted CB1 and CB2 (Caribbean), SA1 and SA2 (South America excluding Brazil), and CA1 (Central America). Thin grey horizontal lines along the bottom of the figure denote sequences from Brazil.
Fig. 4
Fig. 4. Establishment of Am-ZIKV in the Americas
The earliest inferred dates of lineage export to non-Brazilian regions, represented by box-and-whisker plots. Each plot corresponds to the earliest movement between a pair of locations with well-supported virus lineage migration. The first exports to South America outside Brazil (SA1 in Fig. 3), to Central America (CA1) and to the Caribbean (CB1) are shown in panels a–c, respectively. Box and whisker plots were generated in ggplot2, with boxes representing the median and interquartile ranges of the estimated date of earliest movement. In each of ac, dashed lines show the estimated climatic vector suitability score for each recipient region, averaged across the countries for which sequence data is available (see Methods). In each of a-c, the bar plots show available notified ZIKV case data (plots adapted from PAHO) for the countries with the earliest confirmed cases (Colombia in panel a, Mexico in b, and Puerto Rico in c). Coloured arrows indicate the earliest confirmation of ZIKV autochthonous cases in each non-Brazilian region. The vertical dashed line represents the date of ZIKV confirmation in the Americas.

Comment in

References

    1. Kindhauser MK, Allen T, Frank V, Santhana RS, Dye C. Zika: the origin and spread of a mosquito-borne virus. Bulletin of the World Health Organization. 2016;94:675–686C. doi: 10.2471/BLT.16.171082. - DOI - PMC - PubMed
    1. Ministério da Saúde. Boletins Epidemiológicos—Secretaria de Vigilância em Saúde. 2017 http://portalsaude.saude.gov.br/index.php/o-ministerio/principal/secreta....
    1. WHO. Situation Report - Zika virus, microcephaly, Guillain-Brarré syndrome. 2017 Jan 18; ( http://apps.who.int/iris/bitstream/10665/253604/1/zikasitrep20Jan17-eng...., 2017)
    1. Faria NR, et al. Zika virus in the Americas: Early epidemiological and genetic findings. Science. 2016;352:345–349. doi: 10.1126/science.aaf5036. - DOI - PMC - PubMed
    1. Alex Perkins T, Siraj AS, Ruktanonchai CW, Kraemer MU, Tatem AJ. Model-based projections of Zika virus infections in childbearing women in the Americas. Nat Microbiol. 2016;1:16126. doi: 10.1038/nmicrobiol.2016.126. - DOI - PubMed

Publication types