Clinical exome sequencing uncovers genetic disorders in neonates with suspected hypoxic-ischemic encephalopathy: A retrospective analysis.

Hypoxic-ischemic encephalopathy (HIE) occurs in up to 7 out of 1000 births and accounts for almost a quarter of neonatal deaths worldwide. Despite the name, many newborns with HIE have little evidence of perinatal hypoxia. We hypothesized that some infants with HIE have genetic disorders that resemble encephalopathy. We reviewed genetic results for newborns with HIE undergoing exome or genome sequencing at a clinical laboratory (2014-2022). Neonates were included if they had a diagnosis of HIE and were delivered ≥35 weeks. Neonates were excluded for cardiopulmonary pathology resulting in hypoxemia or if neuroimaging suggested postnatal hypoxic-ischemic injury. Of 24 patients meeting inclusion criteria, six (25%) were diagnosed with a genetic condition. Four neonates had variants at loci linked to conditions with phenotypic features resembling HIE, including KIF1A, GBE1, ACTA1, and a 15q13.3 deletion. Two additional neonates had variants in genes not previously associated with encephalopathy, including DUOX2 and PTPN11. Of the six neonates with a molecular diagnosis, two had isolated HIE without apparent comorbidities to suggest a genetic disorder. Genetic diagnoses were identified among neonates with and without sentinel labor events, abnormal umbilical cord gasses, and low Apgar scores. These results suggest that genetic evaluation is clinically relevant for patients with perinatal HIE.

Deep Sequencing to Detect Diversity of Trypanosoma cruzi Infection in Patients Coinfected With Human Immunodeficiency Virus and Chagas Disease.

Chagas disease, caused by Trypanosoma cruzi, can reactivate and cause severe acute disease in immunocompromised patients such as those infected with human immunodeficiency virus (HIV). We conducted amplicon deep sequencing of a 327-bp fragment of the tcscd5 gene using an Ion Torrent PGM directly from clinical samples from HIV patients with high parasitemia. We describe the within-host diversity, both characterizing the discrete typing unit of the infections and confirming the presence of multistrain infections, directly from clinical samples. This method can rapidly provide information on the genetic diversity of T. cruzi infection, which can have direct impacts on clinical disease.

Video education about genetic privacy and patient perspectives about sharing prenatal genetic data: a randomized trial.

BACKGROUND: Laboratories offering cell-free DNA often reserve the right to share prenatal genetic data for research or even commercial purposes, and obtain this permission on the patient consent form. Although it is known that nonpregnant patients are often reluctant to share their genetic data for research, pregnant patients' knowledge of, and opinions about, genetic data privacy are unknown. OBJECTIVE: We investigated whether pregnant patients who had already undergone cell-free DNA screening were aware that genetic data derived from cell-free DNA may be shared for research. Furthermore, we examined whether pregnant patients exposed to video education about the Genetic Information Nondiscrimination Act-a federal law that mandates workplace and health insurance protections against genetic discrimination-were more willing to share cell-free DNA-related genetic data for research than pregnant patients who were unexposed. STUDY DESIGN: In this randomized controlled trial (ClinicalTrials.gov Identifier: NCT04420858), English-speaking patients with singleton pregnancies who underwent cell-free DNA and subsequently presented at 17 0/7 to 23 6/7 weeks of gestation for a detailed anatomy scan were randomized 1:1 to a control or intervention group. Both groups viewed an infographic about cell-free DNA. In addition, the intervention group viewed an educational video about the Genetic Information Nondiscrimination Act. The primary outcomes were knowledge about, and willingness to share, prenatal genetic data from cell-free DNA by commercial laboratories for nonclinical purposes, such as research. The secondary outcomes included knowledge about existing genetic privacy laws, knowledge about the potential for reidentification of anonymized genetic data, and acceptability of various use and sharing scenarios for prenatal genetic data. Eighty-one participants per group were required for 80% power to detect an increase in willingness to share data from 60% to 80% (α=0.05). RESULTS: A total of 747 pregnant patients were screened, and 213 patients were deemed eligible and approached for potential study participation. Of these patients, 163 (76.5%) consented and were randomized; one participant discontinued the intervention, and two participants were excluded from analysis after the intervention when it was discovered that they did not fulfill all eligibility criteria. Overall, 160 (75.1%) of those approached were included in the final analysis. Most patients in the control group (72 [90.0%]) and intervention (76 [97.4%]) group were either unsure about or incorrectly thought that cell-free DNA companies could not share prenatal genetic data for research. Participants in the intervention group were more likely to incorrectly believe that their prenatal genetic data would not be shared for nonclinical purposes than participants in the control group (28.8% in the control group vs 46.2% in the intervention; P=.03). However, video education did not increase participant willingness to share genetic data in multiple scenarios. Non-White participants were less willing than White participants to allow sharing of genetic data specifically for academic research (P<.001). CONCLUSION: Most participants were unaware that their prenatal genetic data may be used for nonclinical purposes. Pregnant patients who were educated about the Genetic Information Nondiscrimination Act were not more willing to share genetic data than those who did not receive this education. Surprisingly, video education about the Genetic Information Nondiscrimination Act led patients to falsely believe that their data would not be shared for research, and participants who identified as racial minorities were less willing to share genetic data. New strategies are needed to improve pregnant patients' understanding of genetic privacy.

Privacy practices using genetic data from cell-free DNA aneuploidy screening.

PURPOSE: Cell-free fetal DNA (cfDNA) analyzes maternal and fetoplacental DNA, generating highly personal genetic information for both mother and fetus. This study aimed to determine how laboratories retain, use, and share genetic information from cfDNA. Other outcomes included laboratories' adherence to American Society of Human Genetics (ASHG) privacy principles, and the readability of privacy policies. METHODS: Laboratories offering cfDNA aneuploidy screening were identified from online searches, curated databases, and a genomics news website. Of 124 laboratories identified, 13 were commercial laboratories offering cfDNA aneuploidy screening in the United States, and were included. Genetic privacy policies from eligible laboratories were identified by reviewing requisition and consent forms, which were obtained online or by direct contact. RESULTS: Most laboratories use prenatal genetic information for research (n = 10, 77%), and more than half (n = 7, 54%) shared genetic information with others. Overall, laboratories inadequately disclosed privacy risks. In a readability analysis, 9 of 11 (82%) laboratories' genetic privacy policies were written at or above a 12th grade reading level. CONCLUSION: Most laboratories allowed for prolonged use and sharing of cfDNA data, demonstrated incomplete adherence to ASHG privacy recommendations, and provided consents written in college-level language. Laboratories should revise their consent forms, and providers should help patients understand these forms.

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.

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

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

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia.

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity.

Background: Amplified copy number in the plasmepsin II/III genes within Plasmodium falciparum has been associated with decreased sensitivity to piperaquine. To examine this association and test whether additional loci might also contribute, we performed a genome-wide association study of ex vivo P. falciparum susceptibility to piperaquine. Methods: Plasmodium falciparum DNA from 183 samples collected primarily from Cambodia was genotyped at 33716 genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models and random forests were used to estimate associations between parasite genotypes and piperaquine susceptibility. Candidate polymorphisms were evaluated for their association with dihydroartemisinin-piperaquine treatment outcomes in an independent dataset. Results: Single nucleotide polymorphisms on multiple chromosomes were associated with piperaquine 90% inhibitory concentrations (IC90) in a genome-wide analysis. Fine-mapping of genomic regions implicated in genome-wide analyses identified multiple SNPs in linkage disequilibrium with each other that were significantly associated with piperaquine IC90, including a novel mutation within the gene encoding the P. falciparum chloroquine resistance transporter, PfCRT. This mutation (F145I) was associated with dihydroartemisinin-piperaquine treatment failure after adjusting for the presence of amplified plasmepsin II/III, which was also associated with decreased piperaquine sensitivity. Conclusions: Our data suggest that, in addition to plasmepsin II/III copy number, other loci, including pfcrt, may also be involved in piperaquine resistance.

A Prolonged Outbreak of KPC-3-Producing Enterobacter cloacae and Klebsiella pneumoniae Driven by Multiple Mechanisms of Resistance Transmission at a Large Academic Burn Center.

Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacter cloacae has been recently recognized in the United States. Whole-genome sequencing (WGS) has become a useful tool for analysis of outbreaks and for determining transmission networks of multidrug-resistant organisms in health care settings, including carbapenem-resistant Enterobacteriaceae (CRE). We experienced a prolonged outbreak of CRE E. cloacae and K. pneumoniae over a 3-year period at a large academic burn center despite rigorous infection control measures. To understand the molecular mechanisms that sustained this outbreak, we investigated the CRE outbreak isolates by using WGS. Twenty-two clinical isolates of CRE, including E. cloacae (n = 15) and K. pneumoniae (n = 7), were sequenced and analyzed genetically. WGS revealed that this outbreak, which seemed epidemiologically unlinked, was in fact genetically linked over a prolonged period. Multiple mechanisms were found to account for the ongoing outbreak of KPC-3-producing E. cloacae and K. pneumoniae This outbreak was primarily maintained by a clonal expansion of E. cloacae sequence type 114 (ST114) with distribution of multiple resistance determinants. Plasmid and transposon analyses suggested that the majority of blaKPC-3 was transmitted via an identical Tn4401b element on part of a common plasmid. WGS analysis demonstrated complex transmission dynamics within the burn center at levels of the strain and/or plasmid in association with a transposon, highlighting the versatility of KPC-producing Enterobacteriaceae in their ability to utilize multiple modes to resistance gene propagation.

Genomic Analysis of Multidrug-Resistant Escherichia coli from North Carolina Community Hospitals: Ongoing Circulation of CTX-M-Producing ST131-H30Rx and ST131-H30R1 Strains.

Escherichia coli sequence type 131 (ST131) predominates globally among multidrug-resistant (MDR) E. coli strains. We used whole-genome sequencing (WGS) to investigate 63 MDR E. coli isolates from 7 North Carolina community hospitals (2010 to 2015). Of these, 39 (62%) represented ST131, including 37 (95%) from the ST131-H30R subclone: 10 (27%) from its H30R1 subset and 27 (69%) from its H30Rx subset. ST131 core genomes differed by a median of 15 (range, 0 to 490) single-nucleotide variants (SNVs) overall versus only 7 within H30R1 (range, 3 to 12 SNVs) and 11 within H30Rx (range, 0 to 21). The four isolates with identical core genomes were all H30Rx. Epidemiological and clinical characteristics did not vary significantly by strain type, but many patients with MDR E. coli or H30Rx infection were critically ill and had poor outcomes. H30Rx isolates characteristically exhibited fluoroquinolone resistance and CTX-M-15 production, had a high prevalence of trimethoprim-sulfamethoxazole resistance (89%), sul1 (89%), and dfrA17 (85%), and were enriched for specific virulence traits, and all qualified as extraintestinal pathogenic E. coli The high overall prevalence of CTX-M-15 appeared to be possibly attributable to its association with the ST131-H30Rx subclone and IncF[F2:A1:B-] plasmids. Some phylogenetically clustered non-ST131 MDR E. coli isolates also had distinctive serotypes/fimH types, fluoroquinolone mutations, CTX-M variants, and IncF types. Thus, WGS analysis of our community hospital source MDR E. coli isolates suggested ongoing circulation and differentiation of E. coli ST131 subclones, with clonal segregation of CTX-M variants, other resistance genes, Inc-type plasmids, and virulence genes.

Microscopic Plasmodium falciparum Gametocytemia and Infectivity to Mosquitoes in Cambodia.

Although gametocytes are essential for malaria transmission, in Africa many falciparum-infected persons without smear-detectable gametocytes still infect mosquitoes. To see whether the same is true in Southeast Asia, we determined the infectiousness of 119 falciparum-infected Cambodian adults to Anopheles dirus mosquitoes by membrane feeding. Just 5.9% of subjects infected mosquitoes. The 8.4% of patients with smear-detectable gametocytes were >20 times more likely to infect mosquitoes than those without and were the source of 96% of all mosquito infections. In low-transmission settings, targeting transmission-blocking interventions to those with microscopic gametocytemia may have an outsized effect on malaria control and elimination.

Absence of putative artemisinin resistance mutations among Plasmodium falciparum in Sub-Saharan Africa: a molecular epidemiologic study.

Plasmodium falciparum parasites that are resistant to artemisinins have been detected in Southeast Asia. Resistance is associated with several polymorphisms in the parasite's K13-propeller gene. The molecular epidemiology of these artemisinin resistance genotypes in African parasite populations is unknown. We developed an assay to quantify rare polymorphisms in parasite populations that uses a pooled deep-sequencing approach to score allele frequencies, validated it by evaluating mixtures of laboratory parasite strains, and then used it to screen P. falciparum parasites from >1100 African infections collected since 2002 from 14 sites across sub-Saharan Africa. We found no mutations in African parasite populations that are associated with artemisinin resistance in Southeast Asian parasites. However, we observed 15 coding mutations, including 12 novel mutations, and limited allele sharing between parasite populations, consistent with a large reservoir of naturally occurring K13-propeller variation. Although polymorphisms associated with artemisinin resistance in P. falciparum in Southeast Asia are not prevalent in sub-Saharan Africa, numerous K13-propeller coding polymorphisms circulate in Africa. Although their distributions do not support a widespread selective sweep for an artemisinin-resistant phenotype, the impact of these mutations on artemisinin susceptibility is unknown and will require further characterization. Rapid, scalable molecular surveillance offers a useful adjunct in tracking and containing artemisinin resistance.

Next-Generation Sequencing and Comparative Analysis of Sequential Outbreaks Caused by Multidrug-Resistant Acinetobacter baumannii at a Large Academic Burn Center.

Next-generation sequencing (NGS) analysis has emerged as a promising molecular epidemiological method for investigating health care-associated outbreaks. Here, we used NGS to investigate a 3-year outbreak of multidrug-resistant Acinetobacter baumannii (MDRAB) at a large academic burn center. A reference genome from the index case was generated using de novo assembly of PacBio reads. Forty-six MDRAB isolates were analyzed by pulsed-field gel electrophoresis (PFGE) and sequenced using an Illumina platform. After mapping to the index case reference genome, four samples were excluded due to low coverage, leaving 42 samples for further analysis. Multilocus sequence types (MLST) and the presence of acquired resistance genes were also determined from the sequencing data. A transmission network was inferred from genomic and epidemiological data using a Bayesian framework. Based on single-nucleotide variant (SNV) differences, this MDRAB outbreak represented three sequential outbreaks caused by distinct clones. The first and second outbreaks were caused by sequence type 2 (ST2), while the third outbreak was caused by ST79. For the second outbreak, the MLST and PFGE results were discordant. However, NGS-based SNV typing detected a recombination event and consequently enabled a more accurate phylogenetic analysis. The distribution of resistance genes varied among the three outbreaks. The first- and second-outbreak strains possessed a blaOXA-23-like group, while the third-outbreak strains harbored a blaOXA-40-like group. NGS-based analysis demonstrated the superior resolution of outbreak transmission networks for MDRAB and provided insight into the mechanisms of strain diversification between sequential outbreaks through recombination.

Variation in the Microbiota of Ixodes Ticks with Regard to Geography, Species, and Sex.

Ixodes scapularis is the principal vector of Lyme disease on the East Coast and in the upper Midwest regions of the United States, yet the tick is also present in the Southeast, where Lyme disease is absent or rare. A closely related species, I. affinis, also carries the pathogen in the South but does not seem to transmit it to humans. In order to better understand the geographic diversity of the tick, we analyzed the microbiota of 104 adult I. scapularis and 13 adult I. affinis ticks captured in 19 locations in South Carolina, North Carolina, Virginia, Connecticut, and New York. Initially, ticks from 4 sites were analyzed by 454 pyrosequencing. Subsequently, ticks from these sites plus 15 others were analyzed by sequencing with an Illumina MiSeq machine. By both analyses, the microbiomes of female ticks were significantly less diverse than those of male ticks. The dissimilarity between tick microbiomes increased with distance between sites, and the state in which a tick was collected could be inferred from its microbiota. The genus Rickettsia was prominent in all locations. Borrelia was also present in most locations and was present at especially high levels in one site in western Virginia. In contrast, members of the family Enterobacteriaceae were very common in North Carolina I. scapularis ticks but uncommon in I. scapularis ticks from other sites and in North Carolina I. affinis ticks. These data suggest substantial variations in the Ixodes microbiota in association with geography, species, and sex.

Genetic Evidence of Importation of Drug-Resistant Plasmodium falciparum to Guatemala from the Democratic Republic of the Congo.

Imported malaria threatens control and elimination efforts in countries that have low rates of transmission. In 2010, an outbreak of Plasmodium falciparum malaria was reported among United Nations peacekeeping soldiers from Guatemala who had recently returned from the Democratic Republic of the Congo (DRC). Epidemiologic evidence suggested that the soldiers were infected in the DRC, but local transmission could not be ruled out in all cases. We used population genetic analyses of neutral microsatellites to determine the outbreak source. Genetic relatedness was compared among parasites found in samples from the soldiers and parasite populations collected in the DRC and Guatemala; parasites identified in the soldiers were more closely related to those from the DRC. A phylogenetic clustering analysis confirms this identification with >99.9% confidence. Thus, results support the hypothesis that the soldiers likely imported malaria from the DRC. This study demonstrates the utility of molecular genotyping in outbreak investigations.

Multilocus microsatellite genotyping array for investigation of genetic epidemiology of Pneumocystis jirovecii.

Pneumocystis jirovecii is a symbiotic respiratory fungus that causes pneumonia (PcP) in immunosuppressed patients. Because P. jirovecii cannot be reliably cultured in vitro, it has proven difficult to study and gaps in our understanding of the organism persist. The release of a draft genome for the organism opens the door for the development of new genotyping approaches for studying its molecular epidemiology and global population structure. We identified and validated 8 putatively neutral microsatellite markers and 1 microsatellite marker linked to the dihydropteroate synthase gene (dhps), the enzymatic target of sulfa drugs used for PcP prevention and treatment. Using these tools, we analyzed P. jirovecii isolates from HIV-infected patients from three geographically distant populations: Uganda, the United States, and Spain. Among the 8 neutral markers, we observed high levels of allelic heterozygosity (average He, 0.586 to 0.842). Consistent with past reports, we observed limited global population structuring, with only the Ugandan isolates showing minor differentiation from the other two populations. In Ugandan isolates that harbored mutations in dhps, the microsatellite locus linked to dhps demonstrated a depressed He, consistent with positive directional selection for sulfa resistance mutations. Using a subset of these microsatellites, analyses of individual and paired samples from infections in San Francisco, CA, showed reliable typeability within a single infection and high discriminatory power between infections. These features suggest that this novel microsatellite typing approach will be an effective tool for molecular-epidemiological investigations into P. jirovecii population structure, transmission, and drug resistance.

Diversity of Rickettsiales in the microbiome of the lone star tick, Amblyomma americanum.

Ticks are important vectors for many emerging pathogens. However, they are also infected with many symbionts and commensals, often competing for the same niches. In this paper, we characterize the microbiome of Amblyomma americanum (Acari: Ixodidae), the lone star tick, in order to better understand the evolutionary relationships between pathogens and nonpathogens. Multitag pyrosequencing of prokaryotic 16S rRNA genes (16S rRNA) was performed on 20 lone star ticks (including males, females, and nymphs). Pyrosequencing of the rickettsial sca0 gene (also known as ompA or rompA) was performed on six ticks. Female ticks had less diverse microbiomes than males and nymphs, with greater population densities of Rickettsiales. The most common members of Rickettsiales were "Candidatus Rickettsia amblyommii" and "Candidatus Midichloria mitochondrii." "Ca. Rickettsia amblyommii" was 2.6-fold more common in females than males, and there was no sequence diversity in the sca0 gene. These results are consistent with a predominantly vertical transmission pattern for "Ca. Rickettsia amblyommii."

Pooled deep sequencing of Plasmodium falciparum isolates: an efficient and scalable tool to quantify prevailing malaria drug-resistance genotypes.

Molecular surveillance for drug-resistant malaria parasites requires reliable, timely, and scalable methods. These data may be efficiently produced by genotyping parasite populations using second-generation sequencing (SGS). We designed and validated a SGS protocol to quantify mutant allele frequencies in the Plasmodium falciparum genes dhfr and dhps in mixed isolates. We applied this new protocol to field isolates from children and compared it to standard genotyping using Sanger sequencing. The SGS protocol accurately quantified dhfr and dhps allele frequencies in a mixture of parasite strains. Using SGS of DNA that was extracted and then pooled from individual isolates, we estimated mutant allele frequencies that were closely correlated to those estimated by Sanger sequencing (correlations, >0.98). The SGS protocol obviated most molecular steps in conventional methods and is cost saving for parasite populations >50. This SGS genotyping method efficiently and reproducibly estimates parasite allele frequencies within populations of P. falciparum for molecular epidemiologic studies.

Hypnozoite depletion in successive Plasmodium vivax relapses.

Genotyping Plasmodium vivax relapses can provide insights into hypnozoite biology. We performed targeted amplicon sequencing of 127 relapses occurring in Indonesian soldiers returning to malaria-free Java after yearlong deployment in malarious Eastern Indonesia. Hepatic carriage of multiple hypnozoite clones was evident in three-quarters of soldiers with two successive relapses, yet the majority of relapse episodes only displayed one clonal population. The number of clones detected in relapse episodes decreased over time and through successive relapses, especially in individuals who received hypnozoiticidal therapy. Interrogating the multiplicity of infection in this P. vivax relapse cohort reveals evidence of independent activation and slow depletion of hypnozoites over many months by multiple possible mechanisms, including parasite senescence and host immunity.