Evolution of Partial Resistance to Artemisinins in Malaria Parasites in Uganda.

BACKGROUND: Partial resistance of Plasmodium falciparum to the artemisinin component of artemisinin-based combination therapies, the most important malaria drugs, emerged in Southeast Asia and now threatens East Africa. Partial resistance, which manifests as delayed clearance after therapy, is mediated principally by mutations in the kelch protein K13 (PfK13). Limited longitudinal data are available on the emergence and spread of artemisinin resistance in Africa. METHODS: We performed annual surveillance among patients who presented with uncomplicated malaria at 10 to 16 sites across Uganda from 2016 through 2022. We sequenced the gene encoding kelch 13 (pfk13) and analyzed relatedness using molecular methods. We assessed malaria metrics longitudinally in eight Ugandan districts from 2014 through 2021. RESULTS: By 2021-2022, the prevalence of parasites with validated or candidate resistance markers reached more than 20% in 11 of the 16 districts where surveillance was conducted. The PfK13 469Y and 675V mutations were seen in far northern Uganda in 2016-2017 and increased and spread thereafter, reaching a combined prevalence of 10 to 54% across much of northern Uganda, with spread to other regions. The 469F mutation reached a prevalence of 38 to 40% in one district in southwestern Uganda in 2021-2022. The 561H mutation, previously described in Rwanda, was first seen in southwestern Uganda in 2021, reaching a prevalence of 23% by 2022. The 441L mutation reached a prevalence of 12 to 23% in three districts in western Uganda in 2022. Genetic analysis indicated local emergence of mutant parasites independent of those in Southeast Asia. The emergence of resistance was observed predominantly in areas where effective malaria control had been discontinued or transmission was unstable. CONCLUSIONS: Data from Uganda showed the emergence of partial resistance to artemisinins in multiple geographic locations, with increasing prevalence and regional spread over time. (Funded by the National Institutes of Health.).

Malaria Species Positivity Rates Among Symptomatic Individuals Across Regions of Differing Transmission Intensities in Mainland Tanzania.

BACKGROUND: Recent data indicate that non-Plasmodium falciparum species may be more prevalent than thought in sub-Saharan Africa. Although Plasmodium malariae, Plasmodium ovale spp., and Plasmodium vivax are less severe than P. falciparum, treatment and control are more challenging, and their geographic distributions are not well characterized. METHODS: We randomly selected 3284 of 12 845 samples collected from cross-sectional surveys in 100 health facilities across 10 regions of Mainland Tanzania and performed quantitative real-time PCR to determine presence and parasitemia of each malaria species. RESULTS: P. falciparum was most prevalent, but P. malariae and P. ovale were found in all but 1 region, with high levels (>5%) of P. ovale in 7 regions. The highest P. malariae positivity rate was 4.5% in Mara and 8 regions had positivity rates ≥1%. We only detected 3 P. vivax infections, all in Kilimanjaro. While most nonfalciparum malaria-positive samples were coinfected with P. falciparum, 23.6% (n = 13 of 55) of P. malariae and 14.7% (n = 24 of 163) of P. ovale spp. were monoinfections. CONCLUSIONS: P. falciparum remains by far the largest threat, but our data indicate that malaria elimination efforts in Tanzania will require increased surveillance and improved understanding of the biology of nonfalciparum species.

Varied prevalence of antimalarial drug resistance markers in different populations of newly arrived refugees in Uganda.

Newly arrived refugees offer insights into malaria epidemiology in their countries of origin. We evaluated asymptomatic refugee children within 7 days of arrival in Uganda from South Sudan and the Democratic Republic of Congo (DRC) in 2022 for parasitemia, parasite species, and Plasmodium falciparum drug resistance markers. Asymptomatic P. falciparum infections were common in both populations. Co-infection with P. malariae was more common in DRC refugees. Prevalences of markers of aminoquinoline resistance (PfCRT K76T, PfMDR1 N86Y) were much higher in South Sudan refugees, of antifolate resistance (PfDHFR C59R and I164L, PfDHPS A437G and K540E) much higher in DRC refugees, and of artemisinin partial resistance (ART-R; PfK13 C469Y and A675V) moderate in both populations. Prevalences of most mutations differed from those seen in Ugandans attending health centers near the refugee centers. Refugee evaluations yielded insights into varied malaria epidemiology and identified markers of ART-R in two previously little-studied countries.

Ex vivo drug susceptibility and resistance mediating genetic polymorphisms of Plasmodium falciparum in Bobo-Dioulasso, Burkina Faso.

Malaria remains a leading cause of morbidity and mortality in Burkina Faso, which utilizes artemether-lumefantrine as the principal therapy to treat uncomplicated malaria and seasonal malaria chemoprevention with monthly sulfadoxine-pyrimethamine plus amodiaquine in children during the transmission season. Monitoring the activities of available antimalarial drugs is a high priority. We assessed the ex vivo susceptibility of Plasmodium falciparum to 11 drugs in isolates from patients presenting with uncomplicated malaria in Bobo-Dioulasso in 2021 and 2022. IC50 values were derived using a standard 72 h growth inhibition assay. Parasite DNA was sequenced to characterize known drug resistance-mediating polymorphisms. Isolates were generally susceptible, with IC50 values in the low-nM range, to chloroquine (median IC5010 nM, IQR 7.9-24), monodesethylamodiaquine (22, 14-46) piperaquine (6.1, 3.6-9.2), pyronaridine (3.0, 1.3-5.5), quinine (50, 30-75), mefloquine (7.1, 3.7-10), lumefantrine (7.1, 4.5-12), dihydroartemisinin (3.7, 2.2-5.5), and atovaquone (0.2, 0.1-0.3) and mostly resistant to cycloguanil (850, 543-1,290) and pyrimethamine (33,200, 18,400-54,200), although a small number of outliers were seen. Considering genetic markers of resistance to aminoquinolines, most samples had wild-type PfCRT K76T (87%) and PfMDR1 N86Y (95%) sequences. For markers of resistance to antifolates, established PfDHFR and PfDHPS mutations were highly prevalent, the PfDHPS A613S mutation was seen in 19% of samples, and key markers of high-level resistance (PfDHFR I164L; PfDHPS K540E) were absent or rare (A581G). Mutations in the PfK13 propeller domain known to mediate artemisinin partial resistance were not detected. Overall, our results suggest excellent susceptibilities to drugs now used to treat malaria and moderate, but stable, resistance to antifolates used to prevent malaria.

coiaf: Directly estimating complexity of infection with allele frequencies.

In malaria, individuals are often infected with different parasite strains. The complexity of infection (COI) is defined as the number of genetically distinct parasite strains in an individual. Changes in the mean COI in a population have been shown to be informative of changes in transmission intensity with a number of probabilistic likelihood and Bayesian models now developed to estimate the COI. However, rapid, direct measures based on heterozygosity or FwS do not properly represent the COI. In this work, we present two new methods that use easily calculated measures to directly estimate the COI from allele frequency data. Using a simulation framework, we show that our methods are computationally efficient and comparably accurate to current approaches in the literature. Through a sensitivity analysis, we characterize how the distribution of parasite densities, the assumed sequencing depth, and the number of sampled loci impact the bias and accuracy of our two methods. Using our developed methods, we further estimate the COI globally from Plasmodium falciparum sequencing data and compare the results against the literature. We show significant differences in the estimated COI globally between continents and a weak relationship between malaria prevalence and COI.

Expansion of artemisinin partial resistance mutations and lack of histidine rich protein-2 and -3 deletions in Plasmodium falciparum infections from Rukara, Rwanda.

BACKGROUND: Emerging artemisinin partial resistance and diagnostic resistance are a threat to malaria control in Africa. Plasmodium falciparum kelch13 (k13) propeller-domain mutations that confer artemisinin partial resistance have emerged in Africa. k13-561H was initially described at a frequency of 7.4% from Masaka in 2014-2015, but not present in nearby Rukara. By 2018, 19.6% of isolates in Masaka and 22% of isolates in Rukara contained the mutation. Longitudinal monitoring is essential to inform control efforts. In Rukara, an assessment was conducted to evaluate recent k13-561H prevalence changes, as well as other key mutations. Prevalence of hrp2/3 deletions was also assessed. METHODS: Samples collected in Rukara in 2021 were genotyped for key artemisinin and partner drug resistance mutations using molecular inversion probe assays and for hrp2/3 deletions using qPCR. RESULTS: Clinically validated k13 artemisinin partial resistance mutations continue to increase in prevalence with the overall level of mutant infections reaching 32% in Rwanda. The increase appears to be due to the rapid emergence of k13-675V (6.4%, 6/94 infections), previously not observed, rather than continued expansion of 561H (23.5% 20/85). Mutations to partner drugs and other anti-malarials were variable, with high levels of multidrug resistance 1 (mdr1) N86 (95.5%) associated with lumefantrine decreased susceptibility and dihydrofolate reductase (dhfr) 164L (24.7%) associated with a high level of antifolate resistance, but low levels of amodiaquine resistance polymorphisms with chloroquine resistance transporter (crt) 76T: at 6.1% prevalence. No hrp2 or hrp3 gene deletions associated with diagnostic resistance were found. CONCLUSIONS: Increasing prevalence of artemisinin partial resistance due to k13-561H and the rapid expansion of k13-675V is concerning for the longevity of artemisinin effectiveness in the region. False negative RDT results do not appear to be an issue with no hrp2 or hpr3 deletions detected. Continued molecular surveillance in this region and surrounding areas is needed to follow artemisinin partial resistance and provide early detection of partner drug resistance, which would likely compromise control and increase malaria morbidity and mortality in East Africa.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania.

BACKGROUND: Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania. METHODS: Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites. RESULTS: Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic. CONCLUSION: Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity.

BACKGROUND: In 2021 and 2023, the World Health Organization approved RTS,S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (PfCSP), but polymorphisms in the gene raise concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission intensities in Mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country. METHODS: The analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of Mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright's inbreeding coefficient (Fws), Wright's fixation index (FST), principal component analysis, nucleotide diversity, and Tajima's D were used to assess within-host parasite diversity, population structure and natural selection. RESULTS: Based on Fws (< 0.95), there was high polyclonality (ranging from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST = 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima's D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences revealed 50 different haplotypes (H_1 to H_50), with only 2% of sequences matching the 3D7 strain haplotype (H_50). Conversely, with the NF54 strain, the Pfcsp C-terminal sequences revealed 49 different haplotypes (H_1 to H_49), with only 0.4% of the sequences matching the NF54 strain (Hap_49). CONCLUSIONS: The findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima's D values, consistent with balancing selection for variants within Th2R and Th3R regions. The study observed differences between the intended haplotypes incorporated into the design of RTS,S and R21 vaccines and those present in natural parasite populations. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

SARS-CoV-2 Serosurveys: How Antigen, Isotype and Threshold Choices Affect the Outcome.

BACKGROUND: Evaluating the performance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) serological assays and clearly articulating the utility of selected antigens, isotypes, and thresholds is crucial to understanding the prevalence of infection within selected communities. METHODS: This cross-sectional study, implemented in 2020, screened PCRconfirmed coronavirus disease 2019 patients (n 86), banked prepandemic and negative samples (n 96), healthcare workers and family members (n 552), and university employees (n 327) for antiSARS-CoV-2 receptor-binding domain, trimeric spike protein, and nucleocapsid protein immunoglobulin (Ig)G and IgA antibodies with a laboratory-developed enzyme-linked immunosorbent assay and tested how antigen, isotype and threshold choices affected the seroprevalence outcomes. The following threshold methods were evaluated: (i) mean 3 standard deviations of the negative controls; (ii) 100 specificity for each antigen-isotype combination; and (iii) the maximal Youden index. RESULTS: We found vastly different seroprevalence estimates depending on selected antigens and isotypes and the applied threshold method, ranging from 0.0 to 85.4. Subsequently, we maximized specificity and reported a seroprevalence, based on more than one antigen, ranging from 9.3 to 25.9. CONCLUSIONS: This study revealed the importance of evaluating serosurvey tools for antigen-, isotype-, and threshold-specific sensitivity and specificity, to interpret qualitative serosurvey outcomes reliably and consistently across studies.

Expression Microdissection for the Analysis of miRNA in a Single-Cell Type.

Cell-specific microRNA (miRNA) expression estimates are important in characterizing the localization of miRNA signaling within tissues. Much of these data are obtained from cultured cells, a process known to significantly alter miRNA expression levels. Thus, our knowledge of in vivo cell miRNA expression estimates is poor. We previously demonstrated expression microdissection-miRNA-sequencing (xMD-miRNA-seq) to acquire in vivo estimates, directly from formalin-fixed tissues, albeit with a limited yield. In this study, we optimized each step of the xMD process, including tissue retrieval, tissue transfer, film preparation, and RNA isolation, to increase RNA yields and ultimately show strong enrichment for in vivo miRNA expression by qPCR array. These method improvements, such as the development of a noncrosslinked ethylene vinyl acetate membrane, resulted in a 23- to 45-fold increase in miRNA yield, depending on the cell type. By qPCR, miR-200a increased by 14-fold in xMD-derived small intestine epithelial cells, with a concurrent 336-fold reduction in miR-143 relative to the matched nondissected duodenal tissue. xMD is now an optimized method to obtain robust in vivo miRNA expression estimates from cells. xMD will allow formalin-fixed tissues from surgical pathology archives to make theragnostic biomarker discoveries.

An optimized GATK4 pipeline for Plasmodium falciparum whole genome sequencing variant calling and analysis.

BACKGROUND: Accurate variant calls from whole genome sequencing (WGS) of Plasmodium falciparum infections are crucial in malaria population genomics. Here a falciparum variant calling pipeline based on GATK version 4 (GATK4) was optimized and applied to 6626 public Illumina WGS samples. METHODS: Control WGS and accurate PacBio assemblies of 10 laboratory strains were leveraged to optimize parameters that control the heterozygosity, local assembly region size, ploidy, mapping and base quality in both GATK HaplotypeCaller and GenotypeGVCFs. From these controls, a high-quality training dataset was generated to recalibrate the raw variant data. RESULTS: On current high-quality samples (read length = 250 bp, insert size = 405-524 bp), the optimized pipeline shows improved sensitivity (86.6 ± 1.7% for SNPs and 82.2 ± 5.9% for indels) compared to the default GATK4 pipeline (77.7 ± 1.3% for SNPs; and 73.1 ± 5.1% for indels, adjusted P < 0.001) and previous variant calling with GATK version 3 (GATK3, 70.3 ± 3.0% for SNPs and 59.7 ± 5.8% for indels, adjusted P < 0.001). Its sensitivity on simulated mixed infection samples (80.8 ± 6.1% for SNPs and 78.3 ± 5.1% for indels) was again improved relative to default GATK4 (68.8 ± 6.0% for SNPs and 38.9 ± 0.7% for indels, adjusted, adjusted P < 0.001). Precision was high and comparable across all pipelines on each type of data tested. The resulting combination of high-quality SNPs and indels increases the resolution of local population population structure detection in sub-Saharan Africa. Finally, increasing ploidy improves the detection of drug resistance mutations and estimation of complexity of infection. CONCLUSIONS: Overall, this study provides an optimized falciparum GATK4 pipeline resource for variant calling which should help improve genomic studies of malaria.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia.

BACKGROUND: Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions. METHODS: In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses. RESULTS: The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure. CONCLUSION: Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

Background malaria incidence and parasitemia during the three-dose RTS,S/AS01 vaccination series do not reduce magnitude of antibody response nor efficacy against the first case of malaria.

BACKGROUND: RTS,S/AS01 has been recommended by WHO for widespread implementation in medium to high malaria transmission settings. Previous analyses have noted lower vaccine efficacies in higher transmission settings, possibly due to the more rapid development of naturally acquired immunity in the control group. METHODS: To investigate a reduced immune response to vaccination as a potential mechanism behind lower efficacy in high transmission areas, we examine initial vaccine antibody (anti-CSP IgG) response and vaccine efficacy against the first case of malaria (to exclude the effect of naturally acquired immunity) using data from three study areas (Kintampo, Ghana; Lilongwe, Malawi; Lambaréné, Gabon) from the 2009-2014 phase III trial (NCT00866619). Our key exposures are parasitemia during the vaccination series and background malaria incidence. We calculate vaccine efficacy (one minus hazard ratio) using a cox-proportional hazards model and allowing for the time-varying effect of RTS,S/AS01. RESULTS: We find that antibody responses to the primary three-dose vaccination series were higher in Ghana than in Malawi and Gabon, but that neither antibody levels nor vaccine efficacy against the first case of malaria varied by background incidence or parasitemia during the primary vaccination series. CONCLUSIONS: We find that vaccine efficacy is unrelated to infections during vaccination. Contributing to a conflicting literature, our results suggest that vaccine efficacy is also unrelated to infections before vaccination, meaning that control-group immunity is likely a major reason for lower efficacy in high transmission settings, not reduced immune responses to RTS,S/AS01. This may be reassuring for implementation in high transmission settings, though further studies are needed.

Distinctive Kaposi Sarcoma-Associated Herpesvirus Serological Profile during Acute Plasmodium falciparum Malaria Episodes.

The seroprevalence of Kaposi sarcoma-associated herpesvirus (KSHV) and the incidence of endemic Kaposi sarcoma (KS) overlap with regions of malaria endemicity in sub-Saharan Africa. Multiple studies have shown an increased risk of KSHV seroconversion in children from high malaria compared to low malaria regions; however, the impact of acute episodes of Plasmodium falciparum (P. falciparum) malaria on KSHV's biphasic life cycle and lytic reactivation has not been determined. Here, we examined KSHV serological profiles and viral loads in 134 children with acute malaria and 221 healthy children from high malaria regions in Kisumu, as well as 77 healthy children from low malaria regions in Nandi. We assayed KSHV, Epstein-Barr virus (EBV), and P. falciparum malaria antibody responses in these three by multiplexed Luminex assay. We confirmed that KSHV seroprevalence was significantly associated with malaria endemicity (OR = 1.95, 1.18-3.24 95% CI, p = 0.01) with 71-77% seropositivity in high-malaria (Kisumu) compared to 28% in low-malaria (Nandi) regions. Furthermore, KSHV serological profiles during acute malaria episodes were distinct from age-matched non-malaria-infected children from the same region. Paired IgG levels also varied after malaria treatment, with significantly higher anti-ORF59 at day 0 but elevated ORF38, ORF73, and K8.1 at day 3. Acute malaria episodes is characterized by perturbation of KSHV latency in seropositive children, providing further evidence that malaria endemicity contributes to the observed increase in endemic KS incidence in sub-Saharan Africa.

Sensitive, Highly Multiplexed Sequencing of Microhaplotypes From the Plasmodium falciparum Heterozygome.

BACKGROUND: Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited. METHODS: A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls. RESULTS: Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes. CONCLUSIONS: The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Amplicon Sequencing as a Potential Surveillance Tool for Complexity of Infection and Drug Resistance Markers in Plasmodium falciparum Asymptomatic Infections.

BACKGROUND: Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers. METHODS: We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker. RESULTS: We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65-2.77] vs 2.22 [95% CI, 2.17-2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased. CONCLUSIONS: This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Malaria Transmission Intensity Likely Modifies RTS, S/AS01 Efficacy Due to a Rebound Effect in Ghana, Malawi, and Gabon.

BACKGROUND: RTS,S/AS01 is the first malaria vaccine to be approved and recommended for widespread implementation by the World Health Organization (WHO). Trials reported lower vaccine efficacies in higher-incidence sites, potentially due to a "rebound" in malaria cases in vaccinated children. When naturally acquired protection in the control group rises and vaccine protection in the vaccinated wanes concurrently, malaria incidence can become greater in the vaccinated than in the control group, resulting in negative vaccine efficacies. METHODS: Using data from the 2009-2014 phase III trial (NCT00866619) in Lilongwe, Malawi; Kintampo, Ghana; and Lambaréné, Gabon, we evaluate this hypothesis by estimating malaria incidence in each vaccine group over time and in varying transmission settings. After estimating transmission intensities using ecological variables, we fit models with 3-way interactions between vaccination, time, and transmission intensity. RESULTS: Over time, incidence decreased in the control group and increased in the vaccine group. Three-dose efficacy in the lowest-transmission-intensity group (0.25 cases per person-year [CPPY]) decreased from 88.2% to 15.0% over 4.5 years, compared with 81.6% to -27.7% in the highest-transmission-intensity group (3 CPPY). CONCLUSIONS: These findings suggest that interventions, including the fourth RTS,S dose, that protect vaccinated individuals during the potential rebound period should be implemented for high-transmission settings.

Decreased Susceptibility to Dihydrofolate Reductase Inhibitors Associated With Genetic Polymorphisms in Ugandan Plasmodium falciparum Isolates.

BACKGROUND: The Plasmodium falciparum dihydrofolate reductase (PfDHFR) inhibitors pyrimethamine and cycloguanil (the active metabolite of proguanil) have important roles in malaria chemoprevention, but drug resistance challenges their efficacies. A new compound, P218, was designed to overcome resistance, but drug-susceptibility data for P falciparum field isolates are limited. METHODS: We studied ex vivo PfDHFR inhibitor susceptibilities of 559 isolates from Tororo and Busia districts, Uganda, from 2016 to 2020, sequenced 383 isolates, and assessed associations between genotypes and drug-susceptibility phenotypes. RESULTS: Median half-maximal inhibitory concentrations (IC50s) were 42 100 nM for pyrimethamine, 1200 nM for cycloguanil, 13000 nM for proguanil, and 0.6 nM for P218. Among sequenced isolates, 3 PfDHFR mutations, 51I (100%), 59R (93.7%), and 108N (100%), were very common, as previously seen in Uganda, and another mutation, 164L (12.8%), had moderate prevalence. Increasing numbers of mutations were associated with decreasing susceptibility to pyrimethamine, cycloguanil, and P218, but not proguanil, which does not act directly against PfDHFR. Differences in P218 susceptibilities were modest, with median IC50s of 1.4 nM for parasites with mixed genotype at position 164 and 5.7 nM for pure quadruple mutant (51I/59R/108N/164L) parasites. CONCLUSIONS: Resistance-mediating PfDHFR mutations were common in Ugandan isolates, but P218 retained excellent activity against mutant parasites.

Whole Blood Transfusion for Severe Malarial Anemia in a High Plasmodium falciparum Transmission Setting.

BACKGROUND: Severe malaria resulting from Plasmodium falciparum infection is the leading parasitic cause of death in children worldwide, and severe malarial anemia (SMA) is the most common clinical presentation. The evidence in support of current blood transfusion guidelines for patients with SMA is limited. METHODS: We conducted a retrospective cohort study of 911 hospitalized children with SMA in a holoendemic region of Zambia to examine the association of whole blood transfusion with in-hospital survival. Data were analyzed in adjusted logistic regression models using multiple imputation for missing data. RESULTS: The median age of patients was 24 months (interquartile range, 16-30) and overall case fatality was 16%. Blood transfusion was associated with 35% reduced odds of death in children with SMA (odds ratio, 0.65; 95% confidence interval, .52-.81; P = .0002) corresponding to a number-needed-to-treat (NNT) of 14 patients. Children with SMA complicated by thrombocytopenia were more likely to benefit from transfusion than those without thrombocytopenia (NNT = 5). Longer storage time of whole blood was negatively associated with survival and with the posttransfusion rise in the platelet count but was not associated with the posttransfusion change in hemoglobin concentration. CONCLUSIONS: Whole blood given to pediatric patients with SMA was associated with improved survival, mainly among those with thrombocytopenia who received whole blood stored for <4 weeks. These findings point to a potential use for incorporating thrombocytopenia into clinical decision making and management of severe malaria, which can be further assessed in prospective studies, and underline the importance of maintaining reliable blood donation networks in areas of high malaria transmission.

mirTarRnaSeq: An R/Bioconductor Statistical Package for miRNA-mRNA Target Identification and Interaction Analysis.

We introduce mirTarRnaSeq, an R/Bioconductor package for quantitative assessment of miRNA-mRNA relationships within sample cohorts. mirTarRnaSeq is a statistical package to explore predicted or pre-hypothesized miRNA-mRNA relationships following target prediction.We present two use cases applying mirTarRnaSeq. First, to identify miRNA targets, we examined EBV miRNAs for interaction with human and virus transcriptomes of stomach adenocarcinoma. This revealed enrichment of mRNA targets highly expressed in CD105+ endothelial cells, monocytes, CD4+ T cells, NK cells, CD19+ B cells, and CD34 cells. Next, to investigate miRNA-mRNA relationships in SARS-CoV-2 (COVID-19) infection across time, we used paired miRNA and RNA sequenced datasets of SARS-CoV-2 infected lung epithelial cells across three time points (4, 12, and 24 hours post-infection). mirTarRnaSeq identified evidence for human miRNAs targeting cytokine signaling and neutrophil regulation immune pathways from 4 to 24 hours after SARS-CoV-2 infection. Confirming the clinical relevance of these predictions, three of the immune specific mRNA-miRNA relationships identified in human lung epithelial cells after SARS-CoV-2 infection were also observed to be differentially expressed in blood from patients with COVID-19. Overall, mirTarRnaSeq is a robust tool that can address a wide-range of biological questions providing improved prediction of miRNA-mRNA interactions.