Classification of lung adenocarcinoma based on stemness scores in bulk and single cell transcriptomes

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

Identifying cancer subtypes based on embryonic and hematopoietic stem cell signatures in pan-cancer

PURPOSE

Cancer cells with stem cell-like properties may contribute to cancer development and therapy resistance. The advancement of multi-omics technology has sparked interest in exploring cancer stemness from a multi-omics perspective. However, there is a limited number of studies that have attempted to subtype cancer by combining different types of stem cell signatures.

METHODS

In this study, 10,323 cancer specimens from 33 TCGA cancer types were clustered based on the enrichment scores of six stemness gene sets, representing two types of stem cell backgrounds: embryonic stem cells (ESCs) and hematopoietic stem cells (HSCs).

RESULTS

We identified four subtypes of pan-cancer, termed StC1, StC2, StC3 and StC4, which displayed distinct molecular and clinical features, including stemness, genome integrity, intratumor heterogeneity, methylation levels, tumor microenvironment, tumor progression, responses to chemotherapy and immunotherapy, and survival prognosis. Importantly, this subtyping method for pan-cancer is reproducible at the protein level.

CONCLUSION

Our findings indicate that the ESC signature is an adverse prognostic factor in cancer, while the HSC signature and ratio of HSC/ESC signatures are positive prognostic factors. The subtyping of cancer based on ESC and HSC signatures may provide insights into cancer biology and clinical implications of cancer.

Stromal Signals Dominate Gene Expression Signature Scores That Aim to Describe Cancer Cell-intrinsic Stemness or Mesenchymality Characteristics

Epithelial-to-mesenchymal transition (EMT) in cancer cells confers migratory abilities, a crucial aspect in the metastasis of tumors that frequently leads to death. In multiple studies, authors proposed gene expression signatures for EMT, stemness, or mesenchymality of tumors based on bulk tumor expression profiling. However, recent studies suggested that noncancerous cells from the microenvironment or macroenvironment heavily influence such signature profiles. Here, we strengthen these findings by investigating 11 published and frequently referenced gene expression signatures that were proposed to describe EMT-related (EMT, mesenchymal, or stemness) characteristics in various cancer types. By analyses of bulk, single-cell, and pseudobulk expression data, we show that the cell type composition of a tumor sample frequently dominates scores of these EMT-related signatures. A comprehensive, integrated analysis of bulk RNA sequencing (RNA-seq) and single-cell RNA-seq data shows that stromal cells, most often fibroblasts, are the main drivers of EMT-related signature scores. We call attention to the risk of false conclusions about tumor properties when interpreting EMT-related signatures, especially in a clinical setting: high patient scores of EMT-related signatures or calls of "stemness subtypes" often result from low cancer cell content in tumor biopsies rather than cancer cell-specific stemness or mesenchymal/EMT characteristics.

SIGNIFICANCE

Cancer self-renewal and migratory abilities are often characterized via gene module expression profiles, also called EMT or stemness gene expression signatures. Using published clinical tumor samples, cancer cell lines, and single cancer cells, we highlight the dominating influence of noncancer cells in low cancer cell content biopsies on their scores. We caution on their application for low cancer cell content clinical cancer samples with the intent to assign such characteristics or subtypes.

Proteomics-based clustering of lung adenocarcinoma identifies three subtypes with significantly different clinical and molecular features

BACKGROUND

Lung adenocarcinoma (LUAD) is a predominant subtype of lung cancer. Although molecular classification of LUAD has been widely explored, proteomics-based subtyping of LUAD remains scarce.

METHODS

We proposed a subtyping method for LUAD based on the expression profiles of 500 proteins with the largest expression variability across LUAD. Furthermore, we comprehensively compared molecular and clinical features among the LUAD subtypes.

RESULTS

Consensus clustering identified three subtypes of LUAD, namely MtE, DrE, and StE. We demonstrated this subtyping method to be reproducible by analyzing two independent LUAD cohorts. MtE was characterized by high enrichment of metabolic pathways, high EGFR mutation rate, low stemness, proliferation, invasion, metastasis and inflammation signatures, favorable prognosis; DrE was characterized by high enrichment of DNA repair pathways, high TP53 mutation rate, and high levels of genomic instability, stemness, proliferation, and intratumor heterogeneity (ITH); and StE was characterized by high enrichment of stroma-related pathways, high KRAS mutation rate, and low levels of genomic instability.

CONCLUSIONS

The proteomics-based clustering analysis identified three LUAD subtypes with significantly different molecular and clinical properties. The novel subtyping method offers new perspectives on the cancer biology and holds promise in improving the clinical management of LUAD.

The Tumor Stemness Indice mRNAsi can Act as Molecular Typing Tool for Lung Adenocarcinoma

Due to the high heterogeneity, lung adenocarcinoma (LUAD) cannot be distinguished into precise molecular subtypes, thereby resulting in poor therapeutic effect and low 5-year survival rate clinically. Although the tumor stemness score (mRNAsi) has been shown to accurately characterize the similarity index of cancer stem cells (CSCs), whether mRNAsi can serve as an effective molecular typing tool for LUAD isn't reported to date. In this study, we first demonstrate that mRNAsi is significantly correlated with the prognosis and disease degree of LUAD patients, i.e., the higher the mRNAsi, the worse the prognosis and the higher the disease degree. Second, we identify 449 mRNAsi-related genes based on both weighted gene co-expression network analysis (WGCNA) and univariate regression analysis. Third, our results display that 449 mRNAsi-related genes can accurately distinguish the LUAD patients into two molecular subtypes: ms-H subtype (with high mRNAsi) and ms-L subtype (with low mRNAsi), particularly the ms-H subtype has a worse prognosis. Remarkably, significant differences in clinical characteristics, immune microenvironment, and somatic mutation exist between the two molecular subtypes, which might lead to the poorer prognosis of the ms-H subtype patients than that of the ms-L subtype ones. Finally, we establish a prognostic model containing 8 mRNAsi-related genes, which can effectively predict the survival rate of LUAD patients. Taken together, our work provides the first molecular subtype related to mRNAsi in LUAD, and reveals that these two molecular subtypes, the prognostic model and marker genes may have important clinical value for effectively monitoring and treating LUAD patients.