Genomic epidemiology demonstrates spatially clustered, local transmission of Plasmodium falciparum in forest-going populations in southern Lao PDR.

While there has been significant progress in controlling falciparum malaria in the Lao People's 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.

Sensitive, Highly Multiplexed Sequencing of Microhaplotypes From the Plasmodium falciparum Heterozygome.

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.

Amplicon Sequencing as a Potential Surveillance Tool for Complexity of Infection and Drug Resistance Markers in Plasmodium falciparum Asymptomatic Infections.

BACKGROUND: Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers. METHODS: We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker. RESULTS: We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65-2.77] vs 2.22 [95% CI, 2.17-2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased. CONCLUSIONS: This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Deep Gene Sequence Cluster Analyses of Multi-Virus-Infected Mucosal Tissue Reveal Enhanced Transmission of Acute HIV-1.

Exposure of the genital mucosa to a genetically diverse viral swarm from the donor HIV-1 can result in breakthrough and systemic infection by a single transmitted/founder (TF) virus in the recipient. The highly diverse HIV-1 envelope (Env) in this inoculating viral swarm may have a critical role in transmission and subsequent immune response. Thus, chronic (Envchronic) and acute (Envacute) Env chimeric HIV-1 were tested using multivirus competition assays in human mucosal penile and cervical tissues. Viral competition analysis revealed that Envchronic viruses resided and replicated mainly in the tissue, while Envacute viruses penetrated the human tissue and established infection of CD4+ T cells more efficiently. Analysis of the replication fitness, as tested in peripheral blood mononuclear cells (PBMCs), showed similar replication fitness of Envacute and Envchronic viruses, which did not correlate with transmission fitness in penile tissue. Further, we observed that chimeric Env viruses with higher replication in genital mucosal tissue (chronic Env viruses) had higher binding affinity to C-type lectins. Data presented herein suggest that the inoculating HIV-1 may be sequestered in the genital mucosal tissue (represented by chronic Env HIV-1) but that a single HIV-1 clone (e.g., acute Env HIV-1) can escape this trapped replication for systemic infection.IMPORTANCE During heterosexual HIV-1 transmission, a genetic bottleneck occurs in the newly infected individual as the virus passes from the mucosa, leading to systemic infection with a single transmitted HIV-1 clone in the recipient. This bottleneck in the recipient has just been described (K. Klein et al., PLoS Pathog 14:e1006754, https://doi.org/10.1371/journal.ppat.1006754), and the mechanisms involved in this selection process have not been elucidated. However, understanding mucosal restriction is of the utmost importance for understanding dynamics of infections and for designing focused vaccines. Using our human penile and cervical mucosal tissue models for mixed HIV infections, we provide evidence that HIV-1 from acute/early infection, compared to that from chronic infection, can more efficiently traverse the mucosal epithelium and be transmitted to T cells, suggesting higher transmission fitness. This study focused on the role of the HIV-1 envelope in transmission and provides strong evidence that HIV transmission may involve breaking the mucosal lectin trap.

Efficient Transmission of Mixed Plasmodium falciparum/vivax Infections From Humans to Mosquitoes.

BACKGROUND: In Southeast Asia, people are often coinfected with different species of malaria (Plasmodium falciparum [Pf] and Plasmodium vivax [Pv]) as well as with multiple clones of the same species. Whether particular species or clones within mixed infections are more readily transmitted to mosquitoes remains unknown. METHODS: Laboratory-reared Anopheles dirus were fed on blood from 119 Pf-infected Cambodian adults, with 5950 dissected to evaluate for transmitted infection. Among 12 persons who infected mosquitoes, polymerase chain reaction and amplicon deep sequencing were used to track species and clone-specific transmission to mosquitoes. RESULTS: Seven of 12 persons that infected mosquitoes harbored mixed Pf/Pv infection. Among these 7 persons, all transmitted Pv with 2 transmitting both Pf and Pv, leading to Pf/Pv coinfection in 21% of infected mosquitoes. Up to 4 clones of each species were detected within persons. Shifts in clone frequency were detected during transmission. However, in general, all parasite clones in humans were transmitted to mosquitoes, with individual mosquitoes frequently carrying multiple transmitted clones. CONCLUSIONS: Malaria diversity in human hosts was maintained in the parasite populations recovered from mosquitoes fed on their blood. However, in persons with mixed Pf/Pv malaria, Pv appears to be transmitted more readily, in association with more prevalent patent gametocytemia.

Antimalarial Drug Resistance Profiling of Plasmodium falciparum Infections in Ghana Using Molecular Inversion Probes and Next-Generation Sequencing.

A key drawback to monitoring the emergence and spread of antimalarial drug resistance in sub-Saharan Africa is early detection and containment. Next-generation sequencing methods offer the resolution, sensitivity, and scale required to fill this gap by surveilling for molecular markers of drug resistance. We performed targeted sequencing using molecular inversion probes to interrogate five Plasmodium falciparum genes (pfcrt, pfmdr1, pfdhps, pfdhfr, and pfk13) implicated in chloroquine, sulfadoxine-pyrimethamine (SP), and artemisinin resistance in two sites in Ghana. A total of 803 dried blood spots from children aged between 6 months and 14 years presenting with uncomplicated P. falciparum malaria at the Begoro District Hospital in Begoro and the Ewim Polyclinic in Cape Coast, Ghana, from 2014 to 2017 were prepared on filter paper. Thirteen years after the removal of drug pressure, chloroquine-sensitive parasite strains with pfcrt K76 have increased nearly to fixation in Begoro, in the forest area (prevalence = 95%), but at a lower rate in Cape Coast, in the coastal region (prevalence = 71%, Z = -3.5, P < 0.001). In addition, pfmdr1 184F-bearing parasites are under strong selection. The pfdhfr/pfdhps quadruple genotype ( IRNG K), associated with SP resistance, is near saturation. Our study identified at a 2 to 10% prevalence pfdhps 581G, which is a sulfadoxine resistance marker that correlates with the failure of SP prophylaxis in pregnancy and which has not been observed in Ghana. The differences in the reexpansion of chloroquine-sensitive strains observed at the two study sites, the stronger SP resistance, and the high prevalence of pfmdr1 184F should be further monitored to inform malaria control strategies in Ghana.

SeekDeep: single-base resolution de novo clustering for amplicon deep sequencing.

PCR amplicon deep sequencing continues to transform the investigation of genetic diversity in viral, bacterial, and eukaryotic populations. In eukaryotic populations such as Plasmodium falciparum infections, it is important to discriminate sequences differing by a single nucleotide polymorphism. In bacterial populations, single-base resolution can provide improved resolution towards species and strains. Here, we introduce the SeekDeep suite built around the qluster algorithm, which is capable of accurately building de novo clusters representing true, biological local haplotypes differing by just a single base. It outperforms current software, particularly at low frequencies and at low input read depths, whether resolving single-base differences or traditional OTUs. SeekDeep is open source and works with all major sequencing technologies, making it broadly useful in a wide variety of applications of amplicon deep sequencing to extract accurate and maximal biologic information.

Plasmodium falciparum genetic variation of var2csa in the Democratic Republic of the Congo.

BACKGROUND: The Democratic Republic of the Congo (DRC) bears a high burden of malaria, which is exacerbated in pregnant women. The VAR2CSA protein plays a crucial role in pregnancy-associated malaria (PAM), and hence quantifying diversity at the var2csa locus in the DRC is important in understanding the basic epidemiology of PAM, and in developing a robust vaccine against PAM. METHODS: Samples were taken from the 2013-14 Demographic and Health Survey conducted in the DRC, focusing on children under 5 years of age. A short subregion of the var2csa gene was sequenced in 115 spatial clusters, giving country-wide estimates of sequence polymorphism and spatial population structure. RESULTS: Results indicate that var2csa is highly polymorphic, and that diversity is being maintained through balancing selection, however, there is no clear signal of phylogenetic or geographic structure to this diversity. Linear modelling demonstrates that the number of var2csa variants in a cluster correlates directly with cluster prevalence, but not with other epidemiological factors such as urbanicity. CONCLUSIONS: Results suggest that the DRC fits within the global pattern of high var2csa diversity and little genetic differentiation between regions. A broad multivalent VAR2CSA vaccine candidate could benefit from targeting stable regions and common variants to address the substantial genetic diversity.

A deep sequencing approach to estimate Plasmodium falciparum complexity of infection (COI) and explore apical membrane antigen 1 diversity.

BACKGROUND: Humans living in regions with high falciparum malaria transmission intensity harbour multi-strain infections comprised of several genetically distinct malaria haplotypes. The number of distinct malaria parasite haplotypes identified from an infected human host at a given time is referred to as the complexity of infection (COI). In this study, an amplicon-based deep sequencing method targeting the Plasmodium falciparum apical membrane antigen 1 (pfama1) was utilized to (1) investigate the relationship between P. falciparum prevalence and COI, (2) to explore the population genetic structure of P. falciparum parasites from malaria asymptomatic individuals participating in the 2007 Demographic and Health Survey (DHS) in the Democratic Republic of Congo (DRC), and (3) to explore selection pressures on geospatially divergent parasite populations by comparing AMA1 amino acid frequencies in the DRC and Mali. RESULTS: A total of 900 P. falciparum infections across 11 DRC provinces were examined. Deep sequencing of both individuals, for COI analysis, and pools of individuals, to examine population structure, identified 77 unique pfama1 haplotypes. The majority of individual infections (64.5%) contained polyclonal (COI > 1) malaria infections based on the presence of genetically distinct pfama1 haplotypes. A minimal correlation between COI and malaria prevalence as determined by sensitive real-time PCR was identified. Population genetic analyses revealed extensive haplotype diversity, the vast majority of which was shared across the sites. AMA1 amino acid frequencies were similar between parasite populations in the DRC and Mali. CONCLUSIONS: Amplicon-based deep sequencing is a useful tool for the detection of multi-strain infections that can aid in the understanding of antigen heterogeneity of potential malaria vaccine candidates, population genetics of malaria parasites, and factors that influence complex, polyclonal malaria infections. While AMA1 and other diverse markers under balancing selection may perform well for understanding COI, they may offer little geographic or temporal discrimination between parasite populations.

Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing.

BACKGROUND: Malaria in pregnancy (MiP) remains a major public health challenge in areas of high malaria transmission. Intermittent preventive treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is recommended to prevent the adverse consequences of MiP. The effectiveness of SP for IPTp may be reduced in areas where the dhps581 mutation (a key marker of high level SP resistance) is found; this mutation was previously reported to be common in the Tanga Region of northern Tanzania, but there are limited data from other areas. The frequency of molecular markers of SP resistance was investigated in malaria parasites from febrile patients at health centres (HC) in seven regions comprising the Lake and Southern Zones of mainland Tanzania as part of the ongoing efforts to generate national-wide data of SP resistance. METHODS: A cross-sectional survey was conducted in the outpatient departments of 14 HCs in seven regions from April to June, 2015. 1750 dried blood spot (DBS) samples were collected (117 to 160 per facility) from consenting patients with positive rapid diagnostic tests for malaria, and no recent (within past 2 months) exposure to SP or related drugs. DNA was extracted from the DBS, pooled by HC, and underwent pooled targeted amplicon deep sequencing to yield estimates of mutated parasite allele frequency at each locus of interest. RESULTS: The dhps540 mutation was common across all 14 sites, ranging from 55 to 98.4% of sequences obtained. Frequency of the dhps581 mutation ranged from 0 to 2.4%, except at Kayanga HC (Kagera Region, Lake Zone) where 24.9% of sequences obtained were mutated. The dhfr164 mutation was detected only at Kanyanga HC (0.06%). CONCLUSION: By pooling DNA extracts, the allele frequency of mutations in 14 sites could be directly determined on a single deep-sequencing run. The dhps540 mutant was very common at all locations. Surprisingly, the dhps581 was common at one health center, but rare in all the others, suggesting that there is geographic micro-heterogeneity in mutant distribution and that accurate surveillance requires inclusion of multiple sites. A better understanding of the effect of the dhps581 mutant on the efficacy of IPTp-SP is needed.