Signature-scoring methods developed for bulk samples are not adequate for cancer single-cell RNA sequencing data

Cancer cell - Fibroblast crosstalk via HB-EGF, EGFR, and MAPK signaling promotes the expression of macrophage chemo-attractants in squamous cell carcinoma

Interactions between cells in the tumor microenvironment (TME) shape cancer progression and patient prognosis. To gain insights into how the TME influences cancer outcomes, we derive gene expression signatures indicative of signaling between stromal fibroblasts and cancer cells, and demonstrate their prognostic significance in multiple and independent squamous cell carcinoma cohorts. By leveraging information within the signatures, we discover that the HB-EGF/EGFR/MAPK axis represents a hub of tumor-stroma crosstalk, promoting the expression of CSF2 and LIF and favoring the recruitment of macrophages. Together, these analyses demonstrate the utility of our approach for interrogating the extent and consequences of TME crosstalk.

Intratumoral radiation dose heterogeneity augments antitumor immunity in mice and primes responses to checkpoint blockade

Radiation therapy (RT) activates multiple immunologic effects in the tumor microenvironment (TME), with diverse dose-response relationships observed. We hypothesized that, in contrast with homogeneous RT, a heterogeneous RT dose would simultaneously optimize activation of multiple immunogenic effects in a single TME, resulting in a more effective antitumor immune response. Using high-dose-rate brachytherapy, we treated mice bearing syngeneic tumors with a single fraction of heterogeneous RT at a dose ranging from 2 to 30 gray. When combined with dual immune checkpoint inhibition in murine models, heterogeneous RT generated more potent antitumor responses in distant, nonirradiated tumors compared with any homogeneous dose. The antitumor effect after heterogeneous RT required CD4 and CD8 T cells and low-dose RT to a portion of the tumor. At the 3-day post-RT time point, dose heterogeneity imprinted the targeted TME with spatial differences in immune-related gene expression, antigen presentation, and susceptibility of tumor cells to immune-mediated destruction. At a later 10-day post-RT time point, high-, moderate-, or low-RT-dose regions demonstrated distinct infiltrating immune cell populations. This was associated with an increase in the expression of effector-associated cytokines in circulating CD8 T cells. Consistent with enhanced adaptive immune priming, heterogeneous RT promoted clonal expansion of effector CD8 T cells. These findings illuminate the breadth of dose-dependent effects of RT on the TME and the capacity of heterogeneous RT to promote antitumor immunity when combined with immune checkpoint inhibitors.

Differential modulation of mutant CALR and JAK2 V617F-driven oncogenesis by HLA genotype in myeloproliferative neoplasms

It has been demonstrated previously that human leukocyte antigen class I (HLA-I) and class II (HLA-II) alleles may modulate JAK2 V617F and CALR mutation (CALRmut)-associated oncogenesis in myeloproliferative neoplasms (MPNs). However, the role of immunogenetic factors in MPNs remains underexplored. We aimed to investigate the potential involvement of HLA genes in CALRmut+ MPNs. High-resolution genotyping of HLA-I and -II loci was conducted in 42 CALRmut+ and 158 JAK2 V617F+ MPN patients and 1,083 healthy controls. A global analysis of the diversity of HLA-I genotypes revealed no significant differences between CALRmut+ patients and controls. However, one HLA-I allele (C*06:02) showed an inverse correlation with presence of CALR mutation. A meta-analysis across independent cohorts and healthy individuals from the 1000 Genomes Project confirmed an inverse correlation between the presentation capabilities of the HLA-I loci for JAK2 V617F and CALRmut-derived peptides in both patients and healthy individuals. scRNA-Seq analysis revealed low expression of TAP1 and CIITA genes in CALRmut+ hematopoietic stem and progenitor cells. In conclusion, the HLA-I genotype differentially restricts JAK2 V617F and CALRmut-driven oncogenesis potentially explaining the mutual exclusivity of the two mutations and differences in their presentation latency. These findings have practical implications for the development of neoantigen-based vaccines in MPNs.

Comparative analysis of single-cell pathway scoring methods and a novel approach

Single-cell gene set analysis (scGSA) provides a useful approach for quantifying molecular functions and pathways in high-throughput transcriptomic data, facilitating the biological interpretation of complex human datasets. However, various factors such as gene set size, quality of the gene sets and the dropouts impact the performance of scGSA. To address these limitations, we present a single-cell Pathway Score (scPS) method to measure gene set activity at single-cell resolution. Furthermore, we benchmark our method with six other methods: AUCell, AddModuleScore, JASMINE, UCell, SCSE and ssGSEA. The comparison across all the methods using two different simulation approaches highlights the effect of cell count, gene set size, noise, condition-specific genes and zero imputation on their performance. The results of our study indicate that the scPS is comparable with other single-cell scoring methods and detects fewer false positives. Importantly, this work reveals critical variables in the scGSA.

Innovative aspects and applications of single cell technology for different diseases

Recent developments in single-cell technologies have provided valuable insights from cancer genomics to complex microbial communities. Single-cell technologies including the RNA-seq, next-generation sequencing (NGS), epigenomics, genomics, and transcriptomics can be used to uncover the single cell nature and molecular characterization of individual cells. These technologies also reveal the cellular transition states, evolutionary relationships between genes, the complex structure of single-cell populations, cell-to-cell interaction leading to biological discoveries and more reliable than traditional bulk technologies. These technologies are becoming the first choice for the early detection of inflammatory biomarkers affecting the proliferation and progression of tumor cells in the tumor microenvironment and improving the clinical efficacy of patients undergoing immunotherapy. These technologies also hold a central position in the detection of checkpoint inhibitors and thus determining the signaling pathways evoked by tumor invasion. This review addressed the emerging approaches of single cell-based technologies in cancer immunotherapies and different human diseases at cellular and molecular levels and the emerging role of sequencing technologies leading to drug discovery. Advancements in these technologies paved for discovering novel diagnostic markers for better understanding the pathological and biochemical mechanisms also for controlling the rate of different diseases.

Polyploid cancer cells reveal signatures of chemotherapy resistance

Therapeutic resistance in cancer significantly contributes to mortality, with many patients eventually experiencing recurrence after initial treatment responses. Recent studies have identified therapy-resistant large polyploid cancer cells in patient tissues, particularly in late-stage prostate cancer, linking them to advanced disease and relapse. Here, we analyzed bone marrow aspirates from 44 advanced prostate cancer patients and found the presence of circulating tumor cells with increased genomic content (CTC-IGC) was significantly associated with poorer progression-free survival. Single cell copy number profiling of CTC-IGC displayed clonal origins with typical CTCs, suggesting complete polyploidization. Induced polyploid cancer cells from PC3 and MDA-MB-231 cell lines treated with docetaxel or cisplatin were examined through single cell DNA sequencing, RNA sequencing, and protein immunofluorescence. Novel RNA and protein markers, including HOMER1, TNFRSF9, and LRP1, were identified as linked to chemotherapy resistance. These markers were also present in a subset of patient CTCs and associated with recurrence in public gene expression data. This study highlights the prognostic significance of large polyploid tumor cells, their role in chemotherapy resistance, and their expression of markers tied to cancer relapse, offering new potential avenues for therapeutic development.

Dysregulated alveolar epithelial cell progenitor function and identity in Hermansky-Pudlak syndrome pulmonary fibrosis

Hermansky-Pudlak syndrome (HPS) is a genetic disorder of endosomal protein trafficking associated with pulmonary fibrosis in specific subtypes, including HPS-1 and HPS-2. Single mutant HPS1 and HPS2 mice display increased fibrotic sensitivity while double mutant HPS1/2 mice exhibit spontaneous fibrosis with aging, which has been attributed to HPS mutations in alveolar epithelial type II (AT2) cells. We utilized HPS mouse models and human lung tissue to investigate mechanisms of AT2 cell dysfunction driving fibrotic remodeling in HPS. Starting at 8 weeks of age, HPS mice exhibited progressive loss of AT2 cell numbers. HPS AT2 cell was impaired ex vivo and in vivo. Incorporating AT2 cell lineage tracing in HPS mice, we observed aberrant differentiation with increased AT2-derived alveolar epithelial type I cells. Transcriptomic analysis of HPS AT2 cells revealed elevated expression of genes associated with aberrant differentiation and p53 activation. Lineage tracing and modeling studies demonstrated that HPS AT2 cells were primed to persist in a Krt8+ reprogrammed transitional state, mediated by p53 activity. Intrinsic AT2 progenitor cell dysfunction and p53 pathway dysregulation are novel mechanisms of disease in HPS-related pulmonary fibrosis, with the potential for early targeted intervention before the onset of fibrotic lung disease.

irGSEA: the integration of single-cell rank-based gene set enrichment analysis

irGSEA is an R package designed to assess the outcomes of various gene set scoring methods when applied to single-cell RNA sequencing data. This package incorporates six distinct scoring methods that rely on the expression ranks of genes, emphasizing relative expression levels over absolute values. The implemented methods include AUCell, UCell, singscore, ssGSEA, JASMINE and Viper. Previous studies have demonstrated the robustness of these methods to variations in dataset size and composition, generating enrichment scores based solely on the relative gene expression of individual cells. By employing the robust rank aggregation algorithm, irGSEA amalgamates results from all six methods to ascertain the statistical significance of target gene sets across diverse scoring methods. The package prioritizes user-friendliness, allowing direct input of expression matrices or seamless interaction with Seurat objects. Furthermore, it facilitates a comprehensive visualization of results. The irGSEA package and its accompanying documentation are accessible on GitHub (https://github.com/chuiqin/irGSEA).

Dynamic Hippo pathway activity underlies mesenchymal differentiation during lung alveolar morphogenesis

Alveologenesis, the final stage in lung development, substantially remodels the distal lung, expanding the alveolar surface area for efficient gas exchange. Secondary crest myofibroblasts (SCMF) exist transiently in the neonatal distal lung and are crucial for alveologenesis. However, the pathways that regulate SCMF function, proliferation and temporal identity remain poorly understood. To address this, we purified SCMFs from reporter mice, performed bulk RNA-seq and found dynamic changes in Hippo-signaling components during alveologenesis. We deleted the Hippo effectors Yap/Taz from Acta2-expressing cells at the onset of alveologenesis, causing a significant arrest in alveolar development. Using single cell RNA-seq, we identified a distinct cluster of cells in mutant lungs with altered expression of marker genes associated with proximal mesenchymal cell types, airway smooth muscle and alveolar duct myofibroblasts. In vitro studies confirmed that Yap/Taz regulates myofibroblast-associated gene signature and contractility. Together, our findings show that Yap/Taz is essential for maintaining functional myofibroblast identity during postnatal alveologenesis.

Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition

Reversible transitions between epithelial and mesenchymal cell states are a crucial form of epithelial plasticity for development and disease progression. Recent experimental data and mechanistic models showed multiple intermediate epithelial-mesenchymal transition (EMT) states as well as trajectories of EMT underpinned by complex gene regulatory networks. In this review, we summarize recent progress in quantifying EMT and characterizing EMT paths with computational methods and quantitative experiments including omics-level measurements. We provide perspectives on how these studies can help relating fundamental cell biology to physiological and pathological outcomes of EMT.

Fallopian tube single cell analysis reveals myeloid cell alterations in high-grade serous ovarian cancer

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

Transcriptional profiling of peripheral blood mononuclear cells identifies inflammatory phenotypes in Ataxia Telangiectasia

INTRODUCTION

Ataxia telangiectasia (A-T) is an autosomal recessive neurodegenerative disease with widespread systemic manifestations and marked variability in clinical phenotypes. In this study, we sought to determine whether transcriptomic profiling of peripheral blood mononuclear cells (PBMCs) defines subsets of individuals with A-T beyond mild and classic phenotypes, enabling identification of novel features for disease classification and treatment response to therapy.

METHODS

Participants with classic A-T (n = 77), mild A-T (n = 13), and unaffected controls (n = 15) were recruited from two outpatient clinics. PBMCs were isolated and bulk RNAseq was performed. Plasma was also isolated in a subset of individuals. Affected individuals were designated mild or classic based on ATM mutations and clinical and laboratory features.

RESULTS

People with classic A-T were more likely to be younger and IgA deficient and to have higher alpha-fetoprotein levels and lower % forced vital capacity compared to individuals with mild A-T. In classic A-T, the expression of genes required for V(D)J recombination was lower, and the expression of genes required for inflammatory activity was higher. We assigned inflammatory scores to study participants and found that inflammatory scores were highly variable among people with classic A-T and that higher scores were associated with lower ATM mRNA levels. Using a cell type deconvolution approach, we inferred that CD4 + T cells and CD8 + T cells were lower in number in people with classic A-T. Finally, we showed that individuals with classic A-T exhibit higher SERPINE1 (PAI-1) mRNA and plasma protein levels, irrespective of age, and higher FLT4 (VEGFR3) and IL6ST (GP130) plasma protein levels compared with mild A-T and controls.

CONCLUSION

Using a transcriptomic approach, we identified novel features and developed an inflammatory score to identify subsets of individuals with different inflammatory phenotypes in A-T. Findings from this study could be used to help direct treatment and to track treatment response to therapy.

Single-cell sequencing analysis within biologically relevant dimensions

The currently predominant approach to transcriptomic and epigenomic single-cell analysis depends on a rigid perspective constrained by reduced dimensions and algorithmically derived and annotated clusters. Here, we developed Seqtometry (sequencing-to-measurement), a single-cell analytical strategy based on biologically relevant dimensions enabled by advanced scoring with multiple gene sets (signatures) for examination of gene expression and accessibility across various organ systems. By utilizing information only in the form of specific signatures, Seqtometry bypasses unsupervised clustering and individual annotations of clusters. Instead, Seqtometry combines qualitative and quantitative cell-type identification with specific characterization of diverse biological processes under experimental or disease conditions. Comprehensive analysis by Seqtometry of various immune cells as well as other cells from different organs and disease-induced states, including multiple myeloma and Alzheimer's disease, surpasses corresponding cluster-based analytical output. We propose Seqtometry as a single-cell sequencing analysis approach applicable for both basic and clinical research.

Single-Cell RNA Sequencing and Transcriptome Analysis Revealed the Immune Microenvironment and Gene Markers of Acute Respiratory Distress Syndrome

Background

Acute respiratory distress syndrome (ARDS) is caused by severe pulmonary inflammation and the leading cause of death in the intensive care unit.

Methods

We used single-cell RNA sequencing to compare peripheral blood mononuclear cells from sepsis-induced ARDS (SEP-ARDS) and pneumonic ARDS (PNE-ARDS) patient. Then, we used the GSE152978 and GSE152979 datasets to identify molecular dysregulation mechanisms at the transcriptional level in ARDS.

Results

Markedly increased CD14 cells were the predominant immune cell type observed in SEP-ARDS and PNE-ARDS patients. Cytotoxic cells and natural killer (NK) T cells were exclusively identified in patients with PNE-ARDS. An enrichment analysis of differentially expressed genes (DEGs) suggested that Th1 cell differentiation and Th2 cell differentiation were enriched in cytotoxic cells, and that the IL-17 signaling pathway, NOD receptor signaling pathway, and complement and coagulation cascades were enriched in CD14 cells. Furthermore, according to GSE152978 and GSE152979, 1939 DEGs were identified in patients with ARDS and controls; they were mainly enriched in the Kyoto Encyclopedia of Genes and Genomes pathways. RBP7 had the highest area under the curve values among the 12 hub genes and was mainly expressed in CD14 cells. Additionally, hub genes were negatively correlated with NK cells and positively correlated with neutrophils, cytotoxic cells, B cells, and macrophages.

Conclusion

A severe imbalance in the proportion of immune cells and immune dysfunction were observed in SEP-ARDS and PNE-ARDS patients. RBP7 may be immunologically associated with CD14 cells and serve as a potential marker of ARDS.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Single-cell gene set enrichment analysis and transfer learning for functional annotation of scRNA-seq data

Although an essential step, cell functional annotation often proves particularly challenging from single-cell transcriptional data. Several methods have been developed to accomplish this task. However, in most cases, these rely on techniques initially developed for bulk RNA sequencing or simply make use of marker genes identified from cell clustering followed by supervised annotation. To overcome these limitations and automatize the process, we have developed two novel methods, the single-cell gene set enrichment analysis (scGSEA) and the single-cell mapper (scMAP). scGSEA combines latent data representations and gene set enrichment scores to detect coordinated gene activity at single-cell resolution. scMAP uses transfer learning techniques to re-purpose and contextualize new cells into a reference cell atlas. Using both simulated and real datasets, we show that scGSEA effectively recapitulates recurrent patterns of pathways' activity shared by cells from different experimental conditions. At the same time, we show that scMAP can reliably map and contextualize new single-cell profiles on a breast cancer atlas we recently released. Both tools are provided in an effective and straightforward workflow providing a framework to determine cell function and significantly improve annotation and interpretation of scRNA-seq data.

Cell-type deconvolution of bulk RNA-Seq from kidney using opensource bioinformatic tools

Traditional bulk RNA-Seq pipelines do not assess cell-type composition within heterogeneous tissues. Therefore, it is difficult to determine whether conflicting findings among samples or datasets are the result of biological differences or technical differences due to variation in sample collections. This report provides a user-friendly, open source method to assess cell-type composition in bulk RNA-Seq datasets for heterogeneous tissues using published single cell (sc)RNA-Seq data as a reference. As an example, we apply the method to analysis of kidney cortex bulk RNA-Seq data from female (N=8) and male (N=9) baboons to assess whether observed transcriptome sex differences are biological or technical, i.e., variation due to ultrasound guided biopsy collections. We found cell-type composition was not statistically different in female versus male transcriptomes based on expression of 274 kidney cell-type specific transcripts, indicating differences in gene expression are not due to sampling differences. This method of cell-type composition analysis is recommended for providing rigor in analysis of bulk RNA-Seq datasets from complex tissues. It is clear that with reduced costs, more analyses will be done using scRNA-Seq; however, the approach described here is relevant for data mining and meta analyses of the thousands of bulk RNA-Seq data archived in the NCBI GEO public database.

Protocol to benchmark gene expression signature scoring techniques for single-cell RNA sequencing data in cancer

Scoring gene signatures is common for bulk and single-cell RNA sequencing (scRNAseq) data. Here, using cancer as a data model, we describe steps to benchmark signature scoring techniques for scRNAseq data in the context of uneven gene dropouts. These steps include identifying and comparing deregulated signatures, generating gold standard signatures for specificity and sensitivity tests, and simulating the impact of dropouts using down sampling. The protocol provides a framework for benchmarking scRNAseq algorithms in such context. For complete details on the use and execution of this protocol, please refer to Noureen et al. (2022).1.

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.

A machine learning framework for scRNA-seq UMI threshold optimization and accurate classification of cell types

Recent advances in single cell RNA sequencing (scRNA-seq) technologies have been invaluable in the study of the diversity of cancer cells and the tumor microenvironment. While scRNA-seq platforms allow processing of a high number of cells, uneven read quality and technical artifacts hinder the ability to identify and classify biologically relevant cells into correct subtypes. This obstructs the analysis of cancer and normal cell diversity, while rare and low expression cell populations may be lost by setting arbitrary high cutoffs for UMIs when filtering out low quality cells. To address these issues, we have developed a novel machine-learning framework that: 1. Trains cell lineage and subtype classifier using a gold standard dataset validated using marker genes 2. Systematically assess the lowest UMI threshold that can be used in a given dataset to accurately classify cells 3. Assign accurate cell lineage and subtype labels to the lower read depth cells recovered by setting the optimal threshold. We demonstrate the application of this framework in a well-curated scRNA-seq dataset of breast cancer patients and two external datasets. We show that the minimum UMI threshold for the breast cancer dataset could be lowered from the original 1500 to 450, thereby increasing the total number of recovered cells by 49%, while achieving a classification accuracy of >0.9. Our framework provides a roadmap for future scRNA-seq studies to determine optimal UMI threshold and accurately classify cells for downstream analyses.

Angiogenesis goes computational - The future way forward to discover new angiogenic targets?

Multi-omics technologies are being increasingly utilized in angiogenesis research. Yet, computational methods have not been widely used for angiogenic target discovery and prioritization in this field, partly because (wet-lab) vascular biologists are insufficiently familiar with computational biology tools and the opportunities they may offer. With this review, written for vascular biologists who lack expertise in computational methods, we aspire to break boundaries between both fields and to illustrate the potential of these tools for future angiogenic target discovery. We provide a comprehensive survey of currently available computational approaches that may be useful in prioritizing candidate genes, predicting associated mechanisms, and identifying their specificity to endothelial cell subtypes. We specifically highlight tools that use flexible, machine learning frameworks for large-scale data integration and gene prioritization. For each purpose-oriented category of tools, we describe underlying conceptual principles, highlight interesting applications and discuss limitations. Finally, we will discuss challenges and recommend some guidelines which can help to optimize the process of accurate target discovery.