Travel-related malaria is regularly encountered in the United States, and the U.S. Centers for Disease Control and Prevention (CDC) characterizes Plasmodium falciparum drug-resistance genotypes routinely for travel-related cases. An important aspect of antimalarial drug resistance is understanding its geographic distribution. However, specimens submitted to CDC laboratories may have missing, incomplete, or inaccurate travel data. To complement genotyping for drug-resistance markers Pfcrt, Pfmdr1, Pfk13, Pfdhps, Pfdhfr, and PfcytB at CDC, amplicons of Pfs47 and Pfcpmp are also sequenced as markers of geographic origin. Here, a bi-allele likelihood (BALK) classifier was trained using Pfs47 and Pfcpmp sequences from published P. falciparum genomes of known geographic origin to classify clinical genotypes to a continent. Among P. falciparum-positive blood samples received at CDC for drug-resistance genotyping from 2018 to 2021 (n = 380), 240 included a travel history with the submission materials, though 6 were excluded due to low sequence quality. Classifications obtained for the remaining 234 were compared to their travel histories. Classification results were over 96% congruent with reported travel for clinical samples, and with collection sites for field isolates. Among travel-related samples, only two incongruent results occurred; a specimen submitted citing Costa Rican travel classified to Africa, and a specimen with travel referencing Sierra Leone classified to Asia. Subsequently, the classifier was applied to specimens with unreported travel histories (n = 140; 5 were excluded due to low sequence quality). For the remaining 135 samples, geographic classification data were paired with results generated using CDC's Malaria Resistance Surveillance (MaRS) protocol, which detects single-nucleotide polymorphisms in and generates haplotypes for Pfcrt, Pfmdr1, Pfk13, Pfdhps, Pfdhfr, and PfcytB. Given the importance of understanding the geographic distribution of antimalarial drug resistance, this work will complement domestic surveillance efforts by expanding knowledge on the geographic origin of drug-resistant P. falciparum entering the USA.

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. © 2024

All errors found in this paper are due to six samples included in the paper incorrectly assigned as being from Queensland, Australia. The report of Strongyloides fuelleborni infections from Australia made in this paper was incorrect, as those samples in fact originated in Guinea-Bissau and Senegal. | | There is an error in Table 4. Specimen Human 333_Au from Queensland (Australia) should be listed as Human 333_GuBi from Guinea-Bissau. Specimen Human 368_16_Au from Queensland (Australia) should be listed as Human 368_16_Se from Senegal. Specimen Human 378_Au from Queensland (Australia) should be listed as Human 378_Bo from Bolivia. Specimen Human 507_Au from Queensland (Australia) should be listed as Human 507_Ni from Nigeria. Specimen Human 524_Au from Queensland (Australia) should be listed as Human 524_Ni from Nigeria. Specimen Human 563_Au from Queensland (Australia) should be listed as Human 563_GuBi from Guinea-Bissau. The authors have provided a corrected Table 4 with the corrected specimens and locations in red.

We thank Rasmussen and Ringwald for further highlighting the importance of routine monitoring of antimalarial drug efficacy in sub-Saharan Africa, including Angola (1). Longitudinal monitoring is critical to identify potential new, persistent, and/or expanding foci of parasite resistance to available drugs. In 3 of the last 4 rounds, artemether-lumefantrine (AL) was estimated to have an efficacy of <90% at one of the three sentinel sites in Angola. To our knowledge, in sub-Saharan Africa, only Angola and Burkina Faso (2) have shown AL efficacy of <90% across multiple therapeutic efficacy study (TES) rounds. Thus, we chose a title to highlight this persistent concern. | | We concur that the significance of the high rates of day 2 slide positivity in Lunda Sul Province is not fully known, and as pointed out, there may be various explanations for this finding. Measuring drug levels is resource intensive and not feasible every year, but this could help rule out underdosing in future studies. However, we believe our study procedures, as described in this and previous studies, are robust and thus make systematic underdosing unlikely. We have always strictly adhered to WHO guidelines, including hemoglobin criteria and analysis of day 1 severe cases, to inform our classifications.

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.

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.

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 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.

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.

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.

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. © 2022 by the authors. Licensee MDPI, Basel, Switzerland.

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.

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: 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.

BACKGROUND: Recent studies show 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 deep amplicon sequencing (TADS) of 116 isolates collected in two areas of Cameroon to genotype the major drug resistance genes k13, crt, mdr1, dhfr, dhps, and the cytochrome b (cytb) in P. 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 non-synonymous mutations was observed across all the genes, except Pfcytb, suggesting continued selection pressure. CONCLUSIONS: The results from this study support the continued use of artemisinin combination therapy (ACT) for treatment and administration of sulphadoxine-pyrimethamine for intermittent preventive therapy in pregnant women and for seasonal chemoprevention in these study sites in Cameroon.

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.

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.

BACKGROUND: Biennial therapeutic efficacy monitoring is a crucial activity for ensuring efficacy of currently used artemisinin-based combination therapy in Angola. METHODS: Children with acute uncomplicated P. falciparum infection in sentinel sites in Benguela, Zaire, and Lunda Sul Provinces were treated with artemether-lumefantrine (AL) or artesunate amodiaquine (ASAQ) and followed for 28 days to assess clinical and parasitological response. Molecular correction was performed using seven microsatellite markers. Samples from treatment failures were genotyped for the pfk13, pfcrt, and pfmdr1 genes. RESULTS: 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-95%) for the AL arm in Lunda Sul, 92.2% (95%CI: 87-98%) for AL in Zaire, 95.6% (95%CI: 91-100%) for ASAQ in Zaire, 98.4% (95%CI: 96-100%) for AL in Benguela, and 100% for ASAQ in Benguela and Lunda Sul. All 103 analyzed samples had wildtype pfk13 sequences. The 76T pfcrt allele was found in most (92%, 11/12) ASAQ late failure samples but only 16% (4/25) of AL failure samples. The N86 pfmdr1 allele was found in 97% (34/35) of treatment failures. CONCLUSION: 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, observed ASAQ efficacy has not been below 95% to date in Angola, including this latest round.

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.

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.

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.

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 Labe 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 Labe. 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 Labe. 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.

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.

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 Sanita (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.

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.

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.

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.