A machine learning approach identifies distinct early-symptom cluster phenotypes which correlate with hospitalization, failure to return to activities, and prolonged COVID-19 symptoms.

BACKGROUND: Accurate COVID-19 prognosis is a critical aspect of acute and long-term clinical management. We identified discrete clusters of early stage-symptoms which may delineate groups with distinct disease severity phenotypes, including risk of developing long-term symptoms and associated inflammatory profiles. METHODS: 1,273 SARS-CoV-2 positive U.S. Military Health System beneficiaries with quantitative symptom scores (FLU-PRO Plus) were included in this analysis. We employed machine-learning approaches to identify symptom clusters and compared risk of hospitalization, long-term symptoms, as well as peak CRP and IL-6 concentrations. RESULTS: We identified three distinct clusters of participants based on their FLU-PRO Plus symptoms: cluster 1 ("Nasal cluster") is highly correlated with reporting runny/stuffy nose and sneezing, cluster 2 ("Sensory cluster") is highly correlated with loss of smell or taste, and cluster 3 ("Respiratory/Systemic cluster") is highly correlated with the respiratory (cough, trouble breathing, among others) and systemic (body aches, chills, among others) domain symptoms. Participants in the Respiratory/Systemic cluster were twice as likely as those in the Nasal cluster to have been hospitalized, and 1.5 times as likely to report that they had not returned-to-activities, which remained significant after controlling for confounding covariates (P < 0.01). Respiratory/Systemic and Sensory clusters were more likely to have symptoms at six-months post-symptom-onset (P = 0.03). We observed higher peak CRP and IL-6 in the Respiratory/Systemic cluster (P < 0.01). CONCLUSIONS: We identified early symptom profiles potentially associated with hospitalization, return-to-activities, long-term symptoms, and inflammatory profiles. These findings may assist in patient prognosis, including prediction of long COVID risk.

The development of broad-spectrum antiviral medical countermeasures to treat viral hemorrhagic fevers caused by natural or weaponized virus infections.

The Joint Program Executive Office for Chemical, Biological, Radiological, and Nuclear Defense (JPEO-CBRND) began development of a broad-spectrum antiviral countermeasure against deliberate use of high-consequence viral hemorrhagic fevers (VHFs) in 2016. The effort featured comprehensive preclinical research, including laboratory testing and rapid advancement of lead molecules into nonhuman primate (NHP) models of Ebola virus disease (EVD). Remdesivir (GS-5734, Veklury, Gilead Sciences) was the first small molecule therapeutic to successfully emerge from this effort. Remdesivir is an inhibitor of RNA-dependent RNA polymerase, a viral enzyme that is essential for viral replication. Its robust potency and broad-spectrum antiviral activity against certain RNA viruses including Ebola virus and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) led to its clinical evaluation in randomized, controlled trials (RCTs) in human patients during the 2018 EVD outbreak in the Democratic Republic of the Congo (DRC) and the ongoing Coronavirus Disease 2019 (COVID-19) pandemic today. Remdesivir was recently approved by the US Food and Drug Administration (FDA) for the treatment of COVID-19 requiring hospitalization. Substantial gaps remain in improving the outcomes of acute viral infections for patients afflicted with both EVD and COVID-19, including how to increase therapeutic breadth and strategies for the prevention and treatment of severe disease. Combination therapy that joins therapeutics with complimentary mechanisms of action appear promising, both preclinically and in RCTs. Importantly, significant programmatic challenges endure pertaining to a clear drug and biological product development pathway for therapeutics targeting biodefense and emerging pathogens when human efficacy studies are not ethical or feasible. For example, remdesivir's clinical development was facilitated by outbreaks of Ebola and SARS-CoV-2; as such, the development pathway employed for remdesivir is likely to be the exception rather than the rule. The current regulatory licensure pathway for therapeutics targeting rare, weaponizable VHF agents is likely to require use of FDA's established Animal Rule (21 CFR 314.600-650 for drugs; 21 CFR 601.90-95 for biologics). The FDA may grant marketing approval based on adequate and well-controlled animal efficacy studies when the results of those studies establish that the drug is safe and likely to produce clinical benefit in humans. In practical terms, this is anticipated to include a series of rigorous, well-documented, animal challenge studies, to include aerosol challenge, combined with human safety data. While small clinical studies against naturally occurring, high-consequence pathogens are typically performed where possible, approval for the therapeutics currently under development against biodefense pathogens will likely require the Animal Rule pathway utilizing studies in NHPs. We review the development of remdesivir as illustrative of the effort that will be needed to field future therapeutics against highly lethal, infectious agents.

Natural History of Aerosol-Induced Ebola Virus Disease in Rhesus Macaques.

Ebola virus disease (EVD) is a serious global health concern because case fatality rates are approximately 50% due to recent widespread outbreaks in Africa. Well-defined nonhuman primate (NHP) models for different routes of Ebola virus exposure are needed to test the efficacy of candidate countermeasures. In this natural history study, four rhesus macaques were challenged via aerosol with a target titer of 1000 plaque-forming units per milliliter of Ebola virus. The course of disease was split into the following stages for descriptive purposes: subclinical, clinical, and decompensated. During the subclinical stage, high levels of venous partial pressure of carbon dioxide led to respiratory acidemia in three of four of the NHPs, and all developed lymphopenia. During the clinical stage, all animals had fever, viremia, and respiratory alkalosis. The decompensatory stage involved coagulopathy, cytokine storm, and liver and renal injury. These events were followed by hypotension, elevated lactate, metabolic acidemia, shock and mortality similar to historic intramuscular challenge studies. Viral loads in the lungs of aerosol-exposed animals were not distinctly different compared to previous intramuscularly challenged studies. Differences in the aerosol model, compared to intramuscular model, include an extended subclinical stage, shortened clinical stage, and general decompensated stage. Therefore, the shortened timeframe for clinical detection of the aerosol-induced disease can impair timely therapeutic administration. In summary, this nonhuman primate model of aerosol-induced EVD characterizes early disease markers and additional details to enable countermeasure development.

Plasmodium falciparum phenotypic and genotypic resistance profile during the emergence of Piperaquine resistance in Northeastern Thailand.

Malaria remains a public health problem in Thailand, especially along its borders where highly mobile populations can contribute to persistent transmission. This study aimed to determine resistant genotypes and phenotypes of 112 Plasmodium falciparum isolates from patients along the Thai-Cambodia border during 2013-2015. The majority of parasites harbored a pfmdr1-Y184F mutation. A single pfmdr1 copy number had CVIET haplotype of amino acids 72-76 of pfcrt and no pfcytb mutations. All isolates had a single pfk13 point mutation (R539T, R539I, or C580Y), and increased % survival in the ring-stage survival assay (except for R539I). Multiple copies of pfpm2 and pfcrt-F145I were detected in 2014 (12.8%) and increased to 30.4% in 2015. Parasites containing either multiple pfpm2 copies with and without pfcrt-F145I or a single pfpm2 copy with pfcrt-F145I exhibited elevated IC90 values of piperaquine. Collectively, the emergence of these resistance patterns in Thailand near Cambodia border mirrored the reports of dihydroartemisinin-piperaquine treatment failures in the adjacent province of Cambodia, Oddar Meanchey, suggesting a migration of parasites across the border. As malaria elimination efforts ramp up in Southeast Asia, host nations militaries and other groups in border regions need to coordinate the proposed interventions.

Safety and immunogenicity of an inactivated eastern equine encephalitis virus vaccine.

BACKGROUND: Eastern equine encephalitis virus (EEEV) is a mosquito borne alphavirus spread primarily in Atlantic and Gulf Coast regions of the United States. EEEV is the causative agent of a devastating meningoencephalitis syndrome, with approximately 30% mortality and significant morbidity. There is no licensed human vaccine against EEEV. An inactivated EEEV vaccine has been offered under investigational new drug (IND) protocols at the United States Army Medical Research Institute of Infectious Diseases (USAMRIID) since 1976. METHODS: Healthy at-risk laboratory personnel received inactivated PE-6 strain EEEV (TSI-GSD 104) vaccine under two separate IND protocols. Protocol FY 99-11 (2002-2008) had a primary series consisting of doses on day 0, 7, and 28. Protocol FY 06-31 (2008-2016) utilized a primary series with doses on day 0 and 28, and month 6. Participants with an inadequate immune response, plaque reduction neutralization test with 80% cut-off (PRNT80) titer < 40, received booster vaccination. Volunteers with prior EEEV vaccination were eligible to enroll for booster doses based on annual titer evaluation. RESULTS: The FY06-31 dosing schema resulted in significantly greater post-primary series immune response (PRNT80 ≥ 40) rates (84% vs 54%) and geometric mean titers (184.1 vs 39.4). The FY 06-31 dosing schema also resulted in significantly greater cumulative annual immune response rates from 1 to up to 7 years post vaccination (75% vs 59%) and geometric mean of titers (60.1 vs 43.0). The majority of probably or definitely related adverse events were mild and local; there were no probably or definitely related serious adverse events. CONCLUSIONS: Inactivated PE-6 EEEV vaccine is safe and immunogenic in at-risk laboratory personnel. A prolonged primary series, with month 6 dose, significantly improved vaccine immunogenicity both post-primary series and longitudinally on annual titers. Despite decades of safe use under IND, full licensure is not planned due to manufacturing constraints, and ongoing development of alternatives.

Distribution and Temporal Dynamics of Plasmodium falciparum Chloroquine Resistance Transporter Mutations Associated With Piperaquine Resistance in Northern Cambodia.

BACKGROUND: Newly emerged mutations within the Plasmodium falciparum chloroquine resistance transporter (PfCRT) can confer piperaquine resistance in the absence of amplified plasmepsin II (pfpm2). In this study, we estimated the prevalence of co-circulating piperaquine resistance mutations in P. falciparum isolates collected in northern Cambodia from 2009 to 2017. METHODS: The sequence of pfcrt was determined for 410 P. falciparum isolates using PacBio amplicon sequencing or whole genome sequencing. Quantitative polymerase chain reaction was used to estimate pfpm2 and pfmdr1 copy number. RESULTS: Newly emerged PfCRT mutations increased in prevalence after the change to dihydroartemisinin-piperaquine in 2010, with >98% of parasites harboring these mutations by 2017. After 2014, the prevalence of PfCRT F145I declined, being outcompeted by parasites with less resistant, but more fit PfCRT alleles. After the change to artesunate-mefloquine, the prevalence of parasites with amplified pfpm2 decreased, with nearly half of piperaquine-resistant PfCRT mutants having single-copy pfpm2. CONCLUSIONS: The large proportion of PfCRT mutants that lack pfpm2 amplification emphasizes the importance of including PfCRT mutations as part of molecular surveillance for piperaquine resistance in this region. Likewise, it is critical to monitor for amplified pfmdr1 in these PfCRT mutants, as increased mefloquine pressure could lead to mutants resistant to both drugs.

Selection of Cytochrome b Mutants Is Rare among Plasmodium falciparum Patients Failing Treatment with Atovaquone-Proguanil in Cambodia.

Atovaquone-proguanil remains effective against multidrug-resistant Plasmodium falciparum in Southeast Asia, but resistance is mediated by a single point mutation in cytochrome b (cytb) that can arise during treatment. Among 14 atovaquone-proguanil treatment failures in a clinical trial in Cambodia, only one recrudescence harbored the cytb mutation Y268C. Deep sequencing did not detect the mutation at baseline or in the first 3 days of treatment, suggesting that it arose de novo Further sequencing across cytb similarly found no low-frequency cytb mutations that were up-selected from baseline to recrudescence. Copy number amplification in dihydroorotate dehydrogenase (DHODH) and cytb as markers of atovaquone tolerance was also absent. Cytb mutation played a minor role in atovaquone-proguanil treatment failures in an active comparator clinical trial.

The risk of Plasmodium vivax parasitaemia after P. falciparum malaria: An individual patient data meta-analysis from the WorldWide Antimalarial Resistance Network.

BACKGROUND: There is a high risk of Plasmodium vivax parasitaemia following treatment of falciparum malaria. Our study aimed to quantify this risk and the associated determinants using an individual patient data meta-analysis in order to identify populations in which a policy of universal radical cure, combining artemisinin-based combination therapy (ACT) with a hypnozoitocidal antimalarial drug, would be beneficial. METHODS AND FINDINGS: A systematic review of Medline, Embase, Web of Science, and the Cochrane Database of Systematic Reviews identified efficacy studies of uncomplicated falciparum malaria treated with ACT that were undertaken in regions coendemic for P. vivax between 1 January 1960 and 5 January 2018. Data from eligible studies were pooled using standardised methodology. The risk of P. vivax parasitaemia at days 42 and 63 and associated risk factors were investigated by multivariable Cox regression analyses. Study quality was assessed using a tool developed by the Joanna Briggs Institute. The study was registered in the International Prospective Register of Systematic Reviews (PROSPERO: CRD42018097400). In total, 42 studies enrolling 15,341 patients were included in the analysis, including 30 randomised controlled trials and 12 cohort studies. Overall, 14,146 (92.2%) patients had P. falciparum monoinfection and 1,195 (7.8%) mixed infection with P. falciparum and P. vivax. The median age was 17.0 years (interquartile range [IQR] = 9.0-29.0 years; range = 0-80 years), with 1,584 (10.3%) patients younger than 5 years. 2,711 (17.7%) patients were treated with artemether-lumefantrine (AL, 13 studies), 651 (4.2%) with artesunate-amodiaquine (AA, 6 studies), 7,340 (47.8%) with artesunate-mefloquine (AM, 25 studies), and 4,639 (30.2%) with dihydroartemisinin-piperaquine (DP, 16 studies). 14,537 patients (94.8%) were enrolled from the Asia-Pacific region, 684 (4.5%) from the Americas, and 120 (0.8%) from Africa. At day 42, the cumulative risk of vivax parasitaemia following treatment of P. falciparum was 31.1% (95% CI 28.9-33.4) after AL, 14.1% (95% CI 10.8-18.3) after AA, 7.4% (95% CI 6.7-8.1) after AM, and 4.5% (95% CI 3.9-5.3) after DP. By day 63, the risks had risen to 39.9% (95% CI 36.6-43.3), 42.4% (95% CI 34.7-51.2), 22.8% (95% CI 21.2-24.4), and 12.8% (95% CI 11.4-14.5), respectively. In multivariable analyses, the highest rate of P. vivax parasitaemia over 42 days of follow-up was in patients residing in areas of short relapse periodicity (adjusted hazard ratio [AHR] = 6.2, 95% CI 2.0-19.5; p = 0.002); patients treated with AL (AHR = 6.2, 95% CI 4.6-8.5; p < 0.001), AA (AHR = 2.3, 95% CI 1.4-3.7; p = 0.001), or AM (AHR = 1.4, 95% CI 1.0-1.9; p = 0.028) compared with DP; and patients who did not clear their initial parasitaemia within 2 days (AHR = 1.8, 95% CI 1.4-2.3; p < 0.001). The analysis was limited by heterogeneity between study populations and lack of data from very low transmission settings. Study quality was high. CONCLUSIONS: In this meta-analysis, we found a high risk of P. vivax parasitaemia after treatment of P. falciparum malaria that varied significantly between studies. These P. vivax infections are likely attributable to relapses that could be prevented with radical cure including a hypnozoitocidal agent; however, the benefits of such a novel strategy will vary considerably between geographical areas.

Correction: Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: A systematic review and meta-analysis.

[This corrects the article DOI: 10.1371/journal.pmed.1003084.].

Piperaquine resistant Cambodian Plasmodium falciparum clinical isolates: in vitro genotypic and phenotypic characterization.

BACKGROUND: High rates of dihydroartemisinin-piperaquine (DHA-PPQ) treatment failures have been documented for uncomplicated Plasmodium falciparum in Cambodia. The genetic markers plasmepsin 2 (pfpm2), exonuclease (pfexo) and chloroquine resistance transporter (pfcrt) genes are associated with PPQ resistance and are used for monitoring the prevalence of drug resistance and guiding malaria drug treatment policy. METHODS: To examine the relative contribution of each marker to PPQ resistance, in vitro culture and the PPQ survival assay were performed on seventeen P. falciparum isolates from northern Cambodia, and the presence of E415G-Exo and pfcrt mutations (T93S, H97Y, F145I, I218F, M343L, C350R, and G353V) as well as pfpm2 copy number polymorphisms were determined. Parasites were then cloned by limiting dilution and the cloned parasites were tested for drug susceptibility. Isobolographic analysis of several drug combinations for standard clones and newly cloned P. falciparum Cambodian isolates was also determined. RESULTS: The characterization of culture-adapted isolates revealed that the presence of novel pfcrt mutations (T93S, H97Y, F145I, and I218F) with E415G-Exo mutation can confer PPQ-resistance, in the absence of pfpm2 amplification. In vitro testing of PPQ resistant parasites demonstrated a bimodal dose-response, the existence of a swollen digestive vacuole phenotype, and an increased susceptibility to quinine, chloroquine, mefloquine and lumefantrine. To further characterize drug sensitivity, parental parasites were cloned in which a clonal line, 14-B5, was identified as sensitive to artemisinin and piperaquine, but resistant to chloroquine. Assessment of the clone against a panel of drug combinations revealed antagonistic activity for six different drug combinations. However, mefloquine-proguanil and atovaquone-proguanil combinations revealed synergistic antimalarial activity. CONCLUSIONS: Surveillance for PPQ resistance in regions relying on DHA-PPQ as the first-line treatment is dependent on the monitoring of molecular markers of drug resistance. P. falciparum harbouring novel pfcrt mutations with E415G-exo mutations displayed PPQ resistant phenotype. The presence of pfpm2 amplification was not required to render parasites PPQ resistant suggesting that the increase in pfpm2 copy number alone is not the sole modulator of PPQ resistance. Genetic background of circulating field isolates appear to play a role in drug susceptibility and biological responses induced by drug combinations. The use of latest field isolates may be necessary for assessment of relevant drug combinations against P. falciparum strains and when down-selecting novel drug candidates.

Recent successes in therapeutics for Ebola virus disease: no time for complacency.

The PALM trial in the Democratic Republic of the Congo identified a statistically significant survival benefit for two monoclonal antibody-based therapeutics in the treatment of acute Ebola virus disease; however, substantial gaps remain in improving the outcomes of acute Ebola virus disease and for the survivors. Ongoing efforts are needed to develop more effective strategies, particularly for individuals with severe disease, for prevention and treatment of viral persistence in immune-privileged sites, for optimisation of post-exposure prophylaxis, and to increase therapeutic breadth. As antibody-based approaches are identified and advanced, promising small-molecule antivirals currently in clinical stage development should continue to be evaluated for filovirus diseases, with consideration of their added value in combination approaches with bundled supportive care, their penetration in tissues of interest, the absence of interaction with glycoprotein-based vaccines, and filoviral breadth.

Natural History of Aerosol Induced Lassa Fever in Non­Human Primates.

Lassa virus (LASV), an arenavirus causing Lassa fever, is endemic to West Africa with up to 300,000 cases and between 5000 and 10,000 deaths per year. Rarely seen in the United States, Lassa virus is a CDC category A biological agent inasmuch deliberate aerosol exposure can have high mortality rates compared to naturally acquired infection. With the need for an animal model, specific countermeasures remain elusive as there is no FDA-approved vaccine. This natural history of aerosolized Lassa virus exposure in Macaca fascicularis was studied under continuous telemetric surveillance. The macaque response to challenge was largely analogous to severe human disease with fever, tachycardia, hypotension, and tachypnea. During initial observations, an increase trend of activated monocytes positive for viral glycoprotein was accompanied by lymphocytopenia. Disease uniformly progressed to high viremia followed by low anion gap, alkalosis, anemia, and thrombocytopenia. Hypoproteinemia occurred late in infection followed by increased levels of white blood cells, cytokines, chemokines, and biochemical markers of liver injury. Viral nucleic acids were detected in tissues of three non­survivors at endpoint, but not in the lone survivor. This study provides useful details to benchmark a pivotal model of Lassa fever in support of medical countermeasure development for both endemic disease and traditional biodefense purposes.

Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: A systematic review and meta-analysis.

BACKGROUND: The radical cure of Plasmodium vivax and P. ovale requires treatment with primaquine or tafenoquine to clear dormant liver stages. Either drug can induce haemolysis in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency, necessitating screening. The reference diagnostic method for G6PD activity is ultraviolet (UV) spectrophotometry; however, a universal G6PD activity threshold above which these drugs can be safely administered is not yet defined. Our study aimed to quantify assay-based variation in G6PD spectrophotometry and to explore the diagnostic implications of applying a universal threshold. METHODS AND FINDINGS: Individual-level data were pooled from studies that used G6PD spectrophotometry. Studies were identified via PubMed search (25 April 2018) and unpublished contributions from contacted authors (PROSPERO: CRD42019121414). Studies were excluded if they assessed only individuals with known haematological conditions, were family studies, or had insufficient details. Studies of malaria patients were included but analysed separately. Included studies were assessed for risk of bias using an adapted form of the Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool. Repeatability and intra- and interlaboratory variability in G6PD activity measurements were compared between studies and pooled across the dataset. A universal threshold for G6PD deficiency was derived, and its diagnostic performance was compared to site-specific thresholds. Study participants (n = 15,811) were aged between 0 and 86 years, and 44.4% (7,083) were women. Median (range) activity of G6PD normal (G6PDn) control samples was 10.0 U/g Hb (6.3-14.0) for the Trinity assay and 8.3 U/g Hb (6.8-15.6) for the Randox assay. G6PD activity distributions varied significantly between studies. For the 13 studies that used the Trinity assay, the adjusted male median (AMM; a standardised metric of 100% G6PD activity) varied from 5.7 to 12.6 U/g Hb (p < 0.001). Assay precision varied between laboratories, as assessed by variance in control measurements (from 0.1 to 1.5 U/g Hb; p < 0.001) and study-wise mean coefficient of variation (CV) of replicate measures (from 1.6% to 14.9%; p < 0.001). A universal threshold of 100% G6PD activity was defined as 9.4 U/g Hb, yielding diagnostic thresholds of 6.6 U/g Hb (70% activity) and 2.8 U/g Hb (30% activity). These thresholds diagnosed individuals with less than 30% G6PD activity with study-wise sensitivity from 89% (95% CI: 81%-94%) to 100% (95% CI: 96%-100%) and specificity from 96% (95% CI: 89%-99%) to 100% (100%-100%). However, when considering intermediate deficiency (<70% G6PD activity), sensitivity fell to a minimum of 64% (95% CI: 52%-75%) and specificity to 35% (95% CI: 24%-46%). Our ability to identify underlying factors associated with study-level heterogeneity was limited by the lack of availability of covariate data and diverse study contexts and methodologies. CONCLUSIONS: Our findings indicate that there is substantial variation in G6PD measurements by spectrophotometry between sites. This is likely due to variability in laboratory methods, with possible contribution of unmeasured population factors. While an assay-specific, universal quantitative threshold offers robust diagnosis at the 30% level, inter-study variability impedes performance of universal thresholds at the 70% level. Caution is advised in comparing findings based on absolute G6PD activity measurements across studies. Novel handheld quantitative G6PD diagnostics may allow greater standardisation in the future.

Prevalence of CYP2D6 Genotypes and Predicted Phenotypes in a Cohort of Cambodians at High Risk for Infections with Plasmodium vivax.

Clinical failure of primaquine (PQ) has been demonstrated in people with CYP450 2D6 genetic polymorphisms that result in reduced or no enzyme activity. The distribution of CYP2D6 genotypes and predicted phenotypes in the Cambodian population is not well described. Surveys in other Asian countries have shown an approximate 50% prevalence of the reduced activity CYP2D6 allele *10, which could translate into increased risk of PQ radical cure failure and repeated relapses, making interruption of transmission and malaria elimination difficult to achieve. We determined CYP2D6 genotypes from 96 volunteers from Oddor Meanchey Province, Cambodia, an area endemic for Plasmodium vivax. We found a 54.2% frequency of the *10 allele, but in approximately half of our subjects, it was paired with a normal activity allele, either *1 or *2. The prevalence of *5, a null allele, was 9.4%. Overall predicted phenotype percentages were normal metabolizers, 46%; intermediate metabolizers, 52%; and poor metabolizers, 1%.

Pharmacogenetic assessment of tafenoquine efficacy in patients with Plasmodium vivax malaria.

Plasmodium vivax has the largest geographic range of human malaria species and is challenging to manage and eradicate due to its ability to establish a dormant liver stage, the hypnozoite, which can reactivate leading to relapse. Until recently, the only treatment approved to kill hypnozoites was the 8-aminoquinoline, primaquine, requiring daily treatment for 14 days. Tafenoquine, an 8-aminoquinoline single-dose treatment with activity against P. vivax hypnozoites, has recently been approved by the US Food and Drug Administration and Australian Therapeutic Goods Administration for the radical cure of P. vivax malaria in patients 16 years and older. We conducted an exploratory pharmacogenetic analysis (GSK Study 208099) to assess the role of host genome-wide variation on tafenoquine efficacy in patients with P. vivax malaria using data from three GSK clinical trials, GATHER and DETECTIVE Part 1 and Part 2. Recurrence-free efficacy at 6 and 4 months and time to recurrence up to 6 months postdosing were analyzed in 438 P. vivax malaria patients treated with tafenoquine. Among the approximately 10.6 million host genetic variants analyzed, two signals reached genome-wide significance (P value ≤ 5 × 10). rs62103056, and variants in a chromosome 12 intergenic region, were associated with recurrence-free efficacy at 6 and 4 months, respectively. Neither of the signals has an obvious biological rationale and would need replication in an independent population. This is the first genome-wide association study to evaluate genetic influence on response to tafenoquine in P. vivax malaria.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces.

BACKGROUND: Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam. METHODS: The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns. RESULTS: Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns. CONCLUSIONS: Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential.

BACKGROUND: Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions. METHODS: Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia. RESULTS: While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region. CONCLUSIONS: These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

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

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

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns.

Estimates of Plasmodium falciparum migration may inform strategies for malaria elimination. Here we elucidate fine-scale parasite population structure and infer recent migration across Southeast Asia using identity-by-descent (IBD) approaches based on genome-wide single nucleotide polymorphisms called in 1722 samples from 54 districts. IBD estimates are consistent with isolation-by-distance. We observe greater sharing of larger IBD segments between artemisinin-resistant parasites versus sensitive parasites, which is consistent with the recent spread of drug resistance. Our IBD analyses reveal actionable patterns, including isolated parasite populations, which may be prioritized for malaria elimination, as well as asymmetrical migration identifying potential sources and sinks of migrating parasites.

Single-Dose Tafenoquine to Prevent Relapse of Plasmodium vivax Malaria.

BACKGROUND: Treatment of Plasmodium vivax malaria requires the clearing of asexual parasites, but relapse can be prevented only if dormant hypnozoites are cleared from the liver (a treatment termed "radical cure"). Tafenoquine is a single-dose 8-aminoquinoline that has recently been registered for the radical cure of P. vivax. METHODS: This multicenter, double-blind, double-dummy, parallel group, randomized, placebo-controlled trial was conducted in Ethiopia, Peru, Brazil, Cambodia, Thailand, and the Philippines. We enrolled 522 patients with microscopically confirmed P. vivax infection (>100 to <100,000 parasites per microliter) and normal glucose-6-phosphate dehydrogenase (G6PD) activity (with normal activity defined as ≥70% of the median value determined at each trial site among 36 healthy male volunteers who were otherwise not involved in the trial). All patients received a 3-day course of chloroquine (total dose of 1500 mg). In addition, patients were assigned to receive a single 300-mg dose of tafenoquine on day 1 or 2 (260 patients), placebo (133 patients), or a 15-mg dose of primaquine once daily for 14 days (129 patients). The primary outcome was the Kaplan-Meier estimated percentage of patients who were free from recurrence at 6 months, defined as P. vivax clearance without recurrent parasitemia. RESULTS: In the intention-to-treat population, the percentage of patients who were free from recurrence at 6 months was 62.4% in the tafenoquine group (95% confidence interval [CI], 54.9 to 69.0), 27.7% in the placebo group (95% CI, 19.6 to 36.6), and 69.6% in the primaquine group (95% CI, 60.2 to 77.1). The hazard ratio for the risk of recurrence was 0.30 (95% CI, 0.22 to 0.40) with tafenoquine as compared with placebo (P<0.001) and 0.26 (95% CI, 0.18 to 0.39) with primaquine as compared with placebo (P<0.001). Tafenoquine was associated with asymptomatic declines in hemoglobin levels, which resolved without intervention. CONCLUSIONS: Single-dose tafenoquine resulted in a significantly lower risk of P. vivax recurrence than placebo in patients with phenotypically normal G6PD activity. (Funded by GlaxoSmithKline and Medicines for Malaria Venture; DETECTIVE ClinicalTrials.gov number, NCT01376167 .).