LuxHMM: DNA methylation analysis with genome segmentation via hidden Markov model.

BACKGROUND: DNA methylation plays an important role in studying the epigenetics of various biological processes including many diseases. Although differential methylation of individual cytosines can be informative, given that methylation of neighboring CpGs are typically correlated, analysis of differentially methylated regions is often of more interest. RESULTS: We have developed a probabilistic method and software, LuxHMM, that uses hidden Markov model (HMM) to segment the genome into regions and a Bayesian regression model, which allows handling of multiple covariates, to infer differential methylation of regions. Moreover, our model includes experimental parameters that describe the underlying biochemistry in bisulfite sequencing and model inference is done using either variational inference for efficient genome-scale analysis or Hamiltonian Monte Carlo (HMC). CONCLUSIONS: Analyses of real and simulated bisulfite sequencing data demonstrate the competitive performance of LuxHMM compared with other published differential methylation analysis methods.

LuxRep: a technical replicate-aware method for bisulfite sequencing data analysis.

BACKGROUND: DNA methylation is commonly measured using bisulfite sequencing (BS-seq). The quality of a BS-seq library is measured by its bisulfite conversion efficiency. Libraries with low conversion rates are typically excluded from analysis resulting in reduced coverage and increased costs. RESULTS: We have developed a probabilistic method and software, LuxRep, that implements a general linear model and simultaneously accounts for technical replicates (libraries from the same biological sample) from different bisulfite-converted DNA libraries. Using simulations and actual DNA methylation data, we show that including technical replicates with low bisulfite conversion rates generates more accurate estimates of methylation levels and differentially methylated sites. Moreover, using variational inference speeds up computation time necessary for whole genome analysis. CONCLUSIONS: In this work we show that taking into account technical replicates (i.e. libraries) of BS-seq data of varying bisulfite conversion rates, with their corresponding experimental parameters, improves methylation level estimation and differential methylation detection.

Peripheral blood DNA methylation differences in twin pairs discordant for Alzheimer's disease.

BACKGROUND: Alzheimer's disease results from a neurodegenerative process that starts well before the diagnosis can be made. New prognostic or diagnostic markers enabling early intervention into the disease process would be highly valuable. Environmental and lifestyle factors largely modulate the disease risk and may influence the pathogenesis through epigenetic mechanisms, such as DNA methylation. As environmental and lifestyle factors may affect multiple tissues of the body, we hypothesized that the disease-associated DNA methylation signatures are detectable in the peripheral blood of discordant twin pairs. RESULTS: Comparison of 23 disease discordant Finnish twin pairs with reduced representation bisulfite sequencing revealed peripheral blood DNA methylation differences in 11 genomic regions with at least 15.0% median methylation difference and FDR adjusted p value ≤ 0.05. Several of the affected genes are primarily associated with neuronal functions and pathologies and do not display disease-associated differences in gene expression in blood. The DNA methylation mark in ADARB2 gene was found to be differentially methylated also in the anterior hippocampus, including entorhinal cortex, of non-twin cases and controls. Targeted bisulfite pyrosequencing of the DNA methylation mark in ADARB2 gene in 62 Finnish and Swedish twin pairs revealed that, in addition to the disease status, DNA methylation of this region is influenced by gender, age, zygosity, APOE genotype, and smoking. Further analysis of 120 Swedish twin pairs indicated that this specific DNA methylation mark is not predictive for Alzheimer's disease and becomes differentially methylated after disease onset. CONCLUSIONS: DNA methylation differences can be detected in the peripheral blood of twin pairs discordant for Alzheimer's disease. These DNA methylation signatures may have value as disease markers and provide insights into the molecular mechanisms of pathogenesis. We found no evidence that the DNA methylation marks would be associated with gene expression in blood. Further studies are needed to elucidate the potential importance of the associated genes in neuronal functions and to validate the prognostic or diagnostic value of the individual marks or marker panels.

RNA Polymerase III Subunit POLR3G Regulates Specific Subsets of PolyA+ and SmallRNA Transcriptomes and Splicing in Human Pluripotent Stem Cells.

POLR3G is expressed at high levels in human pluripotent stem cells (hPSCs) and is required for maintenance of stem cell state through mechanisms not known in detail. To explore how POLR3G regulates stem cell state, we carried out deep-sequencing analysis of polyA+ and smallRNA transcriptomes present in hPSCs and regulated in POLR3G-dependent manner. Our data reveal that POLR3G regulates a specific subset of the hPSC transcriptome, including multiple transcript types, such as protein-coding genes, long intervening non-coding RNAs, microRNAs and small nucleolar RNAs, and affects RNA splicing. The primary function of POLR3G is in the maintenance rather than repression of transcription. The majority of POLR3G polyA+ transcriptome is regulated during differentiation, and the key pluripotency factors bind to the promoters of at least 30% of the POLR3G-regulated transcripts. Among the direct targets of POLR3G, POLG is potentially important in sustaining stem cell status in a POLR3G-dependent manner.

Mature Let-7 miRNAs fine tune expression of LIN28B in pluripotent human embryonic stem cells.

MicroRNAs (miRNA) are central regulators of diverse biological processes and are important in the regulation of stem cell self-renewal. One of the widely studied miRNA-protein regulators is the Lin28-Let-7 pair. In this study, we demonstrate that contrary to the well-established models of mouse ES cells (mESC) and transformed human cancer cells, the pluripotent state of human ES cells (hESC) involves expression of mature Let-7 family miRNAs with concurrent expression of all LIN28 proteins. We show that mature Let-7 miRNAs are regulated during hESC differentiation and have opposite expression profile with LIN28B. Moreover, mature Let-7 miRNAs fine tune the expression levels of LIN28B protein in pluripotent hESCs, whereas silencing of LIN28 proteins have no effect on mature Let-7 levels. These results bring novel information to the highly complex network of human pluripotency and suggest that maintenance of hESC pluripotency differs greatly from the mESCs in regard to LIN28-Let-7 regulation.

Epigenetic Silencing of the Key Antioxidant Enzyme Catalase in Karyotypically Abnormal Human Pluripotent Stem Cells.

Epigenomic regulation is likely to be important in the maintenance of genomic integrity of human pluripotent stem cells, however, the mechanisms are unknown. We explored the epigenomes and transcriptomes of human pluripotent stem cells before and after spontaneous transformation to abnormal karyotypes and in correlation to cancer cells. Our results reveal epigenetic silencing of Catalase, a key regulator of oxidative stress and DNA damage control in abnormal cells. Our findings provide novel insight into the mechanisms associated with spontaneous transformation of human pluripotent stem cells towards malignant fate. The same mechanisms may control the genomic stability of cells in somatic tissues.