An integrated single-cell transcriptomic dataset for non-small cell lung cancer

Development of a prognostic model for NSCLC based on differential genes in tumour stem cells

Non-small cell lung cancer (NSCLC) constitutes a significant portion of lung cancers and cytotoxic drugs (e.g. cisplatin) are currently the first-line treatment. However, NSCLC has developed resistance to this drug, which limits the therapeutic effect and thus affects prognosis. NSCLC sc-RNA-seq data were downloaded from the GEO database and Ku Leuven Laboratory for Functional Epigenetics, and bulk RNA-seq data were obtained from the TCGA database. The "Seurat" package was employed for scRNA-seq data processing, and the uniform manifold approximation and projection (UMAP) were applied for downscaling and cluster identification. Use the FindAllMarkers function to find differential genes (DEGs) for tumor stem cells. Then, we performed univariate regression analyses on the DEGs to identify potential prognostic genes. We created a machine learning framework based on potential prognostic genes, which combines 10 machine learning methods and their 101 combinations to get the optimal prognostic risk model. The model was evaluated in the training set and validation set. A nomogram was developed to provide physicians with a quantitative tool for prognosis prediction. Finally, we evaluated the expression and functionality of SLC2A1. We discovered 22 cell clusters containing 218379 cells by examining single-cell RNA sequencing datasets (GSE148071, KU_lom, GSE131907, GSE136246, GSE127465). Tumour cells were isolated for subpopulation analysis and 162 differential genes from SOX2_cancer were obtained. After univariate Cox analysis, we found 23 genes with prognostic potential prognostic value and utilized them to develop 101‑combination machine learning computational framework. We eventually picked the best performing 'StepCox[both] + RSF', which includes 8 genes. The model has a relatively high prediction accuracy in both TCGA and GEO datasets. In in vitro investigations, targeted suppression of the SLC2A1 gene resulted in significant reductions in proliferation, invasion and migration in A549 cells. In addition, a significant reduction in cisplatin resistance was seen in A549/DDP cells. The outcomes demonstrated the precision and credibility of the prognostic model for NSCLC, highlighting its potential significance in the treatment and prognosis of individuals affected by this disease. SLC2A1 may become a promising prognostic marker and a potential therapeutic target, offering valuable insights to inform clinical treatment decisions.

Single cell transcriptomic analysis reveals tumor immune infiltration by NK cells gene signature in lung adenocarcinoma

Background

Natural Killer (NK) cells are vital components of the innate immune system, crucial for combating infections and tumor growth, making them pivotal in cancer prognosis and immunotherapy. We sought to understand the diverse characteristics of NK cells within lung adenocarcinoma (LUAD) by conducting single-cell RNA sequencing analyses.

Methods

Using the scRNA-seq dataset for multiple primary lung cancers (MPLCs), we examined two major NK cell groups, NK1 and NK2, comparing the expression profiles of 422 differentially expressed NK signature genes. We identified eight genes (SPON2, PLEKHG3, CAMK2N1, RAB27B, CTBP2, EFHD2, GOLM1, and PLOD1) that distinguish NK1 from NK2 cells. A prognostic signature, the NK gene signature (NKGS) score, was established through LASSO Cox regression. High NKGS scores were linked to poorer overall survival in TCGA-LUAD patients and consistently validated in other datasets (GSE31210 and GSE14814).

Results

Functional analysis revealed an enrichment of genes related to the TGF-β signaling pathway in the high NKGS score group. Moreover, a high NKGS score correlated with an immunosuppressive tumor microenvironment (TME) driven by immune evasion mechanisms. We also observed reduced T-cell receptor (TCR) repertoire diversity in the high-risk NKGS group, indicating a negative association between inflammation and risk score.

Conclusion

This study introduced the innovative NKGS score, differentiating NK1 from NK2 cells. High NKGS scores were associated with the TGF-β pathway and provided insights into LUAD prognosis and immune activities.

Obesity-associated microbiomes instigate visceral adipose tissue inflammation by recruitment of distinct neutrophils

Neutrophils are increasingly implicated in chronic inflammation and metabolic disorders. Here, we show that visceral adipose tissue (VAT) from individuals with obesity contains more neutrophils than in those without obesity and is associated with a distinct bacterial community. Exploring the mechanism, we gavaged microbiome-depleted mice with stool from patients with and without obesity during high-fat or normal diet administration. Only mice receiving high-fat diet and stool from subjects with obesity show enrichment of VAT neutrophils, suggesting donor microbiome and recipient diet determine VAT neutrophilia. A rise in pro-inflammatory CD4+ Th1 cells and a drop in immunoregulatory T cells in VAT only follows if there is a transient spike in neutrophils. Human VAT neutrophils exhibit a distinct gene expression pattern that is found in different human tissues, including tumors. VAT neutrophils and bacteria may be a novel therapeutic target for treating inflammatory-driven complications of obesity, including insulin resistance and colon cancer.

Comprehensive single-cell atlas of the mouse retina

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of cellular heterogeneity by characterizing cell types across tissues and species. While several mouse retinal scRNA-seq datasets exist, each dataset is either limited in cell numbers or focused on specific cell classes, thereby hindering comprehensive gene expression analysis across all retina types. To fill the gap, we generated the largest retinal scRNA-seq dataset to date, comprising approximately 190,000 single cells from C57BL/6J mouse retinas, enriched for rare population cells via antibody-based magnetic cell sorting. Integrating this dataset with public datasets, we constructed the Mouse Retina Cell Atlas (MRCA) for wild-type mice, encompassing over 330,000 cells, characterizing 12 major classes and 138 cell types. The MRCA consolidates existing knowledge, identifies new cell types, and is publicly accessible via CELLxGENE, UCSC Cell Browser, and the Broad Single Cell Portal, providing a user-friendly resource for the mouse retina research community.

Fc-Silent Anti-TIGIT Antibodies Potentiate Antitumor Immunity without Depleting Regulatory T Cells

T-cell immunoreceptor with immunoglobulin and immunoreceptor tyrosine-based inhibitory motif domains (TIGIT) is an inhibitory receptor on immune cells that outcompetes an activating receptor, CD226, for shared ligands. Tumor-infiltrating lymphocytes express TIGIT and CD226 on regulatory T cells (Treg) and on CD8+ T cells with tumor-reactive or exhausted phenotypes, supporting the potential of therapeutically targeting TIGIT to enhance antitumor immunity. To optimize the efficacy of therapeutic antibodies against TIGIT, it is necessary to understand IgG Fc (Fcγ) receptor binding for therapeutic benefit. In this study, we showed that combining Fc-enabled (Fce) or Fc-silent (Fcs) anti-TIGIT with antiprogrammed cell death protein 1 in mice resulted in enhanced control of tumors by differential mechanisms: Fce anti-TIGIT promoted the depletion of intratumoral Treg, whereas Fcs anti-TIGIT did not. Despite leaving Treg numbers intact, Fcs anti-TIGIT potentiated the activation of tumor-specific exhausted CD8+ populations in a lymph node-dependent manner. Fce anti-TIGIT induced antibody-dependent cell-mediated cytotoxicity against human Treg in vitro, and significant decreases in Treg were measured in the peripheral blood of patients with phase I solid tumor cancer treated with Fce anti-TIGIT. In contrast, Fcs anti-TIGIT did not deplete human Treg in vitro and was associated with anecdotal objective clinical responses in two patients with phase I solid tumor cancer whose peripheral Treg frequencies remained stable on treatment. Collectively, these data provide evidence for pharmacologic activity and antitumor efficacy of anti-TIGIT antibodies lacking the ability to engage Fcγ receptor.

SIGNIFICANCE

Fcs-silent anti-TIGIT antibodies enhance the activation of tumor-specific pre-exhausted T cells and promote antitumor efficacy without depleting T regulatory cells.

Single-Cell Informatics for Tumor Microenvironment and Immunotherapy

Cancer comprises malignant cells surrounded by the tumor microenvironment (TME), a dynamic ecosystem composed of heterogeneous cell populations that exert unique influences on tumor development. The immune community within the TME plays a substantial role in tumorigenesis and tumor evolution. The innate and adaptive immune cells "talk" to the tumor through ligand-receptor interactions and signaling molecules, forming a complex communication network to influence the cellular and molecular basis of cancer. Such intricate intratumoral immune composition and interactions foster the application of immunotherapies, which empower the immune system against cancer to elicit durable long-term responses in cancer patients. Single-cell technologies have allowed for the dissection and characterization of the TME to an unprecedented level, while recent advancements in bioinformatics tools have expanded the horizon and depth of high-dimensional single-cell data analysis. This review will unravel the intertwined networks between malignancy and immunity, explore the utilization of computational tools for a deeper understanding of tumor-immune communications, and discuss the application of these approaches to aid in diagnosis or treatment decision making in the clinical setting, as well as the current challenges faced by the researchers with their potential future improvements.

Auranofin repurposing for lung and pancreatic cancer: low CA12 expression as a marker of sensitivity in patient-derived organoids, with potentiated efficacy by AKT inhibition

BACKGROUND

This study explores the repurposing of Auranofin (AF), an anti-rheumatic drug, for treating non-small cell lung cancer (NSCLC) adenocarcinoma and pancreatic ductal adenocarcinoma (PDAC). Drug repurposing in oncology offers a cost-effective and time-efficient approach to developing new cancer therapies. Our research focuses on evaluating AF's selective cytotoxicity against cancer cells, identifying RNAseq-based biomarkers to predict AF response, and finding the most effective co-therapeutic agents for combination with AF.

METHODS

Our investigation employed a comprehensive drug screening of AF in combination with eleven anticancer agents in cancerous PDAC and NSCLC patient-derived organoids (n = 7), and non-cancerous pulmonary organoids (n = 2). Additionally, we conducted RNA sequencing to identify potential biomarkers for AF sensitivity and experimented with various drug combinations to optimize AF's therapeutic efficacy.

RESULTS

The results revealed that AF demonstrates a preferential cytotoxic effect on NSCLC and PDAC cancer cells at clinically relevant concentrations below 1 µM, sparing normal epithelial cells. We identified Carbonic Anhydrase 12 (CA12) as a significant RNAseq-based biomarker, closely associated with the NF-κB survival signaling pathway, which is crucial in cancer cell response to oxidative stress. Our findings suggest that cancer cells with low CA12 expression are more susceptible to AF treatment. Furthermore, the combination of AF with the AKT inhibitor MK2206 was found to be particularly effective, exhibiting potent and selective cytotoxic synergy, especially in tumor organoid models classified as intermediate responders to AF, without adverse effects on healthy organoids.

CONCLUSION

Our research offers valuable insights into the use of AF for treating NSCLC and PDAC. It highlights AF's cancer cell selectivity, establishes CA12 as a predictive biomarker for AF sensitivity, and underscores the enhanced efficacy of AF when combined with MK2206 and other therapeutics. These findings pave the way for further exploration of AF in cancer treatment, particularly in identifying patient populations most likely to benefit from its use and in optimizing combination therapies for improved patient outcomes.

Differentiation signals induce APOBEC3A expression via GRHL3 in squamous epithelia and squamous cell carcinoma

Two APOBEC (apolipoprotein-B mRNA editing enzyme catalytic polypeptide-like) DNA cytosine deaminase enzymes (APOBEC3A and APOBEC3B) generate somatic mutations in cancer, driving tumour development and drug resistance. Here we used single cell RNA sequencing to study APOBEC3A and APOBEC3B expression in healthy and malignant mucosal epithelia, validating key observations with immunohistochemistry, spatial transcriptomics and functional experiments. Whereas APOBEC3B is expressed in keratinocytes entering mitosis, we show that APOBEC3A expression is confined largely to terminally differentiating cells and requires Grainyhead-like transcription factor 3 (GRHL3). Thus, in normal tissue, neither deaminase appears to be expressed at high levels during DNA replication, the cell cycle stage associated with APOBEC-mediated mutagenesis. In contrast, we show that in squamous cell carcinoma tissues, there is expansion of GRHL3 expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings indicate a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

Comprehensive single-cell atlas of the mouse retina

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of cellular heterogeneity at the single-cell resolution by classifying and characterizing cell types in multiple tissues and species. While several mouse retinal scRNA-seq reference datasets have been published, each dataset either has a relatively small number of cells or is focused on specific cell classes, and thus is suboptimal for assessing gene expression patterns across all retina types at the same time. To establish a unified and comprehensive reference for the mouse retina, we first generated the largest retinal scRNA-seq dataset to date, comprising approximately 190,000 single cells from C57BL/6J mouse whole retinas. This dataset was generated through the targeted enrichment of rare population cells via antibody-based magnetic cell sorting. By integrating this new dataset with public datasets, we conducted an integrated analysis to construct the Mouse Retina Cell Atlas (MRCA) for wild-type mice, which encompasses over 330,000 single cells. The MRCA characterizes 12 major classes and 138 cell types. It captured consensus cell type characterization from public datasets and identified additional new cell types. To facilitate the public use of the MRCA, we have deposited it in CELLxGENE, UCSC Cell Browser, and the Broad Single Cell Portal for visualization and gene expression exploration. The comprehensive MRCA serves as an easy-to-use, one-stop data resource for the mouse retina communities.

The collagen landscape in cancer: profiling collagens in tumors and in circulation reveals novel markers of cancer-associated fibroblast subtypes

Cancer-associated fibroblasts (CAFs) deposit and remodel collagens in the tumor stroma, impacting cancer progression and efficacy of interventions. CAFs are the focus of new therapeutics with the aim of normalizing the tumor microenvironment. To do this, a better understanding of CAF heterogeneity and collagen composition in cancer is needed. In this study, we sought to profile the expression of collagens at multiple levels with the goal of identifying cancer biomarkers. We investigated the collagen expression pattern in various cell types and CAF subtypes in a publicly available single-cell RNA sequencing (RNA-seq) dataset of pancreatic ductal adenocarcinoma. Next, we investigated the collagen expression profile in tumor samples across cancer types from The Cancer Genome Atlas (TCGA) database and evaluated if specific patterns of collagen expression were associated with prognosis. Finally, we profiled circulating collagen peptides using a panel of immunoassays to measure collagen fragments in the serum of cancer patients. We found that pancreatic stellate cells and fibroblasts were the primary producers of collagens in the pancreas. COL1A1, COL3A1, COL5A1, COL6A1 were expressed in all CAF subtypes, whereas COL8A1, COL10A1, COL11A1, COL12A1 were specific to myofibroblast CAFs (myCAF) and COL14A1 specific to inflammatory CAFs (iCAF). In TCGA database, myCAF collagens COL10A1 and COL11A1 were elevated across solid tumor types, and multiple associations between high expression and worse survival were found. Finally, circulating collagen biomarkers were elevated in the serum of patients with cancer relative to healthy controls with COL11A1 (myCAF) having the best diagnostic accuracy of the markers measured. In conclusion, CAFs express a noncanonical collagen profile with specific collagen subtypes associated with iCAFs and myCAFs in PDAC. These collagens are deregulated at the cellular, tumor, and systemic levels across different solid tumors and associate with survival. These findings could lead to new discoveries such as novel biomarkers and therapeutic targets. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Genomic Landscape of Pleural Mesothelioma and Therapeutic Aftermaths

PURPOSE OF REVIEW

In this article, we provide a comprehensive analysis of recent progress in the genetic characterisation of pleural mesothelioma, and the translation of these findings to clinical practice.

RECENT FINDINGS

Advancements in sequencing technology have allowed the identification of driver mutations and improved our understanding of how these mutations may shape the mesothelioma tumour microenvironment. However, the identification of frequently mutated regions including CDKN2A, BAP1 and NF2 have, to date, not yet yielded targeted therapy options that outperform standard chemo- and immunotherapies. Similarly, the association between mutational profile and the immune microenvironment or immunotherapy response is not well characterised. Further research into the link between tumour mutational profile and response to therapy is critical for identifying targetable vulnerabilities and stratifying patients for therapy.

M2-like tumor-associated macrophages promote epithelial-mesenchymal transition through the transforming growth factor β/Smad/zinc finger e-box binding homeobox pathway with increased metastatic potential and tumor cell proliferation in lung squamous cell carcinoma

Epithelial-mesenchymal transition (EMT) promotes primary tumor progression toward a metastatic state. The role of tumor-associated macrophages (TAMs) in inducing EMT in lung squamous cell carcinoma (LUSC) remains unclear. We aimed to clarify the significance of TAMs in relation to EMT in LUSC. We collected 221 LUSC specimens from patients who had undergone surgery. Immunohistochemistry was performed to evaluate M1-like and M2-like TAM distribution and EMT by E-cadherin and vimentin staining. Human LUSC cell lines (H226 and EBC-1) and a human monocyte cell line (THP-1) were used for in vitro experiments. M2-like polarization of TAMs and EMT marker expression in LUSC cells were evaluated by western blotting. The biological behavior of LUSC cells was evaluated by migration, invasion, and cell proliferation assays. Immunohistochemical analysis showed that 166 (75.1%) tumors were E-cadherin-positive and 44 (19.9%) were vimentin-positive. M2-like TAM density in the tumor stroma was significantly associated with vimentin positivity and worse overall survival. Western blotting demonstrated higher levels of CD163, CD206, vascular endothelial growth factor, and transforming growth factor beta 1 (TGF-β1) in TAMs versus unstimulated macrophages. Furthermore, increased TGF-β1 secretion from TAMs was confirmed by ELISA. TAM-co-cultured H226 and EBC-1 cells exhibited EMT (decreased E-cadherin, increased vimentin). Regarding EMT-activating transcriptional factors, phosphorylated Smad3 and ZEB-family proteins were higher in TAM-co-cultured LUSC cells than in parental cells. TAM-co-cultured H226 and EBC-1 cells demonstrated enhanced migration and invasion capabilities and improved proliferation. Overall, the present study suggests that TAMs can induce EMT with increased metastatic potential and tumor cell proliferation in LUSC.

CD28/PD1 co-expression: dual impact on CD8+ T cells in peripheral blood and tumor tissue, and its significance in NSCLC patients' survival and ICB response

BACKGROUND

Immune checkpoint blockade (ICB) has significantly prolonged survival of non-small cell lung cancer (NSCLC) patients, although most patients develop mechanisms of resistance. Recently single-cell RNA-sequencing (scRNA-Seq) revealed a huge T-cell phenotypic and (dys)functional state variability. Accordingly, T-cell exhaustion is recognized as a functional adaptation, with a dynamic progression from a long-lived "pre-exhausted stem-like progenitor" to a "terminally exhausted" state. In this scenario it is crucial to understand the complex interplay between co-stimulatory and inhibitory molecules in CD8+ T-cell functionality.

METHODS

To gain a baseline landscape of the composition, functional states, and transcriptomic signatures predictive of prognosis, we analyzed CD8+ T-cell subsets characterized by the presence/absence of PD1 and CD28 from periphery, adjacent non-tumor tissue and tumor site of a cohort of treatment-naïve NSCLC patients, by integrated multiparametric flow cytometry, targeted multi-omic scRNA-seq analyses, and computational pipelines.

RESULTS

Despite the increased PD1 levels, an improved PD1+CD28+ T-cell polyfunctionality was observed with the transition from periphery to tumor site, associated with lack of TIGIT, TIM-3 and LAG-3, but not with Ag-experienced-marker CD11a. Differently from CD28+ T cells, the increased PD1 levels in the tumor were associated with reduced functionality in PD1+CD28- T cells. CD11ahigh, although expressed only in a small fraction of this subset, still sustained its functionality. Absence of TIGIT, TIM-3 and CTLA-4, alone or combined, was beneficial to CD28- T cells. Notably, we observed distinct TRM phenotypes in the different districts, with CD28+ T cells more capable of producing TGFβ in the periphery, potentially contributing to elevated CD103 levels. In contrast CD28- TRM mainly produced CXCL13 within the tumor. ScRNA-seq revealed 5 different clusters for each of the two subsets, with distinctive transcriptional profiles in the three districts. By interrogating the TCGA dataset of patients with lung adenocarcinoma (LUAD) and metastatic NSCLC treated with atezolizumab, we found signatures of heterogeneous TRM and "pre-exhausted" long-lived effector memory CD8+ T cells associated with improved response to ICB only in the presence of CD28.

CONCLUSIONS

Our findings identify signatures able to stratify survival of LUAD patients and predict ICB response in advanced NSCLC. CD28 is advocated as a key determinant in the signatures identified, in both periphery and tumor site, thus likely providing feasible biomarkers of ICB response.

Prognostic Factors and Markers in Non-Small Cell Lung Cancer: Recent Progress and Future Challenges

Lung cancer is a highly aggressive neoplasm and, despite the development of recent therapies, tumor progression and recurrence following the initial response remains unsolved. Several questions remain unanswered about non-small cell lung cancer (NSCLC): (1) Which patients will actually benefit from therapy? (2) What are the predictive factors of response to MAbs and TKIs? (3) What are the best combination strategies with conventional treatments or new antineoplastic drugs? To answer these questions, an integrative literature review was carried out, searching articles in PUBMED, NCBI-PMC, Google Academic, and others. Here, we will examine the molecular genetics of lung cancer, emphasizing NSCLC, and delineate the primary categories of inhibitors based on their molecular targets, alongside the main treatment alternatives depending on the type of acquired resistance. We highlighted new therapies based on epigenetic information and a single-cell approach as a potential source of new biomarkers. The current and future of NSCLC management hinges upon genotyping correct prognostic markers, as well as on the evolution of precision medicine, which guarantees a tailored drug combination with precise targeting.

Single-cell RNA sequencing reveals sexual diversity in the human bladder and its prospective impacts on bladder cancer and urinary tract infection

BACKGROUND

Some bladder-related diseases, such as bladder urinary tract infection (UTI) and bladder cancer (BCa), have significant six differences in incidence and prognosis. However, the molecular mechanisms underlying these sex differences are still not fully understood. Understanding the sex-biased differences in gene expression in normal bladder cells can help resolve these problems.

METHODS

We first collected published single-cell RNA sequencing (scRNA-seq) data of normal human bladders from females and males to map the bladder transcriptomic landscape. Then, Gene Ontology (GO) analysis and gene set enrichment analysis (GSEA) were used to determine the significant pathways that changed in the specific cell populations. The Monocle2 package was performed to reconstruct the differentiation trajectories of fibroblasts. In addition, the scMetabolism package was used to analyze the metabolic activity at the single-cell level, and the SCENIC package was used to analyze the regulatory network.

RESULTS

In total, 27,437 cells passed stringent quality control, and eight main cell types in human bladder were identified according to classical markers. Sex-based differential gene expression profiles were mainly observed in human bladder urothelial cells, fibroblasts, B cells, and T cells. We found that urothelial cells in males demonstrated a higher growth rate. Moreover, female fibroblasts produced more extracellular matrix, including seven collagen genes that may mediate BCa progression. Furthermore, the results showed that B cells in female bladders exhibited more B-cell activated signals and a higher expression of immunoglobulin genes. We also found that T cells in female bladders exhibited more T-cell activated signals. These different biological functions and properties of these cell populations may correlate with sex differences in UTI and BCa, and result in different disease processes and outcomes.

CONCLUSIONS

Our study provides reasonable insights for further studies of sex-based physiological and pathological disparities in the human bladder, which will contribute to the understanding of epidemiological differences in UTI and BCa.