Transcriptional control of the Cryptosporidium life cycle.

The parasite Cryptosporidium is a leading agent of diarrhoeal disease in young children, and a cause and consequence of chronic malnutrition1,2. There are no vaccines and only limited treatment options3. The parasite infects enterocytes, in which it engages in asexual and sexual replication4, both of which are essential to continued infection and transmission. However, their molecular mechanisms remain largely unclear5. Here we use single-cell RNA sequencing to reveal the gene expression programme of the entire Cryptosporidium parvum life cycle in culture and in infected animals. Diverging from the prevailing model6, we find support for only three intracellular stages: asexual type-I meronts, male gamonts and female gametes. We reveal a highly organized program for the assembly of components at each stage. Dissecting the underlying regulatory network, we identify the transcription factor Myb-M as the earliest determinant of male fate, in an organism that lacks genetic sex determination. Conditional expression of this factor overrides the developmental program and induces widespread maleness, while conditional deletion ablates male development. Both have a profound impact on the infection. A large set of stage-specific genes now provides the opportunity to understand, engineer and disrupt parasite sex and life cycle progression to advance the development of vaccines and treatments.

Analysis of intestinal epithelial cell responses to Cryptosporidium highlights the temporal effects of IFN-γ on parasite restriction.

The production of IFN-γ is crucial for control of multiple enteric infections, but its impact on intestinal epithelial cells (IEC) is not well understood. Cryptosporidium parasites exclusively infect epithelial cells and the ability of interferons to activate the transcription factor STAT1 in IEC is required for parasite clearance. Here, the use of single cell RNA sequencing to profile IEC during infection revealed an increased proportion of mid-villus enterocytes during infection and induction of IFN-γ-dependent gene signatures that was comparable between uninfected and infected cells. These analyses were complemented by in vivo studies, which demonstrated that IEC expression of the IFN-γ receptor was required for parasite control. Unexpectedly, treatment of Ifng-/- mice with IFN-γ showed the IEC response to this cytokine correlates with a delayed reduction in parasite burden but did not affect parasite development. These data sets provide insight into the impact of IFN-γ on IEC and suggest a model in which IFN-γ signalling to uninfected enterocytes is important for control of Cryptosporidium.

Latent transcriptional variations of individual Plasmodium falciparum uncovered by single-cell RNA-seq and fluorescence imaging.

Malaria parasites follow a complex life cycle that consists of multiple stages that span from the human host to the mosquito vector. Among the species causing malaria, Plasmodium falciparum is the most lethal, with clinical symptoms manifesting during the intraerythrocytic developmental cycle (IDC). During the IDC, P. falciparum progresses through a synchronous and continuous cascade of transcriptional programming previously established using population analyses. While individual parasites are known to exhibit transcriptional variations to evade the host immune system or commit to a sexual fate, such rare expression heterogeneity is largely undetectable on a population level. Therefore, we combined single-cell RNA-sequencing (scRNA-seq) on a microfluidic platform and fluorescence imaging to delineate the transcriptional variations among individual parasites during late asexual and sexual stages. The comparison between asexual and sexual parasites uncovered a set of previously undefined sex-specific genes. Asexual parasites were segregated into three distinct clusters based on the differential expression of genes encoding SERAs, rhoptry proteins, and EXP2 plus transporters. Multiple pseudotime analyses revealed that these stage-specific transitions are distinct. RNA fluorescent in situ hybridization of cluster-specific genes validated distinct stage-specific expression and transitions during the IDC and defined the highly variable transcriptional pattern of EXP2. Additionally, these analyses indicated huge variations in the stage-specific transcript levels among parasites. Overall, scRNA-seq and RNA-FISH of P. falciparum revealed distinct stage transitions and unexpected degrees of heterogeneity with potential impact on transcriptional regulation during the IDC and adaptive responses to the host.

GNA13 loss in germinal center B cells leads to impaired apoptosis and promotes lymphoma in vivo.

GNA13 is the most frequently mutated gene in germinal center (GC)-derived B-cell lymphomas, including nearly a quarter of Burkitt lymphoma and GC-derived diffuse large B-cell lymphoma. These mutations occur in a pattern consistent with loss of function. We have modeled the GNA13-deficient state exclusively in GC B cells by crossing the Gna13 conditional knockout mouse strain with the GC-specific AID-Cre transgenic strain. AID-Cre(+) GNA13-deficient mice demonstrate disordered GC architecture and dark zone/light zone distribution in vivo, and demonstrate altered migration behavior, decreased levels of filamentous actin, and attenuated RhoA activity in vitro. We also found that GNA13-deficient mice have increased numbers of GC B cells that display impaired caspase-mediated cell death and increased frequency of somatic hypermutation in the immunoglobulin VH locus. Lastly, GNA13 deficiency, combined with conditional MYC transgene expression in mouse GC B cells, promotes lymphomagenesis. Thus, GNA13 loss is associated with GC B-cell persistence, in which impaired apoptosis and ongoing somatic hypermutation may lead to an increased risk of lymphoma development.

Trans-kingdom small RNA transfer during host-pathogen interactions: The case of P. falciparum and erythrocytes.

This review focuses on the role of trans-kingdom movement of small RNA (sRNA) molecules between parasites, particularly Plasmodium falciparum, and their respective host cells. While the intercellular transfer of sRNAs within organisms is well recognized, recent studies illustrate many examples of trans-kingdom sRNA exchange within the context of host-parasite interactions. These interactions are predominantly found in the transfer of host sRNAs between erythrocytes and the invading P. falciparum, as well as other host cell types. In addition, parasite-encoded sRNAs can also be transferred to host cells to evade the immune system. The transport of these parasite sRNAs in the body fluids of the host may also offer means to detect and monitor the parasite infection. These isolated examples may only represent the tip of the iceberg in which the transfer of sRNA between host and parasites is a critical aspect of host-pathogen interactions. In addition, the levels of these sRNAs and their speed of transfer may vary dramatically under different contexts to push the biologic equilibrium toward the benefit of hosts vs. parasites. Therefore, these sRNA transfers may offer potential strategies to detect, prevent or treat parasite infections. Here, we review a brief history of the discovery of host erythrocyte sRNAs, their transfers and interactions in the context of P. falciparum infection. We also provide examples and discuss the functional significance of the reciprocal transfer of parasite-encoded sRNAs into hosts. These understandings of sRNA exchanges are put in the context of their implications for parasite pathogenesis, host defenses and the evolution of host polymorphisms driven by host interactions with these parasites.

Hammondia hammondi, an avirulent relative of Toxoplasma gondii, has functional orthologs of known T. gondii virulence genes.

Toxoplasma gondii is a ubiquitous protozoan parasite capable of infecting all warm-blooded animals, including humans. Its closest extant relative, Hammondia hammondi, has never been found to infect humans and, in contrast to T. gondii, is highly attenuated in mice. To better understand the genetic bases for these phenotypic differences, we sequenced the genome of a H. hammondi isolate (HhCatGer041) and found the genomic synteny between H. hammondi and T. gondii to be >95%. We used this genome to determine the H. hammondi primary sequence of two major T. gondii mouse virulence genes, TgROP5 and TgROP18. When we expressed these genes in T. gondii, we found that H. hammondi orthologs of TgROP5 and TgROP18 were functional. Similar to T. gondii, the HhROP5 locus is expanded, and two distinct HhROP5 paralogs increased the virulence of a T. gondii TgROP5 knockout strain. We also identified a 107 base pair promoter region, absent only in type III TgROP18, which is necessary for TgROP18 expression. This result indicates that the ROP18 promoter was active in the most recent common ancestor of these two species and that it was subsequently inactivated in progenitors of the type III lineage. Overall, these data suggest that the virulence differences between these species are not solely due to the functionality of these key virulence factors. This study provides evidence that other mechanisms, such as differences in gene expression or the lack of currently uncharacterized virulence factors, may underlie the phenotypic differences between these species.

Hammondia hammondi harbors functional orthologs of the host-modulating effectors GRA15 and ROP16 but is distinguished from Toxoplasma gondii by a unique transcriptional profile.

Toxoplasma gondii and its nearest extant relative, Hammondia hammondi, are phenotypically distinct despite their remarkable similarity in gene content, synteny, and functionality. To begin to identify genetic differences that might drive distinct infection phenotypes of T. gondii and H. hammondi, in the present study we (i) determined whether two known host-interacting proteins, dense granule protein 15 (GRA15) and rhoptry protein 16 (ROP16), were functionally conserved in H. hammondi and (ii) performed the first comparative transcriptional analysis of H. hammondi and T. gondii sporulated oocysts. We found that GRA15 and ROP16 from H. hammondi (HhGRA15 and HhROP16) modulate the host NF-κB and STAT6 pathways, respectively, when expressed heterologously in T. gondii. We also found the transcriptomes of H. hammondi and T. gondii to be highly distinct. Consistent with the spontaneous conversion of H. hammondi tachyzoites into bradyzoites both in vitro and in vivo, H. hammondi high-abundance transcripts are enriched for genes that are of greater abundance in T. gondii bradyzoites. We also identified genes that are of high transcript abundance in H. hammondi but are poorly expressed in multiple T. gondii life stages, suggesting that these genes are uniquely expressed in H. hammondi. Taken together, these data confirm the functional conservation of known T. gondii virulence effectors in H. hammondi and point to transcriptional differences as a potential source of the phenotypic differences between these species.

Genetic Ablation of a Female-Specific Apetala 2 Transcription Factor Blocks Oocyst Shedding in Cryptosporidium parvum.

The apicomplexan parasite Cryptosporidium is a leading global cause of diarrheal disease, and the infection poses a particularly grave threat to young children and those with weakened immune function. Infection occurs by ingestion of meiotic spores called oocysts, and transmission relies on fecal shedding of new oocysts. The entire life cycle thus occurs in a single host and features asexual as well as sexual forms of replication. Here, we identify and locus tag two Apetala 2-type (AP2) transcription factors and demonstrate that they are exclusively expressed in male and female gametes, respectively. To enable functional studies of essential genes in Cryptosporidium parvum, we develop and validate a small-molecule-inducible gene excision system, which we apply to the female factor AP2-F to achieve conditional gene knockout. Analyzing this mutant, we find the factor to be dispensable for asexual growth and early female fate determination in vitro but to be required for oocyst shedding in infected animals in vivo. Transcriptional analyses conducted in the presence or absence of AP2-F revealed that the factor controls the transcription of genes encoding crystalloid body proteins, which are exclusively expressed in female gametes. In C. parvum, the organelle is restricted to sporozoites, and its loss in other apicomplexan parasites leads to blocked transmission. Overall, our development of conditional gene ablation in C. parvum provides a robust method for genetic analysis in this parasite that enabled us to identify AP2-F as an essential regulator of transcription required for oocyst shedding and transmission. IMPORTANCE The parasite Cryptosporidium infects millions of people worldwide each year, leading to life-threatening diarrheal disease in young children and immunosuppressed individuals. There is no vaccine and only limited treatment. Transmission occurs via the fecal-oral route by an environmentally resilient spore-like oocyst. Infection takes place in the intestinal epithelium, where parasites initially propagate asexually before transitioning to male and female gametes, with sex leading to the formation of new oocysts. The essential role of sexual development for continuous infection and transmission makes it an attractive target for therapy and prevention. To study essential genes and potential drug targets across the life cycle, we established inducible gene excision for C. parvum. We determined that the female-specific transcription factor AP2-F is not required for asexual growth and early female development in vitro but is necessary for oocyst shedding in vivo. This work enhances the genetic tools available to study Cryptosporidium gene function.

Analysis of intestinal epithelial cell responses to Cryptosporidium highlights the temporal effects of IFN-γ on parasite restriction.

The production of IFN-γ is crucial for control of multiple enteric infections, but its impact on intestinal epithelial cells (IEC) is not well understood. Cryptosporidium parasites exclusively infect epithelial cells and the ability of interferons to activate the transcription factor STAT1 in IEC is required for parasite clearance. The use of single cell RNA sequencing to profile IEC during infection revealed induction of IFN-γ-dependent gene signatures that was comparable between uninfected and infected cells, and IEC expression of the IFN-γ receptor was required for parasite control. Unexpectedly, treatment of Ifng-/- mice with IFN-γ demonstrated the IEC response to this cytokine correlates with a delayed reduction in parasite burden but did not affect parasite development. These data sets provide insight into the impact of IFN-γ on IEC and suggest a model in which IFN-γ-mediated bystander activation of uninfected enterocytes is important for control of Cryptosporidium.

Single-Cell Analysis Reveals Distinct Gene Expression and Heterogeneity in Male and Female Plasmodium falciparum Gametocytes.

Sexual reproduction is an obligate step in the Plasmodium falciparum life cycle, with mature gametocytes being the only form of the parasite capable of human-to-mosquito transmission. Development of male and female gametocytes takes 9 to 12 days, and although more than 300 genes are thought to be specific to gametocytes, only a few have been postulated to be male or female specific. Because these genes are often expressed during late gametocyte stages and for some, male- or female-specific transcript expression is debated, the separation of male and female populations is technically challenging. To overcome these challenges, we have developed an unbiased single-cell approach to determine which transcripts are expressed in male versus female gametocytes. Using microfluidic technology, we isolated single mid- to late-stage gametocytes to compare the expression of 91 genes, including 87 gametocyte-specific genes, in 90 cells. Such analysis identified distinct gene clusters whose expression was associated with male, female, or all gametocytes. In addition, a small number of male gametocytes clustered separately from female gametocytes based on sex-specific expression independent of stage. Many female-enriched genes also exhibited stage-specific expression. RNA fluorescent in situ hybridization of male and female markers validated the mutually exclusive expression pattern of male and female transcripts in gametocytes. These analyses uncovered novel male and female markers that are expressed as early as stage III gametocytogenesis, providing further insight into Plasmodium sex-specific differentiation previously masked in population analyses. Our single-cell approach reveals the most robust markers for sex-specific differentiation in Plasmodium gametocytes. Such single-cell expression assays can be generalized to all eukaryotic pathogens.IMPORTANCE Most human deaths that result from malaria are caused by the eukaryotic parasite Plasmodium falciparum The only form of this parasite that is transmitted to the mosquito is the sexual form, called the gametocyte. The production of mature gametocytes can take up to 2 weeks and results in phenotypically distinct males and females, although what causes this gender-specific differentiation remains largely unknown. Here, we demonstrate the first use of microfluidic technology to capture single gametocytes and determine their temporal sex-specific gene expression in an unbiased manner. We were able to determine male or female identity of single cells based on the upregulation of gender-specific genes as early as mid-stage gametocytes. This analysis has revealed strong markers for male and female gametocyte differentiation that were previously concealed in population analyses. Similar single-cell analyses in eukaryotic pathogens using this method may uncover rare cell types and heterogeneity previously masked in population studies.

Automated Detection of P. falciparum Using Machine Learning Algorithms with Quantitative Phase Images of Unstained Cells.

Malaria detection through microscopic examination of stained blood smears is a diagnostic challenge that heavily relies on the expertise of trained microscopists. This paper presents an automated analysis method for detection and staging of red blood cells infected by the malaria parasite Plasmodium falciparum at trophozoite or schizont stage. Unlike previous efforts in this area, this study uses quantitative phase images of unstained cells. Erythrocytes are automatically segmented using thresholds of optical phase and refocused to enable quantitative comparison of phase images. Refocused images are analyzed to extract 23 morphological descriptors based on the phase information. While all individual descriptors are highly statistically different between infected and uninfected cells, each descriptor does not enable separation of populations at a level satisfactory for clinical utility. To improve the diagnostic capacity, we applied various machine learning techniques, including linear discriminant classification (LDC), logistic regression (LR), and k-nearest neighbor classification (NNC), to formulate algorithms that combine all of the calculated physical parameters to distinguish cells more effectively. Results show that LDC provides the highest accuracy of up to 99.7% in detecting schizont stage infected cells compared to uninfected RBCs. NNC showed slightly better accuracy (99.5%) than either LDC (99.0%) or LR (99.1%) for discriminating late trophozoites from uninfected RBCs. However, for early trophozoites, LDC produced the best accuracy of 98%. Discrimination of infection stage was less accurate, producing high specificity (99.8%) but only 45.0%-66.8% sensitivity with early trophozoites most often mistaken for late trophozoite or schizont stage and late trophozoite and schizont stage most often confused for each other. Overall, this methodology points to a significant clinical potential of using quantitative phase imaging to detect and stage malaria infection without staining or expert analysis.

Hemoglobin consumption by P. falciparum in individual erythrocytes imaged via quantitative phase spectroscopy.

Plasmodium falciparum infection causes structural and biochemical changes in red blood cells (RBCs). To quantify these changes, we apply a novel optical technique, quantitative phase spectroscopy (QPS) to characterize individual red blood cells (RBCs) during the intraerythrocytic life cycle of P. falciparum. QPS captures hyperspectral holograms of individual RBCs to measure spectroscopic changes across the visible wavelength range (475-700 nm), providing complex information, i.e. amplitude and phase, about the light field which has interacted with the cell. The complex field provides complimentary information on hemoglobin content and cell mass, which are both found to dramatically change upon infection by P. falciparum. Hb content progressively decreases with parasite life cycle, with an average 72.2% reduction observed for RBCs infected by schizont-stage P. falciparum compared to uninfected cells. Infection also resulted in a 33.1% reduction in RBC's optical volume, a measure of the cells' non-aqueous components. Notably, optical volume is only partially correlated with hemoglobin content, suggesting that changes in other dry mass components such as parasite mass may also be assessed using this technique. The unique ability of QPS to discriminate individual healthy and infected cells using spectroscopic changes indicates that the approach can be used to detect disease.

Methods to Investigate the Regulatory Role of Small RNAs and Ribosomal Occupancy of Plasmodium falciparum.

The genetic variation responsible for the sickle cell allele (HbS) enables erythrocytes to resist infection by the malaria parasite, P. falciparum. The molecular basis of this resistance, which is known to be multifactorial, remains incompletely understood. Recent studies found that the differential expression of erythrocyte microRNAs, once translocated into malaria parasites, affect both gene regulation and parasite growth. These miRNAs were later shown to inhibit mRNA translation by forming a chimeric RNA transcript via 5' RNA fusion with discreet subsets of parasite mRNAs. Here, the techniques that were used to study the functional role and putative mechanism underlying erythrocyte microRNAs on the gene regulation and translational potential of P. falciparum, including the transfection of modified synthetic microRNAs into host erythrocytes, will be detailed.  Finally, a polysome gradient method is used to determine the extent of translation of these transcripts. Together, these techniques allowed us to demonstrate that the dysregulated levels of erythrocyte microRNAs contribute to cell-intrinsic malaria resistance of sickle erythrocytes.

A single chromosome unexpectedly links highly divergent isolates of Toxoplasma gondii.

Toxoplasma gondii is an obligate intracellular parasite that can cause disease in all warm-blooded animals studied to date, including humans. Over a billion people have been infected with this parasite worldwide. In Europe and North America, Toxoplasma has a clonal population structure, where only three lineages are highly dominant (strain types I, II, and III). Khan et al. [mBio 2(6): e00228-11, 2011] have carried out phylogenetic analyses on a large number of diverse strains from outside of these lineages and found evidence for a significant split between the clonal North American/European lineages and those in South America. In contrast to most of the genome, nearly all North American/European strains sampled, and the majority of South American strains sampled, harbored at least portions of a monomorphic chromosome Ia (Ia*). In contrast to previous models, these data suggest that the monomorphic haplotype originated in South America and migrated to the North. These authors propose that South American haplotype 12 was a precursor to modern-day type II, while South American haplotypes 6 and 9 crossed with haplotype 12 to give rise to the type I and III lineages, respectively. However, the findings reported by Khan et al. complicate the origin of chromosome Ia, since there are members of haplotypes 9 and 12 with nearly complete versions of Ia* and members of haplotypes 6 and 12 with over 50% of Ia*. This unexpected finding raises exciting new questions about how an entire common chromosome can be found within strains that are highly divergent at most other genomic loci.