Epigenome-wide association study in hepatocellular carcinoma: Identification of stochastic epigenetic mutations through an innovative statistical approach

Molecular Mechanisms in Tumorigenesis of Hepatocellular Carcinoma and in Target Treatments-An Overview

Hepatocellular carcinoma is the most common primary malignancy of the liver, with hepatocellular differentiation. It is ranked sixth among the most common cancers worldwide and is the third leading cause of cancer-related deaths. The most important etiological factors discussed here are viral infection (HBV, HCV), exposure to aflatoxin B1, metabolic syndrome, and obesity (as an independent factor). Directly or indirectly, they induce chromosomal aberrations, mutations, and epigenetic changes in specific genes involved in intracellular signaling pathways, responsible for synthesis of growth factors, cell proliferation, differentiation, survival, the metastasis process (including the epithelial-mesenchymal transition and the expression of adhesion molecules), and angiogenesis. All these disrupted molecular mechanisms contribute to hepatocarcinogenesis. Furthermore, equally important is the interaction between tumor cells and the components of the tumor microenvironment: inflammatory cells and macrophages-predominantly with a pro-tumoral role-hepatic stellate cells, tumor-associated fibroblasts, cancer stem cells, extracellular vesicles, and the extracellular matrix. In this paper, we reviewed the molecular biology of hepatocellular carcinoma and the intricate mechanisms involved in hepatocarcinogenesis, and we highlighted how certain signaling pathways can be pharmacologically influenced at various levels with specific molecules. Additionally, we mentioned several examples of recent clinical trials and briefly described the current treatment protocol according to the NCCN guidelines.

Genome-wide association study reveals a locus in ADARB2 for complete freedom from headache in Danish Blood Donors

Headache disorders are the most common disorders of the nervous system. The lifetime prevalence of headache disorders show that some individuals never experience headache. The etiology of complete freedom from headache is not known. To assess genetic variants associated with complete freedom from headache, we performed a genome-wide association study of individuals who have never experienced a headache. We included 63,992 individuals (2,998 individuals with complete freedom from headache and 60,994 controls) from the Danish Blood Donor Study Genomic Cohort. Participants were included in two rounds, from 2015 to 2018 and in 2020. We discovered a genome-wide significant association, with the lead variant rs7904615[G] in ADARB2 (EAF = 27%, OR = 1.20 [1.13-1.27], p = 3.92 × 10-9). The genomic locus was replicated in a non-overlapping cohort of 13,032 individuals (539 individuals with complete freedom from headache and 12,493 controls) from the Danish Blood Donor Study Genomic Cohort (p < 0.05, two-sided). Participants for the replication were included from 2015 to 2020. In conclusion, we show that complete freedom from headache has a genetic component, and we suggest that ADARB2 is involved in complete freedom from headache. The genomic locus was specific for complete freedom from headache and was not associated with any primary headache disorders.

A novel BRAF::PTPRN2 fusion in meningioma: a case report

Gene fusion events have been linked to oncogenesis in many cancers. However, gene fusions in meningioma are understudied compared to somatic mutations, chromosomal gains/losses, and epigenetic changes. Fusions involving B-raf proto-oncogene, serine/threonine kinase (BRAF) are subtypes of oncogenic BRAF genetic abnormalities that have been reported in certain cases of brain tumors, such as pilocytic astrocytomas. However, BRAF fusions have not been recognized in meningioma. We present the case of an adult female presenting with episodic partial seizures characterized by déjà vu, confusion, and cognitive changes. Brain imaging revealed a cavernous sinus and sphenoid wing mass and she underwent resection. Histopathology revealed a World Health Organization (WHO) grade 1 meningioma. Genetic profiling with next generation sequencing and microarray analysis revealed an in-frame BRAF::PTPRN2 fusion affecting the BRAF kinase domain as well as chromothripsis of chromosome 7q resulting in multiple segmental gains and losses including amplifications of cyclin dependent kinase 6 (CDK6), tyrosine protein-kinase Met (MET), and smoothened (SMO). Elevated pERK staining in tumor cells provided evidence of activated mitogen-activated protein kinase (MAPK) signaling. This report raises the possibility that gene fusion events may be involved in meningioma pathogenesis and warrant further investigation.

The prognostic value of whole-genome DNA methylation in response to Leflunomide in patients with Rheumatoid Arthritis

Objective

Although Leflunomide (LEF) is effective in treating rheumatoid arthritis (RA), there are still a considerable number of patients who respond poorly to LEF treatment. Till date, few LEF efficacy-predicting biomarkers have been identified. Herein, we explored and developed a DNA methylation-based predictive model for LEF-treated RA patient prognosis.

Methods

Two hundred forty-five RA patients were prospectively enrolled from four participating study centers. A whole-genome DNA methylation profiling was conducted to identify LEF-related response signatures via comparison of 40 samples using Illumina 850k methylation arrays. Furthermore, differentially methylated positions (DMPs) were validated in the 245 RA patients using a targeted bisulfite sequencing assay. Lastly, prognostic models were developed, which included clinical characteristics and DMPs scores, for the prediction of LEF treatment response using machine learning algorithms.

Results

We recognized a seven-DMP signature consisting of cg17330251, cg19814518, cg20124410, cg21109666, cg22572476, cg23403192, and cg24432675, which was effective in predicting RA patient's LEF response status. In the five machine learning algorithms, the support vector machine (SVM) algorithm provided the best predictive model, with the largest discriminative ability, accuracy, and stability. Lastly, the AUC of the complex model(the 7-DMP scores with the lymphocyte and the diagnostic age) was higher than the simple model (the seven-DMP signature, AUC:0.74 vs 0.73 in the test set).

Conclusion

In conclusion, we constructed a prognostic model integrating a 7-DMP scores with the clinical patient profile to predict responses to LEF treatment. Our model will be able to effectively guide clinicians in determining whether a patient is LEF treatment sensitive or not.

Notch, SUMOylation, and ESR-Mediated Signalling Are the Main Molecular Pathways Showing Significantly Different Epimutation Scores between Expressing or Not Oestrogen Receptor Breast Cancer in Three Public EWAS Datasets

Oestrogen receptor expression in breast cancer (BC) cells is a marker of high cellular differentiation and allows the identification of two BC groups (ER-positive and ER-negative) that, although not completely homogeneous, differ in biological characteristics, clinical behaviour, and therapeutic options. The study, based on three publicly available EWAS (Epigenetic Wide Association Study) datasets, focuses on the comparison between these two groups of breast cancer using an epimutation score. The score is calculated not only based on the presence of the epimutation, but also on the deviation amplitude of the methylation outlier value. For each dataset, we performed a functional analysis based first on the functional gene region of each annotated gene (we aggregated the data per gene region TSS1500, TSS200, first-exon, and body-gene identified by the information from the Illumina Data Sheet), and then, we performed a pathway enrichment analysis through the REACTOME database based on the genes with the highest epimutation score. Thus, we blended our results and found common pathways for all three datasets. We found that a higher and significant epimutation score due to hypermethylation in ER-positive BC is present in the promoter region of the genes belonging to the SUMOylation pathway, the Notch pathway, the IFN-γ signalling pathway, and the deubiquitination protease pathway, while a higher and significant level of epimutation due to hypomethylation in ER-positive BC is present in the promoter region of the genes belonging to the ESR-mediated pathway. The presence of this state of promoter hypomethylation in the ESR-mediated signalling genes is consistent and coherent with an active signalling pathway mediated by oestrogen function in the group of ER-positive BC. The SUMOylation and Notch pathways are associated with BC pathogenesis and have been found to play distinct roles in the two BC subgroups. We speculated that the altered methylation profile may play a role in regulating signalling pathways with specific functions in the two subgroups of ER BC.

Genomic and transcriptomic profiling reveals key molecules in metastatic potentials and organ-tropisms of hepatocellular carcinoma

Metastasis is a landmark event for rapid postsurgical relapse and death of HCC patients. Although distinct genomic and transcriptomic profiling of HCC metastasis had been reported previously, the causal relationships of somatic mutants, mRNA levels and metastatic potentials were difficult to be established in clinic. Therefore, 11 human HCC cell lines and 7 monoclonal derivatives with definite metastatic potentials and tropisms were subjected to whole exome sequencing (WES) and whole transcriptome sequencing (WTS). TP53, MYO5A, ROS1 and ARID2 were the prominent mutants of metastatic drivers in HCC cells. During HCC clonal evaluation, TP53, MYO5A and ROS1 mutations occurred in the early stage, EXT2 and NIN in the late stage. NF1 mutant was unique in lung tropistic cell lines, RNF126 mutant in lymphatic tropistic ones. PER1, LMO2, GAS7, NR4A3 expression levels were positively associated with relapse-free survival (RFS) of HCC patients. The integrative analysis revealed 58 genes exhibited both somatic mutation and dysregulated mRNA levels in high metastatic cells. Altogether, metastatic drivers could accumulate gradually at different stages during HCC progression, some drivers might modulate HCC metastatic potentials and the others regulate metastatic tropisms.

Dynamic Network Construction for Identifying Early Warning Signals Based On a Data-Driven Approach: Early Diagnosis Biomarker Discovery for Gastric Cancer

During the development of complex diseases, there is a critical transition from one status to another at a tipping point, which can be an early indicator of disease deterioration. To effectively enhance the performance of early risk identification, a novel dynamic network construction algorithm for identifying early warning signals based on a data-driven approach (EWS-DDA) was proposed. In EWS-DDA, the shrunken centroid was introduced to measure dynamic expression changes in assumed pathway reactions during the progression of complex disease for network construction and to define early warning signals by means of a data-driven approach. We applied EWS-DDA to perform a comprehensive analysis of gene expression profiles of gastric cancer (GC) from The Cancer Genome Atlas database and the Gene Expression Omnibus database. Six crucial genes were selected as potential biomarkers for the early diagnosis of GC. The experimental results of statistical analysis and biological analysis suggested that the six genes play important roles in GC occurrence and development. Then, EWS-DDA was compared with other state-of-the-art network methods to validate its performance. The theoretical analysis and comparison results suggested that EWS-DDA has great potential for a more complete presentation of disease deterioration and effective extraction of early warning information.

Tumor Mutational Burden for Predicting Prognosis and Therapy Outcome of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), the primary hepatic malignancy, represents the second-highest cause of cancer-related death worldwide. Many efforts have been devoted to finding novel biomarkers for predicting both patients' survival and the outcome of pharmacological treatments, with a particular focus on immunotherapy. In this regard, recent studies have focused on unravelling the role of tumor mutational burden (TMB), i.e., the total number of mutations per coding area of a tumor genome, to ascertain whether it can be considered a reliable biomarker to be used either for the stratification of HCC patients in subgroups with different responsiveness to immunotherapy, or for the prediction of disease progression, particularly in relation to the different HCC etiologies. In this review, we summarize the recent advances on the study of TMB and TMB-related biomarkers in the HCC landscape, focusing on their feasibility as guides for therapy decisions and/or predictors of clinical outcome.

Stochastic epigenetic mutations as possible explanation for phenotypical discordance among twins with congenital hypothyroidism

PURPOSE

The elevated frequency of discordance for congenital hypothyroidism (CH) phenotype between monozygotic twins suggests the involvement of non-mendelian mechanisms. The aim of the study was to investigate the role of epigenetics in CH pathogenesis.

METHODS

A genome-wide DNA methylation analysis was performed on the peripheral blood of 23 twin pairs (10 monozygotic and 13 dizygotic), 4 concordant and 19 discordant pairs for CH at birth.

RESULTS

Differential methylation analysis did not show significant differences in methylation levels between CH cases and controls, but a different methylation status of several genes may explain the CH discordance of a monozygotic twin couple carrying a monoallelic nonsense mutation of DUOX2. In addition, the median number of hypo-methylated Stochastic Epigenetic Mutations (SEMs) resulted significantly increased in cases compared to controls. The prioritization analysis for CH performed on the genes epimutated exclusively in the cases identified SLC26A4, FOXI1, NKX2-5 and TSHB as the genes with the highest score. The analysis of significantly SEMs-enriched regions led to the identification of two genes (FAM50B and MEG8) that resulted epigenetically dysregulated in cases.

CONCLUSION

Epigenetic modifications may potentially account for CH pathogenesis and explain discordance among monozygotic twins.

Epigenetics of Autism Spectrum Disorders: A Multi-level Analysis Combining Epi-signature, Age Acceleration, Epigenetic Drift and Rare Epivariations Using Public Datasets

BACKGROUND

Epigenetics of Autism Spectrum Disorders (ASD) is still an understudied field. The majority of the studies on the topic used an approach based on mere classification of cases and controls.

OBJECTIVE

The present study aimed at providing a multi-level approach in which different types of epigenetic analysis (epigenetic drift, age acceleration) are combined.

METHODS

We used publicly available datasets from blood (n = 3) and brain tissues (n = 3), separately. Firstly, we evaluated for each dataset and meta-analyzed the differential methylation profile between cases and controls. Secondly, we analyzed age acceleration, epigenetic drift and rare epigenetic variations.

RESULTS

We observed a significant epi-signature of ASD in blood but not in brain specimens. We did not observe significant age acceleration in ASD, while epigenetic drift was significantly higher compared to controls. We reported the presence of significant rare epigenetic variations in 41 genes, 35 of which were never associated with ASD. Almost all genes were involved in pathways linked to ASD etiopathogenesis (i.e., neuronal development, mitochondrial metabolism, lipid biosynthesis and antigen presentation).

CONCLUSION

Our data support the hypothesis of the use of blood epi-signature as a potential tool for diagnosis and prognosis of ASD. The presence of an enhanced epigenetic drift, especially in brain, which is linked to cellular replication, may suggest that alteration in epigenetics may occur at a very early developmental stage (i.e., fetal) when neuronal replication is still high.

Acetyl-CoA metabolism drives epigenome change and contributes to carcinogenesis risk in fatty liver disease

BACKGROUND

The incidence of non-alcoholic fatty liver disease (NAFLD)-associated hepatocellular carcinoma (HCC) is increasing worldwide, but the steps in precancerous hepatocytes which lead to HCC driver mutations are not well understood. Here we provide evidence that metabolically driven histone hyperacetylation in steatotic hepatocytes can increase DNA damage to initiate carcinogenesis.

METHODS

Global epigenetic state was assessed in liver samples from high-fat diet or high-fructose diet rodent models, as well as in cultured immortalized human hepatocytes (IHH cells). The mechanisms linking steatosis, histone acetylation and DNA damage were investigated by computational metabolic modelling as well as through manipulation of IHH cells with metabolic and epigenetic inhibitors. Chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) and transcriptome (RNA-seq) analyses were performed on IHH cells. Mutation locations and patterns were compared between the IHH cell model and genome sequence data from preneoplastic fatty liver samples from patients with alcohol-related liver disease and NAFLD.

RESULTS

Genome-wide histone acetylation was increased in steatotic livers of rodents fed high-fructose or high-fat diet. In vitro, steatosis relaxed chromatin and increased DNA damage marker γH2AX, which was reversed by inhibiting acetyl-CoA production. Steatosis-associated acetylation and γH2AX were enriched at gene clusters in telomere-proximal regions which contained HCC tumour suppressors in hepatocytes and human fatty livers. Regions of metabolically driven epigenetic change also had increased levels of DNA mutation in non-cancerous tissue from NAFLD and alcohol-related liver disease patients. Finally, genome-scale network modelling indicated that redox balance could be a key contributor to this mechanism.

CONCLUSIONS

Abnormal histone hyperacetylation facilitates DNA damage in steatotic hepatocytes and is a potential initiating event in hepatocellular carcinogenesis.

Epigenetic Regulation in Gastroenteropancreatic Neuroendocrine Tumors

Gastroenteropancreatic neuroendocrine neoplasms are a rare, diverse group of neuroendocrine tumors that form in the pancreatic and gastrointestinal tract, and often present with side effects due to hormone hypersecretion. The pathogenesis of these tumors is known to be linked to several genetic disorders, but sporadic tumors occur due to dysregulation of additional genes that regulate proliferation and metastasis, but also the epigenome. Epigenetic regulation in these tumors includes DNA methylation, chromatin remodeling and regulation by noncoding RNAs. Several large studies demonstrate the identification of epigenetic signatures that may serve as biomarkers, and others identify innovative, epigenetics-based targets that utilize both pharmacological and theranostic approaches towards the development of new treatment approaches.

Stochastic Epigenetic Mutations Influence Parkinson's Disease Risk, Progression, and Mortality

BACKGROUND

Stochastic epigenetic mutations (SEM) reflect a deviation from normal site-specific methylation patterns. Epigenetic mutation load (EML) captures the accumulation of SEMs across an individual's genome and may reflect dysfunction of the epigenetic maintenance system in response to epigenetic challenges.

OBJECTIVE

We investigate whether EML is associated with PD risk and time to events (i.e., death and motor symptom decline).

METHODS

We employed logistic regression and Cox proportional hazards regression to assess the association between EML and several outcomes. Our analyses are based on 568 PD patients and 238 controls from the Parkinson's disease, Environment and Genes (PEG) study, for whom blood-based methylation data was available.

RESULTS

We found an association for PD onset and EML in all genes (OR = 1.90; 95%CI 1.52-2.37) and PD-related genes (OR = 1.87; 95%CI 1.50-2.32). EML was also associated with time to a minimum score of 35 points on the motor UPDRS exam (OR = 1.28; 95%CI 1.06-1.56) and time to death (OR = 1.29, 95%CI 1.11-1.49). An analysis of PD related genes only revealed five intragenic hotspots of high SEM density associated with PD risk.

CONCLUSION

Our findings suggest an enrichment of methylation dysregulation in PD patients in general and specifically in five PD related genes. EML may also be associated with time to death and motor symptom progression in PD patients.

A CpG Methylation Signature as a Potential Marker for Early Diagnosis of Hepatocellular Carcinoma From HBV-Related Liver Disease Using Multiplex Bisulfite Sequencing

Background

Aberrant methylation of CpG sites served as an epigenetic marker for building diagnostic, prognostic, and recurrence models for hepatocellular carcinoma (HCC).

Methods

Using Illumina 450K and EPIC Beadchip, we identified 34 CpG sites in peripheral blood mononuclear cell (PBMC) DNA that were differentially methylated in early HCC versus HBV-related liver diseases (HBVLD). We employed multiplex bisulfite sequencing (MBS) based on next-generation sequencing (NGS) to measure methylation of 34 CpG sites in PBMC DNA from 654 patients that were divided into a training set (n = 442) and a test set (n = 212). Using the training set, we selected and built a six-CpG-scorer (namely, cg14171514, cg07721852, cg05166871, cg18087306, cg05213896, and cg18772205), applying least absolute shrinkage and selection operator (LASSO) regression. We performed multivariable analyses of four candidate risk predictors (namely, six-CpG-scorer, age, sex, and AFP level), using 20 times imputation of missing data, non-linearly transformed, and backwards feature selection with logistic regression. The final model’s regression coefficients were calculated according to “Rubin’s Rules”. The diagnostic accuracy of the model was internally validated with a 10,000 bootstrap validation dataset and then applied to the test set for validation.

Results

The area under the receiver operating characteristic curve (AUROC) of the model was 0.81 (95% CI, 0.77–0.85) and it showed good calibration and decision curve analysis. Using enhanced bootstrap validation, adjusted C-statistics and adjusted Brier score were 0.809 and 0.199, respectively. The model also showed an AUROC value of 0.84 (95% CI 0.79–0.88) of diagnosis for early HCC in the test set.

Conclusions

Our model based on the six-CpG-scorer was a reliable diagnosis tool for early HCC from HBVLD. The usage of the MBS method can realize large-scale detection of CpG sites in clinical diagnosis of early HCC and benefit the majority of patients.

DNA methylation-based biomarkers of aging were slowed down in a two-year diet and physical activity intervention trial: the DAMA study

Several biomarkers of healthy aging have been proposed in recent years, including the epigenetic clocks, based on DNA methylation (DNAm) measures, which are getting increasingly accurate in predicting the individual biological age. The recently developed "next-generation clock" DNAmGrimAge outperforms "first-generation clocks" in predicting longevity and the onset of many age-related pathological conditions and diseases. Additionally, the total number of stochastic epigenetic mutations (SEMs), also known as the epigenetic mutation load (EML), has been proposed as a complementary DNAm-based biomarker of healthy aging. A fundamental biological property of epigenetic, and in particular DNAm modifications, is the potential reversibility of the effect, raising questions about the possible slowdown of epigenetic aging by modifying one's lifestyle. Here, we investigated whether improved dietary habits and increased physical activity have favorable effects on aging biomarkers in healthy postmenopausal women. The study sample consists of 219 women from the "Diet, Physical Activity, and Mammography" (DAMA) study: a 24-month randomized factorial intervention trial with DNAm measured twice, at baseline and the end of the trial. Women who participated in the dietary intervention had a significant slowing of the DNAmGrimAge clock, whereas increasing physical activity led to a significant reduction of SEMs in crucial cancer-related pathways. Our study provides strong evidence of a causal association between lifestyle modification and slowing down of DNAm aging biomarkers. This randomized trial elucidates the causal relationship between lifestyle and healthy aging-related epigenetic mechanisms.

DNA methylome in visceral adipose tissue can discriminate patients with and without colorectal cancer

Adipose tissue dysfunction, particularly the visceral (VAT) compartment, has been proposed to play a relevant role in colorectal cancer (CRC) development and progression. Epigenetic mechanisms could be involved in this association. The current study aimed to evaluate if specific epigenetic marks in VAT are associated with colorectal cancer (CRC) to identify epigenetic hallmarks of adipose tissue-related CRC. Epigenome-wide DNA methylation was evaluated in VAT from 25 healthy participants and 29 CRC patients, using the Infinium HumanMethylation450K BeadChip. The epigenome-wide methylation analysis identified 170,184 sites able to perfectly separate the CRC and healthy samples. The differentially methylated CpG sites (DMCpGs) showed a global trend for increased methylated levels in CRC with respect to healthy group. Most of the genes encoded by the DMCpGs belonged to metabolic pathways and cell cycle, insulin resistance, and adipocytokine signalling, as well as tumoural transformation processes. In gene-specific analyses, involved genes biologically relevant for the development of CRC include PTPRN2, MAD1L1, TNXB, DIP2C, INPP5A, HDCA4, PRDM16, RPTOR, ATP11A, TBCD, PABPC3, and IER2. The methylation level of some of them showed a discriminatory capacity for detecting CRC higher than 90%, showing IER2 to have the highest capacity. This study reveals that a specific methylation pattern of VAT is associated with CRC. Some of the epigenetic marks identified could provide useful tools for the prediction and personalized treatment of CRC connected to excess adiposity.

Epigenetic mapping of the somatotropic axis in Nile tilapia reveals differential DNA hydroxymethylation marks associated with growth

In vertebrates, the somatotropic axis comprising the pituitary gland, liver and muscle plays a major role in myogenesis. Its output in terms of muscle growth is highly affected by nutritional and environmental cues, and thus likely epigenetically regulated. Hydroxymethylation is emerging as a DNA modification that modulates gene expression but a holistic characterization of the hydroxymethylome of the somatotropic axis has not been investigated to date. Using reduced representation 5-hydroxymethylcytosine profiling we demonstrate tissue-specific localization of 5-hydroxymethylcytosines at single nucleotide resolution. Their abundance within gene bodies and promoters of several growth-related genes supports their pertinent role in gene regulation. We propose that cytosine hydroxymethylation may contribute to the phenotypic plasticity of growth through epigenetic regulation of the somatotropic axis.

Genetic Variability in Molecular Pathways Implicated in Alzheimer's Disease: A Comprehensive Review

Alzheimer's disease (AD) is a complex neurodegenerative disease, affecting a significant part of the population. The majority of AD cases occur in the elderly with a typical age of onset of the disease above 65 years. AD presents a major burden for the healthcare system and since population is rapidly aging, the burden of the disease will increase in the future. However, no effective drug treatment for a full-blown disease has been developed to date. The genetic background of AD is extensively studied; numerous genome-wide association studies (GWAS) identified significant genes associated with increased risk of AD development. This review summarizes more than 100 risk loci. Many of them may serve as biomarkers of AD progression, even in the preclinical stage of the disease. Furthermore, we used GWAS data to identify key pathways of AD pathogenesis: cellular processes, metabolic processes, biological regulation, localization, transport, regulation of cellular processes, and neurological system processes. Gene clustering into molecular pathways can provide background for identification of novel molecular targets and may support the development of tailored and personalized treatment of AD.

Molecular and immunological paradigms of hepatocellular carcinoma: Special reference to therapeutic approaches

The development of hepatocellular carcinoma (HCC) is a multistep process with a complex interaction of various genetic backgrounds and the tumor microenvironment. In addition to the development of rational approaches to epidemiologic research, early detection, and diagnosis, considerable progress has been made in systemic treatment with molecular-targeted agents for patients with advanced HCC. Moreover, encouraging reports of recent clinical trials of combination therapy with immune-checkpoint inhibitors (ICIs) has raised high hopes. Each HCC is the result of a unique combination of somatic alterations, including genetic, epigenetic, transcriptomic, and metabolic events, leading to conclusive tumoral heterogeneity. Recent advances in comprehensive genetic analysis have accelerated molecular classification and defined subtypes with specific characteristics, including immune-associated molecular profiles reflecting the immune reactivity in the tumor. In considering the development of therapeutic strategies in combination with immunotherapy, proper interpretation of molecular pathological characterization could lead to effective therapeutic deployment and enable individualization of the management of HCC. Here, we review distinctive molecular alterations in the subtype classification of HCC, current therapies, and representative clinical trials with alternative immune-combination approaches from a molecular pathological point.

Genomic Landscape of HCC

INTRODUCTION

Hepatocellular carcinoma (HCC) is a leading cause of cancer related mortality in the world and it has limited treatment options. Understanding the molecular drivers of HCC is important to develop novel biomarkers and therapeutics.

PURPOSE OF REVIEW

HCC arises in a complex background of chronic hepatitis, fibrosis and liver regeneration which lead to genomic changes. Here, we summarize studies that have expanded our understanding of the molecular landscape of HCC.

RECENT FINDINGS

Recent technological advances in next generation sequencing (NGS) have elucidated specific genetic and molecular programs involved in hepatocarcinogenesis. We summarize the major somatic mutations and epigenetic changes have been identified in NGS-based studies. We also describe promising molecular therapies and immunotherapies which target specific genetic and epigenetic molecular events.

SUMMARY

The genomic landscape of HCC is incredibly complex and heterogeneous. Promising new developments are helping us decipher the molecular drivers of HCC and leading to new therapies.

Epigenetic mutation load is weakly correlated with epigenetic age acceleration

DNA methylation (DNAm) age estimators are widely used to study aging-related conditions. It is not yet known whether DNAm age is associated with the accumulation of stochastic epigenetic mutations (SEMs), which reflect dysfunctions of the epigenetic maintenance system. Here, we defined epigenetic mutation load (EML) as the total number of SEMs per individual. We assessed associations between EML and DNAm age acceleration estimators using biweight midcorrelations in four population-based studies (total n = 6,388). EML was not only positively associated with chronological age (meta r = 0.171), but also with four measures of epigenetic age acceleration: the Horvath pan tissue clock, intrinsic epigenetic age acceleration, the Hannum clock, and the GrimAge clock (meta-analysis correlation ranging from r = 0.109 to 0.179). We further conducted pathway enrichment analyses for each participant's SEMs. The enrichment result demonstrated the stochasticity of epigenetic mutations, meanwhile implicated several pathways: signaling, neurogenesis, neurotransmitter, glucocorticoid, and circadian rhythm pathways may contribute to faster DNAm age acceleration. Finally, investigating genomic-region specific EML, we found that EMLs located within regions of transcriptional repression (TSS1500, TSS200, and 1stExon) were associated with faster age acceleration. Overall, our findings suggest a role for the accumulation of epigenetic mutations in the aging process.

DNA Methylome Distinguishes Head and Neck Cancer from Potentially Malignant Oral Lesions and Healthy Oral Mucosa

There is a strong need to find new, good biomarkers of head and neck squamous cell carcinoma (HNSCC) because of the bad prognoses and high mortality rates. The aim of this study was to identify the potential biomarkers in HNSCC that have differences in their DNA methylome and potentially premalignant oral lesions, in comparison to healthy oral mucosa. In this study, 32 oral samples were tested: nine healthy oral mucosae, 13 HNSCC, and 10 oral lesions for DNA methylation by the Infinium MethylationEPIC BeadChip. Our findings showed that a panel of genes significantly hypermethylated in their promoters or specific sites in HNSCC samples in comparison to healthy oral samples, which are mainly oncogenes, receptor, and transcription factor genes, or genes included in cell cycle, transformation, apoptosis, and autophagy. A group of hypomethylated genes in HNSCC, in comparison to healthy oral mucosa, are mainly involved in the host immune response and transcriptional regulation. The results also showed significant differences in gene methylation between HNSCC and potentially premalignant oral lesions, as well as differently methylated genes that discriminate between oral lesions and healthy mucosa. The given methylation panels point to novel potential biomarkers for early diagnostics of HNSCC, as well as potentially premalignant oral lesions.

A comparative analysis of cell-type adjustment methods for epigenome-wide association studies based on simulated and real data sets

Technological advances and reduced costs of high-density methylation arrays have led to an increasing number of association studies on the possible relationship between human disease and epigenetic variability. DNA samples from peripheral blood or other tissue types are analyzed in epigenome-wide association studies (EWAS) to detect methylation differences related to a particular phenotype. Since information on the cell-type composition of the sample is generally not available and methylation profiles are cell-type specific, statistical methods have been developed for adjustment of cell-type heterogeneity in EWAS. In this study we systematically compared five popular adjustment methods: the factored spectrally transformed linear mixed model (FaST-LMM-EWASher), the sparse principal component analysis algorithm ReFACTor, surrogate variable analysis (SVA), independent SVA (ISVA) and an optimized version of SVA (SmartSVA). We used real data and applied a multilayered simulation framework to assess the type I error rate, the statistical power and the quality of estimated methylation differences according to major study characteristics. While all five adjustment methods improved false-positive rates compared with unadjusted analyses, FaST-LMM-EWASher resulted in the lowest type I error rate at the expense of low statistical power. SVA efficiently corrected for cell-type heterogeneity in EWAS up to 200 cases and 200 controls, but did not control type I error rates in larger studies. Results based on real data sets confirmed simulation findings with the strongest control of type I error rates by FaST-LMM-EWASher and SmartSVA. Overall, ReFACTor, ISVA and SmartSVA showed the best comparable statistical power, quality of estimated methylation differences and runtime.

Stochastic Epigenetic Mutations Are Associated with Risk of Breast Cancer, Lung Cancer, and Mature B-cell Neoplasms

BACKGROUND

Age-related epigenetic dysregulations are associated with several diseases, including cancer. The number of stochastic epigenetic mutations (SEM) has been suggested as a biomarker of life-course accumulation of exposure-related DNA damage; however, the predictive role of SEMs in cancer has seldom been investigated.

METHODS

A SEM, at a given CpG site, was defined as an extreme outlier of DNA methylation value distribution across individuals. We investigated the association of the total number of SEMs with the risk of eight cancers in 4,497 case-control pairs nested in three prospective cohorts. Furthermore, we investigated whether SEMs were randomly distributed across the genome or enriched in functional genomic regions.

RESULTS

In the three-study meta-analysis, the estimated ORs per one-unit increase in log(SEM) from logistic regression models adjusted for age and cancer risk factors were 1.25; 95% confidence interval (CI), 1.11-1.41 for breast cancer, and 1.23; 95% CI, 1.07-1.42 for lung cancer. In the Melbourne Collaborative Cohort Study, the OR for mature B-cell neoplasm was 1.46; 95% CI, 1.25-1.71. Enrichment analyses indicated that SEMs frequently occur in silenced genomic regions and in transcription factor binding sites regulated by EZH2 and SUZ12 (P < 0.0001 and P = 0.0005, respectively): two components of the polycomb repressive complex 2 (PCR2). Finally, we showed that PCR2-specific SEMs are generally more stable over time compared with SEMs occurring in the whole genome.

CONCLUSIONS

The number of SEMs is associated with a higher risk of different cancers in prediagnostic blood samples.

IMPACT

We identified a candidate biomarker for cancer early detection, and we described a carcinogenesis mechanism involving PCR2 complex proteins worthy of further investigations.

Epigenome Wide Association and Stochastic Epigenetic Mutation Analysis on Cord Blood of Preterm Birth

Preterm birth (PTB) can be defined as the endpoint of a complex process that could be influenced by maternal and environmental factors. Epigenetics recently emerged as an interesting field of investigation since it represents an important mechanism of regulation. This study evaluates epigenetic impact of preterm birth on DNA methylation. Genome-wide DNAm was measured using the Illumina 450K array in cord blood samples obtained from 72 full term and 18 preterm newborns. Lymphocyte composition was calculated based on specific epigenetic markers that are present on the 450k array. Differential methylation analysis was performed both at site and region level; moreover, stochastic epigenetic mutations (SEMs) were also evaluated. The study showed significant differences in blood cell composition between the two groups. Moreover, after multiple testing correction, statistically significant differences in DNA methylation levels emerged between the two groups both at site and region levels. Results obtained were compared to those reported by previous EWAS, leading to a list of more consistent genes associated with PTB. Finally, the SEMs analysis revealed that the burden of SEMs resulted significantly higher in the preterm group. In conclusion, PTB resulted associated to specific epigenetic signatures that involve immune system. Moreover, SEMs analysis revealed an increased epigenetic drift at birth in the preterm group.

Activation of Adenosine Monophosphate-Activated Protein Kinase Reduces the Onset of Diet-Induced Hepatocellular Carcinoma in Mice

The worldwide obesity and type 2 diabetes epidemics have led to an increase in nonalcoholic fatty liver disease (NAFLD). NAFLD covers a spectrum of hepatic pathologies ranging from simple steatosis to nonalcoholic steatohepatitis, characterized by fibrosis and hepatic inflammation. Nonalcoholic steatohepatitis predisposes to the onset of hepatocellular carcinoma (HCC). Here, we characterized the effect of a pharmacological activator of the intracellular energy sensor adenosine monophosphate–activated protein kinase (AMPK) on NAFLD progression in a mouse model. The compound stimulated fat oxidation by activating AMPK in both liver and skeletal muscle, as revealed by indirect calorimetry. This translated into an ameliorated hepatic steatosis and reduced fibrosis progression in mice fed a diet high in fat, cholesterol, and fructose for 20 weeks. Feeding mice this diet for 80 weeks caused the onset of HCC. The administration of the AMPK activator for 12 weeks significantly reduced tumor incidence and size. Conclusion: Pharmacological activation of AMPK reduces NAFLD progression to HCC in preclinical models.

Whole genome sequencing analysis identifies recurrent structural alterations in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is the predominant type of esophageal cancer in the Asian region, including Japan. A previous study reported mutational landscape of Japanese ESCCs by using exome sequencing. However, somatic structural alterations were yet to be explored. To provide a comprehensive mutational landscape, we performed whole genome sequencing (WGS) analysis of biopsy specimens from 20 ESCC patients in a Japanese population. WGS analysis identified non-silent coding mutations of TP53, ZNF750 and FAT1 in ESCC. We detected six mutational signatures in ESCC, one of which showed significant association with smoking status. Recurrent structural variations, many of which were chromosomal deletions, affected genes such as LRP1B, TTC28, CSMD1, PDE4D, SDK1 and WWOX in 25%–30% of tumors. Somatic copy number amplifications at 11q13.3 (CCND1), 3q26.33 (TP63/SOX2), and 8p11.23 (FGFR1) and deletions at 9p21.3 (CDKN2A) were identified. Overall, these multi-dimensional view of genomic alterations improve the understanding of the ESCC development at molecular level and provides future prognosis and therapeutic implications for ESCC in Japan.

Socioeconomic position, lifestyle habits and biomarkers of epigenetic aging: a multi-cohort analysis

Differences in health status by socioeconomic position (SEP) tend to be more evident at older ages, suggesting the involvement of a biological mechanism responsive to the accumulation of deleterious exposures across the lifespan. DNA methylation (DNAm) has been proposed as a biomarker of biological aging that conserves memory of endogenous and exogenous stress during life.

We examined the association of education level, as an indicator of SEP, and lifestyle-related variables with four biomarkers of age-dependent DNAm dysregulation: the total number of stochastic epigenetic mutations (SEMs) and three epigenetic clocks (Horvath, Hannum and Levine), in 18 cohorts spanning 12 countries.

The four biological aging biomarkers were associated with education and different sets of risk factors independently, and the magnitude of the effects differed depending on the biomarker and the predictor. On average, the effect of low education on epigenetic aging was comparable with those of other lifestyle-related risk factors (obesity, alcohol intake), with the exception of smoking, which had a significantly stronger effect.

Our study shows that low education is an independent predictor of accelerated biological (epigenetic) aging and that epigenetic clocks appear to be good candidates for disentangling the biological pathways underlying social inequalities in healthy aging and longevity.

RNA binding candidates for human ADAR3 from substrates of a gain of function mutant expressed in neuronal cells

Human ADAR3 is a catalytically inactive member of the Adenosine Deaminase Acting on RNA (ADAR) protein family, whose active members catalyze A-to-I RNA editing in metazoans. Until now, the reasons for the catalytic incapability of ADAR3 has not been defined and its biological function rarely explored. Yet, its exclusive expression in the brain and involvement in learning and memory suggest a central role in the nervous system. Here we describe the engineering of a catalytically active ADAR3 enzyme using a combination of computational design and functional screening. Five mutations (A389V, V485I, E527Q, Q549R and Q733D) engender RNA deaminase in human ADAR3. By way of its catalytic activity, the ADAR3 pentamutant was used to identify potential binding targets for wild type ADAR3 in a human glioblastoma cell line. Novel ADAR3 binding sites discovered in this manner include the 3'-UTRs of the mRNAs encoding early growth response 1 (EGR1) and dual specificity phosphatase 1 (DUSP1); both known to be activity-dependent immediate early genes that respond to stimuli in the brain. Further studies reveal that the wild type ADAR3 protein can regulate transcript levels for DUSP1 and EGR1, suggesting a novel role ADAR3 may play in brain function.

Comprehensive longitudinal study of epigenetic mutations in aging

Background

The role of DNA methylation in aging has been widely studied. However, epigenetic mutations, here defined as aberrant methylation levels compared to the distribution in a population, are less understood. Hence, we investigated longitudinal accumulation of epigenetic mutations, using 994 blood samples collected at up to five time points from 375 individuals in old ages.

Results

We verified earlier cross-sectional evidence on the increase of epigenetic mutations with age, and identified important contributing factors including sex, CD19+ B cells, genetic background, cancer diagnosis, and technical artifacts. We further classified epigenetic mutations into High/Low Methylation Outliers (HMO/LMO) according to their changes in methylation, and specifically studied methylation sites (CpGs) that were prone to mutate (frequently mutated CpGs). We validated four epigenetically mutated CpGs using pyrosequencing in 93 samples. Furthermore, by using twins, we concluded that the age-related accumulation of epigenetic mutations was not related to genetic factors, hence driven by stochastic or environmental effects.

Conclusions

Here we conducted a comprehensive study of epigenetic mutation and highlighted its important role in aging process and cancer development.

DNA methylation of individual repetitive elements in hepatitis C virus infection-induced hepatocellular carcinoma

Background

The two most common repetitive elements (REs) in humans, long interspersed nuclear element-1 (LINE-1) and Alu element (Alu), have been linked to various cancers. Hepatitis C virus (HCV) may cause hepatocellular carcinoma (HCC) by suppressing host defenses, through DNA methylation that controls the mobilization of REs. We aimed to investigate the role of RE methylation in HCV-induced HCC (HCV-HCC).

Results

We studied methylation of over 30,000 locus-specific REs across the genome in HCC, cirrhotic, and healthy liver tissues obtained by surgical resection. Relative to normal liver tissue, we observed the largest number of differentially methylated REs in HCV-HCC followed by alcohol-induced HCC (EtOH-HCC). After excluding EtOH-HCC-associated RE methylation (FDR < 0.001) and those unable to be validated in The Cancer Genome Atlas (TCGA), we identified 13 hypomethylated REs (11 LINE-1 and 2 Alu) and 2 hypermethylated REs (1 LINE-1 and 1 Alu) in HCV-HCC (FDR < 0.001). A majority of these REs were located in non-coding regions, preferentially enriched with chromatin repressive marks H3K27me3, and positively associated with gene expression (median correlation r = 0.32 across REs). We further constructed an HCV-HCC RE methylation score that distinguished HCV-HCC (lowest score), HCV-cirrhosis, and normal liver (highest score) in a dose-responsive manner (p for trend < 0.001). HCV-cirrhosis had a lower score than EtOH-cirrhosis (p = 0.038) and HCV-HCC had a lower score than EtOH-HCC in TCGA (p = 0.024).

Conclusions

Our findings indicate that HCV infection is associated with loss of DNA methylation in specific REs, which could implicate molecular mechanisms in liver cancer development. If our findings are validated in larger sample sizes, methylation of these REs may be useful as an early detection biomarker for HCV-HCC and/or a target for prevention of HCC in HCV-positive individuals.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0733-y) contains supplementary material, which is available to authorized users.

Methylation alteration of SHANK1 as a predictive, diagnostic and prognostic biomarker for chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is a clinically heterogeneous disease characterized by the clonal expansion of malignant B cells. To predict the clinical course of the disease, the identification of diagnostic biomarkers is urgently needed. Aberrant methylation patterns may predict CLL development and its course, being very early changes during carcinogenesis. Our aim was to identify CLL specific methylation patterns and to evaluate whether methylation aberrations in selected genes are associated with changes in gene expression. Here, by performing a genome-wide methylation analysis, we identified several CLL-specific methylation alterations. We focused on the most altered one, at a CpG island located in the body of SHANK1 gene, in our CLL cases compared to healthy controls. This methylation alteration was successfully validated in a larger cohort including 139 CLL and 20 control in silico samples. We also found a positive correlation between SHANK1 methylation level and absolute lymphocyte count, in particular CD19+ B cells, in CLL patients. Moreover, we were able to detect gains of methylation at SHANK1 in blood samples collected years prior to diagnosis. Overall, our results suggest methylation alteration at this SHANK1 CpG island as a biomarker for risk and diagnosis of CLL, and also in the personalized quantification of tumor aggressiveness.

Clustered protocadherins methylation alterations in cancer

Background

Clustered protocadherins (PCDHs) map in tandem at human chromosome 5q31 and comprise three multi-genes clusters: α-, β- and γ-PCDH. The expression of this cluster consists of a complex mechanism involving DNA hub formation through DNA-CCTC binding factor (CTCF) interaction. Methylation alterations can affect this interaction, leading to transcriptional dysregulation. In cancer, clustered PCDHs undergo a mechanism of long-range epigenetic silencing by hypermethylation.

Results

In this study, we detected frequent methylation alterations at CpG islands associated to these clustered PCDHs in all the solid tumours analysed (colorectal, gastric and biliary tract cancers, pilocytic astrocytoma), but not hematologic neoplasms such as chronic lymphocytic leukemia. Importantly, several altered CpG islands were associated with CTCF binding sites. Interestingly, our analysis revealed a hypomethylation event in pilocytic astrocytoma, suggesting that in neuronal tissue, where PCDHs are highly expressed, these genes become hypomethylated in this type of cancer. On the other hand, in tissues where PCDHs are lowly expressed, these CpG islands are targeted by DNA methylation. In fact, PCDH-associated CpG islands resulted hypermethylated in gastrointestinal tumours.

Conclusions

Our study highlighted a strong alteration of the clustered PCDHs methylation pattern in the analysed solid cancers and suggested these methylation aberrations in the CpG islands associated with PCDH genes as powerful diagnostic biomarkers.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0695-0) contains supplementary material, which is available to authorized users.

Transcriptomic profiling of neural stem cell differentiation on graphene substrates

Graphene exhibits excellent mechanical strength, electrical conductivity and good biocompatibility, which make it a suitable candidate as a neural interfacing material in regenerative medicine and tissue engineering. Graphene is reported to promote both of neural stem cells (NSCs) proliferation and differentiation. However, the transcriptomes of 2D graphene-regulated NSC differentiation have not yet been investigated. To identify candidate genes, on which graphene may affect, we used next-generation RNA sequencing to analyze the transcriptome of NSCs differentiated for 21 days on a graphene substrate. These NSCs displayed highly enriched and differentially expressed genes compared with traditional cell culture in vitro. Of these, we identified motor protein genes that might regulate NSC differentiation, including cytoplasmic dynein and axonemal dynein genes, Ccdc108, Dnah5, and Dnah11. Furthermore, we analyzed the cell signaling pathway genes that might regulate NSC differentiation, and we constructed a protein-protein interaction network for the genes that are differentially expressed in NSCs on graphene compared to commercial tissue culture polystyrene substrates. We have identified genes potentially regulating the differentiation and migration of NSCs on graphene substrates, and our findings provide mechanistic evidence for the biological activities of graphene, especially in view of graphene-stem cell interactions.

Exposure to polybrominated biphenyl and stochastic epigenetic mutations: application of a novel epigenetic approach to environmental exposure in the Michigan polybrominated biphenyl registry

Endocrine-disrupting compounds are associated with altered epigenetic regulation and adverse health outcomes, although inconsistent results suggest that people have varied responses to the same exposure. Interpersonal variation in response to environmental exposures is not identified using standard, population-based methods. However, methods that capture an individual's response, such as analyzing stochastic epigenetic mutations (SEMs), may capture currently missed effects of environmental exposure. To test whether polybrominated biphenyl (PBB) was associated with SEMs, DNA methylation was measured using Illumina's MethylationEPIC array in PBB-exposed individuals, and SEMs were identified. Association was tested using a linear regression with robust sandwich variance estimators, controlling for age, sex, lipids, and cell types. The number of SEMs was variable (range: 119-18,309), and positively associated with age (p = 1.23e-17), but not with sex (p = 0.97). PBBs and SEMs were only positively associated in people who were older when they were exposed (p = 0.02 vs. p = 0.91). Many subjects had SEMs enriched in biological pathways, particularly in pathways involved with xenobiotic metabolism and endocrine function. Higher number of SEMs was also associated with higher age acceleration (intrinsic: p = 1.70e-3; extrinsic: p = 3.59e-11), indicating that SEMs may be associated with age-related health problems. Finding an association between environmental contaminants and higher SEMs may provide insight into individual differences in response to environmental contaminants, as well as into the biological mechanism behind SEM formation. Furthermore, these results suggest that people may be particularly vulnerable to epigenetic dysregulation from environmental exposures as they age.

Epigenetics of hepatocellular carcinoma

In recent years, large scale genomics and genome-wide studies using comprehensive genomic tools have reshaped our understanding of cancer evolution and heterogeneity. Hepatocellular carcinoma, being one of the most deadly cancers in the world has been well established as a disease of the genome that harbours a multitude of genetic and epigenetic aberrations during the process of liver carcinogenesis. As such, in depth understanding of the cancer epigenetics in cancer specimens and biopsy can be useful in clinical settings for molecular subclassification, prognosis, and prediction of therapeutic responses. In this review, we present a concise discussion on recent progress in the field of liver cancer epigenetics and some of the current works that contribute to the progress of liver cancer therapeutics.

The association of genetically controlled CpG methylation (cg158269415) of protein tyrosine phosphatase, receptor type N2 (PTPRN2) with childhood obesity

Protein tyrosine phosphatase, receptor type N2 (PTPRN2) encodes a major islet autoantigen in type-1 diabetes. Previous genetic studies have shown its significant association with obesity. PTPRN2 plays an important role in epigenetic regulation of metabolic diseases and cancers. We investigated CpG methylations (cg17429772 and cg158269415) in PTPRN2 by pyrosequencing from blood samples of childhood obesity (n = 638). cg158269415 had significant positive correlations with body mass index (BMI) and waist-hip ratio (WHR). Case-control analysis showed that cg158269415 methylation in blood sample was significantly more hypermethylated in obese cases (n = 252), an average of 2.93% more than that that in controls (n = 386). The cg158269415 methylation has a trimodal distribution pattern with strong dependency on nearby located rs1670344 G > A genotype. Correlations of cg158269415 with BMI and WHR were significant and strong in major G allele carriers (GG + GA). Our study showed that an epigenetic association of PTPRN2 gene with childhood obesity was under certain genetic background. The genetic and epigenetic interplay of PTPRN2 gene may implicate a mechanism of childhood obesity. Whether these small changes in DNA methylation from whole blood are causally or consequently related to childhood obesity outcome and their clinical relevance remains to be determined. However, this study presents a promising obesity risk marker that warrants further investigation.

Distinct genome-wide methylation patterns in sporadic and hereditary nonfunctioning pancreatic neuroendocrine tumors

BACKGROUND

Aberrant methylation is a known cause of cancer initiation and/or progression. There are scant data on the genome-wide methylation pattern of nonfunctioning pancreatic neuroendocrine tumors (NFPanNETs) and sporadic and hereditary NFPanNETs.

METHODS

Thirty-three tissue samples were analyzed: they included samples from sporadic (n = 9), von Hippel-Lindau (VHL)-related (n = 10), and multiple endocrine neoplasia type 1 (MEN1)-related NFPanNETs (n = 10) as well as normal islet cells (n = 4) for comparison. Genome-wide CpG methylation profiling was performed with the Infinium MethylationEPIC BeadChip assay and was analyzed with R-based tools.

RESULTS

In unsupervised hierarchical clustering, sporadic and MEN1-related NFPanNETs clustered together, and the VHL group was in a separate cluster. MEN1-related NFPanNETs had a higher rate of hypermethylated CpG sites in comparison with sporadic and VHL-related tumor groups. Differentially methylated region analysis confirmed the higher rate of hypermethylation in MEN1-related tumors. Moreover, in an integrated analysis of gene expression data for the same tumor samples, downregulated gene expression was found in most genes that were hypermethylated. In a CpG island methylator phenotype analysis, 3 genes were identified and confirmed to have downregulated gene expression: secreted frizzle-related protein 5 (SFRP5) in sporadic NFPanNETs and cell division cycle-associated 7-like (CDCA7L) and RNA binding motif 47 (RBM47) in MEN1-related NFPanNETs.

CONCLUSIONS

MEN1 NFPanNETs have a higher rate of geno me-wide hypermethylation than other NFPanNET subtypes. The similarity between the pathways enriched in a methylation analysis of known genes involved in NFPanNET tumorigenesis suggests a key role for aberrant methylation in the pathogenesis of NFPanNETs.

Molecular Scoring of Hepatocellular Carcinoma for Predicting Metastatic Recurrence and Requirements of Systemic Chemotherapy

Hepatocellular carcinoma (HCC) causes one of the most frequent cancer-related deaths; an HCC subset shows rapid progression that affects survival. We clarify molecular features of aggressive HCC, and establish a molecular scoring system that predicts metastasis after curative treatment. In total, 125 HCCs were examined for TP53, CTNNB1, and TERT promoter mutation, methylation of 8 tumor suppressor genes, and 3 repetitive DNA sequences to estimate promoter hypermethylation and global hypomethylation. A fractional allelic loss (FAL) was calculated to represent chromosomal instability through microsatellite analysis. Molecular subclasses were determined using corresponding and hierarchical clustering analyses. Next, twenty-five HCC patients who underwent liver transplantation were analyzed for associations between molecular characteristics and metastatic recurrence; survival analyses were validated using a publicly available dataset of 376 HCC cases from the Cancer Genome Atlas (TCGA). An HCC subtype characterized by TP53 mutation, high FAL, and global hypomethylation was associated with aggressive tumor characteristics, like vascular invasion; CTNNB1 mutation was a feature of the less-progressive phenotype. A number of molecular risk factors, including TP53 mutation, high FAL, significant global hypomethylation, and absence of CTNNB1 mutation, were noted to predict shorter recurrence-free survival in patients who underwent liver transplantation (p = 0.0090 by log-rank test). These findings were validated in a cohort of resected HCC cases from TCGA (p = 0.0076). We concluded that molecular risks determined by common genetic and epigenetic alterations could predict metastatic recurrence after curative treatments, and could be a marker for considering systemic therapy for HCC patients.

Sonoporation by microbubbles as gene therapy approach against liver cancer

INTRODUCTION

An innovative method, known as sonoporation, was used to induce the expression of silenced genes, such as (but not restricted to) TRAIL and p53, in liver cancer cells (HepG2). The principal aim of the present study was the re-activation of silenced apoptotic pathways in liver cancer models, by using diagnostic synovial microbubble as plasmid gene delivery tools in combination with epigenetic treatments.

MATERIAL AND METHODS

HepG2 cells were used as a liver cancer model. Microbubbles (Sonovue®) were chosen as gene deliver system in combination with the sonoporation approach. Plasmid pEGFP-TRAIL and pEGFP-p53 were selected and propagated in Escherichia coli grown in LB broth, in order to obtain the necessary amount.

RESULTS

Sonoporation was induced by using transducer (Sonitron 2000) and, among the several conditions tested, 3 MHz, 51% Duty Cycle, and 5 W/cm2, 30 s resulted as the best parameters. Data collected showed a dose dependent effect in terms of output energy. A transfection efficacy of 30 - 50% was achieved and recombinant gene expression induced apoptotic effects. In order to increase efficacy, we used the histone deacetylase inhibitor (HDACi, entinostat) MS-275, able to activate TRAIL and thus inducing a stronger pro-apoptotic effect in combination with TRAIL-gene re-expression.

CONCLUSION

For the first time, it was shown the possibility to induce the exogenous expression of the pro-apoptotic gene TRAIL and p53 in a liver cancer HepG2 cells via a sonoporation procedure. The epigenetic treatment using HDACi was able to increase the pro-apoptotic effects of the gene therapy.

Multifactorial analysis of the stochastic epigenetic variability in cord blood confirmed an impact of common behavioral and environmental factors but not of in vitro conception

Background

An increased incidence of imprint-associated disorders has been reported in babies born from assisted reproductive technology (ART). However, previous studies supporting an association between ART and an altered DNA methylation status of the conceived babies have been often conducted on a limited number of methylation sites and without correction for critical potential confounders. Moreover, all the previous studies focused on the identification of methylation changes shared among subjects while an evaluation of stochastic differences has never been conducted. This study aims to evaluate the effect of ART and other common behavioral or environmental factors associated with pregnancy on stochastic epigenetic variability using a multivariate approach.

Results

DNA methylation levels of cord blood from 23 in vitro and 41 naturally conceived children were analyzed using the Infinium HumanMethylation450 BeadChips. After multiple testing correction, no statistically significant difference emerged in the number of cord blood stochastic epigenetic variations or in the methylation levels between in vitro- and in vivo-conceived babies. Conversely, four multiple factor analysis dimensions summarizing common phenotypic, behavioral, or environmental factors (cord blood cell composition, pre or post conception supplementation of folates, birth percentiles, gestational age, cesarean section, pre-gestational mother’s weight, parents’ BMI and obesity status, presence of adverse pregnancy outcomes, mother’s smoking status, and season of birth) were significantly associated with stochastic epigenetic variability. The stochastic epigenetic variation analysis allowed the identification of a rare imprinting defect in the locus GNAS in one of the babies belonging to the control population, which would not have emerged using a classical case-control association analysis.

Conclusions

We confirmed the effect of several common behavioral or environmental factors on the epigenome of newborns and described for the first time an epigenetic effect related to season of birth. Children born after ART did not appear to have an increased risk of genome-wide changes in DNA methylation either at specific loci or randomly scattered throughout the genome. The inability to identify differences between cases and controls suggests that the number of stochastic epigenetic variations potentially induced by ART was not greater than that naturally produced in response to maternal behavior or other common environmental factors.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0510-3) contains supplementary material, which is available to authorized users.

Colorectal cancer early methylation alterations affect the crosstalk between cell and surrounding environment, tracing a biomarker signature specific for this tumor

Colorectal cancer (CRC) develops through the accumulation of both genetic and epigenetic alterations. However, while the former are already used as prognostic and predictive biomarkers, the latter are less well characterized. Here, performing global methylation analysis on both CRCs and adenomas by Illumina Infinium HumanMethylation450 Bead Chips, we identified a panel of 74 altered CpG islands, demonstrating that the earliest methylation alterations affect genes coding for proteins involved in the crosstalk between cell and surrounding environment. The panel discriminates CRCs and adenomas from peritumoral and normal mucosa with very high specificity (100%) and sensitivity (99.9%). Interestingly, over 70% of the hypermethylated islands resulted in downregulation of gene expression. To establish the possible usefulness of these non-invasive markers for detection of colon cancer, we selected three biomarkers and identified the presence of altered methylation in stool DNA and plasma cell-free circulating DNA from CRC patients.

Primary Osteoporosis Is Not Reflected by Disease-Specific DNA Methylation or Accelerated Epigenetic Age in Blood

Osteoporosis is an age-related metabolic bone disease. Hence, osteoporotic patients might suffer from molecular features of accelerated aging, which is generally reflected by specific age-associated DNA methylation (DNAm) changes. In this study, we analyzed genomewide DNAm profiles of peripheral blood from patients with manifest primary osteoporosis and non-osteoporotic controls. Statistical analysis did not reveal any individual CG dinucleotides (CpG sites) with significant aberrant DNAm in osteoporosis. Subsequently, we analyzed if age-associated DNAm patterns are increased in primary osteoporosis (OP). Using three independent age-predictors we did not find any evidence for accelerated epigenetic age in blood of osteoporotic patients. Taken together, osteoporosis is not reflected by characteristic DNAm patterns of peripheral blood that might be used as biomarker for the disease. The prevalence of osteoporosis is age-associated-but it is not associated with premature epigenetic aging in peripheral blood. © 2017 American Society for Bone and Mineral Research.