Higher Nevus Count Exhibits a Distinct DNA Methylation Signature in Healthy Human Skin: Implications for Melanoma

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

Genetic effects on the skin methylome in healthy older twins

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

Map of epigenetic age acceleration: A worldwide analysis

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Development of an epigenetic clock resistant to changes in immune cell composition

Epigenetic clocks are age predictors that use machine-learning models trained on DNA CpG methylation values to predict chronological or biological age. Increases in predicted epigenetic age relative to chronological age (epigenetic age acceleration) are connected to aging-associated pathologies, and changes in epigenetic age are linked to canonical aging hallmarks. However, epigenetic clocks rely on training data from bulk tissues whose cellular composition changes with age. Here, we found that human naive CD8+ T cells, which decrease in frequency during aging, exhibit an epigenetic age 15-20 years younger than effector memory CD8+ T cells from the same individual. Importantly, homogenous naive T cells isolated from individuals of different ages show a progressive increase in epigenetic age, indicating that current epigenetic clocks measure two independent variables, aging and immune cell composition. To isolate the age-associated cell intrinsic changes, we created an epigenetic clock, the IntrinClock, that did not change among 10 immune cell types tested. IntrinClock shows a robust predicted epigenetic age increase in a model of replicative senescence in vitro and age reversal during OSKM-mediated reprogramming.

Epistemic uncertainty challenges aging clock reliability in predicting rejuvenation effects

Epigenetic aging clocks have been widely used to validate rejuvenation effects during cellular reprogramming. However, these predictions are unverifiable because the true biological age of reprogrammed cells remains unknown. We present an analytical framework to consider rejuvenation predictions from the uncertainty perspective. Our analysis reveals that the DNA methylation profiles across reprogramming are poorly represented in the aging data used to train clock models, thus introducing high epistemic uncertainty in age estimations. Moreover, predictions of different published clocks are inconsistent, with some even suggesting zero or negative rejuvenation. While not questioning the possibility of age reversal, we show that the high clock uncertainty challenges the reliability of rejuvenation effects observed during in vitro reprogramming before pluripotency and throughout embryogenesis. Conversely, our method reveals a significant age increase after in vivo reprogramming. We recommend including uncertainty estimation in future aging clock models to avoid the risk of misinterpreting the results of biological age prediction.

Biology and genetics of acquired and congenital melanocytic naevi

Acquired and congenital melanocytic naevi are common benign neoplasms. Understanding their biology and genetics will help clinicians and pathologists correctly diagnose melanocytic tumours, and generate insights into naevus aetiology and melanomagenesis. Genomic data from published studies analysing acquired and congenital melanocytic naevi, including oncogenic driver mutations, common melanoma associated mutations, copy number aberrations, somatic mutation signature patterns, methylation profile, and single nucleotide polymorphisms, were reviewed. Correlation of genomic changes to dermoscopic features, particular anatomic sites and total body naevus counts, was also performed. This review also highlights current scientific theories and evidence concerning naevi growth arrest. Acquired and congenital melanocytic naevi show simple genomes, typically characterised by mutually exclusive single oncogenic driver mutations in either BRAF or NRAS genes. Genomic differences exist between acquired and congenital naevi, common and dysplastic naevi, and by dermoscopic features. Acquired naevi show a higher rate of BRAF hotspot mutations and a lower rate of NRAS hotspot mutations compared to congenital naevi. Dysplastic naevi show upregulation of follicular keratinocyte-related genes compared to common naevi. Anatomical locations and DNA signatures of naevi implicates ultraviolet radiation and non-ultraviolet radiation pathways in naevogenesis. DNA driver point mutations in acquired and congenital melanocytic naevi have been well characterised. Future research is required to better understand transcriptional and epigenetic changes in naevi, as well as those regulating naevus growth arrest and cell environment signalling.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Genome-scale DNA methylation analysis identifies repeat element alterations that modulate the genomic stability of melanocytic nevi

Acquired melanocytic nevi grow and persist in a stable form into adulthood. Using genome-wide methylation profiling, we evaluated 32 histopathologically and dermoscopically characterized nevi, to identify key epigenetic regulatory mechanisms involved in nevogenesis. Benign (69% globular and 31% non-specific dermoscopic pattern) and dysplastic (95% reticular/nonspecific dermoscopic pattern) nevi were dissimilar with only two shared differentially methylated (DM) loci. Benign nevi demonstrated an increase in both genome-scale methylation and methylation of Alu/LINE-1 retrotransposable elements, a marker of genomic stability, as well as global methylation. In contrast, dysplastic nevi showed evidence for genomic instability via hypomethylation of Alu/LINE-1 (Alu; P=0.00019 and LINE-1; P=0.000035). Using dermoscopic classifications, reticular/non-specific patterned nevi had 59,572 CpG DM loci (Q < 0.05), whereas globular nevi had no significant DM loci. In reticular/non-specific patterned nevi, the tumor suppressor PTEN had the greatest proportion of hypermethylated CpG loci in its promoter region compared to all other assayed gene promoters. The relative activity of reticular/non-specific nevi was evidenced by 50,720 hypomethylated loci being enriched for accessible chromatin, and 8,852 hypermethylated loci strongly enriched, for example, marks of active gene promoters, which suggests that gain of DNA methylation observed in these nevus types plays a role in gene regulation.

Melanoma susceptibility: an update on genetic and epigenetic findings

Malignant melanoma is one of the most highly ranked cancers in terms of years of life lost. Hereditary melanoma with its increased familial susceptibility is thought to affect up to 12% of all melanoma patients. In the past, only a few high-penetrance genes associated with familial melanoma, such as CDKN2A and CDK4, have been clinically tested. However, findings now indicate that melanoma is a cancer most likely to develop not only due to high-penetrance variants but also due to polygenic inheritance patterns, leaving no clear division between the hereditary and sporadic development of malignant melanoma. Various pathogenic low-penetrance variants were recently discovered through genome-wide association studies, and are now translated into polygenic risk scores. These can show superior sensitivity rates for the prediction of melanoma susceptibility and related mixed cancer syndromes than risk scores based on phenotypic traits of the patients, with odds ratios of up to 5.7 for patients in risk groups. In addition to describing genetic findings, we also review the first results of epigenetic research showing constitutional methylation changes that alter the susceptibility to cutaneous melanoma and its risk factors.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Investigating DNA methylation as a potential mediator between pigmentation genes, pigmentary traits and skin cancer

Pigmentation characteristics are well-known risk factors for skin cancer. Polymorphisms in pigmentation genes have been associated with these traits and with the risk of malignancy. However, the functional relationship between genetic variation and disease is still unclear. This study aims to assess whether pigmentation SNPs are associated with pigmentary traits and skin cancer via DNA methylation (DNAm). Using a meta-GWAS of whole-blood DNAm from 36 European cohorts (N = 27,750; the Genetics of DNA Methylation Consortium, GoDMC), we found that 19 out of 27 SNPs in 10 pigmentation genes were associated with 391 DNAm sites across 30 genomic regions. We examined the effect of 25 selected DNAm sites on pigmentation traits, sun exposure phenotypes and skin cancer and on gene expression in whole blood. We uncovered an association of DNAm site cg07402062 with red hair in the Avon Longitudinal Study of Parents and Children (ALSPAC). We also found that the expression of ASIP and CDK10 was associated with hair colour, melanoma and basal cell carcinoma. Our results indicate that DNAm and expression of pigmentation genes may play a role as potential mediators of the relationship between genetic variants, pigmentation phenotypes and skin cancer and thus deserve further scrutiny.

Terminal keratinocyte differentiation in vitro is associated with a stable DNA methylome

The epidermal compartment of the skin is regenerated constantly by proliferation of epidermal keratinocytes. Differentiation of a subset of these keratinocytes allows the epidermis to retain its barrier properties. Regulation of keratinocyte fate-whether to remain proliferative or terminally differentiate-is complex and not fully understood. The objective of our study was to assess if DNA methylation changes contribute to the regulation of keratinocyte fate. We employed genome-wide MethylationEPIC beadchip array measuring approximately 850 000 probes combined with RNA sequencing of in vitro cultured non-differentiated and terminally differentiated adult human primary keratinocytes. We did not observe a correlation between methylation status and transcriptome changes. Moreover, only two differentially methylated probes were detected, of which one was located in the TRIM29 gene. Although TRIM29 knock-down resulted in lower expression levels of terminal differentiation genes, these changes were minor. From these results, we conclude that-in our in vitro experimental setup-it is unlikely that changes in DNA methylation have an important regulatory role in terminal keratinocyte differentiation.

Epigenetic regulation of melanogenesis

Melanogenesis is a complex process in which melanin is synthesized in melanocytes and transported to keratinocytes, which involves multiple genes and signaling pathways. Epigenetics refers to the potential genetic changes that affect gene expression without involving changes in the original sequence of DNA nucleotides. DNA methylation regulates the expression of key genes such as tyrosinase (TYR), tyrosinase-related protein 1 (TYRP1), dopachrome tautomerase (DCT) and microphthalmia-associated transcription factor (MITF), as well as paracrine factors such as stem cell factor (SCF) and endothelin-1 (ET-1) in melanogenesis. Potential DNA methylation sites are present in the genes of melanogenesis-related signaling pathways such as "Wnt", "PI3K/Akt/CREB" and "MAPK". H3K27 acetylation is abundant in melanogenesis-related genes. Both the upstream activation and downstream regulation of MITF depend on histone acetyltransferase CBP/p300, and pH-induced H3K27 acetylation may be the amplifying mechanism of MITF's effect. HDAC1 and HDAC10 catalyze histone deacetylation of melanogenesis-related gene promoters. Chromatin remodelers SWI/SNF complex and ISWI complex use the energy of ATP hydrolysis to rearrange nucleosomes, while their active subunits BRG1, BRM and BPTF, act as activators and cofactors of MITF. MicroRNAs (miRNAs) can directly target a large number of melanogenesis-related genes, while long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) regulate melanogenesis in a variety of ways. Interactions exist among the epigenetic mechanisms of melanogenesis. For example, the methyl CpG binding domain protein 2 (MeCP2) links DNA methylation, histone deacetylation, and histone methylation. Epigenetic-based therapy provides novel opportunities for treating dermatoses that are caused by pigmentation disturbances. This review summarizes the epigenetic regulation mechanisms of melanogenesis, and examines the pathogenesis and treatment of epigenetics in pigmentation disorders.

The role of MC1R gene variants and phenotypical features in predicting high nevus count

Variants in the Melanocortin 1 Receptor (MC1R) gene have been associated with an increased risk of melanoma, but the role in nevus count is unclear. We investigated if specific MC1R gene variants or the number of MC1R gene variants and phenotypical features were associated with nevus count. A total of 494 participants of the 'Leiden skin cancer study' were included and the MC1R gene coding sequence was analysed by single-strand conformation polymorphism analysis followed by sequencing of unknown variants. The association between MC1R gene variants and nevus count and the association between age, gender and phenotypical features and nevus count were studied using the Chi-square test. Study of nine frequently occurring MC1R gene variants in participants without skin cancer (n = 203) showed that the 'r' Val60Leu variant was significantly associated with high nevus count (>50 nevi) (P = 0.017). This association was very strong among women (P < 0.001), but not present among men. Having one or two MC1R variants in general did not show a significant difference in the nevus count. Hair colour, skin type, eye colour and age were not significantly associated with nevus count, whereas gender showed a significant association (P = 0.008), with the highest nevus counts in female. The Val60Leu variant of the MC1R gene could be a promising candidate as an independent predictor of high nevus count, particularly in women. This information about the genetic makeup could promote personalized follow-up strategies and might help to prevent skin cancer in the future.

Hard wiring of normal tissue-specific chromosome-wide gene expression levels is an additional factor driving cancer type-specific aneuploidies

Background

Many carcinomas have recurrent chromosomal aneuploidies specific to the tissue of tumor origin. The reason for this specificity is not completely understood.

Methods

In this study, we looked at the frequency of chromosomal arm gains and losses in different cancer types from the The Cancer Genome Atlas (TCGA) and compared them to the mean gene expression of each chromosome arm in corresponding normal tissues of origin from the Genotype-Tissue Expression (GTEx) database, in addition to the distribution of tissue-specific oncogenes and tumor suppressors on different chromosome arms.

Results

This analysis revealed a complex picture of factors driving tumor karyotype evolution in which some recurrent chromosomal copy number reflect the chromosome arm-wide gene expression levels of the their normal tissue of tumor origin.

Conclusions

We conclude that the cancer type-specific distribution of chromosomal arm gains and losses is potentially “hardwiring” gene expression levels characteristic of the normal tissue of tumor origin, in addition to broadly modulating the expression of tissue-specific tumor driver genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00905-y.

Identification of Differentially Methylated CpG Sites in Fibroblasts from Keloid Scars

As a part of an abnormal healing process of dermal injuries and irritation, keloid scars arise on the skin as benign fibroproliferative tumors. Although the etiology of keloid scarring remains unsettled, considerable recent evidence suggested that keloidogenesis may be driven by epigenetic changes, particularly, DNA methylation. Therefore, genome-wide scanning of methylated cytosine-phosphoguanine (CpG) sites in extracted DNA from 12 keloid scar fibroblasts (KF) and 12 control skin fibroblasts (CF) (six normal skin fibroblasts and six normotrophic fibroblasts) was conducted using the Illumina Human Methylation 450K BeadChip in two replicates for each sample. Comparing KF and CF used a Linear Models for Microarray Data (Limma) model revealed 100,000 differentially methylated (DM) CpG sites, 20,695 of which were found to be hypomethylated and 79,305 were hypermethylated. The top DM CpG sites were associated with TNKS2, FAM45B, LOC723972, GAS7, RHBDD2 and CAMKK1. Subsequently, the most functionally enriched genes with the top 100 DM CpG sites were significantly (p ≤ 0.05) associated with SH2 domain binding, regulation of transcription, DNA-templated, nucleus, positive regulation of protein targeting to mitochondrion, nucleoplasm, Swr1 complex, histone exchange, and cellular response to organic substance. In addition, NLK, CAMKK1, LPAR2, CASP1, and NHS showed to be the most common regulators in the signaling network analysis. Taken together, these findings shed light on the methylation status of keloids that could be implicated in the underlying mechanism of keloid scars formation and remission.

TwinsUK: The UK Adult Twin Registry Update

TwinsUK is the largest cohort of community-dwelling adult twins in the UK. The registry comprises over 14,000 volunteer twins (14,838 including mixed, single and triplets); it is predominantly female (82%) and middle-aged (mean age 59). In addition, over 1800 parents and siblings of twins are registered volunteers. During the last 27 years, TwinsUK has collected numerous questionnaire responses, physical/cognitive measures and biological measures on over 8500 subjects. Data were collected alongside four comprehensive phenotyping clinical visits to the Department of Twin Research and Genetic Epidemiology, King’s College London. Such collection methods have resulted in very detailed longitudinal clinical, biochemical, behavioral, dietary and socioeconomic cohort characterization; it provides a multidisciplinary platform for the study of complex disease during the adult life course, including the process of healthy aging. The major strength of TwinsUK is the availability of several ‘omic’ technologies for a range of sample types from participants, which includes genomewide scans of single-nucleotide variants, next-generation sequencing, metabolomic profiles, microbiomics, exome sequencing, epigenetic markers, gene expression arrays, RNA sequencing and telomere length measures. TwinsUK facilitates and actively encourages sharing the ‘TwinsUK’ resource with the scientific community — interested researchers may request data via the TwinsUK website (http://twinsuk.ac.uk/resources-for-researchers/access-our-data/) for their own use or future collaboration with the study team. In addition, further cohort data collection is planned via the Wellcome Open Research gateway (https://wellcomeopenresearch.org/gateways). The current article presents an up-to-date report on the application of technological advances, new study procedures in the cohort and future direction of TwinsUK.

DNA methylation associates with survival in non-metastatic clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is the most common subtype among renal cancer and is associated with poor prognosis if metastasized. Up to one third of patients with local disease at diagnosis will develop metastasis after nephrectomy, and there is a need for new molecular markers to identify patients with high risk of tumor progression. In the present study, we performed genome-wide promoter DNA methylation analysis at diagnosis to identify DNA methylation profiles associated with risk for progress.

METHOD

Diagnostic tissue samples from 115 ccRCC patients were analysed by Illumina HumanMethylation450K arrays and methylation status of 155,931 promoter associated CpGs were related to genetic aberrations, gene expression and clinicopathological parameters.

RESULTS

The ccRCC samples separated into two clusters (cluster A/B) based on genome-wide promoter methylation status. The samples in these clusters differed in tumor diameter (p < 0.001), TNM stage (p < 0.001), morphological grade (p < 0.001), and patients outcome (5 year cancer specific survival (pCSS5yr) p < 0.001 and cumulative incidence of progress (pCIP5yr) p < 0.001. An integrated genomic and epigenomic analysis in the ccRCCs, revealed significant correlations between the total number of genetic aberrations and total number of hypermethylated CpGs (R = 0.435, p < 0.001), and predicted mitotic age (R = 0.407, p < 0.001). We identified a promoter methylation classifier (PMC) panel consisting of 172 differently methylated CpGs accompanying progress of disease. Classifying non-metastatic patients using the PMC panel showed that PMC high tumors had a worse prognosis compared with the PMC low tumors (pCIP5yr 38% vs. 8%, p = 0.001), which was confirmed in non-metastatic ccRCCs in the publically available TCGA-KIRC dataset (pCIP5yr 39% vs. 16%, p < 0.001).

CONCLUSION

DNA methylation analysis at diagnosis in ccRCC has the potential to improve outcome-prediction in non-metastatic patients at diagnosis.

The relationship between naevus count, memory function and telomere length in the Twins UK cohort

The presence of a skin-brain connection whereby alterations in the skin can inform on mechanisms underlying neurodegenerative diseases is increasingly recognized. In this study, we used a discovery (n = 321) and replication (n = 147) sample from the Twins UK population to test the association between naevus count and memory function, and its mediation by telomeres. Memory function was assessed in 1999 and 2009 using the paired associates learning test (PAL), while naevus count and leucocyte telomere length (LTL, assessed by the terminal restriction fragment assay) were measured once. Higher baseline naevus count was significantly associated with fewer errors at the baseline and follow-up PAL, as well as with change in PAL score over 10 years. This association was significantly attenuated after adjustment for LTL. The significant association between naevus count and PAL score was reproduced in the replication sample. These findings suggest that melanocytes might be used as model system to study the biological ageing pathways involved in neurodegeneration.

Smoking induces coordinated DNA methylation and gene expression changes in adipose tissue with consequences for metabolic health

Background

Tobacco smoking is a risk factor for multiple diseases, including cardiovascular disease and diabetes. Many smoking-associated signals have been detected in the blood methylome, but the extent to which these changes are widespread to metabolically relevant tissues, and impact gene expression or metabolic health, remains unclear.

Methods

We investigated smoking-associated DNA methylation and gene expression variation in adipose tissue biopsies from 542 healthy female twins. Replication, tissue specificity, and longitudinal stability of the smoking-associated effects were explored in additional adipose, blood, skin, and lung samples. We characterized the impact of adipose tissue smoking methylation and expression signals on metabolic disease risk phenotypes, including visceral fat.

Results

We identified 42 smoking-methylation and 42 smoking-expression signals, where five genes (AHRR, CYP1A1, CYP1B1, CYTL1, F2RL3) were both hypo-methylated and upregulated in current smokers. CYP1A1 gene expression achieved 95% prediction performance of current smoking status. We validated and replicated a proportion of the signals in additional primary tissue samples, identifying tissue-shared effects. Smoking leaves systemic imprints on DNA methylation after smoking cessation, with stronger but shorter-lived effects on gene expression. Metabolic disease risk traits such as visceral fat and android-to-gynoid ratio showed association with methylation at smoking markers with functional impacts on expression, such as CYP1A1, and at tissue-shared smoking signals, such as NOTCH1. At smoking-signals, BHLHE40 and AHRR DNA methylation and gene expression levels in current smokers were predictive of future gain in visceral fat upon smoking cessation.

Conclusions

Our results provide the first comprehensive characterization of coordinated DNA methylation and gene expression markers of smoking in adipose tissue. The findings relate to human metabolic health and give insights into understanding the widespread health consequence of smoking outside of the lung.

Electronic supplementary material

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

Tissue-specific DNA methylation loss during ageing and carcinogenesis is linked to chromosome structure, replication timing and cell division rates

DNA methylation is an epigenetic mechanism known to affect gene expression and aberrant DNA methylation patterns have been described in cancer. However, only a small fraction of differential methylation events target genes with a defined role in cancer, raising the question of how aberrant DNA methylation contributes to carcinogenesis. As recently a link has been suggested between methylation patterns arising in ageing and those arising in cancer, we asked which aberrations are unique to cancer and which are the product of normal ageing processes. We therefore compared the methylation patterns between ageing and cancer in multiple tissues. We observed that hypermethylation preferentially occurs in regulatory elements, while hypomethylation is associated with structural features of the chromatin. Specifically, we observed consistent hypomethylation of late-replicating, lamina-associated domains. The extent of hypomethylation was stronger in cancer, but in both ageing and cancer it was proportional to the replication timing of the region and the cell division rate of the tissue. Moreover, cancer patients who displayed more hypomethylation in late-replicating, lamina-associated domains had higher expression of cell division genes. These findings suggest that different cell division rates contribute to tissue- and cancer type-specific DNA methylation profiles.