Skin Cancer Risk Is Increased by Somatic Mutations Detected Noninvasively in Healthy-Appearing Sun-Exposed Skin.

Skin cancer risk is increased by exposure to ultraviolet radiation (UVR). Because UVR exposure accumulates over time and lighter skin is more susceptible to UVR, age and skin tone are risk factors for skin cancer. However, measurements of somatic mutations in healthy-appearing skin have not been used to calculate skin cancer risk. In this study, we developed a noninvasive test that quantifies somatic mutations in healthy-appearing sun-exposed skin and applied it to a 1038-subject cohort. Somatic mutations were combined with other known skin cancer risk factors to train a model to calculate risk. The final model (DNA-Skin Cancer Assessment of Risk) was trained to predict personal history of skin cancer from age, family history, skin tone, and mutation count. The addition of mutation count significantly improved model performance (OR = 1.3, 95% confidence interval = 1.14-1.48; P = 5.3 × 10-6) and made a more significant contribution than skin tone. Calculations of skin cancer risk matched the known United States population prevalence, indicating that DNA-Skin Cancer Assessment of Risk was well-calibrated. In conclusion, somatic mutations in healthy-appearing sun-exposed skin increase skin cancer risk, and mutations capture risk information that is not accounted for by other risk factors. Clinical utility is supported by the noninvasive nature of skin sample collection through adhesive patches.

Identification of DNA motifs that regulate DNA methylation.

DNA methylation is an important epigenetic mark but how its locus-specificity is decided in relation to DNA sequence is not fully understood. Here, we have analyzed 34 diverse whole-genome bisulfite sequencing datasets in human and identified 313 motifs, including 92 and 221 associated with methylation (methylation motifs, MMs) and unmethylation (unmethylation motifs, UMs), respectively. The functionality of these motifs is supported by multiple lines of evidence. First, the methylation levels at the MM and UM motifs are respectively higher and lower than the genomic background. Second, these motifs are enriched at the binding sites of methylation modifying enzymes including DNMT3A and TET1, indicating their possible roles of recruiting these enzymes. Third, these motifs significantly overlap with "somatic QTLs" (quantitative trait loci) of methylation and expression. Fourth, disruption of these motifs by mutation is associated with significantly altered methylation level of the CpGs in the neighbor regions. Furthermore, these motifs together with somatic mutations are predictive of cancer subtypes and patient survival. We revealed some of these motifs were also associated with histone modifications, suggesting a possible interplay between the two types of epigenetic modifications. We also found some motifs form feed forward loops to contribute to DNA methylation dynamics.

Epigenomic analysis reveals DNA motifs regulating histone modifications in human and mouse.

Histones are modified by enzymes that act in a locus, cell-type, and developmental stage-specific manner. The recruitment of enzymes to chromatin is regulated at multiple levels, including interaction with sequence-specific DNA-binding factors. However, the DNA-binding specificity of the regulatory factors that orchestrate specific histone modifications has not been broadly mapped. We have analyzed 6 histone marks (H3K4me1, H3K4me3, H3K27ac, H3K27me3, K3H9me3, H3K36me3) across 121 human cell types and tissues from the NIH Roadmap Epigenomics Project as well as 8 histone marks (with addition of H3K4me2 and H3K9ac) from the mouse ENCODE Consortium. We have identified 361 and 369 DNA motifs in human and mouse, respectively, that are the most predictive of each histone mark. Interestingly, 107 human motifs are conserved between the two species. In human embryonic cell line H1, we mutated only the found DNA motifs at particular loci and the significant reduction of H3K27ac levels validated the regulatory roles of the perturbed motifs. The functionality of these motifs was also supported by the evidence that histone-associated motifs, especially H3K4me3 motifs, significantly overlap with the expression of quantitative trait loci SNPs in cancer patients more than the known and random motifs. Furthermore, we observed possible feedbacks to control chromatin dynamics as the found motifs appear in the promoters or enhancers associated with various histone modification enzymes. These results pave the way toward revealing the molecular mechanisms of epigenetic events, such as histone modification dynamics and epigenetic priming.

Comprehensive epigenetic landscape of rheumatoid arthritis fibroblast-like synoviocytes.

Epigenetics contributes to the pathogenesis of immune-mediated diseases like rheumatoid arthritis (RA). Here we show the first comprehensive epigenomic characterization of RA fibroblast-like synoviocytes (FLS), including histone modifications (H3K27ac, H3K4me1, H3K4me3, H3K36me3, H3K27me3, and H3K9me3), open chromatin, RNA expression and whole-genome DNA methylation. To address complex multidimensional relationship and reveal epigenetic regulation of RA, we perform integrative analyses using a novel unbiased method to identify genomic regions with similar profiles. Epigenomically similar regions exist in RA cells and are associated with active enhancers and promoters and specific transcription factor binding motifs. Differentially marked genes are enriched for immunological and unexpected pathways, with "Huntington's Disease Signaling" identified as particularly prominent. We validate the relevance of this pathway to RA by showing that Huntingtin-interacting protein-1 regulates FLS invasion into matrix. This work establishes a high-resolution epigenomic landscape of RA and demonstrates the potential for integrative analyses to identify unanticipated therapeutic targets.

Epipolymorphisms associated with the clinical outcome of autoimmune arthritis affect CD4+ T cell activation pathways.

Multifactorial diseases, including autoimmune juvenile idiopathic arthritis (JIA), result from a complex interplay between genetics and environment. Epigenetic mechanisms are believed to integrate such gene-environment interactions, fine-tuning gene expression, and possibly contributing to immune system dysregulation. Although anti-TNF therapy has strongly increased JIA remission rates, it is not curative and up to 80% of patients flare upon treatment withdrawal. Thus, a crucial unmet medical and scientific need is to understand the immunological mechanisms associated with remission or flare to inform clinical decisions. Here, we explored the CD4+ T-cell DNA methylome of 68 poly-articular and extended oligo-articular JIA patients, before and after anti-TNF therapy withdrawal, to identify features associated with maintenance of inactive disease. Individual CpG sites were clustered in coherent modules without a priori knowledge of their function through network analysis. The methylation level of several CpG modules, specifically those enriched in CpG sites belonging to genes that mediate T-cell activation, uniquely correlated with clinical activity. Differences in DNA methylation were already detectable at the time of therapy discontinuation, suggesting epigenetic predisposition. RNA profiling also detected differences in T-cell activation markers (including HLA-DR) but, overall, its sensitivity was lower than epigenetic profiling. Changes to the T-cell activation signature at the protein level were detectable by flow cytometry, confirming the biological relevance of the observed alterations in methylation. Our work proposes epigenetic discrimination between clinical activity states, and reveals T-cell-related biological functions tied to, and possibly predicting or causing, clinical outcome.

The Transcriptional Response to DNA-Double-Strand Breaks in Physcomitrella patens.

The model bryophyte Physcomitrella patens is unique among plants in supporting the generation of mutant alleles by facile homologous recombination-mediated gene targeting (GT). Reasoning that targeted transgene integration occurs through the capture of transforming DNA by the homology-dependent pathway for DNA double-strand break (DNA-DSB) repair, we analysed the genome-wide transcriptomic response to bleomycin-induced DNA damage and generated mutants in candidate DNA repair genes. Massively parallel (Illumina) cDNA sequencing identified potential participants in gene targeting. Transcripts encoding DNA repair proteins active in multiple repair pathways were significantly up-regulated. These included Rad51, CtIP, DNA ligase 1, Replication protein A and ATR in homology-dependent repair, Xrcc4, DNA ligase 4, Ku70 and Ku80 in non-homologous end-joining and Rad1, Tebichi/polymerase theta, PARP in microhomology-mediated end-joining. Differentially regulated cell-cycle components included up-regulated Rad9 and Hus1 DNA-damage-related checkpoint proteins and down-regulated D-type cyclins and B-type CDKs, commensurate with the imposition of a checkpoint at G2 of the cell cycle characteristic of homology-dependent DNA-DSB repair. Candidate genes, including ATP-dependent chromatin remodelling helicases associated with repair and recombination, were knocked out and analysed for growth defects, hypersensitivity to DNA damage and reduced GT efficiency. Targeted knockout of PpCtIP, a cell-cycle activated mediator of homology-dependent DSB resection, resulted in bleomycin-hypersensitivity and greatly reduced GT efficiency.

Abnormal PTPN11 enhancer methylation promotes rheumatoid arthritis fibroblast-like synoviocyte aggressiveness and joint inflammation.

The PTPN11 gene, encoding the tyrosine phosphatase SHP-2, is overexpressed in rheumatoid arthritis (RA) fibroblast-like synoviocytes (FLS) compared with osteoarthritis (OA) FLS and promotes RA FLS invasiveness. Here, we explored the molecular basis for PTPN11 overexpression in RA FLS and the role of SHP-2 in RA pathogenesis. Using computational methods, we identified a putative enhancer in PTPN11 intron 1, which contained a glucocorticoid receptor- binding (GR-binding) motif. This region displayed enhancer function in RA FLS and contained 2 hypermethylation sites in RA compared with OA FLS. RA FLS stimulation with the glucocorticoid dexamethasone induced GR binding to the enhancer and PTPN11 expression. Glucocorticoid responsiveness of PTPN11 was significantly higher in RA FLS than OA FLS and required the differentially methylated CpGs for full enhancer function. SHP-2 expression was enriched in the RA synovial lining, and heterozygous Ptpn11 deletion in radioresistant or innate immune cells attenuated K/BxN serum transfer arthritis in mice. Treatment with SHP-2 inhibitor 11a-1 reduced RA FLS migration and responsiveness to TNF and IL-1ß stimulation and reduced arthritis severity in mice. Our findings demonstrate how abnormal epigenetic regulation of a pathogenic gene determines FLS behavior and demonstrate that targeting SHP-2 or the SHP-2 pathway could be a therapeutic strategy for RA.

Integrative genomic deconvolution of rheumatoid arthritis GWAS loci into gene and cell type associations.

BACKGROUND: Although genome-wide association studies (GWAS) have identified over 100 genetic loci associated with rheumatoid arthritis (RA), our ability to translate these results into disease understanding and novel therapeutics is limited. Most RA GWAS loci reside outside of protein-coding regions and likely affect distal transcriptional enhancers. Furthermore, GWAS do not identify the cell types where the associated causal gene functions. Thus, mapping the transcriptional regulatory roles of GWAS hits and the relevant cell types will lead to better understanding of RA pathogenesis. RESULTS: We combine the whole-genome sequences and blood transcription profiles of 377 RA patients and identify over 6000 unique genes with expression quantitative trait loci (eQTLs). We demonstrate the quality of the identified eQTLs through comparison to non-RA individuals. We integrate the eQTLs with immune cell epigenome maps, RA GWAS risk loci, and adjustment for linkage disequilibrium to propose target genes of immune cell enhancers that overlap RA risk loci. We examine 20 immune cell epigenomes and perform a focused analysis on primary monocytes, B cells, and T cells. CONCLUSIONS: We highlight cell-specific gene associations with relevance to RA pathogenesis including the identification of FCGR2B in B cells as possessing both intragenic and enhancer regulatory GWAS hits. We show that our RA patient cohort derived eQTL network is more informative for studying RA than that from a healthy cohort. While not experimentally validated here, the reported eQTLs and cell type-specific RA risk associations can prioritize future experiments with the goal of elucidating the regulatory mechanisms behind genetic risk associations.

LBH Gene Transcription Regulation by the Interplay of an Enhancer Risk Allele and DNA Methylation in Rheumatoid Arthritis.

OBJECTIVE: To identify nonobvious therapeutic targets for rheumatoid arthritis (RA), we performed an integrative analysis incorporating multiple "omics" data and the Encyclopedia of DNA Elements (ENCODE) database for potential regulatory regions. This analysis identified the limb bud and heart development (LBH) gene, which has risk alleles associated with RA/celiac disease and lupus, and can regulate cell proliferation in RA. We identified a novel LBH transcription enhancer with an RA risk allele (rs906868 G [Ref]/T) 6 kb upstream of the LBH gene with a differentially methylated locus. The confluence of 3 regulatory elements, rs906868, an RA differentially methylated locus, and a putative enhancer, led us to investigate their effects on LBH regulation in fibroblast-like synoviocytes (FLS). METHODS: We cloned the 1.4-kb putative enhancer with either the rs906868 Ref allele or single-nucleotide polymorphism (SNP) variant into reporter constructs. The constructs were methylated in vitro and transfected into cultured FLS by nucleofection. RESULTS: We found that both variants increased transcription, thereby confirming the region's enhancer function. Unexpectedly, the transcriptional activity of the Ref risk allele was significantly lower than that of the SNP variant and is consistent with low LBH levels as a risk factor for aggressive FLS behavior. Using RA FLS lines with a homozygous Ref or SNP allele, we confirmed that homozygous Ref lines expressed lower LBH messenger RNA levels than did the SNP lines. Methylation significantly reduced enhancer activity for both alleles, indicating that enhancer function is dependent on its methylation status. CONCLUSION: This study shows how the interplay between genetics and epigenetics can affect expression of LBH in RA.

Constructing 3D interaction maps from 1D epigenomes.

The human genome is tightly packaged into chromatin whose functional output depends on both one-dimensional (1D) local chromatin states and three-dimensional (3D) genome organization. Currently, chromatin modifications and 3D genome organization are measured by distinct assays. An emerging question is whether it is possible to deduce 3D interactions by integrative analysis of 1D epigenomic data and associate 3D contacts to functionality of the interacting loci. Here we present EpiTensor, an algorithm to identify 3D spatial associations within topologically associating domains (TADs) from 1D maps of histone modifications, chromatin accessibility and RNA-seq. We demonstrate that active promoter-promoter, promoter-enhancer and enhancer-enhancer associations identified by EpiTensor are highly concordant with those detected by Hi-C, ChIA-PET and eQTL analyses at 200 bp resolution. Moreover, EpiTensor has identified a set of interaction hotspots, characterized by higher chromatin and transcriptional activity as well as enriched TF and ncRNA binding across diverse cell types, which may be critical for stabilizing the local 3D interactions.

Human body epigenome maps reveal noncanonical DNA methylation variation.

Understanding the diversity of human tissues is fundamental to disease and requires linking genetic information, which is identical in most of an individual's cells, with epigenetic mechanisms that could have tissue-specific roles. Surveys of DNA methylation in human tissues have established a complex landscape including both tissue-specific and invariant methylation patterns. Here we report high coverage methylomes that catalogue cytosine methylation in all contexts for the major human organ systems, integrated with matched transcriptomes and genomic sequence. By combining these diverse data types with each individuals' phased genome, we identified widespread tissue-specific differential CG methylation (mCG), partially methylated domains, allele-specific methylation and transcription, and the unexpected presence of non-CG methylation (mCH) in almost all human tissues. mCH correlated with tissue-specific functions, and using this mark, we made novel predictions of genes that escape X-chromosome inactivation in specific tissues. Overall, DNA methylation in several genomic contexts varies substantially among human tissues.

Integrative omics analysis of rheumatoid arthritis identifies non-obvious therapeutic targets.

Identifying novel therapeutic targets for the treatment of disease is challenging. To this end, we developed a genome-wide approach of candidate gene prioritization. We independently collocated sets of genes that were implicated in rheumatoid arthritis (RA) pathogenicity through three genome-wide assays: (i) genome-wide association studies (GWAS), (ii) differentially expression in RA fibroblast-like synoviocytes (FLS), and (iii) differentially methylation in RA FLS. Integrated analysis of these complementary data sets identified a significant enrichment of multi-evidence genes (MEGs) within pathways relating to RA pathogenicity. One MEG is Engulfment and Cell Motility Protein-1 (ELMO1), a gene not previously considered as a therapeutic target in RA FLS. We demonstrated in RA FLS that ELMO1 is: (i) expressed, (ii) promotes cell migration and invasion, and (iii) regulates Rac1 activity. Thus, we created links between ELMO1 and RA pathogenicity, which in turn validates ELMO1 as a potential RA therapeutic target. This study illustrated the power of MEG-based approaches for therapeutic target identification.

Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

The Rheumatoid Arthritis Risk Gene LBH Regulates Growth in Fibroblast-like Synoviocytes.

OBJECTIVE: Fibroblast-like synoviocytes (FLS) are key players in the synovial pathology of rheumatoid arthritis (RA). Currently, there is no treatment that specifically targets these aggressive cells. By combining 3 different "omics" data sets, i.e., 1) risk genes in RA, 2) differentially expressed genes, and 3) differential DNA methylation in RA versus osteoarthritis (OA) FLS, we identified LBH (limb bud and heart development) as a candidate gene in RA. The present study was undertaken to define the role of this gene in FLS. METHODS: Synovial tissue specimens from RA and OA patients were collected at the time of joint replacement surgery. LBH expression was silenced using small interfering RNA or overexpressed using an LBH expression vector in primary FLS. Gene expression profiles were determined by microarray and assessed using Ingenuity Pathway Analysis. Effects of modified LBH expression were investigated in functional assays. RESULTS: LBH was expressed in the synovial lining layer in patients with RA. Transforming growth factor ß1 significantly increased LBH expression in primary FLS, and platelet-derived growth factor BB decreased it. Pathway analysis of the transcriptome of LBH-deficient FLS compared to control FLS identified "cellular growth and proliferation" as the most significantly enriched pathway. In growth assays, LBH deficiency increased FLS proliferation. Conversely, LBH overexpression significantly inhibited cell growth. Cell cycle analysis demonstrated a marked increase in cells entering the cell cycle in LBH-deficient FLS compared to controls. LBH did not alter apoptosis. CONCLUSION: LBH is a candidate gene for synovial pathology in RA. It is regulated by growth factors and modulates cell growth in primary FLS. Our data suggest a novel mechanism for synovial intimal hyperplasia and joint damage in RA.

Computational schemes for the prediction and annotation of enhancers from epigenomic assays.

Identifying and annotating distal regulatory enhancers is critical to understand the mechanisms that control gene expression and cell-type-specific activities. Next-generation sequencing techniques have provided us an exciting toolkit of genome-wide assays that can be used to predict and annotate enhancers. However, each assay comes with its own specific set of analytical needs if enhancer prediction is to be optimal. Furthermore, integration of multiple genome-wide assays allows for different genomic features to be combined, and can improve predictive performance. Herein, we review the genome-wide assays and analysis schemes that are used to predict and annotate enhancers. In particular, we focus on three key computational topics: predicting enhancer locations, determining the cell-type-specific activity of enhancers, and linking enhancers to their target genes.

Predicting the human epigenome from DNA motifs.

The epigenome is established and maintained by the site-specific recruitment of chromatin-modifying enzymes and their cofactors. Identifying the cis elements that regulate epigenomic modification is critical for understanding the regulatory mechanisms that control gene expression patterns. We present Epigram, an analysis pipeline that predicts histone modification and DNA methylation patterns from DNA motifs. The identified cis elements represent interactions with the site-specific DNA-binding factors that establish and maintain epigenomic modifications. We cataloged the cis elements in embryonic stem cells and four derived lineages and found numerous motifs that have location preference, such as at the center of H3K27ac or at the edges of H3K4me3 and H3K9me3, which provides mechanistic insight about the shaping of the epigenome. The Epigram pipeline and predictive motifs are at http://wanglab.ucsd.edu/star/epigram/.