Search

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.

BackgroundInference 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.MethodsThe influence of different combinations of spatial, temporal, and travel-history data on transmission network inferences for Plasmodium falciparum malaria were evaluated.ResultsThe 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.ConclusionsThese results suggest that transmission network inferences made with routine malaria surveillance data should be interpreted with caution.

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. [ABSTRACT FROM AUTHOR]

The filamentous fungus Ashbya gossypii is currently used in industry for the production of vitamin [B.sub.2] (1). It is also an attractive model to study filamentous growth, because of its [...], 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.

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. [ABSTRACT FROM AUTHOR]

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% (R2) 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. [ABSTRACT FROM AUTHOR]

Summary: 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. • The world's first malaria vaccine RTS,S has low efficacy against P. falciparum. • Infections were highly polyclonal each with multiple non-3D7 PfCSP variants. • Genetic differences between 3D7 and non-3D7 PfCSP do not affect binding to HLA. • CSP peptide length significantly affects its binding affinity to the HLA. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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., Competing Interests: The authors declare that no competing interests exist., (Copyright: © 2024 Holzschuh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

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.