Person:

Buckee, Caroline

The var gene family of Plasmodium falciparum encodes the immunodominant variant surface antigens PfEMP1. These highly polymorphic proteins are important virulence factors that mediate cytoadhesion to a variety of host tissues, causing sequestration of parasitized red blood cells in vital organs, including the brain or placenta. Acquisition of variant-specific antibodies correlates with protection against severe malarial infections; however, understanding the relationship between gene expression and infection outcome is complicated by the modular genetic architectures of var genes that encode varying numbers of antigenic domains with differential binding specificities. By analyzing the domain architectures of fully sequenced var gene repertoires we reveal a significant, non-random association between the number of domains comprising a var gene and their sequence conservation. As such, var genes can be grouped into those that are short and diverse and genes that are long and conserved, suggesting gene length as an important characteristic in the classification of var genes. We then use an evolutionary framework to demonstrate how the same evolutionary forces acting on the level of an individual gene may have also shaped the parasite's gene repertoire. The observed associations between sequence conservation, gene architecture and repertoire structure can thus be explained by a trade-off between optimizing within-host fitness and minimizing between-host immune selection pressure. Our results demonstrate how simple evolutionary mechanisms can explain var gene structuring on multiple levels and have important implications for understanding the multifaceted epidemiology of P. falciparum malaria.

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.

Background: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73% of reported poliomyelitis in 2015 and 54% in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. Methods and findings We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95% confidence interval [CI] 0.67–0.84; and OR = 0.75, 95% CI 0.66–0.85, respectively, for each 10% increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95% CI 1.02–1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013–2016 showed good predictive ability (area under the curve range: 0.76–0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July–December 2016 and January–June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. Conclusions: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.

BACKGROUND Quantifying the effect of natural disasters on society is critical for recovery of public health services and infrastructure. The death toll can be difficult to assess in the aftermath of a major disaster. In September 2017, Hurricane Maria caused massive infrastructural damage to Puerto Rico, but its effect on mortality remains contentious. The official death count is 64.

METHODS Using a representative, stratified sample, we surveyed 3299 randomly chosen households across Puerto Rico to produce an independent estimate of all-cause mortality after the hurricane. Respondents were asked about displacement, infrastructure loss, and causes of death. We calculated excess deaths by comparing our estimated post-hurricane mortality rate with official rates for the same period in 2016.

RESULTS From the survey data, we estimated a mortality rate of 14.3 deaths (95% confidence interval [CI], 9.8 to 18.9) per 1000 persons from September 20 through December 31, 2017. This rate yielded a total of 4645 excess deaths during this period (95% CI, 793 to 8498), equivalent to a 62% increase in the mortality rate as compared with the same period in 2016. However, this number is likely to be an underestimate because of survivor bias. The mortality rate remained high through the end of December 2017, and one third of the deaths were attributed to delayed or interrupted health care. Hurricane-related migration was substantial.

CONCLUSIONS This household-based survey suggests that the number of excess deaths related to Hurricane Maria in Puerto Rico is more than 70 times the official estimate. (Funded by the Harvard T.H. Chan School of Public Health and others.)

With the rapidly increasing abundance and accessibility of genomic data, there is a growing interest in using population genetic approaches to characterize fine-scale dispersal of organisms, providing insight into biological processes across a broad range of fields including ecology, evolution and epidemiology. For sexually recombining haploid organisms such as the human malaria parasite P. falciparum, however, there have been no systematic assessments of the type of data and methods required to resolve fine scale connectivity. This analytical gap hinders the use of genomics for understanding local transmission patterns, a crucial goal for policy makers charged with eliminating this important human pathogen. Here we use data collected from four clinics with a catchment area spanning approximately 120 km of the Thai-Myanmar border to compare the ability of divergence (FST) and relatedness based on identity by descent (IBD) to resolve spatial connectivity between malaria parasites collected from proximal clinics. We found no relationship between inter-clinic distance and FST, likely due to sampling of highly related parasites within clinics, but a significant decline in IBD-based relatedness with increasing inter-clinic distance. This association was contingent upon the data set type and size. We estimated that approximately 147 single-infection whole genome sequenced parasite samples or 222 single-infection parasite samples genotyped at 93 single nucleotide polymorphisms (SNPs) were sufficient to recover a robust spatial trend estimate at this scale. In summary, surveillance efforts cannot rely on classical measures of genetic divergence to measure P. falciparum transmission on a local scale. Given adequate sampling, IBD-based relatedness provides a useful alternative, and robust trends can be obtained from parasite samples genotyped at approximately 100 SNPs.

Background: Oral cholera vaccination is an approach to preventing outbreaks in at-risk settings and controlling cholera in endemic settings. However, vaccine-derived herd immunity may be short-lived due to interactions between human mobility and imperfect or waning vaccine efficacy. As the supply and utilization of oral cholera vaccines grows, critical questions related to herd immunity are emerging, including: who should be targeted; when should revaccination be performed; and why have cholera outbreaks occurred in recently vaccinated populations? Methods and findings We use mathematical models to simulate routine and mass oral cholera vaccination in populations with varying degrees of migration, transmission intensity, and vaccine coverage. We show that migration and waning vaccine efficacy strongly influence the duration of herd immunity while birth and death rates have relatively minimal impacts. As compared to either periodic mass vaccination or routine vaccination alone, a community could be protected longer by a blended “Mass and Maintain” strategy. We show that vaccination may be best targeted at populations with intermediate degrees of mobility as compared to communities with very high or very low population turnover. Using a case study of an internally displaced person camp in South Sudan which underwent high-coverage mass vaccination in 2014 and 2015, we show that waning vaccine direct effects and high population turnover rendered the camp over 80% susceptible at the time of the cholera outbreak beginning in October 2016. Conclusions: Oral cholera vaccines can be powerful tools for quickly protecting a population for a period of time that depends critically on vaccine coverage, vaccine efficacy over time, and the rate of population turnover through human mobility. Due to waning herd immunity, epidemics in vaccinated communities are possible but become less likely through complementary interventions or data-driven revaccination strategies.

Seasonal variation in human mobility is globally ubiquitous and affects the spatial spread of infectious diseases, but the ability to measure seasonality in human movement has been limited by data availability. Here, we use mobile phone data to quantify seasonal travel and directional asymmetries in Kenya, Namibia, and Pakistan, across a spectrum from rural nomadic populations to highly urbanized communities. We then model how the geographic spread of several acute pathogens with varying life histories could depend on country-wide connectivity fluctuations through the year. In all three countries, major national holidays are associated with shifts in the scope of travel. Within this broader pattern, the relative importance of particular routes also fluctuates over the course of the year, with increased travel from rural to urban communities after national holidays, for example. These changes in travel impact how fast communities are likely to be reached by an introduced pathogen.

Climate change is likely to profoundly modulate the burden of infectious diseases. However, attributing health impacts to a changing climate requires being able to associate changes in infectious disease incidence with the potentially complex influences of climate. This aim is further complicated by nonlinear feedbacks inherent in the dynamics of many infections, driven by the processes of immunity and transmission. Here, we detail the mechanisms by which climate drivers can shape infectious disease incidence, from direct effects on vector life history to indirect effects on human susceptibility, and detail the scope of variation available with which to probe these mechanisms. We review approaches used to evaluate and quantify associations between climate and infectious disease incidence, discuss the array of data available to tackle this question, and detail remaining challenges in understanding the implications of climate change for infectious disease incidence. We point to areas where synthesis between approaches used in climate science and infectious disease biology provide potential for progress.