Inactive or moderately active human promoters are enriched for inter-individual epialleles

Epialleles and epiallelic heterogeneity in hematological malignancies

DNA methylation has a well-established role in the pathogenesis, prognosis, and response to treatment in all the spectra of hematological malignancies. However, most of the data reported involve average DNA methylation observed in a sample. The emergence of bisulfite sequencing methods such as enhanced reduced representation that permit analyze adjacent CpGs led to exciting findings. Among these are the epialleles shift and the resulting epigenetic heterogeneity observed in leukemias and lymphomas. Epialleles seem to have an influential role as the cause of mutations that characterize leukemias, may stratify groups with different prognosis and response to treatment, and may be redistributed in the genome at different time points of the disease promoting activation of alternate transcriptional networks. Epiallelic shift may be responsible for the intratumor heterogeneity observed within the cells of the same tumor which increases with disease aggressiveness. It may also responsible for the interpatient heterogeneity explaining why blood cancers exhibit different behavior among different patients. Understanding better epiallelic conformation and the consequent chromatin conformational changes and the pathways that may be affected will permit deeper understanding of hematological malignancies pathogenesis and treatment.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Superovulation alters the expression of endometrial genes critical to tissue remodeling and placentation

PURPOSE

Epidemiologic data suggest that in vitro fertilization (IVF) is associated with an increased risk of disorders of placentation including preeclampsia and fetal growth restriction. Specifically, studies have demonstrated that singleton pregnancies conceived following a fresh embryo transfer are at an increased risk of delivering an infant with low birth weight compared to those conceived following a frozen embryo transfer. The mechanism responsible for this association remains unclear. Procedures utilized in IVF have also been linked with epigenetic changes and gene expression changes in both fetal and maternal tissues. Data suggest that modifications in the maternal endometrium can lead to disordered trophoblast invasion and placentation. This study examines the effect of ovarian stimulation on endometrial gene expression and DNA methylation during the window of implantation to examine potential pathways playing a role in the adverse outcomes associated with IVF.

METHODS

Endometrial biopsies were obtained from oocyte donors and age-matched naturally cycling women 11 days following oocyte retrieval in donors or 12 days following luteinizing hormone (LH) surge in naturally cycling women. Global gene expression was analyzed via Affymetrix Human Gene 1.1 ST array and confirmed with RT-qPCR. DNA methylation was assessed with the Infinium DNA methylation 450 K BeadChip.

RESULTS

Analysis of endometrial gene expression from 23 women (11 oocyte donors and 12 controls) demonstrated 165 genes with a greater than twofold change in expression between donors and controls. While there were 785 genes with significant differential methylation in the endometrium of donors when compared with control subjects, none of the genes with altered expression showed significant changes in DNA methylation. Analysis of the differentially expressed genes showed enrichment for genes involved in endometrial remodeling including PLAT, HSPE2, MMP2, and TIMP1. Validation studies using RT-qPCR found a 73% reduction in expression of heparanase 2 (HSPE2) an enzyme associated with both angiogenesis and cell invasion, a greater than twofold increase in tissue-type plasminogen activator (PLAT), a serine protease participating in matrix degradation, and a 70% increase in MMP2, a gelatinase involved in collagen and fibronectin breakdown.

CONCLUSIONS

Superovulation alters expression of genes critical to endometrial remodeling during early implantation. Such changes could lead to altered trophoblast migration and impaired endovascular invasion. These findings offer a potential mechanism for the adverse perinatal outcomes observed following embryo transfer during fresh IVF cycles.

Sixty-five years of the long march in protein secondary structure prediction: the final stretch?

Protein secondary structure prediction began in 1951 when Pauling and Corey predicted helical and sheet conformations for protein polypeptide backbone even before the first protein structure was determined. Sixty-five years later, powerful new methods breathe new life into this field. The highest three-state accuracy without relying on structure templates is now at 82-84%, a number unthinkable just a few years ago. These improvements came from increasingly larger databases of protein sequences and structures for training, the use of template secondary structure information and more powerful deep learning techniques. As we are approaching to the theoretical limit of three-state prediction (88-90%), alternative to secondary structure prediction (prediction of backbone torsion angles and Cα-atom-based angles and torsion angles) not only has more room for further improvement but also allows direct prediction of three-dimensional fragment structures with constantly improved accuracy. About 20% of all 40-residue fragments in a database of 1199 non-redundant proteins have <6 Å root-mean-squared distance from the native conformations by SPIDER2. More powerful deep learning methods with improved capability of capturing long-range interactions begin to emerge as the next generation of techniques for secondary structure prediction. The time has come to finish off the final stretch of the long march towards protein secondary structure prediction.

Neuroimaging Epigenetics: Challenges and Recommendations for Best Practices

Neuroimaging epigenetics is an interdisciplinary application of epigenetics to cognitive neuroscience that seeks to identify molecular and neural predictors of human behavior. This approach can be sensitive to the dynamic interaction between biological predisposition and environmental influences, and is potentially more informative than an approach using static genetic code. Recent work in this field has generated considerable enthusiasm, yet caution is warranted since any novel cross-disciplinary approach lacks a set of established conventions or standards. In this paper we review existing research in the field of imaging epigenetics, outline important caveats and considerations, and suggest a set of guidelines for researchers conducting this work.

Genome-scale characterization of RNA tertiary structures and their functional impact by RNA solvent accessibility prediction

As most RNA structures are elusive to structure determination, obtaining solvent accessible surface areas (ASAs) of nucleotides in an RNA structure is an important first step to characterize potential functional sites and core structural regions. Here, we developed RNAsnap, the first machine-learning method trained on protein-bound RNA structures for solvent accessibility prediction. Built on sequence profiles from multiple sequence alignment (RNAsnap-prof), the method provided robust prediction in fivefold cross-validation and an independent test (Pearson correlation coefficients, r, between predicted and actual ASA values are 0.66 and 0.63, respectively). Application of the method to 6178 mRNAs revealed its positive correlation to mRNA accessibility by dimethyl sulphate (DMS) experimentally measured in vivo (r = 0.37) but not in vitro (r = 0.07), despite the lack of training on mRNAs and the fact that DMS accessibility is only an approximation to solvent accessibility. We further found strong association across coding and noncoding regions between predicted solvent accessibility of the mutation site of a single nucleotide variant (SNV) and the frequency of that variant in the population for 2.2 million SNVs obtained in the 1000 Genomes Project. Moreover, mapping solvent accessibility of RNAs to the human genome indicated that introns, 5' cap of 5' and 3' cap of 3' untranslated regions, are more solvent accessible, consistent with their respective functional roles. These results support conformational selections as the mechanism for the formation of RNA-protein complexes and highlight the utility of genome-scale characterization of RNA tertiary structures by RNAsnap. The server and its stand-alone downloadable version are available at http://sparks-lab.org.

Genetic and epigenetic variation among inbred mouse littermates: identification of inter-individual differentially methylated regions

Background

Phenotypic variability among inbred littermates reared in controlled environments remains poorly understood. Metastable epialleles refer to loci that intrinsically behave in this way and a few examples have been described. They display differential methylation in association with differential expression. For example, inbred mice carrying the agouti viable yellow (A^vy) allele show a range of coat colours associated with different DNA methylation states at the locus. The availability of next-generation sequencing, in particular whole genome sequencing of bisulphite converted DNA, allows us, for the first time, to search for metastable epialleles at base pair resolution.

Results

Using whole genome bisulphite sequencing of DNA from the livers of five mice from the A^vy colony, we searched for sites at which DNA methylation differed among the mice. A small number of loci, 356, were detected and we call these inter-individual Differentially Methylated Regions, iiDMRs, 55 of which overlap with endogenous retroviral elements (ERVs). Whole genome resequencing of two mice from the colony identified very few differences and these did not occur at or near the iiDMRs. Further work suggested that the majority of ERV iiDMRs are metastable epialleles; the level of methylation was maintained in tissue from other germ layers and the level of mRNA from the neighbouring gene inversely correlated with methylation state. Most iiDMRs that were not overlapping ERV insertions occurred at tissue-specific DMRs and it cannot be ruled out that these are driven by changes in the ratio of cell types in the tissues analysed.

Conclusions

Using the most thorough genome-wide profiling technologies for differentially methylated regions, we find very few intrinsically epigenetically variable regions that we term iiDMRs. The most robust of these are at retroviral elements and appear to be metastable epialleles. The non-ERV iiDMRs cannot be described as metastable epialleles at this stage but provide a novel class of variably methylated elements for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0047-z) contains supplementary material, which is available to authorized users.

Genome-wide DNA methylation map of human neutrophils reveals widespread inter-individual epigenetic variation

The extent of variation in DNA methylation patterns in healthy individuals is not yet well documented. Identification of inter-individual epigenetic variation is important for understanding phenotypic variation and disease susceptibility. Using neutrophils from a cohort of healthy individuals, we generated base-resolution DNA methylation maps to document inter-individual epigenetic variation. We identified 12851 autosomal inter-individual variably methylated fragments (iVMFs). Gene promoters were the least variable, whereas gene body and upstream regions showed higher variation in DNA methylation. The iVMFs were relatively enriched in repetitive elements compared to non-iVMFs, and were associated with genome regulation and chromatin function elements. Further, variably methylated genes were disproportionately associated with regulation of transcription, responsive function and signal transduction pathways. Transcriptome analysis indicates that iVMF methylation at differentially expressed exons has a positive correlation and local effect on the inclusion of that exon in the mRNA transcript.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

BACKGROUND

Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

RESULTS

First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

CONCLUSIONS

The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

The meta-epigenomic structure of purified human stem cell populations is defined at cis-regulatory sequences

The mechanism and significance of epigenetic variability in the same cell type between healthy individuals are not clear. Here, we purify human CD34+ hematopoietic stem and progenitor cells (HSPCs) from different individuals and find that there is increased variability of DNA methylation at loci with properties of promoters and enhancers. The variability is especially enriched at candidate enhancers near genes transitioning between silent and expressed states, and encoding proteins with leukocyte differentiation properties. Our findings of increased variability at loci with intermediate DNA methylation values, at candidate “poised” enhancers, and at genes involved in HSPC lineage commitment suggest that CD34+ cell subtype heterogeneity between individuals is a major mechanism for the variability observed. Epigenomic studies performed on cell populations, even when purified, are testing collections of epigenomes, or meta-epigenomes. Our findings show that meta-epigenomic approaches to data analysis can provide insights into cell subpopulation structure.

H3K4me1 marks DNA regions hypomethylated during aging in human stem and differentiated cells

In differentiated cells, aging is associated with hypermethylation of DNA regions enriched in repressive histone post-translational modifications. However, the chromatin marks associated with changes in DNA methylation in adult stem cells during lifetime are still largely unknown. Here, DNA methylation profiling of mesenchymal stem cells (MSCs) obtained from individuals aged 2 to 92 yr identified 18,735 hypermethylated and 45,407 hypomethylated CpG sites associated with aging. As in differentiated cells, hypermethylated sequences were enriched in chromatin repressive marks. Most importantly, hypomethylated CpG sites were strongly enriched in the active chromatin mark H3K4me1 in stem and differentiated cells, suggesting this is a cell type–independent chromatin signature of DNA hypomethylation during aging. Analysis of scedasticity showed that interindividual variability of DNA methylation increased during aging in MSCs and differentiated cells, providing a new avenue for the identification of DNA methylation changes over time. DNA methylation profiling of genetically identical individuals showed that both the tendency of DNA methylation changes and scedasticity depended on nongenetic as well as genetic factors. Our results indicate that the dynamics of DNA methylation during aging depend on a complex mixture of factors that include the DNA sequence, cell type, and chromatin context involved and that, depending on the locus, the changes can be modulated by genetic and/or external factors.

DNA methylation abnormalities at gene promoters are extensive and variable in the elderly and phenocopy cancer cells

Abnormal patterns of DNA methylation are one of the hallmarks of cancer cells. The process of aging has also been associated with similar, albeit less dramatic, changes in methylation patterns, leading to the hypothesis that age-related changes in DNA methylation may partially underlie the increased risk of cancer in the elderly. Here we studied 377 participants aged 85 yr from the Newcastle 85+ Study to investigate the extent of, and interindividual variation in, age-related changes in DNA methylation at specific CpG islands. Using highly quantitative pyrosequencing analysis, we found extensive and highly variable methylation of promoter-associated CpG islands with levels ranging from 4% to 35%, even at known tumor suppressor genes such as TWIST2. Furthermore, the interindividual differences in methylation seen across this elderly population phenocopies multiple features of the altered methylation patterns seen in cancer cells. Both aging- and cancer-related methylation can occur at similar sets of genes, both result in the formation of densely methylated, and likely transcriptionally repressed, alleles, and both exhibit coordinate methylation across multiple loci. In addition, high methylation levels were associated with subsequent diagnosis of leukemia or lymphoma during a 3-yr follow-up period (P=0.00008). These data suggest that the accumulation of age-related changes in promoter-associated CpG islands may contribute to the increased cancer risk seen during aging.-Gautrey, H. E., van Otterdijk, S. D., Cordell, H. J., Newcastle 85+ study core team, Mathers, J. C., Strathdee, G. DNA methylation abnormalities at gene promoters are extensive and variable in the elderly and phenocopy cancer cells.

Age-associated epigenetic drift: implications, and a case of epigenetic thrift?

It is now well established that the genomic landscape of DNA methylation (DNAm) gets altered as a function of age, a process we here call ‘epigenetic drift’. The biological, functional, clinical and evolutionary significance of this epigenetic drift, however, remains unclear. We here provide a brief review of epigenetic drift, focusing on the potential implications for ageing, stem cell biology and disease risk prediction. It has been demonstrated that epigenetic drift affects most of the genome, suggesting a global deregulation of DNAm patterns with age. A component of this drift is tissue-specific, allowing remarkably accurate age-predictive models to be constructed. Another component is tissue-independent, targeting stem cell differentiation pathways and affecting stem cells, which may explain the observed decline of stem cell function with age. Age-associated increases in DNAm target developmental genes, overlapping those associated with environmental disease risk factors and with disease itself, notably cancer. In particular, cancers and precursor cancer lesions exhibit aggravated age DNAm signatures. Epigenetic drift is also influenced by genetic factors. Thus, drift emerges as a promising biomarker for premature or biological ageing, and could potentially be used in geriatrics for disease risk prediction. Finally, we propose, in the context of human evolution, that epigenetic drift may represent a case of epigenetic thrift, or bet-hedging. In summary, this review demonstrates the growing importance of the ‘ageing epigenome’, with potentially far-reaching implications for understanding the effect of age on stem cell function and differentiation, as well as for disease prevention.