Pan-cancer genomic analysis links 3'UTR DNA methylation with increased gene expression in T cells

Prognostic effect of DNA methylation of BTG2 gene in Chinese hepatocellular carcinoma

Background

This study aims to develop a prognostic model for overall survival based on potential methylation sites within B-cell translocation gene 2 (BTG2) in Chinese patients with hepatocellular carcinoma (HCC).

Methods

This is a retrospective study. The beta values of nine CpG sites and RSEM normalized count values of BTG2 gene were extracted from the Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) (TCGA-LIHC) dataset, with the beta value representing the methylation level by indicating the ratio of the intensity of the methylated bead type to the combined locus intensity. Pyrosequencing was performed to determine the range of methylation values surrounding cg01798157 site in BTG2 gene. A weighted linear model was developed to predict the overall survival (OS).

Results

The beta value of cg01798157 was significantly negatively associated with the mRNA expression of BTG2 in the TCGA-LIHC dataset (Spearman's rho = -0.5306, P = 2.27 × 10-27). The methylation level of cg01798157 was significantly associated with OS in the cohort of 51 Chinese HCC patients (Hazard ratio = 0.597, 95% CI: 0.434-0.820, P = 0.001). Multivariate Cox regression analysis identified methylation level of cg01798157, cirrhosis, and microvascular invasion as independent prognostic factors. The prognostic efficiency of death risk score was superior to that of cirrhosis or microvascular invasion alone.

Conclusions

The methylation level of cg01798157 in BTG2 may be an epigenetic biomarker in Chinese patients with resectable HCC.

HOXA1 3'UTR Methylation Is a Potential Prognostic Biomarker in Oral Squamous cell Carcinoma

BACKGROUND

HOXA1 is a prognostic marker and a potential predictive biomarker for radioresistance in head and neck tumors. Its overexpression has been associated with promoter methylation and a worse prognosis in oral squamous cell carcinoma (OSCC) patients. However, opposite outcomes are also described. The effect of the methylation of this gene on different gene regions, other than the promoter, remains uncertain. We investigated the methylation profile at different genomic regions of HOXA1 in OSCC and correlated differentially methylated CpG sites with clinicopathological data.

METHODS

The HOXA1 DNA methylation status was evaluated by analyzing data from The Cancer Genome Atlas and three Gene Expression Omnibus datasets. Significant differentially methylated CpG sites were considered with a |∆β| ≥ 0.10 and a Bonferroni-corrected p-value < 0.01. Differentially methylated CpGs were validated by pyrosequencing using two independent cohorts of 15 and 47 OSCC patients, respectively.

RESULTS

Compared to normal tissues, we found significantly higher DNA methylation levels in the 3'UTR region of HOXA1 in OSCC. Higher methylation levels in tumor samples were positively correlated with smoking habits and patients' overall survival.

CONCLUSIONS

Our findings suggest that HOXA1 gene body methylation is a promising prognostic biomarker for OSCC with potential clinical applications in patient monitoring.

The predictive value of BTG1 for the response of newly diagnosed acute myeloid leukemia to decitabine

BACKGROUND

Decitabine has been widely used to treat acute myeloid leukemia (AML); however as AML is a heterogeneous disease, not all patients benefit from decitabine. This study aimed to identify markers for predicting the response to decitabine.

METHODS

An intersection of in vitro experiments and bioinformatics was performed using a combination of epigenetic and transcriptomic analysis. A tumor-suppressor gene associated with methylation and the response to decitabine was screened. Then the sensitivity and specificity of this marker in predicting the response to decitabine was confirmed in 54 samples from newly diagnosed AML patients treated with decitabine plus IA regimen in a clinical trial (ChiCTR2000037928).

RESULTS

In vitro experiments showed that decitabine caused hypomethylation and upregulation of BTG1, while downregulation of BTG1 attenuated the inhibitory effect of decitabine. In newly diagnosed AML patients who received decitabine plus IA regimen, the predictive value of BTG1 to predict complete remission (CR) was assigned with a sensitivity of 86.7% and a specificity of 100.0% when BTG1 expression was < 0.292 (determined using real-time quantitative PCR), with area under the curve (AUC) = 0.933, P = 0.021. The predictive value of BTG1 to predict measurable residual disease (MRD) negativity was assigned with a sensitivity of 100.0% and a specificity of 80.0% when BTG1 expression was < 0.292 (AUC = 0.892, P = 0.012). Patients were divided into low and high BTG1 expression groups according to a cutoff of 0.292, and the CR rate of the low-expression group was significantly higher than that of the high-expression group (97.5% vs. 50%, P < 0.001).

CONCLUSIONS

Low expression of BTG1 was associated with CR and MRD negativity in newly diagnosed AML patients treated with a decitabine-containing regimen, suggesting that BTG1 is a potential marker for predicting the response to decitabine in newly diagnosed AML.

CLINICAL TRIAL REGISTRATION

ChiCTR2000037928.

Blood-based DNA methylation in advanced Nasopharyngeal Carcinoma exhibited distinct CpG methylation signature

The TNM staging system is currently used to detect cancer stages. Regardless, a small proportion of cancer patients recur even after therapy, suggesting more specific molecular tools are required to justify the stage-specific detection and prompt cancer diagnosis. Thus, we aimed to explore the blood-based DNA methylation signature of metastatic nasopharyngeal carcinoma (NPC) to establish a holistic methylation biomarker panel. For the identification of methylation signature, the EPIC BeadChip-based array was performed. Comparative analysis for identifying unique probes, validation, and functional studies was investigated by analyzing GEO and TCGA datasets. We observed 4093 differentially methylated probes (DMPs), 1232 hydroxymethylated probes, and 25 CpG islands. Gene expression study revealed both upregulated and downregulated genes. Correlation analysis suggested a positive (with a positive r, p ≤ 0.05) and negative (with a negative r, p ≤ 0.05) association with different cancers. TFBS analysis exhibited the binding site for many TFs. Furthermore, gene enrichment analysis indicated the involvement of those identified genes in biological pathways. However, blood-based DNA methylation data uncovered a distinct DNA methylation pattern, which might have an additive role in NPC progression by altering the TFs binding. Moreover, based on tissue-specificity, a variation of correlation between methylation and gene expression was noted in different cancers.

Pan-precancer and cancer DNA methylation profiles revealed significant tissue specificity of interrupted biological processes in tumorigenesis

DNA methylation (DNAme) alterations are known to initiate from the precancerous stage of tumorigenesis. Herein, we investigated the global and local patterns of DNAme perturbations in tumorigenesis by analysing the genome-wide DNAme profiles of the cervix, colorectum, stomach, prostate, and liver at precancerous and cancer stages. We observed global hypomethylation in tissues of both two stages, except for the cervix, whose global DNAme level in normal tissue was lower than that of the other four tumour types. For alterations shared by both stages, there were common hyper-methylation (sHyperMethyl) and hypo-methylation (sHypoMethyl) changes, of which the latter type was more frequently identified in all tissues. Biological pathways interrupted by sHyperMethyl and sHypoMethyl alterations demonstrated significant tissue specificity. DNAme bidirectional chaos indicated by the enrichment of both sHyperMethyl and sHypoMethyl changes in the same pathway was observed in most tissues and was a common phenomenon, particularly in liver lesions. Moreover, for the same enriched pathways, different tissues may be affected by distinct DNAme types. For the PI3K-Akt signalling pathway, sHyperMethyl enrichment was observed in the prostate dataset, but sHypoMethyl enrichment was observed in the colorectum and liver datasets. Nevertheless, they did not show an increased possibility in survival prediction of patients in comparison with other DNAme types. Additionally, our study demonstrated that gene-body DNAme changes of tumour suppressor genes and oncogenes may persist from precancerous lesions to the tumour. Overall, we demonstrate the tissue specificity and commonality of cross-stage alterations in DNA methylation profiles in multi-tissue tumorigenesis.

DNA methylation patterns in CD4+ T-cells separate psoriasis patients from healthy controls, and skin psoriasis from psoriatic arthritis

Background

Psoriasis is an autoimmune/inflammatory disorder primarily affecting the skin. Chronic joint inflammation triggers the diagnosis of psoriatic arthritis (PsA) in approximately one-third of psoriasis patients. Although joint disease typically follows the onset of skin psoriasis, in around 15% of cases it is the initial presentation, which can result in diagnostic delays. The pathophysiological mechanisms underlying psoriasis and PsA are not yet fully understood, but there is evidence pointing towards epigenetic dysregulation involving CD4+ and CD8+ T-cells.

Objectives

The aim of this study was to investigate disease-associated DNA methylation patterns in CD4+ T-cells from psoriasis and PsA patients that may represent potential diagnostic and/or prognostic biomarkers.

Methods

PBMCs were collected from 12 patients with chronic plaque psoriasis and 8 PsA patients, and 8 healthy controls. CD4+ T-cells were separated through FACS sorting, and DNA methylation profiling was performed (Illumina EPIC850K arrays). Bioinformatic analyses, including gene ontology (GO) and KEGG pathway analysis, were performed using R. To identify genes under the control of interferon (IFN), the Interferome database was consulted, and DNA Methylation Scores were calculated.

Results

Numbers and proportions of CD4+ T-cell subsets (naïve, central memory, effector memory, CD45RA re-expressing effector memory cells) did not vary between controls, skin psoriasis and PsA patients. 883 differentially methylated positions (DMPs) affecting 548 genes were identified between controls and "all" psoriasis patients. Principal component and partial least-squares discriminant analysis separated controls from skin psoriasis and PsA patients. GO analysis considering promoter DMPs delivered hypermethylation of genes involved in "regulation of wound healing, spreading of epidermal cells", "negative regulation of cell-substrate junction organization" and "negative regulation of focal adhesion assembly". Comparing controls and "all" psoriasis, a majority of DMPs mapped to IFN-related genes (69.2%). Notably, DNA methylation profiles also distinguished skin psoriasis from PsA patients (2,949 DMPs/1,084 genes) through genes affecting "cAMP-dependent protein kinase inhibitor activity" and "cAMP-dependent protein kinase regulator activity". Treatment with cytokine inhibitors (IL-17/TNF) corrected DNA methylation patterns of IL-17/TNF-associated genes, and methylation scores correlated with skin disease activity scores (PASI).

Conclusion

DNA methylation profiles in CD4+ T-cells discriminate between skin psoriasis and PsA. DNA methylation signatures may be applied for quantification of disease activity and patient stratification towards individualized treatment.

Epigenome-wide development and validation of a prognostic methylation score in intrahepatic cholangiocarcinoma based on machine learning strategies

Background

Clinical parameter-based nomograms and staging systems provide limited information for the prediction of survival in intrahepatic cholangiocarcinoma (ICC) patients. In this study, we developed a methylation signature that precisely predicts overall survival (OS) after surgery.

Methods

An epigenome-wide study of DNA methylation based on whole-genome bisulfite sequencing (WGBS) was conducted for two independent cohorts (discovery cohort, n=164; validation cohort, n=170) from three hepatobiliary centers in China. By referring to differentially methylated regions (DMRs), we proposed the concept of prognostically methylated regions (PMRs), which were composed of consecutive prognostically methylated CpGs (PMCs). Using machine learning strategies (Random Forest and the least absolute shrinkage and selector regression), a prognostic methylation score (PMS) was constructed based on 14 PMRs in the discovery cohort and confirmed in the validation cohort.

Results

The C-indices of the PMS for predicting OS in the discovery and validation cohorts were 0.79 and 0.74, respectively. In the whole cohort, the PMS was an independent predictor of OS [hazard ratio (HR) =8.12; 95% confidence interval (CI): 5.48-12.04; P<0.001], and the C-index (0.78) of the PMS was significantly higher than that of the Johns Hopkins University School of Medicine (JHUSM) nomogram (0.69, P<0.001), the Eastern Hepatobiliary Surgery Hospital (EHBSH) nomogram (0.67, P<0.001), American Joint Committee on Cancer (AJCC) tumor-node-metastasis (TNM) staging system (0.61, P<0.001), and MEGNA prognostic score (0.60, P<0.001). The patients in quartile 4 of PMS could benefit from adjuvant therapy (AT) (HR =0.54; 95% CI: 0.32-0.91; log-rank P=0.043), whereas those in the quartiles 1-3 could not. However, other nomograms and staging system failed to do so. Further analyses of potential mechanisms showed that the PMS was associated with tumor biological behaviors, pathway activation, and immune microenvironment.

Conclusions

The PMS could improve the prognostic accuracy and identify patients who would benefit from AT for ICC patients, and might facilitate decisions in treatment of ICC patients.

Reshaping the tumour immune microenvironment in solid tumours via tumour cell and immune cell DNA methylation: from mechanisms to therapeutics

In recent years, the tumour microenvironment (TME) of solid tumours has attracted more and more attention from researchers, especially those non-tumour components such as immune cells. Infiltration of various immune cells causes tumour immune microenvironment (TIME) heterogeneity, and results in different therapeutic effects. Accumulating evidence showed that DNA methylation plays a crucial role in remodelling TIME and is associated with the response towards immune checkpoint inhibitors (ICIs). During carcinogenesis, DNA methylation profoundly changes, specifically, there is a global loss of DNA methylation and increased DNA methylation at the promoters of suppressor genes. Immune cell differentiation is disturbed, and exclusion of immune cells from the TME occurs at least in part due to DNA methylation reprogramming. Therefore, pharmaceutical interventions targeting DNA methylation are promising. DNA methyltransferase inhibitors (DNMTis) enhance antitumor immunity by inducing transcription of transposable elements and consequent viral mimicry. DNMTis upregulate the expression of tumour antigens, mediate immune cells recruitment and reactivate exhausted immune cells. In preclinical studies, DNMTis have shown synergistic effect when combined with immunotherapies, suggesting new strategies to treat refractory solid tumours.

Genome-Wide DNA Methylation Differences between Bos indicus and Bos taurus

Disease risk is a persistent problem in domestic cattle farming, while economic traits are the main concern. This study aimed to reveal the epigenetic basis for differences between zebu (Bos indicus) and taurine cattle (Bos taurus) in disease, disease resistance, and economic traits, and provide a theoretical basis for the genetic improvement of domestic cattle. In this study, whole genome bisulfite sequencing (WGBS) was used to analyze the whole-genome methylation of spleen and liver samples from Yunnan zebu and Holstein cattle. In the genome-wide methylation pattern analysis, it was found that the methylation pattern of all samples was dominated by the CG type, which accounted for >94.9%. The DNA methylation levels of different functional regions and transcriptional elements in the CG background varied widely. However, the methylation levels of different samples in the same functional regions or transcriptional elements did not differ significantly. In addition, we identified a large number of differentially methylation region (DMR) in both the spleen and liver groups, of which 4713 and 4663 were annotated to functional elements, and most of them were annotated to the intronic and exonic regions of genes. GO and KEGG functional analysis of the same differentially methylation region (DMG) in the spleen and liver groups revealed that significantly enriched pathways were involved in neurological, disease, and growth functions. As a result of the results of DMR localization, we screened six genes (DNM3, INPP4B, PLD, PCYT1B, KCNN2, and SLIT3) that were tissue-specific candidates for economic traits, disease, and disease resistance in Yunnan zebu. In this study, DNA methylation was used to construct links between genotypes and phenotypes in domestic cattle, providing useful information for further screening of epigenetic molecular markers in zebu and taurine cattle.

DNA Methylation Variation Is a Possible Mechanism in the Response of Haemaphysalis longicornis to Low-Temperature Stress

Ticks are notorious ectoparasites and transmit the greatest variety of pathogens than any other arthropods. Cold tolerance is a key determinant of tick abundance and distribution. While studies have shown that DNA methylation is one of the important epigenetic regulations found across many species and plays a significant role in their response to low-temperature stress, its role in the response of ticks to low-temperature stress remains unexplored. Herein, we explored the DNA methylation profile of the tick, Haemaphysalis longicornis, exposed to low-temperature stress (4 °C) using whole-genome bisulfite sequencing (WGBS). We found that approximately 0.95% and 0.94% of the genomic C sites were methylated in the control and low-temperature groups, respectively. Moreover, the methylation level under the CG context was about 3.86% and 3.85% in the control and low-temperature groups, respectively. In addition, a total of 6087 differentially methylated regions (DMRs) were identified between the low-temperature and control groups, including 3288 hypermethylated and 2799 hypomethylated DMRs. Further, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially methylated genes revealed that most of the DMGs were significantly enriched in binding and RNA transport pathways. Taken together, this research confirmed, for the first time, the whole genome DNA methylation profile of H. longicornis and provided new insights into the DNA methylation changes relating to low-temperature stress in H. longicornis, as well as provided a foundation for future studies on the epigenetic mechanism underlying the responses of ticks to abiotic stress.

Differential expression of HAVCR2 gene in pan-cancer: A potential biomarker for survival and immunotherapy

T-cell immunoglobulin mucin 3 (TIM-3) has emerged as a promising immune checkpoint target in cancer therapy. However, the profile of the hepatitis A virus cellular receptor 2 (HAVCR2) gene, encoding TIM-3 expression, is still obscure, along with its role in cancer immunity and prognosis. This study comprehensively analyzed HAVCR2 expression patterns in pan-cancer and underlined its potential value for immune checkpoint inhibitor-based immunotherapy. Our results displayed that HAVCR2 was differentially expressed and closely corresponded to survival status in pan-cancer. More importantly, the HAVCR2 expression level was also significantly related to cancer immune infiltration, immune checkpoint genes, and immune marker genes. Enrichment analyses implicated HAVCR2-associated terms in cancer, including immunity, metabolism, and inflammation. Our study demonstrated that HAVCR2 could participate in differing degrees of immune infiltration in tumorigenesis. The highlights of the HAVCR2 pathway revealed that TIM-3 could function as both a biomarker and clinical target to improve the therapeutic efficacy of immunotherapy.

Regulation of DNA methylation during the testicular development of Shaziling pigs

DNA methylation is one of the key epigenetic regulatory mechanisms in development and spermatogenesis. However, the dynamic regulatory mechanisms of genome-wide DNA methylation during testicular development remain largely unknown. Herein, we generated a single-base resolution DNA methylome and transcriptome atlas of precocious porcine testicular tissues across three developmental stages (1, 75, and 150 days old). The results showed that the dynamic methylation patterns were directly related to the expression of the DNMT3A gene. Conjoint analysis revealed a negative regulatory pattern between promoter methylation and the positive regulation of 3'-untranslated region (3'UTR) methylation. Mechanistically, the decrease in promoter methylation affected the upregulation of meiosis-related genes, such as HORMAD1, SPO11, and SYCE1. Demethylation in the 3'UTR induced the downregulation of the INHBA gene and knockdown of INHBA inhibited cell proliferation by reducing the synthesis of activin A. These findings contribute to exploring the regulatory mechanisms of DNA methylation in testicular development.

FAM81A identified as a stemness-related gene by screening DNA methylation sites based on machine learning-accessed stemness in pancreatic cancer

Aim: We thoroughly discuss the interaction between the stemness index and DNA methylation in pancreatic cancer (PC). Materials & methods: First, the stemness indices of PC (denoted mRNAsi and mDNAsi) were calculated using a one-class logistic regression machine-learning algorithm. Second, we screened the central methylation sites associated with stemness and screened out the key genes. We investigated the DNA methylation regulators associated with the key genes. Finally, using CIBERSORT and TIMER, we assessed the influence of stemness indexes and key genes on PC microenvironment formation. Results: In this study we quantified the stemness indices for PC and screened 20 related central DNA methylation sites. Further analysis of the methylation site cg22687244, located in the 3' UTR, revealed that it promoted the expression of the key gene FAM81A. We show that FAM81A may be regulated by DNA methylation regulators. Furthermore, immune cells were found to be more abundant in PC microenvironments with high expression of FAM81A. Conclusion: We report for the first time that the 3' UTR methylation of FAM81A is closely related to PC stemness and contributes to tumor immune infiltration. Therefore FAM81A may serve as a potential marker to guide the treatment of PC.

Genome-wide DNA methylation patterns reveal clinically relevant predictive and prognostic subtypes in human osteosarcoma

Aberrant methylation of genomic DNA has been reported in many cancers. Specific DNA methylation patterns have been shown to provide clinically useful prognostic information and define molecular disease subtypes with different response to therapy and long-term outcome. Osteosarcoma is an aggressive malignancy for which approximately half of tumors recur following standard combined surgical resection and chemotherapy. No accepted prognostic factor save tumor necrosis in response to adjuvant therapy currently exists, and traditional genomic studies have thus far failed to identify meaningful clinical associations. We studied the genome-wide methylation state of primary tumors and tested how they predict patient outcomes. We discovered relative genomic hypomethylation to be strongly predictive of response to standard chemotherapy. Recurrence and survival were also associated with genomic methylation, but through more site-specific patterns. Furthermore, the methylation patterns were reproducible in three small independent clinical datasets. Downstream transcriptional, in vitro, and pharmacogenomic analysis provides insight into the clinical translation of the methylation patterns. Our findings suggest the assessment of genomic methylation may represent a strategy for stratifying patients for the application of alternative therapies.

Examining age-dependent DNA methylation patterns and gene expression in the male and female mouse hippocampus

DNA methylation is a well-characterized epigenetic modification involved in numerous molecular and cellular functions. Methylation patterns have also been associated with aging mechanisms. However, how DNA methylation patterns change within key brain regions involved in memory formation in an age- and sex-specific manner remains unclear. Here, we performed reduced representation bisulfite sequencing (RRBS) from mouse dorsal hippocampus - which is necessary for the formation and consolidation of specific types of memories - in young and aging mice of both sexes. Overall, our findings demonstrate that methylation levels within the dorsal hippocampus are divergent between sexes during aging in genomic features correlating to mRNA functionality, transcription factor binding sites, and gene regulatory elements. These results define age-related changes in the methylome across genomic features and build a foundation for investigating potential target genes regulated by DNA methylation in an age- and sex-specific manner.

Gene Body Methylation of the Lymphocyte-Specific Gene CARD11 Results in Its Overexpression and Regulates Cancer mTOR Signaling

Investigations into the function of nonpromoter DNA methylation have yielded new insights into epigenetic regulation of gene expression. Previous studies have highlighted the importance of distinguishing between DNA methylation in discrete functional regions; however, integrated nonpromoter DNA methylation and gene expression analyses across a wide number of tumor types and corresponding normal tissues have not been performed. Through integrated analysis of gene expression and DNA methylation profiles, we examined 32 tumor types and identified 57 tumor suppressors and oncogenes out of 260 genes exhibiting a correlation of > 0.5 between gene body methylation and gene expression in at least one tumor type. The lymphocyte-specific gene CARD11 exhibits robust association between gene body methylation and expression across 19 of 32 tumor types examined. It is significantly overexpressed in kidney renal cell carcinoma (KIRC) and lung adenocarcinoma (LUAD) tumor tissues in comparison with respective control samples; and is significantly associated with lower overall survival in KIRC. Contrary to its canonical function in lymphocyte NFκB activation, CARD11 activates the mTOR pathway in KIRC and LUAD, resulting in suppressed autophagy. Furthermore, demethylation of a CpG island within the gene body of CARD11 decreases gene expression. Collectively, our study highlights how DNA methylation outside the promoter region can impact tumor progression. IMPLICATIONS: Our study describes a novel regulatory role of gene body DNA methylation-dependent CARD11 expression on mTOR signaling and its impact on tumor progression.

DNA Methylation as a Marker of Body Shape in Premenopausal Women

Preferential accumulation of fat in the gluteo-femoral (GF) depot (pear shape) rather than in the abdominal (A) depot (apple shape), protects against the development of metabolic diseases but the underlying molecular mechanism is still unknown. Recent data, including our work, suggest that differential epigenetic marking is associated with regulation of genes attributed to distinct fat distribution. Here, we aimed to compare the genomic DNA methylation signatures between apple and pear-shaped premenopausal women. To investigate the contribution of upper and lower body fat, we used paired samples of A-FAT and GF-FAT, analyzed on the BeadChip Methylation Array and quantified the differentially methylated sites between the 2 groups of women. We found unique DNA methylation patterns within both fat depots that are significantly different depending on the body fat distribution. Around 60% of the body shape specific DNA methylation sites identified in adipose tissue are maintained ex vivo in cultured preadipocytes. As it has been reported before in other cell types, we found only a hand full of genes showing coordinated differential methylation and expression levels. Finally, we determined that more than 50% of the body shape specific DNA methylation sites could also be detected in whole blood derived DNA. These data reveal a strong DNA methylation program associated with adipose tissue distribution with the possibility that a simple blood test could be used as a predictive diagnostic indicator of young women who are at increased risk for progressing to the apple body shape with a higher risk of developing obesity related complications. Clinical Trial Registration:https://clinicaltrials.gov/ct2/show/NCT02728635 and https://clinicaltrials.gov/ct2/show/NCT02226640, identifiers NCT02728635 and NCT02226640.

Genome-wide DNA methylation profiles provide insight into epigenetic regulation of red and white muscle development in Chinese perch Siniperca chuatsi

Fish skeletal muscles are composed of spatially well-separated fiber types, namely, red and white muscles with different physiological functions and metabolism. To compare the DNA methylation profiles of the two types of muscle tissues and identify potential candidate genes for the muscle growth and development under epigenetic regulation, genome-wide DNA methylation of the red and white muscle in Chinese perch Siniperca chuatsi were comparatively analyzed using bisulfate sequencing methods. An average of 0.9 billion 150-bp paired-end reads were obtained, of which 86% were uniquely mapped to the genome. Methylation mostly occurred at CG sites at a ratio of 94.43% in the red muscle and 93.16% in the white muscle. The mean methylation levels at C-sites were 5.95% in red muscle and 5.83% in white muscle, whereas the mean methylation levels of CG, CHG, and CHH were 73.23%, 0.62%, and 0.67% in red muscle, and 71.01%, 0.62%, and 0.67% in white muscle, respectively. A total of 4192 differentially methylated genes (DMGs) were identified significantly enriched in cell signaling pathways related to skeletal muscle differentiation and growth. Various muscle-related genes, including myosin gene isoforms and regulatory factors, are differentially methylated in the promoter region between the red and white muscles. Further analysis of the transcriptional expression of these genes showed that the muscle regulatory factors (myf5, myog, pax3, pax7, and twitst2) and myosin genes (myh10, myh16, myo18a, myo7a, myo9a, and myl3) were differentially expressed between the two kinds of muscles, consistent with the DNA methylation analysis results. ELISA assays confirmed that the level of 5mC in red muscle was significantly higher than in white muscle (P < 0.05). The RT-qPCR assays revealed that the expression levels of the three DNA methylation transferase (dnmt) subtypes, dnmt1, dnmt3ab, and dnmt3bb1, were significantly higher in red muscle than in white muscle. The higher DNA methylation levels in the red muscle may result from higher DNA methylation transferase expression in the red muscles. Thus, this study might provide a theoretical foundation to better understand epigenetic regulation in the growth and development of red and white muscles in animals, at least in Chinese perch fish.

Polymorphisms in eicosanoid-related biosynthesis enzymes associated with acute urticaria/angioedema induced by nonsteroidal anti-inflammatory drug hypersensitivity

BACKGROUND

Nonsteroidal anti-inflammatory drugs (NSAIDs) are the main triggers of drug hypersensitivity, with NSAID-induced acute urticaria/angioedema (NIUA) the most frequent phenotype. NSAID hypersensitivity is caused by cyclooxygenase 1 inhibition, which leads to an imbalance in prostaglandin (PG) and cysteinyl leukotriene (CysLT) synthesis. As only susceptible individuals develop NSAID hypersensitivity, genetic factors are believed to be involved; however, no study has assessed the overall genetic variability of key enzymes in PG and CysLT synthesis in NSAID hypersensitivity.

OBJECTIVES

To evaluate simultaneously variants in the main genes involved in PG and CysLT biosynthesis in NIUA.

METHODS

Two independent cohorts of patients were recruited in Spain, alongside NSAID-tolerant controls. The discovery cohort included only patients with NIUA; the replication cohort included patients with NSAID-exacerbated respiratory disease (NERD). A set of tagging single-nucleotide polymorphisms (tagSNPs) in PTGS1, PTGS2, ALOX5 and LTC4S was genotyped using mass spectrometry coupled with endpoint polymerase chain reaction.

RESULTS

The study included 1272 individuals. Thirty-five tagSNPs were successfully genotyped in the discovery cohort, with three being significantly associated after Bonferroni correction (rs10306194 and rs1330344 in PTGS1; rs28395868 in ALOX5). These polymorphisms were genotyped in the replication cohort: rs10306194 and rs28395868 remained associated with NIUA, and rs28395868 was marginally associated with NERD. Odds ratios (ORs) in the combined analysis (discovery and replication NIUA populations) were 1·7 for rs10306194 [95% confidence interval (CI) 1·34-2·14; P_corrected = 2·83 × 10^-4 ) and 2·19 for rs28395868 (95% CI 1·43-3·36; P_corrected = 0·002).

CONCLUSIONS

Variants of PTGS1 and ALOX5 may play a role in NIUA and NERD, supporting the proposed mechanisms of NSAID-hypersensitivity and shedding light on their genetic basis.

Genetic Variants in Cytosolic Phospholipase A2 Associated With Nonsteroidal Anti-Inflammatory Drug-Induced Acute Urticaria/Angioedema

Nonsteroidal anti-inflammatory drugs (NSAIDs) are among the main triggers of drug hypersensitivity reactions, probably due to their high consumption worldwide. The most frequent type of NSAID hypersensitivity is NSAID cross-hypersensitivity, in which patients react to NSAIDs from different chemical groups in the absence of a specific immunological response. The underlying mechanism of NSAID cross-hypersensitivity has been linked to cyclooxygenase (COX)-1 inhibition causing an imbalance in the arachidonic acid pathway. Despite NSAID-induced acute urticaria/angioedema (NIUA) being the most frequent clinical phenotype, most studies have focused on NSAID-exacerbated respiratory disease. As NSAID cross-hypersensitivity reactions are idiosyncratic, only appearing in some subjects, it is believed that individual susceptibility is under the influence of genetic factors. Although associations with polymorphisms in genes from the AA pathway have been described, no previous study has evaluated the potential role of cytosolic phospholipase A2 (cPLA2) variants. This enzyme catalyzes the initial hydrolysis of membrane phospholipids to release AA, which can be subsequently metabolized into eicosanoids. Here, we analyzed for the first time the overall genetic variation in the cPLA2 gene (PLA2G4A) in NIUA patients. For this purpose, a set of tagging single nucleotide polymorphisms (tagSNPs) in PLA2G4A were selected using data from Europeans subjects in the 1,000 Genomes Project, and genotyped with the iPlex Sequenom MassArray technology. Two independent populations, each comprising NIUA patients and NSAID-tolerant controls, were recruited in Spain, for the purposes of discovery and replication, comprising a total of 1,128 individuals. Fifty-eight tagSNPs were successfully genotyped in the discovery cohort, of which four were significantly associated with NIUA after Bonferroni correction (rs2049963, rs2064471, rs12088010, and rs12746200). These polymorphisms were then genotyped in the replication cohort: rs2049963 was associated with increased risk for NIUA after Bonferroni correction under the dominant and additive models, whereas rs12088010 and rs12746200 were protective under these two inheritance models. Our results suggest a role for PLA2G4A polymorphisms in NIUA. However, further studies are required to replicate our findings, elucidate the mechanistic role, and evaluate the participation of PLA2G4A variants in other phenotypes induced by NSAID cross-hypersensitivity.

An Unbiased Predictive Model to Detect DNA Methylation Propensity of CpG Islands in the Human Genome

Background:

Epigenetic repression mechanisms play an important role in gene regulation, specifically in cancer development. In many cases, a CpG island’s (CGI) susceptibility or resistance to methylation is shown to be contributed by local DNA sequence features.

Objective:

To develop unbiased machine learning models–individually and combined for different biological features–that predict the methylation propensity of a CGI.

Methods:

We developed our model consisting of CGI sequence features on a dataset of 75 sequences (28 prone, 47 resistant) representing a genome-wide methylation structure. We tested our model on two independent datasets that are chromosome (132 sequences) and disease (70 sequences) specific.

Results:

We provided improvements in prediction accuracy over previous models. Our results indicate that combined features better predict the methylation propensity of a CGI (area under the curve (AUC) ~0.81). Our global methylation classifier performs well on independent datasets reaching an AUC of ~0.82 for the complete model and an AUC of ~0.88 for the model using select sequences that better represent their classes in the training set. We report certain de novo motifs and transcription factor binding site (TFBS) motifs that are consistently better in separating prone and resistant CGIs.

Conclusion:

Predictive models for the methylation propensity of CGIs lead to a better understanding of disease mechanisms and can be used to classify genes based on their tendency to contain methylation prone CGIs, which may lead to preventative treatment strategies. MATLAB® and Python™ scripts used for model building, prediction, and downstream analyses are available at https://github.com/dicleyalcin/methylProp_predictor.

Exploring Common Therapeutic Targets for Neurodegenerative Disorders Using Transcriptome Study

Alzheimer's disease (AD) and Parkinson's disease (PD) are well-known neuronal degenerative disorders that share common pathological events. Approved medications alleviate symptoms but do not address the root cause of the disease. Energy dysfunction in the neuronal population leads to various pathological events and ultimately results in neuronal death. Identifying common therapeutic targets for these disorders may help in the drug discovery process. The Brodmann area 9 (BA9) region is affected in both the disease conditions and plays an essential role in cognitive, motor, and memory-related functions. Analyzing transcriptome data of BA9 provides deep insights related to common pathological pathways involved in AD and PD. In this work, we map the preprocessed BA9 fastq files generated by RNA-seq for disease and control samples with reference hg38 genomic assembly and identify common variants and differentially expressed genes (DEG). These variants are predominantly located in the 3' UTR (non-promoter) region, affecting the conserved transcription factor (TF) binding motifs involved in the methylation and acetylation process. We have constructed BA9-specific functional interaction networks, which show the relationship between TFs and DEGs. Based on expression signature analysis, we propose that MAPK1, VEGFR1/FLT1, and FGFR1 are promising drug targets to restore blood-brain barrier functionality by reducing neuroinflammation and may save neurons.

Role of Methylation in Pro- and Anti-Cancer Immunity

DNA and RNA methylation play a vital role in the transcriptional regulation of various cell types including the differentiation and function of immune cells involved in pro- and anti-cancer immunity. Interactions of tumor and immune cells in the tumor microenvironment (TME) are complex. TME shapes the fate of tumors by modulating the dynamic DNA (and RNA) methylation patterns of these immune cells to alter their differentiation into pro-cancer (e.g., regulatory T cells) or anti-cancer (e.g., CD8+ T cells) cell types. This review considers the role of DNA and RNA methylation in myeloid and lymphoid cells in the activation, differentiation, and function that control the innate and adaptive immune responses in cancer and non-cancer contexts. Understanding the complex transcriptional regulation modulating differentiation and function of immune cells can help identify and validate therapeutic targets aimed at targeting DNA and RNA methylation to reduce cancer-associated morbidity and mortality.

TBC1D16 predicts chemosensitivity and prognosis in adult acute myeloid leukemia (AML) patients

Acute myeloid leukemia (AML) is a hematopoietic disease with poor survival. Chemotherapy resistance is one of the determinant factors influencing AML prognosis. To identify genes possibly affecting the drug responses in AML, the Illumina Infinium MethylationEPIC (850K) was used to screen for differential DNA methylation loci between patients achieved complete remission (CR) or not (non-CR) after induction therapy in 37 AML patients. Then, 32 differentially methylated sites (DMS) were selected for replication in another 86 AML patients by next-generation sequencing. Nine sites including cg03988660, cg16804603, cg18166936, cg11308319, cg09095403, cg18493214, cg01443536, cg16030878 and cg10143426 were replicated. Analysis of the Gene Expression Omnibus (GEO) database showed that mRNA expression of TBC1D16 and HDAC4 was associated with AML prognosis. Methylation level of the cg16030878 in TBC1D16 3'-UTR correlated positively with TBC1D16 mRNA expression in samples both in the TCGA database and clinically collected in the study. Both higher cg16030878 methylation and higher TBC1D16 mRNA expression were associated with increased risk of non-CR and worse overall survival (OS) in AML patients. In AML cells, knockdown of TBC1D16 decreased cell proliferation and ERK phosphorylation levels, as well as increased sensitivity to mitoxantrone and decitabine indicated by IC₅₀. In patients with combined use of decitabine, those patients with CR showed significantly lower TBC1D16 mRNA expression. On the contrary, knockdown of TBC1D16 resulted in decreased sensitivity to cytarabine in U937 cells. Our findings implicated that TBC1D16 is a potential predictor for chemosensitivity and prognosis in adult AML patients.

Analysis of retrotransposon subfamily DNA methylation reveals novel early epigenetic changes in chronic lymphocytic leukemia

Retrotransposons such as LINE-1 and Alu comprise >25% of the human genome. While global hypomethylation of these elements has been widely reported in solid tumours, their epigenetic dysregulation is yet to be characterised in chronic lymphocytic leukemia (CLL), and there has been scant consideration of their evolutionary history that mediates sensitivity to hypomethylation. Here, we developed an approach for locus- and evolutionary subfamily-specific analysis of retrotransposons using the Illumina Infinium Human Methylation 450K microarray platform, which we applied to publicly-available datasets from CLL and other haematological malignancies. We identified 9,797 microarray probes mapping to 117 LINE-1 subfamilies and 13,130 mapping to 37 Alu subfamilies. Of these, 10,782 were differentially methylated (PFDR<0.05) in CLL patients (n=139) compared with healthy individuals (n=14), with enrichment at enhancers (P=0.002). Differential methylation was associated with evolutionary age of LINE-1 (r²=0.31, P=0.003) and Alu (r²=0.74, P=0.002) elements, with greater hypomethylation of older subfamilies (L1M, AluJ). Locus-specific hypomethylation was associated with differential expression of proximal genes, including DCLK2, HK1, ILRUN, TANK, TBCD, TNFRSF1B and TXNRD2, with higher expression of DCLK2 and TNFRSF1B associated with reduced patient survival. Hypomethylation at nine loci was highly frequent in CLL (>90% patients) but not observed in healthy individuals or other leukaemias, and was detectable in blood samples taken prior to CLL diagnosis in 9 of 82 individuals from the Melbourne Collaborative Cohort Study. Our results demonstrate differential methylation of retrotransposons in CLL by their evolutionary heritage that modulates expression of proximal genes.

Tissue- and Sex-Specific DNA Methylation Changes in Mice Perinatally Exposed to Lead (Pb)

Lead (Pb) is a well-known toxicant that interferes with the development of a child’s nervous and metabolic systems and increases the risk of developing diseases later in life. Although studies have investigated epigenetic effects associated with Pb exposure, knowledge of genome-wide changes with in vivo low dose perinatal Pb exposure in multiple tissues is limited. Within the Toxicant Exposures and Responses by Genomic and Epigenomic Regulators of Transcription (TaRGET II) consortium, we utilized a mouse model to investigate tissue- and sex-specific DNA methylation. Dams were assigned to control or Pb-acetate water, respectively. Exposures started 2 weeks prior to mating and continued until weaning at post-natal day 21 (PND21). Liver and blood were collected from PND21 mice, and the DNA methylome was assessed using enhanced reduced representation bisulfite sequencing (ERRBS). We identified ∼1000 perinatal Pb exposure related differentially methylated cytosines (DMCs) for each tissue- and sex-specific comparison, and hundreds of tissue- and sex-specific differentially methylated regions (DMRs). Several mouse imprinted genes were differentially methylated across both tissues in males and females. Overall, our findings demonstrate that perinatal Pb exposure can induce tissue- and sex-specific DNA methylation changes and provide information for future Pb studies in humans.

Identification of epigenetic mechanisms in paddy crop associated with lowering environmentally related cadmium risks to food safety

Cadmium (Cd) is a toxic metal that contributes to human diseases such as pediatric cancer and cardiovascular dysfunction. Epigenetic modification caused by Cd exposure is the major factor in etiology of environmentally-relevant diseases. However, the underlying epigenetic mechanism for Cd uptake and accumulation in food crops, particularly those growing in Cd-contaminated environments, is largely unknown. This study investigated uncharacterized regulatory mechanisms and biological functions of global DNA hypomethylation at CG sites that are associated with gene expression for Cd detoxification and accumulation in the food crop rice. Mutation of the CG maintenance enzyme OsMET1 confers rice tolerance to Cd exposure. Genome-wide analysis of OsMET1 loss of function mutant Osmet1 and its wild type shows numerous loci differentially methylated and upregulated genes for Cd detoxification, transport and accumulation. We functionally identified a new locus for a putative cadmium tolerance factor (here termed as OsCTF) and demonstrated that Cd-induced DNA demethylation is the drive of OsCTF expression. The 3'-UTR of OsCTF is the primary site of DNA and histone (H3K9me2) demethylation, which is associated with higher levels of OsCTF transcripts detected in the Osmet1 and Ossdg714 mutant lines. Mutation of OsCTF in rice led to hypersensitivity to Cd and the Osctf line accumulated more Cd, whereas transfer of OsCTF back to the Osctf mutant completely restored the normal phenotype. Our work unveiled an important epigenetic mechanism and will help develop breeding crops that contribute to food security and better human health.

Overexpression of eIF3D in Lung Adenocarcinoma Is a New Independent Prognostic Marker of Poor Survival

The eukaryotic initiation factor 3 (eIF3) is the largest and most complex translation initiation factor in mammalian cells. It consists of 13 subunits and among which several were implicated to have significant prognostic effects on multiple human cancer entities. To examine the expression profiles of eIF3 subunits and determine their prognostic value in patients with lung adenocarcinoma (LUAD), the genomic data, survival data, and related clinical information were obtained from The Cancer Genome Atlas (TCGA) project for a secondary analysis. The results showed that among ten aberrantly expressed eIF3 subunits in tumours compared with adjacent normal counterparts (p < 0.05), only upregulated eIF3D could predict poor overall survival (OS) outcome independent of multiple clinicopathological parameters (HR = 2.043, 95% CI: 1.132-3.689, p = 0.018). Chi-square analysis revealed that the highly expressed eIF3D group had larger ratios of patients with advanced pathological stage (68/40 vs. 184/206, p = 0.0046), residual tumour (13/4 vs. 163/176, p = 0.0257), and targeted molecular therapy (85/65 vs. 138/164, p = 0.0357). In silico analysis demonstrated that the altered expression of eIF3D was at least regulated by both copy number alterations (CNAs) and the hypomethylation of cg14297023 site. In conclusion, high eIF3D expression might serve as a valuable independent prognostic indicator of shorter OS in patients with LUAD.