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

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

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.

Ready for renascence in mosquito: The regulation of gene expression in Plasmodium sexual development

Malaria remains one of the most perilous vector-borne infectious diseases for humans globally. Sexual gametocyte represents the exclusive stage at which malaria parasites are transmitted from the vertebrate to the Anopheles host. The feasible and effective approach to prevent malaria transmission is by addressing the sexual developmental processes, that is, gametocytogenesis and gametogenesis. Thus, this review will comprehensively cover advances in the regulation of gene expression surrounding the transmissible stages, including epigenetic, transcriptional, and post-transcriptional control.

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.

Positive-unlabeled learning identifies vaccine candidate antigens in the malaria parasite Plasmodium falciparum

Malaria vaccine development is hampered by extensive antigenic variation and complex life stages of Plasmodium species. Vaccine development has focused on a small number of antigens, many of which were identified without utilizing systematic genome-level approaches. In this study, we implement a machine learning-based reverse vaccinology approach to predict potential new malaria vaccine candidate antigens. We assemble and analyze P. falciparum proteomic, structural, functional, immunological, genomic, and transcriptomic data, and use positive-unlabeled learning to predict potential antigens based on the properties of known antigens and remaining proteins. We prioritize candidate antigens based on model performance on reference antigens with different genetic diversity and quantify the protein properties that contribute most to identifying top candidates. Candidate antigens are characterized by gene essentiality, gene ontology, and gene expression in different life stages to inform future vaccine development. This approach provides a framework for identifying and prioritizing candidate vaccine antigens for a broad range of pathogens.

The malarial blood transcriptome: translational applications

The blood transcriptome of malaria patients has been used extensively to elucidate the pathophysiological mechanisms and host immune responses to disease, identify candidate diagnostic and prognostic biomarkers, and reveal new therapeutic targets for drug discovery. This review gives a high-level overview of the three main translational applications of these studies (diagnostics, prognostics, and therapeutics) by summarising recent literature and outlining the main limitations and future directions of each application. It highlights the need for consistent and accurate definitions of disease states and subject groups and discusses how prognostic studies must distinguish clearly between analyses that attempt to predict future disease states and those which attempt to discriminate between current disease states (classification). Lastly it examines how many promising therapeutics fail due to the choice of imperfect animal models for pre-clinical testing and lack of appropriate validation studies in humans, and how future transcriptional studies may be utilised to overcome some of these limitations.

Unraveling the intricacies of host-pathogen interaction through single-cell genomics

Single-cell genomics provide researchers with tools to assess host-pathogen interactions at a resolution previously inaccessible. Transcriptome analysis, epigenome analysis, and immune profiling techniques allow for a better comprehension of the heterogeneity underlying both the host response and infectious agents. Here, we highlight technological advancements and data analysis workflows that increase our understanding of host-pathogen interactions at the single-cell level. We review various studies that have used these tools to better understand host-pathogen dynamics in a variety of infectious disease contexts, including viral, bacterial, and parasitic diseases. We conclude by discussing how single-cell genomics can advance our understanding of host-pathogen interactions.

Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission

Most Plasmodium vivax infections contain genetically distinct parasites, but the consequences of this polyclonality on the development of asexual parasites, their sexual differentiation, and their transmission remain unknown. We describe infections of Saimiri monkeys with two strains of P. vivax and the analyses of 117,350 parasites characterized by single cell RNA sequencing and individually genotyped. In our model, consecutive inoculations fail to establish polyclonal infections. By contrast, simultaneous inoculations of two strains lead to sustained polyclonal infections, although without detectable differences in parasite regulation or sexual commitment. Analyses of sporozoites dissected from mosquitoes fed on coinfected monkeys show that all genotypes are successfully transmitted to mosquitoes. However, after sporozoite inoculation, not all genotypes contribute to the subsequent blood infections, highlighting an important bottleneck during pre-erythrocytic development. Overall, these studies provide new insights on the mechanisms regulating the establishment of polyclonal P. vivax infections and their consequences for disease transmission.

Single cell transcriptomics shows that malaria promotes unique regulatory responses across multiple immune cell subsets

Plasmodium falciparum malaria drives immunoregulatory responses across multiple cell subsets, which protects from immunopathogenesis, but also hampers the development of effective anti-parasitic immunity. Understanding malaria induced tolerogenic responses in specific cell subsets may inform development of strategies to boost protective immunity during drug treatment and vaccination. Here, we analyse the immune landscape with single cell RNA sequencing during P. falciparum malaria. We identify cell type specific responses in sub-clustered major immune cell types. Malaria is associated with an increase in immunosuppressive monocytes, alongside NK and γδ T cells which up-regulate tolerogenic markers. IL-10-producing Tr1 CD4 T cells and IL-10-producing regulatory B cells are also induced. Type I interferon responses are identified across all cell types, suggesting Type I interferon signalling may be linked to induction of immunoregulatory networks during malaria. These findings provide insights into cell-specific and shared immunoregulatory changes during malaria and provide a data resource for further analysis.

Extracellular vesicles could be a putative posttranscriptional regulatory mechanism that shapes intracellular RNA levels in Plasmodium falciparum

Plasmodium falciparum secretes extracellular vesicles (PfEVs) that contain parasite-derived RNA. However, the significance of the secreted RNA remains unexplored. Here, we compare secreted and intracellular RNA from asexual cultures of six P. falciparum lines. We find that secretion of RNA via extracellular vesicles is not only periodic throughout the asexual intraerythrocytic developmental cycle but is also highly conserved across P. falciparum isolates. We further demonstrate that the phases of RNA secreted via extracellular vesicles are discernibly shifted compared to those of the intracellular RNA within the secreting whole parasite. Finally, transcripts of genes with no known function during the asexual intraerythrocytic developmental cycle are enriched in PfEVs compared to the whole parasite. We conclude that the secretion of extracellular vesicles could be a putative posttranscriptional RNA regulation mechanism that is part of or synergise the classic RNA decay processes to maintain intracellular RNA levels in P. falciparum.

Dual-scRNA-seq analysis reveals rare and uncommon parasitized cell populations in chronic L. donovani infection

Although phagocytic cells are documented targets of Leishmania parasites, it is unclear whether other cell types can be infected. Here, we use unbiased single-cell RNA sequencing (scRNA-seq) to simultaneously analyze host cell and Leishmania donovani transcriptomes to identify and annotate parasitized cells in spleen and bone marrow in chronically infected mice. Our dual-scRNA-seq methodology allows the detection of heterogeneous parasitized populations. In the spleen, monocytes and macrophages are the dominant parasitized cells, while megakaryocytes, basophils, and natural killer (NK) cells are found to be unexpectedly infected. In the bone marrow, the hematopoietic stem cells (HSCs) expressing phagocytic receptors FcγR and CD93 are the main parasitized cells. Additionally, we also detect parasitized cycling basal cells, eosinophils, and macrophages in chronically infected mice. Flow cytometric analysis confirms the presence of parasitized HSCs. Our unbiased dual-scRNA-seq method identifies rare, parasitized cells, potentially implicated in pathogenesis, persistence, and protective immunity, using a non-targeted approach.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

Evaluation of genetic demultiplexing of single-cell sequencing data from model species

Single-cell sequencing (sc-seq) provides a species agnostic tool to study cellular processes. However, these technologies are expensive and require sufficient cell quantities and biological replicates to avoid artifactual results. An option to address these problems is pooling cells from multiple individuals into one sc-seq library. In humans, genotype-based computational separation (i.e., demultiplexing) of pooled sc-seq samples is common. This approach would be instrumental for studying non-isogenic model organisms. We set out to determine whether genotype-based demultiplexing could be more broadly applied among species ranging from zebrafish to non-human primates. Using such non-isogenic species, we benchmark genotype-based demultiplexing of pooled sc-seq datasets against various ground truths. We demonstrate that genotype-based demultiplexing of pooled sc-seq samples can be used with confidence in several non-isogenic model organisms and uncover limitations of this method. Importantly, the only genomic resource required for this approach is sc-seq data and a de novo transcriptome. The incorporation of pooling into sc-seq study designs will decrease cost while simultaneously increasing the reproducibility and experimental options in non-isogenic model organisms.

Identification of gametocyte-associated pir genes in the rodent malaria parasite, Plasmodium chabaudi chabaudi AS

OBJECTIVE

To analyse the transcriptional profiles of the pir multigene family of Plasmodium chabaudi chabaudi in male and female gametocytes isolated from the blood of infected mice.

RESULTS

Infected red blood cells containing female and male P. chabaudi gametocytes transcribe a distinct set of genes encoded by the multigene family pir. The overall patterns are similar to what has been observed in the close relative P. berghei, but here we show that gametocyte-associated pir genes are distinct from those involved in chronic blood-stage infection and highlight a male-associated pir gene which should be the focus of future studies.

A Microtubule-Associated Protein Is Essential for Malaria Parasite Transmission

Mature gametocytes of Plasmodium falciparum display a banana (falciform) shape conferred by a complex array of subpellicular microtubules (SPMT) associated with the inner membrane complex (IMC). Microtubule-associated proteins (MAPs) define MT populations and modulate interaction with pellicular components. Several MAPs have been identified in Toxoplasma gondii, and homologues can be found in the genomes of Plasmodium species, but the function of these proteins for asexual and sexual development of malaria parasites is still unknown. Here, we identified a novel subpellicular MAP, termed SPM3, that is conserved within the genus Plasmodium, especially within the subgenus Laverania, but absent in other Apicomplexa. Conditional knockdown and targeted gene disruption of Pfspm3 in Plasmodium falciparum cause severe morphological defects during gametocytogenesis, leading to round, nonfalciform gametocytes with an aberrant SPMT pattern. In contrast, Pbspm3 knockout in Plasmodium berghei, a species with round gametocytes, caused no defect in gametocytogenesis, but sporozoites displayed an aberrant motility and a dramatic defect in invasion of salivary glands, leading to a decreased efficiency in transmission. Electron microscopy revealed a dissociation of the SPMT from the IMC in Pbspm3 knockout parasites, suggesting a function of SPM3 in anchoring MTs to the IMC. Overall, our results highlight SPM3 as a pellicular component with essential functions for malaria parasite transmission. IMPORTANCE A key structural feature driving the transition between different life cycle stages of the malaria parasite is the unique three-membrane pellicle, consisting of the parasite plasma membrane (PPM) and a double membrane structure underlying the PPM termed the inner membrane complex (IMC). Additionally, there are numerous linearly arranged intramembranous particles (IMPs) linked to the IMC, which likely link the IMC to the subpellicular microtubule cytoskeleton. Here, we identified, localized, and characterized a novel subpellicular microtubule-associated protein unique to the genus Plasmodium. The knockout of this protein in the human-pathogenic species P. falciparum resulted in malformed gametocytes and aberrant microtubules. We confirmed the microtubule association in the P. berghei rodent malaria homologue and show that its knockout results in a perturbed microtubule architecture, aberrant sporozoite motility, and decreased transmission efficiency.

Single-Cell Transcriptomics To Define Plasmodium falciparum Stage Transition in the Mosquito Midgut

Malaria inflicts the highest rate of morbidity and mortality among the vector-borne diseases. The dramatic bottleneck of parasite numbers that occurs in the gut of the obligatory mosquito vector provides a promising target for novel control strategies. Using single-cell transcriptomics, we analyzed Plasmodium falciparum development in the mosquito gut, from unfertilized female gametes through the first 20 h after blood feeding, including the zygote and ookinete stages. This study revealed the temporal gene expression of the ApiAP2 family of transcription factors and of parasite stress genes in response to the harsh environment of the mosquito midgut. Further, employing structural protein prediction analyses, we found several upregulated genes predicted to encode intrinsically disordered proteins (IDPs), a category of proteins known for their importance in regulation of transcription, translation, and protein-protein interactions. IDPs are known for their antigenic properties and may serve as suitable targets for antibody- or peptide-based transmission suppression strategies. In total, this study uncovers the P. falciparum transcriptome from early to late parasite development in the mosquito midgut, inside its natural vector, which provides an important resource for future malaria transmission-blocking initiatives. IMPORTANCE The malaria parasite Plasmodium falciparum causes more than half a million deaths per year. The current treatment regimen targets the symptom-causing blood stage inside the human host. However, recent incentives in the field call for novel interventions to block parasite transmission from humans to the mosquito vector. Therefore, we need to better understand the parasite biology during its development inside the mosquito, including a deeper understanding of the expression of genes controlling parasite progression during these stages. Here, we have generated single-cell transcriptome data, covering P. falciparum's development, from gamete to ookinete inside the mosquito midgut, uncovering previously untapped parasite biology, including a repertoire of novel biomarkers to be explored in future transmission-blocking efforts. We anticipate that our study provides an important resource, which can be further explored to improve our understanding of the parasite biology as well as aid in guiding future malaria intervention strategies.

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.

RefTM: reference-guided topic modeling of single-cell chromatin accessibility data

Single-cell analysis is a valuable approach for dissecting the cellular heterogeneity, and single-cell chromatin accessibility sequencing (scCAS) can profile the epigenetic landscapes for thousands of individual cells. It is challenging to analyze scCAS data, because of its high dimensionality and a higher degree of sparsity compared with scRNA-seq data. Topic modeling in single-cell data analysis can lead to robust identification of the cell types and it can provide insight into the regulatory mechanisms. Reference-guided approach may facilitate the analysis of scCAS data by utilizing the information in existing datasets. We present RefTM (Reference-guided Topic Modeling of single-cell chromatin accessibility data), which not only utilizes the information in existing bulk chromatin accessibility and annotated scCAS data, but also takes advantage of topic models for single-cell data analysis. RefTM simultaneously models: (1) the shared biological variation among reference data and the target scCAS data; (2) the unique biological variation in scCAS data; (3) other variations from known covariates in scCAS data.

scDCCA: deep contrastive clustering for single-cell RNA-seq data based on auto-encoder network

The advances in single-cell ribonucleic acid sequencing (scRNA-seq) allow researchers to explore cellular heterogeneity and human diseases at cell resolution. Cell clustering is a prerequisite in scRNA-seq analysis since it can recognize cell identities. However, the high dimensionality, noises and significant sparsity of scRNA-seq data have made it a big challenge. Although many methods have emerged, they still fail to fully explore the intrinsic properties of cells and the relationship among cells, which seriously affects the downstream clustering performance. Here, we propose a new deep contrastive clustering algorithm called scDCCA. It integrates a denoising auto-encoder and a dual contrastive learning module into a deep clustering framework to extract valuable features and realize cell clustering. Specifically, to better characterize and learn data representations robustly, scDCCA utilizes a denoising Zero-Inflated Negative Binomial model-based auto-encoder to extract low-dimensional features. Meanwhile, scDCCA incorporates a dual contrastive learning module to capture the pairwise proximity of cells. By increasing the similarities between positive pairs and the differences between negative ones, the contrasts at both the instance and the cluster level help the model learn more discriminative features and achieve better cell segregation. Furthermore, scDCCA joins feature learning with clustering, which realizes representation learning and cell clustering in an end-to-end manner. Experimental results of 14 real datasets validate that scDCCA outperforms eight state-of-the-art methods in terms of accuracy, generalizability, scalability and efficiency. Cell visualization and biological analysis demonstrate that scDCCA significantly improves clustering and facilitates downstream analysis for scRNA-seq data. The code is available at https://github.com/WJ319/scDCCA.

Single-cell transcriptomics unveils xylem cell development and evolution

BACKGROUND

Xylem, the most abundant tissue on Earth, is responsible for lateral growth in plants. Typical xylem has a radial system composed of ray parenchyma cells and an axial system of fusiform cells. In most angiosperms, fusiform cells comprise vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed by tracheids in other vascular plants such as gymnosperms. Little is known about the developmental programs and evolutionary relationships of these xylem cell types.

RESULTS

Through both single-cell and laser capture microdissection transcriptomic profiling, we determine the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Based on cross-species analyses of single-cell clusters and overlapping trajectories, we reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages, whereas the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed trait, exhibit strong transcriptomic similarity to vessel elements rather than libriform fibers.

CONCLUSIONS

This evo-devo framework provides a comprehensive understanding of the formation of xylem cell lineages across multiple plant species spanning over a hundred million years of evolutionary history.

Global gene expression of human malaria parasite liver stages throughout intrahepatocytic development

Plasmodium falciparum (Pf) is causing the greatest malaria burden, yet the liver stages (LS) of this most important parasite species have remained poorly studied. Here, we used a human liver-chimeric mouse model in combination with a novel fluorescent PfNF54 parasite line (PfNF54cspGFP) to isolate PfLS-infected hepatocytes and generate transcriptomes that cover the major LS developmental phases in human hepatocytes. RNA-seq analysis of early Pf LS trophozoites two days after infection, revealed a central role of translational regulation in the transformation of the extracellular invasive sporozoite into intracellular LS. The developmental time course gene expression analysis indicated that fatty acid biosynthesis, isoprenoid biosynthesis and iron metabolism are sustaining LS development along with amino acid metabolism and biosynthesis. Countering oxidative stress appears to play an important role during intrahepatic LS development. Furthermore, we observed expression of the variant PfEMP1 antigen-encoding var genes, and we confirmed expression of PfEMP1 protein during LS development. Transcriptome comparison of the late Pf liver stage schizonts with P. vivax (Pv) late liver stages revealed highly conserved gene expression profiles among orthologous genes. A notable difference however was the expression of genes regulating sexual stage commitment. While Pv schizonts expressed markers of sexual commitment, the Pf LS parasites were not sexually committed and showed expression of gametocytogenesis repression factors. Our results provide the first comprehensive gene expression profile of the human malaria parasite Pf LS isolated during in vivo intrahepatocytic development. This data will inform biological studies and the search for effective intervention strategies that can prevent infection.

Long read single cell RNA sequencing reveals the isoform diversity of Plasmodium vivax transcripts

Plasmodium vivax infections often consist of heterogenous populations of parasites at different developmental stages and with distinct transcriptional profiles, which complicates gene expression analyses. The advent of single cell RNA sequencing (scRNA-seq) enabled disentangling this complexity and has provided robust and stage-specific characterization of Plasmodium gene expression. However, scRNA-seq information is typically derived from the end of each mRNA molecule (usually the 3'-end) and therefore fails to capture the diversity in transcript isoforms documented in bulk RNA-seq data. Here, we describe the sequencing of scRNA-seq libraries using Pacific Biosciences (PacBio) chemistry to characterize full-length Plasmodium vivax transcripts from single cell parasites. Our results show that many P. vivax genes are transcribed into multiple isoforms, primarily through variations in untranslated region (UTR) length or splicing, and that the expression of many isoforms is developmentally regulated. Our findings demonstrate that long read sequencing can be used to characterize mRNA molecules at the single cell level and provides an additional resource to better understand the regulation of gene expression throughout the Plasmodium life cycle.

Molecular mechanisms of cellular quiescence in apicomplexan parasites

Quiescence is a reversible nonproliferative cellular state that allows organisms to persist through unfavorable conditions. Quiescence can be stimulated by a wide range of external or intrinsic factors. Cells undergo a coordinated molecular program to enter and exit from the quiescent state, which is governed by signaling, transcriptional and translational changes, epigenetic mechanisms, metabolic switches, and changes in cellular architecture. These mechanisms have been extensively studied in model organisms, and a growing number of studies have identified conserved mechanisms in apicomplexan parasites. Quiescence in the context of a parasitic infection has significant clinical impact: quiescent forms may underlie treatment failures, relapsing infections, and stress tolerance. Here, we review the latest understanding of quiescence in apicomplexa, synthesizing these studies to highlight conserved mechanisms, and identifying technologies to assist in further characterization of quiescence. Understanding conserved mechanisms of quiescence in apicomplexans will provide avenues for transmission prevention and radical cure of infections.

Repurposing of Plasmodium falciparum var genes beyond the blood stage

A commonly observed survival strategy in protozoan parasites is the sequential expression of clonally variant-surface antigens to avoid elimination by the host's immune response. In malaria-causing P. falciparum, the immunovariant erythrocyte-membrane protein-1 (PfEMP1) adhesin family, encoded by var genes, is responsible for both antigenic variation and cytoadherence of infected erythrocytes to the microvasculature. Until recently, the biological function of these variant genes was believed to be restricted to intraerythrocytic developmental stages. With the advent of new technologies, var gene expression has been confirmed in transmission and pre-erythrocytic stages. Here, we discuss how repurposing of var gene expression beyond chronic blood-stage infection may be critical for successful transmission.

A non-coding GWAS variant impacts anthracycline-induced cardiotoxic phenotypes in human iPSC-derived cardiomyocytes

Anthracyclines, widely used to treat breast cancer, have the potential for cardiotoxicity. We have previously identified and validated a germline single nucleotide polymorphism, rs28714259, associated with an increased risk of anthracycline-induced heart failure. We now provide insights into the mechanism by which rs28714259 might confer increased risk of cardiac damage. Using hiPSC-derived cardiomyocyte cell lines with either intrinsic polymorphism or CRISPR-Cas9-mediated deletion of rs28714259 locus, we demonstrate that glucocorticoid receptor signaling activated by dexamethasone pretreatment prior to doxorubicin exposure preserves cardiomyocyte viability and contractility in cardiomyocytes containing the major allele. Homozygous loss of the rs28714259 major allele diminishes dexamethasone's protective effect. We further demonstrate that the risk allele of rs28714259 disrupts glucocorticoid receptor and rs28714259 binding affinity. Finally, we highlight the activation of genes and pathways involved in cardiac hypertrophy signaling that are blocked by the risk allele, suggesting a decreased adaptive survival response to doxorubicin-related stress.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

Naturally Acquired Kelch13 Mutations in Plasmodium falciparum Strains Modulate In Vitro Ring-Stage Artemisinin-Based Drug Tolerance and Parasite Survival in Response to Hyperoxia

The ring-stage survival assay was utilized to assess the impact of physiological hyperoxic stress on dihydroartemisinin (DHA) tolerance for a panel of Plasmodium falciparum strains with and without Kelch13 mutations. Strains without naturally acquired Kelch13 mutations or the postulated genetic background associated with delayed parasite clearance time demonstrated reduced proliferation under hyperoxic conditions in the subsequent proliferation cycle. Dihydroartemisinin tolerance in three isolates with naturally acquired Kelch13 mutations but not two genetically manipulated laboratory strains was modulated by in vitro hyperoxic stress exposure of early-ring-stage parasites in the cycle before drug exposure. Reduced parasite tolerance to additional derivatives, including artemisinin, artesunate, and OZ277, was observed within the second proliferation cycle. OZ439 and epoxomicin completely prevented parasite survival under both hyperoxia and normoxic in vitro culture conditions, highlighting the unique relationship between DHA tolerance and Kelch13 mutation-associated genetic background. IMPORTANCE Artemisinin-based combination therapy (ACT) for treating malaria is under intense scrutiny following treatment failures in the Greater Mekong subregion of Asia. This is further compounded by the potential for extensive loss of life if treatment failures extend to the African continent. Although Plasmodium falciparum has become resistant to all antimalarial drugs, artemisinin "resistance" does not present in the same way as resistance to other antimalarial drugs. Instead, a partial resistance or tolerance is demonstrated, associated with the parasite's genetic profile and linked to a molecular marker referred to as K13. It is suggested that parasites may have adapted to drug treatment, as well as the presence of underlying population health issues such as hemoglobinopathies, and/or environmental pressures, resulting in parasite tolerance to ACT. Understanding parasite evolution and control of artemisinin tolerance will provide innovative approaches to mitigate the development of artemisinin tolerance and thereby artemisinin-based drug treatment failure and loss of life globally to malaria infections.

A novel deep learning-assisted hybrid network for plasmodium falciparum parasite mitochondrial proteins classification

Plasmodium falciparum is a parasitic protozoan that can cause malaria, which is a deadly disease. Therefore, the accurate identification of malaria parasite mitochondrial proteins is essential for understanding their functions and identifying novel drug targets. For classifying protein sequences, several adaptive statistical techniques have been devised. Despite significant gains, prediction performance is still constrained by the lack of appropriate feature descriptors and learning strategies in current systems. Moreover, good ground truth data is important for Artificial Intelligence (AI)-based models but there is a lack of that data in the literature. Therefore, in this work, we propose a novel hybrid network that combines 1D Convolutional Neural Network (CNN) and Bidirectional Gated Recurrent Unit (BGRU) to classify the malaria parasite mitochondrial proteins. Furthermore, we curate a sequential data that are collected from National Center for Biotechnology Information (NCBI) and UniProtKB/Swiss-Prot proteins databanks to prepare a dataset that can be used by the research community for AI-based algorithms evaluation. We obtain 4204 cases after preprocessing of the collected data and denote this set of proteins as PF4204. Finally, we conduct an ablation study on several conventional and deep models using PF4204 and the benchmark PF2095 datasets. The proposed model 'CNN-BGRU' obtains the accuracy values of 0.9096 and 0.9857 on PF4204 and PF2095 datasets, respectively. In addition, the CNN-BGRU is compared with state-of-the-arts, where the results illustrate that it can extract robust features and identify proteins accurately.

A reproducible and sensitive method for generating high-quality transcriptomes from single whitefly salivary glands and other low-input tissues

Transcriptomic studies are an important tool for understanding the molecular pathways underlying host plant use by agricultural pests, including vectors of damaging plant pathogens. Thus far, bulk RNA-Seq has been the main approach for non-model insects. This method relies on pooling large numbers of whole organisms or hundreds of individually dissected organs. The latter approach is logistically challenging, may introduce artifacts of handling and storage, and is not compatible with biological replication. Here, we tested an approach to generate transcriptomes of individual salivary glands and other low-input body tissues from whiteflies (Bemisia tabaci MEAM1), which are major vectors of plant viruses. By comparing our outputs to published bulk RNA-Seq datasets for whole whitefly bodies and pools of salivary glands, we demonstrate that this approach recovers similar numbers of transcripts relative to bulk RNA-Seq in a tissue-specific manner, and for some metrics, exceeds performance of bulk tissue RNA-Seq. Libraries generated from individual salivary glands also yielded additional novel transcripts not identified in pooled salivary gland datasets, and had hundreds of enriched transcripts when compared with whole head tissues. Overall, our study demonstrates that it is feasible to produce high quality, replicated transcriptomes of whitefly salivary glands and other low-input tissues. We anticipate that our approach will expand hypothesis-driven research on salivary glands of whiteflies and other Hemiptera, thus enabling novel control strategies to disrupt feeding and virus transmission.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

Single cell RNA sequencing reveals hemocyte heterogeneity in Biomphalaria glabrata: Plasticity over diversity

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.

Determination of the Stage Composition of Plasmodium Infections from Bulk Gene Expression Data

Malaria symptoms are caused by the development of the parasites within the blood of an infected host. Bulk RNA sequencing (RNA-seq) of infected blood can reveal interactions between parasites and the host immune system during an infection, but because multiple developmental stages with distinct transcriptional profiles are concurrently present in infected blood, it is necessary to correct such analyses for differences in cell composition among samples. Gene expression deconvolution is a statistical approach that has been developed for inferring the cell composition of complex tissues characterized by bulk RNA-seq using gene expression profiles from reference cell types. Here, we describe the evaluation of a species-agnostic reference data set that can be used for efficient and accurate gene expression deconvolution of bulk RNA-seq data generated from any Plasmodium species and for correct gene expression analyses for biases caused by differences in stage composition among samples. IMPORTANCE Differences in cell type proportions among samples can introduce artifacts in gene expression analyses and mask genuine differences in gene regulation. Gene expression deconvolution allows estimation of the proportion of each cell type present in one sample directly from bulk RNA sequencing data, but this approach requires a reference data set with the signature profile of each cell type. Here, we evaluate the suitability of a rodent malaria parasite gene expression data set for estimating the proportions of each parasite developmental stage present in bulk RNA sequencing data generated from blood-stage infections with the human parasites Plasmodium falciparum and Plasmodium vivax. These analyses provide a species-agnostic approach for reliably estimating stage proportions in infected human blood and correcting subsequent gene expression analyses for these variations.

Single-cell RNA profiling of Plasmodium vivax-infected hepatocytes reveals parasite- and host- specific transcriptomic signatures and therapeutic targets

The resilience of Plasmodium vivax, the most widely-distributed malaria-causing parasite in humans, is attributed to its ability to produce dormant liver forms known as hypnozoites, which can activate weeks, months, or even years after an initial mosquito bite. The factors underlying hypnozoite formation and activation are poorly understood, as is the parasite's influence on the host hepatocyte. Here, we shed light on transcriptome-wide signatures of both the parasite and the infected host cell by sequencing over 1,000 P. vivax-infected hepatocytes at single-cell resolution. We distinguish between replicating schizonts and hypnozoites at the transcriptional level, identifying key differences in transcripts encoding for RNA-binding proteins associated with cell fate. In infected hepatocytes, we show that genes associated with energy metabolism and antioxidant stress response are upregulated, and those involved in the host immune response downregulated, suggesting both schizonts and hypnozoites alter the host intracellular environment. The transcriptional markers in schizonts, hypnozoites, and infected hepatocytes revealed here pinpoint potential factors underlying dormancy and can inform therapeutic targets against P. vivax liver-stage infection.

Genomic epidemiological models describe pathogen evolution across fitness valleys

Genomics is fundamentally changing epidemiological research. However, systematically exploring hypotheses in pathogen evolution requires new modeling tools. Models intertwining pathogen epidemiology and genomic evolution can help understand processes such as the emergence of novel pathogen genotypes with higher transmissibility or resistance to treatment. In this work, we present Opqua, a flexible simulation framework that explicitly links epidemiology to sequence evolution and selection. We use Opqua to study determinants of evolution across fitness valleys. We confirm that competition can limit evolution in high-transmission environments and find that low transmission, host mobility, and complex pathogen life cycles facilitate reaching new adaptive peaks through population bottlenecks and decoupling of selective pressures. The results show the potential of genomic epidemiological modeling as a tool in infectious disease research.

Genome-wide functional analysis reveals key roles for kinesins in the mammalian and mosquito stages of the malaria parasite life cycle

Kinesins are microtubule (MT)-based motors important in cell division, motility, polarity, and intracellular transport in many eukaryotes. However, they are poorly studied in the divergent eukaryotic pathogens Plasmodium spp., the causative agents of malaria, which manifest atypical aspects of cell division and plasticity of morphology throughout the life cycle in both mammalian and mosquito hosts. Here, we describe a genome-wide screen of Plasmodium kinesins, revealing diverse subcellular locations and functions in spindle assembly, axoneme formation, and cell morphology. Surprisingly, only kinesin-13 is essential for growth in the mammalian host while the other 8 kinesins are required during the proliferative and invasive stages of parasite transmission through the mosquito vector. In-depth analyses of kinesin-13 and kinesin-20 revealed functions in MT dynamics during apical cell polarity formation, spindle assembly, and axoneme biogenesis. These findings help us to understand the importance of MT motors and may be exploited to discover new therapeutic interventions against malaria.

Differential Trafficking and Expression of PIR Proteins in Acute and Chronic Plasmodium Infections

Plasmodium multigene families are thought to play important roles in the pathogenesis of malaria. Plasmodium interspersed repeat (pir) genes comprise the largest multigene family in many Plasmodium species. However, their expression pattern and localisation remain to be elucidated. Understanding protein subcellular localisation is fundamental to reveal the functional importance and cell-cell interactions of the PIR proteins. Here, we use the rodent malaria parasite, Plasmodium chabaudi chabaudi, as a model to investigate the localisation pattern of this gene family. We found that most PIR proteins are co-expressed in clusters during acute and chronic infection; members of the S7 clade are predominantly expressed during the acute-phase, whereas members of the L1 clade dominate the chronic-phase of infection. Using peptide antisera specific for S7 or L1 PIRS, we show that these PIRs have different localisations within the infected red blood cells. S7 PIRs are exported into the infected red blood cell cytoplasm where they are co-localised with parasite-induced host cell modifications termed Maurer’s clefts, whereas L1 PIRs are localised on or close to the parasitophorous vacuolar membrane. This localisation pattern changes following mosquito transmission and during progression from acute- to chronic-phase of infection. The presence of PIRs in Maurer’s clefts, as seen for Plasmodium falciparum RIFIN and STEVOR proteins, might suggest trafficking of the PIRs on the surface of the infected erythrocytes. However, neither S7 nor L1 PIR proteins detected by the peptide antisera are localised on the surface of infected red blood cells, suggesting that they are unlikely to be targets of surface variant-specific antibodies or to be directly involved in adhesion of infected red blood cells to host cells, as described for Plasmodium falciparum VAR proteins. The differences in subcellular localisation of the two major clades of Plasmodium chabaudi PIRs across the blood cycle, and the apparent lack of expression on the red cell surface strongly suggest that the function(s) of this gene family may differ from those of other multigene families of Plasmodium, such as the var genes of Plasmodium falciparum.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Pervasive sequence-level variation in the transcriptome of Plasmodium falciparum

Single-nucleotide variations (SNVs) in RNA, arising from co- and post-transcriptional phenomena including transcription errors and RNA-editing, are well studied in a range of organisms. In the malaria parasite Plasmodium falciparum, stage-specific and non-specific gene-expression variations accompany the parasite's array of developmental and morphological phenotypes over the course of its complex life cycle. However, the extent, rate and effect of sequence-level variation in the parasite's transcriptome are unknown. Here, we report the presence of pervasive, non-specific SNVs in the P. falciparum transcriptome. SNV rates for a gene were correlated to gene length (r[Formula: see text]0.65-0.7) but not to the AT-content of that gene. Global SNV rates for the P. falciparum lines we used, and for publicly available P. vivax and P. falciparum clinical isolate datasets, were of the order of 10-3 per base, ∼10× higher than rates we calculated for bacterial datasets. These variations may reflect an intrinsic transcriptional error rate in the parasite, and RNA editing may be responsible for a subset of them. This seemingly characteristic property of the parasite may have implications for clinical outcomes and the basic biology and evolution of P. falciparum and parasite biology more broadly. We anticipate that our study will prompt further investigations into the exact sources, consequences and possible adaptive roles of these SNVs.

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.

Genetically identical cells can be phenotypically diverse to allow adaptive flexibility in a given environment. This phenotypic diversity is driven by epigenetic and transcriptional variability. Here, Tripathi et al. perform scRNA-seq of isogenic and non-isogenic Plasmodium falciparum schizont populations to explore transcriptional heterogeneity and stochastic gene expression during the course of development.

Non-coding RNAs in malaria infection

Malaria is one of the most severe infectious diseases affecting humans and it is caused by protozoan pathogens of the species Plasmodium (spp.). The malaria parasite Plasmodium is characterized by a complex, multistage life cycle that requires tight gene regulation which allows for host invasion and defense against host immune responses. Unfortunately, the mechanisms regulating gene expression during Plasmodium infection remain largely elusive, though several lines of evidence implicate a major involvement of non-coding RNAs (ncRNAs). The ncRNAs have been found to play a key role in regulating transcriptional and post-transcriptional events in a broad range of organisms including Plasmodium. In Plasmodium ncRNAs have been shown to regulate key events in the multistage life cycle and virulence ability. Here we review recent progress involving ncRNAs (microRNAs, long non-coding RNAs, and circular RNAs) and their role as regulators of gene expression during Plasmodium infection in human hosts with focus on the possibility of using these molecules as biomarkers for monitoring disease status. We also discuss the surprising function of ncRNAs in mediating the complex interplay between parasite and human host and future perspectives of the field. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

Single-cell transcriptomics reveals expression profiles of Trypanosoma brucei sexual stages

Early diverging lineages such as trypanosomes can provide clues to the evolution of sexual reproduction in eukaryotes. In Trypanosoma brucei, the pathogen that causes Human African Trypanosomiasis, sexual reproduction occurs in the salivary glands of the insect host, but analysis of the molecular signatures that define these sexual forms is complicated because they mingle with more numerous, mitotically-dividing developmental stages. We used single-cell RNA-sequencing (scRNAseq) to profile 388 individual trypanosomes from midgut, proventriculus, and salivary glands of infected tsetse flies allowing us to identify tissue-specific cell types. Further investigation of salivary gland parasite transcriptomes revealed fine-scale changes in gene expression over a developmental progression from putative sexual forms through metacyclics expressing variant surface glycoprotein genes. The cluster of cells potentially containing sexual forms was characterized by high level transcription of the gamete fusion protein HAP2, together with an array of surface proteins and several genes of unknown function. We linked these expression patterns to distinct morphological forms using immunofluorescence assays and reporter gene expression to demonstrate that the kinetoplastid-conserved gene Tb927.10.12080 is exclusively expressed at high levels by meiotic intermediates and gametes. Further experiments are required to establish whether this protein, currently of unknown function, plays a role in gamete formation and/or fusion.

ROSIE: RObust Sparse ensemble for outlIEr detection and gene selection in cancer omics data

The extraction of novel information from omics data is a challenging task, in particular, since the number of features (e.g. genes) often far exceeds the number of samples. In such a setting, conventional parameter estimation leads to ill-posed optimization problems, and regularization may be required. In addition, outliers can largely impact classification accuracy.

Here we introduce ROSIE, an ensemble classification approach, which combines three sparse and robust classification methods for outlier detection and feature selection and further performs a bootstrap-based validity check. Outliers of ROSIE are determined by the rank product test using outlier rankings of all three methods, and important features are selected as features commonly selected by all methods.

We apply ROSIE to RNA-Seq data from The Cancer Genome Atlas (TCGA) to classify observations into Triple-Negative Breast Cancer (TNBC) and non-TNBC tissue samples. The pre-processed dataset consists of 16,600 genes and more than 1,000 samples. We demonstrate that ROSIE selects important features and outliers in a robust way. Identified outliers are concordant with the distribution of the commonly selected genes by the three methods, and results are in line with other independent studies. Furthermore, we discuss the association of some of the selected genes with the TNBC subtype in other investigations. In summary, ROSIE constitutes a robust and sparse procedure to identify outliers and important genes through binary classification. Our approach is ad hoc applicable to other datasets, fulfilling the overall goal of simultaneously identifying outliers and candidate disease biomarkers to the targeted in therapy research and personalized medicine frameworks.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

A transcriptional switch controls sex determination in Plasmodium falciparum

Sexual reproduction and meiotic sex are deeply rooted in the eukaryotic tree of life, but mechanisms determining sex or mating types are extremely varied and are only well characterized in a few model organisms¹. In malaria parasites, sexual reproduction coincides with transmission to the vector host. Sex determination is non-genetic, with each haploid parasite capable of producing either a male or a female gametocyte in the human host². The hierarchy of events and molecular mechanisms that trigger sex determination and maintenance of sexual identity are yet to be elucidated. Here we show that the male development 1 (md1) gene is both necessary and sufficient for male fate determination in the human malaria parasite Plasmodium falciparum. We show that Md1 has a dual function stemming from two separate domains: in sex determination through its N terminus and in male development from its conserved C-terminal LOTUS/OST-HTH domain. We further identify a bistable switch at the md1 locus, which is coupled with sex determination and ensures that the male-determining gene is not expressed in the female lineage. We describe one of only a few known non-genetic mechanisms of sex determination in a eukaryote and highlight Md1 as a potential target for interventions that block malaria transmission.

Analysis of pir gene expression across the Plasmodium life cycle

Background

Plasmodium interspersed repeat (pir) is the largest multigene family in the genomes of most Plasmodium species. A variety of functions for the PIR proteins which they encode have been proposed, including antigenic variation, immune evasion, sequestration and rosetting. However, direct evidence for these is lacking. The repetitive nature of the family has made it difficult to determine function experimentally. However, there has been some success in using gene expression studies to suggest roles for some members in virulence and chronic infection.

Methods

Here pir gene expression was examined across the life cycle of Plasmodium berghei using publicly available RNAseq data-sets, and at high resolution in the intraerythrocytic development cycle using new data from Plasmodium chabaudi.

Results

Expression of pir genes is greatest in stages of the parasite which invade and reside in red blood cells. The marked exception is that liver merozoites and male gametocytes produce a very large number of pir gene transcripts, notably compared to female gametocytes, which produce relatively few. Within the asexual blood stages different subfamilies peak at different times, suggesting further functional distinctions. Representing a subfamily of its own, the highly conserved ancestral pir gene warrants further investigation due to its potential tractability for functional investigation. It is highly transcribed in multiple life cycle stages and across most studied Plasmodium species and thus is likely to play an important role in parasite biology.

Conclusions

The identification of distinct expression patterns for different pir genes and subfamilies is likely to provide a basis for the design of future experiments to uncover their function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03979-6.

A cell-cell atlas approach for understanding symbiotic interactions between microbes

Natural environments are composed of a huge diversity of microorganisms interacting with each other to form complex functional networks. Our understanding of the operative nature of host-symbiont associations is limited because propagating such associations in a laboratory is challenging. The advent of single-cell technologies applied to, for example, animal cells and apicomplexan parasites has revolutionized our understanding of development and disease. Such cell atlas approaches generate maps of cell-specific processes and variations within cellular populations. These methods can now be combined with cellular-imaging so that interaction stage versus transcriptome state can be quantized for microbe-microbe interactions. We predict that the combination of these methods applied to the study of symbioses will transform our understanding of many ecological interactions, including those sampled directly from natural environments.