Title: Early Transmission Dynamics, Spread, and Genomic Characterization of Sars-Cov-2 in Panama
Title: Early Transmission Dynamics, Spread, and Genomic Characterization of Sars-Cov-2 in Panama
Title: Early Transmission Dynamics, Spread, and Genomic Characterization of Sars-Cov-2 in Panama
Title: Early transmission dynamics, spread, and genomic characterization of SARS-CoV-2 in Panama.
Danilo Franco1§, Claudia Gonzalez1,7§, Leyda E. Abrego1,7§, Jean Paul Carrera1,3, Yamilka Diaz1, Yaset
Caisedo4, Ambar Moreno1, Oris Chavarria1, Jessica Gondola1, Marlene Castillo1, Elimelec Valdespino1, Melissa
Gaitán1, Jose Martínez-Mandiche1, Lizbeth Hayer2, Pablo Gonzalez2, Carmen Lange2, Yadira Molto2, Dalis
Mojica1, Ruben Ramos1, Maria Mastelari2, Lizbeth Cerezo2, Lourdes Moreno2, Christl A. Donnelly5,6, Nuno
Rodrigues Faria3,6, Juan Miguel Pascale 1,7, Sandra Lopez-Verges1,7*, Alexander A Martinez1,7*, on behalf of
Gorgas COVID19 team and Panama COVID19 laboratory network#.
§
Contributed equally to this work as co-first authors.
*Co-corresponding authors.
#
Author list can be found in Supplementary information.
Affiliations
1. Gorgas Memorial Institute of Health Studies, Panama City, Panama
2. Ministry of Health of Panama, Panama City, Panama
3. Department of Zoology, University of Oxford, United Kingdom
4. Fundación Valle del Lili. Centro de Investigaciones Clínicas, Cali, Colombia
5. Department of Statistics, University of Oxford, United Kingdom
6. MRC Centre for Global Infectious Disease Analysis (MRC-GIDA), Department of Infectious Disease
Epidemiology, Imperial College London, London, United Kingdom
7. Universidad de Panama, Panama City, Panama
Corresponding authors
Alexander Augusto Martinez Caballero
Department of Research in Genomics and Proteomics.
E-mail: [email protected].
ORCID: 0000-0002-7816-3040
Sandra Lopez-Verges
Department of Research in Virology and Biotechnology.
E-mail: [email protected].
ORCID: 0000-0002-1106-8479
Address:
Gorgas Memorial Institute of Health Studies
Justo Arosemena Street 35th, Panama, Panama.
Postal Code 0816-02593
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
Research in Context
Summary
Background With more than 50000 accumulated cases, Panama has one of the highest incidences of SARS-
CoV-2 in Central America, despite the fast implementation of disease control strategies. We investigated the
early transmission patterns of the virus and the outcomes of mitigation measures in the country.
Methods We collected information from epidemiological surveillance, including contact tracing, and genetic
data from SARS-CoV-2 whole genomes, of the first five weeks of the outbreak. These data were used to
estimate the exponential growth rate, doubling time and the time-varying effective reproductive number (Rt)
using date of symptom onset in a Bayesian framework. The time of most recent ancestor for the introduced and
circulating lineages was estimated by Bayesian analysis.
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
Findings A total of 4210 subjects were SARS-CoV-2 positive during the period evaluated, of them we
sequenced 313 cases, detecting the circulation of 10 SARS-CoV-2 lineages. Whole genomes analysis identified
the local transmission of one cryptic lineage as early as 2 weeks before it was detected by surveillance systems.
Analysis of transmission dynamics showed that lockdown reduced Rt and increased the doubling time, however,
these measures did not stop the circulation of this lineage in the country.
Interpretation These results demonstrate the value of epidemiological modeling and genome surveillance to
assess mitigation strategies. At the same time, an active search for cryptic transmission clusters is crucial to
interrupt local transmission of SARS-CoV-2 in a region.
Funding
Ministry of Health, Contribution from private donors and Secretaría Nacional de Ciencia y Tecnología.
Introduction
A severe respiratory disease called Coronavirus Disease 19 (COVID-19) emerged in early December 2019, in
Wuhan City, Capital of Hubei province in China 1. The etiologic agent was designated as SARS-CoV-2, a single
strand, positive sense enveloped RNA virus of the genus betacoronavirus, family Coronaviridae2. After several
zoonotic viral spillover events2–4, epidemiological investigations confirmed human-to-human transmission of
SARS-CoV-2 among the Wuhan population5.
By January 30th 2020, COVID-19 was declared a public health emergency of international concern by the World
Health Organization (WHO), and by March 11, COVID-19 was declared a pandemic. Previous analysis of
transmission dynamics suggested that SARS-CoV-2 is highly contagious and transmissible among the human
population, initial calculations estimated that the basic reproductive number (R0) was approximately 2.2 1,2,6.
However, others estimates indicated that the R0 could be higher7 with a significant role of asymptomatic
carriers8. As of July 21, SARS-CoV-2 has spread to 188 countries with around 14 763 911 confirmed cases and
611 322 fatalities9. To characterize SARS-CoV-2 local transmission across the globe10, the few mutations
accumulated to date in the viral genome11–13 have allowed the characterization of different lineages or clades (A
and B)14.
In Latin America, SARS-CoV-2 was first reported on February 25th, in Brazil15. In Panama, the first COVID-19
confirmed case was on March 9th in an inbound traveler returning from Spain. However, on March 7th a fatal
case of severe respiratory disease was identified as a COVID-19 suspected case, a tissue sample was collected
the 8th and laboratory confirmation was performed the next day. Although Panama rapidly implemented disease
control strategies, in Latin America it is among the countries with the highest cumulative incidence rate and
fatalities16. We hypothesized that the high incidence of COVID-19 observed in Panama was associated with
SARS-CoV-2 cryptic transmission and early spread, despite the early laboratory analysis of suspected cases and
the rapid implementation of disease control strategies. In this article, we gather data on epidemiological
surveillance, contact tracing, and genetics from whole genomes during the first five weeks of the outbreak, to
evaluate the dynamics of the introduction and spread of SARS-CoV-2 in Panama.
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
interval with 95 % credible intervals (95% Cl), using the EpiEstim packaged implemented in R18. Details are
provided in Supplementary Information.
SARS-CoV-2 Sequencing
Reverse transcription and multiplex PCR were performed according to ARTIC Network protocol
(https://artic.network/), with modifications. Amplicons generated were pooled and prepared for Illumina
Sequencing with the Nextera XT library preparation protocol. All samples were sequenced with MiSeq V2
reagent kits for 500 cycles. Reads obtained were aligned with bowtie-2 and consensus genome was generated
with Quasitools-Hydra. Details are provided in Supplementary Information.
Genome Analysis
To evaluate virus lineage distribution in Panama, lineages were classified using Pangolin14, and their distribution
was assessed per geographic region of origin of the sample and epidemiological information. Time-trees were
estimated using BEAST v1.10.4 19 and BEAGLE v3 20 to enhance computational speed. Posterior estimates of
most recent common ancestor (tMRCA) of each locally transmitted lineage was obtained as described in
Supplementary Information.
Control measures in Latin America
Reports from the WHO, the Organization for Economic Co-operation and Development (OECD) and official
communications from the ministries of health were used to obtain data on control and mitigation strategies at the
country level to reconstruct a timeline representing the strategies adopted by Latin American countries. Only
when no official information was found, press articles from local newspapers were consulted, as described in
Supplementary Information.
Results
Early epidemic dynamics in Panama
Epidemiological surveillance began on January 23 and, by April 16, a total of 18 559 suspected cases of
COVID-19 had been investigated. Of these, 4 210 were confirmed as SARS-CoV-2 infections at that date. On
March 7th, the National Surveillance system detected a fatal case in which the SARS-CoV-2 infection was
confirmed on March 9th. This case was not associated with previous travel history outside the country. The first
case detected by active surveillance at international airports, came from Spain, had mild disease and was
confirmed on March 9th. Most imported cases were associated with travel from Europe, North America, and one
from China.
Epidemiological investigation confirmed that the date of onset of symptoms of one local case in our study goes
back as early as February 15 (Figure 1A). Most local cases had mild disease (Figure 1B). By April 16, a total of
341 hospitalized cases (77 were hospitalized at diagnosis confirmation) and 116 fatalities (31 were dead at the
diagnosis confirmation) were reported in Panama (Figure 1C and 1D). The highest proportion of confirmed
cases were observed in the age group 20-59 years (Supplementary Figure 1A). A higher proportion of females
were tested while more males were positive (Supplementary Figure 1B). A rapid growth rate of 0·13 cases/per
day (Figure 2A) and a short doubling time was also observed during the early stages of the epidemic, with an
increasing doubling time over the analyzed interval (Figure 2B). We estimated an early epidemic potential of
transmission of SARS-CoV-2 in Panama of R0= 2·22 (95% CI, 2·08 - 2·37).
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
Panama detected the first SARS-CoV-2 infection one month after Brazil15, being the 11th country with
confirmed cases in Latin America (Supplementary Figure 2). School closures were implemented the day after
first case confirmation, social distance measures six days later, and stay-at-home 24-hour curfew 14 days later.
Compared to other countries in the Latin American region, Panama was among those that rapidly implemented
control strategies after its first case (Supplementary Figure 2). The effective reproduction number since the
accumulated 25 cases was Rt= 3·01 (95% Crl: 2·54 - 3·52) and rapidly dropped over the course of 17 days to
Rt= 1·08 (95% Crl: 1·00 -1·17) (Figure 2C, Supplementary Table 1). Until April 16, Panama was the country
from Central America with the highest cumulative incidence of cases and fatality (Supplementary Figure 3).
Discussion
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
Panama has the most confirmed SARS-CoV-2 infections and associated fatalities in Central America, although
control strategies were rapidly implemented at the beginning of the outbreak. The reasons for this pattern are
unclear, but they may include, 1) cryptic circulation, 2) rapid early transmission dynamics and 3) actively search
of cases and their contacts, combined with strong molecular capacity of testing in the country.
Epidemiological investigation showed early local transmission events as February 23th, based on reports of
symptoms onset of the first detected local case on March 9th. The potential source of this infection dated back an
onset of symptoms on February 15th. Taken together the epidemiological evidence, along with genetic analysis
of SARS-CoV-2, confirm cryptic circulation in Panama with an interval window of introduction between late
January and early February.
A high early median transmission potential of SARS-CoV-2- was estimated in R0= 2.22 (2.08-2.37), like those
early estimated reported in China and Europe23,24. In addition, the initial time-variant reproduction number
Rt=3.01 (2.54-3.51), estimated since the incidence reached 25 cases, also suggest a high early transmission
dynamic. However, the Rt rapidly dropped after implementation of control strategies.
The genetic analysis of the diverse introductions of lineages and the dynamics of their transmission support
these epidemiological observations. At the beginning of the outbreak, we detected the introduction of ten
lineages, and from them, five were still co-circulating in the country in mid-April. The numbers of detected
lineages could be underestimated because we did not sequence all positive cases, thus there is a possibility of
not-detecting uncommon lineages. Most of the lineages associated with imported cases (A.1, A.3, B, B.1, B.2.1)
were detected by surveillance of suspected cases and their contacts through active contact tracing. However, this
did not happen with the cases associated to the lineage A.2. According to our phylogenetic dating analysis, this
lineage circulated in the country at least three weeks before it was first detected. This lineage was detected
among the school and the police officers clusters, and most of the cases belonging to these groups had their
onset of symptoms between February 15th and March 15th, including the first detected fatal case. The failure to
detect and control the spread of these clusters gave enough time to A.2 lineage to acquire a synonymous
mutation that was fixed in the sequences analyzed and became the A.2 lineage locally transmitted reported in
this study. This descendent lineage (A.2.Pan) was detected in all regions of the country and was not associated
with any epidemic transmission cluster isolated by the public health systems. Moreover, our observations
suggest that in Panama we had a local transmission of lineage A.2.Pan one week before respiratory disease
surveillance detected it. Once this lineage was detected, active contact tracing and testing activities promptly
confirmed more cases. More extensive sampling should be implemented in the country following well-
established pipelines as suggested previously.25
When analyzing the different lineages, the most common mutations (C8782T and T28144C) observed in the
analyzed sequences have been previously reported in South America, United States, Europe and China26,27,
whereas the observed mutations T9477A and C28657T were reported in South America28. These mutations were
not associated with a specific lineage nor with specific transmission or pathogenic effects. The C12815T
mutation detected only in the A.2.Pan lineage has not been reported to date.
The D614G variant was detected in all lineages circulating in Panama during the beginning of SARS-CoV-2
outbreak. This variant was also reported in Europe, Mexico, Brazil, and Wuhan since early February 2020. It
has been related to more efficient virus transmission,29 but with our data do not support the predominance of this
variant in the circulating lineages.
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
The pattern of cryptic transmission in Panama was also observed in other countries which also experienced early
and rapid dissemination of SARS-CoV-2.27,30 We suspect that the tight classification of cases at the beginning of
the surveillance in January, which focused mainly in symptomatic China travelers, precluded the opportunity to
detect these early importations from Europe and USA. Indeed, in these regions the virus was already circulating
at that time.23,27,30 Our forecast of growth rate and R0 showed that mitigation measures undertaken by MoH
helped to reduce the number of new cases of SARS-CoV-2. Implementation of epidemiological control
measurements such as active contact tracing and isolation, social distancing and region-based quarantine where
active transmission clusters are found, would help to effectively control the spread of SARS-CoV-2 as
suggested by results obtained in other countries.24 To implement and maintain these measures until the outbreak
is under control, it is important to improve the communication strategy to provide the public with the
appropriate advice to better prevent infection and subsequent transmission. Improvements in the surveillance
and tracing, data collection and digitalization systems should increase the capacity of the country to respond to
any new transmission waves.
Acknowledgments
Thanks to Gorgas’ administration, general services for their support and extra-work during the COVID19
response. Thanks to all the Health institutions from MINSA, CSS and private hospitals. AAM, JMP, SLV, LEA,
JMP are members of the Sistema Nacional de Investigación (SNI) del SENACYT, Panamá.
Financial support
Laboratory surveillance was done using the Panama National COVID19 response MoH funds and private
donations, Genome surveillance was partially funded by SENACYT COVID19-43 grant and SNI from AAM,
SLV, JMP, and LEA.
References
1 Li Q, Guan X, Wu P, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-
Infected Pneumonia. N Engl J Med 2020; 382: 1199–207.
2 Zhu N, Zhang D, Wang W, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N
Engl J Med 2020; : 1–7.
3 Andersen KG, Rambaut A, Lipkin WI, Holmes EC, Garry RF. The proximal origin of SARS-CoV-2.
Nat Med 2020; 26: 450–2.
4 Zhou P, Yang X-L, Wang X-G, et al. A pneumonia outbreak associated with a new coronavirus of
probable bat origin. Nature 2020; 579: 270–3.
5 Chan JFW, Yuan S, Kok KH, et al. A familial cluster of pneumonia associated with the 2019 novel
coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet 2020; 395:
514–23.
6 Nishiura H, Linton NM, Akhmetzhanov AR. Serial interval of novel coronavirus (COVID-19)
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
April 2020 and its Implications: Data Analysis. JMIR Public Heal Surveill 2020; 6: e19170.
26 Poterico JA, Mestanza O. Genetic variants and source of introduction of SARS-CoV-2 in South
America. J Med Virol 2020; : 0–1.
27 Deng X, Gu W, Federman S, et al. Genomic surveillance reveals multiple introductions of SARS-CoV-2
into Northern California. Science 2020; 9263: 1–11.
28 Castillo AE, Parra B, Tapia P, et al. Phylogenetic analysis of the first four SARS-CoV-2 cases in Chile.
J Med Virol 2020; : 1–5.
29 Korber B, Fischer WM, Gnanakaran S, et al. Spike mutation pipeline reveals the emergence of a more
transmissible form of SARS-CoV-2. bioRxiv 2020; : 2020.04.29.069054.
30 Bedford T, Greninger AL, Roychoudhury P, et al. Cryptic transmission of SARS-CoV-2 in Washington
State. medRxiv Prepr Serv Heal Sci 2020; : 2020.04.02.20051417.
Figure 1. Epidemic Curve of confirmed SARS-CoV-2 infections. Daily incidence of: A) imported and local
confirmed infections, B) ambulatory cases (asymptomatic, pre-symptomatic, mild outpatient cases), C)
hospitalized cases and D) fatal cases, with symptom onset from February 13th to April 13 th. In all four charts, the
y axis represents the daily incidence and the x axis represents the date of symptom onset for each reported case.
The blue vertical dashed lines represent the first recorded onset of symptoms and vertical red dashed lines the
date of first confirmed case by the surveillance system.
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
Figure 2. Epidemiologic parameters: Number of COVID-19 daily confirmed infections (left axis and bars)
overlaid with estimates of A) Fitted exponential growth in cases daily growth rate, B) Doubling time and C)
Time-varying effective reproduction number R, for a time frame of 45 days (x axis). For C), green shaded areas
show 95% credible intervals around the median estimated R. The threshold value Rt = 1 is indicated by the
horizontal red dashed line. For both B) and C), the vertical blue dashed line represents the schools' closure, the
vertical lines red and green represent night curfews and restrictions of movement respectively, and the purple
dashed line the 24-hour curfew.
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
lineages among regions in the country. C) Distribution of the lineages by type of exposure detected by the
surveillance system.
Figure 4. Mutation profile of the sequences obtained in the study. A) Frequency of single nucleotide
polymorphism in the genome (position according to MN908947) for all analyzed sequences (n=313) in the
medRxiv preprint doi: https://doi.org/10.1101/2020.07.31.20160929.this version posted August 4, 2020. The copyright holder for this preprint
(which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
It is made available under a CC-BY-NC 4.0 International license .
study. Distribution of S protein variants (D614 or G614) B) in different regions of the country, and C) among
lineages.