Tumor intrinsic immunity related proteins may be novel tumor suppressors in some types of cancer

Growth differentiation factor-15 may be a novel biomarker in pancreatic cancer: A review

Pancreatic cancer is a highly malignant and invasive gastrointestinal tumor that is often diagnosed at an advanced stage with a poor prognosis and high mortality. Currently, carbohydrate antigen199(CA199) is the only biomarker approved by the FDA for the diagnosis of pancreatic cancer, but it has great limitations. Growth differentiation factor-15 (GDF-15) is expected to be a novel biomarker for the diagnosis, efficacy prediction, and prognosis assessment of pancreatic cancer patients. In this paper, we searched the keywords GDF-15, macrophage inhibitory cytokine-1 (MIC-1), CA199, pancreatic cancer, and tumor markers in PubMed and Web of Science, searched related articles, and read and analyzed the retrieved papers. Finally, we systematically described the characteristics, mechanism of action, and clinical value of GDF-15, aiming to provide help for the detection and treatment of pancreatic cancer.

Integrated transcriptomics and proteomics data analysis identifies CDH17 as a key cell surface target in colorectal cancer

Immunotherapy development against colorectal cancer (CRC) is hindered by the lack of cell surface target highly expressed in cancer cells but with restricted presence in normal tissues to minimize off-tumor toxicities. In this in silico analysis, a longlist of genes (n = 13,488) expressed in CRCs according to the Human Protein Atlas (HPA) database were evaluated to shortlist for potential surface targets based on the following prerequisites: (i) Absent from the brain and lung tissues to minimize the likelihood of neurologic and pulmonary toxicities; (ii) Restricted expression profile in other normal human tissues; (iii) Genes that potentially encode cell surface proteins and; (iv) At least moderately expressed in CRC cases. Fifteen potential targets were shortlisted and subsequently ranked according to the combination of their transcript and protein expression levels in CRCs derived from multiple datasets (i.e. DepMap, TCGA, CPTAC-2, and HPA CRCs). The top-ranked target with the highest and homogenous expression in CRCs was cadherin 17 (CDH17). Downstream analysis of CRC transcriptomics and proteomics datasets showed that CDH17 was significantly correlated with carcinoembryonic antigen expression. Moreover, CDH17 expression was significantly lower in CRC cases with high microsatellite instability, as well as negatively associated with immune response gene sets and the expression of MHC class I and II molecules. CDH17 represents an optimal target for therapeutic development against CRCs, and this study provides a novel framework to identify key cell surface targets for therapeutic development against other malignancies.

PES1 reduces CD8+ T cell infiltration and immunotherapy sensitivity via interrupting ILF3-IL15 complex in esophageal squamous cell carcinoma

BACKGROUND

Although immune checkpoint blockade (ICB) therapy has brought survival benefits to patients with specific cancer types, most of cancer patients remain refractory to the ICB therapy, which is largely attributed to the immunosuppressive tumor microenvironment. Thereby, it is urgent to profile key molecules and signal pathways responsible for modification of tumor microenvironment.

METHODS

Multiple databases of esophageal squamous cell carcinoma (ESCC) were integratively analyzed to screen candidate genes responsible for infiltration of CD8⁺ T cells. Expression of pescadillo ribosomal biogenesis factor 1 (PES1) in clinical ESCC samples was examined by qRT-PCR, western blotting, and immunohistochemistry. The mechanisms of PES1 were investigated via RNA sequencing and mass spectrometry followed by immunoprecipitation and proximity ligation assay. The clinical and therapeutic significance of PES1 in ESCC was comprehensively investigated using ESCC cells and mouse model.

RESULTS

PES1 was significantly upregulated and correlated with poor prognosis in ESCC patients. PES1 knockdown decreased ESCC cell growth in vitro and in vivo and enhanced the efficacy of ICB therapy in mouse model, which was established through subcutaneous inoculation with ESCC cells. Analyses on RNA sequencing and mass spectrometry suggested that PES1 expression was negatively correlated with IL15 and ILF3 was one of the PES1-associated proteins. It has been known that ILF3 interacts with and stabilizes IL15 mRNA to increase IL15 protein level. Our data further indicated that PES1 interfered with the interaction between ILF3 and IL15 mRNA and impaired ILF3-mediated stabilization of IL15 mRNA, which eventually reduced the protein level of IL15. Interestingly, the inhibitory effect of ICB therapy boosted by PES1 knockdown dramatically antagonized by knockdown of IL15, which suppressed the tumor-infiltrated CD8⁺ T cells in ESCC. Finally, we confirmed the relationships among PES1, IL15, and CD8⁺ T cell infiltration in 10 locally advanced ESCC patients receiving ICB neoadjuvant therapy and demonstrated that ICB therapy would be more effective in those with low expression of PES1.

CONCLUSIONS

Altogether, our findings herein provided novel insights on biological function and clinical significance of PES1 and suggested that high expression of PES1 could suppress ILF3-IL15 axis-mediated immunosurveillance and promote resistance to ICB through restraining tumor-infiltrated CD8⁺ T cells.

Identification and development of a novel risk model based on cuproptosis-associated RNA methylation regulators for predicting prognosis and characterizing immune status in hepatocellular carcinoma

BACKGROUND

Cuproptosis, a novel cell death caused by excess copper, is quite obscure in hepatocellular carcinoma (HCC) and needs more investigation.

METHODS

RNA-seq and clinical data of HCC patients TCGA database were analyzed to establish a predictive model through LASSO Cox regression analysis. External dataset ICGC was used for the verification. GSEA and CIBERSORT were applied to investigate the molecular mechanisms and immune microenvironment of HCC. Cuproptosis induced by elesclomol was confirmed via various in vitro experiments. The expression of prognostic genes was verified in HCC tissues using qRT-PCR analysis.

RESULTS

Initially, 18 cuproptosis-associated RNA methylation regulators (CARMRs) were selected for prognostic analysis. A nine-gene signature was created by applying the LASSO Cox regression method. Survival and ROC assays were carried out to validate the model using TCGA and ICGC database. Moreover, there exhibited obvious differences in drug sensitivity in terms of common drugs. A higher tumor mutation burden was shown in the high-risk group. Additionally, significant discrepancies were found between the two groups in metabolic pathways and RNA processing via GSEA analysis. Meanwhile, CIBERSORT analysis indicated different infiltrating levels of various immune cells between the two groups. Elesclomol treatment caused a unique form of programmed cell death accompanied by loss of lipoylated mitochondrial proteins and Fe-S cluster protein. The results of qRT-PCR indicated that most prognostic genes were differentially expressed in the HCC tissues.

CONCLUSION

Overall, our predictive signature displayed potential value in the prediction of overall survival of HCC patients and might provide valuable clues for personalized therapies.

A novel chromatin regulator-related immune checkpoint related gene prognostic signature and potential candidate drugs for endometrial cancer patients

Background

Endometrial cancer (EC) is the most common gynecologic malignancy in developed countries and its prevalence is increasing. As an emerging therapy with a promising efficacy, immunotherapy has been extensively applied in the treatment of solid tumors. In addition, chromatin regulators (CRs), as essential upstream regulators of epigenetics, play a significant role in tumorigenesis and cancer development.

Methods

CRs and immune checkpoint-related genes (ICRGs) were obtained from the previous top research. The Genome Cancer Atlas (TCGA) was utilized to acquire the mRNA expression and clinical information of patients with EC. Correlation analysis was utilized for screen CRs-related ICRGs (CRRICRGs). By Cox regression and least absolute shrinkage and selection operator (LASSO) analysis, prognosis related CRRICRGs were screened out and risk model was constructed. The Kaplan–Meier curve was used to estimate the prognosis between high- and low-risk group. By comparing the IC50 value, the drugs sensitivity difference was explored. We obtained small molecule drugs for the treatment of UCEC patients based on CAMP dataset.

Results

We successfully constructed a 9 CRRICRs-based prognostic signature for patients with UCEC and found the riskscore was an independent prognostic factor. The results of functional analysis suggested that CRRICRGs may be involved in immune processes associated with cancer. Immune characteristics analysis provided further evidence that the CRRICRGs-based model was correlated with immune cells infiltration and immune checkpoint. Eight small molecule drugs that may be effective for the treatment of UCEC patients were screened. Effective drugs identified by drug sensitivity profiling in high- and low-risk groups.

Conclusion

In summary, our study provided novel insights into the function of CRRICRGs in UCEC. We also developed a reliable prognostic panel for the survival of patients with UCEC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00253-w.

Pancancer Analysis of Revealed TDO2 as a Biomarker of Prognosis and Immunotherapy

Background

Tryptophan 2,3-dioxygenase (TDO) encoded by TDO2, a rate-limiting enzyme in the kynurenine pathway, catabolizes tryptophan to kynurenine, evades immune surveillance, and promotes tumor growth. Although accumulating evidence suggests a crucial role of TDO2 during tumor formation and development, systematic evaluation of TDO2 across human cancers has rarely been reported.

Methods

To shed more light on the role of TDO2 in human cancer, we explored the expression profiles of TDO2 and identified its prognostic value in pancancer analysis through TCGA, CCLE, and GTEx databases. We further utilized TCGA data to evaluate the association between TDO2 and tumor immunological features, such as mismatch repair (MMR), tumor immune infiltration, immune checkpoint-related genes, tumor mutational burden (TMB), microsatellite instability (MSI), and DNA methyltransferase (DNMT).

Results

TDO2 exhibited different expression levels in various cancer cell lines. Frequently, TDO2 was detected to be highly expressed in the majority of cancers. In addition, high TDO2 expression was correlated with an unfavorable prognosis for patients in KIRP, LGG, TGCT, and UVM. Moreover, high TDO2 expression level positively correlated with higher immune infiltration, especially dendritic cells. Additionally, there is a close relationship between TDO2 and immune checkpoint-related gene markers, such as LAIR1, CD276, NRP1, CD80, and CD86. Finally, correlation analysis has demonstrated a high-correlation between TDO2 and TMB, MSI, MMR, and DNMT of multiple cancer types.

Conclusion

Therefore, our results suggest that TDO2 can function as a potential prognostic biomarker due to its role in tumor immunity regulation.

Aberrant promoter methylation contributes to LRIG1 silencing in basal/triple-negative breast cancer

BACKGROUND

LRIG1, the founding member of the LRIG (leucine-rich repeat and immunoglobulin-like domain) family of transmembrane proteins, is a negative regulator of receptor tyrosine kinases and a tumour suppressor. Decreased LRIG1 expression is consistently observed in cancer, across diverse tumour types, and is linked to poor patient prognosis. However, mechanisms by which LRIG1 is repressed are not fully understood. Silencing of LRIG1 through promoter CpG island methylation has been reported in colorectal and cervical cancer but studies in breast cancer remain limited.

METHODS

In silico analysis of human breast cancer patient data were used to demonstrate a correlation between DNA methylation and LRIG1 silencing in basal/triple-negative breast cancer, and its impact on patient survival. LRIG1 gene expression, protein abundance, and methylation enrichment were examined by quantitative reverse-transcription PCR, immunoblotting, and methylation immunoprecipitation, respectively, in breast cancer cell lines in vitro. We examined the impact of global demethylation on LRIG1 expression and methylation enrichment using 5-aza-2'-deoxycytidine. We also examined the effects of targeted demethylation of the LRIG1 CpG island, and transcriptional activation of LRIG1 expression, using the RNA guided deadCas9 transactivation system.

RESULTS

Across breast cancer subtypes, LRIG1 expression is lowest in the basal/triple-negative subtype so we investigated whether differential methylation may contribute to this. Indeed, we find that LRIG1 CpG island methylation is most prominent in basal/triple-negative cell lines and patient samples. Use of the global demethylating agent 5-aza-2'-deoxycytidine decreases methylation leading to increased LRIG1 transcript expression in basal/triple-negative cell lines, while having no effect on LRIG1 expression in luminal/ER-positive cell lines. Using a CRISPR/deadCas9 (dCas9)-based targeting approach, we demonstrate that TET1-mediated demethylation (Tet1-dCas9) along with VP64-mediated transcriptional activation (VP64-dCas9) at the CpG island, increased endogenous LRIG1 expression in basal/triple-negative breast cancer cells, without transcriptional upregulation at predicted off-target sites. Activation of LRIG1 by the dCas9 transactivation system significantly increased LRIG1 protein abundance, reduced site-specific methylation, and reduced cancer cell viability. Our findings suggest that CRISPR-mediated targeted activation may be a feasible way to restore LRIG1 expression in cancer.

CONCLUSIONS

Our study contributes novel insight into mechanisms which repress LRIG1 in triple-negative breast cancer and demonstrates for the first time that targeted de-repression of LRIG1 in cancer cells is possible. Understanding the epigenetic mechanisms associated with repression of tumour suppressor genes holds potential for the advancement of therapeutic approaches.

Modulation of Lymphocyte Functions in the Microenvironment by Tumor Oncogenic Pathways

Despite the broad application of different immunotherapeutic strategies for the treatment of solid as well as hematopoietic cancers, the efficacy of these therapies is still limited, with only a minority of patients having a long-term benefit resulting in an improved survival rate. In order to increase the response rates of patients to the currently available immunotherapies, a better understanding of the molecular mechanisms underlying the intrinsic and/or extrinsic resistance to treatment is required. There exist increasing evidences that activation of different oncogenic pathways as well as inactivation of tumor suppressor genes (TSG) in tumor cells inhibit the immune cell recognition and influegnce the composition of the tumor microenvironment (TME), thus leading to an impaired anti-tumoral immune response. A deeper understanding of the link between the tumor milieu and genomic alterations of TSGs and oncogenes is indispensable for the optimization of immunotherapies and to predict the patients’ response to these treatments. This review summarizes the role of different cancer-related, oncogene- and TSG-controlled pathways in the context of anti-tumoral immunity and response to different immunotherapies.

Roles and Regulation of Growth differentiation factor-15 in the Immune and tumor microenvironment

Growth differentiation factor-15 (GDF-15), a member of the TGF-β superfamily, plays multiple roles in a wide variety of cellular processes. It is expressed at low levels under normal conditions but is highly expressed in tumor and tumor microenvironment (TME)-related cells, such as fibroblasts and immune cells. The TME consists of the noncancerous cells present in the tumor, including immune cells, fibroblasts, blood vessel signaling molecules and extracellular matrix, which play a key role in tumor development. GDF-15 affects both stromal cells and immune cells in the TME. It also acts on immune checkpoints, such as PD-1/PDL-1 that regulate stemness of cancer cells, indicating that GDF-15 plays a prominent role in cancer, exhibiting both protumorigenic and antitumorigenic effects, although the latter are reported much less often than the former. The present review addresses novel ideas regarding communication between GDF-15 and stromal cells, immune cells, and cancer cells in the TME. In addition, it discusses the possibility of GDF-15's clinical application as a diagnostic biomarker and therapeutic target in cancer.

Clinical Research on the Mechanisms Underlying Immune Checkpoints and Tumor Metastasis

This study highlights aspects of the latest clinical research conducted on the relationship between immune checkpoints and tumor metastasis. The overview of each immune checkpoint is divided into the following three sections: 1) structure and expression; 2) immune mechanism related to tumor metastasis; and 3) clinical research related to tumor metastasis. This review expands on the immunological mechanisms of 17 immune checkpoints, including TIM-3, CD47, and OX-40L, that mediate tumor metastasis; evidence shows that most of these immune checkpoints are expressed on the surface of T cells, which mainly exert immunomodulatory effects. Additionally, we have summarized the roles of these immune checkpoints in the diagnosis and treatment of metastatic tumors, as these checkpoints are considered common predictors of metastasis in various cancers such as prostate cancer, non-Hodgkin lymphoma, and melanoma. Moreover, certain immune checkpoints can be used in synergy with PD-1 and CTLA-4, along with the implementation of combination therapies such as LIGHT-VTR and anti-PD-1 antibodies. Presently, most monoclonal antibodies generated against immune checkpoints are under investigation as part of ongoing preclinical or clinical trials conducted to evaluate their efficacy and safety to establish a better combination treatment strategy; however, no significant progress has been made regarding monoclonal antibody targeting of CD28, VISTA, or VTCN1. The application of immune checkpoint inhibitors in early stage tumors to prevent tumor metastasis warrants further evidence; the immune-related adverse events should be considered before combination therapy. This review aims to elucidate the mechanisms of immune checkpoint and the clinical progress on their use in metastatic tumors reported over the last 5 years, which may provide insights into the development of novel therapeutic strategies that will assist with the utilization of various immune checkpoint inhibitors.

Cancer Immune Evasion Through Loss of MHC Class I Antigen Presentation

Major histocompatibility class I (MHC I) molecules bind peptides derived from a cell's expressed genes and then transport and display this antigenic information on the cell surface. This allows CD8 T cells to identify pathological cells that are synthesizing abnormal proteins, such as cancers that are expressing mutated proteins. In order for many cancers to arise and progress, they need to evolve mechanisms to avoid elimination by CD8 T cells. MHC I molecules are not essential for cell survival and therefore one mechanism by which cancers can evade immune control is by losing MHC I antigen presentation machinery (APM). Not only will this impair the ability of natural immune responses to control cancers, but also frustrate immunotherapies that work by re-invigorating anti-tumor CD8 T cells, such as checkpoint blockade. Here we review the evidence that loss of MHC I antigen presentation is a frequent occurrence in many cancers. We discuss new insights into some common underlying mechanisms through which some cancers inactivate the MHC I pathway and consider some possible strategies to overcome this limitation in ways that could restore immune control of tumors and improve immunotherapy.

Integrative Genomic and Transcriptomic Analyses of Tumor Suppressor Genes and Their Role on Tumor Microenvironment and Immunity in Lung Squamous Cell Carcinoma

Non-small-cell lung cancers (NSCLCs) are largely classified into lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), which have different therapeutic options according to its molecular profiles and immune checkpoint expression, especially PD-L1, which is a suppressive factor in the tumor microenvironment. The tumor microenvironment can be altered by the genomic mutations on specific innate immune genes as well as tumor suppressor genes, so it is essential to comprehend the association between tumor microenvironment and tumor suppressor genes to discover the promising immunotherapeutic strategy to overcome the resistance of immune check point blockade. In this study, we aimed to analyze how the somatic mutations in tumor suppressor genes affect the tumor immune microenvironment through a comprehensive analysis of mutational profiling on the representative tumor suppressor genes (TP53, CDKN2A, PTEN, RB1, BRCA1, BRCA2) and immune gene expression in The Cancer Genome Atlas (TCGA) 155 lung squamous cell carcinoma (LUSC) and 196 lung adenocarcinoma (LUAD) samples. Several microenvironmental factors, such as the infiltrating immune and stromal cells, were suppressed by the mutated tumor suppressor genes in LUSC, unlike in the LUAD samples. In particular, infiltrating immune cells such as macrophage, neutrophil, and dendritic cells were significantly reduced in tumors with mutated tumor suppressor genes' group. In addition, the gene expressions for interleukin production and lymphocyte differentiation and PGC, C7, HGF, PLA2G2A, IL1RL1, CCR2, ALOX15B, CXCL11, FCN3 were significantly down-regulated, which were key immune genes for the cross-talk between LUSC microenvironment and tumor suppressors. Therefore, we generated evidence that TSG mutations in LUSC have an impact on tumor immune microenvironment, which suggests that TSG non-mutated patients will have the more inflamed tumors and are more likely to respond to immune checkpoint blockade therapy.

Clinical Significance of Down-Regulated CD70 and CD27 Expression in Poor Prognosis of Esophageal Squamous Cell Carcinoma

Introduction

CD27 is a co-stimulatory immune checkpoint molecule in the tumor necrosis factor receptor superfamily. CD27 regulates the generation and maintenance of T cell immunity by binding to CD70 and regulating B-cell activation and immunoglobulin synthesis.

Materials and Methods

CD27 and CD70 expression were assessed in esophageal squamous cell carcinoma (ESCC) compared to normal tissue samples in the GSE53625 dataset of 179 paired cases and in 153 Chinese cases using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry. The correlation was also investigated between CD27 and CD70 expression and immune-related pathways, including CD8⁺ T cell recruitment, function, and other inhibitory immune checkpoints.

Results

Levels of both CD27 and CD70 expression were down-regulated in ESCC compared to the paired normal tissues. CD27 and CD70 expression was mainly present in lymphocytes surrounding and infiltrating the tumor lesions but rarely expressed in tumor cells. Lost expression of CD27 and CD70 was associated with clinicopathological features, including depth of tumor invasion and better patient survival. Furthermore, CD27 expression was significantly associated with levels of CD8A, GZMB, IFNG, the CD8⁺ T cell recruitment-associated chemokines (CXCL9, CXCL10, and CXCL11), and CD8 receptors (CCR5, CXCR6, and CXCR3), while CD70 expression was inversely associated with levels of immunosuppressive checkpoints (PD-L1, PD-L2, and HHLA2).

Conclusion

Detection of CD70/CD27 expression could be further verified as a biomarker for ESCC early detection and prognosis prediction.

Molecular, clinicopathological, and immune correlates of LAG3 promoter DNA methylation in melanoma

Background

The co-receptor lymphocyte-activation gene-3 (LAG3, LAG-3, CD223) is a potential target for immune checkpoint inhibition immunotherapies. However, little is known about the biological and clinical significance of LAG3 DNA methylation in melanoma and its microenvironment.

Methods

We evaluated LAG3 promoter and gene body methylation in a cohort of N = 470 melanoma patients obtained from The Cancer Genome Atlas (TCGA cohort), an independent cohort of N = 120 patients from the University Hospital Bonn, and in subsets of peripheral blood leukocytes, melanocytes, and melanoma cell lines. We validated the association of LAG3 methylation with mRNA expression in vitro in the melanoma cell line A375 treated with the hypomethylating agent 5-azacytidine and stimulated with interferon-γ. Finally, we investigated correlations between LAG3 methylation and progression-free survival in patients treated with immune checkpoint blockade (ICB cohort, N = 118).

Findings

Depending on the analysed locus (promoter, gene body) we found region-dependent significant LAG3 methylation differences between monocytes, B cells, CD8⁺ and CD4⁺ T cells, regulatory T cells, melanocytes, and melanoma cell lines. In tumor tissues, methylation correlated significantly with LAG3 mRNA expression, immune cell infiltrates (histopathologic lymphocyte score and RNA-Seq signatures of distinct immune infiltrates), and an interferon-γ signature. Finally, LAG3 methylation was associated with overall survival in the TCGA cohort and progression-free survival in the ICB cohort. We detected basal LAG3 mRNA expression in the melanoma cell A375 and an interferon-γ inducible expression after demethylation with 5-azacytidine.

Interpretation

Our study points towards an epigenetic regulation of LAG3 via promoter methylation and suggests a prognostic and predictive significance of LAG3 methylation in melanoma. Our results give insight in the tumor cell-intrinsic transcriptional regulation of LAG3 in melanoma. In perspective, our results might pave the way for investigating LAG3 methylation as a predictive biomarker for response to anti-LAG3 immune checkpoint blockage.

Funding

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Genome-Wide Identification and Analysis of Enhancer-Regulated microRNAs Across 31 Human Cancers

Enhancers are cis-regulatory DNA elements that positively regulate the transcription of target genes in a tissue-specific manner, and dysregulation of target genes could lead to various diseases, such as cancer. Recent studies have shown that enhancers can regulate microRNAs (miRNAs) and participate in their biological synthesis. However, the network of enhancer-regulated miRNAs across multiple cancers is still unclear. Here, a total of 2,418 proximal enhancer-miRNA interactions and 1,280 distal enhancer-miRNA interactions were identified through the integration of genomic distance, co-expression, and 3D genome data in 31 cancers. The results showed that both proximal and distal interactions exhibited a significant cancer type-specific feature trend at the tissue level rather than at the single-cell level, and there was a noteworthy positive correlation between the expression of the miRNA and the number of enhancers regulating the same miRNA in most cancers. Furthermore, we found that there was a high correlation between the formation of enhancer-miRNA pairs and the expression of enhancer RNAs (eRNAs) whether in distal or proximal regulation. The characteristics analysis showed that miRes (enhancers that regulated miRNAs) and non-miRes presented significant differences in sequence conservation, guanine-cytosine (GC) content, and histone modification signatures. Notably, GC content, H3K4me1, and H3K36me3 were present differently between distal and proximal regulation, suggesting that they might participate in chromosome looping of enhancer-miRNA interactions. Finally, we introduced a case study, enhancer: chr1:1186391-1186507 ∼ miR-200a was highly relevant to the survival of thyroid cancer patients and a cis-eQTL SNP on the enhancer affected the expression of the TNFRSF18 gene as a tumor suppressor.

Molecular Signature of Subtypes of Non-Small-Cell Lung Cancer by Large-Scale Transcriptional Profiling: Identification of Key Modules and Genes by Weighted Gene Co-Expression Network Analysis (WGCNA)

Non-small-cell lung cancer (NSCLC) represents a heterogeneous group of malignancies consisting essentially of adenocarcinoma (ADC) and squamous cell carcinoma (SCC). Although the diagnosis and treatment of ADC and SCC have been greatly improved in recent decades, there is still an urgent need to identify accurate transcriptome profile associated with the histological subtypes of NSCLC. The present study aims to identify the key dysregulated pathways and genes involved in the development of lung ADC and SCC and to relate them with the clinical traits. The transcriptional changes between tumour and normal lung tissues were investigated by RNA-seq. Gene ontology (GO), canonical pathways analysis with the prediction of upstream regulators, and weighted gene co-expression network analysis (WGCNA) to identify co-expressed modules and hub genes were used to explore the biological functions of the identified dysregulated genes. It was indicated that specific gene signatures differed significantly between ADC and SCC related to the distinct pathways. Of identified modules, four and two modules were the most related to clinical features in ADC and SCC, respectively. CTLA4, MZB1, NIP7, and BUB1B in ADC, as well as GNG11 and CCNB2 in SCC, are novel top hub genes in modules associated with tumour size, SUVmax, and recurrence-free survival. Our research provides a more effective understanding of the importance of biological pathways and the relationships between major genes in NSCLC in the perspective of searching for new molecular targets.