Deciphering the cancer imprintome

AIMER: A SNP-independent software for identifying imprinting-like allelic methylated regions from DNA methylome

Genomic imprinting is essential for mammalian growth and embryogenesis. High-throughput bisulfite sequencing accompanied with parental haplotype-specific information allows analysis of imprinted genes and imprinting control regions (ICRs) on a large scale. Currently, although several allelic methylated regions (AMRs) detection software were developed, methods for detecting imprinted AMRs is still limited. Here, we developed a SNP-independent statistical approach, AIMER, to detect imprinting-like AMRs. By using the mouse frontal cortex methylome as input, we demonstrated that AIMER performs very well in detecting known germline ICRs compared with other methods. Furthermore, we found the putative parental AMRs AIMER detected could be distinguished from sequence-dependent AMRs. Finally, we found a novel germline imprinting-like AMR using WGBS data from 17 distinct mouse tissue samples. The results indicate that AIMER is a good choice for detecting imprinting-like (parent-of-origin-dependent) AMRs. We hope this method will be helpful for future genomic imprinting studies. The Python source code for our project is now publicly available on both GitHub (https://github.com/ZhaoLab-TMU/AIMER) and Gitee (https://gitee.com/zhaolab_tmu/AIMER).

Creation and Validation of the First Infinium DNA Methylation Array for the Human Imprintome

Background

Differentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile these regions, but it is very expensive and it requires not only a high coverage but it is also computationally intensive to assess those regions.

Findings

To address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 9,757 probes map to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites and obtained a mean R 2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array and obtained a mean R 2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799-0.945).

Conclusions

Our custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.

Validation of a DNA methylation microarray for 285,000 CpG sites in the mouse genome

The mouse has been extensively used as the model organism in many studies to characterize biological pathways, drug effects and to mimic human diseases. The similar DNA sequences between both species facilitates these type of experiments. However, much less is known about the mouse epigenome, particularly for DNA methylation. Progress in delivering mouse DNA methylomes has been slow due to the currently available time-consuming and expensive methodologies. Following the great acceptance of the human DNA methylation microarrays, we have herein validated a newly developed DNA methylation microarray (Infinium Mouse Methylation BeadChip) that interrogates 280,754 unique CpG sites within the mouse genome. The CpGs included in the platform cover CpG islands, shores, shelfs and open sea sequences; and loci surrounding transcription start sites and gene bodies. From the functional standpoint, mouse ENCODE representative DNase hypersensitivity sites (rDHSs) and candidate cis-Regulatory Elements (cCREs) are also included. Herein, we show that the profiled mouse DNA methylation microarray provides reliable values among technical replicates; matched results from fresh frozen versus formalin-fixed samples; detects hemimethylated X-chromosome and imprinted CpG sites; and is able to determine CpG methylation changes in mouse cell lines treated with a DNA demethylating agent or upon genetic disruption of a DNA methyltransferase. Most important, using unsupervised hierarchical clustering and t-SNE approaches, the platform is able to classify all types of normal mouse tissues and organs. These data underscore the great features of the assessed microarray to obtain comprehensive DNA methylation profiles of the mouse genome.

Imprinted Genes and Multiple Sclerosis: What Do We Know?

Multiple sclerosis (MS) is a chronic autoimmune neurodegenerative disease of the central nervous system that arises from interplay between non-genetic and genetic risk factors. The epigenetics functions as a link between these factors, affecting gene expression in response to external influence, and therefore should be extensively studied to improve the knowledge of MS molecular mechanisms. Among others, the epigenetic mechanisms underlie the establishment of parent-of-origin effects that appear as phenotypic differences depending on whether the allele was inherited from the mother or father. The most well described manifestation of parent-of-origin effects is genomic imprinting that causes monoallelic gene expression. It becomes more obvious that disturbances in imprinted genes at the least affecting their expression do occur in MS and may be involved in its pathogenesis. In this review we will focus on the potential role of imprinted genes in MS pathogenesis.

ZNF597 is a maternally expressed imprinted gene in the Holstein breed

Genomic imprinting is an epigenetic phenomenon that leads to the preferential expression of genes from either the paternal or maternal allele. Imprinted genes play important roles in mammalian growth and development and a central role in placental function. ZNF597 and NAA60 are two paternally imprinted genes in the human ZNF597-NAA60 imprinted locus, both of which show biallelic expression in the mouse, but their imprinting status in cattle is still unknown. In this study, we examined the allelic expression of ZNF597 and NAA60 in adult bovine placental and somatic tissues. By comparing the mRNA-based genotypes with the genomic DNA-based genotypes, we identified monoallelic expression of ZNF597 in the placenta and in seven other tissues, including the cerebrum, heart, liver, spleen, lung, kidney, and muscle. Nevertheless, analysis revealed biallelic expression of the NAA60 gene in these tissues. Moreover, we tested the imprinting status of ZNF597 and confirmed that the maternal allele is expressed in the bovine placenta. To determine the role of DNA methylation in regulating monoallelic/imprinted expression of bovine ZNF597, the methylation status of two CpG-enriched regions in the bovine ZNF597-NAA60 locus was analyzed using the bisulfite sequencing method. Differentially methylated regions were detected on ten CpG loci in the bovine ZNF597 promoter region. In summary, the bovine ZNF597 gene is a maternally expressed gene, and its expression is regulated by DNA methylation, whereas the NAA60 gene is not imprinted in cattle.

Computational models for lncRNA function prediction and functional similarity calculation

From transcriptional noise to dark matter of biology, the rapidly changing view of long non-coding RNA (lncRNA) leads to deep understanding of human complex diseases induced by abnormal expression of lncRNAs. There is urgent need to discern potential functional roles of lncRNAs for further study of pathology, diagnosis, therapy, prognosis, prevention of human complex disease and disease biomarker detection at lncRNA level. Computational models are anticipated to be an effective way to combine current related databases for predicting most potential lncRNA functions and calculating lncRNA functional similarity on the large scale. In this review, we firstly illustrated the biological function of lncRNAs from five biological processes and briefly depicted the relationship between mutations or dysfunctions of lncRNAs and human complex diseases involving cancers, nervous system disorders and others. Then, 17 publicly available lncRNA function-related databases containing four types of functional information content were introduced. Based on these databases, dozens of developed computational models are emerging to help characterize the functional roles of lncRNAs. We therefore systematically described and classified both 16 lncRNA function prediction models and 9 lncRNA functional similarity calculation models into 8 types for highlighting their core algorithm and process. Finally, we concluded with discussions about the advantages and limitations of these computational models and future directions of lncRNA function prediction and functional similarity calculation. We believe that constructing systematic functional annotation systems is essential to strengthen the prediction accuracy of computational models, which will accelerate the identification process of novel lncRNA functions in the future.

How Does Reprogramming to Pluripotency Affect Genomic Imprinting?

Human induced Pluripotent Stem Cells (hiPSCs) have the capacity to generate a wide range of somatic cells, thus representing an ideal tool for regenerative medicine. Patient-derived hiPSCs are also used for in vitro disease modeling and drug screenings. Several studies focused on the identification of DNA mutations generated, or selected, during the derivation of hiPSCs, some of which are known to drive cancer formation. Avoiding such stable genomic aberrations is paramount for successful use of hiPSCs, but it is equally important to ensure that their epigenetic information is correct, given the critical role of epigenetics in transcriptional regulation and its involvement in a plethora of pathologic conditions. In this review we will focus on genomic imprinting, a prototypical epigenetic mechanism whereby a gene is expressed in a parent-of-origin specific manner, thanks to the differential methylation of specific DNA sequences. Conventional hiPSCs are thought to be in a pluripotent state primed for differentiation. They display a hypermethylated genome with an unexpected loss of DNA methylation at imprinted loci. Several groups recently reported the generation of hiPSCs in a more primitive developmental stage, called naïve pluripotency. Naïve hiPSCs share several features with early human embryos, such as a global genome hypomethylation, which is also accompanied by a widespread loss of DNA methylation at imprinted loci. Given that loss of imprinting has been observed in genetic developmental disorders as well as in a wide range of cancers, it is fundamental to make sure that hiPSCs do not show such epigenetic aberrations. We will discuss what specific imprinted genes, associated with human pathologies, have been found commonly misregulated in hiPSCs and suggest strategies to effectively detect and avoid such undesirable epigenetic abnormalities.

DNA methylation of imprinted gene control regions in the regression of low-grade cervical lesions

The role of host epigenetic mechanisms in the natural history of low-grade cervical intraepithelial neoplasia (CIN1) is not well characterized. We explored differential methylation of imprinted gene regulatory regions as predictors of the risk of CIN1 regression. A total of 164 patients with CIN1 were recruited from 10 Duke University clinics for the CIN Cohort Study. Participants had colposcopies at enrollment and up to five follow-up visits over 3 years. DNA was extracted from exfoliated cervical cells for methylation quantitation at CpG (cytosine-phosphate-guanine) sites and human papillomavirus (HPV) genotyping. Hazard ratios (HR) and 95% confidence intervals (CI) were estimated using Cox regression to quantify the effect of methylation on CIN1 regression over two consecutive visits, compared to non-regression (persistent CIN1; progression to CIN2+; or CIN1 regression at a single time-point), adjusting for age, race, high-risk HPV (hrHPV), parity, oral contraceptive and smoking status. Median participant age was 26.6 years (range: 21.0-64.4 years), 39% were African-American, and 11% were current smokers. Most participants were hrHPV-positive at enrollment (80.5%). Over one-third of cases regressed (n = 53, 35.1%). Median time-to-regression was 12.6 months (range: 4.5-24.0 months). Probability of CIN1 regression was negatively correlated with methylation at IGF2AS CpG 5 (HR = 0.41; 95% CI = 0.23-0.77) and PEG10 DMR (HR = 0.80; 95% CI = 0.65-0.98). Altered methylation of imprinted IGF2AS and PEG10 DMRs may play a role in the natural history of CIN1. If confirmed in larger studies, further research on imprinted gene DMR methylation is warranted to determine its efficacy as a biomarker for cervical cancer screening.

Copy number rather than epigenetic alterations are the major dictator of imprinted methylation in tumors

It has been postulated that imprinting aberrations are common in tumors. To understand the role of imprinting in cancer, we have characterized copy-number and methylation in over 280 cancer cell lines and confirm our observations in primary tumors. Imprinted differentially methylated regions (DMRs) regulate parent-of-origin monoallelic expression of neighboring transcripts in cis. Unlike single-copy CpG islands that may be prone to hypermethylation, imprinted DMRs can either loose or gain methylation during tumorigenesis. Here, we show that methylation profiles at imprinted DMRs often not represent genuine epigenetic changes but simply the accumulation of underlying copy-number aberrations (CNAs), which is independent of the genome methylation state inferred from cancer susceptible loci. Our results reveal that CNAs also influence allelic expression as loci with copy-number neutral loss-of-heterozygosity or amplifications may be expressed from the appropriate parental chromosomes, which is indicative of maintained imprinting, although not observed as a single expression foci by RNA FISH.Altered genomic imprinting is frequently reported in cancer. Here, the authors analyze copy number and methylation in cancer cell lines and primary tumors to show that imprinted methylation profiles represent the accumulation of copy number alteration, rather than epigenetic alterations.

Epigenetic response of imprinted domains during carcinogenesis

Background

Imprinted domains have been identified as targets for aberrant DNA methylation during carcinogenesis, but it remains unclear when these epigenetic alterations occur and how they contribute to tumor progression. Epigenetic instability at key cis-regulatory elements within imprinted domains can concomitantly activate proto-oncogenes and turn off tumor suppressor genes. Thus, to further characterize the epigenetic response of imprinted domains during carcinogenesis, we compared the stability of DNA methylation at a variety of cis-regulatory elements within imprinted domains in two fundamentally different mouse tumors, benign and malignant, induced by the KrasG12D mutation.

Results

We report that imprinted domains remain stable in benign processes but are highly susceptible to epigenetic alterations in infiltrative lesions. The preservation of DNA methylation within imprinted domains in benign tumors throughout their duration suggests that imprinted genes are not involved with the initiation of carcinogenesis or the growth of tumors. However, the frequent detection of DNA methylation changes at imprinting control regions in infiltrative lesions suggest that imprinted genes are associated with tumor cells gaining the ability to defy tissue boundaries.

Conclusion

Overall, this study demonstrates that imprinted domains are targeted for DNA hypermethylation when benign tumor cells transition to malignant. Thus, monitoring DNA methylation within imprinted domains may be useful in evaluating the progression of neoplasms.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0393-8) contains supplementary material, which is available to authorized users.

Epigenetic instability at imprinting control regions in a Kras(G12D)-induced T-cell neoplasm

Although aberrant DNA methylation within imprinted domains has been reported in a variety of neoplastic diseases, it remains largely uncharacterized in the context of carcinogenesis. In this study, we induced T-cell lymphoma in mice by employing a breeding scheme involving mouse strains, LSL-Kras(G12D) and MMTV-Cre. We then systematically surveyed imprinted domains for DNA methylation changes during tumor progression using combined bisulfite restriction analysis and NGS-based bisulfite sequencing. We detected hyper- or hypo-methylation at the imprinting control regions (ICRs) of the Dlk1, Peg10, Peg3, Grb10, and Gnas domains. These DNA methylation changes at ICRs were more prevalent and consistent than those observed at the promoter regions of well-known tumor suppressors, such as Mgmt, Fhit, and Mlh1. Thus, the changes observed at these imprinted domains are the outcome of isolated incidents affecting DNA methylation settings. Within imprinted domains, DNA methylation changes tend to be restricted to ICRs as nearby somatic differentially methylated regions and promoter regions experience no change. Furthermore, detailed analyses revealed that small cis-regulatory elements within ICRs tend to be resistant to DNA methylation changes, suggesting potential protection by unknown trans-factors. Overall, this study demonstrates that DNA methylation changes at ICRs are dynamic during carcinogenesis and advocates that detection of aberrant DNA methylation at ICRs may serve as a biomarker to enhance diagnostic procedures.

Imprinted genes and the environment: links to the toxic metals arsenic, cadmium, lead and mercury

Imprinted genes defy rules of Mendelian genetics with their expression tied to the parent from whom each allele was inherited. They are known to play a role in various diseases/disorders including fetal growth disruption, lower birth weight, obesity, and cancer. There is increasing interest in understanding their influence on environmentally-induced disease. The environment can be thought of broadly as including chemicals present in air, water and soil, as well as food. According to the Agency for Toxic Substances and Disease Registry (ATSDR), some of the highest ranking environmental chemicals of concern include metals/metalloids such as arsenic, cadmium, lead and mercury. The complex relationships between toxic metal exposure, imprinted gene regulation/expression and health outcomes are understudied. Herein we examine trends in imprinted gene biology, including an assessment of the imprinted genes and their known functional roles in the cell, particularly as they relate to toxic metals exposure and disease. The data highlight that many of the imprinted genes have known associations to developmental diseases and are enriched for their role in the TP53 and AhR pathways. Assessment of the promoter regions of the imprinted genes resulted in the identification of an enrichment of binding sites for two transcription factor families, namely the zinc finger family II and PLAG transcription factors. Taken together these data contribute insight into the complex relationships between toxic metals in the environment and imprinted gene biology.

Genome-wide parent-of-origin DNA methylation analysis reveals the intricacies of human imprinting and suggests a germline methylation-independent mechanism of establishment

Differential methylation between the two alleles of a gene has been observed in imprinted regions, where the methylation of one allele occurs on a parent-of-origin basis, the inactive X-chromosome in females, and at those loci whose methylation is driven by genetic variants. We have extensively characterized imprinted methylation in a substantial range of normal human tissues, reciprocal genome-wide uniparental disomies, and hydatidiform moles, using a combination of whole-genome bisulfite sequencing and high-density methylation microarrays. This approach allowed us to define methylation profiles at known imprinted domains at base-pair resolution, as well as to identify 21 novel loci harboring parent-of-origin methylation, 15 of which are restricted to the placenta. We observe that the extent of imprinted differentially methylated regions (DMRs) is extremely similar between tissues, with the exception of the placenta. This extra-embryonic tissue often adopts a different methylation profile compared to somatic tissues. Further, we profiled all imprinted DMRs in sperm and embryonic stem cells derived from parthenogenetically activated oocytes, individual blastomeres, and blastocysts, in order to identify primary DMRs and reveal the extent of reprogramming during preimplantation development. Intriguingly, we find that in contrast to ubiquitous imprints, the majority of placenta-specific imprinted DMRs are unmethylated in sperm and all human embryonic stem cells. Therefore, placental-specific imprinting provides evidence for an inheritable epigenetic state that is independent of DNA methylation and the existence of a novel imprinting mechanism at these loci.

Variable maternal methylation overlapping the nc886/vtRNA2-1 locus is locked between hypermethylated repeats and is frequently altered in cancer

Cancer is as much an epigenetic disease as a genetic one; however, the interplay between these two processes is unclear. Recently, it has been shown that a large proportion of DNA methylation variability can be explained by allele-specific methylation (ASM), either at classical imprinted loci or those regulated by underlying genetic variants. During a recent screen for imprinted differentially methylated regions, we identified the genomic interval overlapping the non-coding nc886 RNA (previously known as vtRNA2-1) as an atypical ASM that shows variable levels of methylation, predominantly on the maternal allele in many tissues. Here we show that the nc886 interval is the first example of a polymorphic imprinted DMR in humans. Further analysis of the region suggests that the interval subjected to ASM is approximately 2 kb in size and somatically acquired. An in depth analysis of this region in primary cancer samples with matching normal adjacent tissue from the Cancer Genome Atlas revealed that aberrant methylation in bladder, breast, colon and lung tumors occurred in approximately 27% of cases. Hypermethylation occurred more frequently than hypomethylation. Using additional normal-tumor paired samples we show that on rare occasions the aberrant methylation profile is due to loss-of-heterozygosity. This work therefore suggests that the nc886 locus is subject to variable allelic methylation that undergoes cancer-associated epigenetic changes in solid tumors.

The human imprintome: regulatory mechanisms, methods of ascertainment, and roles in disease susceptibility

Imprinted genes form a special subset of the genome, exhibiting monoallelic expression in a parent-of-origin-dependent fashion. This monoallelic expression is controlled by parental-specific epigenetic marks, which are established in gametogenesis and early embryonic development and are persistent in all somatic cells throughout life. We define this specific set of cis-acting epigenetic regulatory elements as the imprintome, a distinct and specially tasked subset of the epigenome. Imprintome elements contain DNA methylation and histone modifications that regulate monoallelic expression by affecting promoter accessibility, chromatin structure, and chromatin configuration. Understanding their regulation is critical because a significant proportion of human imprinted genes are implicated in complex diseases. Significant species variation in the repertoire of imprinted genes and their epigenetic regulation, however, will not allow model organisms solely to be used for this crucial purpose. Ultimately, only the human will suffice to accurately define the human imprintome.

Function of lncRNAs and approaches to lncRNA-protein interactions

Long non-coding RNAs (lncRNAs), which represent a new frontier in molecular biology, play important roles in regulating gene expression at epigenetic, transcriptional and post-transcriptional levels. More and more lncRNAs have been found to play important roles in normal cell physiological activities, and participate in the development of varieties of tumors and other diseases. Previously, we have only been able to determine the function of lncRNAs through multiple mechanisms, including genetic imprinting, chromatin remodeling, splicing regulation, mRNA decay, and translational regulation. Application of technological advances to research into the function of lncRNAs is extremely important. The major tools for exploring lncRNAs include microarrays, RNA sequencing (RNA-seq), Northern blotting, real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR), fluorescence in situ hybridization (FISH), RNA interference (RNAi), RNA-binding protein immunoprecipitation (RIP), chromatin isolation by RNA purification (ChIRP), crosslinking-immunopurification (CLIP), and bioinformatic prediction. In this review, we highlight the functions of lncRNAs, and advanced methods to research lncRNA-protein interactions.

Medulloblastoma biology in the post-genomic era

Medulloblastomas, the most common malignant pediatric brain tumors, are comprised of four molecularly distinct subtypes. However, treatment has yet to exploit these molecular vulnerabilities. Three recent studies sequenced a total of 310 primary tumors and identified that two of the four medulloblastoma subtypes are concomitantly associated with subtype-specific mutations as previously characterized. In contrast, the overwhelming majority of mutations occurred only once in the entire cohort and just 12 genes were recurrently mutated with statistical significance. Perturbations in epigenetic regulation are emerging as a unifying theme in cancer and similarly recurring mutations in epigenetic mechanisms were distributed across all subtypes in medulloblastoma. Designing targeted therapies to such a molecularly diverse disease in the post-genomic era presents new challenges. This will require novel methods to link these nonrecurrent mutations into pathways, and preclinical models that faithfully recapitulate patient driver events. Presently, medulloblastoma reinforces epigenetic mechanisms as a tantalizing therapeutic target in cancers.

Leukocyte DNA as surrogate for the evaluation of imprinted Loci methylation in mammary tissue DNA

There is growing interest in identifying surrogate tissues to identify epimutations in cancer patients since primary target tissues are often difficult to obtain. Methylation patterns at imprinted loci are established during gametogenesis and post fertilization and their alterations have been associated with elevated risk of cancer. Methylation at several imprinted differentially methylated regions (GRB10 ICR, H19 ICR, KvDMR, SNRPN/SNURF ICR, IGF2 DMR0, and IGF2 DMR2) were analyzed in DNA from leukocytes and mammary tissue (normal, benign diseases, or malignant tumors) from 87 women with and without breast cancer (average age of cancer patients: 53; range: 31-77). Correlations between genomic variants and DNA methylation at the studied loci could not be assessed, making it impossible to exclude such effects. Methylation levels observed in leukocyte and mammary tissue DNA were close to the 50% expected for monoallellic methylation. While no correlation was observed between leukocyte and mammary tissue DNA methylation for most of the analyzed imprinted genes, Spearman's correlations were statistically significant for IGF2 DMR0 and IGF2 DMR2, although absolute methylation levels differed. Leukocyte DNA methylation levels of selected imprinted genes may therefore serve as surrogate markers of DNA methylation in cancer tissue.

Imprinting at the PLAGL1 domain is contained within a 70-kb CTCF/cohesin-mediated non-allelic chromatin loop

Paternal duplications of chromosome 6q24, a region that contains the imprinted PLAGL1 and HYMAI transcripts, are associated with transient neonatal diabetes mellitus. A common feature of imprinted genes is that they tend to cluster together, presumably as a result of sharing common cis-acting regulatory elements. To determine the extent of this imprinted cluster in human and mouse, we have undertaken a systematic analysis of allelic expression and DNA methylation of the genes mapping within an ∼1.4-Mb region flanking PLAGL1/Plagl1. We confirm that all nine neighbouring genes are biallelically expressed in both species. In human we identify two novel paternally expressed PLAGL1 coding transcripts that originate from unique promoter regions. Chromatin immunoprecipitation for CTCF and the cohesin subunits RAD21 and SMC3 reveals evolutionarily conserved binding sites within unmethylated regions ∼5 kb downstream of the PLAGL1 differentially methylated region and within the PLAGL1 3' untranslated region (UTR). Higher-order chromatin looping occurs between these regions in both expressing and non-expressing tissues, forming a non-allelic chromatin loop around the PLAGL1/Plagl1 gene. In placenta and brain tissues, we identify an additional interaction between the PLAGL1 P3/P4 promoters and the unmethylated element downstream of the PLAGL1 differentially methylated region that we propose facilitates imprinted expression of these alternative isoforms.

Dynamic DNA methylation across diverse human cell lines and tissues

As studies of DNA methylation increase in scope, it has become evident that methylation has a complex relationship with gene expression, plays an important role in defining cell types, and is disrupted in many diseases. We describe large-scale single-base resolution DNA methylation profiling on a diverse collection of 82 human cell lines and tissues using reduced representation bisulfite sequencing (RRBS). Analysis integrating RNA-seq and ChIP-seq data illuminates the functional role of this dynamic mark. Loci that are hypermethylated across cancer types are enriched for sites bound by NANOG in embryonic stem cells, which supports and expands the model of a stem/progenitor cell signature in cancer. CpGs that are hypomethylated across cancer types are concentrated in megabase-scale domains that occur near the telomeres and centromeres of chromosomes, are depleted of genes, and are enriched for cancer-specific EZH2 binding and H3K27me3 (repressive chromatin). In noncancer samples, there are cell-type specific methylation signatures preserved in primary cell lines and tissues as well as methylation differences induced by cell culture. The relationship between methylation and expression is context-dependent, and we find that CpG-rich enhancers bound by EP300 in the bodies of expressed genes are unmethylated despite the dense gene-body methylation surrounding them. Non-CpG cytosine methylation occurs in human somatic tissue, is particularly prevalent in brain tissue, and is reproducible across many individuals. This study provides an atlas of DNA methylation across diverse and well-characterized samples and enables new discoveries about DNA methylation and its role in gene regulation and disease.

Characterization of novel paternal ncRNAs at the Plagl1 locus, including Hymai, predicted to interact with regulators of active chromatin

Genomic imprinting is a complex epigenetic mechanism of transcriptional control that utilizes DNA methylation and histone modifications to bring about parent-of-origin specific monoallelic expression in mammals. Genes subject to imprinting are often organised in clusters associated with large non-coding RNAs (ncRNAs), some of which have cis-regulatory functions. Here we have undertaken a detailed allelic expression analysis of an imprinted domain on mouse proximal chromosome 10 comprising the paternally expressed Plagl1 gene. We identified three novel Plagl1 transcripts, only one of which contains protein-coding exons. In addition, we characterised two unspliced ncRNAs, Hymai, the mouse orthologue of HYMAI, and Plagl1it (Plagl1 intronic transcript), a transcript located in intron 5 of Plagl1. Imprinted expression of these novel ncRNAs requires DNMT3L-mediated maternal DNA methylation, which is also indispensable for establishing the correct chromatin profile at the Plagl1 DMR. Significantly, the two ncRNAs are retained in the nucleus, consistent with a potential regulatory function at the imprinted domain. Analysis with catRAPID, a protein-ncRNA association prediction algorithm, suggests that Hymai and Plagl1it RNAs both have potentially high affinity for Trithorax chromatin regulators. The two ncRNAs could therefore help to protect the paternal allele from DNA methylation by attracting Trithorax proteins that mediate H3 lysine-4 methylation.

Submitted GenBank nucleotides sequences:

Plagl1it: JN595789

Hymai: JN595790

Genomic landscape of human allele-specific DNA methylation

DNA methylation mediates imprinted gene expression by passing an epigenomic state across generations and differentially marking specific regulatory regions on maternal and paternal alleles. Imprinting has been tied to the evolution of the placenta in mammals and defects of imprinting have been associated with human diseases. Although recent advances in genome sequencing have revolutionized the study of DNA methylation, existing methylome data remain largely untapped in the study of imprinting. We present a statistical model to describe allele-specific methylation (ASM) in data from high-throughput short-read bisulfite sequencing. Simulation results indicate technical specifications of existing methylome data, such as read length and coverage, are sufficient for full-genome ASM profiling based on our model. We used our model to analyze methylomes for a diverse set of human cell types, including cultured and uncultured differentiated cells, embryonic stem cells and induced pluripotent stem cells. Regions of ASM identified most consistently across methylomes are tightly connected with known imprinted genes and precisely delineate the boundaries of several known imprinting control regions. Predicted regions of ASM common to multiple cell types frequently mark noncoding RNA promoters and represent promising starting points for targeted validation. More generally, our model provides the analytical complement to cutting-edge experimental technologies for surveying ASM in specific cell types and across species.

The evolution of mammalian genomic imprinting was accompanied by the acquisition of novel CpG islands

Parent-of-origin–dependent expression of imprinted genes is mostly associated with allele-specific DNA methylation of the CpG islands (CGIs) called germ line differentially methylated regions (gDMRs). Although the essential role of gDMRs for genomic imprinting has been well established, little is known about how they evolved. In several imprinted loci, the CGIs forming gDMRs may have emerged with the insertion of a retrotransposon or retrogene. To examine the generality of the hypothesis that the CGIs forming gDMRs were novel CGIs recently acquired during mammalian evolution, we reviewed the time of novel CGI emergence for all the maternal gDMR loci using the novel data analyzed in this study combined with the data from previous reports. The comparative sequence analyses using mouse, human, dog, cow, elephant, tammar, opossum, platypus, and chicken genomic sequences were carried out for Peg13, Meg1/Grb10, Plagl1/Zac1, Gnas, and Slc38a4 imprinted loci to obtain comprehensive results. The combined data showed that emergence of novel CGIs occurred universally in the maternal gDMR loci at various time points during mammalian evolution. Furthermore, the analysis of Meg1/Grb10 locus provided evidence that gradual base pair–wise sequence change was involved in the accumulation of CpG sequence, suggesting the mechanism of novel CGI emergence is more complex than the suggestion that CpG sequences originated solely by insertion of CpG-rich transposable elements. We propose that acquisition of novel CGIs was a key genomic change for the evolution of imprinting and that it usually occurred in the maternal gDMR loci.

Pathology from evolutionary conflict, with a theory of X chromosome versus autosome conflict over sexually antagonistic traits

Evolutionary conflicts cause opponents to push increasingly hard and in opposite directions on the regulation of traits. One can see only the intermediate outcome from the balance of the exaggerated and opposed forces. Intermediate expression hides the underlying conflict, potentially misleading one to conclude that trait regulation is designed to achieve efficient and robust expression, rather than arising by the precarious resolution of conflict. Perturbation often reveals the underlying nature of evolutionary conflict. Upon mutation or knockout of one side in the conflict, the other previously hidden and exaggerated push on the trait may cause extreme, pathological expression. In this regard, pathology reveals hidden evolutionary design. We first review several evolutionary conflicts between males and females, including conflicts over mating, fertilization, and the growth rate of offspring. Perturbations of these conflicts lead to infertility, misregulated growth, cancer, behavioral abnormalities, and psychiatric diseases. We then turn to antagonism between the sexes over traits present in both males and females. For many traits, the different sexes favor different phenotypic values, and constraints prevent completely distinct expression in the sexes. In this case of sexual antagonism, we present a theory of conflict between X-linked genes and autosomal genes. We suggest that dysregulation of the exaggerated conflicting forces between the X chromosome and the autosomes may be associated with various pathologies caused by extreme expression along the male-female axis. Rapid evolution of conflicting X-linked and autosomal genes may cause divergence between populations and speciation.

A DNA methylation fingerprint of 1628 human samples

Most of the studies characterizing DNA methylation patterns have been restricted to particular genomic loci in a limited number of human samples and pathological conditions. Herein, we present a compromise between an extremely comprehensive study of a human sample population with an intermediate level of resolution of CpGs at the genomic level. We obtained a DNA methylation fingerprint of 1628 human samples in which we interrogated 1505 CpG sites. The DNA methylation patterns revealed show this epigenetic mark to be critical in tissue-type definition and stemness, particularly around transcription start sites that are not within a CpG island. For disease, the generated DNA methylation fingerprints show that, during tumorigenesis, human cancer cells underwent a progressive gain of promoter CpG-island hypermethylation and a loss of CpG methylation in non-CpG-island promoters. Although transformed cells are those in which DNA methylation disruption is more obvious, we observed that other common human diseases, such as neurological and autoimmune disorders, had their own distinct DNA methylation profiles. Most importantly, we provide proof of principle that the DNA methylation fingerprints obtained might be useful for translational purposes by showing that we are able to identify the tumor type origin of cancers of unknown primary origin (CUPs). Thus, the DNA methylation patterns identified across the largest spectrum of samples, tissues, and diseases reported to date constitute a baseline for developing higher-resolution DNA methylation maps and provide important clues concerning the contribution of CpG methylation to tissue identity and its changes in the most prevalent human diseases.

Methylation screening of reciprocal genome-wide UPDs identifies novel human-specific imprinted genes

Nuclear transfer experiments undertaken in the mid-80's revealed that both maternal and paternal genomes are necessary for normal development. This is due to genomic imprinting, an epigenetic mechanism that results in parent-of-origin monoallelic expression of genes regulated by germline-derived allelic methylation. To date, ∼100 imprinted transcripts have been identified in mouse, with approximately two-thirds showing conservation in humans. It is currently unknown how many imprinted genes are present in humans, and to what extent these transcripts exhibit human-specific imprinted expression. This is mainly due to the fact that the majority of screens for imprinted genes have been undertaken in mouse, with subsequent analysis of the human orthologues. Utilizing extremely rare reciprocal genome-wide uniparental disomy samples presenting with Beckwith-Wiedemann and Silver-Russell syndrome-like phenotypes, we analyzed ∼0.1% of CpG dinculeotides present in the human genome for imprinted differentially methylated regions (DMRs) using the Illumina Infinium methylation27 BeadChip microarray. This approach identified 15 imprinted DMRs associated with characterized imprinted domains, and confirmed the maternal methylation of the RB1 DMR. In addition, we discovered two novel DMRs, first, one maternally methylated region overlapping the FAM50B promoter CpG island, which results in paternal expression of this retrotransposon. Secondly, we found a paternally methylated, bidirectional repressor located between maternally expressed ZNF597 and NAT15 genes. These three genes are biallelically expressed in mice due to lack of differential methylation, suggesting that these genes have become imprinted after the divergence of mouse and humans.

The functional role of long non-coding RNA in human carcinomas

Long non-coding RNAs (lncRNAs) are emerging as new players in the cancer paradigm demonstrating potential roles in both oncogenic and tumor suppressive pathways. These novel genes are frequently aberrantly expressed in a variety of human cancers, however the biological functions of the vast majority remain unknown. Recently, evidence has begun to accumulate describing the molecular mechanisms by which these RNA species function, providing insight into the functional roles they may play in tumorigenesis. In this review, we highlight the emerging functional role of lncRNAs in human cancer.

Disruption of genomic neighbourhood at the imprinted IGF2-H19 locus in Beckwith-Wiedemann syndrome and Silver-Russell syndrome

Hyper- and hypomethylation at the IGF2-H19 imprinting control region (ICR) result in reciprocal changes in IGF2-H19 expression and the two contrasting growth disorders, Beckwith-Wiedemann syndrome (BWS) and Silver-Russell syndrome (SRS). DNA methylation of the ICR controls the reciprocal imprinting of IGF2 and H19 by preventing the binding of the insulator protein, CTCF. We here show that local changes in histone modifications and CTCF--cohesin binding at the ICR in BWS and SRS together with DNA methylation correlate with the higher order chromatin structure at the locus. In lymphoblastoid cells from control individuals, we found the repressive histone H3K9me3 and H4K20me3 marks associated with the methylated paternal ICR allele and the bivalent H3K4me2/H3K27me3 mark together with H3K9ac and CTCF--cohesin associated with the non-methylated maternal allele. In patient-derived cell lines, the mat/pat asymmetric distribution of these epigenetic marks was lost with H3K9me3 and H4K20me3 becoming biallelic in the BWS and H3K4me2, H3K27me3 and H3K9ac together with CTCF-cohesin becoming biallelic in the SRS. We further show that in BWS and SRS cells, there is opposing chromatin looping conformation mediated by CTCF--cohesin binding sites surrounding the locus. In normal cells, lack of CTCF--cohesin binding at the paternal ICR is associated with monoallelic interaction between two CTCF sites flanking the locus. CTCF--cohesin binding at the maternal ICR blocks this interaction by associating with the CTCF site downstream of the enhancers. The two alternative chromatin conformations are differently favoured in BWS and SRS likely predisposing the locus to the activation of IGF2 or H19, respectively.