Bridging epigenomics and complex disease: the basics

A CNN based m5c RNA methylation predictor

Post-transcriptional modifications of RNA play a key role in performing a variety of biological processes, such as stability and immune tolerance, RNA splicing, protein translation and RNA degradation. One of these RNA modifications is m5c which participates in various cellular functions like RNA structural stability and translation efficiency, got popularity among biologists. By applying biological experiments to detect RNA m5c methylation sites would require much more efforts, time and money. Most of the researchers are using pre-processed RNA sequences of 41 nucleotides where the methylated cytosine is in the center. Therefore, it is possible that some of the information around these motif may have lost. The conventional methods are unable to process the RNA sequence directly due to high dimensionality and thus need optimized techniques for better features extraction. To handle the above challenges the goal of this study is to employ an end-to-end, 1D CNN based model to classify and interpret m5c methylated data sites. Moreover, our aim is to analyze the sequence in its full length where the methylated cytosine may not be in the center. The evaluation of the proposed architecture showed a promising results by outperforming state-of-the-art techniques in terms of sensitivity and accuracy. Our model achieve 96.70% sensitivity and 96.21% accuracy for 41 nucleotides sequences while 96.10% accuracy for full length sequences.

Signaling pathways in cancer metabolism: mechanisms and therapeutic targets

A wide spectrum of metabolites (mainly, the three major nutrients and their derivatives) can be sensed by specific sensors, then trigger a series of signal transduction pathways and affect the expression levels of genes in epigenetics, which is called metabolite sensing. Life body regulates metabolism, immunity, and inflammation by metabolite sensing, coordinating the pathophysiology of the host to achieve balance with the external environment. Metabolic reprogramming in cancers cause different phenotypic characteristics of cancer cell from normal cell, including cell proliferation, migration, invasion, angiogenesis, etc. Metabolic disorders in cancer cells further create a microenvironment including many kinds of oncometabolites that are conducive to the growth of cancer, thus forming a vicious circle. At the same time, exogenous metabolites can also affect the biological behavior of tumors. Here, we discuss the metabolite sensing mechanisms of the three major nutrients and their derivatives, as well as their abnormalities in the development of various cancers, and discuss the potential therapeutic targets based on metabolite-sensing signaling pathways to prevent the progression of cancer.

DNA Methylation Research in Autologous Hematopoietic Stem Cell Transplant Population

Despite increased sophistication in DNA methylation (DNAm) measurement and methods, conducting studies in specific populations such as the hematopoietic stem cell transplant (HCT) population, presents unique challenges and study design considerations. In this article, we explain the motivation for investigating DNAm in the HCT population, highlighting important study design features and key findings in a longitudinal prospective pilot study of DNAm in 32 patients undergoing autologous HCT in Central Virginia, USA. We also discuss limitations and challenges to generating robust results. We observed that HCT does not prevent high-quality DNA from being extracted from whole blood for DNAm research and that longitudinal prospective studies that span pre- and 2-months post-HCT are feasible. Critically, we did not observe significant impacts of cancer diagnosis, time since transplant, age, or chromosomal sex on overall DNAm data dimensionality. These observations demonstrate that while extreme care is required to ensure generalizable, accurate, and interpretable results, researchers should not avoid HCT-DNAm research simply for fear that the transplant procedure or presence of a cancer diagnosis will prevent meaningful conclusions from being drawn. DNAm is an attractive biomarker that is understudied in patients undergoing HCT and needs to expand to improve precise prediction of HCT outcomes.

Epigenome-Wide Study of Posttraumatic Stress Disorder Symptom Severity in a Treatment-Seeking Adolescent Sample

Emerging research has demonstrated that psychosocial trauma exposure may elicit epigenetic changes, with downstream effects on the transcriptional regulation of genes. Epigenome-wide association studies (EWAS) offer an agnostic approach to examine DNA methylation (DNAm) associations and are a valuable tool to aid in the identification of biological pathways involved in posttraumatic stress disorder (PTSD). This study represents the first EWAS of PTSD in an adolescent sample, an important group given the significance of this developmental period regarding both DNAm changes and PTSD risk. The sample (n = 39, M age = 15.41 years, SD = 1.27, 84.6% female) comprised adolescents who experienced interpersonal trauma and were enrolled in a treatment study. Participants were assessed using the UCLA PTSD Reaction Index for DSM-IV-Adolescent Version and provided a blood sample at baseline. Genomic DNA was isolated from whole blood and assayed using the Illumina Infinium MethylationEPIC BeadChip. The primary analysis estimated the associations among individual CpG sites and PTSD symptom scores. Of the 793,575 screened probes tested, two were significant at a false discovery rate (FDR) < 10%. Hypomethylation of both sites was associated with increased PTSD symptom scores. Analysis of differentially methylated regions (DMR) identified a DMR associated with PTSD symptom scores at an FDR < 10%. Results from follow-up models are also discussed. Findings from this preliminary investigation suggest the importance of further research conducted in adolescent samples. The analytic pipeline and results are documented for use in future meta-analytic work as more such samples become available.

Stress resets ancestral heritable small RNA responses

Transgenerational inheritance of small RNAs challenges basic concepts of heredity. In Caenorhabditis elegans nematodes, small RNAs are transmitted across generations to establish a transgenerational memory trace of ancestral environments and distinguish self-genes from non-self-elements. Carryover of aberrant heritable small RNA responses was shown to be maladaptive and to lead to sterility. Here, we show that various types of stress (starvation, high temperatures, and high osmolarity) induce resetting of ancestral small RNA responses and a genome-wide reduction in heritable small RNA levels. We found that mutants that are defective in various stress pathways exhibit irregular RNAi inheritance dynamics even in the absence of stress. Moreover, we discovered that resetting of ancestral RNAi responses is specifically orchestrated by factors that function in the p38 MAPK pathway and the transcription factor SKN-1/Nrf2. Stress-dependent termination of small RNA inheritance could protect from run-on of environment-irrelevant heritable gene regulation.

Power analysis of transcriptome-wide association study: Implications for practical protocol choice

The transcriptome-wide association study (TWAS) has emerged as one of several promising techniques for integrating multi-scale 'omics' data into traditional genome-wide association studies (GWAS). Unlike GWAS, which associates phenotypic variance directly with genetic variants, TWAS uses a reference dataset to train a predictive model for gene expressions, which allows it to associate phenotype with variants through the mediating effect of expressions. Although effective, this core innovation of TWAS is poorly understood, since the predictive accuracy of the genotype-expression model is generally low and further bounded by expression heritability. This raises the question: to what degree does the accuracy of the expression model affect the power of TWAS? Furthermore, would replacing predictions with actual, experimentally determined expressions improve power? To answer these questions, we compared the power of GWAS, TWAS, and a hypothetical protocol utilizing real expression data. We derived non-centrality parameters (NCPs) for linear mixed models (LMMs) to enable closed-form calculations of statistical power that do not rely on specific protocol implementations. We examined two representative scenarios: causality (genotype contributes to phenotype through expression) and pleiotropy (genotype contributes directly to both phenotype and expression), and also tested the effects of various properties including expression heritability. Our analysis reveals two main outcomes: (1) Under pleiotropy, the use of predicted expressions in TWAS is superior to actual expressions. This explains why TWAS can function with weak expression models, and shows that TWAS remains relevant even when real expressions are available. (2) GWAS outperforms TWAS when expression heritability is below a threshold of 0.04 under causality, or 0.06 under pleiotropy. Analysis of existing publications suggests that TWAS has been misapplied in place of GWAS, in situations where expression heritability is low.

IMI - Myopia Genetics Report

The knowledge on the genetic background of refractive error and myopia has expanded dramatically in the past few years. This white paper aims to provide a concise summary of current genetic findings and defines the direction where development is needed. We performed an extensive literature search and conducted informal discussions with key stakeholders. Specific topics reviewed included common refractive error, any and high myopia, and myopia related to syndromes. To date, almost 200 genetic loci have been identified for refractive error and myopia, and risk variants mostly carry low risk but are highly prevalent in the general population. Several genes for secondary syndromic myopia overlap with those for common myopia. Polygenic risk scores show overrepresentation of high myopia in the higher deciles of risk. Annotated genes have a wide variety of functions, and all retinal layers appear to be sites of expression. The current genetic findings offer a world of new molecules involved in myopiagenesis. As the missing heritability is still large, further genetic advances are needed. This Committee recommends expanding large-scale, in-depth genetic studies using complementary big data analytics, consideration of gene-environment effects by thorough measurement of environmental exposures, and focus on subgroups with extreme phenotypes and high familial occurrence. Functional characterization of associated variants is simultaneously needed to bridge the knowledge gap between sequence variance and consequence for eye growth.

Mammalian SWI/SNF Enzymes and the Epigenetics of Tumor Cell Metabolic Reprogramming

Tumor cells reprogram their metabolism to survive and grow in a challenging microenvironment. Some of this reprogramming is performed by epigenetic mechanisms. Epigenetics is in turn affected by metabolism; chromatin modifying enzymes are dependent on substrates that are also key metabolic intermediates. We have shown that the chromatin remodeling enzyme Brahma-related gene 1 (BRG1), an epigenetic regulator, is necessary for rapid breast cancer cell proliferation. The mechanism for this requirement is the BRG1-dependent transcription of key lipogenic enzymes and regulators. Reduction in lipid synthesis lowers proliferation rates, which can be restored by palmitate supplementation. This work has established BRG1 as an attractive target for breast cancer therapy. Unlike genetic alterations, epigenetic mechanisms are reversible, promising gentler therapies without permanent off-target effects at distant sites.

Genomics era and complex disorders: Implications of GWAS with special reference to coronary artery disease, type 2 diabetes mellitus, and cancers

The Human Genome Project (HGP) has identified millions of single nucleotide polymorphisms (SNPs) and their association with several diseases, apart from successfully characterizing the Mendelian/monogenic diseases. However, the dissection of precise etiology of complex genetic disorders still poses a challenge for human geneticists. This review outlines the landmark results of genome-wide association studies (GWAS) with respect to major complex diseases - Coronary artery disease (CAD), type 2 diabetes mellitus (T2DM), and predominant cancers. A brief account on the current Indian scenario is also given. All the relevant publications till mid-2015 were accessed through web databases such as PubMed and Google. Several databases providing genetic information related to these diseases were tabulated and in particular, the list of the most significant SNPs identified through GWAS was made, which may be useful for designing studies in functional validation. Post-GWAS implications and emerging concepts such as epigenomics and pharmacogenomics were also discussed.

The significance of the increased expression of phosphorylated MeCP2 in the membranes from patients with proliferative diabetic retinopathy

The purpose of this study was to evaluate the correlation of expression of phosphorylated methyl-CpG binding protein 2-Ser421 (MeCP2-S421) and VEGF in the membranes of patients with PDR. We examined the expression of phospho-MeCP2-S80, S421, VEGF and PEDF in surgically excised PDR membranes from 33 patients with diabetes, and idiopathic epiretinal membranes from 11 patients without diabetes, using immunohistochemistry and western blot. The colocalization of MeCP2-S421 with VEGF, PEDF, CD31, GFAP and αSMA was revealed by fluorescent double labeling. The effect of CoCl2 and knock down MeCP2 using specific siRNA on the expression of MeCP2 and VEGF were analyzed in HUCAC cells by Western blot. We found that phospho-MeCP2-S421 was significantly increased in the membranes from the patients with PDR compared with the specimens from patients without diabetes (P < 0.01). The expression of phospho-MeCP2-S421 was much stronger than that of phospho-MeCP2-S80 in the PDR membranes. Double labeling showed that the high phospho-MeCP2-S421 expression was associated with strong expression of VEGF, but not PEDF. Further, phospho-MeCP2-S421 and VEGF were increased by the stimulation of CoCl2 and knock down MeCP2 inhibited the expression of VEGF. Our result suggests that phospho-MeCP2-S421 might involve in the pathogenesis of PDR.

Integrating Epigenomic Elements and GWASs Identifies BDNF Gene Affecting Bone Mineral Density and Osteoporotic Fracture Risk

To identify susceptibility genes for osteoporosis, we conducted an integrative analysis that combined epigenomic elements and previous genome-wide association studies (GWASs) data, followed by validation at population and functional levels, which could identify common regulatory elements and predict new susceptibility genes that are biologically meaningful to osteoporosis. By this approach, we found a set of distinct epigenomic elements significantly enriched or depleted in the promoters of osteoporosis-associated genes, including 4 transcription factor binding sites, 27 histone marks, and 21 chromatin states segmentation types. Using these epigenomic marks, we performed reverse prediction analysis to prioritize the discovery of new candidate genes. Functional enrichment analysis of all the prioritized genes revealed several key osteoporosis related pathways, including Wnt signaling. Genes with high priority were further subjected to validation using available GWASs datasets. Three genes were significantly associated with spine bone mineral density, including BDNF, PDE4D, and SATB2, which all closely related to bone metabolism. The most significant gene BDNF was also associated with osteoporotic fractures. RNA interference revealed that BDNF knockdown can suppress osteoblast differentiation. Our results demonstrated that epigenomic data could be used to indicate common epigenomic marks to discover additional loci with biological functions for osteoporosis.

Epigenetic regulation in heart failure: part II DNA and chromatin

Epigenetic regulatory mechanisms play key roles in cardiac development, differentiation, homeostasis, response to stress and injury, and disease. Human heart failure (HF) epigenetic regulatory mechanisms have not been deciphered to date. This 2-part review distills the rapidly evolving research focused on human HF epigenetic regulatory mechanisms. Part I, which was published in the September/October issue, focused on epigenetic regulatory mechanisms involving RNA, specifically the role of short, intermediate, and long noncoding RNAs (lncRNAs) and endogenous competing RNA regulatory networks. Part II, now in the November/December issue, focuses on the epigenetic regulatory mechanisms involving DNA, including DNA methylation, histone modifications, and chromatin conformational changes. Part II concludes with 2 examples of well-studied integrated epigenetic regulatory mechanisms: the structural and functional roles of the Mediator complex in regulating transcription and the epigenetic networked "cross-talk" regulating atrial natriuretic peptide and brain natriuretic peptide promoter activation.

Characterization of genome-methylome interactions in 22 nuclear pedigrees

Genetic polymorphisms can shape the global landscape of DNA methylation, by either changing substrates for DNA methyltransferases or altering the DNA binding affinity of cis-regulatory proteins. The interactions between CpG methylation and genetic polymorphisms have been previously investigated by methylation quantitative trait loci (mQTL) and allele-specific methylation (ASM) analysis. However, it remains unclear whether these approaches can effectively and comprehensively identify all genetic variants that contribute to the inter-individual variation of DNA methylation levels. Here we used three independent approaches to systematically investigate the influence of genetic polymorphisms on variability in DNA methylation by characterizing the methylation state of 96 whole blood samples in 52 parent-child trios from 22 nuclear pedigrees. We performed targeted bisulfite sequencing with padlock probes to quantify the absolute DNA methylation levels at a set of 411,800 CpG sites in the human genome. With mid-parent offspring analysis (MPO), we identified 10,593 CpG sites that exhibited heritable methylation patterns, among which 70.1% were SNPs directly present in methylated CpG dinucleotides. We determined the mQTL analysis identified 49.9% of heritable CpG sites for which regulation occurred in a distal cis-regulatory manner, and that ASM analysis was only able to identify 5%. Finally, we identified hundreds of clusters in the human genome for which the degree of variation of CpG methylation, as opposed to whether or not CpG sites were methylated, was associated with genetic polymorphisms, supporting a recent hypothesis on the genetic influence of phenotypic plasticity. These results show that cis-regulatory SNPs identified by mQTL do not comprise the full extent of heritable CpG methylation, and that ASM appears overall unreliable. Overall, the extent of genome-methylome interactions is well beyond what is detectible with the commonly used mQTL and ASM approaches, and is likely to include effects on plasticity.

Epigenetic modification of DRG neuronal gene expression subsequent to nerve injury: etiological contribution to complex regional pain syndromes (Part I)

DRG is of importance in relaying painful stimulation to the higher pain centers and therefore could be a crucial target for early intervention aimed at suppressing primary afferent stimulation. Complex regional pain syndrome (CRPS) is a common pain condition with an unknown etiology. Recently added new information enriches our understanding of CRPS pathophysiology. Researches on genetics, biogenic amines, neurotransmitters, and mechanisms of pain modulation, central sensitization, and autonomic functions in CRPS revealed various abnormalities indicating that multiple factors and mechanisms are involved in the pathogenesis of CRPS. Epigenetics refers to mitotically and meiotically heritable changes in gene expression that do not affect the DNA sequence. As epigenetic modifications potentially play an important role in inflammatory cytokine metabolism, neurotransmitter responsiveness, and analgesic sensitivity, they are likely key factors in the development of chronic pain. In this dyad review series, we systematically examine the nerve injury-related changes in the neurological system and their contribution to CRPS. In this part, we first reviewed and summarized the role of neural sensitization in DRG neurons in performing function in the context of pain processing. Particular emphasis is placed on the cellular and molecular changes after nerve injury as well as different models of inflammatory and neuropathic pain. These were considered as the potential molecular bases that underlie nerve injury-associated pathogenesis of CRPS.

DNA methylation is globally disrupted and associated with expression changes in chronic obstructive pulmonary disease small airways

DNA methylation is an epigenetic modification that is highly disrupted in response to cigarette smoke and involved in a wide spectrum of malignant and nonmalignant diseases, but surprisingly not previously assessed in small airways of patients with chronic obstructive pulmonary disease (COPD). Small airways are the primary sites of airflow obstruction in COPD. We sought to determine whether DNA methylation patterns are disrupted in small airway epithelia of patients with COPD, and evaluate whether changes in gene expression are associated with these disruptions. Genome-wide methylation and gene expression analysis were performed on small airway epithelial DNA and RNA obtained from the same patient during bronchoscopy, using Illumina's Infinium HM27 and Affymetrix's Genechip Human Gene 1.0 ST arrays. To control for known effects of cigarette smoking on DNA methylation, methylation and gene expression profiles were compared between former smokers with and without COPD matched for age, pack-years, and years of smoking cessation. Our results indicate that aberrant DNA methylation is (1) a genome-wide phenomenon in small airways of patients with COPD, and (2) associated with altered expression of genes and pathways important to COPD, such as the NF-E2-related factor 2 oxidative response pathway. DNA methylation is likely an important mechanism contributing to modulation of genes important to COPD pathology. Because these methylation events may underlie disease-specific gene expression changes, their characterization is a critical first step toward the development of epigenetic markers and an opportunity for developing novel epigenetic therapeutic interventions for COPD.

Epigenomic elements enriched in the promoters of autoimmunity susceptibility genes

Genome-wide association studies have identified a number of autoimmune disease-susceptibility genes. Whether or not these loci share any regulatory or functional elements, however, is an open question. Finding such common regulators is of considerable research interest in order to define systemic therapeutic targets. The growing amount of experimental genomic annotations, particularly those from the ENCODE project, provide a wealth of opportunities to search for such commonalities. We hypothesized that regulatory commonalities might not only delineate a regulatory landscape predisposing to autoimmune diseases, but also define functional elements distinguishing specific diseases. We further investigated if, and how, disease-specific epigenomic elements can identify novel genes yet to be associated with the diseases. We evaluated transcription factors, histone modifications, and chromatin state data obtained from the ENCODE project for statistically significant over- or under-representation in the promoters of genes associated with Systemic Lupus Erythematosus (SLE), Rheumatoid Arthritis (RA), and Systemic Sclerosis (SSc). We identified BATF, BCL11A, IRF4, NFkB, PAX5, and PU.1 as transcription factors over-represented in SLE- and RA-susceptibility gene promoters. H3K4me1 and H3K4me2 epigenomic marks were associated with SLE susceptibility genes, and H3K9me3 was common to both SLE and RA. In contrast to a transcriptionally active signature in SLE and RA, SSc-susceptibility genes were depleted in activating epigenomic elements. Using epigenomic elements enriched in SLE and RA, we identified additional immune and B cell signaling-related genes with the same elements in their promoters. Our analysis suggests common and disease-specific epigenomic elements that may define novel therapeutic targets for controlling aberrant activation of autoimmune susceptibility genes.

Autoantibody repertoires, natural biomarkers, and system controllers

The immune system is composed of networks of interacting cells and molecules; therefore, to understand and control immune behavior we need to adopt the thinking and tools of systems immunology. This review describes the use of an antigen microarray device and informatics to profile the repertoires of autoantibodies in health and disease. Autoantibody profiling provides an insight into the biomarkers used by the immune system in its dialog with the body. Heat shock protein 60 (HSP60) and HSP70 are cited as examples of key hubs in physiological regulatory networks; HSP molecules and peptides can be viewed as natural regulators because the immune system itself deploys them to modulate inflammatory reactions. The discovery of such natural biomarkers paves the way towards natural control.