BACKGROUND: Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited. METHODS: A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls. RESULTS: Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes. CONCLUSIONS: The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

While there has been significant progress in controlling falciparum malaria in the Lao Peoples Democratic Republic (PDR), sporadic cases persist in southern provinces where the extent and patterns of transmission remain largely unknown. To assess parasite transmission in this area, 53 Plasmodium falciparum (Pf) positive cases detected through active test and treat campaigns from December 2017 to November 2018 were sequenced, targeting 204 highly polymorphic amplicons. Two R packages, MOIRE and Dcifer, were applied to assess the multiplicity of infections (MOI), effective MOI (eMOI), within-host parasite relatedness, and between-host parasite relatedness ([Formula: see text]). Genomic data were integrated with survey data to characterize the temporal and spatial structures of identified clusters. The positive cases were mainly captured during the focal test and treat campaign conducted in 2018, and in the Pathoomphone area, which had the highest test positivity and forest activity. About 30% of the cases were polyclonal infections, with over half of theses (63%) showing within-host relatedness greater than 0.6, suggesting that cotransmission rather than superinfection was primarily responsible for maintaining polyclonality. A large majority of cases (81%) were infected by parasites genetically linked to one or more other cases. We identified five genetically distinct clusters in forest fringe villages within the Pathoomphone district, characterized by a high degree of genetic relatedness between parasites (mean [Formula: see text] = 0.8). Four smaller clusters of 2-3 cases linked Moonlapamok and Pathoomphone districts, with an average [Formula: see text] of 0.6, suggesting cross-district transmission. Most of the clustered cases occurred within 20 km and 2 months of each other, consistent with focal transmission. Transmission clusters identified in this study confirm the role of ongoing focal parasite transmission occurring within the forest or forest-fringe in the highly mobile population.

Multiple studies have reported a male bias in incidence and/or prevalence of malaria infection in males compared to females. To test the hypothesis that sex-based differences in host-parasite interactions affect the epidemiology of malaria, we intensively followed Plasmodium falciparum infections in a cohort in a malaria endemic area of eastern Uganda and estimated both force of infection (FOI) and rate of clearance using amplicon deep-sequencing. We found no evidence of differences in behavioral risk factors, incidence of malaria, or FOI by sex. In contrast, females cleared asymptomatic infections at a faster rate than males (hazard ratio [HR]=1.82, 95% CI 1.20 to 2.75 by clone and HR = 2.07, 95% CI 1.24 to 3.47 by infection event) in multivariate models adjusted for age, timing of infection onset, and parasite density. These findings implicate biological sex-based differences as an important factor in the host response to this globally important pathogen.

Ixodes scapularis is the principal vector of Lyme disease on the East Coast and in the upper Midwest regions of the United States, yet the tick is also present in the Southeast, where Lyme disease is absent or rare. A closely related species, I. affinis, also carries the pathogen in the South but does not seem to transmit it to humans. In order to better understand the geographic diversity of the tick, we analyzed the microbiota of 104 adult I. scapularis and 13 adult I. affinis ticks captured in 19 locations in South Carolina, North Carolina, Virginia, Connecticut, and New York. Initially, ticks from 4 sites were analyzed by 454 pyrosequencing. Subsequently, ticks from these sites plus 15 others were analyzed by sequencing with an Illumina MiSeq machine. By both analyses, the microbiomes of female ticks were significantly less diverse than those of male ticks. The dissimilarity between tick microbiomes increased with distance between sites, and the state in which a tick was collected could be inferred from its microbiota. The genus Rickettsia was prominent in all locations. Borrelia was also present in most locations and was present at especially high levels in one site in western Virginia. In contrast, members of the family Enterobacteriaceae were very common in North Carolina I. scapularis ticks but uncommon in I. scapularis ticks from other sites and in North Carolina I. affinis ticks. These data suggest substantial variations in the Ixodes microbiota in association with geography, species, and sex.

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.