Plasmodium falciparum transmission in the highlands of Ethiopia is driven by closely related and clonal parasites.

Malaria cases are frequently recorded in the Ethiopian highlands even at altitudes above 2000 m. The epidemiology of malaria in the Ethiopian highlands, and, in particular, the role of importation by human migration from the highly endemic lowlands is not well understood. We sequenced 187 Plasmodium falciparum samples from two sites in the Ethiopian highlands, Gondar (n = 159) and Ziway (n = 28), using a multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug resistance loci. Here, we characterize the parasite population structure and genetic relatedness. We identify moderate parasite diversity (mean HE : 0.54) and low infection complexity (74.9% monoclonal). A significant percentage of infections share microhaplotypes, even across transmission seasons and sites, indicating persistent local transmission. We identify multiple clusters of clonal or near-clonal infections, highlighting high genetic relatedness. Only 6.3% of individuals diagnosed with P. falciparum reported recent travel. Yet, in clonal or near-clonal clusters, infections of travellers were frequently observed first in time, suggesting that parasites may have been imported and then transmitted locally. 31.1% of infections are pfhrp2-deleted and 84.4% pfhrp3-deleted, and 28.7% have pfhrp2/3 double deletions. Parasites with pfhrp2/3 deletions and wild-type parasites are genetically distinct. Mutations associated with resistance to sulphadoxine-pyrimethamine or suggested to reduce sensitivity to lumefantrine are observed at near-fixation. In conclusion, genomic data corroborate local transmission and the importance of intensified control in the Ethiopian highlands.

Using a mobile nanopore sequencing lab for end-to-end genomic surveillance of Plasmodium falciparum: A feasibility study.

Genomic epidemiology holds promise for malaria control and elimination efforts, for example by informing on Plasmodium falciparum genetic diversity and prevalence of mutations conferring anti-malarial drug resistance. Limited sequencing infrastructure in many malaria-endemic areas prevents the rapid generation of genomic data. To address these issues, we developed and validated assays for P. falciparum nanopore sequencing in endemic sites using a mobile laboratory, targeting key antimalarial drug resistance markers and microhaplotypes. Using two multiplexed PCR reactions, we amplified six highly polymorphic microhaplotypes and ten drug resistance markers. We developed a bioinformatics workflow that allows genotyping of polyclonal malaria infections, including minority clones. We validated the panels on mock dried blood spot (DBS) and rapid diagnostic test (RDT) samples and archived DBS, demonstrating even, high read coverage across amplicons (range: 580x to 3,212x median coverage), high haplotype calling accuracy, and the ability to explore within-sample diversity of polyclonal infections. We field-tested the feasibility of rapid genotyping in Zanzibar in close collaboration with the local malaria elimination program using DBS and routinely collected RDTs as sample inputs. Our assay identified haplotypes known to confer resistance to known antimalarials in the dhfr, dhps and mdr1 genes, but no evidence of artemisinin partial resistance. Most infections (60%) were polyclonal, with high microhaplotype diversity (median HE = 0.94). In conclusion, our assays generated actionable data within a few days, and we identified current challenges for implementing nanopore sequencing in endemic countries to accelerate malaria control and elimination.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania.

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

Genetic variations of Plasmodium falciparum circumsporozoite protein and the impact on interactions with human immunoproteins and malaria vaccine efficacy.

In October 2021, the world's first malaria vaccine RTS,S was endorsed by WHO for broad use in children, despite its low efficacy. This study examined polyclonal infections and the associations of parasite genetic variations with binding affinity to human leukocyte antigen (HLA). Multiplicity of infection was determined by amplicon deep sequencing of PfMSP1. Genetic variations in PfCSP were examined across 88 samples from Ghana and analyzed together with 1655 PfCSP sequences from other African and non-African isolates. Binding interactions of PfCSP peptide variants and HLA were predicted using NetChop and HADDOCK. High polyclonality was detected among infections, with each infection harboring multiple non-3D7 PfCSP variants. Twenty-seven PfCSP haplotypes were detected in the Ghanaian samples, and they broadly represented PfCSP diversity across Africa. The number of genetic differences between 3D7 and non-3D7 PfCSP variants does not influence binding to HLA. However, CSP peptide length after proteolytic degradation significantly affects its molecular weight and binding affinity to HLA. Despite the high diversity of HLA, the majority of the HLAI and II alleles interacted/bound with all Ghana CSP peptides. Multiple non-3D7 strains among P. falciparum infections could impact the effectiveness of RTS,S. Longer peptides of the Th2R/Th3R CSP regions should be considered in future versions of RTS,S.

Projecting vaccine demand and impact for emerging zoonotic pathogens.

BACKGROUND: Despite large outbreaks in humans seeming improbable for a number of zoonotic pathogens, several pose a concern due to their epidemiological characteristics and evolutionary potential. To enable effective responses to these pathogens in the event that they undergo future emergence, the Coalition for Epidemic Preparedness Innovations is advancing the development of vaccines for several pathogens prioritized by the World Health Organization. A major challenge in this pursuit is anticipating demand for a vaccine stockpile to support outbreak response. METHODS: We developed a modeling framework for outbreak response for emerging zoonoses under three reactive vaccination strategies to assess sustainable vaccine manufacturing needs, vaccine stockpile requirements, and the potential impact of the outbreak response. This framework incorporates geographically variable zoonotic spillover rates, human-to-human transmission, and the implementation of reactive vaccination campaigns in response to disease outbreaks. As proof of concept, we applied the framework to four priority pathogens: Lassa virus, Nipah virus, MERS coronavirus, and Rift Valley virus. RESULTS: Annual vaccine regimen requirements for a population-wide strategy ranged from > 670,000 (95% prediction interval 0-3,630,000) regimens for Lassa virus to 1,190,000 (95% PrI 0-8,480,000) regimens for Rift Valley fever virus, while the regimens required for ring vaccination or targeting healthcare workers (HCWs) were several orders of magnitude lower (between 1/25 and 1/700) than those required by a population-wide strategy. For each pathogen and vaccination strategy, reactive vaccination typically prevented fewer than 10% of cases, because of their presently low R0 values. Targeting HCWs had a higher per-regimen impact than population-wide vaccination. CONCLUSIONS: Our framework provides a flexible methodology for estimating vaccine stockpile needs and the geographic distribution of demand under a range of outbreak response scenarios. Uncertainties in our model estimates highlight several knowledge gaps that need to be addressed to target vulnerable populations more accurately. These include surveillance gaps that mask the true geographic distribution of each pathogen, details of key routes of spillover from animal reservoirs to humans, and the role of human-to-human transmission outside of healthcare settings. In addition, our estimates are based on the current epidemiology of each pathogen, but pathogen evolution could alter vaccine stockpile requirements.

Inferring person-to-person networks of Plasmodium falciparum transmission: are analyses of routine surveillance data up to the task?

BACKGROUND: Inference of person-to-person transmission networks using surveillance data is increasingly used to estimate spatiotemporal patterns of pathogen transmission. Several data types can be used to inform transmission network inferences, yet the sensitivity of those inferences to different data types is not routinely evaluated. METHODS: The influence of different combinations of spatial, temporal, and travel-history data on transmission network inferences for Plasmodium falciparum malaria were evaluated. RESULTS: The information content of these data types may be limited for inferring person-to-person transmission networks and may lead to an overestimate of transmission. Only when outbreaks were temporally focal or travel histories were accurate was the algorithm able to accurately estimate the reproduction number under control, Rc. Applying this approach to data from Eswatini indicated that inferences of Rc and spatiotemporal patterns therein depend upon the choice of data types and assumptions about travel-history data. CONCLUSIONS: These results suggest that transmission network inferences made with routine malaria surveillance data should be interpreted with caution.

Impacts of K-12 school reopening on the COVID-19 epidemic in Indiana, USA.

In the United States, schools closed in March 2020 due to COVID-19 and began reopening in August 2020, despite continuing transmission of SARS-CoV-2. In states where in-person instruction resumed at that time, two major unknowns were the capacity at which schools would operate, which depended on the proportion of families opting for remote instruction, and adherence to face-mask requirements in schools, which depended on cooperation from students and enforcement by schools. To determine the impact of these conditions on the statewide burden of COVID-19 in Indiana, we used an agent-based model calibrated to and validated against multiple data types. Using this model, we quantified the burden of COVID-19 on K-12 students, teachers, their families, and the general population under alternative scenarios spanning three levels of school operating capacity (50 %, 75 %, and 100 %) and three levels of face-mask adherence in schools (50 %, 75 %, and 100 %). Under a scenario in which schools operated remotely, we projected 45,579 (95 % CrI: 14,109-132,546) infections and 790 (95 % CrI: 176-1680) deaths statewide between August 24 and December 31. Reopening at 100 % capacity with 50 % face-mask adherence in schools resulted in a proportional increase of 42.9 (95 % CrI: 41.3-44.3) and 9.2 (95 % CrI: 8.9-9.5) times that number of infections and deaths, respectively. In contrast, our results showed that at 50 % capacity with 100 % face-mask adherence, the number of infections and deaths were 22 % (95 % CrI: 16 %-28 %) and 11 % (95 % CrI: 5 %-18 %) higher than the scenario in which schools operated remotely. Within this range of possibilities, we found that high levels of school operating capacity (80-95 %) and intermediate levels of face-mask adherence (40-70 %) resulted in model behavior most consistent with observed data. Together, these results underscore the importance of precautions taken in schools for the benefit of their communities.

Over 100 Years of Rift Valley Fever: A Patchwork of Data on Pathogen Spread and Spillover.

During the past 100 years, Rift Valley fever virus (RVFV), a mosquito-borne virus, has caused potentially lethal disease in livestock, and has been associated with significant economic losses and trade bans. Spillover to humans occurs and can be fatal. Here, we combined data on RVF disease in humans (22 countries) and animals (37 countries) from 1931 to 2020 with seroprevalence studies from 1950 to 2020 (n = 228) from publicly available databases and publications to draw a more complete picture of the past and current RVFV epidemiology. RVFV has spread from its original locus in Kenya throughout Africa and into the Arabian Peninsula. Throughout the study period seroprevalence increased in both humans and animals, suggesting potentially increased RVFV exposure. In 24 countries, animals or humans tested positive for RVFV antibodies even though outbreaks had never been reported there, suggesting RVFV transmission may well go unnoticed. Among ruminants, sheep were the most likely to be exposed during RVF outbreaks, but not during periods of cryptic spread. We discuss critical data gaps and highlight the need for detailed study descriptions, and long-term studies using a one health approach to further convert the patchwork of data to the tale of RFV epidemiology.

Should Deep-Sequenced Amplicons Become the New Gold Standard for Analyzing Malaria Drug Clinical Trials?

Regulatory clinical trials are required to ensure the continued supply and deployment of effective antimalarial drugs. Patient follow-up in such trials typically lasts several weeks, as the drugs have long half-lives and new infections often occur during this period. "Molecular correction" is therefore used to distinguish drug failures from new infections. The current WHO-recommended method for molecular correction uses length-polymorphic alleles at highly diverse loci but is inherently poor at detecting low-density clones in polyclonal infections. This likely leads to substantial underestimates of failure rates, delaying the replacement of failing drugs with potentially lethal consequences. Deep-sequenced amplicons (AmpSeq) substantially increase the detectability of low-density clones and may offer a new "gold standard" for molecular correction. Pharmacological simulation of clinical trials was used to evaluate the suitability of AmpSeq for molecular correction. We investigated the impact of factors such as the number of amplicon loci analyzed, the informatics criteria used to distinguish genotyping "noise" from real low-density signals, the local epidemiology of malaria transmission, and the potential impact of genetic signals from gametocytes. AmpSeq greatly improved molecular correction and provided accurate drug failure rate estimates. The use of 3 to 5 amplicons was sufficient, and simple, nonstatistical criteria could be used to classify recurrent infections as drug failures or new infections. These results suggest AmpSeq is strongly placed to become the new standard for molecular correction in regulatory trials, with potential extension into routine surveillance once the requisite technical support becomes established.

Estimating unobserved SARS-CoV-2 infections in the United States.

By March 2020, COVID-19 led to thousands of deaths and disrupted economic activity worldwide. As a result of narrow case definitions and limited capacity for testing, the number of unobserved severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections during its initial invasion of the United States remains unknown. We developed an approach for estimating the number of unobserved infections based on data that are commonly available shortly after the emergence of a new infectious disease. The logic of our approach is, in essence, that there are bounds on the amount of exponential growth of new infections that can occur during the first few weeks after imported cases start appearing. Applying that logic to data on imported cases and local deaths in the United States through 12 March, we estimated that 108,689 (95% posterior predictive interval [95% PPI]: 1,023 to 14,182,310) infections occurred in the United States by this date. By comparing the model's predictions of symptomatic infections with local cases reported over time, we obtained daily estimates of the proportion of symptomatic infections detected by surveillance. This revealed that detection of symptomatic infections decreased throughout February as exponential growth of infections outpaced increases in testing. Between 24 February and 12 March, we estimated an increase in detection of symptomatic infections, which was strongly correlated (median: 0.98; 95% PPI: 0.66 to 0.98) with increases in testing. These results suggest that testing was a major limiting factor in assessing the extent of SARS-CoV-2 transmission during its initial invasion of the United States.

Amplicon deep sequencing improves Plasmodium falciparum genotyping in clinical trials of antimalarial drugs.

Clinical trials monitoring malaria drug resistance require genotyping of recurrent Plasmodium falciparum parasites to distinguish between treatment failure and new infection occurring during the trial follow up period. Because trial participants usually harbour multi-clonal P. falciparum infections, deep amplicon sequencing (AmpSeq) was employed to improve sensitivity and reliability of minority clone detection. Paired samples from 32 drug trial participants were Illumina deep-sequenced for five molecular markers. Reads were analysed by custom-made software HaplotypR and trial outcomes compared to results from the previous standard genotyping method based on length-polymorphic markers. Diversity of AmpSeq markers in pre-treatment samples was comparable or higher than length-polymorphic markers. AmpSeq was highly reproducible with consistent quantification of co-infecting parasite clones within a host. Outcomes of the three best-performing markers, cpmp, cpp and ama1-D3, agreed in 26/32 (81%) of patients. Discordance between the three markers performed per sample was much lower by AmpSeq (six patients) compared to length-polymorphic markers (eleven patients). Using AmpSeq for discrimination of recrudescence and new infection in antimalarial drug trials provides highly reproducible and robust characterization of clone dynamics during trial follow-up. AmpSeq overcomes limitations inherent to length-polymorphic markers. Regulatory clinical trials of antimalarial drugs will greatly benefit from this unbiased typing method.

Longitudinal tracking and quantification of individual Plasmodium falciparum clones in complex infections.

Longitudinal tracking of individual Plasmodium falciparum strains in multi-clonal infections is essential for investigating infection dynamics of malaria. The traditional genotyping techniques did not permit tracking changes in individual clone density during persistent natural infections. Amplicon deep sequencing (Amp-Seq) offers a tool to address this knowledge gap. The sensitivity of Amp-Seq for relative quantification of clones was investigated using three molecular markers, ama1-D2, ama1-D3, and cpmp. Amp-Seq and length-polymorphism based genotyping were compared for their performance in following minority clones in longitudinal samples from Papua New Guinea. Amp-Seq markers were superior to length-polymorphic marker msp2 in detecting minority clones (sensitivity Amp-Seq: 95%, msp2: 85%). Multiplicity of infection (MOI) by Amp-Seq was 2.32 versus 1.73 for msp2. The higher sensitivity had no effect on estimates of force of infection because missed minority clones were detected in preceding or succeeding bleeds. Individual clone densities were tracked longitudinally by Amp-Seq despite MOI > 1, thus providing an additional parameter for investigating malaria infection dynamics. Amp-Seq based genotyping of longitudinal samples improves detection of minority clones and estimates of MOI. Amp-Seq permits tracking of clone density over time to study clone competition or the dynamics of specific, i.e. resistance-associated genotypes.

Somatic mutation footprinting reveals a unique tetranucleotide signature associated with intron-exon boundaries in lung cancer.

Cigarette smoke comprises a large number of carcinogenic substances that can increase DNA mutation load in epithelial cells of the mouth, throat and lungs. While a strong C:A substitution preference is abundant in tobacco-related cancer genomes, detection of complex or less abundant somatic mutation signatures may be confounded by the heterogeneity of carcinogens present in smoke. Trinucleotide signatures are defined for a variety of somatic mutation processes, yet the extent to which this configuration optimally defines and discriminates between mutational processes is not clear. Here, we describe a method that determines whether trinucleotide patterns do a good job at encapsulating a mutation signature or whether they mask underlying heterogeneity that alternative pattern structures would better define. The approach works by mapping the dependency of trinucleotide signatures in relation to sequence context to establish a 'footprint' of context dependency. Applying this technique to smoke-associated cancers, we show that a robust tetranucleotide substitution is prevalent in 17% of lung squamous cell carcinoma genomes. The signature is dominated by the substitution CT(C:A)G and is strongly associated with gene expression level and intron-exon junctions. Intriguingly, its distribution across the genome is biased towards 5' splice junctions, suggesting a novel mechanism of mutation.

Development of amplicon deep sequencing markers and data analysis pipeline for genotyping multi-clonal malaria infections.

BACKGROUND: Amplicon deep sequencing permits sensitive detection of minority clones and improves discriminatory power for genotyping multi-clone Plasmodium falciparum infections. New amplicon sequencing and data analysis protocols are needed for genotyping in epidemiological studies and drug efficacy trials of P. falciparum. METHODS: Targeted sequencing of molecular marker csp and novel marker cpmp was conducted in duplicate on mixtures of parasite culture strains and 37 field samples. A protocol allowing to multiplex up to 384 samples in a single sequencing run was applied. Software "HaplotypR" was developed for data analysis. RESULTS: Cpmp was highly diverse (He = 0.96) in contrast to csp (He = 0.57). Minority clones were robustly detected if their frequency was >1%. False haplotype calls owing to sequencing errors were observed below that threshold. CONCLUSIONS: To reliably detect haplotypes at very low frequencies, experiments are best performed in duplicate and should aim for coverage of >10'000 reads/amplicon. When compared to length polymorphic marker msp2, highly multiplexed amplicon sequencing displayed greater sensitivity in detecting minority clones.

QuasR: quantification and annotation of short reads in R.

UNLABELLED: QuasR is a package for the integrated analysis of high-throughput sequencing data in R, covering all steps from read preprocessing, alignment and quality control to quantification. QuasR supports different experiment types (including RNA-seq, ChIP-seq and Bis-seq) and analysis variants (e.g. paired-end, stranded, spliced and allele-specific), and is integrated in Bioconductor so that its output can be directly processed for statistical analysis and visualization. AVAILABILITY AND IMPLEMENTATION: QuasR is implemented in R and C/C++. Source code and binaries for major platforms (Linux, OS X and MS Windows) are available from Bioconductor (www.bioconductor.org/packages/release/bioc/html/QuasR.html). The package includes a 'vignette' with step-by-step examples for typical work flows. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DNA sequence explains seemingly disordered methylation levels in partially methylated domains of Mammalian genomes.

For the most part metazoan genomes are highly methylated and harbor only small regions with low or absent methylation. In contrast, partially methylated domains (PMDs), recently discovered in a variety of cell lines and tissues, do not fit this paradigm as they show partial methylation for large portions (20%-40%) of the genome. While in PMDs methylation levels are reduced on average, we found that at single CpG resolution, they show extensive variability along the genome outside of CpG islands and DNase I hypersensitive sites (DHS). Methylation levels range from 0% to 100% in a roughly uniform fashion with only little similarity between neighboring CpGs. A comparison of various PMD-containing methylomes showed that these seemingly disordered states of methylation are strongly conserved across cell types for virtually every PMD. Comparative sequence analysis suggests that DNA sequence is a major determinant of these methylation states. This is further substantiated by a purely sequence based model which can predict 31% (R(2)) of the variation in methylation. The model revealed CpG density as the main driving feature promoting methylation, opposite to what has been shown for CpG islands, followed by various dinucleotides immediately flanking the CpG and a minor contribution from sequence preferences reflecting nucleosome positioning. Taken together we provide a reinterpretation for the nucleotide-specific methylation levels observed in PMDs, demonstrate their conservation across tissues and suggest that they are mainly determined by specific DNA sequence features.

The Ashbya gossypii genome as a tool for mapping the ancient Saccharomyces cerevisiae genome.

We have sequenced and annotated the genome of the filamentous ascomycete Ashbya gossypii. With a size of only 9.2 megabases, encoding 4718 protein-coding genes, it is the smallest genome of a free-living eukaryote yet characterized. More than 90% of A. gossypii genes show both homology and a particular pattern of synteny with Saccharomyces cerevisiae. Analysis of this pattern revealed 300 inversions and translocations that have occurred since divergence of these two species. It also provided compelling evidence that the evolution of S. cerevisiae included a whole genome duplication or fusion of two related species and showed, through inferred ancient gene orders, which of the duplicated genes lost one copy and which retained both copies.

Reinvestigation of the Saccharomyces cerevisiae genome annotation by comparison to the genome of a related fungus: Ashbya gossypii.

BACKGROUND: The recently sequenced genome of the filamentous fungus Ashbya gossypii revealed remarkable similarities to that of the budding yeast Saccharomyces cerevisiae both at the level of homology and synteny (conservation of gene order). Thus, it became possible to reinvestigate the S. cerevisiae genome in the syntenic regions leading to an improved annotation. RESULTS: We have identified 23 novel S. cerevisiae open reading frames (ORFs) as syntenic homologs of A. gossypii genes; for all but one, homologs are present in other eukaryotes including humans. Other comparisons identified 13 overlooked introns and suggested 69 potential sequence corrections resulting in ORF extensions or ORF fusions with improved homology to the syntenic A. gossypii homologs. Of the proposed corrections, 25 were tested and confirmed by resequencing. In addition, homologs of nearly 1,000 S. cerevisiae ORFs, presently annotated as hypothetical, were found in A. gossypii at syntenic positions and can therefore be considered as authentic genes. Finally, we suggest that over 400 S. cerevisiae ORFs that overlap other ORFs in S. cerevisiae and for which no homolog can be detected in A. gossypii should be regarded as spurious. CONCLUSIONS: Although, the S. cerevisiae genome is rightly considered as one of the most accurately sequenced and annotated eukaryotic genomes, we have shown that it still benefits substantially from comparison to the completed sequence and syntenic gene map of A. gossypii, an evolutionarily related fungus. This type of approach will strongly support the annotation of more complex genomes such as the human and murine genomes.