Parasite Subpopulations Research Articles

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Genomic and epidemiological evidence for the emergence of a L. infantum/L. donovani hybrid with unusual epidemiology in northern Italy.

Leishmania (L.) infantum is one of the main causative agents of animal and human leishmaniasis across many endemic areas in South America, Europe, North Africa, and Asia. Despite its clinical significance, little is known about the genetic diversity of L. infantum circulating in a given endemic area. Here, we investigate this important open question by applying a comparative genomics approach to seven L. infantum isolates from different hosts and Italian regions, including the northern part of the country (Emilia-Romagna, RER), Sicily, and Sardinia, as an initial attempt to explore the breadth of parasite genetic heterogeneity in Italy. Additionally, microsatellite analysis was carried out to compare the isolates from RER with other 70 L. infantum strains from the same region as well as 65 strains belonging to the L. donovani complex from other countries. We revealed important karyotypic instability and identified strain-specific changes in gene dosage, which affected important virulence factors such as amastins and surface antigen-like proteins. Single nucleotide polymorphism-based clustering analysis of these genomes together with over 80 publicly available L. infantum and L. donovani genomes placed the Italian isolates into three geographically distinct clusters within the Mediterranean basin and uncovered three isolates clustering with putative L. infantum/L. donovani hybrids isolated in Cyprus. As judged by microsatellite profiling, these hybrid isolates are representative of a sub-population of parasites circulating in northern Italy that preferentially infect humans but not dogs. Our results place Italy at the crossroads of L. infantum infection in the Mediterranean and call attention to the public health risk represented by the introduction of non-European Leishmania species.IMPORTANCEThis study closes important knowledge gaps with respect to Leishmania (L.) infantum genetic heterogeneity in a given endemic country, as exemplified here for Italy, and reveals genetic hybridization as a main cause for re-emerging human leishmaniasis in northern Italy. The observed high diversity of Leishmania parasites on the Italian peninsula suggests different geographical origins, with genomic adaptation to various ecologies affecting both pathogenicity and transmission potential. This is documented by the discovery of a putative L. infantum/L. donovani hybrid strain, which has been shown to preferentially infect humans but not dogs. Our results provide important information to health authorities, which need to consider the public health risk represented by the introduction of new Leishmania species into EU countries due to population displacement or travel from countries where exotic/allochthonous parasite species are endemic.

Natural resistance to meglumine antimoniate is associated with treatment failure in cutaneous leishmaniasis caused by Leishmania (Viannia) panamensis.

The multifactorial basis of therapeutic response can obscure the relation between antimicrobial drug susceptibility and clinical outcome. To discern the relationship between parasite susceptibility to meglumine antimoniate (SbV) and therapeutic outcome of cutaneous leishmaniasis, risk factors for treatment failure were considered in evaluating this relationship in ninety-one cutaneous leishmaniasis patients and corresponding clinical strains of Leishmania (Viannia) panamensis. Parasite susceptibility to 32 μg SbV/mL (plasma Cmax) was evaluated in primary human macrophages, PBMCs, and U937 macrophages. Early parasitological response to treatment was determined in lesions of a subgroup of patients, and pathogenicity of Sb-resistant and sensitive clinical strains was compared in BALB/c mice. Parasite survival in cell models and patient lesions was determined by qRT-PCR of Leishmania 7SLRNA transcript. Parasite loads in BALB/c mice were quantified by limiting dilution analysis. The disparate Sb-susceptibility of parasite subpopulations distinguished by isoenzyme profiles (zymodemes) was manifest in all cell models. Notably, Sb-resistance defined by parasite survival, was most effectively discerned in U937 macrophages compared with primary human host cells, significantly higher among strains from patients who failed treatment than cured and, significantly associated with treatment failure. Each unit increase in transformed survival rate corresponded to a 10.6-fold rise in the odds of treatment failure. Furthermore, treatment failure was significantly associated with naturally Sb-resistant zymodeme 2.3 strains, which also produced larger lesions and parasite burdens in BALB/c mice than Sb-sensitive zymodeme 2.2 strains. The confounding effect of host risk factors for treatment failure in discerning this association was evidenced in comparing strains from patients with and without the defined risk factors for treatment failure. These results establish the association of natural resistance to meglumine antimoniate with treatment failure, the importance of host risk factors in evaluating drug susceptibility and treatment outcome, and the clinical and epidemiological relevance of natural Sb-resistance in L. (V.) panamensis subpopulations.

PbARID-associated chromatin remodeling events are essential for gametocyte development in Plasmodium.

Gametocyte development of the Plasmodium parasite is a key step for transmission of the parasite. Male and female gametocytes are produced from a subpopulation of asexual blood-stage parasites, but the mechanisms that regulate the differentiation of sexual stages are still under investigation. In this study, we investigated the role of PbARID, a putative subunit of a SWI/SNF chromatin remodeling complex, in transcriptional regulation during the gametocyte development of P. berghei. PbARID expression starts in early gametocytes before the manifestation of male and female-specific features, and disruption of its gene results in the complete loss of gametocytes with detectable male features and the production of abnormal female gametocytes. ChIP-seq analysis of PbARID showed that it forms a complex with gSNF2, an ATPase subunit of the SWI/SNF chromatin remodeling complex, associating with the male cis-regulatory element, TGTCT. Further ChIP-seq of PbARID in gsnf2-knockout parasites revealed an association of PbARID with another cis-regulatory element, TGCACA. RIME and DNA-binding assays suggested that HDP1 is the transcription factor that recruits PbARID to the TGCACA motif. Our results indicated that PbARID could function in two chromatin remodeling events and paly essential roles in both male and female gametocyte development.

Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique.

Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigate transcriptional heterogeneity and stochastic gene expression in Plasmodium falciparum by performing the quasilinear multiple annealing and looping based amplification cycles (MALBAC) based amplification and single cell RNA sequencing of blood stage schizonts. Our data reveals significant transcriptional variations in the schizont stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, the data reflects several diversification processes including putative developmental “checkpoint”; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability are preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we apply quantitative RT-PCR and RNA-FISH approach and confirm stochastic expression of key invasion genes, such as, msp1, msp3, msp7, eba181 and ama1 which represent prime candidates for invasion-blocking vaccines.

Impact of Drug Pressure versus Limited Access to Drug in Malaria Control: The Dilemma.

Malaria burden has severe impact on the world. Several arsenals, including the use of antimalarials, are in place to curb the malaria burden. However, the application of these antimalarials has two extremes, limited access to drug and drug pressure, which may have similar impact on malaria control, leading to treatment failure through divergent mechanisms. Limited access to drugs ensures that patients do not get the right doses of the antimalarials in order to have an effective plasma concentration to kill the malaria parasites, which leads to treatment failure and overall reduction in malaria control via increased transmission rate. On the other hand, drug pressure can lead to the selection of drug resistance phenotypes in a subpopulation of the malaria parasites as they mutate in order to adapt. This also leads to a reduction in malaria control. Addressing these extremes in antimalarial application can be essential in maintaining the relevance of the conventional antimalarials in winning the war against malaria.

Chemical Derivatization and Characterization of Novel Antitrypanosomals for African Trypanosomiasis.

The search for novel antitrypanosomals and the investigation into their mode of action remain crucial due to the toxicity and resistance of commercially available antitrypanosomal drugs. In this study, two novel antitrypanosomals, tortodofuordioxamide (compound 2) and tortodofuorpyramide (compound 3), were chemically derived from the natural N-alkylamide tortozanthoxylamide (compound 1) through structural modification. The chemical structures of these compounds were confirmed through spectrometric and spectroscopic analysis, and their in vitro efficacy and possible mechanisms of action were, subsequently, investigated in Trypanosoma brucei (T. brucei), one of the causative species of African trypanosomiasis (AT). The novel compounds 2 and 3 displayed significant antitrypanosomal potencies in terms of half-maximal effective concentrations (EC50) and selectivity indices (SI) (compound 1, EC50 = 7.3 μM, SI = 29.5; compound 2, EC50 = 3.2 μM, SI = 91.3; compound 3, EC50 = 4.5 μM, SI = 69.9). Microscopic analysis indicated that at the EC50 values, the compounds resulted in the coiling and clumping of parasite subpopulations without significantly affecting the normal ratio of nuclei to kinetoplasts. In contrast to the animal antitrypanosomal drug diminazene, compounds 1, 2 and 3 exhibited antioxidant absorbance properties comparable to the standard antioxidant Trolox (Trolox, 0.11 A; diminazene, 0.50 A; compound 1, 0.10 A; compound 2, 0.09 A; compound 3, 0.11 A). The analysis of growth kinetics suggested that the compounds exhibited a relatively gradual but consistent growth inhibition of T. brucei at different concentrations. The results suggest that further pharmacological optimization of compounds 2 and 3 may facilitate their development into novel AT chemotherapy.

Deep amplicon sequencing highlights low intra-host genetic variability of Echinococcus multilocularis and high prevalence of the European-type haplotypes in coyotes and red foxes in Alberta, Canada.

Echinococcus multilocularis (Em) is a zoonotic parasite considered a global emergent pathogen. Recent findings indicate that the parasite is expanding its range in North America and that European-type haplotypes are circulating in western Canada. However, genetic analyses are usually conducted only on a few parasites out of thousands of individuals within each definitive host, likely underestimating the prevalence of less common haplotypes. Moreover, mixed infections with several mtDNA haplotypes in the same host have been reported, but their relative abundance within the host was never estimated. We aimed to 1) estimate the frequency of co-infections of different Em haplotypes in coyotes (Canis latrans) and red foxes (Vulpes vulpes) from western Canada and their relative abundance within the definitive hosts, 2) detect less prevalent haplotypes by sampling a larger proportion of the parasite subpopulation per host, and 3) investigate differences in the distribution of Em haplotypes in these main definitive hosts; foxes and coyotes. We extracted DNA from ~10% of the worm subpopulation per host (20 foxes and 47 coyotes) and used deep amplicon sequencing (NGS technology) on four loci, targeting the most polymorphic regions from the mitochondrial genes cox1 (814 bp), nad1 (344 bp), and cob (387 bp). We detected the presence of mixed infections with multiple Em haplotypes and with different Echinococcus species including Em and E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the European-type haplotypes in both hosts. Our results suggest a population expansion of the European over the North American strain in Alberta and a limited distribution of some European-type haplotypes. Our findings indicate that deep amplicon sequencing represents a valuable tool to characterize Em in multiple hosts, to assess the current distribution and possible origins of the European strain in North America. The potential use of next-generation sequencing technologies is particularly important to understand the patterns of geographic expansion of this parasite.

Characterization and Follow-Up of Trypanosoma cruzi Natural Populations Refractory to Etiological Chemotherapy in Oral Chagas Disease Patients.

We aimed to characterize the genetic constitution of natural T. cruzi populations involved in an Oral Chagas Disease (OCD) outbreak at a rural school of the community of Chichiriviche de la Costa, Venezuela, which affected patients did not respond to the etiological treatment. Peripheral blood samples and/or hemocultures were obtained from twenty-nine OCD patients at time of diagnosis or along nine years of Post-treatment (Tx) follow-up. The IgG serology, T. cruzi discrete typing units (DTU), satellite DNA-qPCR parasitic loads, and minicircle signatures were determined at Pre-Tx and after Tx. The serological titles and parasitic loads changed after treatment, with a significant decrease of IgG titers (Spearman’s r value= -0.961) and median parasite loads from 2.869 [IQR = 2.113 to 3.720] to 0.105 [IQR = -1.147 to 1.761] log10 par eq. /mL at Pre-Tx and Post-Tx, respectively, suggesting infection evolution from acute to chronic phase, without seroconversion or parasitological eradication, which was indicative of treatment failure. All patients were infected with T. cruzi DTU I populations. At Pre-Tx their median Jaccard genetic distances were 0.775 [IQR = 0.708 to 0.882], decreasing in genetic variability towards the end of follow-up (Mann-Whitney U test p= 0.0031). Interestingly, no Post-Tx minicircle signature was identical to its Pre-Tx counterpart population in a same patient, revealing selection of parasite subpopulations between the primary infection and Post-Tx. The parasitic populations isolated from hemocultures showed a lower number of bands in the minicircle signatures with respect to the signatures obtained directly from the patients’ blood samples, demonstrating a process of parasitic selection and reduction of the population variability that initially infected the patients. Decrease of parasitic loads after treatment as well as Pre- and Post-Tx intra-TcI diversity might be a consequence of both, natural evolution of the acute infection to the chronic phase and persistence of refractory populations due to Tx selection.

Limited diversity in the CD8+ antigen-coding loci in Theileria parva parasites from cattle from southern and eastern Africa

Theileria parva infections in cattle causes huge economic losses in the affected African countries, directly impacting the livelihood of the poor small-holder farmers. The current immunization protocol using live sporozoites in eastern Africa, is among the control measures designed to limit T. parva infections in cattle. However, the ability of the immune protection induced by this immunization to protect against field parasites has been compromised by the diversity of the parasite involving the schizont antigen genes. Previous studies have reported on the antigenic diversity of T. parva parasites from southern and eastern Africa, however, similar reports on T. parva parasites particularly from cattle from southern Africa remains scanty, due to the self-limiting nature of Corridor disease. Thus, we evaluated the diversity of CD8+ T-cell regions of ten schizont antigen genes in T. parva parasites associated with Corridor disease and East Coast fever (ECF) from southern and eastern Africa respectively.Regions of schizont antigen (TpAg) genes containing the CD8+ T-cell epitopes (CTL determinants) were amplified from genomic DNA extracted from blood of T. parva positive samples, cloned and sequenced. The results revealed limited diversity between the two parasite groups from cattle from southern and eastern Africa, defying the widely accepted notion that antigen-encoding loci in cattle-derived parasites are conserved, while in buffalo-derived parasites, they are extensively variable. This suggests that only a sub-population of parasites is successfully transmitted from buffalo to cattle, resulting in the limited antigenic diversity in Corridor disease parasites. Tp4, Tp5, Tp7 and Tp8 showed limited to absence of diversity in both parasite groups, suggesting the need to further investigate their immunogenic properties for consideration as candidates for a subunit vaccine. Distinct and common variants of Tp2 were detected among the ECF parasites from eastern Africa indicating evidence of parasite mixing following immunization. This study provides additional information on the comparative diversity of TpAg genes in buffalo- and cattle-derived T. parva parasites from cattle from southern and eastern Africa.

Genomic surveillance of Plasmodium falciparum and Plasmodium vivax cases at the University Hospital in Tegucigalpa, Honduras

Malaria continues to be an important health problem in Honduras despite major progress achieved reducing its incidence in the last two decades. In a context of case reduction, continuing surveillance of parasite diversity and drug resistance is an important component to assist effective malaria control strategies and support risk assessments. In this study, we employed next generation sequencing on collected Plasmodium vivax and P. falciparum samples from the Hospital Escuela (University Hospital) in Honduras between 2005 and 2017. Hospital Escuela is the main public health hospital in Honduras and receives suspected malaria cases from endemic regions within the country. The resulting sequencing data was used to assess complexity of infections, parasite population structure, parasite diversity and drug resistance profiling. All P. vivax samples and all autochtonous P. falciparum samples were monoclonal and presented a low intra population diversity (π = 0.25 and 0.07, respectively). Genotyping of drug resistance markers showed that three P. falciparum samples presented the chloroquine resistant haplotype SVMNT on pfcrtr (positions 72–76). Epidemiological data suggested that two of these samples were imported cases from Africa whereas the third one was a local case. Three suspected imported cases (two of which were also pfcrt mutants) presented the pfmdr1 86Y mutation that further enhances the CQ resistant genotype. No evidence was found for kelch13 artemisinin resistance associated mutations nor parasite genetic background mutations. Discriminant analysis of principal components and phylogenetic analysis showed two P. vivax and two P. falciparum parasite sub-populations with limited recombination between them. It also confirmed the closer relationship of the three imported cases with African strains. Our findings showed that local Honduras P. falciparum strains do not hold CQ resistance polymorphisms which aligns with clinical data reported by the country and supports the continuity of CQ based treatment in Honduras. In addition, our findings highlight the need of using genomic approaches to provide key information about parasite biology including drug resistance, population structure and HRP2/HRP3 deletions which are becoming relevant as the country move towards elimination.

Efficient Surveillance of Plasmodium knowlesi Genetic Subpopulations, Malaysian Borneo, 2000-2018.

Population genetic analysis revealed that Plasmodium knowlesi infections in Malaysian Borneo are caused by 2 divergent parasites associated with long-tailed (cluster 1) and pig-tailed (cluster 2) macaques. Because the transmission ecology is likely to differ for each macaque species, we developed a simple genotyping PCR to efficiently distinguish between and survey the 2 parasite subpopulations. This assay confirmed differences in the relative proportions in areas of Kapit division of Sarawak state, consistent with multilocus microsatellite analyses. Analyses of 1,204 human infections at Kapit Hospital showed that cluster 1 caused approximately two thirds of cases with no significant temporal changes from 2000 to 2018. We observed an apparent increase in overall numbers in the most recent 2 years studied, driven mainly by increased cluster 1 parasite infections. Continued monitoring of the frequency of different parasite subpopulations and correlation with environmental alterations are necessary to determine whether the epidemiology will change substantially.

Plasma-derived extracellular vesicles from Plasmodium vivax patients signal spleen fibroblasts via NF-kB facilitating parasite cytoadherence

Plasmodium vivax is the most widely distributed human malaria parasite. Previous studies have shown that circulating microparticles during P. vivax acute attacks are indirectly associated with severity. Extracellular vesicles (EVs) are therefore major components of circulating plasma holding insights into pathological processes. Here, we demonstrate that plasma-derived EVs from Plasmodium vivax patients (PvEVs) are preferentially uptaken by human spleen fibroblasts (hSFs) as compared to the uptake of EVs from healthy individuals. Moreover, this uptake induces specific upregulation of ICAM-1 associated with the translocation of NF-kB to the nucleus. After this uptake, P. vivax-infected reticulocytes obtained from patients show specific adhesion properties to hSFs, reversed by inhibiting NF-kB translocation to the nucleus. Together, these data provide physiological EV-based insights into the mechanisms of human malaria pathology and support the existence of P. vivax-adherent parasite subpopulations in the microvasculature of the human spleen.

Whole genome sequencing of amplified Plasmodium knowlesi DNA from unprocessed blood reveals genetic exchange events between Malaysian Peninsular and Borneo subpopulations

The zoonotic Plasmodium knowlesi parasite is the most common cause of human malaria in Malaysia. Genetic analysis has shown that the parasites are divided into three subpopulations according to their geographic origin (Peninsular or Borneo) and, in Borneo, their macaque host (Macaca fascicularis or M. nemestrina). Whilst evidence suggests that genetic exchange events have occurred between the two Borneo subpopulations, the picture is unclear in less studied Peninsular strains. One difficulty is that P. knowlesi infected individuals tend to present with low parasitaemia leading to samples with insufficient DNA for whole genome sequencing. Here, using a parasite selective whole genome amplification approach on unprocessed blood samples, we were able to analyse recent genomes sourced from both Peninsular Malaysia and Borneo. The analysis provides evidence that recombination events are present in the Peninsular Malaysia parasite subpopulation, which have acquired fragments of the M. nemestrina associated subpopulation genotype, including the DBPβ and NBPXa erythrocyte invasion genes. The NBPXb invasion gene has also been exchanged within the macaque host-associated subpopulations of Malaysian Borneo. Our work provides strong evidence that exchange events are far more ubiquitous than expected and should be taken into consideration when studying the highly complex P. knowlesi population structure.

High genetic diversity and differentiation of theBabesia ovispopulation in Turkey

Babesia ovis is a tick-transmitted protozoan haemoparasite causing ovine babesiosis in sheep and goats leading to considerable economic loss in Turkey and neighbouring countries. There are no vaccines available, therapeutic drugs leave toxic residues in meat and milk, and tick vector control entails environmental risks. A panel of eight mini- and micro-satellite marker loci was developed and applied to study genetic diversity and substructuring of B.ovis from western, central and eastern Turkey. A high genetic diversity (He =0.799) was found for the sample of overall B.ovis population (n=107) analyzed. Principle component analysis (PCoA) revealed the existence of three parasite subpopulations: (a) a small subpopulation of isolates from Aydin, western Turkey; (b) a second cluster predominantly generated by isolates from western Turkey; and (c) a third cluster predominantly formed by isolates from central and eastern Turkey. Two B.ovis isolates from Israel included in the analysis clustered with isolates from central and eastern Turkey. This finding strongly suggests substructuring of a major Turkish population into western versus central-eastern subpopulations, while the additional smaller B.ovis population found in Aydin could have been introduced, more recently, to Turkey. STRUCTURE analysis suggests a limited exchange of parasite strains between the western and the central-eastern regions and vice versa, possibly due to limited trading of sheep. Importantly, evidence for recombinant genotypes was obtained in regionally interchanged parasite isolates. Important climatic differences between the western and the central/eastern region, with average yearly temperatures of 21°C versus 15°C, correspond with the identified geographical substructuring. We hypothesize that the different climatic conditions may result in variation in the activity of subpopulations of Rhipicephalus spp. tick vectors, which, in turn, could selectively maintain and transmit different parasite populations. These findings may have important implications for vaccine development and the spread of drug resistance.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil.

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

Rare but evolutionarily consequential outcrossing in a highly inbred zoonotic parasite

Recurrent self-mating can result in nearly clonal propagation of biological lineages, but even occasional outcrossing can serve to redistribute variation in future generations, providing cohesion among regional populations. The zoonotic parasite Trichinella spiralis has been suspected to undergo frequent inbreeding, resulting in genetically uniform larval cohorts which differ markedly from one another. Here, we explored the extent of inbreeding for this parasite by determining how genetic variation (at variable microsatellite markers) is distributed among 1379 larvae derived from 41 wild boars in Extremadura, Spain. In particular, we sought to determine how much of the genetic variation in this region’s parasites occurs among the larvae of any given wild boar, and whether each derives from one, or more, parental lineages. We found strong evidence for inbreeding, resulting in genetically distinct parasite subpopulations among the parasites derived from many pairs of wild boar. Fully two-thirds of these parasite cohorts appear to derive from inbred parents; in 10% of the wild boars, parasites were so inbred as to become absolutely fixed in all of the assayed genetic loci. In spite of this, more than one pair of parents appear to have given rise to the infections in one-third of the sampled wild boars, resulting in mixed infections. These mixed infections should slow losses of heterozygosity and multi-locus polymorphism in any given parasite lineage. Such outcrossing should limit distinctions that would otherwise accumulate among transmission chains, thereby enforcing cohesion through the region’s population in spite of its marked departure from panmixia. Conditions of transmission may differ in other regions, where such epidemiological features may engender different evolutionary outcomes.

Genome-wide mosaicism in divergence between zoonotic malaria parasite subpopulations with separate sympatric transmission cycles.

Plasmodium knowlesi is a significant cause of human malaria transmitted as a zoonosis from macaque reservoir hosts in South‐East Asia. Microsatellite genotyping has indicated that human infections in Malaysian Borneo are an admixture of two highly divergent sympatric parasite subpopulations that are, respectively, associated with long‐tailed macaques (Cluster 1) and pig‐tailed macaques (Cluster 2). Whole‐genome sequences of clinical isolates subsequently confirmed the separate clusters, although fewer of the less common Cluster 2 type were sequenced. Here, to analyse population structure and genomic divergence in subpopulation samples of comparable depth, genome sequences were generated from 21 new clinical infections identified as Cluster 2 by microsatellite analysis, yielding a cumulative sample size for this subpopulation similar to that for Cluster 1. Profound heterogeneity in the level of intercluster divergence was distributed across the genome, with long contiguous chromosomal blocks having high or low divergence. Different mitochondrial genome clades were associated with the two major subpopulations, but limited exchange of haplotypes from one to the other was evident, as was also the case for the maternally inherited apicoplast genome. These findings indicate deep divergence of the two sympatric P. knowlesi subpopulations, with introgression likely to have occurred recently. There is no evidence yet of specific adaptation at any introgressed locus, but the recombinant mosaic types offer enhanced diversity on which selection may operate in a currently changing landscape and human environment. Loci responsible for maintaining genetic isolation of the sympatric subpopulations need to be identified in the chromosomal regions showing fixed differences.

Single-cell RNA sequencing reveals a signature of sexual commitment in malaria parasites.

Pathogens have to balance transmission with persistence. For Plasmodium falciparum, the most widespread and virulent malaria parasite, persistence within its human host requires continuous asexual replication within red blood cells, while its mosquito-borne transmission depends on intra-erythrocytic differentiation into non-replicating sexual stages called gametocytes. Commitment to either fate is determined during the preceding cell cycle that begins with invasion by a single, asexually committed merozoite and ends, 48 hours later, with a schizont releasing newly formed merozoites, all committed to either continued asexual replication or differentiation into gametocytes. Sexual commitment requires the transcriptional activation of ap2-g (PF3D7_1222600), the master regulator of sexual development, from an epigenetically silenced state during asexual replication. AP2-G expression during this 'commitment cycle' prepares gene expression in nascent merozoites to initiate sexual development through a hitherto unknown mechanism. To maintain a persistent infection, the expression of ap2-g is limited to a sub-population of parasites (1-30%, depending on genetic background and growth conditions). As sexually committed schizonts comprise only a sub-population and are morphologically indistinguishable from their asexually committed counterparts, defining their characteristic gene expression has been difficult using traditional, bulk transcriptome profiling. Here we use highly parallel, single-cell RNA sequencing of malaria cultures undergoing sexual commitment to determine the transcriptional changes induced by AP2-G within this sub-population. By analysing more than 18,000 single parasite transcriptomes from a conditional AP2-G knockdown line and NF54 wild-type parasites at multiple stages of development, we show that sexually committed, AP2-G+ mature schizonts specifically upregulate additional regulators of gene expression, including other AP2 transcription factors, histone-modifying enzymes, and regulators of nucleosome positioning. These epigenetic regulators may act to facilitate the expression and/or repression of genes that are necessary for the initiation of gametocyte development in the subsequent cell cycle.