Representing and extracting knowledge from single-cell data.

Single-cell analysis is currently one of the most high-resolution techniques to study biology. The large complex datasets that have been generated have spurred numerous developments in computational biology, in particular the use of advanced statistics and machine learning. This review attempts to explain the deeper theoretical concepts that underpin current state-of-the-art analysis methods. Single-cell analysis is covered from cell, through instruments, to current and upcoming models. The aim of this review is to spread concepts which are not yet in common use, especially from topology and generative processes, and how new statistical models can be developed to capture more of biology. This opens epistemological questions regarding our ontology and models, and some pointers will be given to how natural language processing (NLP) may help overcome our cognitive limitations for understanding single-cell data.

Dissecting the genetic landscape of GPCR signaling through phenotypic profiling in C. elegans.

G protein-coupled receptors (GPCRs) mediate responses to various extracellular and intracellular cues. However, the large number of GPCR genes and their substantial functional redundancy make it challenging to systematically dissect GPCR functions in vivo. Here, we employ a CRISPR/Cas9-based approach, disrupting 1654 GPCR-encoding genes in 284 strains and mutating 152 neuropeptide-encoding genes in 38 strains in C. elegans. These two mutant libraries enable effective deorphanization of chemoreceptors, and characterization of receptors for neuropeptides in various cellular processes. Mutating a set of closely related GPCRs in a single strain permits the assignment of functions to GPCRs with functional redundancy. Our analyses identify a neuropeptide that interacts with three receptors in hypoxia-evoked locomotory responses, unveil a collection of regulators in pathogen-induced immune responses, and define receptors for the volatile food-related odorants. These results establish our GPCR and neuropeptide mutant libraries as valuable resources for the C. elegans community to expedite studies of GPCR signaling in multiple contexts.

scDual-Seq of Toxoplasma gondii-infected mouse BMDCs reveals heterogeneity and differential infection dynamics.

Dendritic cells and macrophages are integral parts of the innate immune system and gatekeepers against infection. The protozoan pathogen, Toxoplasma gondii, is known to hijack host immune cells and modulate their immune response, making it a compelling model to study host-pathogen interactions. Here we utilize single cell Dual RNA-seq to parse out heterogeneous transcription of mouse bone marrow-derived dendritic cells (BMDCs) infected with two distinct genotypes of T. gondii parasites, over multiple time points post infection. We show that the BMDCs elicit differential responses towards T. gondii infection and that the two parasite lineages distinctly manipulate subpopulations of infected BMDCs. Co-expression networks define host and parasite genes, with implications for modulation of host immunity. Integrative analysis validates previously established immune pathways and additionally, suggests novel candidate genes involved in host-pathogen interactions. Altogether, this study provides a comprehensive resource for characterizing host-pathogen interplay at high-resolution.

Type I interferon shapes brain distribution and tropism of tick-borne flavivirus.

Viral tropism within the brain and the role(s) of vertebrate immune response to neurotropic flaviviruses infection is largely understudied. We combine multimodal imaging (cm-nm scale) with single nuclei RNA-sequencing to study Langat virus in wildtype and interferon alpha/beta receptor knockout (Ifnar-/-) mice to visualize viral pathogenesis and define molecular mechanisms. Whole brain viral infection is imaged by Optical Projection Tomography coregistered to ex vivo MRI. Infection is limited to grey matter of sensory systems in wildtype mice, but extends into white matter, meninges and choroid plexus in Ifnar-/- mice. Cells in wildtype display strong type I and II IFN responses, likely due to Ifnb expressing astrocytes, infiltration of macrophages and Ifng-expressing CD8+ NK cells, whereas in Ifnar-/-, the absence of this response contributes to a shift in cellular tropism towards non-activated resident microglia. Multimodal imaging-transcriptomics exemplifies a powerful way to characterize mechanisms of viral pathogenesis and tropism.

Immune-Activated B Cells Are Dominant in Prostate Cancer.

B cells are multifaceted immune cells responding robustly during immune surveillance against tumor antigens by presentation to T cells and switched immunoglobulin production. However, B cells are unstudied in prostate cancer (PCa). We used flow cytometry to analyze B-cell subpopulations in peripheral blood and lymph nodes from intermediate-high risk PCa patients. B-cell subpopulations were related to clinicopathological factors. B-cell-receptor single-cell sequencing and VDJ analysis identified clonal B-cell expansion in blood and lymph nodes. Pathological staging was pT2 in 16%, pT3a in 48%, and pT3b in 36%. Lymph node metastases occurred in 5/25 patients (20%). Compared to healthy donors, the peripheral blood CD19+ B-cell compartment was significantly decreased in PCa patients and dominated by naïve B cells. The nodal B-cell compartment had significantly increased fractions of CD19+ B cells and switched memory B cells. Plasmablasts were observed in tumor-draining sentinel lymph nodes (SNs). VDJ analysis revealed clonal expansion in lymph nodes. Thus, activated B cells are increased in SNs from PCa patients. The increased fraction of switched memory cells and plasmablasts together with the presence of clonally expanded B cells indicate tumor-specific T-cell-dependent responses from B cells, supporting an important role for B cells in the protection against tumors.

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.

Tuned activation of MSLN-CAR T cells induces superior antitumor responses in ovarian cancer models.

BACKGROUND: Limited persistence of functional CAR T cells in the immunosuppressive solid tumor microenvironment remains a major hurdle in the successful translation of CAR T cell therapy to treat solid tumors. Fine-tuning of CAR T cell activation by mutating CD3ζ chain immunoreceptor tyrosine-based activation motifs (ITAMs) in CD19-CAR T cells (containing the CD28 costimulatory domain) has proven to extend functional CAR T cell persistence in preclinical models of B cell malignancies. METHODS: In this study, two conventional second-generation MSLN-CAR T cell constructs encoding for either a CD28 co-stimulatory (M28z) or 4-1BB costimulatory (MBBz) domain and a novel mesothelin (MSLN)-directed CAR T cell construct encoding for the CD28 costimulatory domain and CD3ζ chain containing a single ITAM (M1xx) were evaluated using in vitro and in vivo preclinical models of ovarian cancer. Two ovarian cancer cell lines and two orthotopic models of ovarian cancer in NSG mice were used: SKOV-3 cells inoculated through microsurgery in the ovary and to mimic a disseminated model of advanced ovarian cancer, OVCAR-4 cells injected intraperitoneally. MSLN-CAR T cell treatment efficacy was evaluated by survival analysis and the characterization and quantification of the different MSLN-CAR T cells were performed by flow cytometry, quantitative PCR and gene expression analysis. RESULTS: M1xx CAR T cells elicited superior antitumor potency and persistence, as compared with the conventional second generation M28z and MBBz CAR T cells. Ex vivo M28z and MBBz CAR T cells displayed a more exhausted phenotype than M1xx CAR T cells as determined by co-expression of PD-1, LAG-3 and TIM-3. Furthermore, M1xx CAR T cells showed superior ex vivo IFNy, TNF and GzB production and were characterized by a self-renewal gene signature. CONCLUSIONS: Altogether, our study demonstrates the enhanced therapeutic potential of MSLN-CAR T cells expressing a mutated CD3ζ chain containing a single ITAM for the treatment of ovarian cancer. CAR T cells armored with calibrated activation potential may improve the clinical responses in solid tumors.

Serine Protease Inhibitors Restrict Host Susceptibility to SARS-CoV-2 Infections.

The coronavirus disease 2019, COVID-19, is a complex disease with a wide range of symptoms from asymptomatic infections to severe acute respiratory syndrome with lethal outcome. Individual factors such as age, sex, and comorbidities increase the risk for severe infections, but other aspects, such as genetic variations, are also likely to affect the susceptibility to SARS-CoV-2 infection and disease severity. Here, we used a human 3D lung cell model based on primary cells derived from multiple donors to identity host factors that regulate SARS-CoV-2 infection. With a transcriptomics-based approach, we found that less susceptible donors show a higher expression level of serine protease inhibitors SERPINA1, SERPINE1, and SERPINE2, identifying variation in cellular serpin levels as restricting host factors for SARS-CoV-2 infection. We pinpoint their antiviral mechanism of action to inhibition of the cellular serine protease, TMPRSS2, thereby preventing cleavage of the viral spike protein and TMPRSS2-mediated entry into the target cells. By means of single-cell RNA sequencing, we further locate the expression of the individual serpins to basal, ciliated, club, and goblet cells. Our results add to the importance of genetic variations as determinants for SARS-CoV-2 susceptibility and suggest that genetic deficiencies of cellular serpins might represent risk factors for severe COVID-19. Our study further highlights TMPRSS2 as a promising target for antiviral intervention and opens the door for the usage of locally administered serpins as a treatment against COVID-19. IMPORTANCE Identification of host factors affecting individual SARS-CoV-2 susceptibility will provide a better understanding of the large variations in disease severity and will identify potential factors that can be used, or targeted, in antiviral drug development. With the use of an advanced lung cell model established from several human donors, we identified cellular protease inhibitors, serpins, as host factors that restrict SARS-CoV-2 infection. The antiviral mechanism was found to be mediated by the inhibition of a serine protease, TMPRSS2, which results in a blockage of viral entry into target cells. Potential treatments with these serpins would not only reduce the overall viral burden in the patients, but also block the infection at an early time point, reducing the risk for the hyperactive immune response common in patients with severe COVID-19.

Mapping Rora expression in resting and activated CD4+ T cells.

The transcription factor Rora has been shown to be important for the development of ILC2 and the regulation of ILC3, macrophages and Treg cells. Here we investigate the role of Rora across CD4+ T cells in general, but with an emphasis on Th2 cells, both in vitro as well as in the context of several in vivo type 2 infection models. We dissect the function of Rora using overexpression and a CD4-conditional Rora-knockout mouse, as well as a RORA-reporter mouse. We establish the importance of Rora in CD4+ T cells for controlling lung inflammation induced by Nippostrongylus brasiliensis infection, and have measured the effect on downstream genes using RNA-seq. Using a systematic stimulation screen of CD4+ T cells, coupled with RNA-seq, we identify upstream regulators of Rora, most importantly IL-33 and CCL7. Our data suggest that Rora is a negative regulator of the immune system, possibly through several downstream pathways, and is under control of the local microenvironment.

Stress and odorant receptor feedback during a critical period after hatching regulates olfactory sensory neuron differentiation in Drosophila.

Here, we reveal that the regulation of Drosophila odorant receptor (OR) expression during the pupal stage is permissive and imprecise. We found that directly after hatching an OR feedback mechanism both directs and refines OR expression. We demonstrate that, as in mice, dLsd1 and Su(var)3-9 balance heterochromatin formation to direct OR expression. We show that the expressed OR induces dLsd1 and Su(var)3-9 expression, linking OR level and possibly function to OR expression. OR expression refinement shows a restricted duration, suggesting that a gene regulatory critical period brings olfactory sensory neuron differentiation to an end. Consistent with a change in differentiation, stress during the critical period represses dLsd1 and Su(var)3-9 expression and makes the early permissive OR expression permanent. This induced permissive gene regulatory state makes OR expression resilient to stress later in life. Hence, during a critical period OR feedback, similar to in mouse OR selection, defines adult OR expression in Drosophila.

Meta-Analysis of Gene Popularity: Less Than Half of Gene Citations Stem from Gene Regulatory Networks.

The reasons for selecting a gene for further study might vary from historical momentum to funding availability, thus leading to unequal attention distribution among all genes. However, certain biological features tend to be overlooked in evaluating a gene's popularity. Here we present a meta-analysis of the reasons why different genes have been studied and to what extent, with a focus on the gene-specific biological features. From unbiased datasets we can define biological properties of genes that reasonably may affect their perceived importance. We make use of both linear and nonlinear computational approaches for estimating gene popularity to then compare their relative importance. We find that roughly 25% of the studies are the result of a historical positive feedback, which we may think of as social reinforcement. Of the remaining features, gene family membership is the most indicative followed by disease relevance and finally regulatory pathway association. Disease relevance has been an important driver until the 1990s, after which the focus shifted to exploring every single gene. We also present a resource that allows one to study the impact of reinforcement, which may guide our research toward genes that have not yet received proportional attention.

Regulation of Decay Accelerating Factor Primes Human Germinal Center B Cells for Phagocytosis.

Germinal centers (GC) are sites for extensive B cell proliferation and homeostasis is maintained by programmed cell death. The complement regulatory protein Decay Accelerating Factor (DAF) blocks complement deposition on host cells and therefore also phagocytosis of cells. Here, we show that B cells downregulate DAF upon BCR engagement and that T cell-dependent stimuli preferentially led to activation of DAFlo B cells. Consistent with this, a majority of light and dark zone GC B cells were DAFlo and susceptible to complement-dependent phagocytosis, as compared with DAFhi GC B cells. We could also show that the DAFhi GC B cell subset had increased expression of the plasma cell marker Blimp-1. DAF expression was also modulated during B cell hematopoiesis in the human bone marrow. Collectively, our results reveal a novel role of DAF to pre-prime activated human B cells for phagocytosis prior to apoptosis.

Insertional mutagenesis in the zoonotic pathogen Chlamydia caviae.

The ability to introduce targeted genetic modifications in microbial genomes has revolutionized our ability to study the role and mode of action of individual bacterial virulence factors. Although the fastidious lifestyle of obligate intracellular bacterial pathogens poses a technical challenge to such manipulations, the last decade has produced significant advances in our ability to conduct molecular genetic analysis in Chlamydia trachomatis, a major bacterial agent of infertility and blindness. Similar approaches have not been established for the closely related veterinary Chlamydia spp., which cause significant economic damage, as well as rare but potentially life-threatening infections in humans. Here we demonstrate the feasibility of conducting site-specific mutagenesis for disrupting virulence genes in C. caviae, an agent of guinea pig inclusion conjunctivitis that was recently identified as a zoonotic agent in cases of severe community-acquired pneumonia. Using this approach, we generated C. caviae mutants deficient for the secreted effector proteins IncA and SinC. We demonstrate that C. caviae IncA plays a role in mediating fusion of the bacteria-containing vacuoles inhabited by C. caviae. Moreover, using a chicken embryo infection model, we provide first evidence for a role of SinC in C. caviae virulence in vivo.

CRISPR Screening in Single Cells.

The combination of single-cell RNA-seq and CRISPR allows for efficient interrogation of possibly any number of genes, only limited by the sequencing capability. Here we describe the current protocols for CRISPR screening in single cells, from cloning and virus production to generating sequencing data.

Genome-wide CRISPR Screens in T Helper Cells Reveal Pervasive Crosstalk between Activation and Differentiation.

T helper type 2 (Th2) cells are important regulators of mammalian adaptive immunity and have relevance for infection, autoimmunity, and tumor immunology. Using a newly developed, genome-wide retroviral CRISPR knockout (KO) library, combined with RNA-seq, ATAC-seq, and ChIP-seq, we have dissected the regulatory circuitry governing activation and differentiation of these cells. Our experiments distinguish cell activation versus differentiation in a quantitative framework. We demonstrate that these two processes are tightly coupled and are jointly controlled by many transcription factors, metabolic genes, and cytokine/receptor pairs. There are only a small number of genes regulating differentiation without any role in activation. By combining biochemical and genetic data, we provide an atlas for Th2 differentiation, validating known regulators and identifying factors, such as Pparg and Bhlhe40, as part of the core regulatory network governing Th2 helper cell fates.

Single-cell reconstruction of the early maternal-fetal interface in humans.

During early human pregnancy the uterine mucosa transforms into the decidua, into which the fetal placenta implants and where placental trophoblast cells intermingle and communicate with maternal cells. Trophoblast-decidual interactions underlie common diseases of pregnancy, including pre-eclampsia and stillbirth. Here we profile the transcriptomes of about 70,000 single cells from first-trimester placentas with matched maternal blood and decidual cells. The cellular composition of human decidua reveals subsets of perivascular and stromal cells that are located in distinct decidual layers. There are three major subsets of decidual natural killer cells that have distinctive immunomodulatory and chemokine profiles. We develop a repository of ligand-receptor complexes and a statistical tool to predict the cell-type specificity of cell-cell communication via these molecular interactions. Our data identify many regulatory interactions that prevent harmful innate or adaptive immune responses in this environment. Our single-cell atlas of the maternal-fetal interface reveals the cellular organization of the decidua and placenta, and the interactions that are critical for placentation and reproductive success.

Gene expression variability across cells and species shapes innate immunity.

As the first line of defence against pathogens, cells mount an innate immune response, which varies widely from cell to cell. The response must be potent but carefully controlled to avoid self-damage. How these constraints have shaped the evolution of innate immunity remains poorly understood. Here we characterize the innate immune response's transcriptional divergence between species and variability in expression among cells. Using bulk and single-cell transcriptomics in fibroblasts and mononuclear phagocytes from different species, challenged with immune stimuli, we map the architecture of the innate immune response. Transcriptionally diverging genes, including those that encode cytokines and chemokines, vary across cells and have distinct promoter structures. Conversely, genes that are involved in the regulation of this response, such as those that encode transcription factors and kinases, are conserved between species and display low cell-to-cell variability in expression. We suggest that this expression pattern, which is observed across species and conditions, has evolved as a mechanism for fine-tuned regulation to achieve an effective but balanced response.

Genome-wide analyses reveal the IRE1a-XBP1 pathway promotes T helper cell differentiation by resolving secretory stress and accelerating proliferation.

BACKGROUND: The IRE1a-XBP1 pathway is a conserved adaptive mediator of the unfolded protein response. The pathway is indispensable for the development of secretory cells by facilitating protein folding and enhancing secretory capacity. In the immune system, it is known to function in dendritic cells, plasma cells, and eosinophil development and differentiation, while its role in T helper cell is unexplored. Here, we investigated the role of the IRE1a-XBP1 pathway in regulating activation and differentiation of type-2 T helper cell (Th2), a major T helper cell type involved in allergy, asthma, helminth infection, pregnancy, and tumor immunosuppression. METHODS: We perturbed the IRE1a-XBP1 pathway and interrogated its role in Th2 cell differentiation. We performed genome-wide transcriptomic analysis of differential gene expression to reveal IRE1a-XBP1 pathway-regulated genes and predict their biological role. To identify direct target genes of XBP1 and define XBP1's regulatory network, we performed XBP1 ChIPmentation (ChIP-seq). We validated our predictions by flow cytometry, ELISA, and qPCR. We also used a fluorescent ubiquitin cell cycle indicator mouse to demonstrate the role of XBP1 in the cell cycle. RESULTS: We show that Th2 lymphocytes induce the IRE1a-XBP1 pathway during in vitro and in vivo activation. Genome-wide transcriptomic analysis of differential gene expression by perturbing the IRE1a-XBP1 pathway reveals XBP1-controlled genes and biological pathways. Performing XBP1 ChIPmentation (ChIP-seq) and integrating with transcriptomic data, we identify XBP1-controlled direct target genes and its transcriptional regulatory network. We observed that the IRE1a-XBP1 pathway controls cytokine secretion and the expression of two Th2 signature cytokines, IL13 and IL5. We also discovered that the IRE1a-XBP1 pathway facilitates activation-dependent Th2 cell proliferation by facilitating cell cycle progression through S and G2/M phase. CONCLUSIONS: We confirm and detail the critical role of the IRE1a-XBP1 pathway during Th2 lymphocyte activation in regulating cytokine expression, secretion, and cell proliferation. Our high-quality genome-wide XBP1 ChIP and gene expression data provide a rich resource for investigating XBP1-regulated genes. We provide a browsable online database available at http://data.teichlab.org .

Unbiased classification of mosquito blood cells by single-cell genomics and high-content imaging.

Mosquito blood cells are immune cells that help control infection by vector-borne pathogens. Despite their importance, little is known about mosquito blood cell biology beyond morphological and functional criteria used for their classification. Here, we combined the power of single-cell RNA sequencing, high-content imaging flow cytometry, and single-molecule RNA hybridization to analyze a subset of blood cells of the malaria mosquito Anopheles gambiae By demonstrating that blood cells express nearly half of the mosquito transcriptome, our dataset represents an unprecedented view into their transcriptional program. Analyses of differentially expressed genes identified transcriptional signatures of two cell types and provide insights into the current classification of these cells. We further demonstrate the active transfer of a cellular marker between blood cells that may confound their identification. We propose that cell-to-cell exchange may contribute to cellular diversity and functional plasticity seen across biological systems.

Mixed-species RNAseq analysis of human lymphoma cells adhering to mouse stromal cells identifies a core gene set that is also differentially expressed in the lymph node microenvironment of mantle cell lymphoma and chronic lymphocytic leukemia patients.

A subset of hematologic cancer patients is refractory to treatment or suffers relapse, due in part to minimal residual disease, whereby some cancer cells survive treatment. Cell-adhesion-mediated drug resistance is an important mechanism, whereby cancer cells receive survival signals via interaction with e.g. stromal cells. No genome-wide studies of in vitro systems have yet been performed to compare gene expression in different cell subsets within a co-culture and cells grown separately. Using RNA sequencing and species-specific read mapping, we compared transcript levels in human Jeko-1 mantle cell lymphoma cells stably adhered to mouse MS-5 stromal cells or in suspension within a co-culture or cultured separately as well as in stromal cells in co-culture or in separate culture. From 1050 differentially expressed transcripts in adherent mantle cell lymphoma cells, we identified 24 functional categories that together represent four main functional themes, anti-apoptosis, B-cell signaling, cell adhesion/migration and early mitosis. A comparison with previous mantle cell lymphoma and chronic lymphocytic leukemia studies, of gene expression differences between lymph node and blood, identified 116 genes that are differentially expressed in all three studies. From these genes, we suggest a core set of genes (CCL3, CCL4, DUSP4, ETV5, ICAM1, IL15RA, IL21R, IL4I1, MFSD2A, NFKB1, NFKBIE, SEMA7A, TMEM2) characteristic of cells undergoing cell-adhesion-mediated microenvironment signaling in mantle cell lymphoma/chronic lymphocytic leukemia. The model system developed and characterized here together with the core gene set will be useful for future studies of pathways that mediate increased cancer cell survival and drug resistance mechanisms.