Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes.

Methods to deconvolve single-cell RNA-sequencing (scRNA-seq) data are necessary for samples containing a mixture of genotypes, whether they are natural or experimentally combined. Multiplexing across donors is a popular experimental design that can avoid batch effects, reduce costs and improve doublet detection. By using variants detected in scRNA-seq reads, it is possible to assign cells to their donor of origin and identify cross-genotype doublets that may have highly similar transcriptional profiles, precluding detection by transcriptional profile. More subtle cross-genotype variant contamination can be used to estimate the amount of ambient RNA. Ambient RNA is caused by cell lysis before droplet partitioning and is an important confounder of scRNA-seq analysis. Here we develop souporcell, a method to cluster cells using the genetic variants detected within the scRNA-seq reads. We show that it achieves high accuracy on genotype clustering, doublet detection and ambient RNA estimation, as demonstrated across a range of challenging scenarios.

Influence of the sickle cell trait on Plasmodium falciparum infectivity from naturally infected gametocyte carriers.

BACKGROUND: Sickle cell trait (SCT) refers to the carriage of one abnormal copy of the ß-globin gene, the HbS allele. SCT offers protection against malaria, controlling parasite density and preventing progression to symptomatic malaria. However, it remains unclear whether SCT also affects transmission stages and mosquito infection parameters. Deciphering the impact of the SCT on human to mosquito malaria transmission is key to understanding mechanisms that maintain the trait in malaria endemic areas. METHODS: The study was conducted from June to July 2017 among asymptomatic children living in the locality of Mfou, Cameroon. Blood samples were collected from asymptomatic children to perform malaria diagnosis by microscopy, Plasmodium species by PCR and hemoglobin typing by RFLP. Infectiousness of gametocytes to mosquitoes was assessed by membrane feeding assays using blood from gametocyte carriers of HbAA and HbAS genotypes. A zero-inflated model was fitted to predict distribution of oocysts in mosquitoes according to hemoglobin genotype of the gametocyte source. RESULTS: Among the 1557 children enrolled in the study, 314 (20.16%) were of the HbAS genotype. The prevalence of children with P. falciparum gametocytes was 18.47% in HbAS individuals and 13.57% in HbAA, and the difference is significant (χ2 = 4.61, P = 0.032). Multiplicity of infection was lower in HbAS gametocyte carriers (median = 2 genotypes/carrier in HbAS versus 3.5 genotypes/carrier in HbAA, Wilcoxon sum rank test = 188, P = 0.032). Gametocyte densities in the blood donor significantly influenced mosquito infection prevalence in both HbAS and HbAA individuals. The HbAS genotype had no significant effect on mosquito infection outcomes when using immune or naïve serum in feeding assays. In AB replacement feeding experiments, the odds ratio of mosquito infection for HbAA blood as compared to HbAS was 0.56 (95% CI 0.29-1.10), indicating a twice higher risk of infection in mosquitoes fed on gametocyte-containing blood of HbAS genotype. CONCLUSION: Plasmodium transmission stages were more prevalent in SCT individuals. This may reflect the parasite's enhanced investment in the sexual stage to increase their survival rate when asexual replication is impeded. The public health impact of our results points the need for intensive malaria control interventions in areas with high prevalence of HbAS. The similar infection parameters in feeding experiments where mosquitoes received the original serum from the blood donor indicated that immune responses to gametocyte surface proteins occur in both HbAS and HbAA individuals. The higher risk of infection in mosquitoes fed on HbAS blood depleted of immune factors suggests that changes in the membrane properties in HbAS erythrocytes may impact on the maturation process of gametocytes within circulating red blood cells.

Genetic Diversity of Plasmodium falciparum and Distribution of Antimalarial Drug Resistance Mutations in Symptomatic and Asymptomatic Infections.

Malaria control relies on passive case detection, and this strategy fails detecting asymptomatic infections. In addition, infections in endemic areas harbor multiple parasite genotypes that could affect case management and malaria epidemiology. Here, we performed AmpSeq genotyping to capture polymorphisms associated with antimalarial resistance and the genetic diversity within natural Plasmodium falciparum infections. Known genetic polymorphisms associated with altered drug susceptibility were screened for the five most common marker genes, pfdhfr, pfdhps, pfmdr1, pfcrt, and pfK13, and genetic diversity was established from two known AmpSeq markers, cpmp and csp. Relative abundance of the different genotypes within mixed infections was calculated from the number of reads per genotype. Genotyping was performed on 117 samples, 63 from asymptomatic and 54 from symptomatic individuals. We identified up to 15 genotypes within an infection, and the median multiplicity of infection was higher in asymptomatic infections (median MOI = 5 in asymptomatics versus median MOI = 2 in symptomatics, P < 0.001). No genetic differentiation on parasites from asymptomatic and symptomatic individuals was found. No mutation associated with ART resistance was identified. Prevalence of the P. falciparum chloroquine resistance wild-type genotype (CVMNK) reached 80%, confirming a return to chloroquine (CQ) sensitive parasites in Cameroon. In addition, the CQ-associated resistant genotype (CVIET) was present at very low density in polyclonal infections. Persistence of low-density chloroquine resistant parasites indicates competition-survival trade-offs may contribute to maintaining genetic diversity in natura. Thus, monitoring the expansion of these low-density genotypes in different immune backgrounds will be critical to evaluate drug policy changes.

Actin network disassembly powers dissemination of Listeria monocytogenes.

Several bacterial pathogens hijack the actin assembly machinery and display intracellular motility in the cytosol of infected cells. At the cell cortex, intracellular motility leads to bacterial dissemination through formation of plasma membrane protrusions that resolve into vacuoles in adjacent cells. Here, we uncover a crucial role for actin network disassembly in dissemination of Listeria monocytogenes. We found that defects in the disassembly machinery decreased the rate of actin tail turnover but did not affect the velocity of the bacteria in the cytosol. By contrast, defects in the disassembly machinery had a dramatic impact on bacterial dissemination. Our results suggest a model of L. monocytogenes dissemination in which the disassembly machinery, through local recycling of the actin network in protrusions, fuels continuous actin assembly at the bacterial pole and concurrently exhausts cytoskeleton components from the network distal to the bacterium, which enables membrane apposition and resolution of protrusions into vacuoles.

An essential role of the basal body protein SAS-6 in Plasmodium male gamete development and malaria transmission.

Gametocytes are the sole Plasmodium parasite stages that infect mosquitoes; therefore development of functional gametes is required for malaria transmission. Flagellum assembly of the Plasmodium male gamete differs from that of most other eukaryotes in that it is intracytoplasmic but retains a key conserved feature: axonemes assemble from basal bodies. The centriole/basal body protein SAS-6 normally regulates assembly and duplication of these organelles and its depletion causes severe flagellar/ciliary abnormalities in a diverse array of eukaryotes. Since basal body and flagellum assembly are intimately coupled to male gamete development in Plasmodium, we hypothesized that SAS-6 disruption may cause gametogenesis defects and perturb transmission. We show that Plasmodium berghei sas6 knockouts display severely abnormal male gametogenesis presenting reduced basal body numbers, axonemal assembly defects and abnormal nuclear allocation. The defects in gametogenesis reduce fertilization and render Pbsas6 knockouts less infectious to mosquitoes. Additionally, we show that lack of Pbsas6 blocks transmission from mosquito to vertebrate host, revealing an additional yet undefined role in ookinete to sporulating oocysts transition. These findings underscore the vulnerability of the basal body/SAS-6 to malaria transmission blocking interventions.

Proteomic analysis of the Plasmodium male gamete reveals the key role for glycolysis in flagellar motility.

BACKGROUND: Gametogenesis and fertilization play crucial roles in malaria transmission. While male gametes are thought to be amongst the simplest eukaryotic cells and are proven targets of transmission blocking immunity, little is known about their molecular organization. For example, the pathway of energy metabolism that power motility, a feature that facilitates gamete encounter and fertilization, is unknown. METHODS: Plasmodium berghei microgametes were purified and analysed by whole-cell proteomic analysis for the first time. Data are available via ProteomeXchange with identifier PXD001163. RESULTS: 615 proteins were recovered, they included all male gamete proteins described thus far. Amongst them were the 11 enzymes of the glycolytic pathway. The hexose transporter was localized to the gamete plasma membrane and it was shown that microgamete motility can be suppressed effectively by inhibitors of this transporter and of the glycolytic pathway. CONCLUSIONS: This study describes the first whole-cell proteomic analysis of the malaria male gamete. It identifies glycolysis as the likely exclusive source of energy for flagellar beat, and provides new insights in original features of Plasmodium flagellar organization.

A single cell atlas of sexual development in Plasmodium falciparum.

The developmental decision made by malaria parasites to become sexual underlies all malaria transmission. Here, we describe a rich atlas of short- and long-read single-cell transcriptomes of over 37,000 Plasmodium falciparum cells across intraerythrocytic asexual and sexual development. We used the atlas to explore transcriptional modules and exon usage along sexual development and expanded it to include malaria parasites collected from four Malian individuals naturally infected with multiple P. falciparum strains. We investigated genotypic and transcriptional heterogeneity within and among these wild strains at the single-cell level, finding differential expression between different strains even within the same host. These data are a key addition to the Malaria Cell Atlas interactive data resource, enabling a deeper understanding of the biology and diversity of transmission stages.

Bruton's tyrosine kinase regulates Shigella flexneri dissemination in HT-29 intestinal cells.

Shigella flexneri is a Gram-negative intracellular pathogen that infects the intestinal epithelium and utilizes actin-based motility to spread from cell to cell. S. flexneri actin-based motility has been characterized in various cell lines, but studies in intestinal cells are limited. Here we characterized S. flexneri actin-based motility in HT-29 intestinal cells. In agreement with studies conducted in various cell lines, we showed that S. flexneri relies on neural Wiskott-Aldrich Syndrome protein (N-WASP) in HT-29 cells. We tested the potential role of various tyrosine kinases involved in N-WASP activation and uncovered a previously unappreciated role for Bruton's tyrosine kinase (Btk) in actin tail formation in intestinal cells. We showed that Btk depletion led to a decrease in N-WASP phosphorylation which affected N-WASP recruitment to the bacterial surface, decreased the number of bacteria displaying actin-based motility, and ultimately affected the efficiency of spread from cell to cell. Finally, we showed that the levels of N-WASP phosphorylation and Btk expression were increased in response to infection, which suggests that S. flexneri infection not only triggers the production of proinflammatory factors as previously described but also manipulates cellular processes required for dissemination in intestinal cells.

The Prevalence of Human Plasmodium Species during Peak Transmission Seasons from 2016 to 2021 in the Rural Commune of Ntjiba, Mali.

Up-to-date knowledge of key epidemiological aspects of each Plasmodium species is necessary for making informed decisions on targeted interventions and control strategies to eliminate each of them. This study aims to describe the epidemiology of plasmodial species in Mali, where malaria is hyperendemic and seasonal. Data reports collected during high-transmission season over six consecutive years were analyzed to summarize malaria epidemiology. Malaria species and density were from blood smear microscopy. Data from 6870 symptomatic and 1740 asymptomatic participants were analyzed. The median age of participants was 12 years, and the sex ratio (male/female) was 0.81. Malaria prevalence from all Plasmodium species was 65.20% (95% CI: 60.10-69.89%) and 22.41% (CI: 16.60-28.79%) for passive and active screening, respectively. P. falciparum was the most prevalent species encountered in active and passive screening (59.33%, 19.31%). This prevalence was followed by P. malariae (1.50%, 1.15%) and P. ovale (0.32%, 0.06%). Regarding frequency, P. falciparum was more frequent in symptomatic individuals (96.77% vs. 93.24%, p = 0.014). In contrast, P. malariae was more frequent in asymptomatic individuals (5.64% vs. 2.45%, p < 0.001). P. ovale remained the least frequent species (less than 1%), and no P. vivax was detected. The most frequent coinfections were P. falciparum and P. malariae (0.56%). Children aged 5-9 presented the highest frequency of P. falciparum infections (41.91%). Non-falciparum species were primarily detected in adolescents (10-14 years) with frequencies above 50%. Only P. falciparum infections had parasitemias greater than 100,000 parasites per µL of blood. P. falciparum gametocytes were found with variable prevalence across age groups. Our data highlight that P. falciparum represented the first burden, but other non-falciparum species were also important. Increasing attention to P. malariae and P. ovale is essential if malaria elimination is to be achieved.

Systemic host inflammation induces stage-specific transcriptomic modification and slower maturation in malaria parasites.

Maturation rates of malaria parasites within red blood cells (RBCs) can be influenced by host nutrient status and circadian rhythm; whether host inflammatory responses can also influence maturation remains less clear. Here, we observed that systemic host inflammation induced in mice by an innate immune stimulus, lipopolysaccharide (LPS), or by ongoing acute Plasmodium infection, slowed the progression of a single cohort of parasites from one generation of RBC to the next. Importantly, plasma from LPS-conditioned or acutely infected mice directly inhibited parasite maturation during in vitro culture, which was not rescued by supplementation, suggesting the emergence of inhibitory factors in plasma. Metabolomic assessments confirmed substantial alterations to the plasma of LPS-conditioned and acutely infected mice, and identified a small number of candidate inhibitory metabolites. Finally, we confirmed rapid parasite responses to systemic host inflammation in vivo using parasite scRNA-seq, noting broad impairment in transcriptional activity and translational capacity specifically in trophozoites but not rings or schizonts. Thus, we provide evidence that systemic host inflammation rapidly triggered transcriptional alterations in circulating blood-stage Plasmodium trophozoites and predict candidate inhibitory metabolites in the plasma that may impair parasite maturation in vivo. IMPORTANCE Malaria parasites cyclically invade, multiply, and burst out of red blood cells. We found that a strong inflammatory response can cause changes to the composition of host plasma, which directly slows down parasite maturation. Thus, our work highlights a new mechanism that limits malaria parasite growth in the bloodstream.

Regulators of male and female sexual development are critical for the transmission of a malaria parasite.

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite.

PbGEST mediates malaria transmission to both mosquito and vertebrate host.

The malaria life cycle relies on the successful transfer of the parasite between its human and mosquito hosts. We identified a Plasmodium berghei secreted protein (PBANKA_131270) that plays distinct roles in both the mammal-to-mosquito and the mosquito-to-mammal transitions. This protein, here named gamete egress and sporozoite traversal (GEST), plays an important role in the egress of male and female gametes from the vertebrate red blood cell. Interestingly, GEST is also required following the bite of the infected mosquito, for sporozoite progression through the skin. We found PbGEST to be secreted shortly after activation of the intraerythrocytic gametocyte, and during sporozoite migration. These findings indicate that a single malaria protein may have pleiotropic roles in different parasites stages mediating transmission between its insect and mammalian hosts.

Use of a selective inhibitor to define the chemotherapeutic potential of the plasmodial hexose transporter in different stages of the parasite's life cycle.

During blood infection, malarial parasites use D-glucose as their main energy source. The Plasmodium falciparum hexose transporter (PfHT), which mediates the uptake of D-glucose into parasites, is essential for survival of asexual blood-stage parasites. Recently, genetic studies in the rodent malaria model, Plasmodium berghei, found that the orthologous hexose transporter (PbHT) is expressed throughout the parasite's development within the mosquito vector, in addition to being essential during intraerythrocytic development. Here, using a D-glucose-derived specific inhibitor of plasmodial hexose transporters, compound 3361, we have investigated the importance of D-glucose uptake during liver and transmission stages of P. berghei. Initially, we confirmed the expression of PbHT during liver stage development, using a green fluorescent protein (GFP) tagging strategy. Compound 3361 inhibited liver-stage parasite development, with a 50% inhibitory concentration (IC50) of 11 µM. This process was insensitive to the external D-glucose concentration. In addition, compound 3361 inhibited ookinete development and microgametogenesis, with IC50s in the region of 250 µM (the latter in a D-glucose-sensitive manner). Consistent with our findings for the effect of compound 3361 on vector parasite stages, 1 mM compound 3361 demonstrated transmission blocking activity. These data indicate that novel chemotherapeutic interventions that target PfHT may be active against liver and, to a lesser extent, transmission stages, in addition to blood stages.

Using scRNA-seq to Identify Transcriptional Variation in the Malaria Parasite Ookinete Stage.

The crossing of the mosquito midgut epithelium by the malaria parasite motile ookinete form represents the most extreme population bottleneck in the parasite life cycle and is a prime target for transmission blocking strategies. However, we have little understanding of the clonal variation that exists in a population of ookinetes in the vector, partially because the parasites are difficult to access and are found in low numbers. Within a vector, variation may result as a response to specific environmental cues or may exist independent of those cues as a potential bet-hedging strategy. Here we use single-cell RNA-seq to profile transcriptional variation in Plasmodium berghei ookinetes across different vector species, and between and within individual midguts. We then compare our results to low-input transcriptomes from individual Anopheles coluzzii midguts infected with the human malaria parasite Plasmodium falciparum. Although the vast majority of transcriptional changes in ookinetes are driven by development, we have identified candidate genes that may be responding to environmental cues or are clonally variant within a population. Our results illustrate the value of single-cell and low-input technologies in understanding clonal variation of parasite populations.

Uptake of Plasmodium falciparum Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding.

Plasmodium falciparum remains one of the leading causes of child mortality, and nearly half of the world's population is at risk of contracting malaria. While pathogenesis results from replication of asexual forms in human red blood cells, it is the sexually differentiated forms, gametocytes, which are responsible for the spread of the disease. For transmission to succeed, both mature male and female gametocytes must be taken up by a female Anopheles mosquito during its blood meal for subsequent differentiation into gametes and mating inside the mosquito gut. Observed circulating numbers of gametocytes in the human host are often surprisingly low. A pre-fertilization behavior, such as skin sequestration, has been hypothesized to explain the efficiency of human-to-mosquito transmission but has not been sufficiently tested due to a lack of appropriate tools. In this study, we describe the optimization of a qPCR tool that enables the relative quantification of gametocytes within very small input samples. Such a tool allows for the quantification of gametocytes in different compartments of the host and the vector that could potentially unravel mechanisms that enable highly efficient malaria transmission. We demonstrate the use of our gametocyte quantification method in mosquito blood meals from both direct skin feeding on Plasmodium gametocyte carriers and standard membrane feeding assay. Relative gametocyte abundance was not different between mosquitoes fed through a membrane or directly on the skin suggesting that there is no systematic enrichment of gametocytes picked up in the skin.

Artemisinin Bioactivity and Resistance in Malaria Parasites.

Artemisinin is the most widely-used compound against malaria and plays a critical role in the treatment of malaria worldwide. Resistance to artemisinin emerged about a decade ago in Southeast Asia and it is paramount to prevent its spread or emergence in Africa. Artemisinin has a complex mode of action and can cause widespread injury to many components of the parasite. In this review, we outline the different metabolic pathways affected by artemisinin, including the unfolded protein response, protein polyubiquitination, proteasome, phosphatidylinositol-3-kinase, and the eukaryotic translation initiation factor 2α. Based on recently published data, we present a model of how these different pathways interplay and how mutations in K13, the main identified resistance marker, may help parasites survive under artemisinin pressure.

The Malaria Cell Atlas: Single parasite transcriptomes across the complete Plasmodium life cycle.

Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites.

Single-cell RNA-seq reveals hidden transcriptional variation in malaria parasites.

Single-cell RNA-sequencing is revolutionising our understanding of seemingly homogeneous cell populations but has not yet been widely applied to single-celled organisms. Transcriptional variation in unicellular malaria parasites from the Plasmodium genus is associated with critical phenotypes including red blood cell invasion and immune evasion, yet transcriptional variation at an individual parasite level has not been examined in depth. Here, we describe the adaptation of a single-cell RNA-sequencing (scRNA-seq) protocol to deconvolute transcriptional variation for more than 500 individual parasites of both rodent and human malaria comprising asexual and sexual life-cycle stages. We uncover previously hidden discrete transcriptional signatures during the pathogenic part of the life cycle, suggesting that expression over development is not as continuous as commonly thought. In transmission stages, we find novel, sex-specific roles for differential expression of contingency gene families that are usually associated with immune evasion and pathogenesis.

Evaluation of the intra- and inter-specific genetic variability of Plasmodium lactate dehydrogenase.

BACKGROUND: Malaria diagnosis is vital to efficient control programmes and the recent advent of malaria rapid diagnostic tests (RDTs) provides a reliable and simple diagnostic method. However a characterization of the efficiency of these tests and the proteins they detect is needed to maximize RDT sensitivity. METHODS: Plasmodial lactate dehydrogenase (pLDH) gene of wild isolates of the four human species of Plasmodium from a variety of malaria endemic settings were sequenced and analysed. RESULTS: No variation in nucleotide was found within Plasmodium falciparum, synonymous mutations were found for Plasmodium malariae and Plasmodium. vivax; and three different types of amino acid sequence were found for Plasmodium ovale. Conserved and variable regions were identified within each species. CONCLUSION: The results indicate that antigen variability is unlikely to explain variability in performance of RDTs detecting pLDH from cases of P. falciparum, P. vivax or P. malariae malaria, but may contribute to poor detection of P. ovale.