Triple-layer dissection of the lung adenocarcinoma transcriptome: regulation at the gene, transcript, and exon levels

TEDC2 plays an oncogenic role and serves as a therapeutic target of hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common malignancies and tends to have a poor prognosis due to its insidious onset, difficulty in early diagnosis, and limited treatment options. Tubulin epsilon and delta complex 2 (TEDC2), also known as C16orf59, is implicated in maintaining centriole stability, but the involvement of TEDC2 in HCC remains unknown. This study aimed to investigate the expression profile and potential mechanisms of TEDC2 in HCC.

METHODS

Multiple RNA sequencing datasets were screened for differentially expressed genes in HCC, and the prognosis-related gene, TEDC2, was further screened as a target gene in this study. The expression of TEDC2 in public datasets and clinical specimens was analyzed, and the involvement of TEDC2 in HCC was investigated by bioinformatic analysis and in vitro experiments.

RESULTS

TEDC2 levels were elevated in HCC compared to healthy livers. Overexpression of TEDC2 was positively correlated with pathologic stage and histologic grade. In addition, TEDC2 was found to be an independent prognostic predictor. An excellent prognostic model of HCC was successfully constructed with TEDC2 in combination with the TNM stage. Bioinformatic analysis revealed that overexpression of TEDC2 might be associated with impaired tumor immunity in HCC, as evidenced by increased infiltration of T helper 2 (Th2) cells and reduced infiltration of cytotoxic cells. Further studies showed that TP53 mutations regulated TEDC2 expression, and TEDC2 was significantly associated with drug sensitivity. Moreover, overexpression of TEDC2 promoted cell metastasis and proliferation in vitro.

CONCLUSION

These findings initially suggested a crucial effect of TEDC2 overexpression on HCC tumor progression, suggesting its potential as a novel prognostic and therapeutic target in HCC.

A pan-cancer analysis of the prognostic implication and oncogenic role of tubulin epsilon and delta complex 2 (TEDC2) in human tumors

Introduction

Tubulin epsilon and delta complex 2 (TEDC2) is widely expressed in various human tissues and primarily governs centriole stability. However, the biological significance of TEDC2 in pan-cancer is unclear.

Methods

In this study, we employed R software and various online bioinformatics analysis tools to investigate the functional attributes of TEDC2 in human tumours and its potential involvement in immune response. The status of TEDC2 expression was evaluated in samples from the TCGA and GEO datasets, as well as in tumour and corresponding normal samples from the TCGA database. Subsequently, Kaplan-Meier estimates, clinical correlations, and univariate Cox regressions were used to analyze the 33 types of tumors from TCGA and determine the prognostic significance of TEDC2. Moreover, nomogram models were formulated using three distinct tumours, namely kidney renal clear cell carcinoma (KIRC), lung adenocarcinoma (LUAD), and liver hepatocellular carcinoma (LIHC), to evaluate the prognostic significance of TEDC2 in tumours. Furthermore, TEDC2 was investigated for its correlation with the levels of immune cell infiltration, and a functional enrichment analysis was conducted to identify potential signalling pathways involving TEDC2.

Results

Differential analysis revealed that 16 tumour types expressed TEDC2 to a greater extent than normal tissues. The abnormal expression of TEDC2 can predict survival outcomes in patients with adrenocortical carcinoma (ACC), KIRC, kidney renal papillary cell carcinoma (KIRP), LUAD, LIHC, lower grade glioma (LGG), and thymoma (THYM). Subsequent results indicated that TEDC2 has the ability to influence ECM regulators, cell cycle, and Immune checkpoint-associated signalling pathways, which could potentially lead to a poor prognosis and tumour progression.

Discussion

TEDC2 has been identified as a potential therapeutic target that could predict the prognosis of multiple tumour types, making it a promising target for reversing tumour development.

TEDC2 correlated with prognosis and immune microenvironment in lung adenocarcinoma

Tubulin epsilon and delta complex 2 (TEDC2) is a protein coding gene whose functions are poorly identified yet. This study aimed to identify the role of TEDC2 in prognosis and immune microenvironment of lung adenocarcinoma (LUAD). Through The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, the mRNA expression of TEDC2 was upregulated in LUAD tissues compared to normal tissues. The protein level of TEDC2 was also higher in LUAD in the Human Protein Atlas. The receiver operating characteristic (ROC) curve showed that high TEDC2 level could distinguish LUAD patients from normal subjects. In addition, the impact of TEDC2 expression on prognosis was evaluated by Kaplan-Meier and Cox regression analyses, and the results suggested that high TEDC2 expression was significantly associated with poor prognosis and was the independent prognostic factor in LUAD. GO and KEGG pathway analyses indicated the co-expressed genes of TEDC2 were mainly related to mitotic cell cycle processes. Importantly, high expression of TEDC2 indicated low infiltration of immune cells, especially dendritic cells and B cells. TEDC2 was also positively correlated with immune checkpoints such as PDCD1, LAG3 and CD276. Taken together, this study preliminarily revealed the clinical significance of TEDC2 in LUAD and provided novel insights into the role of TEDC2 in immune microenvironment.

Identification of Novel Prognostic Markers Associated With Laryngeal Squamous Cell Carcinoma Using Comprehensive Analysis

Background

Laryngeal squamous cell carcinoma (LSCC) is a leading malignant cancer of the head and neck. Patients with LSCC, in which the cancer has infiltrated and metastasized, have a poor prognosis. Therefore, there is an urgent need to identify more potential targets for drugs and biomarkers for early diagnosis.

Methods

RNA sequence data from LSCC and patients’ clinical traits were obtained from the Gene Expression Omnibus (GEO) (GSE142083) and The Cancer Genome Atlas (TCGA) database. Differentially expressed gene (DEG) analysis and weighted gene co-expression network analysis (WGCNA) were performed to identify hub genes. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, prognostic value analysis, receiver operating characteristic (ROC) curve analysis, gene mutation analysis, tumor-infiltrating immune cell abundance profile estimation, gene set variation analysis (GSVA), and gene set enrichment analysis (GSEA) were performed. Single-gene RNA sequencing data were obtained from the GSE150321 dataset. Cell proliferation and viability were confirmed by the CCK-8 assay and real-time PCR.

Results

A total of 701 DEGs, including 329 upregulated and 372 downregulated genes, were screened in the GSE142083 dataset. Using WGCNA, three modules were identified to be closely related to LSCC. After intersecting the DEGs and performing univariate and multivariate Cox analyses, a novel prognostic model based on three genes (SLC35C1, HOXB7, and TEDC2) for LSCC was established. Interfering TEDC2 expression inhibited tumor cell proliferation and migration.

Conclusions

Our results show that SLC35C1, HOXB7, and TEDC2 have the potential to become new therapeutic targets and prognostic biomarkers for LSCC.

Key Candidate Prognostic Biomarkers Correlated with Immune Infiltration in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common subtype of primary liver cancer, which causes ~800,000 deaths annually world-wide. Immune checkpoint inhibitor (ICI) has reformed cancer therapy and achieved unprecedented results in various malignancies, including HCC. However, the response rate of immunotherapy is very low in HCC. Considereing the complicated and unique immune status in liver, we hypothesize that critical molecules will affect prognosis and correlate with immune context in the tumor microenvironment of HCC.

Methods

Using Kaplan–Meier plotter, GEPIA2 and Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB), survival genes and their prognostic value were estimated in HCC. Based on Tumor Immune Estimation Resource (TIMER), association between survival genes and immune infiltration was examined in HCC. FunRich and STRING were used to analyze gene ontology and protein–protein interaction (PPI) Network, qRT-PCR was used to measure mRNA level of candidates; and a Cell Counting Kit-8 was used to measure proliferation of HCC cell line.

Results

Using multiple databases, we identified 36 key prognostic genes highly expressed in HCC and associated with poor survival of patients. Meanwhile, the 36 gene signatures correlated with immune infiltration in HCC. Moreover, these genes were significantly associated with exhausted T cells and polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) in HCC. Among the 36 key genes, SKA3, SGOL2, SPINDOC, TEDC2, TMCO3 and NUP205 were highly expressed in tumor samples compared with adjacent normal tissues in our HCC cohort (n=22). Additionally, proliferation of SMMC7721 cell line was inhibited when it interfered with SiRNA of each gene.

Conclusion

The 36 genes may serve as potential prognostic biomarkers and molecular targets to ameliorate tumor immune microenvironment (TIME) in HCC and therefore represent a novel avenue for individualized immunotherapy in HCC.

The impact of single-nucleotide polymorphisms on intravesical recurrence after bacillus Calmette-Guérin therapy for non-muscle invasive bladder cancer in a genome-wide association study

OBJECTIVE

Bacillus Calmette-Guérin (BCG) instillation therapy is widely used to reduce intravesical recurrence in non-muscle invasive bladder cancer (NMIBC). In this study, we aimed to reveal the genetic variations associated with intravesical recurrence after BCG therapy for NMIBC in a genome-wide association study (GWAS).

MATERIALS AND METHODS

This study included Japanese patients with NMIBC, in whom genomic DNA was obtained from whole blood samples. The association between genetic variation and treatment failure was analyzed by GWAS in 44 patients treated with BCG instillation as a discovery cohort. Candidate single-nucleotide polymorphisms (SNPs) were examined separately in 47 patients treated with BCG instillation and in 62 patients treated with chemotherapeutic agent instillation as validation studies.

RESULTS

Among the 44 patients in the discovery cohort, 14 cases experienced intravesical recurrent diseases. GWAS identified 12 candidate SNPs (rs9374832, rs35176001, rs363765, rs2127120, rs4277759, rs73664140, rs1607282, rs12141654, rs4541358, rs6986852, rs12373386, and rs17637903). In the validation study, a genetic risk stratification model using the number of risk alleles in rs363765 and rs6986852 discriminated the risk of intravesical recurrence after BCG therapy, but not after non-BCG therapy.

CONCLUSION

This study suggested that several SNPs were associated with intravesical recurrence after BCG therapy for NMIBC. A genetic risk model may be useful to predict intravesical recurrence after BCG therapy, warranting further research and development for clinical application.

Genome-wide associations between alcohol consumption and blood DNA methylation: evidence from twin study

Aim: Alcohol intake alters DNA methylation profiles and methylation might mediate the association between alcohol and disease, but limited number of positive CpG sites repeatedly replicated. Materials & methods: In total, 57 monozygotic (MZ) twin pairs discordant for alcohol drinking from the Chinese National Twin Registry and 158 MZ and dizygotic twin pairs in the Swedish Adoption/Twin Study of Aging were evaluated. DNA methylation was detected using the Infinium HumanMethylation450 BeadChip. Results: Among candidate CpG sites, cg07326074 was significantly correlated with drinking after adjusting for covariates in MZ twins in both datasets but not in the entire sample or dizygotic twins. Conclusion: The hypermethylation of cg07326074, located in the tumor-promoting gene C16orf59, was associated with alcohol consumption.

Identification of Key Genes With Differential Correlations in Lung Adenocarcinoma

Background

With the advent of large-scale molecular profiling, an increasing number of oncogenic drivers contributing to precise medicine and reshaping classification of lung adenocarcinoma (LUAD) have been identified. However, only a minority of patients archived improved outcome under current standard therapies because of the dynamic mutational spectrum, which required expanding susceptible gene libraries. Accumulating evidence has witnessed that understanding gene regulatory networks as well as their changing processes was helpful in identifying core genes which acted as master regulators during carcinogenesis. The present study aimed at identifying key genes with differential correlations between normal and tumor status.

Methods

Weighted gene co-expression network analysis (WGCNA) was employed to build a gene interaction network using the expression profile of LUAD from The Cancer Genome Atlas (TCGA). R package DiffCorr was implemented for the identification of differential correlations between tumor and adjacent normal tissues. STRING and Cytoscape were used for the construction and visualization of biological networks.

Results

A total of 176 modules were detected in the network, among which yellow and medium orchid modules showed the most significant associations with LUAD. Then genes in these two modules were further chosen to evaluate their differential correlations. Finally, dozens of novel genes with opposite correlations including ATP13A4-AS1, HIGD1B, DAP3, and ISG20L2 were identified. Further biological and survival analyses highlighted their potential values in the diagnosis and treatment of LUAD. Moreover, real-time qPCR confirmed the expression patterns of ATP13A4-AS1, HIGD1B, DAP3, and ISG20L2 in LUAD tissues and cell lines.

Conclusion

Our study provided new insights into the gene regulatory mechanisms during transition from normal to tumor, pioneering a network-based algorithm in the application of tumor etiology.

Genome wide methylation profiling of selected matched soft tissue sarcomas identifies methylation changes in metastatic and recurrent disease

In this study we used the Illumina Infinium Methylation array to investigate in a cohort of matched archival human tissue samples (n = 32) from 14 individuals with soft tissue sarcomas if genome-wide methylation changes occur during metastatic and recurrent (Met/Rec) disease. A range of sarcoma types were selected for this study: leiomyosarcoma (LMS), myxofibrosarcoma (MFS), rhabdomyosarcoma (RMS) and synovial sarcoma (SS). We identified differential methylation in all Met/Rec matched samples, demonstrating that epigenomic differences develop during the clonal evolution of sarcomas. Differentially methylated regions and genes were detected, not been previously implicated in sarcoma progression, including at PTPRN2 and DAXX in LMS, WT1-AS and TNXB in SS, VENTX and NTRK3 in pleomorphic RMS and MEST and the C14MC / miR-379/miR-656 in MFS. Our overall findings indicate the presence of objective epigenetic differences across primary and Met/Rec human tissue samples not previously reported.

Molecular Pathology of Primary Non-small Cell Lung Cancer

Lung carcinoma is one of the most common human cancers and is estimated to have an incidence of approximately 2 million new cases per year worldwide with a 20% mortality rate. Lung cancer represents one of the leading causes of cancer related death in the world. Of all cancer types to affect the pulmonary system, non-small cell lung carcinoma comprises approximately 80-85% of all tumors. In the past few decades cytogenetic and advanced molecular techniques have helped define the genomic landscape of lung cancer, and in the process, revolutionized the clinical management and treatment of patients with advanced non-small cell lung cancer. The discovery of specific, recurrent genetic abnormalities has led to the development of targeted therapies that have extended the life expectancy of patients who develop carcinoma of the lungs. Patients are now routinely treated with targeted therapies based on identifiable molecular alterations or other predictive biomarkers which has led to a revolution in the field of pulmonary pathology and oncology. Numerous different testing modalities, with various strengths and limitations now exist which complicate diagnostic algorithms, however recently emerging consensus guidelines and recommendations have begun to standardize the way to approach diagnostic testing of lung carcinoma. Herein we provide an overview of the molecular genetic landscape of non-small cell lung carcinoma, with attention to those clinically relevant alterations which drive management, as well as review current recommendations for molecular testing.

A functional variant in TNXB promoter associates with the risk of esophageal squamous-cell carcinoma

Our previous study identified a tag single-nucleotide polymorphism (SNP) rs204900 in TNXB associated with risk of esophageal squamous-cell carcinoma (ESCC) in the Chinese population. However, the functional role of TNXB and causal variants had not been interrogated in that study. In the present study, we explored the effects of TNXB expression in the development of ESCC and searched for functional variants in this gene. We found TNXB was downregulated in ESCC tumors. Using small interfering RNAs and CRISPR-Cas9 methods, we identified that both knockdown and knockout of TNXB significantly promoted ESCC cell growth in vitro, suggesting a tumor suppressor role of this gene in ESCC. Through further fine-mapping analysis, we identified that a noncoding variant in the promoter of TNXB, rs411337, predisposed to ESCC risk (odds ratio = 1.36, 95% confidence interval: 1.22-1.51, P = 9.10 × 10^-9 ). These findings revealed the functional mechanism of TNXB in the development of ESCC and may contribute to the prevention and treatment of this disease in the future.

Integrative analysis with expanded DNA methylation data reveals common key regulators and pathways in cancers

The integration of genomic and DNA methylation data has been demonstrated as a powerful strategy in understanding cancer mechanisms and identifying therapeutic targets. The TCGA consortium has mapped DNA methylation in thousands of cancer samples using Illumina Infinium Human Methylation 450 K BeadChip (Illumina 450 K array) that only covers about 1.5% of CpGs in the human genome. Therefore, increasing the coverage of the DNA methylome would significantly leverage the usage of the TCGA data. Here, we present a new model called EAGLING that can expand the Illumina 450 K array data 18 times to cover about 30% of the CpGs in the human genome. We applied it to analyze 13 cancers in TCGA. By integrating the expanded methylation, gene expression, and somatic mutation data, we identified the genes showing differential patterns in each of the 13 cancers. Many of the triple-evidenced genes identified in majority of the cancers are biomarkers or potential biomarkers. Pan-cancer analysis also revealed the pathways in which the triple-evidenced genes are enriched, which include well known ones as well as new ones, such as axonal guidance signaling pathway and pathways related to inflammatory processing or inflammation response. Triple-evidenced genes, particularly TNXB, RRM2, CELSR3, SLC16A3, FANCI, MMP9, MMP11, SIK1, and TRIM59 showed superior predictive power in both tumor diagnosis and prognosis. These results have demonstrated that the integrative analysis using the expanded methylation data is powerful in identifying critical genes/pathways that may serve as new therapeutic targets.

Identification of DNA methylation associated gene signatures in endometrial cancer via integrated analysis of DNA methylation and gene expression systematically

Objective

Endometrial cancer (EC) is a common gynecologic cancer worldwide. However, the pathogenesis of EC has not been epigenetically elucidated. Here, this study aims to describe the DNA methylation profile and identify favorable gene signatures highly associated with aberrant DNA methylation changes in EC.

Methods

The data regarding DNA methylation and gene expression were downloaded from The Cancer Genome Atlas (TCGA) database. Differentially methylated CpG sites (DMCs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs) were identified, and the relationship between the 2 omics was further analyzed. In addition, weighted CpG site co-methylation network (WCCN) was constructed followed by an integrated analysis of DNA methylation and gene expression data.

Results

Four hundred thirty-one tumor tissues and 46 tissues adjacent tumor of EC patients were analyzed. One thousand one hundred thirty-five DMCs (merging to 10 DMRs), and 1,488 DEGs were obtained between tumor and normal groups, respectively. One hundred forty-eight DMCs-DEGs correlated pairs and 13 regional DMCs-DEGs pairs were obtained. Interestingly, we found that some hub genes in 2 modules among 8 modules of WCCN analysis were down-regulated in tumor samples. Furthermore, protocadherins (PCDHs) clusters, DDP6, TNXB, and ZNF154 were identified as novel deregulated genes with altered methylation in EC.

Conclusion

Based on the analysis of DNA methylation in a systematic view, the potential long-range epigenetic silencing (LRES) composed of PCDHs was reported in ECs for the first time. PCDHs clusters, DDP6, and TNXB were firstly found to be associated with tumorigenesis, and may be novel candidate biomarkers for EC.

Global DNA methylation patterns in Barrett's esophagus, dysplastic Barrett's, and esophageal adenocarcinoma are associated with BMI, gender, and tobacco use

BACKGROUND

The risk of developing Barrett's esophagus (BE) and/or esophageal adenocarcinoma (EAC) is associated with specific demographic and behavioral factors, including gender, obesity/elevated body mass index (BMI), and tobacco use. Alterations in DNA methylation, an epigenetic modification that can affect gene expression and that can be influenced by environmental factors, is frequently present in both BE and EAC and is believed to play a role in the formation of BE and its progression to EAC. It is currently unknown whether obesity or tobacco smoking influences the risk of developing BE/EAC via the induction of alterations in DNA methylation. To investigate this possibility, we assessed the genome-wide methylation status of 81 esophageal tissues, including BE, dysplastic BE, and EAC epithelia using HumanMethylation450 BeadChips (Illumina).

RESULTS

We found numerous differentially methylated loci in the esophagus tissues when comparing males to females, obese to lean individuals, and smokers to nonsmokers. Differences in DNA methylation between these groups were seen in a variety of functional genomic regions and both within and outside of CpG islands. Several cancer-related pathways were found to have differentially methylated genes between these comparison groups.

CONCLUSIONS

Our findings suggest obesity and tobacco smoking may influence DNA methylation in the esophagus and raise the possibility that these risk factors affect the development of BE, dysplastic BE, and EAC through influencing the epigenetic status of specific loci that have a biologically plausible role in cancer formation.