Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Histone deacetylases and inhibitors in diabetes mellitus and its complications

Diabetes mellitus (DM) is a metabolic disorder characterized by hyperglycemia, with its prevalence linked to both genetic predisposition and environmental factors. Epigenetic modifications, particularly through histone deacetylases (HDACs), have been recognized for their significant influence on DM pathogenesis. This review focuses on the classification of HDACs, their role in DM and its complications, and the potential therapeutic applications of HDAC inhibitors. HDACs, which modulate gene expression without altering DNA sequences, are categorized into four classes with distinct functions and tissue specificity. HDAC inhibitors (HDACi) have shown efficacy in various diseases, including DM, by targeting these enzymes. The review highlights how HDACs regulate β-cell function, insulin sensitivity, and hepatic gluconeogenesis in DM, as well as their impact on diabetic cardiomyopathy, nephropathy, and retinopathy. Finally, we suggest that targeted histone modification is expected to become a key method for the treatment of diabetes and its complications. The study of HDACi offers insights into new treatment strategies for DM and its associated complications.

Association of ZBTB38 gene polymorphism (rs724016) with height and fetal hemoglobin in individuals with sickle cell anemia

Objectives

Our study evaluated the association of the polymorphism rs724016 in the ZBTB38 gene, previously associated with height in other populations, with predictors of height, clinical outcomes, and laboratory parameters in sickle cell anemia (SCA).

Methods

Cross-sectional study with individuals with SCA and aged between 3 and 20 years. Clinical, laboratory, molecular, and bone age (BA) data were evaluated. Levels of IGF-1 and IGFBP-3 were adjusted for BA, target height (TH) was calculated as the mean parental height standard deviation score (SDS), and predicted adult height (PAH) SDS was calculated using BA.

Results

We evaluated 80 individuals with SCA. The homozygous genotype of the G allele of rs724016 was associated with a lower height SDS (p < 0.001) and, in a additive genetic model, was negatively associated with HbF levels (p = 0.016). Lower adjusted IGF-1 levels were associated with co-inheritance of alpha-thalassemia and with the absence of HU therapy. Elevated HbF levels were associated with a lower deficit in adjusted growth potential (TH minus PAH).

Conclusion

Our analysis shows that SNP rs724016 in the ZBTB38 is associated with shorter height and lower HbF levels, an important modifier of SCA.

Systematic epigenome editing captures the context-dependent instructive function of chromatin modifications

Chromatin modifications are linked with regulating patterns of gene expression, but their causal role and context-dependent impact on transcription remains unresolved. Here we develop a modular epigenome editing platform that programs nine key chromatin modifications, or combinations thereof, to precise loci in living cells. We couple this with single-cell readouts to systematically quantitate the magnitude and heterogeneity of transcriptional responses elicited by each specific chromatin modification. Among these, we show that installing histone H3 lysine 4 trimethylation (H3K4me3) at promoters can causally instruct transcription by hierarchically remodeling the chromatin landscape. We further dissect how DNA sequence motifs influence the transcriptional impact of chromatin marks, identifying switch-like and attenuative effects within distinct cis contexts. Finally, we examine the interplay of combinatorial modifications, revealing that co-targeted H3K27 trimethylation (H3K27me3) and H2AK119 monoubiquitination (H2AK119ub) maximizes silencing penetrance across single cells. Our precision-perturbation strategy unveils the causal principles of how chromatin modification(s) influence transcription and dissects how quantitative responses are calibrated by contextual interactions.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men

BACKGROUND

American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity.

METHODS

To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays.

RESULTS

Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis.

CONCLUSIONS

These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue.

IMPACT

Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.

Is developmental plasticity triggered by DNA methylation changes in the invasive cane toad (Rhinella marina)?

Many organisms can adjust their development according to environmental conditions, including the presence of conspecifics. Although this developmental plasticity is common in amphibians, its underlying molecular mechanisms remain largely unknown. Exposure during development to either 'cannibal cues' from older conspecifics, or 'alarm cues' from injured conspecifics, causes reduced growth and survival in cane toad (Rhinella marina) tadpoles. Epigenetic modifications, such as changes in DNA methylation patterns, are a plausible mechanism underlying these developmental plastic responses. Here we tested this hypothesis, and asked whether cannibal cues and alarm cues trigger the same DNA methylation changes in developing cane toads. We found that exposure to both cannibal cues and alarm cues was associated with local changes in DNA methylation patterns. These DNA methylation changes affected genes putatively involved in developmental processes, but in different genomic regions for different conspecific-derived cues. Genetic background explains most of the epigenetic variation among individuals. Overall, the molecular mechanisms triggered by exposure to cannibal cues seem to differ from those triggered by alarm cues. Studies linking epigenetic modifications to transcriptional activity are needed to clarify the proximate mechanisms that regulate developmental plasticity in cane toads.

Sex-dependent placental methylation quantitative trait loci provide insight into the prenatal origins of childhood onset traits and conditions

Molecular quantitative trait loci (QTLs) allow us to understand the biology captured in genome-wide association studies (GWASs). The placenta regulates fetal development and shows sex differences in DNA methylation. We therefore hypothesized that placental methylation QTL (mQTL) explain variation in genetic risk for childhood onset traits, and that effects differ by sex. We analyzed 411 term placentas from two studies and found 49,252 methylation (CpG) sites with mQTL and 2,489 CpG sites with sex-dependent mQTL. All mQTL were enriched in regions that typically affect gene expression in prenatal tissues. All mQTL were also enriched in GWAS results for growth- and immune-related traits, but male- and female-specific mQTL were more enriched than cross-sex mQTL. mQTL colocalized with trait loci at 777 CpG sites, with 216 (28%) specific to males or females. Overall, mQTL specific to male and female placenta capture otherwise overlooked variation in childhood traits.

Pan-cancer analysis of promoter activity quantitative trait loci

Altered promoter activity has been generally observed in diverse biological processes, including tumorigenesis. Accumulating evidence suggests that employing a quantitative trait locus mapping approach is effective in comprehending the genetic basis of promoter activity. By utilizing genotype data from The Cancer Genome Atlas and calculating corresponding promoter activity values using proActiv, we systematically evaluated the impact of genetic variants on promoter activity and identified >1.0 million promoter activity quantitative trait loci (paQTLs) as both cis- and trans-acting. Additionally, leveraging data from the genome-wide association study (GWAS) catalog, we discovered >1.3 million paQTLs that overlap with known GWAS linkage disequilibrium regions. Remarkably, ∼9324 paQTLs exhibited significant associations with patient prognosis. Moreover, investigating the impact of promoter activity on >1000 imputed antitumor therapy responses among pan-cancer patients revealed >43 000 million significant associations. Furthermore, ∼25 000 significant associations were identified between promoter activity and immune cell abundance. Finally, a user-friendly data portal, Pancan-paQTL (https://www.hbpding.com/PancanPaQTL/), was constructed for users to browse, search and download data of interest. Pancan-paQTL serves as a comprehensive multidimensional database, enabling functional and clinical investigations into genetic variants associated with promoter activity, drug responses and immune infiltration across multiple cancer types.

Endoplasmic reticulum stress related genome-wide Mendelian randomization identifies therapeutic genes for ulcerative colitis and Crohn's disease

Background: Endoplasmic reticulum stress (ERS) is an important pathophysiological mechanism in ulcerative colitis (UC) and Crohn's disease (CD). ERS-related genes may be influenced by genetic factors and intestinal inflammation. However, the role of ERS as a trigger or potential etiological factor for UC and CD is unclear, as the expression of ERS-related genes in UC and CD may be the cause or subsequent changes in intestinal inflammation. Here, we used a three-step summary data-based Mendelian randomization (SMR) approach integrating multi-omics data to identify putative causal effects of ERS-related genes in UC and CD. Methods: Genome-wide association study (GWAS) summary data for UC (6,968 cases and 20,464 controls) and CD (5,956 cases and 14,927 controls) were extracted as outcome, and DNA methylation quantitative trait loci (mQTL, 1,980 participants) data and expression QTL data (eQTL, 31,684 participants) from the blood were obtained as exposure. The ERS-related genes were extracted from the GeneCards database, and then the GWAS summary data were integrated with the mQTL and eQTL data associated with ERS genes by SMR. Sensitivity analysis included two-sample MR analysis, power calculations, Bayesian co-localization analysis, and phenotype scanning were performed to evaluate the robustness of the results. Results: A total of 1,193 ERS-related genes were obtained. The three-step SMR analysis showed that cg24011261 CpG site regulating GPX1 expression was associated with a low risk of UC, whereas GPX1 expression regulated by a combination of cg05055782, cg24011261, and cg05551922 CpG sites was associated with a low risk of CD. Sensitivity analysis further supports these findings. Conclusion: This multi-omics integration study identifies a causal relationship between the role of ERS in UC and CD and suggests potential new therapeutic targets for clinical practice.

Simultaneous sequencing of genetic and epigenetic bases in DNA

DNA comprises molecular information stored in genetic and epigenetic bases, both of which are vital to our understanding of biology. Most DNA sequencing approaches address either genetics or epigenetics and thus capture incomplete information. Methods widely used to detect epigenetic DNA bases fail to capture common C-to-T mutations or distinguish 5-methylcytosine from 5-hydroxymethylcytosine. We present a single base-resolution sequencing methodology that sequences complete genetics and the two most common cytosine modifications in a single workflow. DNA is copied and bases are enzymatically converted. Coupled decoding of bases across the original and copy strand provides a phased digital readout. Methods are demonstrated on human genomic DNA and cell-free DNA from a blood sample of a patient with cancer. The approach is accurate, requires low DNA input and has a simple workflow and analysis pipeline. Simultaneous, phased reading of genetic and epigenetic bases provides a more complete picture of the information stored in genomes and has applications throughout biomedicine.

Exploring molecular markers and drug candidates for colorectal cancer through comprehensive bioinformatics analysis

Colorectal cancer (CRC) often has a poor prognosis and identifying useful and novel agents for treating CRC is urgently required. This study aimed to examine molecular markers associated with CRC prognosis and to identify potential drug candidates. The differentially expressed genes (DEGs) of CRC in TCGA were identified. The genes associated with CRC, summarized from NCBI-gene, OMIM, and the DEGs, were used to construct a co-expression network by WGCNA. Moreover, the co-expression genes from modules of interest were used to carry out functional enrichment. A total of 2742 DEGs, including 1674 upregulated and 1068 downregulated genes, were identified. Thirteen co-expression modules were constructed with WGCNA. Brown and blue co-expression modules with significant differences in disease phenotype were found. Functional enrichment analysis showed that genes in the brown module were mainly related to cell cycle, cell proliferation, DNA replication, and RNA transport. The genes in the blue module were mainly associated with fatty acid degradation, sulfur metabolism, PPAR signaling pathway and bile secretion. In addition, both the genes in brown and blue were associated with tumor staging. Some prognostic markers and candidate small molecules drugs for CRC treatment were identified. In conclusion, we revealed molecular biomarker profiles in CRC by systematic bioinformatics analysis, constructed regulatory networks of mRNA, ncRNA and transcriptional regulators (TFs), and identified potential drugs targeting hub proteins and TFs.

New genetic and epigenetic insights into the chemokine system: the latest discoveries aiding progression toward precision medicine

Over the past thirty years, the importance of chemokines and their seven-transmembrane G protein-coupled receptors (GPCRs) has been increasingly recognized. Chemokine interactions with receptors trigger signaling pathway activity to form a network fundamental to diverse immune processes, including host homeostasis and responses to disease. Genetic and nongenetic regulation of both the expression and structure of chemokines and receptors conveys chemokine functional heterogeneity. Imbalances and defects in the system contribute to the pathogenesis of a variety of diseases, including cancer, immune and inflammatory diseases, and metabolic and neurological disorders, which render the system a focus of studies aiming to discover therapies and important biomarkers. The integrated view of chemokine biology underpinning divergence and plasticity has provided insights into immune dysfunction in disease states, including, among others, coronavirus disease 2019 (COVID-19). In this review, by reporting the latest advances in chemokine biology and results from analyses of a plethora of sequencing-based datasets, we outline recent advances in the understanding of the genetic variations and nongenetic heterogeneity of chemokines and receptors and provide an updated view of their contribution to the pathophysiological network, focusing on chemokine-mediated inflammation and cancer. Clarification of the molecular basis of dynamic chemokine-receptor interactions will help advance the understanding of chemokine biology to achieve precision medicine application in the clinic.

Local CpG density affects the trajectory and variance of age-associated DNA methylation changes

BACKGROUND

DNA methylation is an epigenetic mark associated with the repression of gene promoters. Its pattern in the genome is disrupted with age and these changes can be used to statistically predict age with epigenetic clocks. Altered rates of aging inferred from these clocks are observed in human disease. However, the molecular mechanisms underpinning age-associated DNA methylation changes remain unknown. Local DNA sequence can program steady-state DNA methylation levels, but how it influences age-associated methylation changes is unknown.

RESULTS

We analyze longitudinal human DNA methylation trajectories at 345,895 CpGs from 600 individuals aged between 67 and 80 to understand the factors responsible for age-associated epigenetic changes at individual CpGs. We show that changes in methylation with age occur at 182,760 loci largely independently of variation in cell type proportions. These changes are especially apparent at 8322 low CpG density loci. Using SNP data from the same individuals, we demonstrate that methylation trajectories are affected by local sequence polymorphisms at 1487 low CpG density loci. More generally, we find that low CpG density regions are particularly prone to change and do so variably between individuals in people aged over 65. This differs from the behavior of these regions in younger individuals where they predominantly lose methylation.

CONCLUSIONS

Our results, which we reproduce in two independent groups of individuals, demonstrate that local DNA sequence influences age-associated DNA methylation changes in humans in vivo. We suggest that this occurs because interactions between CpGs reinforce maintenance of methylation patterns in CpG dense regions.

Kidney omics in hypertension: from statistical associations to biological mechanisms and clinical applications

Hypertension is a major cardiovascular disease risk factor and contributor to premature death globally. Family-based investigations confirmed a significant heritable component of blood pressure (BP), whereas genome-wide association studies revealed >1000 common and rare genetic variants associated with BP and/or hypertension. The kidney is not only an organ of key relevance to BP regulation and the development of hypertension, but it also acts as the tissue mediator of genetic predisposition to hypertension. The identity of kidney genes, pathways, and related mechanisms underlying the genetic associations with BP has started to emerge through integration of genomics with kidney transcriptomics, epigenomics, and other omics as well as through applications of causal inference, such as Mendelian randomization. Single-cell methods further enabled mapping of BP-associated kidney genes to cell types, and in conjunction with other omics, started to illuminate the biological mechanisms underpinning associations of BP-associated genetic variants and kidney genes. Polygenic risk scores derived from genome-wide association studies and refined on kidney omics hold the promise of enhanced diagnostic prediction, whereas kidney omics-informed drug discovery is likely to contribute new therapeutic opportunities for hypertension and hypertension-mediated kidney damage.

Increased 5-hydroxymethylcytosine is a favorable prognostic factor of Helicobacter pylori-negative gastric cancer patients

BACKGROUND

Most gastric cancer (GC) patients are diagnosed at middle or late stage because the symptoms in early stage are obscure, which causes higher mortality rates of GC. Helicobacter pylori (H. pylori) was identified as a class I carcinogen and leads to aberrant DNA methylation/hydroxymethylation. 5-hydroxymethylcytosine (5-hmC) plays complex roles in gene regulation of tumorigenesis and can be considered as an activating epigenetic mark of hydroxymethylation.

AIM

To explore the association between 5-hmC levels and the progression and prognosis of GC patients with or without H. pylori infection.

METHODS

A retrospective cohort study was conducted to estimate the predicted value of 5-hmC level in the progression and prognosis of GC patients with different H. pylori infection status. A total of 144 GC patients were recruited.

RESULTS

The levels of 5-hmC were significantly decreased in tumor tissues (0.076 ± 0.048) compared with the matched control tissues (0.110 ± 0.057, P = 0.001). A high level of 5-hmC was an independent significant favorable predictor of overall survival in GC patients (hazard ratio = 0.61, 95% confidence interval: 0.38-0.98, P = 0.040), the H. pylori-negative GC subgroup (hazard ratio = 0.30, 95% confidence interval: 0.13-0.68, P = 0.004) and the GC patients with TNM stage Ⅰ or Ⅱ (hazard ratio = 0.32, 95% confidence interval: 0.13-0.77, P = 0.011).

CONCLUSION

Increased 5-hmC is a favorable prognostic factor in GC, especially for H. pylori-negative subgroups.

PTBP2 - a gene with relevance for both Anorexia nervosa and body weight regulation

Genetic factors are relevant for both eating disorders and body weight regulation. A recent genome-wide association study (GWAS) for anorexia nervosa (AN) detected eight genome-wide significant chromosomal loci. One of these loci, rs10747478, was also genome-wide and significantly associated with body mass index (BMI). The nearest coding gene is the Polypyrimidine Tract Binding Protein 2 gene (PTBP2). To detect mutations in PTBP2, Sanger sequencing of the coding region was performed in 192 female patients with AN (acute or recovered) and 191 children or adolescents with (extreme) obesity. Twenty-five variants were identified. Twenty-three of these were predicted to be pathogenic or functionally relevant in at least one in silico tool. Two novel synonymous variants (p.Ala77Ala and p.Asp195Asp), one intronic SNP (rs188987764), and the intronic deletion (rs561340981) located in the highly conserved region of PTBP2 may have functional consequences. Ten of 20 genes interacting with PTBP2 were studied for their impact on body weight regulation based on either previous functional studies or GWAS hits for body weight or BMI. In a GWAS for BMI (Pulit et al. 2018), the number of genome-wide significant associations at the PTBP2 locus was different between males (60 variants) and females (two variants, one of these also significant in males). More than 65% of these 61 variants showed differences in the effect size pertaining to BMI between sexes (absolute value of Z-score >2, two-sided p < 0.05). One LD block overlapping 5'UTR and all coding regions of PTBP2 comprises 56 significant variants in males. The analysis based on sex-stratified BMI GWAS summary statistics implies that PTBP2 may have a more pronounced effect on body weight regulation in males than in females.

ezQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci

Genome-wide association studies (GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common trait-associated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus (eQTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loci through statistical colocalization methods. While QTL colocalization is becoming a standard analysis in post-GWAS investigation, an easy web tool for users to perform formal colocalization analyses with either user-provided or public GWAS and eQTL datasets has been lacking. Here, we present ezQTL, a web-based bioinformatic application to interactively visualize and analyze genetic association data such as GWAS loci and molecular QTLs under different linkage disequilibrium (LD) patterns (1000 Genomes Project, UK Biobank, or user-provided data). This application allows users to perform data quality control for variants matched between different datasets, LD visualization, and two-trait colocalization analyses using two state-of-the-art methodologies (eCAVIAR and HyPrColoc), including batch processing. ezQTL is a free and publicly available cross-platform web tool, which can be accessed online at https://analysistools.cancer.gov/ezqtl.

Lysosomal dysfunction in neurodegeneration: emerging concepts and methods

The understanding of lysosomes has come a long way since the initial discovery of their role in degrading cellular waste. The lysosome is now recognized as a highly dynamic organelle positioned at the crossroads of cell signaling, transcription, and metabolism. Underscoring its importance is the observation that, in addition to rare monogenic lysosomal storage disorders, genes regulating lysosomal function are implicated in common sporadic neurodegenerative diseases, including Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS). Developing therapies for these disorders is particularly challenging, largely due to gaps in knowledge of the underlying molecular and cellular processes. In this review, we discuss technological advances that have propelled deeper understanding of the lysosome in neurodegeneration, from elucidating the functions of lysosome-related disease risk variants at the level of the organelle, cell, and tissue, to the development of disease-specific biological models that recapitulate disease manifestations. Finally, we identify key questions to be addressed to successfully bridge the gap to the clinic.

Smoking during pregnancy is associated with child overweight independent of maternal pre-pregnancy BMI and genetic predisposition to adiposity

High maternal body mass index (BMI) and smoking during pregnancy are risk factors for child overweight. Maternal smoking tends to reduce her BMI and the association of smoking with child overweight may be confounded by or interacting with maternal genetic predisposition to adiposity. In the Danish National Birth Cohort, we investigated whether smoking during pregnancy is associated with child BMI/overweight independent of pre-pregnancy BMI and maternal genetic predisposition to adiposity estimated as total, transmitted and non-transmitted genetic risk scores (GRSs) based on 941 common genetic variants associated with BMI. Smoking during pregnancy was associated with higher child BMI and higher odds of child overweight in a dose–response relationship. The odds ratio (95% CI) for smoking 11 + cigarettes in third trimester versus no smoking was 2.42 (1.30; 4.50), irrespective of maternal BMI and maternal GRSs (total, transmitted or non-transmitted). There were no statistically significant interactions between maternal GRSs and smoking (all p-values for interactions > 0.05). In conclusion, in this study, smoking during pregnancy exhibits a dose–response association with increased child BMI/overweight, independent of maternal pre-pregnancy BMI, maternal transmitted, and non-transmitted genetic predisposition to adiposity. Avoidance of smoking during pregnancy may help prevent childhood obesity irrespective of the mother–child genetic predisposition.

Environmentally sensitive hotspots in the methylome of the early human embryo

In humans, DNA methylation marks inherited from gametes are largely erased following fertilisation, prior to construction of the embryonic methylome. Exploiting a natural experiment of seasonal variation including changes in diet and nutritional status in rural Gambia, we analysed three datasets covering two independent child cohorts and identified 259 CpGs showing consistent associations between season of conception (SoC) and DNA methylation. SoC effects were most apparent in early infancy, with evidence of attenuation by mid-childhood. SoC-associated CpGs were enriched for metastable epialleles, parent-of-origin-specific methylation and germline differentially methylated regions, supporting a periconceptional environmental influence. Many SoC-associated CpGs overlapped enhancers or sites of active transcription in H1 embryonic stem cells and fetal tissues. Half were influenced but not determined by measured genetic variants that were independent of SoC. Environmental 'hotspots' providing a record of environmental influence at periconception constitute a valuable resource for investigating epigenetic mechanisms linking early exposures to lifelong health and disease.

Cis-regulatory hubs: a new 3D model of complex disease genetics with an application to schizophrenia

Genes and their regulatory elements are organized in neurons within 3D networks which model functional structures and explain schizophrenia genetic etiology.

The 3D conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases such as schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs, and schizophrenia heritability compared with equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared with promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Identifying causal models between genetically regulated methylation patterns and gene expression in healthy colon tissue

BACKGROUND

DNA methylation is involved in the regulation of gene expression and phenotypic variation, but the inter-relationship between genetic variation, DNA methylation and gene expression remains poorly understood. Here we combine the analysis of genetic variants related to methylation markers (methylation quantitative trait loci: mQTLs) and gene expression (expression quantitative trait loci: eQTLs) with methylation markers related to gene expression (expression quantitative trait methylation: eQTMs), to provide novel insights into the genetic/epigenetic architecture of colocalizing molecular markers.

RESULTS

Normal mucosa from 100 patients with colon cancer and 50 healthy donors included in the Colonomics project have been analyzed. Linear models have been used to find mQTLs and eQTMs within 1 Mb of the target gene. From 32,446 eQTLs previously detected, we found a total of 6850 SNPs, 114 CpGs and 52 genes interrelated, generating 13,987 significant combinations of co-occurring associations (meQTLs) after Bonferromi correction. Non-redundant meQTLs were 54, enriched in genes involved in metabolism of glucose and xenobiotics and immune system. SNPs in meQTLs were enriched in regulatory elements (enhancers and promoters) compared to random SNPs within 1 Mb of genes. Three colorectal cancer GWAS SNPs were related to methylation changes, and four SNPs were related to chemerin levels. Bayesian networks have been used to identify putative causal relationships among associated SNPs, CpG and gene expression triads. We identified that most of these combinations showed the canonical pathway of methylation markers causes gene expression variation (60.1%) or non-causal relationship between methylation and gene expression (33.9%); however, in up to 6% of these combinations, gene expression was causing variation in methylation markers.

CONCLUSIONS

In this study we provided a characterization of the regulation between genetic variants and inter-dependent methylation markers and gene expression in a set of 150 healthy colon tissue samples. This is an important finding for the understanding of molecular susceptibility on colon-related complex diseases.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Mapping methylation quantitative trait loci in cardiac tissues nominates risk loci and biological pathways in congenital heart disease

Background

Most congenital heart defects (CHDs) result from complex interactions among genetic susceptibilities, epigenetic modifications, and maternal environmental exposures. Characterizing the complex relationship between genetic, epigenetic, and transcriptomic variation will enhance our understanding of pathogenesis in this important type of congenital disorder. We investigated cis-acting effects of genetic single nucleotide polymorphisms (SNPs) on local DNA methylation patterns within 83 cardiac tissue samples and prioritized their contributions to CHD risk by leveraging results of CHD genome-wide association studies (GWAS) and their effects on cardiac gene expression.

Results

We identified 13,901 potential methylation quantitative trait loci (mQTLs) with a false discovery threshold of 5%. Further co-localization analyses and Mendelian randomization indicated that genetic variants near the HLA-DRB6 gene on chromosome 6 may contribute to CHD risk by regulating the methylation status of nearby CpG sites. Additional SNPs in genomic regions on chromosome 10 (TNKS2-AS1 gene) and chromosome 14 (LINC01629 gene) may simultaneously influence epigenetic and transcriptomic variations within cardiac tissues.

Conclusions

Our results support the hypothesis that genetic variants may influence the risk of CHDs through regulating the changes of DNA methylation and gene expression. Our results can serve as an important source of information that can be integrated with other genetic studies of heart diseases, especially CHDs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00975-2.

Downregulation of 5-hydroxymethylcytosine is an early event in pancreatic tumorigenesis

Epigenetic alterations are increasingly recognized as important contributors to the development and progression of pancreatic ductal adenocarcinoma. 5-hydroxymethylcytosine (5hmC) is an epigenetic DNA mark generated through the ten-eleven translocation (TET) enzyme-mediated pathway and is closely linked to gene activation. However, the timing of alterations in epigenetic regulation in the progression of pancreatic neoplasia is not well understood. In this study, we hypothesized that aberrant expression of ten-eleven translocation methylcytosine dioxygenase 1 (TET1) and subsequent global 5hmC alteration are linked to early tumorigenesis in the pancreas. Therefore, we evaluated alterations of 5hmC and TET1 levels using immunohistochemistry in pancreatic neoplasms (n = 380) and normal ducts (n = 118). The study cohort included representation of the full spectrum of precancerous lesions from low- and high-grade pancreatic intraepithelial neoplasia (n = 95), intraductal papillary mucinous neoplasms (all subtypes, n = 129), intraductal oncocytic papillary neoplasms (n = 12), and mucinous cystic neoplasms (n = 144). 5hmC and TET1 were significantly downregulated in all types of precancerous lesion and associated invasive pancreatic ductal adenocarcinomas compared with normal ductal epithelium (all p < 0.001), and expression of 5hmC positively correlated with expression of TET1. Importantly, downregulation of both 5hmC and TET1 was observed in most low-grade precancerous lesions. There were no clear associations between 5hmC levels and clinicopathological factors, thereby suggesting a common epigenetic abnormality across precancerous lesions. We conclude that downregulation of 5hmC and TET1 is an early event in pancreatic tumorigenesis. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Uncovering genetic mechanisms of hypertension through multi-omic analysis of the kidney

The kidney is an organ of key relevance to blood pressure (BP) regulation, hypertension and antihypertensive treatment. However, genetically mediated renal mechanisms underlying susceptibility to hypertension remain poorly understood. We integrated genotype, gene expression, alternative splicing and DNA methylation profiles of up to 430 human kidneys to characterize the effects of BP index variants from genome-wide association studies (GWASs) on renal transcriptome and epigenome. We uncovered kidney targets for 479 (58.3%) BP-GWAS variants and paired 49 BP-GWAS kidney genes with 210 licensed drugs. Our colocalization and Mendelian randomization analyses identified 179 unique kidney genes with evidence of putatively causal effects on BP. Through Mendelian randomization, we also uncovered effects of BP on renal outcomes commonly affecting patients with hypertension. Collectively, our studies identified genetic variants, kidney genes, molecular mechanisms and biological pathways of key relevance to the genetic regulation of BP and inherited susceptibility to hypertension.

Genetic impacts on DNA methylation: research findings and future perspectives

Multiple recent studies highlight that genetic variants can have strong impacts on a significant proportion of the human DNA methylome. Methylation quantitative trait loci, or meQTLs, allow for the exploration of biological mechanisms that underlie complex human phenotypes, with potential insights for human disease onset and progression. In this review, we summarize recent milestones in characterizing the human genetic basis of DNA methylation variation over the last decade, including heritability findings and genome-wide identification of meQTLs. We also discuss challenges in this field and future areas of research geared to generate insights into molecular processes underlying human complex traits.

Interpretation of allele-specific chromatin accessibility using cell state-aware deep learning

Genomic sequence variation within enhancers and promoters can have a significant impact on the cellular state and phenotype. However, sifting through the millions of candidate variants in a personal genome or a cancer genome, to identify those that impact cis-regulatory function, remains a major challenge. Interpretation of noncoding genome variation benefits from explainable artificial intelligence to predict and interpret the impact of a mutation on gene regulation. Here we generate phased whole genomes with matched chromatin accessibility, histone modifications, and gene expression for 10 melanoma cell lines. We find that training a specialized deep learning model, called DeepMEL2, on melanoma chromatin accessibility data can capture the various regulatory programs of the melanocytic and mesenchymal-like melanoma cell states. This model outperforms motif-based variant scoring, as well as more generic deep learning models. We detect hundreds to thousands of allele-specific chromatin accessibility variants (ASCAVs) in each melanoma genome, of which 15%-20% can be explained by gains or losses of transcription factor binding sites. A considerable fraction of ASCAVs are caused by changes in AP-1 binding, as confirmed by matched ChIP-seq data to identify allele-specific binding of JUN and FOSL1. Finally, by augmenting the DeepMEL2 model with ChIP-seq data for GABPA, the TERT promoter mutation, as well as additional ETS motif gains, can be identified with high confidence. In conclusion, we present a new integrative genomics approach and a deep learning model to identify and interpret functional enhancer mutations with allelic imbalance of chromatin accessibility and gene expression.

Methylation quantitative trait locus analysis of chronic postsurgical pain uncovers epigenetic mediators of genetic risk

Background: Overlap of pathways enriched by single nucleotide polymorphisms and DNA-methylation underlying chronic postsurgical pain (CPSP), prompted pilot study of CPSP-associated methylation quantitative trait loci (meQTL). Materials & methods: Children undergoing spine-fusion were recruited prospectively. Logistic-regression for genome- and epigenome-wide CPSP association and DNA-methylation-single nucleotide polymorphism association/mediation analyses to identify meQTLs were followed by functional genomics analyses. Results: CPSP (n = 20/58) and non-CPSP groups differed in pain-measures. Of 2753 meQTLs, DNA-methylation at 127 cytosine-guanine dinucleotides mediated association of 470 meQTLs with CPSP (p < 0.05). At PARK16 locus, CPSP risk meQTLs were associated with decreased DNA-methylation at RAB7L1 and increased DNA-methylation at PM20D1. Corresponding RAB7L1/PM20D1 blood eQTLs (GTEx) and cytosine-guanine dinucleotide-loci enrichment for histone marks, transcription factor binding sites and ATAC-seq peaks suggest altered transcription factor-binding. Conclusion: CPSP-associated meQTLs indicate epigenetic mechanisms mediate genetic risk. Clinical trial registration: NCT01839461, NCT01731873 (ClinicalTrials.gov).

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1119 cases (including 589 cocaine-dependent subjects) and 1092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

Chromatin architecture reveals cell type-specific target genes for kidney disease risk variants

BACKGROUND

Cell type-specific transcriptional programming results from the combinatorial interplay between the repertoire of active regulatory elements. Disease-associated variants disrupt such programming, leading to altered expression of downstream regulated genes and the onset of pathological states. However, due to the non-linear regulatory properties of non-coding elements such as enhancers, which can activate transcription at long distances and in a non-directional way, the identification of causal variants and their target genes remains challenging. Here, we provide a multi-omics analysis to identify regulatory elements associated with functional kidney disease variants, and downstream regulated genes.

RESULTS

In order to understand the genetic risk of kidney diseases, we generated a comprehensive dataset of the chromatin landscape of human kidney tubule cells, including transcription-centered 3D chromatin organization, histone modifications distribution and transcriptome with HiChIP, ChIP-seq and RNA-seq. We identified genome-wide functional elements and thousands of interactions between the distal elements and target genes. The results revealed that risk variants for renal tumor and chronic kidney disease were enriched in kidney tubule cells. We further pinpointed the target genes for the variants and validated two target genes by CRISPR/Cas9 genome editing techniques in zebrafish, demonstrating that SLC34A1 and MTX1 were indispensable genes to maintain kidney function.

CONCLUSIONS

Our results provide a valuable multi-omics resource on the chromatin landscape of human kidney tubule cells and establish a bioinformatic pipeline in dissecting functions of kidney disease-associated variants based on cell type-specific epigenome.

Neurogenetic and Neuroepigenetic Mechanisms in Cognitive Health and Disease

Over the last two decades, the explosion of experimental, computational, and high-throughput technologies has led to critical insights into how the brain functions in health and disease. It has become increasingly clear that the vast majority of brain activities result from the complex entanglement of genetic factors, epigenetic changes, and environmental stimuli, which, when altered, can lead to neurodegenerative and neuropsychiatric disorders. Nevertheless, a complete understanding of the molecular mechanisms underlying neuronal activities and higher-order cognitive processes continues to elude neuroscientists. Here, we provide a concise overview of how the interaction between the environment and genetic as well as epigenetic mechanisms shapes complex neuronal processes such as learning, memory, and synaptic plasticity. We then consider how this interaction contributes to the development of neurodegenerative and psychiatric disorders, and how it can be modeled to predict phenotypic variability and disease risk. Finally, we outline new frontiers in neurogenetic and neuroepigenetic research and highlight the challenges these fields will face in their quest to decipher the molecular mechanisms governing brain functioning.

Comparison of DNA Methylation Profiles of Hemostatic Genes between Liver Tissue and Peripheral Blood within Individuals

DNA methylation has become increasingly recognized in the etiology of complex diseases, including thrombotic disorders. Blood is often collected in epidemiological studies for genotyping and has recently also been used to examine DNA methylation in epigenome-wide association studies. DNA methylation patterns are often tissue-specific, thus, peripheral blood may not accurately reflect the methylation pattern in the tissue of relevance. Here, we collected paired liver and blood samples concurrently from 27 individuals undergoing liver surgery. We performed targeted bisulfite sequencing for a set of 35 hemostatic genes primarily expressed in liver to analyze DNA methylation levels of >10,000 cytosine-phosphate-guanine (CpG) dinucleotides. We evaluated whether DNA methylation in blood could serve as a proxy for DNA methylation in liver at individual CpGs. Approximately 30% of CpGs were nonvariable and were predominantly hypo- (<25%) or hypermethylated (>70%) in both tissues. While blood can serve as a proxy for liver at these CpGs, the low variability renders these unlikely to explain phenotypic differences. We therefore focused on CpG sites with variable methylation levels in liver. The level of blood-liver tissue correlation varied widely across these variable CpGs; moderate correlations (0.5 ≤ r < 0.75) were detected for 6% and strong correlations (r ≥ 0.75) for a further 4%. Our findings indicate that it is essential to study the concordance of DNA methylation between blood and liver at individual CpGs. This paired blood-liver dataset is intended as a resource to aid interpretation of blood-based DNA methylation results.

Integrative analyses prioritize GNL3 as a risk gene for bipolar disorder

Genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms (SNPs) associated with bipolar disorder (BD), but what the causal variants are and how they contribute to BD is largely unknown. In this study, we used FUMA, a GWAS annotation tool, to pinpoint potential causal variants and genes from the latest BD GWAS findings, and performed integrative analyses, including brain expression quantitative trait loci (eQTL), gene coexpression network, differential gene expression, protein-protein interaction, and brain intermediate phenotype association analysis to identify the functions of a prioritized gene and its connection to BD. Convergent lines of evidence prioritized protein-coding gene G Protein Nucleolar 3 (GNL3) as a BD risk gene, with integrative analyses revealing GNL3's roles in cell proliferation, neuronal functions, and brain phenotypes. We experimentally revealed that BD-related eQTL SNPs rs10865973, rs12635140, and rs4687644 regulate GNL3 expression using dual luciferase reporter assay and CRISPR interference experiment in human neural progenitor cells. We further identified that GNL3 knockdown and overexpression led to aberrant neuronal proliferation and differentiation, using two-dimensional human neural cell cultures and three-dimensional forebrain organoid model. This study gathers evidence that BD-related genetic variants regulate GNL3 expression which subsequently affects neuronal proliferation and differentiation.

Effect of maternal preconceptional and pregnancy micronutrient interventions on children's DNA methylation: Findings from the EMPHASIS study

BACKGROUND

Maternal nutrition in pregnancy has been linked to offspring health in early and later life, with changes to DNA methylation (DNAm) proposed as a mediating mechanism.

OBJECTIVE

We investigated intervention-associated DNAm changes in children whose mothers participated in 2 randomized controlled trials of micronutrient supplementation before and during pregnancy, as part of the EMPHASIS (Epigenetic Mechanisms linking Preconceptional nutrition and Health Assessed in India and sub-Saharan Africa) study (ISRCTN14266771).

DESIGN

We conducted epigenome-wide association studies with blood samples from Indian (n = 698) and Gambian (n = 293) children using the Illumina EPIC array and a targeted study of selected loci not on the array. The Indian micronutrient intervention was food based, whereas the Gambian intervention was a micronutrient tablet.

RESULTS

We identified 6 differentially methylated CpGs in Gambians [2.5-5.0% reduction in intervention group, all false discovery rate (FDR) <5%], the majority mapping to ESM1, which also represented a strong signal in regional analysis. One CpG passed FDR <5% in the Indian cohort, but overall effect sizes were small (<1%) and did not have the characteristics of a robust signature. We also found strong evidence for enrichment of metastable epialleles among subthreshold signals in the Gambian analysis. This supports the notion that multiple methylation loci are influenced by micronutrient supplementation in the early embryo.

CONCLUSIONS

Maternal preconceptional and pregnancy micronutrient supplementation may alter DNAm in children measured at 7-9 y. Multiple factors, including differences between the nature of the intervention, participants, and settings, are likely to have contributed to the lack of replication in the Indian cohort. Potential links to phenotypic outcomes will be explored in the next stage of the EMPHASIS study.

ReQTL: identifying correlations between expressed SNVs and gene expression using RNA-sequencing data

MOTIVATION

By testing for associations between DNA genotypes and gene expression levels, expression quantitative trait locus (eQTL) analyses have been instrumental in understanding how thousands of single nucleotide variants (SNVs) may affect gene expression. As compared to DNA genotypes, RNA genetic variation represents a phenotypic trait that reflects the actual allele content of the studied system. RNA genetic variation at expressed SNV loci can be estimated using the proportion of alleles bearing the variant nucleotide (variant allele fraction, VAFRNA). VAFRNA is a continuous measure which allows for precise allele quantitation in loci where the RNA alleles do not scale with the genotype count. We describe a method to correlate VAFRNA with gene expression and assess its ability to identify genetically regulated expression solely from RNA-sequencing (RNA-seq) datasets.

RESULTS

We introduce ReQTL, an eQTL modification which substitutes the DNA allele count for the variant allele fraction at expressed SNV loci in the transcriptome (VAFRNA). We exemplify the method on sets of RNA-seq data from human tissues obtained though the Genotype-Tissue Expression (GTEx) project and demonstrate that ReQTL analyses are computationally feasible and can identify a subset of expressed eQTL loci.

AVAILABILITY AND IMPLEMENTATION

A toolkit to perform ReQTL analyses is available at https://github.com/HorvathLab/ReQTL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Enhancer of zeste 2 polycomb repressive complex 2 subunit polymorphisms in melanoma skin cancer risk

Melanoma is the most deadly skin cancer, and its incidence is growing. EZH2, a member of the Polycomb Group (PcGs) proteins family, plays an important biological role in the occurrence and development of melanoma. EZH2 germline genetic polymorphisms have not been yet evaluated in melanoma predisposition. Three hundred thirty sporadic Italian melanoma patients and 333 healthy volunteers were genotyped to analyse the association between EZH2 variants rs6950683, rs2302427, rs3757441, rs2072408 and melanoma risk. The functionality of rs6950683 alleles was investigated in keratinocytes (HaCat), melanoma cells (A375) and human embryonic kidney cells (HEK293), using promoter-reporter assays. Genotype distribution of SNPs showed that rs6950683T and rs3757441C alleles were positively associated with melanoma risk (P = .003 and .004, respectively). Haplotype analysis revealed that TCCA and CCCG haplotypes were associated with a higher risk of melanoma (P = .02 and .04, respectively). Functional assays demonstrated that allele rs6950683T reduce promoter activity in the three cell lines analysed compared to C allele. rs6950683T and rs3757441C alleles in the EZH2 gene appear positively associated with melanoma risk in the analysed population. In addition, we demonstrated for the first time the functional role of rs6950683 upstream polymorphism on EZH2 gene expression regulation.

Epigenetic mechanisms in schizophrenia and other psychotic disorders: a systematic review of empirical human findings

Schizophrenia and other psychotic disorders are highly debilitating psychiatric conditions that lack a clear etiology and exhibit polygenic inheritance underlain by pleiotropic genes. The prevailing explanation points to the interplay between predisposing genes and environmental exposure. Accumulated evidence suggests that epigenetic regulation of the genome may mediate dynamic gene-environment interactions at the molecular level by modulating the expression of psychiatric phenotypes through transcription factors. This systematic review summarizes the current knowledge linking schizophrenia and other psychotic disorders to epigenetics, based on PubMed and Web of Science database searches conducted according to the PRISMA guidelines. Three groups of mechanisms in case-control studies  of human tissue (i.e., postmortem brain and bio-fluids) were considered: DNA methylation, histone modifications, and non-coding miRNAs. From the initial pool of 3,204 records, 152 studies met our inclusion criteria (11,815/11,528, 233/219, and 2,091/1,827 cases/controls for each group, respectively). Many of the findings revealed associations with epigenetic modulations of genes regulating neurotransmission, neurodevelopment, and immune function, as well as differential miRNA expression (e.g., upregulated miR-34a, miR-7, and miR-181b). Overall, actual evidence moderately supports an association between epigenetics and schizophrenia and other psychotic disorders. However, heterogeneous results and cross-tissue extrapolations call for future work. Integrating epigenetics into systems biology may critically enhance research on psychosis and thus our understanding of the disorder. This may have implications for psychiatry in risk stratification, early recognition, diagnostics, precision medicine, and other interventional approaches targeting epigenetic fingerprints.

EpigenCentral: Portal for DNA methylation data analysis and classification in rare diseases

Epigenetic processes play a key role in regulating gene expression. Genetic variants that disrupt chromatin-modifying proteins are associated with a broad range of diseases, some of which have specific epigenetic patterns, such as aberrant DNA methylation (DNAm), which may be used as disease biomarkers. While much of the epigenetic research has focused on cancer, there is a paucity of resources devoted to neurodevelopmental disorders (NDDs), which include autism spectrum disorder and many rare, clinically overlapping syndromes. To address this challenge, we created EpigenCentral, a free web resource for biomedical researchers, molecular diagnostic laboratories, and clinical practitioners to perform the interactive classification and analysis of DNAm data related to NDDs. It allows users to search for known disease-associated patterns in their DNAm data, classify genetic variants as pathogenic or benign to assist in molecular diagnostics, or analyze patterns of differential methylation in their data through a simple web form. EpigenCentral is freely available at http://epigen.ccm.sickkids.ca/.

Allele-specific DNA methylation is increased in cancers and its dense mapping in normal plus neoplastic cells increases the yield of disease-associated regulatory SNPs

BACKGROUND

Mapping of allele-specific DNA methylation (ASM) can be a post-GWAS strategy for localizing regulatory sequence polymorphisms (rSNPs). The advantages of this approach, and the mechanisms underlying ASM in normal and neoplastic cells, remain to be clarified.

RESULTS

We perform whole genome methyl-seq on diverse normal cells and tissues and three cancer types. After excluding imprinting, the data pinpoint 15,112 high-confidence ASM differentially methylated regions, of which 1838 contain SNPs in strong linkage disequilibrium or coinciding with GWAS peaks. ASM frequencies are increased in cancers versus matched normal tissues, due to widespread allele-specific hypomethylation and focal allele-specific hypermethylation in poised chromatin. Cancer cells show increased allele switching at ASM loci, but disruptive SNPs in specific classes of CTCF and transcription factor binding motifs are similarly correlated with ASM in cancer and non-cancer. Rare somatic mutations affecting these same motif classes track with de novo ASM. Allele-specific transcription factor binding from ChIP-seq is enriched among ASM loci, but most ASM differentially methylated regions lack such annotations, and some are found in otherwise uninformative "chromatin deserts."

CONCLUSIONS

ASM is increased in cancers but occurs by a shared mechanism involving disruptive SNPs in CTCF and transcription factor binding sites in both normal and neoplastic cells. Dense ASM mapping in normal plus cancer samples reveals candidate rSNPs that are difficult to find by other approaches. Together with GWAS data, these rSNPs can nominate specific transcriptional pathways in susceptibility to autoimmune, cardiometabolic, neuropsychiatric, and neoplastic diseases.

CloudASM: an ultra-efficient cloud-based pipeline for mapping allele-specific DNA methylation

SUMMARY

Methods for quantifying the imbalance in CpG methylation between alleles genome-wide have been described but their algorithmic time complexity is quadratic and their practical use requires painstaking attention to infrastructure choice, implementation and execution. To solve this problem, we developed CloudASM, a scalable, ultra-efficient, turn-key, portable pipeline on Google Cloud Platform (GCP) that uses a novel pipeline manager and GCP's serverless enterprise data warehouse.

AVAILABILITY AND IMPLEMENTATION

CloudASM is freely available in the GitHub repository https://github.com/TyckoLab/CloudASM and a sample dataset and its results are also freely available at https://console.cloud.google.com/storage/browser/cloudasm.

CONTACT

emmanuel.dumont@hmh-cdi.org.

DNA methylation near the INS gene is associated with INS genetic variation (rs689) and type 1 diabetes in the Diabetes Autoimmunity Study in the Young

OBJECTIVE

Mechanisms underlying the role of non-human leukocyte antigen (HLA) genetic risk variants in type 1 diabetes (T1D) are poorly understood. We aimed to test the association between methylation and non-HLA genetic risk.

METHODS

We conducted a methylation quantitative trait loci (mQTL) analysis in a nested case-control study from the Dietary Autoimmunity Study in the Young. Controls (n = 83) were frequency-matched to T1D cases (n = 83) based on age, race/ethnicity, and sample availability. We evaluated 13 non-HLA genetic markers known be associated with T1D. Genome-wide methylation profiling was performed on peripheral blood samples collected prior to T1D using the Illumina 450 K (discovery set) and infinium methylation EPIC beadchip (EPIC validation) platforms. Linear regression models, adjusting for age and sex, were used to test to each single nucleotide polymorphism (SNP) -probe combination. Logistic regression models were used to test the association between T1D and methylation levels among probes with a significant mQTL. A meta-analysis was used to combine odds ratios from the two platforms.

RESULTS

We identified 10 SNP-methylation probe pairs (false discovery rate (FDR) adjusted P < .05 and validation P < .05). Probes were associated with the GSDMB, C1QTNF6, IL27, and INS genes. The cg03366382 (OR: 1.9, meta-P = .0495), cg21574853 (OR: 2.5, meta-P = .0232), and cg25336198 (odds ratio: 6.6, meta-P = .0081) probes were significantly associated with T1D. The three probes were located upstream from the INS transcription start site.

CONCLUSIONS

We confirmed an association between DNA methylation and rs689 that has been identified in related studies. Measurements in our study preceded the onset of T1D suggesting methylation may have a role in the relationship between INS variation and T1D development.

SNP2APA: a database for evaluating effects of genetic variants on alternative polyadenylation in human cancers

Alternative polyadenylation (APA) is an important post-transcriptional regulation that recognizes different polyadenylation signals (PASs), resulting in transcripts with different 3' untranslated regions, thereby influencing a series of biological processes and functions. Recent studies have revealed that some single nucleotide polymorphisms (SNPs) could contribute to tumorigenesis and development through dysregulating APA. However, the associations between SNPs and APA in human cancers remain largely unknown. Here, using genotype and APA data of 9082 samples from The Cancer Genome Atlas (TCGA) and The Cancer 3'UTR Altas (TC3A), we systematically identified SNPs affecting APA events across 32 cancer types and defined them as APA quantitative trait loci (apaQTLs). As a result, a total of 467 942 cis-apaQTLs and 30 721 trans-apaQTLs were identified. By integrating apaQTLs with survival and genome-wide association studies (GWAS) data, we further identified 2154 apaQTLs associated with patient survival time and 151 342 apaQTLs located in GWAS loci. In addition, we designed an online tool to predict the effects of SNPs on PASs by utilizing PAS motif prediction tool. Finally, we developed SNP2APA, a user-friendly and intuitive database (http://gong_lab.hzau.edu.cn/SNP2APA/) for data browsing, searching, and downloading. SNP2APA will significantly improve our understanding of genetic variants and APA in human cancers.

Epigenetic Regulation of Alternative mRNA Splicing in Dilated Cardiomyopathy

In recent years, the genetic architecture of dilated cardiomyopathy (DCM) has been more thoroughly elucidated. However, there is still insufficient knowledge on the modifiers and regulatory principles that lead to the failure of myocardial function. The current study investigates the association of epigenome-wide DNA methylation and alternative splicing, both of which are important regulatory principles in DCM. We analyzed screening and replication cohorts of cases and controls and identified distinct transcriptomic patterns in the myocardium that differ significantly, and we identified a strong association of intronic DNA methylation and flanking exons usage (p < 2 × 10^-16). By combining differential exon usage (DEU) and differential methylation regions (DMR), we found a significant change of regulation in important sarcomeric and other DCM-associated pathways. Interestingly, inverse regulation of Titin antisense non-coding RNA transcript splicing and DNA methylation of a locus reciprocal to TTN substantiate these findings and indicate an additional role for non-protein-coding transcripts. In summary, this study highlights for the first time the close interrelationship between genetic imprinting by DNA methylation and the transport of this epigenetic information towards the dynamic mRNA splicing landscape. This expands our knowledge of the genome-environment interaction in DCM besides simple gene expression regulation.

The sperm epigenome does not display recurrent epimutations in patients with severely impaired spermatogenesis

BACKGROUND

In the past 15 years, numerous studies have described aberrant DNA methylation of imprinted genes (e.g. MEST and H19) in sperm of oligozoospermic men, but the prevalence and genomic extent of abnormal methylation patterns have remained unknown.

RESULTS

Using deep bisulfite sequencing (DBS), we screened swim-up sperm samples from 40 normozoospermic and 93 patients diagnosed as oligoasthenoteratozoospermic, oligoteratozoospermic or oligozoospermic, which are termed OATs throughout the manuscript, for H19 and MEST methylation. Based on this screening, we defined three patient groups: normal controls (NC), abnormally methylated oligozoospermic (AMO; n = 7) and normally methylated oligozoospermic (NMO; n = 86). Whole-genome bisulfite sequencing (WGBS) of five NC and five AMO samples revealed abnormal methylation levels of all 50 imprinting control regions in each AMO sample. To investigate whether this finding reflected epigenetic germline mosaicism or the presence of residual somatic DNA, we made a genome-wide inventory of soma-germ cell-specific DNA methylation. We found that > 2000 germ cell-specific genes are promoter-methylated in blood and that AMO samples had abnormal methylation levels at these genes, consistent with the presence of somatic cell DNA. The comparison between the five NC and six NMO samples revealed 19 differentially methylated regions (DMRs), none of which could be validated in an independent cohort of 40 men. Previous studies reported a higher incidence of epimutations at single CpG sites in the CTCF-binding region 6 of H19 in infertile patients. DBS analysis of this locus, however, revealed an association between DNA methylation levels and genotype (rs2071094), but not fertility phenotype.

CONCLUSIONS

Our results suggest that somatic DNA contamination and genetic variation confound methylation studies in sperm of infertile men. While we cannot exclude the existence of rare patients with slightly abnormal sperm methylation at non-recurrent CpG sites, the prevalence of aberrant methylation in swim-up purified sperm of infertile men has likely been overestimated, which is reassuring for patients undergoing assisted reproduction.

Missing heritability in Parkinson's disease: the emerging role of non-coding genetic variation

Parkinson's disease (PD) is a neurodegenerative disorder caused by a complex interplay of genetic and environmental factors. For the stratification of PD patients and the development of advanced clinical trials, including causative treatments, a better understanding of the underlying genetic architecture of PD is required. Despite substantial efforts, genome-wide association studies have not been able to explain most of the observed heritability. The majority of PD-associated genetic variants are located in non-coding regions of the genome. A systematic assessment of their functional role is hampered by our incomplete understanding of genotype-phenotype correlations, for example through differential regulation of gene expression. Here, the recent progress and remaining challenges for the elucidation of the role of non-coding genetic variants is reviewed with a focus on PD as a complex disease with multifactorial origins. The function of gene regulatory elements and the impact of non-coding variants on them, and the means to map these elements on a genome-wide level, will be delineated. Moreover, examples of how the integration of functional genomic annotations can serve to identify disease-associated pathways and to prioritize disease- and cell type-specific regulatory variants will be given. Finally, strategies for functional validation and considerations for suitable model systems are outlined. Together this emphasizes the contribution of rare and common genetic variants to the complex pathogenesis of PD and points to remaining challenges for the dissection of genetic complexity that may allow for better stratification, improved diagnostics and more targeted treatments for PD in the future.

DNA methylation in bull spermatozoa: evolutionary impacts, interindividual variability, and contribution to the embryo

The DNA methylome of spermatozoa results from a unique epigenetic reprogramming crucial for chromatin compaction and the protection of the paternal genetic heritage. Although bull semen is widely used for artificial insemination (AI), little is known about the sperm epigenome in cattle. The purpose of this review is to synthetize recent work on the bull sperm methylome in light of the knowledge accumulated in humans and model species. We will address sperm-specific DNA methylation features and their potential evolutionary impacts, with particular emphasis on hypomethylated regions and repetitive elements. We will review recent examples of interindividual variability and intra-individual plasticity of the bull sperm methylome as related to fertility and age, respectively. Finally, we will address paternal methylome reprogramming after fertilization, as well as the mechanisms potentially involved in epigenetic inheritance, and provide some examples of disturbances that alter the dynamics of reprogramming in cattle. Because the selection of AI bulls is closely based on their genotypes, we will also discuss the complex interplay between sequence polymorphism and DNA methylation, which represents both a difficulty in addressing the role of DNA methylation in shaping phenotypes and an opportunity to better understand genome plasticity.