Analysis of the interplay between methylation and expression reveals its potential role in cancer aetiology

Exon 1 methylation status of CDH13 is associated with decreased overall survival and distant metastasis in patients with postoperative colorectal cancer

BACKGROUND

Cadherin 13 (CDH13) is a member of the cadherin superfamily that exerts tumor-suppressive effects on cancers derived from epithelial cells. Although hypermethylation of CDH13 promoter has been reported in various cancers, its prognostic value for colorectal cancer (CRC) is still controversial. The methylation alterations of CDH13 within exon 1 have not yet been investigated.

METHODS

A total of 49 CRC patients were recruited for the prospective study. The methylation status of CpG sites was quantified by Bisulfite Amplicon Sequencing (BSAS) in malignant tissues and adjacent normal tissues. The primary endpoint of the study was overall survival (OS) after surgery. The relationship between methylation level with pathological stage and OS was also evaluated.

RESULTS

Compared with adjacent normal tissues, the overall average methylation level within exon 1 was significantly increased in tumor tissues (p < 0.001). The association study showed that the hypermethylation status of the CpG1 site was non-significantly associated with the presence of distant metastasis (p = 0.032). Moreover, the hypermethylation of two CpG sites, including CpG1 (p = 0.003) and CpG5 (p = 0.032), was associated with worse OS in CRC. Co-hypermethylation of CpG1 and CpG5 sites was significantly associated with a worse clinical outcome (HR: 4.43 [95% CI 1.27-15.46]; p = 0.019) in multivariate Cox regression analysis.

CONCLUSION

The methylation level of CDH13 exon 1 in CRC tissue was significantly higher than in adjacent normal tissues. Hypermethylation at the CpG1 site suggests a risk of distant metastasis in CRC. The hypermethylation of the CpG1 site and CpG5 site, including the co-hypermethylation of these two sites, may serve as a valuable prognostic biomarker.

Analysis of the Gene Networks and Pathways Correlated with Tissue Differentiation in Prostate Cancer

Prostate cancer (PCa) is the most prevalent non-cutaneous cancer in men. Early PCa detection has been made possible by the adoption of screening methods based on the serum prostate-specific antigen and Gleason score (GS). The aim of this study was to correlate gene expression with the differentiation level of prostate adenocarcinomas, as indicated by GS. We used data from The Cancer Genome Atlas (TCGA) and included 497 prostate cancer patients, 52 of which also had normal tissue sample sequencing data. Gene ontology analysis revealed that higher GSs were associated with greater responses to DNA damage, telomere lengthening, and cell division. Positive correlation was found with transcription factor activator of the adenovirus gene E2 (E2F) and avian myelocytomatosis viral homolog (MYC) targets, G2M checkpoints, DNA repair, and mitotic spindles. Immune cell deconvolution revealed high M0 macrophage counts and an increase in M2 macrophages dependent on the GS. The molecular pathways most correlated with GSs were cell cycle, RNA transport, and calcium signaling (depleted). A combinatorial approach identified a set of eight genes able to differentiate by k-Nearest Neighbors (kNN) between normal tissues, low-Gleason tissues, and high-Gleason tissues with high accuracy. In conclusion, our study could be a step forward to better understanding the link between gene expression and PCa progression and aggressiveness.

Comprehensive multiomics and in silico approach uncovers prognostic, immunological, and therapeutic roles of ANLN in lung adenocarcinoma

The anillin actin-binding protein (ANLN) is immensely overexpressed in cancers, including lung cancer (LC). Phytocompounds have gained interest due to their broader potential and reduced unwanted effects. Screening numerous compounds presents a challenge, but in silico molecular docking is pragmatic. The present study aims to identify the role of ANLN in lung adenocarcinoma (LUAD), along with identification and interaction analysis of anticancer and ANLN inhibitory phytocompounds followed by molecular dynamics (MD) simulation. Using a systematic approach, we found that ANLN is significantly overexpressed in LUAD and mutated with a frequency of 3.73%. It is linked with advanced stages, clinicopathological parameters, worsening of relapse-free survival (RFS), and overall survival (OS), pinpointing its oncogenic and prognostic potential. High-throughput screening and molecular docking of phytocompounds revealed that kaempferol (flavonoid aglycone) interacts strongly with the active site of ANLN protein via hydrogen bonds, Vander Waals interactions, and acts as a potent inhibitor. Furthermore, we discovered that ANLN expression was found to be significantly higher (p) in LC cells compared to normal cells. This is a propitious and first study to demonstrate ANLN and kaempferol interactions, which might eventually lead to removal of rout from cell cycle regulation posed by ANLN overexpression and allow it to resume normal processes of proliferation. Overall, this approach suggested a plausible biomarker role of ANLN and the combination of molecular docking subsequently led to the identification of contemporary phytocompounds, bearing symbolic anticancer effects. The findings would be advantageous for pharmaceutics but require validation using in vitro and in vivo methods. HIGHLIGHTS: • ANLN is significantly overexpressed in LUAD. • ANLN is implicated in the infiltration of TAMs and altering plasticity of TME. • Kaempferol (potential ANLN inhibitor) shows important interactions with ANLN which could remove the alterations in cell cycle regulation, imposed by ANLN overexpression eventually leading to normal process of cell proliferation.

Clinical significance and functional role of adhesion G-protein-coupled receptors in human pancreatic ductal adenocarcinoma

BACKGROUND

The adhesion G-protein-coupled receptors (GPCRs) play crucial roles in tumour pathogenesis, however, their clinical significance in pancreatic ductal adenocarcinoma (PDAC) remains unclear.

METHODS

We analysed 796 PDAC patients, including 331 from public data sets (TCGA, ICGC and GSE57495) and 465 from independent cohorts (training: n = 321, validation: n = 144). Using in-vitro studies, we confirmed the biological function of the candidate GPCRs.

RESULTS

Analysis of all 33 adhesion GPCRs, led to identify GPR115, as the only significant prognostic factor in all public data sets. The patients with high GPR115 expression exhibited significantly poorer prognosis for OS and RFS, in training (P < 0.01, P < 0.01) and validation cohort (P < 0.01, P = 0.04). Multivariate analysis indicated that GPR115 high expression was an independent prognostic factor in both cohorts (HR = 1.43; P = 0.01, HR = 2.55; P < 0.01). A risk-prediction model using Cox regression by incorporating GPR115 and clinicopathological factors accurately predicted 5-year survival following surgery. In addition, GPR115 silencing inhibited cell proliferation and migration in PDAC cells.

CONCLUSION

We demonstrated that GPR115 has important prognostic significance and functional role in tumour progression; providing a rationale that this may be a potential therapeutic target in patients with PDAC.

The impact of Piscirickettsia salmonis infection on genome-wide DNA methylation profile in Atlantic Salmon

Salmon rickettsial septicaemia (SRS), caused by the bacteria Piscirickettsia salmonis (P. salmonis), is responsible for significant mortality in farmed Atlantic salmon in Chile. Currently there are no effective treatments or preventive measures for this disease, although genetic selection or genome engineering to increase salmon resistance to SRS are promising strategies. The accuracy and efficiency of these strategies are usually influenced by the available biological background knowledge of the disease. The aim of this study was to investigate DNA methylation changes in response to P. salmonis infection in the head kidney and liver tissue of Atlantic salmon, and the interaction between gene expression and DNA methylation in the same tissues. The head kidney and liver methylomes of 66 juvenile salmon were profiled using reduced representation bisulphite sequencing (RRBS), and compared between P. salmonis infected animals (3 and 9 days post infection) and uninfected controls, and between SRS resistant and susceptible fish. Methylation was correlated with matching RNA-Seq data from the same animals, revealing that methylation in the first exon leads to an important repression of gene expression. Head kidney methylation showed a clear response to the infection, associated with immunological processes such as actin cytoskeleton regulation, phagocytosis, endocytosis and pathogen associated pattern receptor signaling. Our results contribute to the growing understanding of the role of methylation in regulation of gene expression and response to infectious diseases and could inform the incorporation of epigenetic markers into genomic selection for disease resistant and the design of diagnostic epigenetic markers to better manage fish health in salmon aquaculture.

Correlation Between the RNA Methylation Genes and Immune Infiltration and Prognosis of Patients with Hepatocellular Carcinoma: A Pan-Cancer Analysis

Background

RNA methylation is one of the most common RNA modifications and is dynamic and reversible. The enzymes and downstream effectors associated with RNA methylation modifications can be targeted to regulate RNA methylation levels. This mechanism can affect RNA processing, metabolism, cell proliferation and migration, and regulation of physiological or pathological processes. The aim of this study was to investigate the role of RNA methylation-related genes in hepatocellular carcinoma (HCC).

Methods

Baseline RNA methylation data were extracted from The Cancer Genome Atlas database. The expression pattern, predictive value, mutational profile, and interaction network of RNA methylation genes in pancancer were examined. Then, the association between the expression of RNA methylation genes and immune infiltration was investigated. In addition, a risk score model for HCC was developed and analyzed.

Results

Cancer cells had a higher expression of RNA methylation genes than normal cells in some cancer cells, and a higher expression of RNA methylation genes could negatively affect patient prognosis. Enrichment analysis revealed that RNA methylation genes are involved in the mRNA surveillance pathway and RNA degradation and transport. A 4-gene (ALYREF, NSUN4, TRMT6, YTHDF1) prognostic signature was established to predict HCC prognosis based on RNA methylation-related genes. Finally, the role of prognostic models in HCC was validated.

Conclusion

RNA methylation genes can be an indicator of oncogenicity in relation to HCC prognosis and are associated with immune infiltration in the tumour microenvironment. This finding could provide clinicians with the opportunity to explore new strategic approaches.

The effect of key DNA methylation in different regions on gene expression in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a common cancer with high morbidity and mortality. As we all know, the alteration of DNA methylation has a crucial impact on the occurrence of HCC. However, the mechanism of the effect of DNA methylation in different regions on gene expression is still unclear. Here, by computing and analyzing the distribution of differential methylation in 12 different regions in HCC tissues and adjacent normal tissues, not only the hypermethylation of CpG islands and global hypomethylation were found, but also a stable distribution pattern of differential methylation in HCC was found. Then the correlations between DNA methylations in different regions and gene expressions were calculated, and the diversity of correlations in different regions was determined. The key genes of differential methylation and differential expression related to the survival of HCC patients were obtained by using Cox regression analysis, a four-gene prognostic risk scoring model was constructed, and the prognostic performance was well verified. The regions of the differentially methylated CpG sites corresponding to the four key genes were located and their influences on the expression were analyzed. The results indicate that the promoter, first exon, 5'UTR, sixth exon, N_Shore, and S_Shore hypomethylation promotes the expression of key oncogenes, which together lead to the occurrence of HCC. These results might help to study the role of DNA methylation in HCC and provide potential biomarkers for the diagnosis of HCC.

Identifying molecular targets for reverse aging using integrated network analysis of transcriptomic and epigenomic changes during aging

Aging is associated with widespread physiological changes, including skeletal muscle weakening, neuron system degeneration, hair loss, and skin wrinkling. Previous studies have identified numerous molecular biomarkers involved in these changes, but their regulatory mechanisms and functional repercussions remain elusive. In this study, we conducted next-generation sequencing of DNA methylation and RNA sequencing of blood samples from 51 healthy adults between 20 and 74 years of age and identified aging-related epigenetic and transcriptomic biomarkers. We also identified candidate molecular targets that can reversely regulate the transcriptomic biomarkers of aging by reconstructing a gene regulatory network model and performing signal flow analysis. For validation, we screened public experimental data including gene expression profiles in response to thousands of chemical perturbagens. Despite insufficient data on the binding targets of perturbagens and their modes of action, curcumin, which reversely regulated the biomarkers in the experimental dataset, was found to bind and inhibit JUN, which was identified as a candidate target via signal flow analysis. Collectively, our results demonstrate the utility of a network model for integrative analysis of omics data, which can help elucidate inter-omics regulatory mechanisms and develop therapeutic strategies against aging.

GPR115 Contributes to Lung Adenocarcinoma Metastasis Associated With LAMC2 and Predicts a Poor Prognosis

GPR115, a member of the adhesion G protein-coupled receptor family, is dysregulated in many cancers. However, the expression and function of GRP115 in non-small cell lung cancer (NSCLC) is not clear. Here, we examined the expression pattern, clinical significance, and function of GPR115 in NSCLC by analysis of clinical specimens and human cell lines and bioinformatics analysis. Immunohistochemical analysis of clinical samples showed that GPR115 was significantly upregulated in NSCLC tissues compares with normal lung epithelial tissue (P < 0.05). And GPR115 overexpression is an independent prognostic factor for 5-year overall survival of NSCLC patients [hazard ratio (HR)=1.625, P = 0.008]. Interestingly, higher expression of GPR115 was strongly correlation with differentiation level (P = 0.027), tumor size (P = 0.010), lymph node metastasis (P = 0.022), tumor-node-metastasis stage (P = 0.008), and poor prognosis of lung adenocarcinoma (LUAD, all P = 0.039), but not lung squamous cell carcinoma (LUSC, P > 0.05). Moreover, downregulation of GPR115 by RNA interference in human lung cancer lines inhibited cell proliferation, migration, and invasion. Preliminary bioinformatic analysis confirmed that GPR115 was closely associated with LAMC2 (Spearman correlation coefficient=0.67, P < 0.05), which was accumulated in ECM-receptor interaction and focal adhesion. Consistent with these findings, deceased of GPR115 was associated with E-cadherin, N-cadherin and Vimentin confirmed by western blot. In conclusion, these data suggest that GPR115 may play a role in the tumor growth and metastasis and may have utility as a diagnostic and prognostic marker for LUAD, but not LUSC.

5'UTR methylation in different genes from workers exposed to volatile organic compounds: A new insight for considering an epigenetic mark as a functional correlate

Gene-specific methylation has been related with transcriptional/translational consequences in different cells; also, this epigenetic modification is affected by environmental exposures. In previous studies, CYP2E1 activity in toluene-exposed workers was decreased compared to controls, however, CYP2E1 promoter methylation levels did not show significant differences. Here, we compared gene-specific methylation levels at the 5'UTR region, in a subset of workers whom already participated in two former studies, compared to controls.

METHODS

DNA was obtained from whole blood in five different groups: occupationally exposed to a mixture of volatile organic compounds (VOC): high levels (n = 19); low levels (n = 19) and very low levels (n = 17), toluene-exposed workers (n = 19) and control group (n = 19). We performed PCR-pyrosequencing at the 5'UTR region from four genes: CYP2E1, IL-6, SOD1 and TNF-α.

RESULTS

In participants exposed to high levels of a VOC mixture, we found significant differences: lower methylation levels for IL-6, and higher methylation levels for TNF-α compared to controls. In toluene-exposed workers, we found significant, lower methylation levels for CYP2E1 compared to controls.

CONCLUSION

Lower methylation levels at the 5'UTR region from CYP2E1 in toluene exposed-workers, suggests that this epigenetic modification could represent a functional correlate regarding enzymatic activity, as a response to toluene biotransformation.

Multiple correlation analyses revealed complex relationship between DNA methylation and mRNA expression in human peripheral blood mononuclear cells

DNA methylation is an important regulator on the mRNA expression. However, a genome-wide correlation pattern between DNA methylation and mRNA expression in human peripheral blood mononuclear cells (PBMCs) is largely unknown. The comprehensive relationship between mRNA and DNA methylation was explored by using four types of correlation analyses and a genome-wide methylation-mRNA expression quantitative trait locus (eQTL) analysis in PBMCs in 46 unrelated female subjects. An enrichment analysis was performed to detect biological function for the detected genes. Single pair correlation coefficient (r _T1) between methylation level and mRNA is moderate (-0.63-0.62) in intensity, and the negative and positive correlations are nearly equal in quantity. Correlation analysis on each gene (T4) found 60.1% genes showed correlations between mRNA and gene-based methylation at P < 0.05 and more than 5.96% genes presented very strong correlation (R _T4 > 0.8). Methylation sites have regulation effects on mRNA expression in eQTL analysis, with more often observations in region of transcription start site (TSS). The genes under significant methylation regulation both in correlation analysis and eQTL analysis tend to cluster to the categories (e.g., transcription, translation, regulation of transcription) that are essential for maintaining the basic life activities of cells. Our findings indicated that DNA methylation has predictive regulation effect on mRNA with a very complex pattern in PBMCs. The results increased our understanding on correlation of methylation and mRNA and also provided useful clues for future epigenetic studies in exploring biological and disease-related regulatory mechanisms in PBMC.