Polymorphisms in the Angiogenesis-Related Genes EFNB2, MMP2 and JAG1 Are Associated with Survival of Colorectal Cancer Patients

PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

A DNA methylation signature for the prediction of tumour recurrence in stage II colorectal cancer

BACKGROUND

A major challenge in stage II colorectal carcinoma is to identify patients with increased risk of recurrence. Biomarkers that distinguish patients with poor prognosis from patients without recurrence are currently lacking. This study aims to develop a robust DNA methylation classifier that allows the prediction of recurrence and chemotherapy benefit in patients with stage II colorectal cancer. We performed a genome-wide DNA methylation capture sequencing in 243 stage II colorectal carcinoma samples and identified a relapse-specific DNA methylation signature consisting of eight CpG sites.

METHODS

Two hundred and forty-three patients with stage II CRC were enrolled in this study. In order to select differential methylation sites among recurrence and non-recurrence stage II CRC samples, DNA methylation profiles of 62 tumour samples including 31 recurrence and 31 nonrecurrence samples were analysed using the Agilent SureSelectXT Human Methyl-Seq, a comprehensive target enrichment system to analyse CpG methylation. Pyrosequencing was applied to quantify the methylation level of candidate DNA methylation sites in 243 patients. Least absolute shrinkage and selection operator (LASSO) method was employed to build the disease recurrence prediction classifier.

RESULTS

We identified a relapse-related DNA methylation signature consisting of eight CpG sites in stage II CRC by DNA methylation capture sequencing. The classifier showed significantly higher prognostic accuracy than any clinicopathological risk factors. The Kaplan-Meier survival curve showed an association of high-risk score with poor prognosis. In multivariate analysis, the signature was the most significant prognosis factor, with an HR of 2.80 (95% CI, 1.71-4.58, P < 0.001). The signature could identify patients who are suitable candidates for adjuvant chemotherapy.

CONCLUSIONS

An eight-CpG DNA methylation signature is a reliable prognostic and predictive tool for disease recurrence in patients with stage II CRC.

Classical Angiogenic Signaling Pathways and Novel Anti-Angiogenic Strategies for Colorectal Cancer

Although productive progress has been made in colorectal cancer (CRC) researchs, CRC is the second most frequent type of malignancy and the major cause of cancer-related death among gastrointestinal cancers. As angiogenesis constitutes an important point in the control of CRC progression and metastasis, understanding the key signaling pathways that regulate CRC angiogenesis is critical in elucidating ways to inhibit CRC. Herein, we comprehensively summarized the angiogenesis-related pathways of CRC, including vascular endothelial growth factor (VEGF), nuclear factor-kappa B (NF-κB), Janus kinase (JAK)/signal transducer and activator of transcription (STAT), Wingless and int-1 (Wnt), and Notch signaling pathways. We divided the factors influencing the specific pathway into promoters and inhibitors. Among these, some drugs or natural compounds that have antiangiogenic effects were emphasized. Furthermore, the interactions of these pathways in angiogenesis were discussed. The current review provides a comprehensive overview of the key signaling pathways that are involved in the angiogenesis of CRC and contributes to the new anti-angiogenic strategies for CRC.

Identification of Key Genes and Pathways in Genotoxic Stress Induced Endothelial Dysfunction: Results of Whole Transcriptome Sequencing

Atherosclerosis is a leading cause of cardiovascular morbidity and mortality worldwide. Endothelial disfunction underlying the atherogenesis can be triggered by genotoxic stress in endothelial cells. In the presented research whole transcriptome sequencing (RNA-seq) of human coronary artery (HCAEC) and internal thoracic artery (HITAEC) endothelial cells in vitro exposed to 500 ng/mL mitomycin C (treatment group) or 0.9% NaCl (control group) was performed. Resulting to bioinformatic analysis, 56 upregulated differentially expressed genes (DEGs) and 6 downregulated DEGs with absolute fold change ≥ 2 and FDR p-value < 0.05 were selected in HCAEC exposed to mitomycin C compared to the control group; in HITAEC only one upregulated DEG was found. According to Gene Ontology enrichment analysis, DEGs in HCAEC were classified into 25 functional groups of biological processes, while in HITAEC we found no statistically significant (FDR p-value < 0.05) groups. The four largest groups containing more than 50% DEGs (“signal transduction”, “response to stimulus”, “biological regulation”, and “regulation of biological process”) were identified. Finally, candidate DEGs and pathways underlying the genotoxic stress induced endothelial disfunction have been discovered that could improve our understanding of fundamental basis of atherogenesis and help to justification of genotoxic stress as a novel risk factor for atherosclerosis.

[JAG1 promotes migration, invasion, and adhesion of triple-negative breast cancer cells by promoting angiogenesis]

OBJECTIVE

To investigate the effect of JAG1 on the malignant phenotype of triple-negative breast cancer (TNBC) and its role in angiogenesis in breast cancer microenvironment.

METHODS

The expressions of Notch molecules were detected in human TNBC 231 and 231B cells using RT-qPCR. Five female nude mice were inoculated with 231 cells and another 5 with 231B cells into the mammary fat pads, and 4-6 weeks later, the tumors were collected for immunohistochemical and immunofluorescence tests. 231 cells and 231B cells were treated with recombinant JAG (rJAG) protein and DAPT, respectively, and changes in their malignant phenotypes were assessed using CCK-8 assay, Hoechst 33258 staining, wound healing assay, Transwell chamber assay and endothelial cell adhesion assay. Western blotting was used to detect the changes in the expressions of proteins related with the malignant phenotypes of 231 and 231B cells. The effects of conditioned medium (CM) derived from untreated 231 and 231 B cells, rJAG1-treated 231 cells and DAPT-treated 231B cells on proliferation and tube formation ability of cultured human umbilical vein endothelial cells (HUVECs) were evaluated using CCK-8 assay and tube-forming assay.

RESULTS

The expression of JAG1 was higher in 231B cells than in 231 cells (P < 0.05). Tumor 231B showed higher expression of VEGFA and CD31. Compared with 231-Blank group, the migration, invasion and adhesion of 231 cells in 231-rJAG1 were significantly enhanced (P < 0.05). Protein levels of Twist1 and Snail increased (P < 0.01), anti-apoptotic protein Bcl-2 increased (P < 0.05), while DAPT inhibited the related phenomena and indicators of 231B. The 231-rJAG1-CM increased the cell number and tubule number of HUVEC (P < 0.05).

CONCLUSION

JAG1 may affect the malignant phenotype of TNBC and promote angiogenesis in the tumor microenvironment.

The EPH/Ephrin System in Colorectal Cancer

The EPH/ephrin system constitutes a bidirectional signaling pathway comprised of a family of tyrosine kinase receptors in tandem with their plasma membrane-bound ligand (ephrins). Its significance in a wide variety of physiologic and pathologic processes has been recognized during the past decades. In carcinogenesis, EPH/ephrins coordinate a wide spectrum of pathologic processes, such as angiogenesis, vessel infiltration, and metastasis. Despite the recent advances in colorectal cancer (CRC) diagnosis and treatment, it remains a leading cause of death globally, accounting for 9.2% of all cancer deaths. A growing body of literature has been published lately revitalizing our scientific interest towards the role of EPH/ephrins in pathogenesis and the treatment of CRC. The aim of the present review is to present the recent CRC data which might lead to clinical practice changes in the future.

Integrative analysis and experiments to explore angiogenesis regulators correlated with poor prognosis, immune infiltration and cancer progression in lung adenocarcinoma

Angiogenesis is the process of capillary sprouting from pre-existing vessels and it plays a critical role in the carcinogenic process of lung adenocarcinoma (LUAD). However, the association of angiogenesis regulators with the prognosis and progression of LUAD needs to be further elucidated. In this study, we adopted differential expression analysis, Cox proportional hazards (PH) regression analysis and experimental validation to identify angiogenesis regulators correlated with a poor prognosis, immune infiltration and cancer progression in LUAD. These results showed that the diagnostic and prognostic models based on COL5A2 and EPHB2 served as independent biomarkers with superior predictive ability. The patients in the high-risk group exhibited a worse prognosis in the TCGA cohort (P < 0.001, HR = 1.72, 95% CI 1.28-2.30), GSE310210 cohort (P = 0.005, HR = 2.87, 95% CI 1.46-5.61), and GSE31019 cohort (P = 0.01, HR = 2.14, 95% CI 1.19-3.86) than patients in the low-risk group. The high prognostic risk patients had a higher TMB (P < 0.001); higher fractions of M0 macrophages, neutrophils, NK cells resting, and T cells CD4 memory activated (P < 0.05); and higher expression of immune checkpoints PD-1, PDL-1, PDL-2, and B7H3 (P < 0.001). Patients in the high-risk group were more sensitive to chemotherapeutic drugs and molecular targeted drugs such as cisplatin, doxorubicin, gefitinib, and bosutinib (P < 0.0001). In addition, inhibition of COL5A2 and EPHB2 effectively suppressed the proliferation and migration of LUAD cells. The current study identified angiogenesis regulators as potential biomarkers and therapeutic targets for LUAD and may help to further optimize cancer therapy.

Genome instability-related long non-coding RNA in clear renal cell carcinoma determined using computational biology

Background

There is evidence that long non-coding RNA (lncRNA) is related to genetic stability. However, the complex biological functions of these lncRNAs are unclear.

Method

TCGA - KIRC lncRNAs expression matrix and somatic mutation information data were obtained from TCGA database. “GSVA” package was applied to evaluate the genomic related pathway in each samples. GO and KEGG analysis were performed to show the biological function of lncRNAs-mRNAs. “Survival” package was applied to determine the prognostic significance of lncRNAs. Multivariate Cox proportional hazard regression analysis was applied to conduct lncRNA prognosis model.

Results

In the present study, we applied computational biology to identify genome-related long noncoding RNA and identified 26 novel genomic instability-associated lncRNAs in clear cell renal cell carcinoma. We identified a genome instability-derived six lncRNA-based gene signature that significantly divided clear renal cell samples into high- and low-risk groups. We validated it in test cohorts. To further elucidate the role of the six lncRNAs in the model’s genome stability, we performed a gene set variation analysis (GSVA) on the matrix. We performed Pearson correlation analysis between the GSVA scores of genomic stability-related pathways and lncRNA. It was determined that LINC00460 and LINC01234 could be used as critical factors in this study. They may influence the genome stability of clear cell carcinoma by participating in mediating critical targets in the base excision repair pathway, the DNA replication pathway, homologous recombination, mismatch repair pathway, and the P53 signaling pathway.

Conclusion subsections

These data suggest that LINC00460 and LINC01234 are crucial for the stability of the clear cell renal cell carcinoma genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08356-9.

ATP5A1 Participates in Transcriptional and Posttranscriptional Regulation of Cancer-Associated Genes by Modulating Their Expression and Alternative Splicing Profiles in HeLa Cells

Background: Aberrant expression and alternative splicing of oncogenes are the driving events in tumor initiation and development. But how these events are regulated in cancer cells is largely unknown. Functions of ATP5A1, an important mitochondrial ATP synthase gene, in transcriptional and posttranscriptional regulation were explored in this study. Methods: ATP5A1 was overexpressed using plasmid-transformed HeLa cells, and its influence on cell apoptosis and proliferation is evaluated. Transcriptome sequencing was then performed using RNA-seq to study the changes in gene expression and regulation of alternative splicing events. Validation of the implicated genes was achieved using RT-qPCR analysis. Results: It was found that ATP5A1 could significantly promote cellular apoptosis, but it had no influence on cell proliferation. ATP5A1 overexpression significantly increased the expression levels of genes associated with the innate immune response, angiogenesis, and collagen catabolic processes. This included enrichment of MMP2 and MMP19. It was also found that ATP5A1 could interfere with the alternative splicing of hundreds of genes associated with glucose homeostasis, HIF-1 signaling activation, and several pathways associated with cancers. Eight ATP5A1-regulated differentially expressed genes and 3 genes altered by splicing were selected and validated using RT-qPCR analysis. Conclusions: In summary, we illustrate the regulatory functions of ATP5A1 on the transcriptome of HeLa cells by exploring its influence on gene expression and alternative splicing. The results suggest that ATP5A1 may play an important regulatory role in cervical cancer cells by regulating expression and alternative splicing of cancer-associated genes. This study provides novel insights into the current understanding of the mechanisms of ATP5A1 on carcinogenesis and cancer progression.

Genotype-Based Gene Expression in Colon Tissue-Prediction Accuracy and Relationship with the Prognosis of Colorectal Cancer Patients

Colorectal cancer (CRC) survival has environmental and inherited components. The expression of specific genes can be inferred based on individual genotypes-so called expression quantitative trait loci. In this study, we used the PrediXcan method to predict gene expression in normal colon tissue using individual genotype data from 91 CRC patients and examined the correlation ρ between predicted and measured gene expression levels. Out of 5434 predicted genes, 58% showed a negative ρ value and only 16% presented a ρ higher than 0.10. We subsequently investigated the association between genotype-based gene expression in colon tissue for genes with ρ > 0.10 and survival of 4436 CRC patients. We identified an inverse association between the predicted expression of ARID3B and CRC-specific survival for patients with a body mass index greater than or equal to 30 kg/m2 (HR (hazard ratio) = 0.66 for an expression higher vs. lower than the median, p = 0.005). This association was validated using genotype and clinical data from the UK Biobank (HR = 0.74, p = 0.04). In addition to the identification of ARID3B expression in normal colon tissue as a candidate prognostic biomarker for obese CRC patients, our study illustrates the challenges of genotype-based prediction of gene expression, and the advantage of reassessing the prediction accuracy in a subset of the study population using measured gene expression data.