Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1.

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Cannabinoid enhancement of lncRNA MMP25-AS1/MMP25 interaction reduces neutrophil infiltration and intestinal epithelial injury in HIV/SIV infection.

Intestinal epithelial barrier dysfunction, a hallmark of HIV/SIV infection, persists despite viral suppression by combination antiretroviral therapy (cART). Emerging evidence suggests a critical role for long noncoding RNAs (lncRNAs) in maintaining epithelial homeostasis. We simultaneously profiled lncRNA/mRNA expression exclusively in colonic epithelium (CE) of SIV-infected rhesus macaques (RMs) administered vehicle (VEH) or Δ-9-tetrahydrocannabinol (THC). Relative to controls, fewer lncRNAs were up- or downregulated in CE of THC/SIV compared with VEH/SIV RMs. Importantly, reciprocal expression of the natural antisense lncRNA MMP25-AS1 (up 2.3-fold) and its associated protein-coding gene MMP25 (attracts neutrophils by inactivating alpha-1 anti-trypsin/SERPINA1) (down 2.2-fold) was detected in CE of THC/SIV RMs. Computational analysis verified 2 perfectly matched complementary regions and an energetically stable (normalized binding free energy = -0.2626) MMP25-AS1/MMP25 duplex structure. MMP25-AS1 overexpression blocked IFN-γ-induced MMP25 mRNA and protein expression in vitro. Elevated MMP25 protein expression in CE of VEH/SIV but not THC/SIV RMs was associated with increased infiltration by myeloperoxidase/CD11b++ neutrophils (transendothelial migration) and epithelial CD47 (transepithelial migration) expression. Interestingly, THC administered in combination with cART increased MMP25-AS1 and reduced MMP25 mRNA/protein expression in jejunal epithelium of SIV-infected RMs. Our findings demonstrate that MMP25-AS1 is a potentially unique epigenetic regulator of MMP25 and that low-dose THC can reduce neutrophil infiltration and intestinal epithelial injury potentially by downregulating MMP25 expression through modulation of MMP25-AS1.

Alterations in Abundance and Compartmentalization of miRNAs in Blood Plasma Extracellular Vesicles and Extracellular Condensates during HIV/SIV Infection and Its Modulation by Antiretroviral Therapy (ART) and Delta-9-Tetrahydrocannabinol (Δ9-THC).

In this follow-up study, we investigated the abundance and compartmentalization of blood plasma extracellular miRNA (exmiRNA) into lipid-based carriers-blood plasma extracellular vesicles (EVs) and non-lipid-based carriers-extracellular condensates (ECs) during SIV infection. We also assessed how combination antiretroviral therapy (cART), administered in conjunction with phytocannabinoid delta-9-tetrahydrocannabinol (THC), altered the abundance and compartmentalization of exmiRNAs in the EVs and ECs of SIV-infected rhesus macaques (RMs). Unlike cellular miRNAs, exmiRNAs in blood plasma may serve as minimally invasive disease indicators because they are readily detected in stable forms. The stability of exmiRNAs in cell culture fluids and body fluids (urine, saliva, tears, cerebrospinal fluid (CSF), semen, blood) is based on their association with different carriers (lipoproteins, EVs, and ECs) that protect them from the activities of endogenous RNases. Here, we showed that in the blood plasma of uninfected control RMs, significantly less exmiRNAs were associated with EVs compared to the level (30% higher) associated with ECs, and that SIV infection altered the profile of EVs and ECs miRNAome (Manuscript 1). In people living with HIV (PLWH), host-encoded miRNAs regulate both host and viral gene expression, which may serve as indicators of disease or treatment biomarkers. The profile of miRNAs in blood plasma of PLWH (elite controllers versus viremic patients) are different, indicating that HIV may alter host miRNAome. However, there are no studies assessing the effect of cART or other substances used by PLWH, such as THC, on the abundance of exmiRNA and their association with EVs and ECs. Moreover, longitudinal exmiRNA profiles following SIV infection, treatment with THC, cART, or THC+cART remains unclear. Here, we serially analyzed miRNAs associated with blood plasma derived EVs and ECs. Methods: Paired EVs and ECs were separated from EDTA blood plasma of male Indian rhesus macaques (RMs) in five treatment groups, including VEH/SIV, VEH/SIV/cART, THC/SIV, THC/SIV/cART, or THC alone. Separation of EVs and ECs was achieved with the unparalleled nano-particle purification tool ─PPLC, a state-of-the-art, innovative technology equipped with gradient agarose bead sizes and a fast fraction collector that allows high resolution separation and retrieval of preparative quantities of sub-populations of extracellular structures. Global miRNA profiles of the paired EVs and ECs were determined with RealSeq Biosciences (Santa Cruz, CA) custom sequencing platform by conducting small RNA (sRNA)-seq. The sRNA-seq data were analyzed using various bioinformatic tools. Validation of key exmiRNA was performed using specific TaqMan microRNA stem-loop RT-qPCR assays. Results: We investigated the effect of cART, THC, or both cART and THC together on the abundance and compartmentalization of blood plasma exmiRNA in EVs and ECs in SIV-infected RMs. As shown in Manuscript 1 of this series, were in uninfected RMs, ~30% of exmiRNAs were associated with ECs, we confirmed in this follow up manuscript that exmiRNAs were present in both lipid-based carriers-EVs and non-lipid-based carriers-ECs, with 29.5 to 35.6% and 64.2 to 70.5 % being associated with EVs and ECs, respectively. Remarkably, the different treatments (cART, THC) have distinct effects on the enrichment and compartmentalization pattern of exmiRNAs. In the VEH/SIV/cART group, 12 EV-associated and 15 EC-associated miRNAs were significantly downregulated. EV-associated miR-206, a muscle-specific miRNA that is present in blood, was higher in the VEH/SIV/ART compared to the VEH/SIV group. ExmiR-139-5p that was implicated in endocrine resistance, focal adhesion, lipid and atherosclerosis, apoptosis, and breast cancer by miRNA-target enrichment analysis was significantly lower in VEH/SIV/cART compared to VEH/SIV, irrespective of the compartment. With respect to THC treatment, 5 EV-associated and 21 EC-associated miRNAs were significantly lower in the VEH/THC/SIV. EV-associated miR-99a-5p was higher in VEH/THC/SIV compared to VEH/SIV, while miR-335-5p counts were significantly lower in both EVs and ECs of THC/SIV compared to VEH/SIV. EVs from SIV/cART/THC combined treatment group have significant increases in the count of eight (miR-186-5p, miR-382-5p, miR-139-5p and miR-652, miR-10a-5p, miR-657, miR-140-5p, miR-29c-3p) miRNAs, all of which were lower in VEH/SIV/cART group. Analysis of miRNA-target enrichment showed that this set of eight miRNAs were implicated in endocrine resistance, focal adhesions, lipid and atherosclerosis, apoptosis, and breast cancer as well as cocaine and amphetamine addiction. In ECs and EVs, combined THC and cART treatment significantly increased miR-139-5p counts compared to VEH/SIV group. Significant alterations in these host miRNAs in both EVs and ECs in the untreated and treated (cART, THC, or both) RMs indicate the persistence of host responses to infection or treatments, and this is despite cART suppression of viral load and THC suppression of inflammation. To gain further insight into the pattern of miRNA alterations in EVs and ECs and to assess potential cause-and-effect relationships, we performed longitudinal miRNA profile analysis, measured in terms of months (1 and 5) post-infection (MPI). We uncovered miRNA signatures associated with THC or cART treatment of SIV-infected macaques in both EVs and ECs. While the number of miRNAs was significantly higher in ECs relative to EVs for all groups (VEH/SIV, SIV/cART, THC/SIV, THC/SIV/cART, and THC) longitudinally from 1 MPI to 5 MPI, treatment with cART and THC have longitudinal effects on the abundance and compartmentalization pattern of exmiRNAs in the two carriers. As shown in Manuscript 1 where SIV infection led to longitudinal suppression of EV-associated miRNA-128-3p, administration of cART to SIV-infected RMs did not increase miR-128-3p but resulted in longitudinal increases in six EV-associated miRNAs (miR-484, miR-107, miR-206, miR-184, miR-1260b, miR-6132). Furthermore, administration of cART to THC treated SIV-infected RMs resulted in a longitudinal decrease in three EV-associated miRNAs (miR-342-3p, miR-100-5p, miR181b-5p) and a longitudinal increase in three EC-associated miRNAs (miR-676-3p, miR-574-3p, miR-505-5p). The longitudinally altered miRNAs in SIV-infected RMs may indicate disease progression, while in the cART Group and the THC Group, the longitudinally altered miRNAs may serve as biomarkers of response to treatment. Conclusions: This paired EVs and ECs miRNAome analyses provided a comprehensive cross-sectional and longitudinal summary of the host exmiRNA responses to SIV infection and the impact of THC, cART, or THC and cART together on the miRNAome during SIV infection. Overall, our data point to previously unrecognized alterations in the exmiRNA profile in blood plasma following SIV infection. Our data also indicate that cART and THC treatment independently and in combination may alter both the abundance and the compartmentalization of several exmiRNA related to various disease and biological processes.

SIV Infection Regulates Compartmentalization of Circulating Blood Plasma miRNAs within Extracellular Vesicles (EVs) and Extracellular Condensates (ECs) and Decreases EV-Associated miRNA-128.

Background: This is Manuscript 1 of a two-part Manuscript of the same series. Here, we present findings from our first set of studies on the abundance and compartmentalization of blood plasma extracellular microRNAs (exmiRNAs) into extracellular particles, including blood plasma extracellular vesicles (EVs) and extracellular condensates (ECs) in the setting of untreated HIV/SIV infection. The goals of the study presented in this Manuscript 1 are to (i) assess the abundance and compartmentalization of exmiRNAs in EVs versus ECs in the healthy uninfected state, and (ii) evaluate how SIV infection may affect exmiRNA abundance and compartmentalization in these particles. Considerable effort has been devoted to studying the epigenetic control of viral infection, particularly in understanding the role of exmiRNAs as key regulators of viral pathogenesis. MicroRNA (miRNAs) are small (~20-22 nts) non-coding RNAs that regulate cellular processes through targeted mRNA degradation and/or repression of protein translation. Originally associated with the cellular microenvironment, circulating miRNAs are now known to be present in various extracellular environments, including blood serum and plasma. While in circulation, miRNAs are protected from degradation by ribonucleases through their association with lipid and protein carriers, such as lipoproteins and other extracellular particles-EVs and ECs. Functionally, miRNAs play important roles in diverse biological processes and diseases (cell proliferation, differentiation, apoptosis, stress responses, inflammation, cardiovascular diseases, cancer, aging, neurological diseases, and HIV/SIV pathogenesis). While lipoproteins and EV-associated exmiRNAs have been characterized and linked to various disease processes, the association of exmiRNAs with ECs is yet to be made. Likewise, the effect of SIV infection on the abundance and compartmentalization of exmiRNAs within extracellular particles is unclear. Literature in the EV field has suggested that most circulating miRNAs may not be associated with EVs. However, a systematic analysis of the carriers of exmiRNAs has not been conducted due to the inefficient separation of EVs from other extracellular particles, including ECs. Methods: Paired EVs and ECs were separated from EDTA blood plasma of SIV-uninfected male Indian rhesus macaques (RMs, n = 15). Additionally, paired EVs and ECs were isolated from EDTA blood plasma of combination anti-retroviral therapy (cART) naïve SIV-infected (SIV+, n = 3) RMs at two time points (1- and 5-months post infection, 1 MPI and 5 MPI). Separation of EVs and ECs was achieved with PPLC, a state-of-the-art, innovative technology equipped with gradient agarose bead sizes and a fast fraction collector that allows high-resolution separation and retrieval of preparative quantities of sub-populations of extracellular particles. Global miRNA profiles of the paired EVs and ECs were determined with RealSeq Biosciences (Santa Cruz, CA) custom sequencing platform by conducting small RNA (sRNA)-seq. The sRNA-seq data were analyzed using various bioinformatic tools. Validation of key exmiRNAs was performed using specific TaqMan microRNA stem-loop RT-qPCR assays. Results: We showed that exmiRNAs in blood plasma are not restricted to any type of extracellular particles but are associated with lipid-based carriers-EVs and non-lipid-based carriers-ECs, with a significant (~30%) proportion of the exmiRNAs being associated with ECs. In the blood plasma of uninfected RMs, a total of 315 miRNAs were associated with EVs, while 410 miRNAs were associated with ECs. A comparison of detectable miRNAs within paired EVs and ECs revealed 19 and 114 common miRNAs, respectively, detected in all 15 RMs. Let-7a-5p, Let-7c-5p, miR-26a-5p, miR-191-5p, and let-7f-5p were among the top 5 detectable miRNAs associated with EVs in that order. In ECs, miR-16-5p, miR-451, miR-191-5p, miR-27a-3p, and miR-27b-3p, in that order, were the top detectable miRNAs in ECs. miRNA-target enrichment analysis of the top 10 detected common EV and EC miRNAs identified MYC and TNPO1 as top target genes, respectively. Functional enrichment analysis of top EV- and EC-associated miRNAs identified common and distinct gene-network signatures associated with various biological and disease processes. Top EV-associated miRNAs were implicated in cytokine-cytokine receptor interactions, Th17 cell differentiation, IL-17 signaling, inflammatory bowel disease, and glioma. On the other hand, top EC-associated miRNAs were implicated in lipid and atherosclerosis, Th1 and Th2 cell differentiation, Th17 cell differentiation, and glioma. Interestingly, infection of RMs with SIV revealed that the brain-enriched miR-128-3p was longitudinally and significantly downregulated in EVs, but not ECs. This SIV-mediated decrease in miR-128-3p counts was validated by specific TaqMan microRNA stem-loop RT-qPCR assay. Remarkably, the observed SIV-mediated decrease in miR-128-3p levels in EVs from RMs agrees with publicly available EV miRNAome data by Kaddour et al., 2021, which showed that miR-128-3p levels were significantly lower in semen-derived EVs from HIV-infected men who used or did not use cocaine compared to HIV-uninfected individuals. These findings confirmed our previously reported finding and suggested that miR-128 may be a target of HIV/SIV. Conclusions: In the present study, we used sRNA sequencing to provide a holistic understanding of the repertoire of circulating exmiRNAs and their association with extracellular particles, such as EVs and ECs. Our data also showed that SIV infection altered the profile of the miRNAome of EVs and revealed that miR-128-3p may be a potential target of HIV/SIV. The significant decrease in miR-128-3p in HIV-infected humans and in SIV-infected RMs may indicate disease progression. Our study has important implications for the development of biomarker approaches for various types of cancer, cardiovascular diseases, organ injury, and HIV based on the capture and analysis of circulating exmiRNAs.

Cannabinoids modulate the microbiota-gut-brain axis in HIV/SIV infection by reducing neuroinflammation and dysbiosis while concurrently elevating endocannabinoid and indole-3-propionate levels.

BACKGROUND: Although the advent of combination anti-retroviral therapy (cART) has transformed HIV into a manageable chronic disease, an estimated 30-50% of people living with HIV (PLWH) exhibit cognitive and motor deficits collectively known as HIV-associated neurocognitive disorders (HAND). A key driver of HAND neuropathology is chronic neuroinflammation, where proinflammatory mediators produced by activated microglia and macrophages are thought to inflict neuronal injury and loss. Moreover, the dysregulation of the microbiota-gut-brain axis (MGBA) in PLWH, consequent to gastrointestinal dysfunction and dysbiosis, can lead to neuroinflammation and persistent cognitive impairment, which underscores the need for new interventions. METHODS: We performed RNA-seq and microRNA profiling in basal ganglia (BG), metabolomics (plasma) and shotgun metagenomic sequencing (colon contents) in uninfected and SIV-infected rhesus macaques (RMs) administered vehicle (VEH/SIV) or delta-9-tetrahydrocannabinol (THC) (THC/SIV). RESULTS: Long-term, low-dose THC reduced neuroinflammation and dysbiosis and significantly increased plasma endocannabinoid, endocannabinoid-like, glycerophospholipid and indole-3-propionate levels in chronically SIV-infected RMs. Chronic THC potently blocked the upregulation of genes associated with type-I interferon responses (NLRC5, CCL2, CXCL10, IRF1, IRF7, STAT2, BST2), excitotoxicity (SLC7A11), and enhanced protein expression of WFS1 (endoplasmic reticulum stress) and CRYM (oxidative stress) in BG. Additionally, THC successfully countered miR-142-3p-mediated suppression of WFS1 protein expression via a cannabinoid receptor-1-mediated mechanism in HCN2 neuronal cells. Most importantly, THC significantly increased the relative abundance of Firmicutes and Clostridia including indole-3-propionate (C. botulinum, C. paraputrificum, and C. cadaveris) and butyrate (C. butyricum, Faecalibacterium prausnitzii and Butyricicoccus pullicaecorum) producers in colonic contents. CONCLUSION: This study demonstrates the potential of long-term, low-dose THC to positively modulate the MGBA by reducing neuroinflammation, enhancing endocannabinoid levels and promoting the growth of gut bacterial species that produce neuroprotective metabolites, like indole-3-propionate. The findings from this study may benefit not only PLWH on cART, but also those with no access to cART and more importantly, those who fail to suppress the virus under cART.

Antigenic diversity in malaria parasites is maintained on extrachromosomal DNA.

Sequence variation among antigenic var genes enables Plasmodium falciparum malaria parasites to evade host immunity. Using long sequence reads from haploid clones from a mutation accumulation experiment, we detect var diversity inconsistent with simple chromosomal inheritance. We discover putatively circular DNA that is strongly enriched for var genes, which exist in multiple alleles per locus separated by recombination and indel events. Extrachromosomal DNA likely contributes to rapid antigenic diversification in P. falciparum.

Chronic delta-9-tetrahydrocannabinol (THC) treatment counteracts SIV-induced modulation of proinflammatory microRNA cargo in basal ganglia-derived extracellular vesicles.

BACKGROUND: Early invasion of the central nervous system (CNS) by human immunodeficiency virus (HIV) (Gray et al. in Brain Pathol 6:1-15, 1996; An et al. in Ann Neurol 40:611-6172, 1996), results in neuroinflammation, potentially through extracellular vesicles (EVs) and their micro RNAs (miRNA) cargoes (Sharma et al. in FASEB J 32:5174-5185, 2018; Hu et al. in Cell Death Dis 3:e381, 2012). Although the basal ganglia (BG) is a major target and reservoir of HIV in the CNS (Chaganti et al. in Aids 33:1843-1852, 2019; Mintzopoulos et al. in Magn Reson Med 81:2896-2904, 2019), whether BG produces EVs and the effect of HIV and/or the phytocannabinoid-delta-9-tetrahydrocannabinol (THC) on BG-EVs and HIV neuropathogenesis remain unknown. METHODS: We used the simian immunodeficiency virus (SIV) model of HIV and THC treatment in rhesus macaques (Molina et al. in AIDS Res Hum Retroviruses 27:585-592, 2011) to demonstrate for the first time that BG contains EVs (BG-EVs), and that BG-EVs cargo and function are modulated by SIV and THC. We also used primary astrocytes from the brains of wild type (WT) and CX3CR1+/GFP mice to investigate the significance of BG-EVs in CNS cells. RESULTS: Significant changes in BG-EV-associated miRNA specific to SIV infection and THC treatment were observed. BG-EVs from SIV-infected rhesus macaques (SIV EVs) contained 11 significantly downregulated miRNAs. Remarkably, intervention with THC led to significant upregulation of 37 miRNAs in BG-EVs (SIV-THC EVs). Most of these miRNAs are predicted to regulate pathways related to inflammation/immune regulation, TLR signaling, Neurotrophin TRK receptor signaling, and cell death/response. BG-EVs activated WT and CX3CR1+/GFP astrocytes and altered the expression of CD40, TNFα, MMP-2, and MMP-2 gene products in primary mouse astrocytes in an EV and CX3CR1 dependent manners. CONCLUSIONS: Our findings reveal a role for BG-EVs as a vehicle with potential to disseminate HIV- and THC-induced changes within the CNS.

A Malaria Parasite Cross Reveals Genetic Determinants of Plasmodium falciparum Growth in Different Culture Media.

What genes determine in vitro growth and nutrient utilization in asexual blood-stage malaria parasites? Competition experiments between NF54, clone 3D7, a lab-adapted African parasite, and a recently isolated Asian parasite (NHP4026) reveal contrasting outcomes in different media: 3D7 outcompetes NHP4026 in media containing human serum, while NHP4026 outcompetes 3D7 in media containing AlbuMAX, a commercial lipid-rich bovine serum formulation. To determine the basis for this polymorphism, we conducted parasite genetic crosses using humanized mice and compared genome-wide allele frequency changes in three independent progeny populations cultured in media containing human serum or AlbuMAX. This bulk segregant analysis detected three quantitative trait loci (QTL) regions [on chromosome (chr) 2 containing aspartate transaminase AST; chr 13 containing EBA-140; and chr 14 containing cysteine protease ATG4] linked with differential growth in serum or AlbuMAX in each of the three independent progeny pools. Selection driving differential growth was strong (s = 0.10 - 0.23 per 48-hour lifecycle). We conducted validation experiments for the strongest QTL on chr 13: competition experiments between ΔEBA-140 and 3D7 wildtype parasites showed fitness reversals in the two medium types as seen in the parental parasites, validating this locus as the causative gene. These results (i) demonstrate the effectiveness of bulk segregant analysis for dissecting fitness traits in P. falciparum genetic crosses, and (ii) reveal intimate links between red blood cell invasion and nutrient composition of growth media. Use of parasite crosses combined with bulk segregant analysis will allow systematic dissection of key nutrient acquisition/metabolism and red blood cell invasion pathways in P. falciparum.

Nutrient Limitation Magnifies Fitness Costs of Antimalarial Drug Resistance Mutations.

Drug resistance mutations tend to disrupt key physiological processes and frequently carry fitness costs, which are a central determinant of the rate of spread of these mutations in natural populations. Head-to-head competition assays provide a standard approach to measuring fitness for malaria parasites. These assays typically use a standardized culture medium containing RPMI 1640, which has a 1.4- to 5.5-fold higher concentration of amino acids than human blood. In this rich medium, we predict that fitness costs will be underestimated because resource competition is weak. We tested this prediction using an artemisinin-sensitive parasite edited to contain kelch-C580Y or R561H mutations conferring resistance to artemisinin or synonymous control mutations. We examined the impact of these single amino acid mutations on fitness, using replicated head-to-head competition experiments conducted in media containing (i) normal RPMI, (ii) modified RPMI with reduced amino acid concentration, (iii) RPMI containing only isoleucine, or (iv) 3-fold diluted RPMI. We found a significant 1.3- to 1.4-fold increase in fitness costs measured in modified and isoleucine-only media relative to normal media, while fitness costs were 2.5-fold higher in diluted media. We conclude that fitness costs are strongly affected by media composition and will be significantly underestimated in normal RPMI. Several components differed between media, including pABA and sodium bicarbonate concentrations, so we cannot directly determine which is responsible. Elevated fitness costs in nature will limit spread of artemisinin (ART) resistance but will also promote evolution of compensatory mutations that restore fitness and can be exploited to maximize selection in laboratory experiments.

Genomic analysis of a parasite invasion: Colonization of the Americas by the blood fluke Schistosoma mansoni.

Schistosoma mansoni, a snail-borne, blood fluke that infects humans, was introduced into the Americas from Africa during the Trans-Atlantic slave trade. As this parasite shows strong specificity to the snail intermediate host, we expected that adaptation to South American Biomphalaria spp. snails would result in population bottlenecks and strong signatures of selection. We scored 475,081 single nucleotide variants in 143 S. mansoni from the Americas (Brazil, Guadeloupe and Puerto Rico) and Africa (Cameroon, Niger, Senegal, Tanzania, and Uganda), and used these data to ask: (i) Was there a population bottleneck during colonization? (ii) Can we identify signatures of selection associated with colonization? (iii) What were the source populations for colonizing parasites? We found a 2.4- to 2.9-fold reduction in diversity and much slower decay in linkage disequilibrium (LD) in parasites from East to West Africa. However, we observed similar nuclear diversity and LD in West Africa and Brazil, suggesting no strong bottlenecks and limited barriers to colonization. We identified five genome regions showing selection in the Americas, compared with three in West Africa and none in East Africa, which we speculate may reflect adaptation during colonization. Finally, we infer that unsampled populations from central African regions between Benin and Angola, with contributions from Niger, are probably the major source(s) for Brazilian S. mansoni. The absence of a bottleneck suggests that this is a rare case of a serendipitous invasion, where S. mansoni parasites were pre-adapted to the Americas and able to establish with relative ease.

Cannabinoid control of gingival immune activation in chronically SIV-infected rhesus macaques involves modulation of the indoleamine-2,3-dioxygenase-1 pathway and salivary microbiome.

BACKGROUND: HIV/SIV-associated periodontal disease (gingivitis/periodontitis) (PD) represents a major comorbidity affecting people living with HIV (PLWH) on combination anti-retroviral therapy (cART). PD is characterized by chronic inflammation and dysbiosis. Nevertheless, the molecular mechanisms and use of feasible therapeutic strategies to reduce/reverse inflammation and dysbiosis remain understudied and unaddressed. METHODS: Employing a systems biology approach, we report molecular, metabolome and microbiome changes underlying PD and its modulation by phytocannabinoids [delta-9-tetrahydrocannabinol (Δ9-THC)] in uninfected and SIV-infected rhesus macaques (RMs) untreated (VEH-untreated/SIV) or treated with vehicle (VEH/SIV) or Δ9-THC (THC/SIV). FINDINGS: VEH- untreated/SIV but not THC/SIV RMs showed significant enrichment of genes linked to anti-viral defense, interferon-ß, NFκB, RIG-1, and JAK-STAT signaling. We focused on the anti-microbial DUOX1 and immune activation marker IDO1 that were reciprocally regulated in the gingiva of VEH-untreated/SIV RMs. Both proteins localized to the gingival epithelium and CD163+ macrophages, and showed differential expression in the gingiva of THC/SIV and VEH/SIV RMs. Additionally, inflammation-associated miR-21, miR-142-3p, miR-223, and miR-125a-5p showed significantly higher expression in the gingiva of VEH/SIV RMs. In human primary gingival epithelial cells, miR-125a-5p post-transcriptionally downregulated DUOX1 and THC inhibited IDO1 protein expression through a cannabinoid receptor-2 mediated mechanism. Interestingly, THC/SIV RMs showed relatively reduced plasma levels of kynurenine, kynurenate, and the neurotoxic quinolinate compared to VEH/SIV RMs at 5 months post SIV infection (MPI). Most importantly, THC blocked HIV/SIV-induced depletion of Firmicutes and Bacteroidetes, and reduced Gammaproteobacteria abundance in saliva. Reduced IDO1 protein expression was associated with significantly (p<0.05) higher abundance of Prevotella, Lactobacillus (L. salivarius, L. buchneri, L. fermentum, L. paracasei, L. rhamnosus, L. johnsonii) and Bifidobacteria and reduced abundance of the pathogenic Porphyromonas cangingivalis and Porphyromonas macacae at 5MPI. INTERPRETATION: The data provides deeper insights into the molecular mechanisms underlying HIV/SIV-induced PD and more importantly, the anti-inflammatory and anti-dysbiotic properties of THC in the oral cavity. Overall, these translational findings suggest that phytocannabinoids may help reduce gingival/systemic inflammation, salivary dysbiosis and potentially metabolic disease/syndrome in PLWH on cART and those with no access to cART or do not suppress the virus under cART. FUNDING: Research reported in this publication was supported by the National Institutes of Health Award Numbers R01DA052845 (MM and SNB), R01DA050169 (MM and CO), R01DA042524 and R56DE026930 (MM), and P51OD011104 and P51OD011133. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Genetic analysis of praziquantel response in schistosome parasites implicates a transient receptor potential channel.

Mass drug administration with praziquantel (PZQ) monotherapy is considered the mainstay for control and elimination of the parasites causing schistosomiasis in humans. This drug shows imperfect cure rates in the field, and parasites showing reduced PZQ response can be selected in the laboratory, but the extent of resistance in Schistosoma mansoni populations is unknown. We examined the genetic basis of the variation in response in a PZQ-selected S. mansoni population (SmLE-PZQ-R) in which 35% of the parasitic worms survive high-dose PZQ (73 micrograms per milliliter) treatment. We used genome-wide association to map loci underlying PZQ response and identified a transient receptor potential (Sm.TRPMPZQ) channel (Smp_246790) within the major chromosome 3 peak that is activated by nanomolar concentrations of PZQ. The PZQ response showed recessive inheritance and marker-assisted selection of parasites at a single Sm.TRPMPZQ SNP that produced populations of PZQ-enriched resistant (PZQ-ER) and PZQ-enriched sensitive (PZQ-ES) parasites, exhibiting >377-fold difference in PZQ response. The PZQ-ER parasites survived treatment in rodents at higher frequencies compared with PZQ-ES, and resistant parasites exhibited 2.25-fold lower expression of Sm.TRPMPZQ relative to sensitive parasites. Specific chemical blockers of Sm.TRPMPZQ enhanced PZQ resistance, whereas Sm.TRPMPZQ activators increased sensitivity. We surveyed Sm.TRPMPZQ sequence variations in 259 parasites from different global sites and identified one nonsense mutation that resulted in a truncated protein with no PZQ binding site. Our results demonstrate that Sm.TRPMPZQ underlies variation in PZQ responses in S. mansoni and provides an approach for monitoring emerging PZQ-resistant alleles in schistosome elimination programs.

Genetic architecture of transmission stage production and virulence in schistosome parasites.

Both theory and experimental data from pathogens suggest that the production of transmission stages should be strongly associated with virulence, but the genetic bases of parasite transmission/virulence traits are poorly understood. The blood fluke Schistosoma mansoni shows extensive variation in numbers of cercariae larvae shed and in their virulence to infected snail hosts, consistent with expected trade-offs between parasite transmission and virulence. We crossed schistosomes from two populations that differ 8-fold in cercarial shedding and in their virulence to Biomphalaria glabrata snail hosts, and determined four-week cercarial shedding profiles in F0 parents, F1 parents and 376 F2 progeny from two independent crosses in inbred snails. Sequencing and linkage analysis revealed that cercarial production is polygenic and controlled by five QTLs (i.e. Quantitative Trait Loci). These QTLs act additively, explaining 28.56% of the phenotypic variation. These results demonstrate that the genetic architecture of key traits relevant to schistosome ecology can be dissected using classical linkage mapping approaches.

The power and promise of genetic mapping from Plasmodium falciparum crosses utilizing human liver-chimeric mice.

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.

The hemolymph of Biomphalaria snail vectors of schistosomiasis supports a diverse microbiome.

The microbiome - the microorganism community that is found on or within an organism's body - is increasingly recognized to shape many aspects of its host biology and is a key determinant of health and disease. Microbiomes modulate the capacity of insect disease vectors (mosquitoes, tsetse flies, sandflies) to transmit parasites and disease. We investigate the diversity and abundance of microorganisms within the hemolymph (i.e. blood) of Biomphalaria snails, the intermediate host for Schistosoma mansoni, using Illumina MiSeq sequencing of the bacterial 16S V4 rDNA. We sampled hemolymph from five snails from six different laboratory populations of B. glabrata and one population of B. alexandrina. We observed 279.84 ± 0.79 amplicon sequence variants per snail. There were significant differences in microbiome composition at the level of individual snails, snail populations and species. Snail microbiomes were dominated by Proteobacteria and Bacteroidetes while water microbiomes from snail tank were dominated by Actinobacteria. We investigated the absolute bacterial load using qPCR: hemolymph samples contained 2784 ± 339 bacteria/µl. We speculate that the microbiome may represent a critical, but unexplored intermediary in the snail-schistosome interaction as hemolymph is in very close contact with the parasite at each step of its development.

The extended recovery ring-stage survival assay provides a superior association with patient clearance half-life and increases throughput.

BACKGROUND: Tracking and understanding artemisinin resistance is key for preventing global setbacks in malaria eradication efforts. The ring-stage survival assay (RSA) is the current gold standard for in vitro artemisinin resistance phenotyping. However, the RSA has several drawbacks: it is relatively low throughput, has high variance due to microscopy readout, and correlates poorly with the current benchmark for in vivo resistance, patient clearance half-life post-artemisinin treatment. Here a modified RSA is presented, the extended Recovery Ring-stage Survival Assay (eRRSA), using 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives, including parasite isolates with and without kelch13 mutations. METHODS: Plasmodium falciparum cultures were synchronized with single layer Percoll during the schizont stage of the intraerythrocytic development cycle. Cultures were left to reinvade to early ring-stage and parasitaemia was quantified using flow cytometry. Cultures were diluted to 2% haematocrit and 0.5% parasitaemia in a 96-well plate to start the assay, allowing for increased throughput and decreased variability between biological replicates. Parasites were treated with 700 nM of dihydroartemisinin or 0.02% dimethyl sulfoxide (DMSO) for 6 h, washed three times in drug-free media, and incubated for 66 or 114 h, when samples were collected and frozen for PCR amplification. A SYBR Green-based quantitative PCR method was used to quantify the fold-change between treated and untreated samples. RESULTS: 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives were assayed using the eRRSA. Due to the large number of pyknotic and dying parasites at 66 h post-exposure (72 h sample), parasites were grown for an additional cell cycle (114 h post-exposure, 120 h sample), which drastically improved correlation with patient clearance half-life compared to the 66 h post-exposure sample. A Spearman correlation of - 0.8393 between fold change and patient clearance half-life was identified in these 15 isolates from Southeast Asia, which is the strongest correlation reported to date. CONCLUSIONS: eRRSA drastically increases the efficiency and accuracy of in vitro artemisinin resistance phenotyping compared to the traditional RSA, which paves the way for extensive in vitro phenotyping of hundreds of artemisinin resistant parasites.

Oxamniquine resistance alleles are widespread in Old World Schistosoma mansoni and predate drug deployment.

Do mutations required for adaptation occur de novo, or are they segregating within populations as standing genetic variation? This question is key to understanding adaptive change in nature, and has important practical consequences for the evolution of drug resistance. We provide evidence that alleles conferring resistance to oxamniquine (OXA), an antischistosomal drug, are widespread in natural parasite populations under minimal drug pressure and predate OXA deployment. OXA has been used since the 1970s to treat Schistosoma mansoni infections in the New World where S. mansoni established during the slave trade. Recessive loss-of-function mutations within a parasite sulfotransferase (SmSULT-OR) underlie resistance, and several verified resistance mutations, including a deletion (p.E142del), have been identified in the New World. Here we investigate sequence variation in SmSULT-OR in S. mansoni from the Old World, where OXA has seen minimal usage. We sequenced exomes of 204 S. mansoni parasites from West Africa, East Africa and the Middle East, and scored variants in SmSULT-OR and flanking regions. We identified 39 non-synonymous SNPs, 4 deletions, 1 duplication and 1 premature stop codon in the SmSULT-OR coding sequence, including one confirmed resistance deletion (p.E142del). We expressed recombinant proteins and used an in vitro OXA activation assay to functionally validate the OXA-resistance phenotype for four predicted OXA-resistance mutations. Three aspects of the data are of particular interest: (i) segregating OXA-resistance alleles are widespread in Old World populations (4.29-14.91% frequency), despite minimal OXA usage, (ii) two OXA-resistance mutations (p.W120R, p.N171IfsX28) are particularly common (>5%) in East African and Middle-Eastern populations, (iii) the p.E142del allele has identical flanking SNPs in both West Africa and Puerto Rico, suggesting that parasites bearing this allele colonized the New World during the slave trade and therefore predate OXA deployment. We conclude that standing variation for OXA resistance is widespread in S. mansoni.

Striking differences in virulence, transmission and sporocyst growth dynamics between two schistosome populations.

BACKGROUND: Parasite traits associated with transmission success, such as the number of infective stages released from the host, are expected to be optimized by natural selection. However, in the trematode parasite Schistosoma mansoni, a key transmission trait, i.e. the number of cercariae larvae shed from infected Biomphalaria spp. snails, varies significantly within and between different parasite populations and selection experiments demonstrate that this variation has a strong genetic basis. In this study, we compared the transmission strategies of two laboratory schistosome population and their consequences for their snail host. METHODS: We infected inbred Biomphalaria glabrata snails using two S. mansoni parasite populations (SmBRE and SmLE), both isolated from Brazil and maintained in the laboratory for decades. We compared life history traits of these two parasite populations by quantifying sporocyst growth within infected snails (assayed using qPCR), output of cercaria larvae and impact on snail host physiological response (i.e. hemoglobin rate, laccase-like activity) and survival. RESULTS: We identified striking differences in virulence and transmission between the two studied parasite populations. SmBRE (low shedder (LS) parasite population) sheds very low numbers of cercariae and causes minimal impact on the snail physiological response (i.e. laccase-like activity, hemoglobin rate and snail survival). In contrast, SmLE (high shedder (HS) parasite population) sheds 8-fold more cercariae (mean ± SE cercariae per shedding: 284 ± 19 vs 2352 ± 113), causes high snail mortality and has strong impact on snail physiology. We found that HS sporocysts grow more rapidly inside the snail host, comprising up to 60% of cells within infected snails, compared to LS sporocysts, which comprised up to 31%. Cercarial production is strongly correlated to the number of S. mansoni sporocyst cells present within the snail host tissue, although the proportion of sporocyst cells alone does not explain the low cercarial shedding of SmBRE. CONCLUSIONS: We demonstrated the existence of alternative transmission strategies in the S. mansoni parasite consistent with trade-offs between parasite transmission and host survival: a "boom-bust" strategy characterized by high virulence, high transmission and short duration infections and a "slow and steady" strategy with low virulence, low transmission but long duration of snail host infections.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle.

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

Pairwise growth competitions identify relative fitness relationships among artemisinin resistant Plasmodium falciparum field isolates.

BACKGROUND: Competitive outcomes between co-infecting malaria parasite lines can reveal fitness disparities in blood stage growth. Blood stage fitness costs often accompany the evolution of drug resistance, with the expectation that relatively fitter parasites will be more likely to spread in populations. With the recent emergence of artemisinin resistance, it is important to understand the relative competitive fitness of the metabolically active asexual blood stage parasites. Genetically distinct drug resistant parasite clones with independently evolved sets of mutations are likely to vary in asexual proliferation rate, contributing to their chance of transmission to the mosquito vector. METHODS: An optimized in vitro 96-well plate-based protocol was used to quantitatively measure-head-to-head competitive fitness during blood stage development between seven genetically distinct field isolates from a hotspot of emerging artemisinin resistance and the laboratory strain, NF54. These field isolates were isolated from patients in Southeast Asia carrying different alleles of kelch13 and included both artemisinin-sensitive and artemisinin-resistant isolates. Fluorescent labeled microsatellite markers were used to track the relative densities of each parasite throughout the co-growth period of 14-60 days. All-on-all competitions were conducted for the panel of eight parasite lines (28 pairwise competitions) to determine their quantitative competitive fitness relationships. RESULTS: Twenty-eight pairwise competitive growth outcomes allowed for an unambiguous ranking among a set of seven genetically distinct parasite lines isolated from patients in Southeast Asia displaying a range of both kelch13 alleles and clinical clearance times and a laboratory strain, NF54. This comprehensive series of assays established the growth relationships among the eight parasite lines. Interestingly, a clinically artemisinin resistant parasite line that carries the wild-type form of kelch13 outcompeted all other parasites in this study. Furthermore, a kelch13 mutant line (E252Q) was competitively more fit without drug than lines with other resistance-associated kelch13 alleles, including the C580Y allele that has expanded to high frequencies under drug pressure in Southeast Asian resistant populations. CONCLUSIONS: This optimized competitive growth assay can be employed for assessment of relative growth as an index of fitness during the asexual blood stage growth between natural lines carrying different genetic variants associated with artemisinin resistance. Improved understanding of the fitness costs of different parasites proliferating in human blood and the role different resistance mutations play in the context of specific genetic backgrounds will contribute to an understanding of the potential for specific mutations to spread in populations, with the potential to inform targeted strategies for malaria therapy.