Search

The 2013–2016 epidemic of Ebola virus disease was of unprecedented magnitude, duration and impact. Analysing 1610 Ebola virus genomes, representing over 5% of known cases, we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic ‘gravity’ model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already set the seeds for an international epidemic, rendering these measures ineffective in curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing they were susceptible to significant outbreaks but at lower risk of introductions. Finally, we reveal this large epidemic to be a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help inform interventions in future epidemics., Organismic and Evolutionary Biology, Accepted Manuscript

Background;Limited clinical and laboratory data are available on patients with Ebola virus disease (EVD). The Kenema Government Hospital in Sierra Leone, which had an existing infrastructure for research regarding viral hemorrhagic fever, has received and cared for patients with EVD since the beginning of the outbreak in Sierra Leone in May 2014. Methods;We reviewed available epidemiologic, clinical, and laboratory records of patients in whom EVD was diagnosed between May 25 and June 18, 2014. We used quantitative reverse-transcriptase–polymerase-chain-reaction assays to assess the load of Ebola virus (EBOV, Zaire species) in a subgroup of patients. Results;Of 106 patients in whom EVD was diagnosed, 87 had a known outcome, and 44 had detailed clinical information available. The incubation period was estimated to be 6 to 12 days, and the case fatality rate was 74%. Common findings at presentation included fever (in 89% of the patients), headache (in 80%), weakness (in 66%), dizziness(in 60%), diarrhea (in 51%), abdominal pain (in 40%), and vomiting (in 34%). Clinical and laboratory factors at presentation that were associated with a fatal outcome included fever, weakness, dizziness, diarrhea, and elevated levels of blood urea nitrogen, aspartate aminotransferase, and creatinine. Exploratory analyses indicated that patients under the age of 21 years had a lower case fatality rate than those over the age of 45 years (57% vs. 94%, P=0.03), and patients presenting with fewer than 100,000 EBOV copies per milliliter had a lower case fatality rate than those with 10 million EBOV copies per milliliter or more (33% vs. 94%, P=0.003). Bleeding occurred in only 1 patient.Conclusions; The incubation period and case fatality rate among patients with EVD in Sierra Leone are similar to those observed elsewhere in the 2014 outbreak and in previous outbreaks. Although bleeding was an infrequent finding, diarrhea and other gastrointestinal manifestations were common. (Funded by the National Institutes of Health and others.), Organismic and Evolutionary Biology, Version of Record

Analysis of genome sequences of 159 isolates of Plasmodium falciparum from Senegal yields an extraordinarily high proportion (26.85%) of protein-coding genes with the ratio of nonsynonymous to synonymous polymorphism greater than one. This proportion is much greater than observed in other organisms. Also unusual is that the site-frequency spectra of synonymous and nonsynonymous polymorphisms are virtually indistinguishable. We hypothesized that the complicated life cycle of malaria parasites might lead to qualitatively different population genetics from that predicted from the classical Wright-Fisher (WF) model, which assumes a single random-mating population with a finite and constant population size in an organism with nonoverlapping generations. This paper summarizes simulation studies of random genetic drift and selection in malaria parasites that take into account their unusual life history. Our results show that random genetic drift in the malaria life cycle is more pronounced than under the WF model. Paradoxically, the efficiency of purifying selection in the malaria life cycle is also greater than under WF, and the relative efficiency of positive selection varies according to conditions. Additionally, the site-frequency spectrum under neutrality is also more skewed toward low-frequency alleles than expected with WF. These results highlight the importance of considering the malaria life cycle when applying existing population genetic tools based on the WF model. The same caveat applies to other species with similarly complex life cycles., Organismic and Evolutionary Biology, Accepted Manuscript

Through rapid genetic adaptation and natural selection, the Plasmodium falciparum parasite—the deadliest of those that cause malaria—is able to develop resistance to antimalarial drugs, thwarting present efforts to control it. Genome-wide association studies (GWAS) provide a critical hypothesis-generating tool for understanding how this occurs. However, in P. falciparum, the limited amount of linkage disequilibrium hinders the power of traditional array-based GWAS. Here, we demonstrate the feasibility and power improvements gained by using whole-genome sequencing for association studies. We analyzed data from 45 Senegalese parasites and identified genetic changes associated with the parasites’ in vitro response to 12 different antimalarials. To further increase statistical power, we adapted a common test for natural selection, XP-EHH (cross-population extended haplotype homozygosity), and used it to identify genomic regions associated with resistance to drugs. Using this sequence-based approach and the combination of association and selection-based tests, we detected several loci associated with drug resistance. These loci included the previously known signals at pfcrt, dhfr, and pfmdr1, as well as many genes not previously implicated in drug-resistance roles, including genes in the ubiquitination pathway. Based on the success of the analysis presented in this study, and on the demonstrated shortcomings of array-based approaches, we argue for a complete transition to sequence-based GWAS for small, low linkage-disequilibrium genomes like that of P. falciparum., Molecular and Cellular Biology, Organismic and Evolutionary Biology, Accepted Manuscript

We used a high-density single-nucleotide polymorphism array to genotype 75 Plasmodium falciparum isolates recently collected from Senegal and The Gambia to search for signals of selection in this malaria endemic region. We found little geographic or temporal stratification of the genetic diversity among the sampled parasites. Through application of the iHS and REHH haplotype-based tests for positive selection, we found evidence of recent selective sweeps at a known drug resistance locus, at several known antigenic loci, and at several genomic regions not previously identified as sites of recent selection. We discuss the value of deep population-specific genomic analyses for identifying selection signals within sampled endemic populations of parasites, which may correspond to local selection pressures such as distinctive therapeutic regimes or mosquito vectors., Version of Record

Background: Cerebral malaria, a severe form of Plasmodium falciparum infection, is an important cause of mortality in sub-Saharan African children. A Taqman 24 Single Nucleotide Polymorphisms (SNP) molecular barcode assay was developed for use in laboratory parasites which estimates genotype number and identifies the predominant genotype. Methods The 24 SNP assay was used to determine predominant genotypes in blood and tissues from autopsy and clinical patients with cerebral malaria. Results: Single genotypes were shared between the peripheral blood, the brain, and other tissues of cerebral malaria patients, while malaria-infected patients who died of non-malarial causes had mixed genetic signatures in tissues examined. Children with retinopathy-positive cerebral malaria had significantly less complex infections than those without retinopathy (OR = 3.7, 95% CI [1.51-9.10]).The complexity of infections significantly decreased over the malaria season in retinopathy-positive patients compared to retinopathy-negative patients. Conclusions: Cerebral malaria patients harbour a single or small set of predominant parasites; patients with incidental parasitaemia sustain infections involving diverse genotypes. Limited diversity in the peripheral blood of cerebral malaria patients and correlation with tissues supports peripheral blood samples as appropriate for genome-wide association studies of parasite determinants of pathogenicity., Organismic and Evolutionary Biology, Version of Record

Background: Mounting evidence suggests a major role for epigenetic feedback in Plasmodium falciparum transcriptional regulation. Long non-coding RNAs (lncRNAs) have recently emerged as a new paradigm in epigenetic remodeling. We therefore set out to investigate putative roles for lncRNAs in P. falciparum transcriptional regulation. Results: We used a high-resolution DNA tiling microarray to survey transcriptional activity across 22.6% of the P. falciparum strain 3D7 genome. We identified 872 protein-coding genes and 60 putative P. falciparum lncRNAs under developmental regulation during the parasite's pathogenic human blood stage. Further characterization of lncRNA candidates led to the discovery of an intriguing family of lncRNA telomere-associated repetitive element transcripts, termed lncRNA-TARE. We have quantified lncRNA-TARE expression at 15 distinct chromosome ends and mapped putative transcriptional start and termination sites of lncRNA-TARE loci. Remarkably, we observed coordinated and stage-specific expression of lncRNA-TARE on all chromosome ends tested, and two dominant transcripts of approximately 1.5 kb and 3.1 kb transcribed towards the telomere. Conclusions: We have characterized a family of 22 telomere-associated lncRNAs in P. falciparum. Homologous lncRNA-TARE loci are coordinately expressed after parasite DNA replication, and are poised to play an important role in P. falciparum telomere maintenance, virulence gene regulation, and potentially other processes of parasite chromosome end biology. Further study of lncRNA-TARE and other promising lncRNA candidates may provide mechanistic insight into P. falciparum transcriptional regulation., Organismic and Evolutionary Biology, Stem Cell and Regenerative Biology, Other Research Unit, Accepted Manuscript

The Plasmodium falciparum parasite's ability to adapt to environmental pressures, such as the human immune system and antimalarial drugs, makes malaria an enduring burden to public health. Understanding the genetic basis of these adaptations is critical to intervening successfully against malaria. To that end, we created a high-density genotyping array that assays over 17,000 single nucleotide polymorphisms (~1 SNP/kb), and applied it to 57 culture-adapted parasites from three continents. We characterized genome-wide genetic diversity within and between populations and identified numerous loci with signals of natural selection, suggesting their role in recent adaptation. In addition, we performed a genome-wide association study (GWAS), searching for loci correlated with resistance to thirteen antimalarials; we detected both known and novel resistance loci, including a new halofantrine resistance locus, PF10_0355. Through functional testing we demonstrated that PF10_0355 overexpression decreases sensitivity to halofantrine, mefloquine, and lumefantrine, but not to structurally unrelated antimalarials, and that increased gene copy number mediates resistance. Our GWAS and follow-on functional validation demonstrate the potential of genome-wide studies to elucidate functionally important loci in the malaria parasite genome., Organismic and Evolutionary Biology, Other Research Unit, Accepted Manuscript

Background: The identification of genetic changes that confer drug resistance or other phenotypic changes in pathogens can help optimize treatment strategies, support the development of new therapeutic agents, and provide information about the likely function of genes. Elucidating mechanisms of phenotypic drug resistance can also assist in identifying the mode of action of uncharacterized but potent antimalarial compounds identified in high-throughput chemical screening campaigns against Plasmodium falciparum. Results: Here we show that tiling microarrays can detect de novo a large proportion of the genetic changes that differentiate one genome from another. We show that we detect most single nucleotide polymorphisms or small insertion deletion events and all known copy number variations that distinguish three laboratory isolates using readily accessible methods. We used the approach to discover mutations that occur during the selection process after transfection. We also elucidated a mechanism by which parasites acquire resistance to the antimalarial fosmidomycin, which targets the parasite isoprenoid synthesis pathway. Our microarray-based approach allowed us to attribute in vitro derived fosmidomycin resistance to a copy number variation event in the pfdxr gene, which enables the parasite to overcome fosmidomycin-mediated inhibition of isoprenoid biosynthesis. Conclusions: We show that newly emerged single nucleotide polymorphisms can readily be detected and that malaria parasites can rapidly acquire gene amplifications in response to in vitro drug pressure. The ability to define comprehensively genetic variability in P. falciparum with a single overnight hybridization creates new opportunities to study parasite evolution and improve the treatment and control of malaria., Other Research Unit, Version of Record

Background: Single nucleotide polymorphism (SNP) genotyping provides the means to develop a practical, rapid, inexpensive assay that will uniquely identify any Plasmodium falciparum parasite using a small amount of DNA. Such an assay could be used to distinguish recrudescence from re-infection in drug trials, to monitor the frequency and distribution of specific parasites in a patient population undergoing drug treatment or vaccine challenge, or for tracking samples and determining purity of isolates in the laboratory during culture adaptation and sub-cloning, as well as routine passage. Methods: A panel of twenty-four SNP markers has been identified that exhibit a high minor allele frequency (average MAF > 35%), for which robust TaqMan genotyping assays were constructed. All SNPs were identified through whole genome sequencing and MAF was estimated through Affymetrix array-based genotyping of a worldwide collection of parasites. These assays create a "molecular barcode" to uniquely identify a parasite genome.Results Using 24 such markers no two parasites known to be of independent origin have yet been found to have the same allele signature. The TaqMan genotyping assays can be performed on a variety of samples including cultured parasites, frozen whole blood, or whole blood spotted onto filter paper with a success rate > 99%. Less than 5 ng of parasite DNA is needed to complete a panel of 24 markers. The ability of this SNP panel to detect and identify parasites was compared to the standard molecular methods, MSP-1 and MSP-2 typing. Conclusion: This work provides a facile field-deployable genotyping tool that can be used without special skills with standard lab equipment, and at reasonable cost that will unambiguously identify and track P. falciparum parasites both from patient samples and in the laboratory., Engineering and Applied Sciences, Organismic and Evolutionary Biology, Version of Record