The spread of drug resistance to antimalarial treatments poses a serious public health risk globally. To combat this risk, molecular surveillance of drug resistance is imperative. We report the prevalence of mutations in the Plasmodium falciparum kelch 13 propeller domain associated with partial artemisinin resistance, which we determined by using Sanger sequencing samples from patients enrolled in therapeutic efficacy studies from 9 sub-Saharan countries during 2014-2018. Of the 2,865 samples successfully sequenced before treatment (day of enrollment) and on the day of treatment failure, 29 (1.0%) samples contained 11 unique nonsynonymous mutations and 83 (2.9%) samples contained 27 unique synonymous mutations. Two samples from Kenya contained the S522C mutation, which has been associated with delayed parasite clearance; however, no samples contained validated or candidate artemisinin-resistance mutations.

Well-studied innate immune systems exist throughout bacteria and archaea, but a more recently discovered genomic locus may offer prokaryotes surprising immunological adaptability. Mediated by a cassette-like genomic locus termed Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), the microbial adaptive immune system differs from its eukaryotic immune analogues by incorporating new immunities unidirectionally. CRISPR thus stores genomically recoverable timelines of virus-host coevolution in natural organisms refractory to laboratory cultivation. Here we combined a population genetic mathematical model of CRISPR-virus coevolution with six years of metagenomic sequencing to link the recoverable genomic dynamics of CRISPR loci to the unknown population dynamics of virus and host in natural communities. Metagenomic reconstructions in an acid-mine drainage system document CRISPR loci conserving ancestral immune elements to the base-pair across thousands of microbial generations. This 'trailer-end conservation' occurs despite rapid viral mutation and despite rapid prokaryotic genomic deletion. The trailer-ends of many reconstructed CRISPR loci are also largely identical across a population. 'Trailer-end clonality' occurs despite predictions of host immunological diversity due to negative frequency dependent selection (kill the winner dynamics). Statistical clustering and model simulations explain this lack of diversity by capturing rapid selective sweeps by highly immune CRISPR lineages. Potentially explaining 'trailer-end conservation,' we record the first example of a viral bloom overwhelming a CRISPR system. The polyclonal viruses bloom even though they share sequences previously targeted by host CRISPR loci. Simulations show how increasing random genomic deletions in CRISPR loci purges immunological controls on long-lived viral sequences, allowing polyclonal viruses to bloom and depressing host fitness. Our results thus link documented patterns of genomic conservation in CRISPR loci to an evolutionary advantage against persistent viruses. By maintaining old immunities, selection may be tuning CRISPR-mediated immunity against viruses reemerging from lysogeny or migration.

The ability to identify mixed-species infections and track the origin of Plasmodium parasites can further enhance the development of treatment and prevention recommendations as well as outbreak investigations. Here, we explore the utility of using the full Plasmodium mitochondrial genome to classify Plasmodium species, detect mixed infections, and infer the geographical origin of imported P. falciparum parasites to the United States (U.S.). Using the recently developed standardized, high-throughput Malaria Resistance Surveillance (MaRS) protocol, the full Plasmodium mitochondrial genomes of 265 malaria cases imported to the U.S. from 2014-2017 were sequenced and analyzed. P. falciparum infections were found in 94.7% (251/265) of samples. Five percent (14/265) of samples were identified as mixed- Plasmodium species or non-P. falciparum, including P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. P. falciparum mitochondrial haplotypes analysis revealed greater than eighteen percent of samples to have at least two P. falciparum mitochondrial genome haplotypes, indicating either heteroplasmy or multi-clonal infections. Maximum-likelihood phylogenies of 912 P. falciparum mitochondrial genomes with known country origin were used to infer the geographical origin of thirteen samples from persons with unknown travel histories as: Africa (country unspecified) (n = 10), Ghana (n = 1), Southeast Asia (n = 1), and the Philippines (n = 1). We demonstrate the utility and current limitations of using the Plasmodium mitochondrial genome to classify samples with mixed-infections and infer the geographical origin of imported P. falciparum malaria cases to the U.S. with unknown travel history.

Anti-malarial resistance is a threat to recent gains in malaria control. This study aimed to assess the efficacy and safety of artesunate-amodiaquine (ASAQ) and artemether-lumefantrine (AL) in the management of uncomplicated malaria and to measure the prevalence of molecular markers of resistance of Plasmodium falciparum in sentinel sites in Maferinyah and Labé Health Districts in Guinea in 2016.

Outbreaks of cyclosporiasis, a food-borne illness caused by the coccidian parasite Cyclospora cayetanensis have increased in the USA in recent years, with approximately 2300 laboratory-confirmed cases reported in 2018. Genotyping tools are needed to inform epidemiological investigations, yet genotyping Cyclospora has proven challenging due to its sexual reproductive cycle which produces complex infections characterized by high genetic heterogeneity. We used targeted amplicon deep sequencing and a recently described ensemble-based distance statistic that accommodates heterogeneous (mixed) genotypes and specimens with partial genotyping data, to genotype and cluster 648 C. cayetanensis samples submitted to CDC in 2018. The performance of the ensemble was assessed by comparing ensemble-identified genetic clusters to analogous clusters identified independently based on common food exposures. Using these epidemiologic clusters as a gold standard, the ensemble facilitated genetic clustering with 93.8% sensitivity and 99.7% specificity. Hence, we anticipate that this procedure will greatly complement epidemiologic investigations of cyclosporiasis.

Biennial therapeutic efficacy monitoring is a crucial activity for ensuring the efficacy of currently used artemisinin-based combination therapy in Angola. Children with acute uncomplicated Plasmodium falciparum infection in sentinel sites in the Benguela, Zaire, and Lunda Sul Provinces were treated with artemether-lumefantrine (AL) or artesunate-amodiaquine (ASAQ) and monitored for 28 days to assess clinical and parasitological responses. Molecular correction was performed using seven microsatellite markers. Samples from treatment failures were genotyped for the pfk13, pfcrt, and pfmdr1 genes. Day 3 clearance rates were ≥95% in all arms. Uncorrected day 28 Kaplan-Meier efficacy estimates ranged from 84.2 to 90.1% for the AL arms and 84.7 to 100% for the ASAQ arms. Corrected day 28 estimates were 87.6% (95% confidence interval [CI], 81 to 95%) for the AL arm in Lunda Sul, 92.2% (95% CI, 87 to 98%) for AL in Zaire, 95.6% (95% CI, 91 to 100%) for ASAQ in Zaire, 98.4% (95% CI, 96 to 100%) for AL in Benguela, and 100% for ASAQ in Benguela and Lunda Sul. All 103 analyzed samples had wild-type pfk13 sequences. The 76T pfcrt allele was found in most (92%; 11/12) ASAQ late-failure samples but in only 16% (4/25) of AL failure samples. The N86 pfmdr1 allele was found in 97% (34/35) of treatment failures. The AL efficacy in Lunda Sul was below the 90% World Health Organization threshold, the third time in four rounds that this threshold was crossed for an AL arm in Angola. In contrast, the observed ASAQ efficacy has not been below 95% to date in Angola, including this latest round.

Routine molecular surveillance for imported drug-resistant malaria parasites to the USA and European Union is an important public health activity. The obtained molecular data are used to help keep chemoprophylaxis and treatment guidelines up to date for persons traveling to malaria endemic countries. Recent advances in next-generation sequencing (NGS) technologies provide a new and effective way of tracking malaria drug-resistant parasites.

Sulfadoxine-pyrimethamine (SP) is used for prevention of malaria in pregnant women in Angola. We sequenced the Plasmodium falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) genes, implicated in SP resistance, in samples collected during a 2019 study of artemisinin-based combination therapy efficacy in Benguela, Lunda Sul, and Zaire provinces. A total of 90 day 0 and day of failure samples were individually sequenced, while 508 day 0 samples from participants without recurrent parasitemia were pooled after DNA extraction into 61 pools. The N51I, C59R, and S108N pfdhfr mutations and A437G pfdhps mutations were present at high proportions in all provinces (weighted allele frequencies, 62% to 100%). The K540E pfdhps mutation was present at lower proportions (10% to 14%). The A581G pfdhps mutation was only observed in Zaire, at a 4.6% estimated prevalence. The I431V and A613S mutations were also only observed in Zaire, at a prevalence of 2.8% to 2.9%. The most common (27% to 66%) reconstructed haplotype in all three provinces was the canonical quadruple pfdhfr pfdhps mutant. The canonical quintuple mutant was absent in Lunda Sul and Benguela and present in 7.9% of samples in Zaire. A single canonical sextuple (2.6%) mutant was observed in Zaire Province. Proportions of the pfdhps K540E and A581G mutations were well below the World Health Organization thresholds for meaningful SP resistance (prevalence of 95% for K540E and 10% for A581G). Samples from therapeutic efficacy studies represent a convenient source of samples for monitoring SP resistance markers.