Genetic variants shape rheumatoid arthritis-specific transcriptomic features in CD4+ T cells through differential DNA methylation, explaining a substantial proportion of heritability.

OBJECTIVE: CD4+ T cells have been suggested as the most disease-relevant cell type in rheumatoid arthritis (RA) in which RA-risk non-coding variants exhibit allele-specific effects on regulation of RA-driving genes. This study aimed to understand RA-specific signatures in CD4+ T cells using multi-omics data, interpreting inter-omics relationships in shaping the RA transcriptomic landscape. METHODS: We profiled genome-wide variants, gene expression and DNA methylation in CD4+ T cells from 82 patients with RA and 40 healthy controls using high-throughput technologies. We investigated differentially expressed genes (DEGs) and differential methylated regions (DMRs) in RA and localised quantitative trait loci (QTLs) for expression and methylation. We then integrated these based on individual-level correlations to inspect DEG-regulating sources and investigated the potential regulatory roles of RA-risk variants by a partitioned-heritability enrichment analysis with RA genome-wide association summary statistics. RESULTS: A large number of RA-specific DEGs were identified (n=2575), highlighting T cell differentiation and activation pathways. RA-specific DMRs, preferentially located in T cell regulatory regions, were correlated with the expression levels of 548 DEGs mostly in the same topologically associating domains. In addition, expressional variances in 771 and 83 DEGs were partially explained by expression QTLs for DEGs and methylation QTLs (meQTLs) for DEG-correlated DMRs, respectively. A large number of RA variants were moderately to strongly correlated with meQTLs. DEG-correlated DMRs, enriched with meQTLs, had strongly enriched heritability of RA. CONCLUSION: Our findings revealed that the methylomic changes, driven by RA heritability-explaining variants, shape the differential expression of a substantial fraction of DEGs in CD4+ T cells in patients with RA, reinforcing the importance of a multidimensional approach in disease-relevant tissues.

Epigenetic analysis in rheumatoid arthritis synoviocytes.

Rheumatoid arthritis (RA) is a complex chronic systematic disease with progressive destruction of the joints by invasive synoviocytes. To characterize the key regulators involved in the development of RA, we obtained multilayer epigenomics data including DNA methylation by whole-genome bisulfite sequencing, miRNA profiles, genetic variations by whole-exome sequencing, and mRNA profiles from synoviocytes of RA and osteoarthritis (OA) patients. The overall DNA methylation patterns were not much different between RA and OA, but 523 low-methylated regions (LMRs) were specific to RA. The LMRs were preferentially localized at the 5' introns and overlapped with transcription factor binding motifs for GLI1, RUNX2, and TFAP2A/C. Single base-scale differentially methylated CpGs were linked with several networks related to wound response, tissue development, collagen fibril organization, and the TGF-ß receptor signaling pathway. Further, the DNA methylation of 201 CpGs was significantly correlated with 27 expressed miRNA genes. Our interpretation of epigenomic data of the synoviocytes from RA and OA patients is an informative resource to further investigate regulatory elements and biomarkers responsible for the pathophysiology of RA and OA.

Genome-wide methylation analysis identifies ELOVL5 as an epigenetic biomarker for the risk of type 2 diabetes mellitus.

Genome-wide DNA methylation has been implicated in complex human diseases. Here, we identified epigenetic biomarkers for type 2 diabetes (T2D) underlying obesogenic environments. In a blood-based DNA methylation analysis of 11 monozygotic twins (MZTW) discordant for T2D, we discovered genetically independent candidate methylation sites. In a follow-up replication study (17 MZTW pairs) for external validation, we replicated the T2D-association at a novel CpG signal in the ELOVL fatty acid elongase 5 (ELOVL5) gene specific to T2D-discordant MZTW. For concordant DNA methylation signatures in tissues, we further confirmed that a CpG site (cg18681426) was associated with adipogenic differentiation between human preadipocytes and adipocytes isolated from the same biopsy sample. In addition, the ELOVL5 gene was significantly differentially expressed in adipose tissues from unrelated T2D patients and in human pancreatic islets. Our results demonstrate that blood-derived DNA methylation is associated with T2D risk as a proxy for cumulative epigenetic status in human adipose and pancreatic tissues. Moreover, ELOVL5 expression was increased in cellular and mouse models of induced obesity-related diabetes. These findings may provide new insights into epigenetic architecture by uncovering methylation-based biomarkers.

Data of methylome and transcriptome derived from human dilated cardiomyopathy.

Alterations in DNA methylation and gene expression have been implicated in the development of human dilated cardiomyopathy (DCM). Differentially methylated probes (DMPs) and differentially expressed genes (DEGs) were identified between the left ventricle (LV, a pathological locus for DCM) and the right ventricle (RV, a proxy for normal hearts). The data in this DiB are for supporting our report entitled "Methylome analysis reveals alterations in DNA methylation in the regulatory regions of left ventricle development genes in human dilated cardiomyopathy" (Bong-Seok Jo, In-Uk Koh, Jae-Bum Bae, Ho-Yeong Yu, Eun-Seok Jeon, Hae-Young Lee, Jae-Joong Kim, Murim Choi, Sun Shim Choi, 2016) [1].

Methylome analysis reveals alterations in DNA methylation in the regulatory regions of left ventricle development genes in human dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is one of the main causes of heart failure (called cardiomyopathies) in adults. Alterations in epigenetic regulation (i.e., DNA methylation) have been implicated in the development of DCM. Here, we identified a total of 1828 differentially methylated probes (DMPs) using the Infinium 450K HumanMethylation Bead chip by comparing the methylomes between 18 left ventricles and 9 right ventricles. Alterations in DNA methylation levels were observed mainly in lowly methylated regions corresponding to promoter-proximal regions, which become hypermethylated in severely affected left ventricles. Subsequent mRNA microarray analysis showed that the effect of DNA methylation on gene expression regulation is not unidirectional but is controlled by the functional sub-network context. DMPs were significantly enriched in the transcription factor binding sites (TFBSs) we tested. Alterations in DNA methylation were specifically enriched in the cis-regulatory regions of cardiac development genes, the majority of which are involved in ventricular development (e.g., TBX5 and HAND1).

Genome-wide association meta-analysis identifies novel variants associated with fasting plasma glucose in East Asians.

Fasting plasma glucose (FPG) has been recognized as an important indicator for the overall glycemic state preceding the onset of metabolic diseases. So far, most indentified genome-wide association loci for FPG were derived from populations with European ancestry, with a few exceptions. To extend a thorough catalog for FPG loci, we conducted meta-analyses of 13 genome-wide association studies in up to 24,740 nondiabetic subjects with East Asian ancestry. Follow-up replication analyses in up to an additional 21,345 participants identified three new FPG loci reaching genome-wide significance in or near PDK1-RAPGEF4, KANK1, and IGF1R. Our results could provide additional insight into the genetic variation implicated in fasting glucose regulation.

HBx induces hypomethylation of distal intragenic CpG islands required for active expression of developmental regulators.

Epigenetic alterations caused by viral oncoproteins are strong initiation factors for cancer development, but their mechanisms are largely unknown. To identify the epigenetic effects of viral hepatitis B virus X (HBx) that lead to hepatocellular carcinoma (HCC), we profiled the DNA methylomes of normal and HBx transgenic mouse liver. Intriguingly, severe hypomethylation of intragenic CpG islands (CGIs) was observed in HBx liver before the full development of HCC. Normally, these CGIs were highly methylated (mCGIs) by the DNMT3L complex and marked with epigenetic signatures associated with active expression, such as H3K36me3. Hypomethylation of mCGI was caused by the downregulation of Dnmt3L and Dnmt3a due to HBx bound to their promoters, along with HDAC1. These events lead to the downregulation of many developmental regulators that could facilitate tumorigenesis. Here we provide an intriguing epigenetic regulation mediated by mCGI that is required for cell differentiation and describe a previously unidentified epigenetic role for HBx in promoting HCC development.

Genome-wide methylation profiling of ADPKD identified epigenetically regulated genes associated with renal cyst development.

Autosomal dominant polycystic kidney disease (ADPKD) is a common human genetic disease characterized by the formation of multiple fluid-filled cysts in bilateral kidneys. Although mutations in polycystic kidney disease 1 (PKD1) are predominantly responsible for ADPKD, the focal and sporadic property of individual cystogenesis suggests another molecular mechanism such as epigenetic alterations. To determine the epigenomic alterations in ADPKD and their functional relevance, ADPKD and non-ADPKD individuals were analyzed by unbiased methylation profiling genome-wide and compared with their expression data. Intriguingly, PKD1 and other genes related to ion transport and cell adhesion were hypermethylated in gene-body regions, and their expressions were downregulated in ADPKD, implicating epigenetic silencing as the key mechanism underlying cystogenesis. Especially, in patients with ADPKD, PKD1 was hypermethylated in gene-body region and it was associated with recruitment of methyl-CpG-binding domain 2 proteins. Moreover, treatment with DNA methylation inhibitors retarded cyst formation of Madin-Darby Canine Kidney cells, accompanied with the upregulation of Pkd1 expression. These results are consistent with previous studies that knock-down of PKD1 was sufficient for cystogenesis. Therefore, our results reveal a critical role for hypermethylation of PKD1 and cystogenesis-related regulatory genes in cyst development, suggesting epigenetic therapy as a potential treatment for ADPKD.

Perspectives of international human epigenome consortium.

As the International Human Epigenome Consortium (IHEC) launched officially at the 2010 Washington meeting, a giant step toward the conquest of unexplored regions of the human genome has begun. IHEC aims at the production of 1,000 reference epigenomes to the international scientific community for next 7-10 years. Seven member institutions, including South Korea, Korea National Institute of Health (KNIH), will produce 25-200 reference epigenomes individually, and the produced data will be publically available by using a data center. Epigenome data will cover from whole genome bisulfite sequencing, histone modification, and chromatin access information to miRNA-seq. The final goal of IHEC is the production of reference maps of human epigenomes for key cellular status relevant to health and disease.

The characteristics of genome-wide DNA methylation in naïve CD4+ T cells of patients with psoriasis or atopic dermatitis.

Psoriasis and atopic dermatitis (AD) are skin diseases that are characterized by polarized CD4+ T cell responses. During the polarization of naïve CD4+ T cells, DNA methylation plays an important role in the regulation of gene transcription. In this study, we profiled the genome-wide DNA methylation status of naïve CD4+ T cells in patients with psoriasis or AD and healthy controls using a ChIP-seq method. Only psoriasis patient T cells, not those of AD patients, showed distinct hypomethylation (>4-fold) compared to healthy control T cells in twenty-six regions of the genome ranging in size from 10 to 70 kb. These regions were mostly pericentromeric on 10 different chromosomes and incidentally coincided with various strong epigenomic signals, such as histone modifications and transcription factor binding sites, that had been observed in the ENCODE project implying the potential epigenetic regulation in psoriasis development. The gene-centric analysis indicated that the promoter regions of 121 genes on the X chromosome had dramatically elevated methylation levels in psoriasis patient T-cells compared to those from healthy controls (>4-fold). Moreover, immune-related genes on the X chromosome had higher hypermethylation than other genes (P=0.046). No such patterns were observed with AD patient T cells. These findings imply that methylation changes in naïve CD4+ T cells may affect CD4+ T cell polarization, especially in the pathogenesis of psoriasis.

Genome-wide analysis of DNA methylation and the gene expression change in lung cancer.

INTRODUCTION: The recent DNA methylation studies on cancers have revealed the necessity of profiling an entire human genome and not to restrict the profiling to specific regions of the human genome. It has been suggested that genome-wide DNA methylation analysis enables us to identify the genes that are regulated by DNA methylation in carcinogenesis. METHODS: So, we performed whole-genome DNA methylation analysis for human lung squamous cell carcinoma (SCC), which is strongly related with smoking. We also performed microarrays using 21 pairs of normal lung tissues and tumors from patients with SCC. By combining these data, 30 hypermethylated and down-regulated genes, and 22 hypomethylated and up-regulated genes were selected. The gene expression level and DNA methylation pattern were confirmed by semiquantitative reverse-transcriptase polymerase chain reaction and pyrosequencing, respectively. RESULTS: By these validations, we selected five hypermethylated and down-regulated genes and one hypomethylated and up-regulated gene. Moreover, these six genes were proven to be actually regulated by DNA methylation by confirming the recovery of their DNA methylation pattern and gene expression level using a demethylating agent. The DNA methylation pattern of the CYTL1 promoter region was significantly different between early and advanced stages of SCC. CONCLUSION: In conclusion, by combining the whole-genome DNA methylation pattern and the gene expression profile, we identified the six genes (CCDC37, CYTL1, CDO1, SLIT2, LMO3, and SERPINB5) that are regulated by DNA methylation, and we suggest their value as target molecules for further study of SCC.

Identification of DNA methylation changes associated with human gastric cancer.

BACKGROUND: Epigenetic alteration of gene expression is a common event in human cancer. DNA methylation is a well-known epigenetic process, but verifying the exact nature of epigenetic changes associated with cancer remains difficult. METHODS: We profiled the methylome of human gastric cancer tissue at 50-bp resolution using a methylated DNA enrichment technique (methylated CpG island recovery assay) in combination with a genome analyzer and a new normalization algorithm. RESULTS: We were able to gain a comprehensive view of promoters with various CpG densities, including CpG Islands (CGIs), transcript bodies, and various repeat classes. We found that gastric cancer was associated with hypermethylation of 5' CGIs and the 5'-end of coding exons as well as hypomethylation of repeat elements, such as short interspersed nuclear elements and the composite element SVA. Hypermethylation of 5' CGIs was significantly correlated with downregulation of associated genes, such as those in the HOX and histone gene families. We also discovered long-range epigenetic silencing (LRES) regions in gastric cancer tissue and identified several hypermethylated genes (MDM2, DYRK2, and LYZ) within these regions. The methylation status of CGIs and gene annotation elements in metastatic lymph nodes was intermediate between normal and cancerous tissue, indicating that methylation of specific genes is gradually increased in cancerous tissue. CONCLUSIONS: Our findings will provide valuable data for future analysis of CpG methylation patterns, useful markers for the diagnosis of stomach cancer, as well as a new analysis method for clinical epigenomics investigations.

Nucleosome deposition and DNA methylation at coding region boundaries.

BACKGROUND: Nucleosome deposition downstream of transcription initiation and DNA methylation in the gene body suggest that control of transcription elongation is a key aspect of epigenetic regulation. RESULTS: Here we report a genome-wide observation of distinct peaks of nucleosomes and methylation at both ends of a protein coding unit. Elongating polymerases tend to pause near both coding ends immediately upstream of the epigenetic peaks, causing a significant reduction in elongation efficiency. Conserved features in underlying protein coding sequences seem to dictate their evolutionary conservation across multiple species. The nucleosomal and methylation marks are commonly associated with high sequence-encoded DNA-bending propensity but differentially with CpG density. As the gene grows longer, the epigenetic codes seem to be shifted from variable inner sequences toward boundary regions, rendering the peaks more prominent in higher organisms. CONCLUSIONS: Recent studies suggest that epigenetic inhibition of transcription elongation facilitates the inclusion of constitutive exons during RNA splicing. The epigenetic marks we identified here seem to secure the first and last coding exons from exon skipping as they are indispensable for accurate translation.

SNM1B/Apollo interacts with astrin and is required for the prophase cell cycle checkpoint.

Previously, we have shown that SNM1A is a multifunctional gene involved in both the DNA damage response and in an early mitotic checkpoint in response to spindle stress. Another member of the SNM1 gene family, SNM1B/Apollo, has been shown to have roles in both the response to DNA interstrand cross-linking agents and in telomere protection during S phase. Here, we demonstrate a novel role for SNM1B/Apollo in mitosis in response to spindle stress. SNM1B-deficient cells exhibit a defect in the prophase checkpoint. Loss of the prophase checkpoint induces an extended mitotic delay, which is due to prolonged activation of the spindle checkpoint. In addition, we show that SNM1B/Apollo interacts with the essential microtubule binding protein Astrin. SNM1B/Apollo interacts with Astrin through its conserved metallo-beta-lactamase domain, and disruption of this interaction by point mutations results in a deficient prophase checkpoint. These findings suggest that SNM1B/Apollo and Astrin function together to enforce the prophase checkpoint in response to spindle stress.

Redox-dependent changes in RsrA, an anti-sigma factor in Streptomyces coelicolor: zinc release and disulfide bond formation.

sigmaR is a sigma factor for transcribing genes to defend cells against oxidative stresses in the antibiotic-producing bacterium Streptomyces coelicolor. The availability of sigmaR is regulated by RsrA, an anti-sigma factor, whose sigmaR-binding activity is regulated by redox changes in the environment, via thiol-disulfide exchange. We found that reduced RsrA contains zinc in a stoichiometric amount, whereas oxidized form has very little: 1 mol of zinc per mol of RsrA was released upon oxidation as monitored by a chromogenic Zn-chelator, 4-(2-pyridylazo)-resorcinol (PAR). Measurement of zinc bound in several RsrA mutants of various cysteine and histidine substitutions suggested that C3, H7, C41, and C44 serve as zinc-binding sites. The zinc-binding and sigmaR-binding activities of mutant proteins did not coincide, suggesting that zinc might not be absolutely required for the anti-sigma activity of RsrA. Zn-free apo-RsrA bound sigmaR and inhibited sigmaR-dependent transcription in vitro. Compared with Zn-RsrA, the anti-transcription activity of apo-RsrA was about threefold lower and its sigmaR-binding affinity decreased by about ninefold when measured by surface plasmon resonance analysis. Apo-RsrA was more sensitive to protease, suggesting that zinc allows RsrA to maintain a more compact structure, optimized for binding sigmaR. The cysteine pairs that form disulfide bonds were determined by MALDI-TOF mass spectrometry, revealing formation of the critical disulfide bond between C11 and one of the essential cysteine residues C41 or 44, most likely C44. An improved model for the mechanism of redox-modulation of RsrA was presented.