Close Proximity to Mining Is Associated with Increased Prevalence of the Drug Resistance-Associated Mutation dhps540E in Eastern Democratic Republic of the Congo.

Increasing sulfadoxine-pyrimethamine (SP) resistance in the Democratic Republic of the Congo (DRC) has threatened its use for prevention of malaria in one of the most malarious countries in the world. Using geographic information on mining operations in the DRC and genetic data on SP drug resistance markers from the 2013-2014 Demographic and Health Surveys, we evaluated associations between close residence to mining and the presence of mutations conferring resistance to sulfadoxine. Close residential proximity to mining was associated with increased prevalence odds ratio (POR) of the dhps540E mutation (POR: 2.11, 95% uncertainty interval: 1.15-3.96) with adjustments for confounding variables and space. Our findings indicate that exposure to mining is associated with increased presence of an antimalarial drug resistance haplotype that threatens effective use of SP for vulnerable populations. Areas actively engaged in mining could be considered for interventions to reduce the spread of emerging drug resistance in the DRC.

Strong isolation by distance and evidence of population microstructure reflect ongoing Plasmodium falciparum transmission in Zanzibar.

Background: The Zanzibar archipelago of Tanzania has become a low-transmission area for Plasmodium falciparum. Despite being considered an area of pre-elimination for years, achieving elimination has been difficult, likely due to a combination of imported infections from mainland Tanzania and continued local transmission. Methods: To shed light on these sources of transmission, we applied highly multiplexed genotyping utilizing molecular inversion probes to characterize the genetic relatedness of 282 P. falciparum isolates collected across Zanzibar and in Bagamoyo district on the coastal mainland from 2016 to 2018. Results: Overall, parasite populations on the coastal mainland and Zanzibar archipelago remain highly related. However, parasite isolates from Zanzibar exhibit population microstructure due to the rapid decay of parasite relatedness over very short distances. This, along with highly related pairs within shehias, suggests ongoing low-level local transmission. We also identified highly related parasites across shehias that reflect human mobility on the main island of Unguja and identified a cluster of highly related parasites, suggestive of an outbreak, in the Micheweni district on Pemba island. Parasites in asymptomatic infections demonstrated higher complexity of infection than those in symptomatic infections, but have similar core genomes. Conclusions: Our data support importation as a main source of genetic diversity and contribution to the parasite population in Zanzibar, but they also show local outbreak clusters where targeted interventions are essential to block local transmission. These results highlight the need for preventive measures against imported malaria and enhanced control measures in areas that remain receptive to malaria reemergence due to susceptible hosts and competent vectors. Funding: This research was funded by the National Institutes of Health, grants R01AI121558, R01AI137395, R01AI155730, F30AI143172, and K24AI134990. Funding was also contributed from the Swedish Research Council, Erling-Persson Family Foundation, and the Yang Fund. RV acknowledges funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 program supported by the European Union. RV also acknowledges funding by Community Jameel.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing.

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Strong isolation by distance and evidence of population microstructure reflect ongoing Plasmodium falciparum transmission in Zanzibar.

The Zanzibar archipelago of Tanzania has become a low-transmission area for Plasmodium falciparum. Despite being considered an area of pre-elimination for years, achieving elimination has been difficult, likely due to a combination of imported infections from mainland Tanzania, and continued local transmission. To shed light on these sources of transmission, we applied highly multiplexed genotyping utilizing molecular inversion probes to characterize the genetic relatedness of 282 P. falciparum isolates collected across Zanzibar and in Bagamoyo District on the coastal mainland from 2016-2018. Overall, parasite populations on the coastal mainland and Zanzibar archipelago remain highly related. However, parasite isolates from Zanzibar exhibit population microstructure due to rapid decay of parasite relatedness over very short distances. This, along with highly related pairs within shehias, suggests ongoing low level local transmission. We also identified highly related parasites across shehias that reflect human mobility on the main island of Unguja and identified a cluster of highly related parasites, suggestive of an outbreak, in the Micheweni district on Pemba island. Parasites in asymptomatic infections demonstrated higher complexity of infection than those in symptomatic infections, but have similar core genomes. Our data support importation as a main source of genetic diversity and contribution to the parasite population on Zanzibar, but they also show local outbreak clusters where targeted interventions are essential to block local transmission. These results highlight the need for preventive measures against imported malaria and enhanced control measures in areas that remain receptive for malaria reemergence due to susceptible hosts and competent vectors.

Sampling for malaria molecular surveillance.

Strategic use of Plasmodium falciparum genetic variation has great potential to inform public health actions for malaria control and elimination. Malaria molecular surveillance (MMS) begins with a strategy to identify and collect parasite samples, guided by public-health priorities. In this review we discuss sampling design practices for MMS and point out epidemiological, biological, and statistical factors that need to be considered. We present examples for different use cases, including detecting emergence and spread of rare variants, establishing transmission sources and inferring changes in malaria transmission intensity. This review will potentially guide the collection of samples and data, serve as a starting point for further methodological innovation, and enhance utilization of MMS to support malaria elimination.

Epidemiological drivers of transmissibility and severity of SARS-CoV-2 in England.

As the SARS-CoV-2 pandemic progressed, distinct variants emerged and dominated in England. These variants, Wildtype, Alpha, Delta, and Omicron were characterized by variations in transmissibility and severity. We used a robust mathematical model and Bayesian inference framework to analyse epidemiological surveillance data from England. We quantified the impact of non-pharmaceutical interventions (NPIs), therapeutics, and vaccination on virus transmission and severity. Each successive variant had a higher intrinsic transmissibility. Omicron (BA.1) had the highest basic reproduction number at 8.3 (95% credible interval (CrI) 7.7-8.8). Varying levels of NPIs were crucial in controlling virus transmission until population immunity accumulated. Immune escape properties of Omicron decreased effective levels of immunity in the population by a third. Furthermore, in contrast to previous studies, we found Alpha had the highest basic infection fatality ratio (2.9%, 95% CrI 2.7-3.2), followed by Delta (2.2%, 95% CrI 2.0-2.4), Wildtype (1.2%, 95% CrI 1.1-1.2), and Omicron (0.7%, 95% CrI 0.6-0.8). Our findings highlight the importance of continued surveillance. Long-term strategies for monitoring and maintaining effective immunity against SARS-CoV-2 are critical to inform the role of NPIs to effectively manage future variants with potentially higher intrinsic transmissibility and severe outcomes.

coiaf: Directly estimating complexity of infection with allele frequencies.

In malaria, individuals are often infected with different parasite strains. The complexity of infection (COI) is defined as the number of genetically distinct parasite strains in an individual. Changes in the mean COI in a population have been shown to be informative of changes in transmission intensity with a number of probabilistic likelihood and Bayesian models now developed to estimate the COI. However, rapid, direct measures based on heterozygosity or FwS do not properly represent the COI. In this work, we present two new methods that use easily calculated measures to directly estimate the COI from allele frequency data. Using a simulation framework, we show that our methods are computationally efficient and comparably accurate to current approaches in the literature. Through a sensitivity analysis, we characterize how the distribution of parasite densities, the assumed sequencing depth, and the number of sampled loci impact the bias and accuracy of our two methods. Using our developed methods, we further estimate the COI globally from Plasmodium falciparum sequencing data and compare the results against the literature. We show significant differences in the estimated COI globally between continents and a weak relationship between malaria prevalence and COI.

Estimating the COVID-19 infection fatality ratio accounting for seroreversion using statistical modelling.

Background: The infection fatality ratio (IFR) is a key statistic for estimating the burden of coronavirus disease 2019 (COVID-19) and has been continuously debated throughout the COVID-19 pandemic. The age-specific IFR can be quantified using antibody surveys to estimate total infections, but requires consideration of delay-distributions from time from infection to seroconversion, time to death, and time to seroreversion (i.e. antibody waning) alongside serologic test sensitivity and specificity. Previous IFR estimates have not fully propagated uncertainty or accounted for these potential biases, particularly seroreversion. Methods: We built a Bayesian statistical model that incorporates these factors and applied this model to simulated data and 10 serologic studies from different countries. Results: We demonstrate that seroreversion becomes a crucial factor as time accrues but is less important during first-wave, short-term dynamics. We additionally show that disaggregating surveys by regions with higher versus lower disease burden can inform serologic test specificity estimates. The overall IFR in each setting was estimated at 0.49-2.53%. Conclusion: We developed a robust statistical framework to account for full uncertainties in the parameters determining IFR. We provide code for others to apply these methods to further datasets and future epidemics.

Reassessing Reported Deaths and Estimated Infection Attack Rate during the First 6 Months of the COVID-19 Epidemic, Delhi, India.

India reported >10 million coronavirus disease (COVID-19) cases and 149,000 deaths in 2020. To reassess reported deaths and estimate incidence rates during the first 6 months of the epidemic, we used a severe acute respiratory syndrome coronavirus 2 transmission model fit to data from 3 serosurveys in Delhi and time-series documentation of reported deaths. We estimated 48.7% (95% credible interval 22.1%-76.8%) cumulative infection in the population through the end of September 2020. Using an age-adjusted overall infection fatality ratio based on age-specific estimates from mostly high-income countries, we estimated that just 15.0% (95% credible interval 9.3%-34.0%) of COVID-19 deaths had been reported, indicating either substantial underreporting or lower age-specific infection-fatality ratios in India than in high-income countries. Despite the estimated high attack rate, additional epidemic waves occurred in late 2020 and April-May 2021. Future dynamics will depend on the duration of natural and vaccine-induced immunity and their effectiveness against new variants.

The epidemiology of Plasmodium vivax among adults in the Democratic Republic of the Congo.

Reports of P. vivax infections among Duffy-negative hosts have accumulated throughout sub-Saharan Africa. Despite this growing body of evidence, no nationally representative epidemiological surveys of P. vivax in sub-Saharan Africa have been performed. To overcome this gap in knowledge, we screened over 17,000 adults in the Democratic Republic of the Congo (DRC) for P. vivax using samples from the 2013-2014 Demographic Health Survey. Overall, we found a 2.97% (95% CI: 2.28%, 3.65%) prevalence of P. vivax infections across the DRC. Infections were associated with few risk-factors and demonstrated a relatively flat distribution of prevalence across space with focal regions of relatively higher prevalence in the north and northeast. Mitochondrial genomes suggested that DRC P. vivax were distinct from circulating non-human ape strains and an ancestral European P. vivax strain, and instead may be part of a separate contemporary clade. Our findings suggest P. vivax is diffusely spread across the DRC at a low prevalence, which may be associated with long-term carriage of low parasitemia, frequent relapses, or a general pool of infections with limited forward propagation.

Key epidemiological drivers and impact of interventions in the 2020 SARS-CoV-2 epidemic in England.

We fitted a model of SARS-CoV-2 transmission in care homes and the community to regional surveillance data for England. Compared with other approaches, our model provides a synthesis of multiple surveillance data streams into a single coherent modeling framework, allowing transmission and severity to be disentangled from features of the surveillance system. Of the control measures implemented, only national lockdown brought the reproduction number (Rt eff) below 1 consistently; if introduced 1 week earlier, it could have reduced deaths in the first wave from an estimated 48,600 to 25,600 [95% credible interval (CrI): 15,900 to 38,400]. The infection fatality ratio decreased from 1.00% (95% CrI: 0.85 to 1.21%) to 0.79% (95% CrI: 0.63 to 0.99%), suggesting improved clinical care. The infection fatality ratio was higher in the elderly residing in care homes (23.3%, 95% CrI: 14.7 to 35.2%) than those residing in the community (7.9%, 95% CrI: 5.9 to 10.3%). On 2 December 2020, England was still far from herd immunity, with regional cumulative infection incidence between 7.6% (95% CrI: 5.4 to 10.2%) and 22.3% (95% CrI: 19.4 to 25.4%) of the population. Therefore, any vaccination campaign will need to achieve high coverage and a high degree of protection in vaccinated individuals to allow nonpharmaceutical interventions to be lifted without a resurgence of transmission.

Using syndromic measures of mortality to capture the dynamics of COVID-19 in Java, Indonesia, in the context of vaccination rollout.

BACKGROUND: As in many countries, quantifying COVID-19 spread in Indonesia remains challenging due to testing limitations. In Java, non-pharmaceutical interventions (NPIs) were implemented throughout 2020. However, as a vaccination campaign launches, cases and deaths are rising across the island. METHODS: We used modelling to explore the extent to which data on burials in Jakarta using strict COVID-19 protocols (C19P) provide additional insight into the transmissibility of the disease, epidemic trajectory, and the impact of NPIs. We assess how implementation of NPIs in early 2021 will shape the epidemic during the period of likely vaccine rollout. RESULTS: C19P burial data in Jakarta suggest a death toll approximately 3.3 times higher than reported. Transmission estimates using these data suggest earlier, larger, and more sustained impact of NPIs. Measures to reduce sub-national spread, particularly during Ramadan, substantially mitigated spread to more vulnerable rural areas. Given current trajectory, daily cases and deaths are likely to increase in most regions as the vaccine is rolled out. Transmission may peak in early 2021 in Jakarta if current levels of control are maintained. However, relaxation of control measures is likely to lead to a subsequent resurgence in the absence of an effective vaccination campaign. CONCLUSIONS: Syndromic measures of mortality provide a more complete picture of COVID-19 severity upon which to base decision-making. The high potential impact of the vaccine in Java is attributable to reductions in transmission to date and dependent on these being maintained. Increases in control in the relatively short-term will likely yield large, synergistic increases in vaccine impact.

Within-country age-based prioritisation, global allocation, and public health impact of a vaccine against SARS-CoV-2: A mathematical modelling analysis.

The worldwide endeavour to develop safe and effective COVID-19 vaccines has been extraordinary, and vaccination is now underway in many countries. However, the doses available in 2021 are likely to be limited. We extend a mathematical model of SARS-CoV-2 transmission across different country settings to evaluate the public health impact of potential vaccines using WHO-developed target product profiles. We identify optimal vaccine allocation strategies within- and between-countries to maximise averted deaths under constraints on dose supply. We find that the health impact of SARS-CoV-2 vaccination depends on the cumulative population-level infection incidence when vaccination begins, the duration of natural immunity, the trajectory of the epidemic prior to vaccination, and the level of healthcare available to effectively treat those with disease. Within a country we find that for a limited supply (doses for < 20% of the population) the optimal strategy is to target the elderly. However, with a larger supply, if vaccination can occur while other interventions are maintained, the optimal strategy switches to targeting key transmitters to indirectly protect the vulnerable. As supply increases, vaccines that reduce or block infection have a greater impact than those that prevent disease alone due to the indirect protection provided to high-risk groups. Given a 2 billion global dose supply in 2021, we find that a strategy in which doses are allocated to countries proportional to population size is close to optimal in averting deaths and aligns with the ethical principles agreed in pandemic preparedness planning.

Spatially clustered count data provide more efficient search strategies in invasion biology and disease control.

Geographic profiling, a mathematical model originally developed in criminology, is increasingly being used in ecology and epidemiology. Geographic profiling boasts a wide range of applications, such as finding source populations of invasive species or breeding sites of vectors of infectious disease. The model provides a cost-effective approach for prioritizing search strategies for source locations and does so via simple data in the form of the positions of each observation, such as individual sightings of invasive species or cases of a disease. In doing so, however, classic geographic profiling approaches fail to make the distinction between those areas containing observed absences and those areas where no data were recorded. Absence data are generated via spatial sampling protocols but are often discarded during the inference process. Here we construct a geographic profiling model that resolves these issues by making inferences via count data, analyzing a set of discrete sentinel locations at which the number of encounters has been recorded. Crucially, in our model this number can be zero. We verify the ability of this new model to estimate source locations and other parameters of practical interest via a Bayesian power analysis. We also measure model performance via real-world data in which the model infers breeding locations of mosquitoes in bromeliads in Miami-Dade County, Florida, USA. In both cases, our novel model produces more efficient search strategies by shifting focus from those areas containing observed absences to those with no data, an improvement over existing models that treat these areas equally. Our model makes important improvements upon classic geographic profiling methods, which will significantly enhance real-world efforts to develop conservation management plans and targeted interventions.

Reduction in mobility and COVID-19 transmission.

In response to the COVID-19 pandemic, countries have sought to control SARS-CoV-2 transmission by restricting population movement through social distancing interventions, thus reducing the number of contacts. Mobility data represent an important proxy measure of social distancing, and here, we characterise the relationship between transmission and mobility for 52 countries around the world. Transmission significantly decreased with the initial reduction in mobility in 73% of the countries analysed, but we found evidence of decoupling of transmission and mobility following the relaxation of strict control measures for 80% of countries. For the majority of countries, mobility explained a substantial proportion of the variation in transmissibility (median adjusted R-squared: 48%, interquartile range - IQR - across countries [27-77%]). Where a change in the relationship occurred, predictive ability decreased after the relaxation; from a median adjusted R-squared of 74% (IQR across countries [49-91%]) pre-relaxation, to a median adjusted R-squared of 30% (IQR across countries [12-48%]) post-relaxation. In countries with a clear relationship between mobility and transmission both before and after strict control measures were relaxed, mobility was associated with lower transmission rates after control measures were relaxed indicating that the beneficial effects of ongoing social distancing behaviours were substantial.

Plasmodium interspecies interactions during a period of increasing prevalence of Plasmodium ovale in symptomatic individuals seeking treatment: an observational study.

BACKGROUND: The epidemiology and severity of non-falciparum malaria in endemic settings has garnered little attention. We aimed to characterise the prevalence, interaction, clinical risk factors, and temporal trends of non-falciparum Plasmodium species among symptomatic individuals presenting at health-care facilities in endemic settings of Kenya. METHODS: We diagnosed and analysed infecting malaria species (Plasmodium falciparum, Plasmodium ovale curtisi, Plasmodium ovale wallikeri, and Plasmodium malariae) via PCR in clinical samples collected between March 1, 2008, and Dec 31, 2016, from six hospitals located in different regions of Kenya. We recruited patients aged 6 months or older who presented at outpatient departments with symptoms of malaria or tested positive for uncomplicated malaria by malaria rapid diagnostic test. Descriptive statistics were used to describe the prevalence and distribution of Plasmodium species. A statistical model was designed and used for estimating the frequency of Plasmodium species and assessing interspecies interactions. Mixed-effect linear regression models with random slopes for each location were used to test for change in prevalence over time. FINDINGS: Samples from 2027 symptomatic participants presenting at care facilities were successfully analysed for all Plasmodium species. 1469 (72·5%) of the samples were P falciparum single-species infections, 523 (25·8%) were mixed infections, and only 35 (1·7%) were single non-falciparum species infections. 452 (22·3%) were mixed infections containing P ovale spp. A likelihood-based model calculation of the population frequency of each species estimated a significant within-host interference between P falciparum and P ovale curtisi. Mixed-effect logistic regression models identified a significant increase in P ovale wallikeri (2·1% per year; p=0·043) and P ovale curtisi (0·7% per year; p=0·0002) species over time, with a reciprocal decrease in P falciparum single-species infections (2·5% per year; p=0·0065). The frequency of P malariae infections did not significantly change over time. Risk of P falciparum infections presenting with fever was lower if co-infected with P malariae (adjusted odds ratio 0·43, 95% CI 0·25-0·74; p=0·0023). INTERPRETATION: Our results show a prevalence of non-falciparum species infections of 27·5% among symptomatic individuals presenting at care facilities, which is higher than expected from previous cross-sectional surveys. The proportion of infections with P ovale wallikeri and P ovale curtisi was observed to significantly increase over the period of study, which could be due to attenuated responsiveness of these species to malaria drug treatment. The increase in frequency of P ovale spp could threaten the malaria control efforts in Kenya and pose increased risk of malaria to travellers. FUNDING: Armed Forces Health Surveillance Branch and its Global Emerging Infections Surveillance Section.

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes.

High-throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up-to-date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high-burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome-wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13-R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region-level complexity-of-infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.

Evaluating the Performance of Malaria Genetics for Inferring Changes in Transmission Intensity Using Transmission Modeling.

Substantial progress has been made globally to control malaria, however there is a growing need for innovative new tools to ensure continued progress. One approach is to harness genetic sequencing and accompanying methodological approaches as have been used in the control of other infectious diseases. However, to utilize these methodologies for malaria, we first need to extend the methods to capture the complex interactions between parasites, human and vector hosts, and environment, which all impact the level of genetic diversity and relatedness of malaria parasites. We develop an individual-based transmission model to simulate malaria parasite genetics parameterized using estimated relationships between complexity of infection and age from five regions in Uganda and Kenya. We predict that cotransmission and superinfection contribute equally to within-host parasite genetic diversity at 11.5% PCR prevalence, above which superinfections dominate. Finally, we characterize the predictive power of six metrics of parasite genetics for detecting changes in transmission intensity, before grouping them in an ensemble statistical model. The model predicted malaria prevalence with a mean absolute error of 0.055. Different assumptions about the availability of sample metadata were considered, with the most accurate predictions of malaria prevalence made when the clinical status and age of sampled individuals is known. Parasite genetics may provide a novel surveillance tool for estimating the prevalence of malaria in areas in which prevalence surveys are not feasible. However, the findings presented here reinforce the need for patient metadata to be recorded and made available within all future attempts to use parasite genetics for surveillance.