Mother-child transmission of epigenetic information by tunable polymorphic imprinting

Regulation of nc886 (vtRNA2-1) RNAs is associated with cardiometabolic risk factors and diseases

Non-coding 886 (nc886, vtRNA2-1) is a polymorphically imprinted gene. The methylation status of this locus has been shown to be associated with periconceptional conditions, and both the methylation status and the levels of nc886 RNAs have been shown to associate with later-life health traits. We have previously shown that nc886 RNA levels are associated not only with the methylation status of the locus, but also with a genetic polymorphism upstream from the locus. In this study, we describe the genetic and epigenetic regulators that predict lifelong nc886 RNA levels, as well as their association with cardiometabolic disease (CMD) risk factors and events. We utilised six population cohorts and one CMD cohort comprising 9058 individuals in total. The association of nc886 RNA levels, as predicted by epigenetic and genetic regulators, with CMD phenotypes was analysed using regression models, with a meta-analysis of the results. The meta-analysis showed that individuals with upregulated nc886 RNA levels have higher diastolic blood pressure (β = 0.07, p = 0.008), lower HDL levels (β =  - 0.07, p = 0.006) and an increased incidence of type 2 diabetes (OR = 1.260, p = 0.013). Moreover, CMD patients with upregulated nc886 RNA levels have an increased incidence of stroke (OR = 1.581, p = 0.006) and death (OR = 1.290, p = 0.046). In conclusion, we show that individuals who are predicted to present elevated nc886 RNA levels have poorer cardiovascular health and are at an elevated risk of complications in secondary prevention. This unique mechanism yields metabolic variation in human populations, constituting a CMD risk factor that cannot be modified through lifestyle choices.

Multi-Omics Study Reveals Nc886/vtRNA2-1 as a Positive Regulator of Prostate Cancer Cell Immunity

Noncoding RNA 886 has emerged as a pivotal regulatory RNA with distinct functions across tissues, acting as a regulator of protein activity by directly binding to protein partners. While it is well recognized as a tumor suppressor in prostate cancer, the underlying molecular mechanisms remain elusive. To discover the principal pathways regulated by nc886 in prostate cancer, we used a transcriptomic and proteomic approach analyzing malignant DU145, LNCaP, PC3, and benign RWPE-1 prostate cell line models transiently transfected with in vitro transcribed nc886 or antisense oligonucleotides. Multiomics revelead a significant enrichment of immune system-related pathways across the cell lines, including cytokines and interferon signaling. The interferon response provoked by nc886 was validated by functional assays. The invariability of PKR phosphorylation and NF-κB activity in the gain/loss of nc886 function experiments and the positive regulation of innate immunity suggest a PKR-independent mechanism of nc886 action. Accordingly, the GSEA of the PRAD-TCGA data set revealed immune stimulation as the nc886 most associated node also in the clinical setting. Our study showed that the reduction of nc886 leads to a blunted immune response in prostate cancer.

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

The mystique of epigenetic regulation: the remarkable case of a human noncoding RNA, nc886

nc886 is a regulatory noncoding RNA that is transcribed by RNA polymerase III (Pol III), is variably expressed in different biological contexts, and plays roles in inflammation and cancer. Epigenetic mechanisms play an intriguing role in regulating nc886 expression. As a maternally imprinted gene and metastable epiallele, nc866 exhibits polymorphic imprinting, with a methylation status that is influenced by environmental and biological factors. Consequently, the promoter DNA methylation status and the different resulting RNA expression levels of nc886 are associated with physiological and pathological conditions. In this review, we summarize the literature and explore the significance in relation to diverse roles of nc886.

Expression of the non-coding RNA nc886 facilitates the development of tyrosine kinase inhibitor resistance in EGFR-mutated non-small-cell lung cancer cells

Treatment of non-small-cell lung cancer (NSCLC) patients possessing EGFR-activating mutations with tyrosine kinase inhibitors (TKIs) can confer an initial promising response. However, TKI resistance inevitably arises. Numerous TKI resistance mechanisms are identified including EGFR secondary mutations, bypass receptor tyrosine kinase (RTK) signaling, and cellular transition e.g. epithelial-mesenchymal transition (EMT). To increase the knowledge of TKI resistance we performed an epigenetic screen to identify small non-coding (nc) genes with DNA methylation alterations in HCC827 NSCLC EGFR-mutated cells with acquired TKI resistance. We analyzed Infinium Methylation EPIC 850K Array data for DNA methylation changes present in both TKI-resistant HCC827 cells with EMT and MET-amplification. Hereby, we identified that the polymorphic maternal imprinted gene nc886 (vtRNA2-1) has a decrease in promoter DNA methylation in TKI-resistant cells. This epigenetic change was associated with an increase in the expression of nc886. The induction of EMT did not affect nc886 expression. CRISPR/Cas9-mediated distortion of the nc886 sequence increased the sensitivity of HCC827 cells towards TKI. Finally, nc886 sequence distortion hindered MET RTK activation and instead was EMT the endpoint TKI resistance mechanism. In conclusion, the expression of nc886 contributes to TKI resistance in the HCC827 NSCLC cell line by supporting cell survival and selection of the endpoint TKI resistance mechanism. We propose DNA methylation and expression changes for nc886 to constitute a novel TKI resistance contributing mechanism in NSCLC.

Epigenome-wide analysis across the development span of pediatric acute lymphoblastic leukemia: backtracking to birth

BACKGROUND

Cancer is the leading cause of disease-related mortality in children. Causes of leukemia, the most common form, are largely unknown. Growing evidence points to an origin in-utero, when global redistribution of DNA methylation occurs driving tissue differentiation.

METHODS

Epigenome-wide DNA methylation was profiled in surrogate (blood) and target (bone marrow) tissues at birth, diagnosis, remission and relapse of pediatric pre-B acute lymphoblastic leukemia (pre-B ALL) patients. Double-blinded analyses was performed between prospective cohorts extending from birth to diagnosis and retrospective studies backtracking from clinical disease to birth. Validation was carried out using independent technologies and populations.

RESULTS

The imprinted and immuno-modulating VTRNA2-1 was hypermethylated (FDR<0.05) at birth in nested cases relative to controls in all tested populations (totaling 317 cases and 483 controls), including European and Hispanic ancestries. VTRNA2-1 methylation was stable over follow-up years after birth and across surrogate, target and other tissues (n=5,023 tissues; 30 types). When profiled in leukemic tissues from two clinical cohorts (totaling 644 cases), VTRNA2-1 methylation exhibited higher levels at diagnosis relative to controls, it reset back to normal levels at remission, and then re-increased to above control levels at relapse. Hypermethylation was significantly associated with worse pre-B ALL patient survival and with reduced VTRNA2-1 expression (n=2,294 tissues; 26 types), supporting a functional and translational role for VTRNA2-1 methylation.

CONCLUSION

This study provides proof-of-concept to detect at birth epigenetic precursors of pediatric pre-B ALL. These alterations were reproducible with different technologies, in three continents and in two ethnicities, and can offer biomarkers for early detection and prognosis as well as actionable targets for therapy.

Epigenetics of the non-coding RNA nc886 across blood, adipose tissue and skeletal muscle in offspring exposed to diabetes in pregnancy

BACKGROUND

Diabetes in pregnancy is associated with increased risk of long-term metabolic disease in the offspring, potentially mediated by in utero epigenetic variation. Previously, we identified multiple differentially methylated single CpG sites in offspring of women with gestational diabetes mellitus (GDM), but whether stretches of differentially methylated regions (DMRs) can also be identified in adolescent GDM offspring is unknown. Here, we investigate which DNA regions in adolescent offspring are differentially methylated in blood by exposure to diabetes in pregnancy. The secondary aim was to characterize the RNA expression of the identified DMR, which contained the nc886 non-coding RNA.

METHODS

To identify DMRs, we employed the bump hunter method in samples from young (9-16 yr, n = 92) offspring of women with GDM (O-GDM) and control offspring (n = 94). Validation by pyrosequencing was performed in an adult offspring cohort (age 28-33 years) consisting of O-GDM (n = 82), offspring exposed to maternal type 1 diabetes (O-T1D, n = 67) and control offspring (O-BP, n = 57). RNA-expression was measured using RT-qPCR in subcutaneous adipose tissue and skeletal muscle.

RESULTS

One significant DMR represented by 10 CpGs with a bimodal methylation pattern was identified, located in the nc886/VTRNA2-1 non-coding RNA gene. Low methylation status across all CpGs of the nc886 in the young offspring was associated with maternal GDM. While low methylation degree in adult offspring in blood, adipose tissue, and skeletal muscle was not associated with maternal GDM, adipose tissue nc886 expression was increased in O-GDM compared to O-BP, but not in O-T1D. In addition, adipose tissue nc886 expression levels were positively associated with maternal pre-pregnancy BMI (p = 0.006), but not with the offspring's own adiposity.

CONCLUSIONS

Our results highlight that nc886 is a metastable epiallele, whose methylation in young offspring is negatively correlated with maternal obesity and GDM status. The physiological effect of nc886 may be more important in adipose tissue than in skeletal muscle. Further research should aim to investigate how nc886 regulation in adipose tissue by exposure to GDM may contribute to development of metabolic disease.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans

BACKGROUND

The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging.

RESULTS

We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks.

CONCLUSIONS

These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Passing down maternal dietary memories through lncRNAs

Parental diet is known to influence the offspring in an intergenerational manner, and this has been implicated in species adaptation and general health. Recent studies highlight the role of maternal long noncoding RNAs (lncRNAs) in serving as one of the 'memories' of maternal diet in regulating offspring development and predisposition to metabolic disease.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Batch-effect detection, correction and characterisation in Illumina HumanMethylation450 and MethylationEPIC BeadChip array data

Background

Genomic technologies can be subject to significant batch-effects which are known to reduce experimental power and to potentially create false positive results. The Illumina Infinium Methylation BeadChip is a popular technology choice for epigenome-wide association studies (EWAS), but presently, little is known about the nature of batch-effects on these designs. Given the subtlety of biological phenotypes in many EWAS, control for batch-effects should be a consideration.

Results

Using the batch-effect removal approaches in the ComBat and Harman software, we examined two in-house datasets and compared results with three large publicly available datasets, (1214 HumanMethylation450 and 1094 MethylationEPIC BeadChips in total), and find that despite various forms of preprocessing, some batch-effects persist. This residual batch-effect is associated with the day of processing, the individual glass slide and the position of the array on the slide. Consistently across all datasets, 4649 probes required high amounts of correction. To understand the impact of this set to EWAS studies, we explored the literature and found three instances where persistently batch-effect prone probes have been reported in abstracts as key sites of differential methylation. As well as batch-effect susceptible probes, we also discover a set of probes which are erroneously corrected. We provide batch-effect workflows for Infinium Methylation data and provide reference matrices of batch-effect prone and erroneously corrected features across the five datasets spanning regionally diverse populations and three commonly collected biosamples (blood, buccal and saliva).

Conclusions

Batch-effects are ever present, even in high-quality data, and a strategy to deal with them should be part of experimental design, particularly for EWAS. Batch-effect removal tools are useful to reduce technical variance in Infinium Methylation data, but they need to be applied with care and make use of post hoc diagnostic measures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01277-9.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia

Background

In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.

Data

We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.

Results

In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.

Conclusions

Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

Stage-differentiated ensemble modeling of DNA methylation landscapes uncovers salient biomarkers and prognostic signatures in colorectal cancer progression

Background

Aberrant DNA methylation acts epigenetically to skew the gene transcription rate up or down, contributing to cancer etiology. A gap in our understanding concerns the epigenomics of stagewise cancer progression. In this study, we have developed a comprehensive computational framework for the stage-differentiated modelling of DNA methylation landscapes in colorectal cancer (CRC).

Methods

The methylation β-matrix was derived from the public-domain TCGA data, converted into M-value matrix, annotated with AJCC stages, and analysed for stage-salient genes using an ensemble of approaches involving stage-differentiated modelling of methylation patterns and/or expression patterns. Differentially methylated genes (DMGs) were identified using a contrast against controls (adjusted p-value <0.001 and |log fold-change of M-value| >2), and then filtered using a series of all possible pairwise stage contrasts (p-value <0.05) to obtain stage-salient DMGs. These were then subjected to a consensus analysis, followed by matching with clinical data and performing Kaplan–Meier survival analysis to evaluate the impact of methylation patterns of consensus stage-salient biomarkers on disease prognosis.

Results

We found significant genome-wide changes in methylation patterns in cancer cases relative to controls agnostic of stage. The stage-differentiated models yielded the following consensus salient genes: one stage-I gene (FBN1), one stage-II gene (FOXG1), one stage-III gene (HCN1) and four stage-IV genes (NELL1, ZNF135, FAM123A, LAMA1). All the biomarkers were significantly hypermethylated in the promoter regions, indicating down-regulation of expression and implying a putative CpG island Methylator Phenotype (CIMP) manifestation. A prognostic signature consisting of FBN1 and FOXG1 survived all the analytical filters, and represents a novel early-stage epigenetic biomarker / target.

Conclusions

We have designed and executed a workflow for stage-differentiated epigenomic analysis of colorectal cancer progression, and identified several stage-salient diagnostic biomarkers, and an early-stage prognostic biomarker panel. The study has led to the discovery of an alternative CIMP-like signature in colorectal cancer, reinforcing the role of CIMP drivers in tumor pathophysiology.

Environmentally sensitive hotspots in the methylome of the early human embryo

In humans, DNA methylation marks inherited from gametes are largely erased following fertilisation, prior to construction of the embryonic methylome. Exploiting a natural experiment of seasonal variation including changes in diet and nutritional status in rural Gambia, we analysed three datasets covering two independent child cohorts and identified 259 CpGs showing consistent associations between season of conception (SoC) and DNA methylation. SoC effects were most apparent in early infancy, with evidence of attenuation by mid-childhood. SoC-associated CpGs were enriched for metastable epialleles, parent-of-origin-specific methylation and germline differentially methylated regions, supporting a periconceptional environmental influence. Many SoC-associated CpGs overlapped enhancers or sites of active transcription in H1 embryonic stem cells and fetal tissues. Half were influenced but not determined by measured genetic variants that were independent of SoC. Environmental 'hotspots' providing a record of environmental influence at periconception constitute a valuable resource for investigating epigenetic mechanisms linking early exposures to lifelong health and disease.

Association between parental unhealthy behaviors and offspring's cardiovascular health status: Results from a cross-sectional analysis of parent-offspring pairs in China

BACKGROUND

Cardiovascular health (CVH) in children and adolescents, which might be largely influenced by parental behaviors, may affect the incidence of cardiovascular diseases in adulthood. However, few studies have been conducted to explore the associations between parental behaviors and CVH status of offspring in China.

METHODS

Data were obtained from a cross-sectional survey conducted in Chinese children and adolescents aged 7-18 years old, with a total of 10,043 parent-offspring pairs included. Parental behaviors included moderate to vigorous physical activity (MVPA), dietary behaviors, and weight status. The CVH status of offspring was consulted by The American Heart Association, including seven factors. The associations between parental behaviors and CVH status of offspring were evaluated by multilevel logistic regression. Stratified analyses were conducted to explore the potential modifying influence of sociodemographic factors.

RESULTS

Most of the offspring had five ideal CVH factors; only 21.04% had six to seven ideal CVH factors. Parental unhealthy behaviors were associated with high odds of nonideal CVH status of offspring. Parental overweight/obesity, insufficient MVPA, and unhealthy dietary behaviors could increase the odds of owning one to three ideal CVH factors in offspring, with corresponding odds ratios (ORs) (95% confidence interval) of 1.61 (1.32-1.96), 1.31 (1.10-1.56), and 2.05 (1.43-2.94), respectively. There was a dose-response relationship between parental single unhealthy behavior and the odds of nonideal CVH status in offspring (P-trend < 0.001). Offspring with overweight parents had ORs of 1.25 for nonideal CVH status, compared to offspring with normal-weight parents. Among offspring who had the same number of ideal CVH factors, the cumulative association between unhealthy behaviors of parents and offspring's nonideal CVH status increased if parents had more unhealthy behaviors (P-trend < 0.001).

CONCLUSIONS

Parental overweight/obesity, insufficient MVPA, and unhealthy dietary behaviors were strongly associated with CVH status in offspring. With a cumulative association, more unhealthy parental behaviors were associated with higher odds of offspring's nonideal CVH status, suggesting that targeting parental behaviors might facilitate attainment of improving CVH status of children and adolescents.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome

BACKGROUND

parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist.

METHODS

in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315).

FINDINGS

a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs.

INTERPRETATION

the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases.

FUNDING

NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

A Review on Epigenetic Inheritance of Experiences in Humans

If genetics defines the inheritance of DNA, epigenetics aims to regulate and make it adaptable. Epigenetic alterations include DNA methylation, chromatin remodelling, post-translational modifications of histone proteins and activity of non-coding RNAs. Several studies, especially in animal models, have reported transgenerational inheritance of epigenetic marks. However, evidence of transgenerational inheritance in humans via germline in the absence of any direct exposure to the driving external stimulus remains controversial. Most of the epimutations exist in relation with genetic variants. The present review looks at intergenerational and transgenerational inheritance in humans, (both father and mother) in response to diet, exposure to chemicals, stress, exercise, and disease status. If not transgenerational, at least intergenerational human studies could help to understand early processes of inheritance. In humans, female and male germline development follow separate paths of epigenetic events and both oocyte and sperm possess their own unique epigenomes. While DNA methylation alterations are reset during epigenetic reprogramming, non-coding RNAs via human sperm provide evidence of being reliable carriers for transgenerational inheritance. Human studies reveal that one mechanism of epigenetic inheritance cannot be applied to the complete human genome. Multiple factors including time, type, and tissue of exposure determine if the modified epigenetic mark could be transmissible and till which generation. Population-specific differences should also be taken into consideration while associating inheritance to an environmental exposure. A longitudinal study targeting one environmental factor, but different population groups should be conducted at a specific geographical location to pinpoint heritable epigenetic changes.

Harnessing the power of multi-omics data for predicting climate change response

Predicting how species will respond to future climate change is of central importance in the midst of the global biodiversity crisis, and recent work has demonstrated the utility of population genomics for improving these predictions. Here, we suggest a broadening of the approach to include other types of genomic variants that play an important role in adaptation, like structural (e.g. copy number variants) and epigenetic variants (e.g. DNA methylation). These data could provide additional power for forecasting response, especially in weakly structured or panmictic species. Incorporating structural and epigenetic variation into estimates of climate change vulnerability, or maladaptation, may not only improve prediction power but also provide insight into the molecular mechanisms underpinning species' response to climate change.

Oocyte age and preconceptual alcohol use are highly correlated with epigenetic imprinting of a noncoding RNA (nc886)

Genomic imprinting occurs before fertilization, impacts every cell of the developing child, and may be sensitive to environmental perturbations. The noncoding RNA, nc886 (also called VTRNA2-1) is the only known example of the ∼100 human genes imprinted by DNA methylation, that shows polymorphic imprinting in the population. The nc886 gene is part of an ∼1.6-kb differentially methylated region (DMR) that is methylated in the oocyte and silenced on the maternal allele in about 75% of humans worldwide. Here, we show that the presence or absence of imprinting at the nc886 DMR in an individual is consistent across different tissues, confirming that the imprint is established before cellular differentiation and is maintained into adulthood. We investigated the relationships between the frequency of imprinting in newborns and maternal age, alcohol consumption and cigarette smoking before conception in more than 1,100 mother/child pairs from South Africa. The probability of imprinting in newborns was increased in older mothers and decreased in mothers who drank alcohol before conception. On the other hand, cigarette smoking had no apparent relationship with the frequency of imprinting. These data show an epigenetic change during oocyte maturation which is potentially subject to environmental influence. Much focus has been placed on avoiding alcohol consumption during pregnancy, but our data suggest that drinking before conception may affect the epigenome of the newborn.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Integrative DNA Methylation and Gene Expression Analysis of Cognitive Behavioral Therapy Response in Children and Adolescents with Obsessive-Compulsive Disorder; a Pilot Study

Purpose

Here, we propose an integrative analysis of genome-wide methylation and gene expression to provide new insight into the biological mechanisms of Cognitive behavioral therapy (CBT) in pediatric obsessive-compulsive disorder (OCD).

Patients and Methods

Twelve children and adolescents with OCD receiving CBT for the first time were classified as responders or non-responders after eight weeks of CBT. Differentially methylated positions (DMPs) and gene co-expression modules were identified using specific R software packages. Correlations between the DMPs and gene co-expression modules were investigated.

Results

Two genes were enriched with significant DMPs (Δβ > ± 0.2, FDR-adjusted p-value < 0.05): PIWIL1 and MIR886. The yellowgreen module of co-expressed genes was associated with CBT response (FDR-adjusted p-value = 0.0003). Significant correlations were observed between the yellowgreen module and the CpGs in PIWIL1 and MIR886 (p < 0.008). Patients showing hypermethylation in these CpGs presented an upregulation in the genes in the yellowgreen module.

Conclusion

Taken together, the preliminary results of this systems-level approach, despite the study limitations, provide evidence that the epigenetic regulation of ncRNAs could be a predictor of CBT response.

Limitations

The sample size limited the statistical power, and given that the study was hypothesis-driven, our results should be seen as preliminary.

Pan-cancer chromatin analysis of the human vtRNA genes uncovers their association with cancer biology

Background: The vault RNAs (vtRNAs) are a class of 84-141-nt eukaryotic non-coding RNAs transcribed by RNA polymerase III, associated to the ribonucleoprotein complex known as vault particle. Of the four human vtRNA genes, vtRNA1-1, vtRNA1-2 and vtRNA1-3, clustered at locus 1, are integral components of the vault particle, while vtRNA2-1 is a more divergent homologue located in a second locus. Gene expression studies of vtRNAs in large cohorts have been hindered by their unsuccessful sequencing using conventional transcriptomic approaches.

Methods: VtRNA expression in The Cancer Genome Atlas (TCGA) Pan-Cancer cohort was estimated using the genome-wide DNA methylation and chromatin accessibility data (ATAC-seq) of their genes as surrogate variables. The association between vtRNA expression and patient clinical outcome, immune subtypes and transcriptionally co-regulated gene programs was analyzed in the dataset.

Results: VtRNAs promoters are enriched in transcription factors related to viral infection. VtRNA2-1 is likely the most independently regulated homologue. VtRNA1-1 has the most accessible chromatin, followed by vtRNA1-2, vtRNA2-1 and vtRNA1-3. VtRNA1-1 and vtRNA1-3 chromatin status does not significantly change in cancer tissues. Meanwhile, vtRNA2-1 and vtRNA1-2 expression is widely deregulated in neoplastic tissues and its alteration is compatible with a broad oncogenic role for vtRNA1-2, and both tumor suppressor and oncogenic functions for vtRNA2-1. Yet, vtRNA1-1, vtRNA1-2 and vtRNA2-1 promoter DNA methylation predicts a shorter patient overall survival cancer-wide. In addition, gene ontology analyses of vtRNAs co-regulated genes identify a chromosome regulatory domain, epithelial differentiation, immune and thyroid cancer gene sets for specific vtRNAs. Furthermore, vtRNA expression patterns are associated with cancer immune subtypes and vtRNA1-2 expression is positively associated with cell proliferation and wound healing.

Conclusions: Our study presents the landscape of vtRNA chromatin status cancer-wide, identifying co-regulated gene networks and ontological pathways associated with the different vtRNA genes that may account for their diverse roles in cancer.

VTRNA2-1: Genetic Variation, Heritable Methylation and Disease Association

VTRNA2-1 is a metastable epiallele with accumulating evidence that methylation at this region is heritable, modifiable and associated with disease including risk and progression of cancer. This study investigated the influence of genetic variation and other factors such as age and adult lifestyle on blood DNA methylation in this region. We first sequenced the VTRNA2-1 gene region in multiple-case breast cancer families in which VTRNA2-1 methylation was identified as heritable and associated with breast cancer risk. Methylation quantitative trait loci (mQTL) were investigated using a prospective cohort study (4500 participants with genotyping and methylation data). The cis-mQTL analysis (334 variants ± 50 kb of the most heritable CpG site) identified 43 variants associated with VTRNA2-1 methylation (p < 1.5 × 10⁻⁴); however, these explained little of the methylation variation (R² < 0.5% for each of these variants). No genetic variants elsewhere in the genome were found to strongly influence VTRNA2-1 methylation. SNP-based heritability estimates were consistent with the mQTL findings (h² = 0, 95%CI: −0.14 to 0.14). We found no evidence that age, sex, country of birth, smoking, body mass index, alcohol consumption or diet influenced blood DNA methylation at VTRNA2-1. Genetic factors and adult lifestyle play a minimal role in explaining methylation variability at the heritable VTRNA2-1 cluster.

Methylation status of vault RNA 2-1 promoter is a predictor of glycemic response to glucagon-like peptide-1 analog therapy in type 2 diabetes mellitus

INTRODUCTION

Therapeutic efficiency of glucagon-like peptide-1 (GLP-1) analog is about 50%-70% in type 2 diabetes mellitus (T2DM). Discovery of potential genetic biomarkers for prediction of treatment efficiency of GLP-1 analog before therapy is still necessary. We assess whether DNA methylation was associated with glycemic response to GLP-1 analog therapy in patients with poorly controlled T2DM.

RESEARCH DESIGN AND METHODS

Genomic DNA was extracted from the peripheral blood of training (n=10) and validation (n=128) groups of patients with T2DM receiving GLP-1 analogs. DNA methylome was analyzed using Infinium Human Methylation EPIC Bead Chip in the training group. The candidate genes were examined using a pyrosequencing platform in the validation group. The association between DNA methylation status and glycemic response to GLP-1 was analyzed in these patients.

RESULTS

The most differential methylation region between those with a good (responsive) and poor (unresponsive) glycemic response to GLP-1 analog therapy was located on chromosome 5q31.1 (135415693 to 135416613), the promoter of VTRNA2-1 in the training group. The methylation status of the VTRNA2-1 promoter was examined in the validation group via pyrosequencing reaction, and the hypomethylation of VTRNA2-1 (<40% methylation) was significantly associated with poor glycemic response to GLP-1 treatment (OR 2.757, 95% CI 1.240 to 6.130, p=0.011). Since the VTRNA2-1 promoter region was previously reported maternal imprinting extended to the adjacent centromeric CCCTC-binding factor site that contained an A/C polymorphism (rs2346018), which was associated with methylation density of VTRNA2-1, this A/C polymorphism was also integrated to analyze association with glycemic response to GLP-1 analog therapy. In patients with the A allele of rs2346018 and hypomethylation (<40%) on the VTRNA2-1 promoter, the OR increased to 4.048 (95% CI 1.438 to 11.389, p=0.007).

CONCLUSIONS

The glycemic response to GLP-1 analog treatment is associated with the methylation status of the VTRNA2-1 promoter and polymorphism of rs2346018.

Parental cardiovascular health predicts time to onset of cardiovascular disease in offspring

BACKGROUND

Cardiovascular disease (CVD) risk factors are transmitted from parents to children. We prospectively examined the association between parental cardiovascular health (CVH) and time to onset of CVD in the offspring.

METHODS AND RESULTS

The study consisted of a total of 5967 offspring-mother-father trios derived from the Framingham Heart Study. Cardiovascular health score was defined using the seven American Heart Association's CVH metrics attained at ideal levels: poor (0-2), intermediate (3-4), and ideal CVH (5-7). Multivariable-adjusted Cox proportional hazards regression models, Kaplan-Meier plots, and Irwin's restricted mean were used to examine the association and sex-specific differences between parental CVH and offspring's CVD-free survival. In a total of 71 974 person-years of follow-up among the offspring, 718 incident CVD events occurred. The overall CVD incidence rate was 10 per 1000 person-years [95% confidence interval (CI) 9.3-10.7]. Offspring of mothers with ideal CVH lived 9 more years free of CVD than offspring of mothers with poor CVH (P < 0.001). Maternal poor CVH was associated with twice as high hazard of early onset of CVD compared with maternal ideal CVH (adjusted Hazard Ratio 2.09, 95% CI 1.50-2.92). No statistically significant association was observed in the hazards of CVD-free survival by paternal CVH categories.

CONCLUSIONS

We found that offspring of parents with ideal CVH had a greater CVD-free survival. Maternal CVH was a more robust predictor of offspring's CVD-free survival than paternal CVH, underscoring the need for clinical and policy interventions that involve mothers to break the intergenerational cycle of CVD-related morbidity and mortality.

Loss of Imprinting in Human Placentas Is Widespread, Coordinated, and Predicts Birth Phenotypes

Genomic imprinting leads to mono-allelic expression of genes based on parent of origin. Therian mammals and angiosperms evolved this mechanism in nutritive tissues, the placenta, and endosperm, where maternal and paternal genomes are in conflict with respect to resource allocation. We used RNA-seq to analyze allelic bias in the expression of 91 known imprinted genes in term human placentas from a prospective cohort study in Mali. A large fraction of the imprinted exons (39%) deviated from mono-allelic expression. Loss of imprinting (LOI) occurred in genes with either maternal or paternal expression bias, albeit more frequently in the former. We characterized LOI using binomial generalized linear mixed models. Variation in LOI was predominantly at the gene as opposed to the exon level, consistent with a single promoter driving the expression of most exons in a gene. Some genes were less prone to LOI than others, particularly lncRNA genes were rarely expressed from the repressed allele. Further, some individuals had more LOI than others and, within a person, the expression bias of maternally and paternally imprinted genes was correlated. We hypothesize that trans-acting maternal effect genes mediate correlated LOI and provide the mother with an additional lever to control fetal growth by extending her influence to LOI of the paternally imprinted genes. Limited evidence exists to support associations between LOI and offspring phenotypes. We show that birth length and placental weight were associated with allelic bias, making this the first comprehensive report of an association between LOI and a birth phenotype.

Protein kinase R and its cellular regulators in cancer: An active player or a surveillant?

Protein kinase R (PKR), originally known as an antiviral protein, senses various stresses as well as pathogen-driven double-stranded RNAs. Thereby activated PKR provokes diverse downstream events, including eIF2α phosphorylation and nuclear factor kappa-light-chain-enhancer of activated B cells activation. Consequently, PKR induces apoptosis and inflammation, both of which are highly important in cancer as much as its original antiviral role. Therefore, cellular proteins and RNAs should tightly control PKR activity. PKR and its regulators are often dysregulated in cancer and it is undoubted that such dysregulation contributes to tumorigenesis. However, PKR's precise role in cancer is still in debate, due to incomprehensible and even contradictory data. In this review, we introduce important cellular PKR regulators and discuss about their roles in cancer. Among them, we pay particular attention to nc886, a PKR repressor noncoding RNA that has been identified relatively recently, because its expression pattern in cancer can explain interesting yet obscure oncologic aspects of PKR. Based on nc886 and its regulation of PKR, we have proposed a tumor surveillance model, which reconciles contradictory data about PKR in cancer. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

A genomic atlas of systemic interindividual epigenetic variation in humans

BACKGROUND

DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific.

RESULTS

For each of 10 donors from the NIH Genotype-Tissue Expression (GTEx) program, CpG methylation is measured by deep whole-genome bisulfite sequencing of genomic DNA from tissues representing the three germ layer lineages: thyroid (endoderm), heart (mesoderm), and brain (ectoderm). We develop a computational algorithm to identify genomic regions at which interindividual variation in DNA methylation is consistent across all three lineages. This approach identifies 9926 correlated regions of systemic interindividual variation (CoRSIVs). These regions, comprising just 0.1% of the human genome, are inter-correlated over long genomic distances, associated with transposable elements and subtelomeric regions, conserved across diverse human ethnic groups, sensitive to periconceptional environment, and associated with genes implicated in a broad range of human disorders and phenotypes. CoRSIV methylation in one tissue can predict expression of associated genes in other tissues.

CONCLUSIONS

In addition to charting a previously unexplored molecular level of human individuality, this atlas of human CoRSIVs provides a resource for future population-based investigations into how interindividual epigenetic variation modulates risk of disease.