Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis

A hybrid demultiplexing strategy that improves performance and robustness of cell hashing

Cell hashing, a nucleotide barcode-based method that allows users to pool multiple samples and demultiplex in downstream analysis, has gained widespread popularity in single-cell sequencing due to its compatibility, simplicity, and cost-effectiveness. Despite these advantages, the performance of this method remains unsatisfactory under certain circumstances, especially in experiments that have imbalanced sample sizes or use many hashtag antibodies. Here, we introduce a hybrid demultiplexing strategy that increases accuracy and cell recovery in multi-sample single-cell experiments. This approach correlates the results of cell hashing and genetic variant clustering, enabling precise and efficient cell identity determination without additional experimental costs or efforts. In addition, we developed HTOreader, a demultiplexing tool for cell hashing that improves the accuracy of cut-off calling by avoiding the dominance of negative signals in experiments with many hashtags or imbalanced sample sizes. When compared to existing methods using real-world datasets, this hybrid approach and HTOreader consistently generate reliable results with increased accuracy and cell recovery.

Meiotic transcriptional reprogramming mediated by cell-cell communications in humans and mice revealed by scATAC-seq and scRNA-seq

Meiosis is a highly complex process significantly influenced by transcriptional regulation. However, studies on the mechanisms that govern transcriptomic changes during meiosis, especially in prophase I, are limited. Here, we performed single-cell ATAC-seq of human testis tissues and observed reprogramming during the transition from zygotene to pachytene spermatocytes. This event, conserved in mice, involved the deactivation of genes associated with meiosis after reprogramming and the activation of those related to spermatogenesis before their functional onset. Furthermore, we identified 282 transcriptional regulators (TRs) that underwent activation or deactivation subsequent to this process. Evidence suggested that physical contact signals from Sertoli cells may regulate these TRs in spermatocytes, while secreted ENHO signals may alter metabolic patterns in these cells. Our results further indicated that defective transcriptional reprogramming may be associated with non-obstructive azoospermia (NOA). This study revealed the importance of both physical contact and secreted signals between Sertoli cells and germ cells in meiotic progression.

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

BACKGROUND

Normalization is a critical step in the analysis of single-cell RNA-sequencing (scRNA-seq) datasets. Its main goal is to make gene counts comparable within and between cells. To do so, normalization methods must account for technical and biological variability. Numerous normalization methods have been developed addressing different sources of dispersion and making specific assumptions about the count data.

MAIN BODY

The selection of a normalization method has a direct impact on downstream analysis, for example differential gene expression and cluster identification. Thus, the objective of this review is to guide the reader in making an informed decision on the most appropriate normalization method to use. To this aim, we first give an overview of the different single cell sequencing platforms and methods commonly used including isolation and library preparation protocols. Next, we discuss the inherent sources of variability of scRNA-seq datasets. We describe the categories of normalization methods and include examples of each. We also delineate imputation and batch-effect correction methods. Furthermore, we describe data-driven metrics commonly used to evaluate the performance of normalization methods. We also discuss common scRNA-seq methods and toolkits used for integrated data analysis.

CONCLUSIONS

According to the correction performed, normalization methods can be broadly classified as within and between-sample algorithms. Moreover, with respect to the mathematical model used, normalization methods can further be classified into: global scaling methods, generalized linear models, mixed methods, and machine learning-based methods. Each of these methods depict pros and cons and make different statistical assumptions. However, there is no better performing normalization method. Instead, metrics such as silhouette width, K-nearest neighbor batch-effect test, or Highly Variable Genes are recommended to assess the performance of normalization methods.

Reconstruction Set Test (RESET): A computationally efficient method for single sample gene set testing based on randomized reduced rank reconstruction error

We have developed a new, and analytically novel, single sample gene set testing method called Reconstruction Set Test (RESET). RESET quantifies gene set importance based on the ability of set genes to reconstruct values for all measured genes. RESET is realized using a computationally efficient randomized reduced rank reconstruction algorithm (available via the RESET R package on CRAN) that can effectively detect patterns of differential abundance and differential correlation for self-contained and competitive scenarios. As demonstrated using real and simulated scRNA-seq data, RESET provides superior performance at a lower computational cost relative to other single sample approaches.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

Single-cell analysis of immune and stroma cell remodeling in clear cell renal cell carcinoma primary tumors and bone metastatic lesions

BACKGROUND

Despite therapeutic advances, once a cancer has metastasized to the bone, it represents a highly morbid and lethal disease. One third of patients with advanced clear cell renal cell carcinoma (ccRCC) present with bone metastasis at the time of diagnosis. However, the bone metastatic niche in humans, including the immune and stromal microenvironments, has not been well-defined, hindering progress towards identification of therapeutic targets.

METHODS

We collected fresh patient samples and performed single-cell transcriptomic profiling of solid metastatic tissue (Bone Met), liquid bone marrow at the vertebral level of spinal cord compression (Involved), and liquid bone marrow from a different vertebral body distant from the tumor site but within the surgical field (Distal), as well as bone marrow from patients undergoing hip replacement surgery (Benign). In addition, we incorporated single-cell data from primary ccRCC tumors (ccRCC Primary) for comparative analysis.

RESULTS

The bone marrow of metastatic patients is immune-suppressive, featuring increased, exhausted CD8 + cytotoxic T cells, T regulatory cells, and tumor-associated macrophages (TAM) with distinct transcriptional states in metastatic lesions. Bone marrow stroma from tumor samples demonstrated a tumor-associated mesenchymal stromal cell population (TA-MSC) that appears to be supportive of epithelial-to mesenchymal transition (EMT), bone remodeling, and a cancer-associated fibroblast (CAFs) phenotype. This stromal subset is associated with poor progression-free and overall survival and also markedly upregulates bone remodeling through the dysregulation of RANK/RANKL/OPG signaling activity in bone cells, ultimately leading to bone resorption.

CONCLUSIONS

These results provide a comprehensive analysis of the bone marrow niche in the setting of human metastatic cancer and highlight potential therapeutic targets for both cell populations and communication channels.

Representation and quantification of module activity from omics data with rROMA

The efficiency of analyzing high-throughput data in systems biology has been demonstrated in numerous studies, where molecular data, such as transcriptomics and proteomics, offers great opportunities for understanding the complexity of biological processes. One important aspect of data analysis in systems biology is the shift from a reductionist approach that focuses on individual components to a more integrative perspective that considers the system as a whole, where the emphasis shifted from differential expression of individual genes to determining the activity of gene sets. Here, we present the rROMA software package for fast and accurate computation of the activity of gene sets with coordinated expression. The rROMA package incorporates significant improvements in the calculation algorithm, along with the implementation of several functions for statistical analysis and visualizing results. These additions greatly expand the package's capabilities and offer valuable tools for data analysis and interpretation. It is an open-source package available on github at: www.github.com/sysbio-curie/rROMA . Based on publicly available transcriptomic datasets, we applied rROMA to cystic fibrosis, highlighting biological mechanisms potentially involved in the establishment and progression of the disease and the associated genes. Results indicate that rROMA can detect disease-related active signaling pathways using transcriptomic and proteomic data. The results notably identified a significant mechanism relevant to cystic fibrosis, raised awareness of a possible bias related to cell culture, and uncovered an intriguing gene that warrants further investigation.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

scPML: pathway-based multi-view learning for cell type annotation from single-cell RNA-seq data

Recent developments in single-cell technology have enabled the exploration of cellular heterogeneity at an unprecedented level, providing invaluable insights into various fields, including medicine and disease research. Cell type annotation is an essential step in its omics research. The mainstream approach is to utilize well-annotated single-cell data to supervised learning for cell type annotation of new singlecell data. However, existing methods lack good generalization and robustness in cell annotation tasks, partially due to difficulties in dealing with technical differences between datasets, as well as not considering the heterogeneous associations of genes in regulatory mechanism levels. Here, we propose the scPML model, which utilizes various gene signaling pathway data to partition the genetic features of cells, thus characterizing different interaction maps between cells. Extensive experiments demonstrate that scPML performs better in cell type annotation and detection of unknown cell types from different species, platforms, and tissues.

A multifunctional Wnt regulator underlies the evolution of rodent stripe patterns

Animal pigment patterns are excellent models to elucidate mechanisms of biological organization. Although theoretical simulations, such as Turing reaction-diffusion systems, recapitulate many animal patterns, they are insufficient to account for those showing a high degree of spatial organization and reproducibility. Here, we study the coat of the African striped mouse (Rhabdomys pumilio) to uncover how periodic stripes form. Combining transcriptomics, mathematical modelling and mouse transgenics, we show that the Wnt modulator Sfrp2 regulates the distribution of hair follicles and establishes an embryonic prepattern that foreshadows pigment stripes. Moreover, by developing in vivo gene editing in striped mice, we find that Sfrp2 knockout is sufficient to alter the stripe pattern. Strikingly, mutants exhibited changes in pigmentation, revealing that Sfrp2 also regulates hair colour. Lastly, through evolutionary analyses, we find that striped mice have evolved lineage-specific changes in regulatory elements surrounding Sfrp2, many of which may be implicated in modulating the expression of this gene. Altogether, our results show that a single factor controls coat pattern formation by acting both as an orienting signalling mechanism and a modulator of pigmentation. More broadly, our work provides insights into how spatial patterns are established in developing embryos and the mechanisms by which phenotypic novelty originates.

A hybrid demultiplexing strategy that improves performance and robustness of cell hashing

Recent advances in single cell RNA sequencing allow users to pool multiple samples and demultiplex in downstream analysis, which greatly increase experimental efficiency and cost-effectiveness. Among all the demultiplexing methods, nucleotide barcode-based cell hashing has gained widespread popularity due to its compatibility and simplicity. Despite these advantages, certain issues of this technic remain to be solved, such as challenges in distinguishing true positive from background, high reagent cost for samples with large cell numbers, and unpredictable false negative and false doublet rates. Here, we propose a hybrid demultiplexing strategy that increases calling accuracy and cell recovery of cell hashing without adding experimental cost. In this approach, we computationally cluster all single cells based on their natural genetic variations and assign donor identity by finding the dominant hashtag in each genotype cluster. This hybrid strategy assigns donor identity to any cell that is identified as singlet by either genotype clustering or cell hashing, which allows us to demultiplex most majority of cells even if only a small fraction of cells are labeled with hashtags. When comparing its performance with cell hashing on multiple real-world datasets, this hybrid approach consistently generates reliable demultiplexing results with increased cell recovery and accuracy.

CellGO: a novel deep learning-based framework and webserver for cell-type-specific gene function interpretation

Interpreting the function of genes and gene sets identified from omics experiments remains a challenge, as current pathway analysis tools often fail to consider the critical biological context, such as tissue or cell-type specificity. To address this limitation, we introduced CellGO. CellGO tackles this challenge by leveraging the visible neural network (VNN) and single-cell gene expressions to mimic cell-type-specific signaling propagation along the Gene Ontology tree within a cell. This design enables a novel scoring system to calculate the cell-type-specific gene-pathway paired active scores, based on which, CellGO is able to identify cell-type-specific active pathways associated with single genes. In addition, by aggregating the activities of single genes, CellGO extends its capability to identify cell-type-specific active pathways for a given gene set. To enhance biological interpretation, CellGO offers additional features, including the identification of significantly active cell types and driver genes and community analysis of pathways. To validate its performance, CellGO was assessed using a gene set comprising mixed cell-type markers, confirming its ability to discern active pathways across distinct cell types. Subsequent benchmarking analyses demonstrated CellGO's superiority in effectively identifying cell types and their corresponding cell-type-specific pathways affected by gene knockouts, using either single genes or sets of genes differentially expressed between knockout and control samples. Moreover, CellGO demonstrated its ability to infer cell-type-specific pathogenesis for disease risk genes. Accessible as a Python package, CellGO also provides a user-friendly web interface, making it a versatile and accessible tool for researchers in the field.

Computational single cell oncology: state of the art

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

ZNF683 marks a CD8+ T cell population associated with anti-tumor immunity following anti-PD-1 therapy for Richter syndrome

Unlike many other hematologic malignancies, Richter syndrome (RS), an aggressive B cell lymphoma originating from indolent chronic lymphocytic leukemia, is responsive to PD-1 blockade. To discover the determinants of response, we analyze single-cell transcriptome data generated from 17 bone marrow samples longitudinally collected from 6 patients with RS. Response is associated with intermediate exhausted CD8 effector/effector memory T cells marked by high expression of the transcription factor ZNF683, determined to be evolving from stem-like memory cells and divergent from terminally exhausted cells. This signature overlaps with that of tumor-infiltrating populations from anti-PD-1 responsive solid tumors. ZNF683 is found to directly target key T cell genes (TCF7, LMO2, CD69) and impact pathways of T cell cytotoxicity and activation. Analysis of pre-treatment peripheral blood from 10 independent patients with RS treated with anti-PD-1, as well as patients with solid tumors treated with anti-PD-1, supports an association of ZNF683high T cells with response.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

PanomiR: a systems biology framework for analysis of multi-pathway targeting by miRNAs

Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA-miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

C-type lectin receptor expression is a hallmark of neutrophils infiltrating the skin in epidermolysis bullosa acquisita

Introduction

Inflammatory epidermolysis bullosa acquisita (EBA) is characterized by a neutrophilic response to anti-type VII collagen (COL7) antibodies resulting in the development of skin inflammation and blistering. The antibody transfer model of EBA closely mirrors this EBA phenotype.

Methods

To better understand the changes induced in neutrophils upon recruitment from peripheral blood into lesional skin in EBA, we performed single-cell RNA-sequencing of whole blood and skin dissociate to capture minimally perturbed neutrophils and characterize their transcriptome.

Results

Through this approach, we identified clear distinctions between circulating activated neutrophils and intradermal neutrophils. Most strikingly, the gene expression of multiple C-type lectin receptors, which have previously been reported to orchestrate host defense against fungi and select bacteria, were markedly dysregulated. After confirming the upregulation of Clec4n, Clec4d, and Clec4e in experimental EBA as well as in lesional skin from patients with inflammatory EBA, we performed functional studies in globally deficient Clec4e-/- and Clec4d-/- mice as well as in neutrophil-specific Clec4n-/- mice. Deficiency in these genes did not reduce disease in the EBA model.

Discussion

Collectively, our results suggest that while the upregulation of Clec4n, Clec4d, and Clec4e is a hallmark of activated dermal neutrophil populations, their individual contribution to the pathogenesis of EBA is dispensable.

Target Genes of c-MYC and MYCN with Prognostic Power in Neuroblastoma Exhibit Different Expressions during Sympathoadrenal Development

Deregulation of the MYC family of transcription factors c-MYC (encoded by MYC), MYCN, and MYCL is prevalent in most human cancers, with an impact on tumor initiation and progression, as well as response to therapy. In neuroblastoma (NB), amplification of the MYCN oncogene and over-expression of MYC characterize approximately 40% and 10% of all high-risk NB cases, respectively. However, the mechanism and stage of neural crest development in which MYCN and c-MYC contribute to the onset and/or progression of NB are not yet fully understood. Here, we hypothesized that subtle differences in the expression of MYCN and/or c-MYC targets could more accurately stratify NB patients in different risk groups rather than using the expression of either MYC gene alone. We employed an integrative approach using the transcriptome of 498 NB patients from the SEQC cohort and previously defined c-MYC and MYCN target genes to model a multigene transcriptional risk score. Our findings demonstrate that defined sets of c-MYC and MYCN targets with significant prognostic value, effectively stratify NB patients into different groups with varying overall survival probabilities. In particular, patients exhibiting a high-risk signature score present unfavorable clinical parameters, including increased clinical risk, higher INSS stage, MYCN amplification, and disease progression. Notably, target genes with prognostic value differ between c-MYC and MYCN, exhibiting distinct expression patterns in the developing sympathoadrenal system. Genes associated with poor outcomes are mainly found in sympathoblasts rather than in chromaffin cells during the sympathoadrenal development.

Polygenic regression uncovers trait-relevant cellular contexts through pathway activation transformation of single-cell RNA sequencing data

Advances in single-cell RNA sequencing (scRNA-seq) techniques have accelerated functional interpretation of disease-associated variants discovered from genome-wide association studies (GWASs). However, identification of trait-relevant cell populations is often impeded by inherent technical noise and high sparsity in scRNA-seq data. Here, we developed scPagwas, a computational approach that uncovers trait-relevant cellular context by integrating pathway activation transformation of scRNA-seq data and GWAS summary statistics. scPagwas effectively prioritizes trait-relevant genes, which facilitates identification of trait-relevant cell types/populations with high accuracy in extensive simulated and real datasets. Cellular-level association results identified a novel subpopulation of naive CD8+ T cells related to COVID-19 severity and oligodendrocyte progenitor cell and microglia subsets with critical pathways by which genetic variants influence Alzheimer's disease. Overall, our approach provides new insights for the discovery of trait-relevant cell types and improves the mechanistic understanding of disease variants from a pathway perspective.

Heterogeneity in the gene regulatory landscape of leiomyosarcoma

Characterizing inter-tumor heterogeneity is crucial for selecting suitable cancer therapy, as the presence of diverse molecular subgroups of patients can be associated with disease outcome or response to treatment. While cancer subtypes are often characterized by differences in gene expression, the mechanisms driving these differences are generally unknown. We set out to model the regulatory mechanisms driving sarcoma heterogeneity based on patient-specific, genome-wide gene regulatory networks. We developed a new computational framework, PORCUPINE, which combines knowledge on biological pathways with permutation-based network analysis to identify pathways that exhibit significant regulatory heterogeneity across a patient population. We applied PORCUPINE to patient-specific leiomyosarcoma networks modeled on data from The Cancer Genome Atlas and validated our results in an independent dataset from the German Cancer Research Center. PORCUPINE identified 37 heterogeneously regulated pathways, including pathways representing potential targets for treatment of subgroups of leiomyosarcoma patients, such as FGFR and CTLA4 inhibitory signaling. We validated the detected regulatory heterogeneity through analysis of networks and chromatin states in leiomyosarcoma cell lines. We showed that the heterogeneity identified with PORCUPINE is not associated with methylation profiles or clinical features, thereby suggesting an independent mechanism of patient heterogeneity driven by the complex landscape of gene regulatory interactions.

Endothelial sensing of AHR ligands regulates intestinal homeostasis

Endothelial cells line the blood and lymphatic vasculature, and act as an essential physical barrier, control nutrient transport, facilitate tissue immunosurveillance and coordinate angiogenesis and lymphangiogenesis1,2. In the intestine, dietary and microbial cues are particularly important in the regulation of organ homeostasis. However, whether enteric endothelial cells actively sense and integrate such signals is currently unknown. Here we show that the aryl hydrocarbon receptor (AHR) acts as a critical node for endothelial cell sensing of dietary metabolites in adult mice and human primary endothelial cells. We first established a comprehensive single-cell endothelial atlas of the mouse small intestine, uncovering the cellular complexity and functional heterogeneity of blood and lymphatic endothelial cells. Analyses of AHR-mediated responses at single-cell resolution identified tissue-protective transcriptional signatures and regulatory networks promoting cellular quiescence and vascular normalcy at steady state. Endothelial AHR deficiency in adult mice resulted in dysregulated inflammatory responses and the initiation of proliferative pathways. Furthermore, endothelial sensing of dietary AHR ligands was required for optimal protection against enteric infection. In human endothelial cells, AHR signalling promoted quiescence and restrained activation by inflammatory mediators. Together, our data provide a comprehensive dissection of the effect of environmental sensing across the spectrum of enteric endothelia, demonstrating that endothelial AHR signalling integrates dietary cues to maintain tissue homeostasis by promoting endothelial cell quiescence and vascular normalcy.

Single-Cell Multiomics

Single-cell RNA sequencing methods have led to improved understanding of the heterogeneity and transcriptomic states present in complex biological systems. Recently, the development of novel single-cell technologies for assaying additional modalities, specifically genomic, epigenomic, proteomic, and spatial data, allows for unprecedented insight into cellular biology. While certain technologies collect multiple measurements from the same cells simultaneously, even when modalities are separately assayed in different cells, we can apply novel computational methods to integrate these data. The application of computational integration methods to multimodal paired and unpaired data results in rich information about the identities of the cells present and the interactions between different levels of biology, such as between genetic variation and transcription. In this review, we both discuss the single-cell technologies for measuring these modalities and describe and characterize a variety of computational integration methods for combining the resulting data to leverage multimodal information toward greater biological insight.

Regulation of Photosensitivity by the Hippo Pathway in Lupus Skin

OBJECTIVE

Photosensitivity is one of the most common manifestations of systemic lupus erythematosus (SLE), yet its pathogenesis is not well understood. The normal-appearing epidermis of patients with SLE exhibits increased ultraviolet B (UVB)-driven cell death that persists in cell culture. Here, we investigated the role of epigenetic modification and Hippo signaling in enhanced UVB-induced apoptosis seen in SLE keratinocytes.

METHODS

We analyzed DNA methylation in cultured keratinocytes from SLE patients compared to keratinocytes from healthy controls (n = 6/group). Protein expression was validated in cultured keratinocytes using immunoblotting and immunofluorescence. An immortalized keratinocyte line overexpressing WWC1 was generated via lentiviral vector. WWC1-driven changes were inhibited using a large tumor suppressor kinase 1/2 (LATS1/2) inhibitor (TRULI) and small interfering RNA (siRNA). The interaction between the Yes-associated protein (YAP) and the transcriptional enhancer associate domain (TEAD) was inhibited by overexpression of an N/TERT cell line expressing a tetracycline-inducible green fluorescent protein-tagged protein that inhibits YAP-TEAD binding (TEADi). Apoptosis was assessed using cleaved caspase 3/7 and TUNEL staining.

RESULTS

Hippo signaling was the top differentially methylated pathway in SLE versus control keratinocytes. SLE keratinocytes (n = 6) showed significant hypomethylation (Δβ = -0.153) and thus overexpression of the Hippo regulator WWC1 (P = 0.002). WWC1 overexpression increased LATS1/2 kinase activation, leading to YAP cytoplasmic retention and altered proapoptotic transcription in SLE keratinocytes. Accordingly, UVB-mediated apoptosis in keratinocytes could be enhanced by WWC1 overexpression or YAP-TEAD inhibition, mimicking SLE keratinocytes. Importantly, inhibition of LATS1/2 with either the chemical inhibitor TRULI or siRNA effectively eliminated enhanced UVB-apoptosis in SLE keratinocytes.

CONCLUSION

Our work unravels a novel driver of photosensitivity in SLE: overactive Hippo signaling in SLE keratinocytes restricts YAP transcriptional activity, leading to shifts that promote UVB apoptosis.

IndepthPathway: an integrated tool for in-depth pathway enrichment analysis based on single-cell sequencing data

MOTIVATION

Single-cell sequencing enables exploring the pathways and processes of cells, and cell populations. However, there is a paucity of pathway enrichment methods designed to tolerate the high noise and low gene coverage of this technology. When gene expression data are noisy and signals are sparse, testing pathway enrichment based on the genes expression may not yield statistically significant results, which is particularly problematic when detecting the pathways enriched in less abundant cells that are vulnerable to disturbances.

RESULTS

In this project, we developed a Weighted Concept Signature Enrichment Analysis specialized for pathway enrichment analysis from single-cell transcriptomics (scRNA-seq). Weighted Concept Signature Enrichment Analysis took a broader approach for assessing the functional relations of pathway gene sets to differentially expressed genes, and leverage the cumulative signature of molecular concepts characteristic of the highly differentially expressed genes, which we termed as the universal concept signature, to tolerate the high noise and low coverage of this technology. We then incorporated Weighted Concept Signature Enrichment Analysis into an R package called "IndepthPathway" for biologists to broadly leverage this method for pathway analysis based on bulk and single-cell sequencing data. Through simulating technical variability and dropouts in gene expression characteristic of scRNA-seq as well as benchmarking on a real dataset of matched single-cell and bulk RNAseq data, we demonstrate that IndepthPathway presents outstanding stability and depth in pathway enrichment results under stochasticity of the data, thus will substantially improve the scientific rigor of the pathway analysis for single-cell sequencing data.

AVAILABILITY AND IMPLEMENTATION

The IndepthPathway R package is available through: https://github.com/wangxlab/IndepthPathway.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Applications of single-cell RNA sequencing in drug discovery and development

Single-cell technologies, particularly single-cell RNA sequencing (scRNA-seq) methods, together with associated computational tools and the growing availability of public data resources, are transforming drug discovery and development. New opportunities are emerging in target identification owing to improved disease understanding through cell subtyping, and highly multiplexed functional genomics screens incorporating scRNA-seq are enhancing target credentialling and prioritization. ScRNA-seq is also aiding the selection of relevant preclinical disease models and providing new insights into drug mechanisms of action. In clinical development, scRNA-seq can inform decision-making via improved biomarker identification for patient stratification and more precise monitoring of drug response and disease progression. Here, we illustrate how scRNA-seq methods are being applied in key steps in drug discovery and development, and discuss ongoing challenges for their implementation in the pharmaceutical industry.

Directionality of developing skeletal muscles is set by mechanical forces

Formation of oriented myofibrils is a key event in musculoskeletal development. However, the mechanisms that drive myocyte orientation and fusion to control muscle directionality in adults remain enigmatic. Here, we demonstrate that the developing skeleton instructs the directional outgrowth of skeletal muscle and other soft tissues during limb and facial morphogenesis in zebrafish and mouse. Time-lapse live imaging reveals that during early craniofacial development, myoblasts condense into round clusters corresponding to future muscle groups. These clusters undergo oriented stretch and alignment during embryonic growth. Genetic perturbation of cartilage patterning or size disrupts the directionality and number of myofibrils in vivo. Laser ablation of musculoskeletal attachment points reveals tension imposed by cartilage expansion on the forming myofibers. Application of continuous tension using artificial attachment points, or stretchable membrane substrates, is sufficient to drive polarization of myocyte populations in vitro. Overall, this work outlines a biomechanical guidance mechanism that is potentially useful for engineering functional skeletal muscle.

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets

Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

Reconstruction Set Test (RESET): a computationally efficient method for single sample gene set testing based on randomized reduced rank reconstruction error

We have developed a new, and analytically novel, single sample gene set testing method called Reconstruction Set Test (RESET). RESET quantifies gene set importance at both the sample-level and for the entire dataset based on the ability of set genes to reconstruct values for all measured genes. RESET addresses four important limitations of current techniques: 1) existing single sample methods are designed to detect mean differences and struggle to identify differential correlation patterns, 2) computationally efficient techniques are self-contained methods and cannot directly detect competitive scenarios where set genes differ from non-set genes in the same sample, 3) the scores generated by current methods can only be accurately compared across samples for a single set and not between sets, and 4) the computational performance of even the fastest existing methods be significant on very large datasets. RESET is realized using a computationally efficient randomized reduced rank reconstruction algorithm (available via the RESET R package on CRAN) that can effectively detect patterns of differential abundance and differential correlation for self-contained and competitive scenarios. As demonstrated using real and simulated scRNA-seq data, RESET provides superior accuracy at a lower computational cost relative to other single sample approaches.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 hours post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Multi-modal quantification of pathway activity with MAYA

Signaling pathways can be activated through various cascades of genes depending on cell identity and biological context. Single-cell atlases now provide the opportunity to inspect such complexity in health and disease. Yet, existing reference tools for pathway scoring resume activity of each pathway to one unique common metric across cell types. Here, we present MAYA, a computational method that enables the automatic detection and scoring of the diverse modes of activation of biological pathways across cell populations. MAYA improves the granularity of pathway analysis by detecting subgroups of genes within reference pathways, each characteristic of a cell population and how it activates a pathway. Using multiple single-cell datasets, we demonstrate the biological relevance of identified modes of activation, the robustness of MAYA to noisy pathway lists and batch effect. MAYA can also predict cell types starting from lists of reference markers in a cluster-free manner. Finally, we show that MAYA reveals common modes of pathway activation in tumor cells across patients, opening the perspective to discover shared therapeutic vulnerabilities.

Haplotype-aware analysis of somatic copy number variations from single-cell transcriptomes

Genome instability and aberrant alterations of transcriptional programs both play important roles in cancer. Single-cell RNA sequencing (scRNA-seq) has the potential to investigate both genetic and nongenetic sources of tumor heterogeneity in a single assay. Here we present a computational method, Numbat, that integrates haplotype information obtained from population-based phasing with allele and expression signals to enhance detection of copy number variations from scRNA-seq. Numbat exploits the evolutionary relationships between subclones to iteratively infer single-cell copy number profiles and tumor clonal phylogeny. Analysis of 22 tumor samples, including multiple myeloma, gastric, breast and thyroid cancers, shows that Numbat can reconstruct the tumor copy number profile and precisely identify malignant cells in the tumor microenvironment. We identify genetic subpopulations with transcriptional signatures relevant to tumor progression and therapy resistance. Numbat requires neither sample-matched DNA data nor a priori genotyping, and is applicable to a wide range of experimental settings and cancer types.

Transcriptomic profiling of the developing brain revealed cell-type and brain-region specificity in a mouse model of prenatal stress

BACKGROUND

Prenatal stress (PS) is considered as a risk factor for many mental disorders. PS-induced transcriptomic alterations may contribute to the functional dysregulation during brain development. Here, we used RNA-seq to explore changes of gene expression in the mouse fetal brain after prenatal exposure to chronic unpredictable mild stress (CUMS).

RESULTS

We compared the stressed brains to the controls and identified groups of significantly differentially expressed genes (DEGs). GO analysis on up-regulated DEGs revealed enrichment for the cell cycle pathways, while down-regulated DEGs were mostly enriched in the neuronal pathways related to synaptic transmission. We further performed cell-type enrichment analysis using published scRNA-seq data from the fetal mouse brain and revealed cell-type-specificity for up- and down-regulated DEGs, respectively. The up-regulated DEGs were highly enriched in the radial glia, while down-regulated DEGs were enriched in different types of neurons. Cell deconvolution analysis further showed altered cell fractions in the stressed brain, indicating accumulation of neuroblast and impaired neurogenesis. Moreover, we also observed distinct brain-region expression pattern when mapping DEGs onto the developing Allen brain atlas. The up-regulated DEGs were primarily enriched in the dorsal forebrain regions including the cortical plate and hippocampal formation. Surprisingly, down-regulated DEGs were found excluded from the cortical region, but highly expressed on various regions in the ventral forebrain, midbrain and hindbrain.

CONCLUSION

Taken together, we provided an unbiased data source for transcriptomic alterations of the whole fetal brain after chronic PS, and reported differential cell-type and brain-region vulnerability of the developing brain in response to environmental insults during the pregnancy.

Assessment of Diagnosis, Prognosis and Immune Infiltration Response to the Expression of the Ferroptosis-Related Molecule HAMP in Clear Cell Renal Cell Carcinoma

BACKGROUND

Hepcidin antimicrobial peptide (HAMP) is a key factor in maintaining iron metabolism, which may induce ferroptosis when upregulated. However, its prognostic value and relation to immune infiltrating cells remains unclear.

METHODS

This study analyzed the expression levels of HAMP in the Oncomine, Timer and Ualcan databases, and examined its prognostic potential in KIRC with R programming. The Timer and GEPIA databases were used to estimate the correlations between HAMP and immune infiltration and the markers of immune cells. The intersection genes and the co-expression PPI network were constructed via STRING, R programming and GeneMANIA, and the hub genes were selected with Cytoscape. In addition, we analyzed the gene set enrichment and GO/KEGG pathways by GSEA.

RESULTS

Our study revealed higher HAMP expression levels in tumor tissues including KIRC, which were related to poor prognosis in terms of OS, DSS and PFI. The expression of HAMP was positively related to the immune infiltration level of macrophages, Tregs, etc., corresponding with the immune biomarkers. Based on the intersection genes, we constructed the PPI network and used the 10 top hub genes. Further, we performed a pathway enrichment analysis of the gene sets, including Huntington's disease, the JAK-STAT signaling pathway, ammonium ion metabolic process, and so on.

CONCLUSION

In summary, our study gave an insight into the potential prognosis of HAMP, which may act as a diagnostic biomarker and therapeutic target related to immune infiltration in KIRC.

Standardization of Single-Cell RNA-Sequencing Analysis Workflow to Study Drosophila Ovary

Developments in single-cell technology have considerably changed the way we study biology. Significant efforts have been made over the last few years to build comprehensive cell-type-specific transcriptomic atlases for a wide range of tissues in several model organisms in order to discover cell-type-specific markers and drivers of gene expression. One such tissue is the ovary of the fruit-fly Drosophila melanogaster, which is a popular model system with wide-ranging applications in the study of both development and disease. Three independent studies have recently produced comprehensive maps of cell-type-specific gene expression that describe both spatiotemporal regulation of the process of oogenesis and unique transcriptomic profiles of different cell types that constitute the ovary. In this chapter, we outlined the wet-lab protocol that was followed in our recent study for sample preparation and reanalyze the resultant dataset to discuss the benchmarks in data analysis, which are fundamental to comprehensive curation of the single-cell dataset representing the fly ovary.

Improving the performance of single-cell RNA-seq data mining based on relative expression orderings

The advent of single-cell RNA-sequencing (scRNA-seq) provides an unprecedented opportunity to explore gene expression profiles at the single-cell level. However, gene expression values vary over time and under different conditions even within the same cell. There is an urgent need for more stable and reliable feature variables at the single-cell level to depict cell heterogeneity. Thus, we construct a new feature matrix called the delta rank matrix (DRM) from scRNA-seq data by integrating an a priori gene interaction network, which transforms the unreliable gene expression value into a stable gene interaction/edge value on a single-cell basis. This is the first time that a gene-level feature has been transformed into an interaction/edge-level for scRNA-seq data analysis based on relative expression orderings. Experiments on various scRNA-seq datasets have demonstrated that DRM performs better than the original gene expression matrix in cell clustering, cell identification and pseudo-trajectory reconstruction. More importantly, the DRM really achieves the fusion of gene expressions and gene interactions and provides a method of measuring gene interactions at the single-cell level. Thus, the DRM can be used to find changes in gene interactions among different cell types, which may open up a new way to analyze scRNA-seq data from an interaction perspective. In addition, DRM provides a new method to construct a cell-specific network for each single cell instead of a group of cells as in traditional network construction methods. DRM's exceptional performance is due to its extraction of rich gene-association information on biological systems and stable characterization of cells.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

The landscape of tumor cell states and spatial organization in H3-K27M mutant diffuse midline glioma across age and location

Histone 3 lysine27-to-methionine (H3-K27M) mutations most frequently occur in diffuse midline gliomas (DMGs) of the childhood pons but are also increasingly recognized in adults. Their potential heterogeneity at different ages and midline locations is vastly understudied. Here, through dissecting the single-cell transcriptomic, epigenomic and spatial architectures of a comprehensive cohort of patient H3-K27M DMGs, we delineate how age and anatomical location shape glioma cell-intrinsic and -extrinsic features in light of the shared driver mutation. We show that stem-like oligodendroglial precursor-like cells, present across all clinico-anatomical groups, display varying levels of maturation dependent on location. We reveal a previously underappreciated relationship between mesenchymal cancer cell states and age, linked to age-dependent differences in the immune microenvironment. Further, we resolve the spatial organization of H3-K27M DMG cell populations and identify a mitotic oligodendroglial-lineage niche. Collectively, our study provides a powerful framework for rational modeling and therapeutic interventions.

Molecular Computational Anatomy: Unifying the Particle to Tissue Continuum via Measure Representations of the Brain

OBJECTIVE

The objective of this research is to unify the molecular representations of spatial transcriptomics and cellular scale histology with the tissue scales of computational anatomy for brain mapping.

IMPACT STATEMENT

We present a unified representation theory for brain mapping based on geometric varifold measures of the microscale deterministic structure and function with the statistical ensembles of the spatially aggregated tissue scales.

INTRODUCTION

Mapping across coordinate systems in computational anatomy allows us to understand structural and functional properties of the brain at the millimeter scale. New measurement technologies in digital pathology and spatial transcriptomics allow us to measure the brain molecule by molecule and cell by cell based on protein and transcriptomic functional identity. We currently have no mathematical representations for integrating consistently the tissue limits with the molecular particle descriptions. The formalism derived here demonstrates the methodology for transitioning consistently from the molecular scale of quantized particles-using mathematical structures as first introduced by Dirac as the class of generalized functions-to the tissue scales with methods originally introduced by Euler for fluids.

METHODS

We introduce two mathematical methods based on notions of generalized functions and statistical mechanics. We use geometric varifolds, a product measure on space and function, to represent functional states at the micro-scales-electrophysiology, molecular histology-integrated with a Boltzmann-like program to pass from deterministic particle descriptions to empirical probabilities on the functional states at the tissue scales.

RESULTS

Our space-function varifold representation provides a recipe for traversing from molecular to tissue scales in terms of a cascade of linear space scaling composed with nonlinear functional feature mapping. Following the cascade implies every scale is a geometric measure so that a universal family of measure norms can be introduced which quantifies the geodesic connection between brains in the orbit independent of the probing technology, whether it be RNA identities, Tau or amyloid histology, spike trains, or dense MR imagery.

CONCLUSIONS

We demonstrate a unified brain mapping theory for molecular and tissue scales based on geometric measure representations. We call the consistent aggregation of tissue scales from particle and cellular scales, molecular computational anatomy.

Polygenic enrichment distinguishes disease associations of individual cells in single-cell RNA-seq data

Single-cell RNA sequencing (scRNA-seq) provides unique insights into the pathology and cellular origin of disease. We introduce single-cell disease relevance score (scDRS), an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). We applied scDRS to 74 diseases/traits and 1.3 million single-cell gene-expression profiles across 31 tissues/organs. Cell-type-level results broadly recapitulated known cell-type-disease associations. Individual-cell-level results identified subpopulations of disease-associated cells not captured by existing cell-type labels, including T cell subpopulations associated with inflammatory bowel disease, partially characterized by their effector-like states; neuron subpopulations associated with schizophrenia, partially characterized by their spatial locations; and hepatocyte subpopulations associated with triglyceride levels, partially characterized by their higher ploidy levels. Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

Machine learning for cell type classification from single nucleus RNA sequencing data

With the advent of single cell/nucleus RNA sequencing (sc/snRNA-seq), the field of cell phenotyping is now a data-driven exercise providing statistical evidence to support cell type/state categorization. However, the task of classifying cells into specific, well-defined categories with the empirical data provided by sc/snRNA-seq remains nontrivial due to the difficulty in determining specific differences between related cell types with close transcriptional similarities, resulting in challenges with matching cell types identified in separate experiments. To investigate possible approaches to overcome these obstacles, we explored the use of supervised machine learning methods-logistic regression, support vector machines, random forests, neural networks, and light gradient boosting machine (LightGBM)-as approaches to classify cell types using snRNA-seq datasets from human brain middle temporal gyrus (MTG) and human kidney. Classification accuracy was evaluated using an F-beta score weighted in favor of precision to account for technical artifacts of gene expression dropout. We examined the impact of hyperparameter optimization and feature selection methods on F-beta score performance. We found that the best performing model for granular cell type classification in both datasets is a multinomial logistic regression classifier and that an effective feature selection step was the most influential factor in optimizing the performance of the machine learning pipelines.

A transcriptional cross species map of pancreatic islet cells

OBJECTIVE

Pancreatic islets of Langerhans secrete hormones to regulate systemic glucose levels. Emerging evidence suggests that islet cells are functionally heterogeneous to allow a fine-tuned and efficient endocrine response to physiological changes. A precise description of the molecular basis of this heterogeneity, in particular linking animal models to human islets, is an important step towards identifying the factors critical for endocrine cell function in physiological and pathophysiological conditions.

METHODS

In this study, we used single-cell RNA sequencing to profile more than 50'000 endocrine cells isolated from healthy human, pig and mouse pancreatic islets and characterize transcriptional heterogeneity and evolutionary conservation of those cells across the three species. We systematically delineated endocrine cell types and α- and β-cell heterogeneity through prior knowledge- and data-driven gene sets shared across species, which altogether capture common and differential cellular properties, transcriptional dynamics and putative driving factors of state transitions.

RESULTS

We showed that global endocrine expression profiles correlate, and that critical identity and functional markers are shared between species, while only approximately 20% of cell type enriched expression is conserved. We resolved distinct human α- and β-cell states that form continuous transcriptional landscapes. These states differentially activate maturation and hormone secretion programs, which are related to regulatory hormone receptor expression, signaling pathways and different types of cellular stress responses. Finally, we mapped mouse and pig cells to the human reference and observed that the spectrum of human α- and β-cell heterogeneity and aspects of such functional gene expression are better recapitulated in the pig than mouse data.

CONCLUSIONS

Here, we provide a high-resolution transcriptional map of healthy human islet cells and their murine and porcine counterparts, which is easily queryable via an online interface. This comprehensive resource informs future efforts that focus on pancreatic endocrine function, failure and regeneration, and enables to assess molecular conservation in islet biology across species for translational purposes.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Schwann cell precursors represent a neural crest-like state with biased multipotency

Schwann cell precursors (SCPs) are nerve-associated progenitors that can generate myelinating and non-myelinating Schwann cells but also are multipotent like the neural crest cells from which they originate. SCPs are omnipresent along outgrowing peripheral nerves throughout the body of vertebrate embryos. By using single-cell transcriptomics to generate a gene expression atlas of the entire neural crest lineage, we show that early SCPs and late migratory crest cells have similar transcriptional profiles characterised by a multipotent "hub" state containing cells biased towards traditional neural crest fates. SCPs keep diverging from the neural crest after being primed towards terminal Schwann cells and other fates, with different subtypes residing in distinct anatomical locations. Functional experiments using CRISPR-Cas9 loss-of-function further show that knockout of the common "hub" gene Sox8 causes defects in neural crest-derived cells along peripheral nerves by facilitating differentiation of SCPs towards sympathoadrenal fates. Finally, specific tumour populations found in melanoma, neurofibroma and neuroblastoma map to different stages of SCP/Schwann cell development. Overall, SCPs resemble migrating neural crest cells that maintain multipotency and become transcriptionally primed towards distinct lineages.

Regulatory Programs of B-cell Activation and Germinal Center Reaction Allow B-ALL Escape from CD19 CAR T-cell Therapy

Chimeric antigen receptor (CAR) T-cell therapy has led to tremendous successes in the treatment of B-cell malignancies. However, a large fraction of treated patients relapse, often with disease expressing reduced levels of the target antigen. Here, we report that exposing CD19+ B-cell acute lymphoblastic leukemia (B-ALL) cells to CD19 CAR T cells reduced CD19 expression within hours. Initially, CD19 CAR T cells caused clustering of CD19 at the T cell-leukemia cell interface followed by CD19 internalization and decreased CD19 surface expression on the B-ALL cells. CD19 expression was then repressed by transcriptional rewiring. Using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing, we demonstrated that a subset of refractory CD19low cells sustained decreased CD19 expression through transcriptional programs of physiologic B-cell activation and germinal center reaction. Inhibiting B-cell activation programs with the Bruton's tyrosine kinase inhibitor ibrutinib increased the cytotoxicity of CD19 CAR T cells without affecting CAR T-cell viability. These results demonstrate transcriptional plasticity as an underlying mechanism of escape from CAR T cells and highlight the importance of combining CAR T-cell therapy with targeted therapies that aim to overcome this plasticity. See related Spotlight by Zhao and Melenhorst, p. 1040.

Temporally restricted activation of IFNβ signaling underlies response to immune checkpoint therapy in mice

The biological determinants of the response to immune checkpoint blockade (ICB) in cancer remain incompletely understood. Little is known about dynamic biological events that underpin therapeutic efficacy due to the inability to frequently sample tumours in patients. Here, we map the transcriptional profiles of 144 responding and non-responding tumours within two mouse models at four time points during ICB. We find that responding tumours display on/fast-off kinetics of type-I-interferon (IFN) signaling. Phenocopying of this kinetics using time-dependent sequential dosing of recombinant IFNs and neutralizing antibodies markedly improves ICB efficacy, but only when IFNβ is targeted, not IFNα. We identify Ly6C⁺/CD11b⁺ inflammatory monocytes as the primary source of IFNβ and find that active type-I-IFN signaling in tumour-infiltrating inflammatory monocytes is associated with T cell expansion in patients treated with ICB. Together, our results suggest that on/fast-off modulation of IFNβ signaling is critical to the therapeutic response to ICB, which can be exploited to drive clinical outcomes towards response.

Immune checkpoint blockade (ICB) is partially successful as a cancer therapy. Here using mouse models, the authors transcriptionally monitor responding and non-responding tumours showing that responding tumours were associated with transient IFN-β signalling which could promote the anti-tumour response.

ASURAT: functional annotation-driven unsupervised clustering of single-cell transcriptomes

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) analysis reveals heterogeneity and dynamic cell transitions. However, conventional gene-based analyses require intensive manual curation to interpret biological implications of computational results. Hence, a theory for efficiently annotating individual cells remains warranted.

RESULTS

We present ASURAT, a computational tool for simultaneously performing unsupervised clustering and functional annotation of disease, cell type, biological process and signaling pathway activity for single-cell transcriptomic data, using a correlation graph decomposition for genes in database-derived functional terms. We validated the usability and clustering performance of ASURAT using scRNA-seq datasets for human peripheral blood mononuclear cells, which required fewer manual curations than existing methods. Moreover, we applied ASURAT to scRNA-seq and spatial transcriptome datasets for human small cell lung cancer and pancreatic ductal adenocarcinoma, respectively, identifying previously overlooked subpopulations and differentially expressed genes. ASURAT is a powerful tool for dissecting cell subpopulations and improving biological interpretability of complex and noisy transcriptomic data.

AVAILABILITY AND IMPLEMENTATION

ASURAT is published on Bioconductor (https://doi.org/10.18129/B9.bioc.ASURAT). The codes for analyzing data in this article are available at Github (https://github.com/keita-iida/ASURATBI) and figshare (https://doi.org/10.6084/m9.figshare.19200254.v4).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.