key: cord-0949827-bu9eawah
authors: Cowley, Lauren A.; Afrad, Mokibul Hassan; Rahman, Sadia Isfat Ara; Mamun, Md Mahfuz Al; Chin, Taylor; Mahmud, Ayesha; Rahman, Mohammed Ziaur; Billah, Mallick Masum; Khan, Manjur Hossain; Sultana, Sharmin; Khondaker, Tilovatul; Baker, Stephen; Banik, Nandita; Alam, Ahmed Nawsher; Mannoor, Kaiissar; Banu, Sayera; Chowdhury, Anir; Flora, Meerjady Sabrina; Thomson, Nicholas R.; Buckee, Caroline O.; Qadri, Firdausi; Shirin, Tahmina
title: Genomics, social media and mobile phone data enable mapping of SARS-CoV-2 lineages to inform health policy in Bangladesh
date: 2021-09-08
journal: Nat Microbiol
DOI: 10.1038/s41564-021-00955-3
sha: bf49ff6ea1821bde5b264d4555e0a3f6e0041121
doc_id: 949827
cord_uid: bu9eawah

Genomics, combined with population mobility data, used to map importation and spatial spread of SARS-CoV-2 in high-income countries has enabled the implementation of local control measures. Here, to track the spread of SARS-CoV-2 lineages in Bangladesh at the national level, we analysed outbreak trajectory and variant emergence using genomics, Facebook ‘Data for Good’ and data from three mobile phone operators. We sequenced the complete genomes of 67 SARS-CoV-2 samples (collected by the IEDCR in Bangladesh between March and July 2020) and combined these data with 324 publicly available Global Initiative on Sharing All Influenza Data (GISAID) SARS-CoV-2 genomes from Bangladesh at that time. We found that most (85%) of the sequenced isolates were Pango lineage B.1.1.25 (58%), B.1.1 (19%) or B.1.36 (8%) in early-mid 2020. Bayesian time-scaled phylogenetic analysis predicted that SARS-CoV-2 first emerged during mid-February in Bangladesh, from abroad, with the first case of coronavirus disease 2019 (COVID-19) reported on 8 March 2020. At the end of March 2020, three discrete lineages expanded and spread clonally across Bangladesh. The shifting pattern of viral diversity in Bangladesh, combined with the mobility data, revealed that the mass migration of people from cities to rural areas at the end of March, followed by frequent travel between Dhaka (the capital of Bangladesh) and the rest of the country, disseminated three dominant viral lineages. Further analysis of an additional 85 genomes (November 2020 to April 2021) found that importation of variant of concern Beta (B.1.351) had occurred and that Beta had become dominant in Dhaka. Our interpretation that population mobility out of Dhaka, and travel from urban hotspots to rural areas, disseminated lineages in Bangladesh in the first wave continues to inform government policies to control national case numbers by limiting within-country travel.

surveillance, particularly outside the capital Dhaka, remained challenging. As a result, it has been unclear how the epidemic has spread in Bangladesh and what this means for the potential spread of SARS-CoV-2 variants, and the best use of interventions, including therapeutics and vaccines. Viral genomics analyses have been used to track the SARS-CoV-2 epidemic, enabling fine-scaled transmission mapping and analysis of how changes in population behaviour impact patterns of transmission 4 . For example, genomics analysis has been used to disentangle the timeline of spread in Europe 5 and to understand community transmission versus international importations in New Zealand 6 . Despite sequencing of >2,000 SARS-CoV-2 genomes from Bangladesh, only a handful of studies report the phylogenetic placement of those strains [7] [8] [9] [10] .These phylogenetic approaches provide a window into the biology of transmission that is independent of the capacity of the health system to detect cases. Similarly, mobile phone data have been used extensively as a way to monitor the population behavioural response to the epidemic in real time, and to understand the human drivers of transmission 11, 12 . These new data streams may be particularly powerful tools for monitoring the pandemic in LMICs, in which RT-PCR testing capacity is often highly constrained.

Worldwide genomic surveillance has been rapidly scaled up in response to the emergence of VOCs in late 2020. Pango lineages 4 B.1.1.7, B.1.351, P.1 and B.1.617.2 have been designated as VOCs on the basis of their observed heightened risk in terms of transmissibility, severity, vaccine or immune evasion, and have been renamed Alpha, Beta, Gamma and Delta, respectively. It is essential that these variants and any newly emerging variants are observed and monitored by genomic surveillance in LMICs as well as in high-income countries (HICs). This is partly enabled through international data sharing on global databases such as GISAID.

To inform health policies for SARS-CoV-2 control at the national level, we combined viral genomics and population mobility data to analyse the emergence and outbreak trajectory of SARS-CoV-2 in Bangladesh.

July 2020. SARS-CoV-2-positive samples were collected by the IEDCR from patients in six administrative areas between March and July 2020 (Supplementary Table 1 ). Samples were tested for SARS-CoV-2 using RT-PCR and 67 positive samples were transferred to the Institute for Developing Science and Health Initiatives for sequencing. We combined our data with 324 publicly available genomes from Bangladesh that were sampled during the first wave (Supplementary Table 2 , acknowledged in the Supplementary Information) and added 68,870 international genomes from GISAID for phylogeographical context (Fig. 1a) .

The 391 isolates from Bangladesh fell into 19 lineages, assigned by the Pangolin lineage software 4 Table 3 . Bayesian phylodynamic analysis estimated that the mutation rate of the isolates from Bangladesh is 0.7 × 10 −3 substitutions per site per year (~20 mutations per genome per year), which is consistent with global estimates 5 . Phylodynamic analysis of the most recent common ancestor (MRCA) predicted that SARS-CoV-2 first appeared in mid-February 2020 in Bangladesh (Fig. 1b) , consistent with other global estimates 5 . The first positive case was detected in early March, within two weeks of the predicted first introduction. (Fig. 2) . None of these lineages shared an MRCA with other lineages present only in Bangladesh, rather, they were all derived from lineages that were established and circulating outside Bangladesh and that therefore represent separate introductions from the global outbreak pool (Fig. 1) . B.1.1.25 has also been detected in the United Kingdom and Australia. Our data suggest that lineage B.1.1.25 was imported into Bangladesh at least twice before the cessation of international air travel on 21 March (Fig. 1a,b) . B.1.1 seems to have been imported into Bangladesh at least five times, from both the United States and Europe (Fig. 1a) .

In contrast to lineages B. Lineage A, which was present at the beginning of the pandemic in China in 2019 before the introduction into Europe in 2020, was not detected in Bangladesh until April 2020 (Fig. 1b) , did not expand or disseminate widely, and may have been introduced by widespread dispersal in India before being detected in Bangladesh (Fig. 1a) .

Mobility data from Facebook users were available from 22 March, and from mobile phone operators from 27 April. An analysis of the mobility patterns among approximately 6,600 Facebook users (Fig. 3a) shows that around 14.2% of users left Dhaka between 23 and 26 March, indicating a mass migration out of Dhaka to all areas of the country. The displacement of individuals from the city to other parts of the country is illustrated by the increase in the population relative to the baseline in rural areas such as Barisal. Note the large spike of movement associated with the evacuation of Barisal due to Cyclone Amphan in May 13 . Interestingly, from Fig. 3a , in areas such as Gazipur and Narayanganj, it is apparent that, after initial displacement, people returned in late April 14 ; our data show that this was followed by a return to Dhaka in late May (Fig. 3a) . Figure 3b illustrates the distribution of journey lengths, which are a proxy for movement patterns, taken by approximately 100 million mobile phone subscribers from three of four Bangladeshi operators from May to August. The impact of Eid-a national holiday during which many people travel to visit family and friends-is visible in July, with more long-distance trips occurring compared with other months. The percentage of long-distance trips (>50 km) travelled during the month of Eid (July) was significantly different compared with that of other months-the percentage of long-distance trips in July was 71.1% (95% confidence interval (CI) = 71.0-71.2%); the percentage of long-distance trips in August, the next highest month, was 57.9% (95% CI = 57.9-58.0%). In general, the distances travelled are considerable due to the country's geography (Bangladesh spans 600 km east to west) and are consistent with Bangladesh's highly mobile populations. Figure 3c shows the large number of subscribers travelling to and from Narayanganj, one of the first hotspots of SARS-CoV-2 and a continuing driver of transmission, throughout the summer. Together, these mobility data are consistent with the rapid dissemination of SARS-CoV-2 out of Dhaka to the rest of the country, as people left the city at the end of March, and frequent travel around Bangladesh sustaining transmission subsequently.

Mass migration was the main driver of countrywide spread in 2020. Combining two data streams (population mobility and genomics data) revealed a link between the expansion of three dominant lineages and the mass migration that occurred at the end of March at the beginning of the general holiday. This is evidenced by the maximum clade credibility tree in Fig. 1b that dates the MRCA for all three dominant lineages directly before the observed mass migration when the mobility data usage swapped from cities to regional towns (Fig. 3) . The resulting expansion of lineages B.1.1, B.1.1.25 and B.1.36 is shown in Fig. 2 , which shows that these dominant lineages outcompeted other lineages, and reveals the rapid dissemination of these lineages to most areas of the country after the mass migration event at the end of March 2020.

Given the daily migration of the Bangladeshi population into and out of major cities, the risk for sustained transmission and expansion of SARS-CoV-2 remains extremely high and limits the effectiveness of behavioural interventions. This highly mobile population seems to have rapidly transported lineages across the country. This occurred mainly when workers living and working in cities returned home to rural areas because schools, offices and other working places closed.

To establish whether there had been recent changes in lineage dynamics or international importations, we sequenced a further 85 SARS-CoV-2 samples (Supplementary To investigate the emergence of lineages of concern Alpha (B.1.1.7) and Beta (B.1.351) in Bangladesh, a Nextstrain phylogenetic tree was built (Fig. 4) . This revealed continuous transmission of lineage B.1.1.25 from the first wave and showed that B.1.1 and B.1.36 had been eliminated. Multiple introductions of Alpha B.1.1.7 occurred in Bangladesh after December 2020, but most were not established. More recently, in February 2021, a transmission chain became established in the community that resulted in non-travel-associated cases. Most concerning is the recent dominance of Beta (B.1.351), which has been assigned to a large number of cases that were not associated with travel and is responsible for 47% of our randomly sampled and recently sequenced cases. More than 85% of our Beta (B.1.351) sequences were sampled in Dhaka. Several other researchers who are sequencing viruses sampled in Bangladesh have deposited Beta (B.1.351) lineage sequences in GISAID (384 as of July 2021), indicating a lineage expansion and confirming that the dominance in our random sample is not a result of sampling bias. Further surveillance is necessary to investigate how this imported lineage has established and dominated the epidemic in Bangladesh so rapidly and what effect this lineage is having on transmissibility, mortality and immunity there.

In 2020, we sequenced 67 viral genomes from six administrative areas (divisions) of Bangladesh using a nanopore MinION device and combined these sequencing data with already-deposited global SARS-CoV-2 genomes (including 324 from Bangladesh COVID-19 cases) to show how early repeated international introductions into Bangladesh were replaced by endemic spread of three dominant lineages that dispersed the country in late March. Population mobility patterns analysed from digital trace and mobile phone data showed that the switch in the dynamics of SARS-CoV-2 genomic epidemiology coincided with migration out of cities to the rest of the country. In 2021, we again used nanopore sequencing to provide complete genome sequences of an additional 85 samples that were isolated between November 2020 and April 2021. The 2021 results revealed the presence of VOC Beta (B.1.351) in Bangladesh and its dominance in Dhaka, the capital of Bangladesh.

Our analyses of the genomic epidemiology of SARS-CoV-2 and population mobility in Bangladesh during the first wave of 2020 indicate that repeated international importations until late March were followed by a period of sustained community transmission that was consistent with a mass exodus from urban areas. A ban on international and domestic flights (21 March 2020) was probably the end point of international importations and 23 March marked the beginning of the expansion and dispersal of three dominant lineages in Bangladesh. The predicted introduction of SARS-CoV-2 into Bangladesh during mid-February is consistent with other global estimates that predicted that lineage B.1 was imported into Italy in late January but not detected until mid-February 4,5 .

Recent studies have linked a ~50 kb genomic locus inherited by modern humans from Neanderthals to a heightened risk of severe disease from COVID- 19 (refs. 15,16 ) . Importantly, the highest carrier frequency of this high-risk haplotype is found in Bangladesh, underlining the need for the continued strengthening of existing surveillance, monitoring systems and interventions to reduce transmission. To date, published studies reporting the effect of non-pharmaceutical interventions on lineage diversity have focused on HIC settings 17 . Although only a few studies have investigated the intricacies and specificities of differences in the dissemination of SARS-CoV-2 in low-income settings, those that have revealed that transmission dynamics can be completely different from HIC settings 18 .

Evidence from our study has revealed that stay-at-home orders can exacerbate transmission in LMICs with similar demographic specificities to Bangladesh by inducing mass migrations out of cities and return of transient workers to rural villages once work opportunities shut. Our finding, that non-pharmaceutical interventions (such as stay at home orders) that have been successful in HIC settings can't always be applied identically in LMICs, is important and will continue to be relevant in light of the emergence and dominance of VOCs.

It is important to note that the Oxford-AstraZeneca vaccine was rolled out in Bangladesh from and reduce the possibility of countrywide spread. Our study exemplifies the advantages of combining mobility and genomic data to untangle outbreak dynamics to shape policies and interventions as the outbreak spreads. Here we provide evidence that policies that limit intracountry travel could reduce the incidence of VOCs in more rural areas of Bangladesh and should be swiftly implemented to prevent a rapid rise in cases outside of Dhaka.

Throughout the Bangladesh epidemic, our genomic epidemiology consortium has been in constant contact with key policy makers in Bangladesh. Since March 2020, we have been actively collaborating with national and international institutes on SARS-CoV-2 genomics and have provided genomic information to the Directorate General of Health Services (DGHS) and Ministry of Health and Family Welfare (MoFHW), Bangladesh. This has led to important policy decisions in response to the detection of the Alpha variant in early 2021 and Beta variant a few months later. For example, the government immediately quarantined patients who were positive for the SARS-CoV-2 Alpha variant (B. 1.1.7) , including travellers and their contacts, and enforced mandatory institutional quarantine for 14 days from 29 December 2020, for any incoming passengers from high-risk countries to Bangladesh.

Then, a few months later, when the Beta variant (B.1.351) emerged in Bangladesh, the government immediately quarantined infected individuals and their contacts, imposed restrictions on intercity movement and successfully limited the transmission across Bangladesh.

Moreover, in an important differentiation from interventions that were in place during the first wave, the government of Bangladesh also imposed a restriction on intercity movement (beginning June 2021) to prevent mass migration disseminating SARS-CoV-2 and VOCs across Bangladesh.

The global incidence of VOCs has continued to rise linearly throughout 2021, and lessons learnt during initial waves in 2020 (such as spread through mass migration) must be applied to stem the curve of the ongoing pandemic and increasingly transmissible variants 20 . Our study shows what can be achieved locally as a result of international collaboration using continuous genomic surveillance. Our consortium is committed to continuing to provide information on SARS-CoV-2 as the situation develops to give the required support to the Government of Bangladesh as well as to any other similar country settings.

We hope that our findings on outbreak dynamics in Bangladesh will help countries with similar healthcare challenges to develop research capabilities to monitor outbreaks and inform national social and healthcare policies to suppress the transmission of SARS-CoV-2.

Sample collection. As of 12 July 2021, the total number of COVID-19 testing facilities across Bangladesh was 613 (ref. 3 ). IEDCR was the first institute in the country to start testing for SARS-CoV-2 using RT-PCR. Thus, samples received at the IEDCR were available for sequencing from the start of the outbreak. In 2020, we sequenced 67 SARS-CoV-2-positive samples from between 8 March and 5 July. Furthermore, for our analysis, we included 324 SARS-CoV-2 strains from Bangladesh and 68,870 global SARS-CoV-2 sequences that had been deposited at the GISAID as of 31 July 2020, acknowledged in the Supplementary Information. In 2021, we sequenced an additional 85 SARS-CoV-2-positive samples that were obtained between 11 November 2020 and 15 April 2021.

Whole-genome sequencing. Patient samples that were positive for SARS-CoV-2 on the basis of RT-PCR analysis were selected for sequencing. Stored nasopharyngeal swabs were re-extracted using the QIAamp Viral RNA Mini Kit (QIAGEN) and confirmed with RT-PCR using the WHO recommended primers and probes targeting the E and N genes. Samples with C t values of less than 31 were retained for further analysis. The initial 67 samples were sequenced according to the Arctic nCoV-2019 sequencing protocol v2 (GunIt) (https://www.protocols. io/view/ncov-2019-sequencing-protocol-v2-bdp7i5rn?version_warning=no), while the later 85 samples were sequenced according to the Arctic nCoV-2019 sequencing protocol v3 (LoCost) (https://www.protocols.io/view/ncov-2019-sequ encing-protocol-v3-locost-bh42j8ye). In brief, viral cDNA was synthesized using either SuperScript IV (Thermo Fisher Scientific) or the LunaScript RT SuperMix Kit (New England BioLabs) followed by second-strand synthesis using the Q5 high-fidelity DNA polymerase (New England BioLabs). Sequencing libraries were then constructed using the Oxford Nanopore ligation sequencing kit (SQK-LSK109). Libraries were sequenced using R9.4.1 MinION flow cells. Sequenced genomes were recovered by mapping pass quality reads against the reference Wuhan genome (GenBank: MN908947.3) using the ARTIC medaka (ARTIC-nCoV-bioinformaticsSOP-v.1.1.0) pipeline for sequence correction, single-nucleotide polymorphism calling and generation of the consensus sequence. We used high coverage (200×) to account for random sequencing errors, and regions with low coverage or low quality were masked using Ns. Lineages were assigned according to the proposed nomenclature of Pangolin lineage assignment software v.2.0.7 and lineage version 2020-08-29 (https://github.com/hCoV-2019/ pangolin).

Phylogenetic analysis. SARS-CoV-2 sequences from Bangladesh, deposited at the GISAID before 31 July 2020 (n = 391), were phylogenetically placed onto a global tree of 68,870 SARS-CoV-2 global sequences available from GISAID (31 July 2020). The total available (78,448) was filtered for those with short/ambiguous sequences to leave 68,870 sequences. The global tree was provided by R. Lanfear in his 31 July 2020 update 21 . The alignment of 68,870 sequences and the global tree was used with Llama software v.0.1 to place the 324 Bangladeshi sequenced strains in the global phylogenetic context 22 . Llama parameters, including a selection of five lineage representatives and an extraction of six from the larger tree radius, were used to investigate the phylogeography of the Bangladesh isolates. We used TempEst v.1.5.3 to investigate the evolutionary tempo of the 2020 samples 23 (Supplementary Fig. 1 ). We found a positive correlation with the temporal signal and a low to moderate association between genetic distances and sampling dates (R 2 = 0.2471), reflecting the low mutation rate of SARS-CoV-2, which is consistent with findings reported elsewhere. Given the positive correlation, we determined that these data were suitable for phylogenetic molecular clock analysis in BEAST 24 . We used a time-aware coalescent Bayesian exponential growth model that is available in BEAST (v. 1.10.4 ). The HKY+Γ model of nucleotide substitution was used with a strict molecular clock. Parameters were estimated using the Bayesian Markov Chain Monte Carlo framework, with 100,000,000-step-long chains, sampling every 1,000 steps and removing the initial 10% as burn-in. Sufficient sampling was assessed using Tracer (v.1.7.1), by verifying that every parameter had effective sampling sizes of more than 100. The resulting phylogenies were visualized with a maximum clade credibility tree in FigTree (v. 1.4.4) .

To update our phylogenetic analysis with our recent 2021 samples, we constructed a Nextstrain 25 build using their global subsampling scheme across 1,045,491 sequences downloaded from GISAID on 20 April 20 2021. We selected one representative sequence from each country, year and month from the downloaded global sample. This resulted in 1,317 background global strains to 166 Bangladesh sequences, 152 provided by this study and 14 provided for background context and acknowledged in the Supplementary Information.

Anonymized and aggregated daily population location data of Facebook users in Bangladesh were provided by Facebook Data for Good 26 . These data capture Facebook users who provide location information through the Facebook app by having location services enabled. The population of Facebook users is determined by the modal location for each individual during every eight-hour time window. The data analysed here were aggregated temporally (daily) and spatially (districts) and represent the daily count of Facebook users in each district from 22 March 2020 onwards. We compared these counts to the average number of users in each location during a 45-day baseline period preceding 22 March. The baseline averages were calculated for each unique day of the week and time of day combination. For each district, we calculated the daily percentage change in Facebook users compared with the corresponding baseline average.

We derived population mobility estimates from mobile phone call detail records provided by three of the four telecommunication operators (Grameenphone, Banglalink and Robi Axiata Limited) in Bangladesh. Trips were calculated based on changes in a subscriber's assigned tower location from the previous day. All data were aggregated temporally (daily) and spatially (Upazilas; sub-district) on the basis of tower locations according to previously described methods 11, 27, 28 . The aggregated data consist of the daily number of subscribers for each Upazila and the total number of trips between all pairs of Upazilas from 27 April 2020 onwards. Although the Facebook data provide greater temporal coverage, they represent a much smaller percentage of the Bangladeshi population. The mobile phone data represent around 100 million subscribers across Bangladesh and enable us to estimate the daily number of trips each month, including the Eid holiday period, at the end of May 2020.

No statistical method was used to predetermine sample size and no data were excluded from the analyses unless C t > 31.

Reporting Summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

All sequencing data used in this study are available on GISAID, as described in the Supplementary Information and in the GISAID acknowledgements (Supplementary Data 1). Sequencing reads generated in this study have been made available on the NCBI Sequence Read Archive service (BioProject ID: PRJNA737194).

Corresponding author(s): Firdausi Qadri Last updated by author(s): Jul 23, 2021 Reporting Summary Nature Portfolio wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Portfolio policies, see our Editorial Policies and the Editorial Policy Checklist.

For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement A statement on whether measurements were taken from distinct samples or whether the same sample was measured repeatedly

The statistical test(s) used AND whether they are one-or two-sided Only common tests should be described solely by name; describe more complex techniques in the Methods section.

A description of all covariates tested A description of any assumptions or corrections, such as tests of normality and adjustment for multiple comparisons A full description of the statistical parameters including central tendency (e.g. means) or other basic estimates (e.g. regression coefficient) AND variation (e.g. standard deviation) or associated estimates of uncertainty (e.g. confidence intervals)

For null hypothesis testing, the test statistic (e.g. F, t, r) with confidence intervals, effect sizes, degrees of freedom and P value noted For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio guidelines for submitting code & software for further information.

Policy information about availability of data All manuscripts must include a data availability statement. This statement should provide the following information, where applicable:

-Accession codes, unique identifiers, or web links for publicly available datasets -A description of any restrictions on data availability -For clinical datasets or third party data, please ensure that the statement adheres to our policy Data availability: All consensus genomes are available from the GISAID database and sequencing reads from the NCBI Sequence Read Archive service (BioProject ID: PRJNA737194).

Field-specific reporting Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection. All studies must disclose on these points even when the disclosure is negative.

In Bangladesh genomics, social media and mobile phone data streams are integrated to map the spread of SARS-CoV-2 lineages and inform country-level policies to curb infection rates

Rural population (% of total population)-Bangladesh

ZS?locations=BD 2. Poverty data: Bangladesh (Asian Development Bank

WHO Bangladesh COVID-19 Morbidity and Mortality Weekly Update No. 72

A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology

The emergence of SARS-CoV-2 in Europe and North America

Genomic epidemiology reveals transmission patterns and dynamics of SARS-CoV-2 in Aotearoa New Zealand

Complete genome sequence of a novel coronavirus (SARS-CoV-2) isolate from

Coding-complete genome sequences of three SARS-CoV-2 strains from

Genome sequence of a SARS-CoV-2 strain from Bangladesh that is nearly identical to United Kingdom SARS-CoV-2 variant B.1.1.7. Microbiol

COVID-19 rise in Bangladesh correlates with increasing detection of B.1.351 variant

Mobile phone data for public health: towards data-sharing solutions that protect individual privacy and national security

Association between mobility patterns and COVID-19 transmission in the USA: a mathematical modelling study

Amphan: India and Bangladesh evacuate millions ahead of super cyclone

Bangladesh reopens 600 apparel factories despite virus risk

Genomewide association study of severe Covid-19 with respiratory failure

The major genetic risk factor for severe COVID-19 is inherited from Neanderthals

Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK

Epidemiology and transmission dynamics of COVID-19 in two Indian states

Cellular immune responses are preserved and may contribute to Chadox1 ChAdOx1 nCoV-19 vaccine effectiveness against infection due to SARS-CoV-2 B·1·617·2 Delta variant despite reduced virus neutralisation

SARS-CoV-2 Lineages

A global phylogeny of SARS-CoV-2 from GISAID data, including sequences deposited up to 31

Llama Software: local lineage and monophyly assessment v0

Exploring the temporal structure of heterochronous sequences using TempEst

EAST 2.5: an advanced software platform for Bayesian evolutionary analysis

Nextstrain: real-time tracking of pathogen evolution

Facebook Disaster Maps: aggregate insights for crisis response & recovery

Multinational patterns of seasonal asymmetry in human movement influence infectious disease dynamics

Impact of human mobility on the emergence of dengue epidemics in Pakistan

We dedicate this work to Professor Dr Salehin Qadri for all his effort and support in making this study, and studies like this, possible. The Bill and Melinda Gates foundation, Wellcome Trust (grant 206194) and the Government of Bangladesh supported the in-country sequencing of SARS-CoV-2 Bangladesh samples. We would like to acknowledge Sally Forrest for her efforts in shipping reagents during the pandemic, which enabled the sequencing detailed in this report. We would also like to acknowledge diagnostic testing staff at IEDCR and ideSHI who provided initial diagnostics of SARS-CoV-2 samples. icddr,b is grateful to the Governments of Bangladesh, Canada, Sweden and the UK for providing core/unrestricted support.

We provide 152 whole genome sequences of SARS-CoV-2 sampled in Bangladesh between March 2020 and April 2021. All sequences are indepent and sampled from a clinical case of COVID-19.

No statistical method was used to predetermine sample size and no data were excluded from the analyses unless ct value <31.

All data was collected in Bangladesh using Oxford Nanopore Minion sequencing equipment at ideSHi.Timing and spatial scale Data was collected during March 2020-April 2021.

Samples which had ct values <31 were excluded from sequencing and the study.

All methods are fully reproducible according to the detailed instructions in methods and fully accesible sequencing data in Bioproject PRJNA737194. 

had complete access to all data in the study, and T.S. and F.Q. were responsible for the decision to submit the manuscript for publication.

The authors declare no competing interests.

The online version contains supplementary material available at https://doi.org/10.1038/s41564-021-00955-3.Correspondence and requests for materials should be addressed to F.Q.Peer review information Nature Microbiology thanks Francesco Beghini and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.