Search

Cyclospora cayetanensis is an intestinal parasite responsible for the diarrheal illness, cyclosporiasis. Molecular genotyping, using targeted amplicon sequencing, provides a complementary tool for outbreak investigations, especially when epidemiological data are insufficient for linking cases and identifying clusters. The goal of this study was to identify candidate genotyping markers using a novel workflow for detection of segregating single nucleotide polymorphisms (SNPs) in C. cayetanensis genomes. Four whole C. cayetanensis genomes were compared using this workflow and four candidate markers were selected for evaluation of their genotyping utility by PCR and Sanger sequencing. These four markers covered 13 SNPs and resolved parasites from 57 stool specimens, differentiating C. cayetanensis into 19 new unique genotypes.

Abstract Background Because biotechnological uses of bacteriophage gene products as alternatives to conventional antibiotics will require a thorough understanding of their genomic context, we sequenced and analyzed the genomes of four closely related phages isolated from Clostridium perfringens, an important agricultural and human pathogen. Results Phage whole-genome tetra-nucleotide signatures and proteomic tree topologies correlated closely with host phylogeny. Comparisons of our phage genomes to 26 others revealed three shared COGs; of particular interest within this core genome was an endolysin (PF01520, an N-acetylmuramoyl-L-alanine amidase) and a holin (PF04531). Comparative analyses of the evolutionary history and genomic context of these common phage proteins revealed two important results: 1) strongly significant host-specific sequence variation within the endolysin, and 2) a protein domain architecture apparently unique to our phage genomes in which the endolysin is located upstream of its associated holin. Endolysin sequences from our phages were one of two very distinct genotypes distinguished by variability within the putative enzymatically-active domain. The shared or core genome was comprised of genes with multiple sequence types belonging to five pfam families, and genes belonging to 12 pfam families, including the holin genes, which were nearly identical. Conclusions Significant genomic diversity exists even among closely-related bacteriophages. Holins and endolysins represent conserved functions across divergent phage genomes and, as we demonstrate here, endolysins can have significant variability and host-specificity even among closely-related genomes. Endolysins in our phage genomes may be subject to different selective pressures than the rest of the genome. These findings may have important implications for potential biotechnological applications of phage gene products.

Background: In the United States (U.S.), hantavirus pulmonary syndrome (HPS) and non-HPS hantavirus infection are nationally notifiable diseases. Criteria for identifying human cases are based on clinical symptoms (HPS or non-HPS) and acute diagnostic results (IgM+, rising IgG+ titers, RT-PCR+, or immunohistochemistry (IHC)+). Here we provide an overview of diagnostic testing and summarize human Hantavirus disease occurrence and genotype distribution in the U.S. from 2008 to 2020., Methods: Epidemiological data from the national hantavirus registry was merged with laboratory diagnostic testing results performed at the CDC. Residual hantavirus-positive specimens were sequenced, and the available epidemiological and genetic data sets were linked to conduct a genomic epidemiological study of hantavirus disease in the U.S., Findings: From 1993 to 2020, 833 human hantavirus cases have been identified, and from 2008 to 2020, 335 human cases have occurred. Among New World (NW) hantavirus cases detected at the CDC diagnostic laboratory (representing 29.2% of total cases), most (85.0%) were detected during acute disease, however, some convalescent cases were detected in states not traditionally associated with hantavirus infections (Connecticut, Missouri, New Jersey, Pennsylvania, Tennessee, and Vermont). From 1993 to 2020, 94.9% (745/785) of U.S. hantaviruses cases were detected west of the Mississippi with 45.7% (359/785) in the Four Corners region of the U.S. From 2008 to 2020, 67.7% of NW hantavirus cases were detected between the months of March and August. Sequencing of RT-PCR-positive cases demonstrates a geographic separation of Orthohantavirus sinnombreense species [Sin Nombre virus (SNV), New York virus, and Monongahela virus]; however, there is a large gap in viral sequence data from the Northwestern and Central U.S. Finally, these data indicate that commercial IgM assays are not concordant with CDC-developed assays, and that "concordant positive" (i.e., commercial IgM+ and CDC IgM+ results) specimens exhibit clinical characteristics of hantavirus disease., Interpretation: Hantaviral disease is broadly distributed in the contiguous U.S, viral variants are localised to specific geographic regions, and hantaviral disease infrequently detected in most Southeastern states. Discordant results between two diagnostic detection methods highlight the need for an improved standardised testing plan in the U.S. Hantavirus surveillance and detection will continue to improve with clearly defined, systematic reporting methods, as well as explicit guidelines for clinical characterization and diagnostic criteria., Funding: This work was funded by core funds provided to the Viral Special Pathogens Branch at CDC., Competing Interests: The authors declare no competing interests.

Background: Therapeutic efficacy studies (TESs) and detection of molecular markers of drug resistance are recommended by the World Health Organization (WHO) to monitor the efficacy of artemisinin-based combination therapy (ACT). This study assessed the trends of molecular markers of artemisinin resistance and/or reduced susceptibility to lumefantrine using samples collected in TES conducted in Mainland Tanzania from 2016 to 2021., Methods: A total of 2,015 samples were collected during TES of artemether-lumefantrine at eight sentinel sites (in Kigoma, Mbeya, Morogoro, Mtwara, Mwanza, Pwani, Tabora, and Tanga regions) between 2016 and 2021. Photo-induced electron transfer polymerase chain reaction (PET-PCR) was used to confirm presence of malaria parasites before capillary sequencing, which targeted two genes: Plasmodium falciparum kelch 13 propeller domain (k13) and P. falciparum multidrug resistance 1 (pfmdr1)., Results: Sequencing success was ≥ 87.8%, and 1,724/1,769 (97.5%) k13 wild-type samples were detected. Thirty-seven (2.1%) samples had synonymous mutations and only eight (0.4%) had non-synonymous mutations in the k13 gene; seven of these were not validated by the WHO as molecular markers of resistance. One sample from Morogoro in 2020 had a k13 R622I mutation, which is a validated marker of artemisinin partial resistance. For pfmdr1, all except two samples carried N86 (wild-type), while mutations at Y184F increased from 33.9% in 2016 to about 60.5% in 2021, and only four samples (0.2%) had D1246Y mutations. pfmdr1 haplotypes were reported in 1,711 samples, with 985 (57.6%) NYD, 720 (42.1%) NFD, and six (0.4%) carrying minor haplotypes (three with NYY, 0.2%; YFD in two, 0.1%; and NFY in one sample, 0.1%). Between 2016 and 2021, NYD decreased from 66.1% to 45.2%, while NFD increased from 38.5% to 54.7%., Conclusion: This is the first report of the R622I (k13 validated mutation) in Tanzania. N86 and D1246 were nearly fixed, while increases in Y184F mutations and NFD haplotype were observed between 2016 and 2021. Despite the reports of artemisinin partial resistance in Rwanda and Uganda, this study did not report any other validated mutations in these study sites in Tanzania apart from R622I suggesting that intensified surveillance is urgently needed to monitor trends of drug resistance markers and their impact on the performance of ACT., (© 2024. The Author(s).)

Background: The Democratic Republic of the Congo has had 15 Ebola virus disease (EVD) outbreaks, from 1976 to 2023. On June 1, 2020, the Democratic Republic of the Congo declared an outbreak of EVD in the western Équateur Province (11th outbreak), proximal to the 2018 Tumba and Bikoro outbreak and concurrent with an outbreak in the eastern Nord Kivu Province. In this Article, we assessed whether the 11th outbreak was genetically related to previous or concurrent EVD outbreaks and connected available epidemiological and genetic data to identify sources of possible zoonotic spillover, uncover additional unreported cases of nosocomial transmission, and provide a deeper investigation into the 11th outbreak., Methods: We analysed epidemiological factors from the 11th EVD outbreak to identify patient characteristics, epidemiological links, and transmission modes to explore virus spread through space, time, and age groups in the Équateur Province, Democratic Republic of the Congo. Trained field investigators and health professionals recorded data on suspected, probable, and confirmed cases, including demographic characteristics, possible exposures, symptom onset and signs and symptoms, and potentially exposed contacts. We used blood samples from individuals who were live suspected cases and oral swabs from individuals who were deceased to diagnose EVD. We applied whole-genome sequencing of 87 available Ebola virus genomes (from 130 individuals with EVD between May 19 and Sept 16, 2020), phylogenetic divergence versus time, and Bayesian reconstruction of phylogenetic trees to calculate viral substitution rates and study viral evolution. We linked the available epidemiological and genetic datasets to conduct a genomic and epidemiological study of the 11th EVD outbreak., Findings: Between May 19 and Sept 16, 2020, 130 EVD (119 confirmed and 11 probable) cases were reported across 13 Équateur Province health zones. The individual identified as the index case reported frequent consumption of bat meat, suggesting the outbreak started due to zoonotic spillover. Sequencing revealed two circulating Ebola virus variants associated with this outbreak-a Mbandaka variant associated with the majority (97%) of cases and a Tumba-like variant with similarity to the ninth EVD outbreak in 2018. The Tumba-like variant exhibited a reduced substitution rate, suggesting transmission from a previous survivor of EVD., Interpretation: Integrating genetic and epidemiological data allowed for investigative fact-checking and verified patient-reported sources of possible zoonotic spillover. These results demonstrate that rapid genetic sequencing combined with epidemiological data can inform responders of the mechanisms of viral spread, uncover novel transmission modes, and provide a deeper understanding of the outbreak, which is ultimately needed for infection prevention and control during outbreaks., Funding: WHO and US Centers for Disease Control and Prevention., Competing Interests: Declaration of interests We declare no competing interests., (Published by Elsevier Ltd.)

We recently described a targeted amplicon deep sequencing (TADS) strategy that utilizes a nested PCR targeting the 18S rDNA gene of blood-borne parasites. The assay facilitates selective digestion of host DNA by targeting enzyme restriction sites present in vertebrates but absent in parasites. This enriching of parasite-derived amplicon drastically reduces the proportion of host-derived reads during sequencing and results in the sensitive detection of several clinically important blood parasites including Plasmodium spp., Babesia spp., kinetoplastids, and filarial nematodes. Despite these promising results, high costs and the laborious nature of metagenomics sequencing are prohibitive to the routine use of this assay in most laboratories. We describe and evaluate a new metagenomic approach that utilizes a set of primers modified from our original assay that incorporates Illumina barcodes and adapters during the PCR steps. This modification makes amplicons immediately compatible with sequencing on the Illumina MiSeq platform, removing the need for a separate library preparation, which is expensive and time-consuming. We compared this modified assay to our previous nested TADS assay in terms of preparation speed, limit of detection (LOD), and cost. Our modifications reduced assay turnaround times from 7 to 5 days. The cost decreased from approximately $40 per sample to $11 per sample. The modified assay displayed comparable performance in the detection and differentiation of human-infecting Plasmodium spp., Babesia spp., kinetoplastids, and filarial nematodes in clinical samples. The LOD of this modified approach was determined for malaria parasites and remained similar to that previously reported for our earlier assay (0.58 Plasmodium falciparum parasites/µL of blood). These modifications markedly reduced costs and turnaround times, making the assay more amenable to routine diagnostic applications., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Hantaviruses zoonotically infect humans worldwide with pathogenic consequences and are mainly spread by rodents that shed aerosolized virus particles in urine and feces. Bioinformatics methods for hantavirus diagnostics, genomic surveillance and epidemiology are currently lacking a comprehensive approach for data sharing, integration, visualization, analytics and reporting. With the possibility of hantavirus cases going undetected and spreading over international borders, a significant reporting delay can miss linked transmission events and impedes timely, targeted public health interventions. To overcome these challenges, we built HantaNet, a standalone visualization engine for hantavirus genomes that facilitates viral surveillance and classification for early outbreak detection and response. HantaNet is powered by MicrobeTrace, a browser-based multitool originally developed at the Centers for Disease Control and Prevention (CDC) to visualize HIV clusters and transmission networks. HantaNet integrates coding gene sequences and standardized metadata from hantavirus reference genomes into three separate gene modules for dashboard visualization of phylogenetic trees, viral strain clusters for classification, epidemiological networks and spatiotemporal analysis. We used 85 hantavirus reference datasets from GenBank to validate HantaNet as a classification and enhanced visualization tool, and as a public repository to download standardized sequence data and metadata for building analytic datasets. HantaNet is a model on how to deploy MicrobeTrace-specific tools to advance pathogen surveillance, epidemiology and public health globally.

We identified 2 fatal cases of persons infected with hantavirus in Arizona, USA, 2020; 1 person was co-infected with SARS-CoV-2. Delayed identification of the cause of death led to a public health investigation that lasted ≈9 months after their deaths, which complicated the identification of a vector or exposure.

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.

Rift Valley fever (RVF) is a zoonotic disease of public health and economic importance. Uganda has reported sporadic outbreaks of RVF in both humans and animals across the country, especially in the southwestern part of the "cattle corridor" through an established viral hemorrhagic fever surveillance system. We report 52 human cases of laboratory-confirmed RVF from 2017 to 2020. The case fatality rate was 42%. Among those infected, 92% were males and 90% were adults (≥ 18 years). Clinical symptoms were characterized by fever (69%), unexplained bleeding (69%), headache (51%), abdominal pain (49%), and nausea and vomiting (46%). Most of the cases (95%) originated from central and western districts that are part of the cattle corridor of Uganda, where the main risk factor was direct contact with livestock (P = 0.009). Other predictors of RVF positivity were determined to be male gender (P = 0.001) and being a butcher (P = 0.04). Next-generation sequencing identified the predominant Ugandan clade as Kenya-2, observed previously across East Africa. There is need for further investigation and research into the effect and spread of this neglected tropical disease in Uganda and the rest of Africa. Control measures such as promoting vaccination and limiting animal-human transmission could be explored to reduce the impact of RVF in Uganda and globally.

The prevalence of Plasmodium falciparum hrp2 (pfhrp2)-deleted parasites threatens the efficacy of the most used and sensitive malaria rapid diagnostic tests and highlights the need for continued surveillance for this gene deletion. While PCR methods are adequate for determining pfhrp2 presence or absence, they offer a limited view of its genetic diversity. Here, we present a portable sequencing method using the MinION. Pfhrp2 amplicons were generated from individual samples, barcoded, and pooled for sequencing. To overcome potential crosstalk between barcodes, we implemented a coverage-based threshold for pfhrp2 deletion confirmation. Amino acid repeat types were then counted and visualized with custom Python scripts following de novo assembly. We evaluated this assay using well-characterized reference strains and 152 field isolates with and without pfhrp2 deletions, of which 38 were also sequenced on the PacBio platform to provide a standard for comparison. Of 152 field samples, 93 surpassed the positivity threshold, and of those samples, 62/93 had a dominant pfhrp2 repeat type. PacBio-sequenced samples with a dominant repeat-type profile from the MinION sequencing data matched the PacBio profile. This field-deployable assay can be used alone for surveilling pfhrp2 diversity or as a sequencing-based addition to the World Health Organization's existing deletion surveillance protocol., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Naturally occurring human infections by zoonotic Plasmodium species have been documented for P. knowlesi, P. cynomolgi, P. simium, P. simiovale, P. inui, P. inui-like, P. coatneyi, and P. brasilianum. Accurate detection of each species is complicated by their morphological similarities with other Plasmodium species. PCR-based assays offer a solution but require prior knowledge of adequate genomic targets that can distinguish the species. While whole genomes have been published for P. knowlesi, P. cynomolgi, P. simium, and P. inui, no complete genome for P. brasilianum has been available. Previously, we reported a draft genome for P. brasilianum, and here we report the completed genome for P. brasilianum. The genome is 31.4 Mb in size and comprises 14 chromosomes, the mitochondrial genome, the apicoplast genome, and 29 unplaced contigs. The chromosomes consist of 98.4% nucleotide sites that are identical to the P. malariae genome, the closest evolutionarily related species hypothesized to be the same species as P. brasilianum, with 41,125 non-synonymous SNPs (0.0722% of genome) identified between the two genomes. Furthermore, P. brasilianum had 4864 (82.1%) genes that share 80% or higher sequence similarity with 4970 (75.5%) P. malariae genes. This was demonstrated by the nearly identical genomic organization and multiple sequence alignments for the merozoite surface proteins msp3 and msp7. We observed a distinction in the repeat lengths of the circumsporozoite protein (CSP) gene sequences between P. brasilianum and P. malariae. Our results demonstrate a 97.3% pairwise identity between the P. brasilianum and the P. malariae genomes. These findings highlight the phylogenetic proximity of these two species, suggesting that P. malariae and P. brasilianum are strains of the same species, but this could not be fully evaluated with only a single genomic sequence for each species., (© 2022. The Author(s).)

The Plasmodium falciparum protein VAR2CSA allows infected erythrocytes to accumulate within the placenta, inducing pathology and poor birth outcomes. Multiple exposures to placental malaria (PM) induce partial immunity against VAR2CSA, making it a promising vaccine candidate. However, the extent to which VAR2CSA genetic diversity contributes to immune evasion and virulence remains poorly understood. The deep sequencing of the var2csa DBL3X domain in placental blood from forty-nine primigravid and multigravid women living in malaria-endemic western Kenya revealed numerous unique sequences within individuals in association with chronic PM but not gravidity. Additional analysis unveiled four distinct sequence types that were variably present in mixed proportions amongst the study population. An analysis of the abundance of each of these sequence types revealed that one was inversely related to infant gestational age, another was inversely related to placental parasitemia, and a third was associated with chronic PM. The categorization of women according to the type to which their dominant sequence belonged resulted in the segregation of types as a function of gravidity: two types predominated in multigravidae whereas the other two predominated in primigravidae. The univariate logistic regression analysis of sequence type dominance further revealed that gravidity, maternal age, placental parasitemia, and hemozoin burden (within maternal leukocytes), reported a lack of antimalarial drug use, and infant gestational age and birth weight influenced the odds of membership in one or more of these sequence predominance groups. Cumulatively, these results show that unique var2csa sequences differentially appear in women with different PM exposure histories and segregate to types independently associated with maternal factors, infection parameters, and birth outcomes. The association of some var2csa sequence types with indicators of pathogenesis should motivate vaccine efforts to further identify and target VAR2CSA epitopes associated with maternal morbidity and poor birth outcomes.

Sequencing large numbers of individual samples is often needed for countrywide antimalarial drug resistance surveillance. Pooling DNA from several individual samples is an alternative cost and time saving approach for providing allele frequency (AF) estimates at a population level. Using 100 individual patient DNA samples of dried blood spots from a 2017 nationwide drug resistance surveillance study in Haiti, we compared codon coverage of drug resistance-conferring mutations in four Plasmodium falciparum genes (crt, dhps, dhfr, and mdr1), for the same deep sequenced samples run individually and pooled. Samples with similar real-time PCR cycle threshold (Ct) values (+/- 1.0 Ct value) were combined with ten samples per pool. The sequencing success for samples in pools were higher at a lower parasite density than the individual samples sequence method. The median codon coverage for drug resistance-associated mutations in all four genes were greater than 3-fold higher in the pooled samples than in individual samples. The overall codon coverage distribution for pooled samples was wider than the individual samples. The sample pools with < 40 parasites/μL blood showed more discordance in AF calls for dhfr and mdr1 between the individual and pooled samples. This discordance in AF estimation may be due to low amounts of parasite DNA, which could lead to variable PCR amplification efficiencies. Grouping samples with an estimated ≥ 40 parasites/μL blood prior to pooling and deep sequencing yielded the expected population level AF. Pooling DNA samples based on estimates of > 40 parasites/μL prior to deep sequencing can be used for rapid genotyping of a large number of samples for these four genes and possibly other drug resistant markers in population-based studies. As Haiti is a low malaria transmission country with very few mixed infections and continued chloroquine sensitivity, the pooled sequencing approach can be used for routine national molecular surveillance of resistant parasites., Competing Interests: The authors have declared that no competing interests exist.

Background: Partial artemisinin resistance is suspected if delayed parasite clearance (ie, persistence of parasitaemia on day 3 after treatment initiation) is observed. Validated markers of artemisinin partial resistance in southeast Asia, Plasmodium falciparum kelch13 (Pfkelch13) R561H and P574L, have been reported in Rwanda but no association with parasite clearance has been observed. We aimed to establish the efficacy of artemether-lumefantrine and genetic characterisation of Pfkelch13 alleles and their association with treatment outcomes., Methods: This open-label, single-arm, multicentre, therapeutic efficacy study was done in 2018 in three Rwandan sites: Masaka, Rukara, and Bugarama. Children aged 6-59 months with P falciparum monoinfection and fever were eligible and treated with a 3-day course of artemether-lumefantrine. Treatment response was monitored for 28 days using weekly microscopy screenings of blood samples for P falciparum. Mutations in Pfkelch13 and P falciparum multidrug resistance-1 (Pfmdr1) genes were characterised in parasites collected from enrolled participants. Analysis of flanking microsatellites surrounding Pfkelch13 was done to define the origins of the R561H mutations. The primary endpoint was PCR-corrected parasitological cure on day 28, as per WHO protocol., Findings: 228 participants were enrolled and 224 (98·2%) reached the study endpoint. PCR-corrected efficacies were 97·0% (95% CI 88-100) in Masaka, 93·8% (85-98) in Rukara, and 97·2% (91-100) in Bugarama. Pfkelch13 R561H mutations were present in 28 (13%) of 218 pre-treatment samples and P574L mutations were present in two (1%) pre-treatment samples. 217 (90%) of the 240 Pfmdr1 haplotypes observed in the pretreatment samples, had either the NFD (N86Y, Y184F, D1246Y) or NYD haplotype. Eight (16%) of 51 participants in Masaka and 12 (15%) of 82 participants in Rukara were microscopically positive 3 days after treatment initiation, which was associated with pre-treatment presence of Pfkelch13 R561H in Masaka (p=0·0005). Genetic analysis of Pfkelch13 R561H mutations suggest their common ancestry and local origin in Rwanda., Interpretation: We confirm evidence of emerging artemisinin partial resistance in Rwanda. Although artemether-lumefantrine remains efficacious, vigilance for decreasing efficacy, further characterisation of artemisinin partial resistance, and evaluation of additional antimalarials in Rwanda should be considered., Funding: The US President's Malaria Initiative., Translation: For the French translation of the abstract see Supplementary Materials section., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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.

Background: Recent studies showed the first emergence of the R561H artemisinin-associated resistance marker in Africa, which highlights the importance of continued molecular surveillance to assess the selection and spread of this and other drug resistance markers in the region., Method: In this study, we used targeted amplicon deep sequencing of 116 isolates collected in two areas of Cameroon to genotype the major drug resistance genes, k13, crt, mdr1, dhfr, and dhps, and the cytochrome b gene (cytb) in Plasmodium falciparum., Results: No confirmed or associated artemisinin resistance markers were observed in Pfk13. In comparison, both major and minor alleles associated with drug resistance were found in Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps. Notably, a high frequency of other nonsynonymous mutations was observed across all the genes, except for Pfcytb, suggesting continued selection pressure., Conclusions: The results from this study supported the continued use of artemisinin-based combination therapy and administration of sulfadoxine-pyrimethamine for intermittent preventive therapy in pregnant women, and for seasonal chemoprevention in these study sites in Cameroon., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: The current first-line treatments for uncomplicated malaria recommended by the National Malaria Control Programme in Mali are artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ). From 2015 to 2016, an in vivo study was carried out to assess the clinical and parasitological responses to AL and ASAQ in Sélingué, Mali., Methods: Children between 6 and 59 months of age with uncomplicated Plasmodium falciparum infection and 2000-200,000 asexual parasites/μL of blood were enrolled, randomly assigned to either AL or ASAQ, and followed up for 42 days. Uncorrected and PCR-corrected efficacy results at days 28 and 42. were calculated. Known markers of resistance in the Pfk13, Pfmdr1, and Pfcrt genes were assessed using Sanger sequencing., Results: A total of 449 patients were enrolled: 225 in the AL group and 224 in the ASAQ group. Uncorrected efficacy at day 28 was 83.4% (95% CI 78.5-88.4%) in the AL arm and 93.1% (95% CI 89.7-96.5%) in the ASAQ arm. The per protocol PCR-corrected efficacy at day 28 was 91.0% (86.0-95.9%) in the AL arm and 97.1% (93.6-100%) in the ASAQ arm. ASAQ was significantly (p < 0.05) better than AL for each of the aforementioned efficacy outcomes. No mutations associated with artemisinin resistance were identified in the Pfk13 gene. Overall, for Pfmdr1, the N86 allele and the NFD haplotype were the most common. The NFD haplotype was significantly more prevalent in the post-treatment than in the pre-treatment isolates in the AL arm (p < 0.01) but not in the ASAQ arm. For Pfcrt, the CVIET haplotype was the most common., Conclusions: The findings indicate that both AL and ASAQ remain effective for the treatment of uncomplicated malaria in Sélingué, Mali.

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., (Copyright © 2021 Dimbu et al.)

Background: Targeted amplicon deep sequencing (TADS) has enabled characterization of diverse bacterial communities, yet the application of TADS to communities of parasites has been relatively slow to advance. The greatest obstacle to this has been the genetic diversity of parasitic agents, which include helminths, protozoa, arthropods, and some acanthocephalans. Meanwhile, universal amplification of conserved loci from all parasites without amplifying host DNA has proven challenging. Pan-eukaryotic PCRs preferentially amplify the more abundant host DNA, obscuring parasite-derived reads following TADS. Flaherty et al. (2018) described a pan-parasitic TADS method involving amplification of eukaryotic 18S rDNA regions possessing restriction sites only in vertebrates. Using this method, host DNA in total DNA extracts could be selectively digested prior to PCR using restriction enzymes, thereby increasing the number of parasite-derived reads obtained following NGS. This approach showed promise though was only as sensitive as conventional PCR., Results: Here, we expand on this work by designing a second set of pan-eukaryotic primers flanking the priming sites already described, enabling nested PCR amplification of the established 18S rDNA target. This nested approach facilitated introduction of a second restriction digestion between the first and second PCR, reducing the proportional mass of amplifiable host-derived DNA while increasing the number of PCR amplification cycles. We applied this method to blood specimens containing Babesia, Plasmodium, various kinetoplastids, and filarial nematodes and confirmed its limit of detection (LOD) to be approximately 10-fold lower than previously described, falling within the range of most qPCR methods., Conclusions: The assay detects and differentiates the major malaria parasites of humans, along with several other clinically important blood parasites. This represents an important step towards a TADS-based universal parasite diagnostic (UPDx) test with a sufficient LOD for routine applications. Video Abstract.

The number of Asian migrants working in sub-Saharan developing countries like Angola has been increasing. Their malaria risk, prevention, and care-seeking practices have not been characterized. A cross-sectional survey was conducted in 733 Chinese and Southeast Asian migrants in Angola. Respondents were interviewed and provided blood samples. Samples were analyzed to detect Plasmodium antigen and characterize host anti- Plasmodium response. Positive samples were genotyped using the pfs47 marker. Most respondents (72%; 95% CI: 68-75) reported using bed nets, but less than 1% reported using chemoprophylaxis. Depending on the assay, 1-4% of respondents had evidence of active malaria infection. By contrast, 55% (95% CI: 52-59) were seropositive for Plasmodium antibodies. Most infections were Plasmodium falciparum , but infection and/or exposure to Plasmodium vivax and Plasmodium malariae was also detected. Seroprevalence by time in Angola showed most exposure occurred locally. One respondent had sufficiently high parasitemia for pfs47 genotyping, which showed that the infection was likely locally acquired despite recent travel to home country. Asian migrants to Angola are at substantial risk of malaria. Employers should consider enhanced malaria prevention programs, including chemoprophylaxis; embassies should encourage prevention practices. Angolan healthcare workers should be aware of high malaria exposure in Asian migrants.

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.

Histidine-rich protein 2 (HRP2) detecting rapid diagnostic tests (RDTs) have played an important role in enabling prompt malaria diagnosis in remote locations. However, emergence of pfhrp2 deleted parasites is threatening the efficacy of RDTs, and the World Health Organization (WHO) has highlighted surveillance of these deletions as a priority. Nested PCR is used to confirm pfhrp2 deletion but is costly and laborious. Due to spurious amplification of paralogue pfhrp3, the identity of nested exon 1 PCR product must be confirmed by sequencing. Here we describe a new one-step PCR method for detection of pfhrp2. To determine sensitivity and specificity, all PCRs were performed in triplicate. Using photo-induced electron transfer (PET) PCR detecting 18srRNA as true positive, one-step had comparable sensitivity of 95.0% (88.7-98.4%) to nested exon 1, 99.0% (94.6-99.9%) and nested exon 2, 98.0% (93.0-99.8%), and comparable specificity 93.8% (69.8-99.8%) to nested exon 1 100.0% (79.4-100.0%) and nested exon 2, 100.0% (74.4-100.0%). Sequencing revealed that one step PCR does not amplify pfhrp3. Logistic regression models applied to measure the 95% level of detection of the one-step PCR in clinical isolates provided estimates of 133p/μL (95% confidence interval (CI): 3-793p/μL) for whole blood (WB) samples and 385p/μL (95% CI: 31-2133 p/μL) for dried blood spots (DBSs). When considering protocol attributes, the one-step PCR is less expensive, faster and more suitable for high throughput. In summary, we have developed a more accurate PCR method that may be ideal for the application of the WHO protocol for investigating pfhrp2 deletions in symptomatic individuals presenting to health care facilities., Competing Interests: Authors Gireesh Subramaniam and Jasmine Padilla are employed by the Oak Ridge Institute for Science and Education (ORISE). ORISE provided support in the form of salary for authors GS and JP, but did not have any additional role role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests statement: Authors Sophie Jones and Dhruviben Patel are employed by Williams Consulting LLC. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials. Authors Gireesh Subramaniam and Jasmine Padilla are employed by the Oak Ridge Institute for Science and Education (ORISE). There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Background: 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., Methods: This was a two-arm randomised controlled trial of the efficacy of AL and ASAQ among children aged 6-59 months with uncomplicated Plasmodium falciparum malaria in two sites. Children were followed for 28 days to assess clinical and parasitological response. The primary outcome was the Kaplan-Meier estimate of Day 28 (D28) efficacy after correction by microsatellite-genotyping. Pre-treatment (D0) and day of failure samples were assayed for molecular markers of resistance in the pfk13 and pfmdr1 genes., Results: A total of 421 participants were included with 211 participants in the Maferinyah site and 210 in Labé. No early treatment failure was observed in any study arms. However, 22 (5.3%) participants developed a late treatment failure (8 in the ASAQ arm and 14 in the AL arm), which were further classified as 2 recrudescences and 20 reinfections. The Kaplan-Meier estimate of the corrected efficacy at D28 was 100% for both AL and ASAQ in Maferinyah site and 99% (95% Confidence Interval: 97.2-100%) for ASAQ and 99% (97.1-100%) for AL in Labé. The majority of successfully analysed D0 (98%, 380/389) and all day of failure (100%, 22/22) samples were wild type for pfk13. All 9 observed pfk13 mutations were polymorphisms not associated with artemisinin resistance. The NFD haplotype was the predominant haplotype in both D0 (197/362, 54%) and day of failure samples (11/18, 61%) successfully analysed for pfmdr1., Conclusion: This study observed high efficacy and safety of both ASAQ and AL in Guinea, providing evidence for their continued use to treat uncomplicated malaria. Continued monitoring of ACT efficacy and safety and molecular makers of resistance in Guinea is important to detect emergence of parasite resistance and to inform evidence-based malaria treatment policies.

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics., (Published by Oxford University Press for the Infectious Diseases Society of America 2019.)

Background: 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., Methods: As part of a technology transfer arrangement between the CDC Malaria Branch and the Istituto Superiore di Sanità (ISS), Rome, Italy, the recently described Malaria Resistance Surveillance (MaRS) protocol was used to genotype 148 Plasmodium falciparum isolates from Eritrea for kelch 13 (k13) and cytochrome b (cytb) genes, molecular markers associated with resistance to artemisinin (ART) and atovaquone/proguanil (AP), respectively., Results: Spanning the full-length k13 gene, seven non-synonymous single nucleotide polymorphisms (SNPs) were found (K189N, K189T, E208K, D281V, E401Q, R622I and T535M), of which none have been associated with artemisinin resistance. No mutations were found in cytochrome b., Conclusion: All patients successfully genotyped carried parasites susceptible to ART and AP treatment. Future studies between CDC Malaria Branch and ISS are planned to expand the MaRS system, including data sharing, in an effort to maintain up to date treatment guidelines for travelers to malaria endemic countries.

Strongyloidiasis is a neglected tropical disease caused by the human infective nematodes Strongyloides stercoralis, Strongyloides fuelleborni fuelleborni and Strongyloides fuelleborni kellyi. Previous large-scale studies exploring the genetic diversity of this important genus have focused on Southeast Asia, with a small number of isolates from the USA, Switzerland, Australia and several African countries having been genotyped. Consequently, little is known about the global distribution of geographic sub-variants of these nematodes and the genetic diversity that exists within the genus Strongyloides generally. We extracted DNA from human, dog and primate feces containing Strongyloides, collected from several countries representing all inhabited continents. Using a genotyping assay adapted for deep amplicon sequencing on the Illumina MiSeq platform, we sequenced the hyper-variable I and hyper-variable IV regions of the Strongyloides 18S rRNA gene and a fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1) gene from these specimens. We report several novel findings including unique S. stercoralis and S. fuelleborni genotypes, and the first identifications of a previously unknown S. fuelleborni infecting humans within Australia. We expand on an existing Strongyloides genotyping scheme to accommodate S. fuelleborni and these novel genotypes. In doing so, we compare our data to all 18S and cox1 sequences of S. fuelleborni and S. stercoralis available in GenBank (to our knowledge), that overlap with the sequences generated using our approach. As this analysis represents more than 1,000 sequences collected from diverse hosts and locations, representing all inhabited continents, it allows a truly global understanding of the population genetic structure of the Strongyloides species infecting humans, non-human primates, and domestic dogs., Competing Interests: RB is a co-author on the following International Patent: #PCT/US2018/052469. Talundzic E, Bradbury R, Olsen C, Flaherty B. Removing Interfering Host Nucleic Acids for Molecular Parasite Detection. Meredith Lane is employed by Synergy America as a contractor to provide laboratory technical service to the Division of Parasitic Diseases and Malaria at the Centers for Disease Control and Prevention. In this capacity, there are no competing Interests in Ms. Lane's employment at Synergy and her authorship of this paper. Meredith Lane is correctly affiliated to this company (Synergy America). Meredith Lane and Synergy America hold no consultancies or patents and have no products in development or marketed products that are competing interests with this publication.

Sexually reproducing pathogens such as Cyclospora cayetanensis often produce genetically heterogeneous infections where the number of unique sequence types detected at any given locus varies depending on which locus is sequenced. The genotypes assigned to these infections quickly become complex when additional loci are analysed. This genetic heterogeneity confounds the utility of traditional sequence-typing and phylogenetic approaches for aiding epidemiological trace-back, and requires new methods to address this complexity. Here, we describe an ensemble of two similarity-based classification algorithms, including a Bayesian and heuristic component that infer the relatedness of C. cayetanensis infections. The ensemble requires a set of haplotypes as input and assigns arbitrary distances to specimen pairs reflecting their most likely relationships. The approach was applied to data generated from a test cohort of 88 human fecal specimens containing C. cayetanensis, including 30 from patients whose infections were associated with epidemiologically defined outbreak clusters of cyclosporiasis. The ensemble assigned specimens to plausible clusters of genetically related infections despite their complex haplotype composition. These relationships were corroborated by a significant number of epidemiological linkages (P < 0.0001) suggesting the ensemble's utility for aiding epidemiological trace-back investigations of cyclosporiasis.

Mutations in the Plasmodium falciparum k13 ( Pfk13 ) gene are linked to delayed parasite clearance in response to artemisinin-based combination therapies (ACTs) in Southeast Asia. To explore the evolutionary rate and constraints acting on this gene, k13 orthologs from species sharing a recent common ancestor with P. falciparum and Plasmodium vivax were analyzed. These comparative studies were followed by genetic polymorphism analyses within P. falciparum using 982 complete Pfk13 sequences from public databases and new data obtained by next-generation sequencing from African and Haitian isolates. Although k13 orthologs evolve at heterogeneous rates, the gene was conserved across the genus, with only synonymous substitutions being found at residues where mutations linked to the delayed parasite clearance phenotype have been reported. This suggests that those residues were under constraint from undergoing nonsynonymous changes during evolution of the genus. No fixed nonsynonymous differences were found between Pfk13 and its orthologs in closely related species found in African apes. This indicates that all nonsynonymous substitutions currently found in Pfk13 are younger than the time of divergence between P. falciparum and its closely related species. At the population level, no mutations linked to delayed parasite clearance were found in our samples from Africa and Haiti. However, there is a high number of single Pfk13 mutations segregating in P. falciparum populations, and two predominant alleles are distributed worldwide. This pattern is discussed in terms of how changes in the efficacy of natural selection, affected by population expansion, may have allowed for the emergence of mutations tolerant to ACTs., (Copyright © 2019 American Society for Microbiology.)

Cyclosporiasis is an infection caused by Cyclospora cayetanensis, which is acquired by consumption of contaminated fresh food or water. In the United States, cases of cyclosporiasis are often associated with foodborne outbreaks linked to imported fresh produce or travel to disease-endemic countries. Epidemiologic investigation has been the primary method for linking outbreak cases. A molecular typing marker that can identify genetically related samples would be helpful in tracking outbreaks. We evaluated the mitochondrial junction region as a potential genotyping marker. We tested stool samples from 134 laboratory-confirmed cases in the United States by using PCR and Sanger sequencing. All but 2 samples were successfully typed and divided into 14 sequence types. Typing results were identical among samples within each epidemiologically defined case cluster for 7 of 10 clusters. These findings suggest that this marker can distinguish between distinct case clusters and might be helpful during cyclosporiasis outbreak investigations.

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., Competing Interests: Author Dhruviben Patel is employed by Williams Consulting LLC. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Background: The World Health Organization recommends regular therapeutic efficacy studies (TES) to monitor the performance of first and second-line anti-malarials. In 2016, efficacy and safety of artemether-lumefantrine (AL) for the treatment of uncomplicated falciparum malaria were assessed through a TES conducted between April and October 2016 at four sentinel sites of Kibaha, Mkuzi, Mlimba, and Ujiji in Tanzania. The study also assessed molecular markers of artemisinin and lumefantrine (partner drug) resistance., Methods: Eligible patients were enrolled at the four sites, treated with standard doses of AL, and monitored for 28 days with clinical and laboratory assessments. The main outcomes were PCR corrected cure rates, day 3 positivity rates, safety of AL, and prevalence of single nucleotide polymorphisms in Plasmodium falciparum kelch 13 (Pfk13) (codon positions: 440-600) and P. falciparum multi-drug resistance 1 (Pfmdr1) genes (codons: N86Y, Y184F and D1246Y), markers of artemisinin and lumefantrine resistance, respectively., Results: Of 344 patients enrolled, three withdrew, six were lost to follow-up; and results were analysed for 335 (97.4%) patients. Two patients had treatment failure (one early treatment failure and one recrudescent infection) after PCR correction, yielding an adequate clinical and parasitological response of > 98%. Day 3 positivity rates ranged from 0 to 5.7%. Common adverse events included cough, abdominal pain, vomiting, and diarrhoea. Two patients had serious adverse events; one died after the first dose of AL and another required hospitalization after the second dose of AL (on day 0) but recovered completely. Of 344 samples collected at enrolment (day 0), 92.7% and 100% were successfully sequenced for Pfk13 and Pfmdr1 genes, respectively. Six (1.9%) had non-synonymous mutations in Pfk13, none of which had been previously associated with artemisinin resistance. For Pfmdr1, the NFD haplotype (codons N86, 184F and D1246) was detected in 134 (39.0%) samples; ranging from 33.0% in Mlimba to 45.5% at Mkuzi. The difference among the four sites was not significant (p = 0.578). All samples had a single copy of the Pfmdr1 gene., Conclusion: The study indicated high efficacy of AL and the safety profile was consistent with previous reports. There were no known artemisinin-resistance Pfk13 mutations, but there was a high prevalence of a Pfmdr1 haplotype associated with reduced sensitivity to lumefantrine (but no reduced efficacy was observed in the subjects). Continued TES and monitoring of markers of resistance to artemisinin and partner drugs is critical for early detection of resistant parasites and to inform evidence-based malaria treatment policies. Trial Registration ClinicalTrials.gov NCT03387631.

Background: Cyclospora cayetanensis is a food-borne intestinal human parasite that causes outbreaks of diarrhea. There is a need for efficient laboratory methods for strain-level characterization to assist in outbreak investigations. By using next generation sequencing, genomic sequences can be obtained and compared to identify potential genotyping markers. However, there is no method available to propagate this parasite in the laboratory. Therefore, genomic DNA must be extracted from oocysts purified from human stool. The objective of this study was to apply optimized methods to purify C. cayetanensis oocysts and extract DNA in order to obtain high-quality whole genome sequences with minimum contamination of DNA from other organisms., Results: Oocysts from 21 human stool specimens were separated from other stool components using discontinuous density gradient centrifugation and purified further by flow cytometry. Genomic DNA was used to construct Ovation Ultralow libraries for Illumina sequencing. MiSeq sequencing reads were taxonomically profiled for contamination, de novo assembled, and mapped to a draft genome available in GenBank to assess the quality of the resulting genomic sequences. Following all purification steps, the majority (81-99%) of sequencing reads were from C. cayetanensis . They could be assembled into draft genomes of around 45 MB in length with GC-content of 52%., Conclusions: Density gradients performed in the presence of a detergent followed by flow cytometry sorting of oocysts yielded sufficient genomic DNA largely free from contamination and suitable for whole genome sequencing of C. cayetanensis. The methods described here will facilitate the accumulation of genomic sequences from various samples, which is a prerequisite for the development of typing tools to aid in outbreak investigations.

Background: Targeted amplicon deep sequencing (TADS) of the 16S rRNA gene is commonly used to explore and characterize bacterial microbiomes. Meanwhile, attempts to apply TADS to the detection and characterization of entire parasitic communities have been hampered since conserved regions of many conserved parasite genes, such as the 18S rRNA gene, are also conserved in their eukaryotic hosts. As a result, targeted amplification of 18S rRNA from clinical samples using universal primers frequently results in competitive priming and preferential amplification of host DNA. Here, we describe a novel method that employs a single pair of universal primers to capture all blood-borne parasites while reducing host 18S rRNA template and enhancing the amplification of parasite 18S rRNA for TADS. This was achieved using restriction enzymes to digest the 18S rRNA gene at cut sites present only in the host sequence prior to PCR amplification., Results: This method was validated against 16 species of blood-borne helminths and protozoa. Enzyme digestion prior to PCR enrichment and Illumina amplicon deep sequencing led to a substantial reduction in human reads and a corresponding 5- to 10-fold increase in parasite reads relative to undigested samples. This method allowed for discrimination of all common parasitic agents found in human blood, even in cases of multi-parasite infection, and markedly reduced the limit of detection in digested versus undigested samples., Conclusions: The results herein provide a novel methodology for the reduction of host DNA prior to TADS and establish the validity of a next-generation sequencing-based platform for universal parasite detection.

Background: The Angolan government recommends three artemisinin-based combinations for the treatment of uncomplicated Plasmodium falciparum malaria: artemether-lumefantrine (AL), artesunate-amodiaquine (ASAQ), and dihydroartemisinin-piperaquine (DP). Due to the threat of emerging anti-malarial drug resistance, it is important to periodically monitor the efficacy of artemisinin-based combination therapy (ACT). This study evaluated these medications' therapeutic efficacy in Benguela, Lunda Sul, and Zaire Provinces., Methods: Enrollment occurred between March and July 2017. Study participants were children with P. falciparum monoinfection from each provincial capital. Participants received a 3-day course of a quality-assured artemisinin-based combination and were monitored for 28 (AL and ASAQ arms) or 42 days (DP arm). Each ACT was assessed in two provinces. The primary study endpoints were: (1) follow-up without complications and (2) failure to respond to treatment or development of recurrent P. falciparum infection. Parasites from each patient experiencing recurrent infection were genotyped to differentiate new infection from recrudescence of persistent parasitaemia. These parasites were also analysed for molecular markers associated with ACT resistance., Results: Of 608 children enrolled in the study, 540 (89%) reached a primary study endpoint. Parasitaemia was cleared within 3 days of medication administration in all participants, and no early treatment failures were observed. After exclusion of reinfections, the corrected efficacy of AL was 96% (91-100%, 95% confidence interval) in Zaire and 97% (93-100%) in Lunda Sul. The corrected efficacy of ASAQ was 100% (97-100%) in Benguela and 93% (88-99%) in Zaire. The corrected efficacy of DP was 100% (96-100%) in Benguela and 100% in Lunda Sul. No mutations associated with artemisinin resistance were identified in the pfk13 gene in the 38 cases of recurrent P. falciparum infection. All 33 treatment failures in the AL and ASAQ arms carried pfmdr1 or pfcrt mutations associated with lumefantrine and amodiaquine resistance, respectively, on day of failure., Conclusions: AL, ASAQ, and DP continue to be efficacious against P. falciparum malaria in these provinces of Angola. Rapid parasite clearance and the absence of genetic evidence of artemisinin resistance are consistent with full susceptibility to artemisinin derivatives. Periodic monitoring of in vivo drug efficacy remains a priority routine activity for Angola.

The recent advances in next-generation sequencing technologies provide a new and effective way of tracking malaria drug-resistant parasites. To take advantage of this technology, an end-to-end Illumina targeted amplicon deep sequencing (TADS) and bioinformatics pipeline for molecular surveillance of drug resistance in P. falciparum , called ma laria r esistance s urveillance (MaRS), was developed. TADS relies on PCR enriching genomic regions, specifically target genes of interest, prior to deep sequencing. MaRS enables researchers to simultaneously collect data on allele frequencies of multiple full-length P. falciparum drug resistance genes ( crt , mdr1 , k13 , dhfr , dhps , and the cytochrome b gene), as well as the mitochondrial genome. Information is captured at the individual patient level for both known and potential new single nucleotide polymorphisms associated with drug resistance. The MaRS pipeline was validated using 245 imported malaria cases that were reported to the Centers for Disease Control and Prevention (CDC). The chloroquine resistance crt CV IET genotype (mutations underlined) was observed in 42% of samples, the highly pyrimethamine-resistant dhpsIRN triple mutant in 92% of samples, and the sulfadoxine resistance dhps mutation S GE AA in 26% of samples. The mdr1 N F SND genotype was found in 40% of samples. With the exception of two cases imported from Cambodia, no artemisinin resistance k13 alleles were identified, and 99% of patients carried parasites susceptible to atovaquone-proguanil. Our goal is to implement MaRS at the CDC for routine surveillance of imported malaria cases in the United States and to aid in the adoption of this system at participating state public health laboratories, as well as by global partners., (Copyright © 2018 American Society for Microbiology.)

Artemisinin-based combination therapy (ACT) is the most effective and widely used treatment for uncomplicated Plasmodium falciparum malaria and is a cornerstone for malaria control and prevention globally. Resistance to artemisinin derivatives has been confirmed in the Greater Mekong Subregion (GMS) and manifests as slow parasite clearance in patients and reduced ring stage susceptibility to artemisinins in survival assays. The P. falciparum kelch13 gene mutations associated with artemisinin-resistant parasites are now widespread in the GMS. We genotyped 277 samples collected during an observational study from 2012 to 2016 from eight provinces in Thailand to identify P. falciparum kelch13 mutations. The results were combined with previously reported genotyping results from Thailand to construct a map illustrating the evolution of P. falciparum kelch13 mutations from 2007 to 2016 in that country. Different mutant alleles were found in strains with different geographical origins. The artemisinin resistance-conferring Y493H and R539T mutations were detected mainly in eastern Thailand (bordering Cambodia), while P574L was found only in western Thailand and R561H only in northwestern Thailand. The C580Y mutation was found across the entire country and was nearing fixation along the Thai-Cambodia border. Overall, the prevalence of artemisinin resistance mutations increased over the last 10 years across Thailand, especially along the Thai-Cambodia border. Molecular surveillance and therapeutic efficacy monitoring should be intensified in the region to further assess the extent and spread of artemisinin resistance., (Copyright © 2018 American Society for Microbiology.)