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

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.

Systemic interindividual DNA methylation variants in cattle share major hallmarks with those in humans

BACKGROUND

We recently identified ~ 10,000 correlated regions of systemic interindividual epigenetic variation (CoRSIVs) in the human genome. These methylation variants are amenable to population studies, as DNA methylation measurements in blood provide information on epigenetic regulation throughout the body. Moreover, establishment of DNA methylation at human CoRSIVs is labile to periconceptional influences such as nutrition. Here, we analyze publicly available whole-genome bisulfite sequencing data on multiple tissues of each of two Holstein cows to determine whether CoRSIVs exist in cattle.

RESULTS

Focusing on genomic blocks with ≥ 5 CpGs and a systemic interindividual variation index of at least 20, our approach identifies 217 cattle CoRSIVs, a subset of which we independently validate by bisulfite pyrosequencing. Similar to human CoRSIVs, those in cattle are strongly associated with genetic variation. Also as in humans, we show that establishment of DNA methylation at cattle CoRSIVs is particularly sensitive to early embryonic environment, in the context of embryo culture during assisted reproduction.

CONCLUSIONS

Our data indicate that CoRSIVs exist in cattle, as in humans, suggesting these systemic epigenetic variants may be common to mammals in general. To the extent that individual epigenetic variation at cattle CoRSIVs affects phenotypic outcomes, assessment of CoRSIV methylation at birth may become an important tool for optimizing agriculturally important traits. Moreover, adjusting embryo culture conditions during assisted reproduction may provide opportunities to tailor agricultural outcomes by engineering CoRSIV methylation profiles.

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.

Crosstalk between vault RNAs and innate immunity

PURPOSE

Vault (vt) RNAs are noncoding (nc) RNAs transcribed by RNA polymerase III (RNA Pol III) with 5'-triphosphate (5'-PPP) termini that play significant roles and are recognized by innate immune sensors, including retinoic acid-inducible protein 1 (RIG-I). In addition, vtRNAs adopt secondary structures that can be targets of interferon-inducible protein kinase R (PKR) and the oligoadenylate synthetase (OAS)/RNase L system, both of which are important for activating antiviral defenses. However, changes in the expression of vtRNAs have been associated with pathological processes that activate proinflammatory pathways, which influence cellular events such as differentiation, aging, autophagy, apoptosis, and drug resistance in cancer cells.

RESULTS

In this review, we summarized the biology of vtRNAs and focused on their interactions with the innate immune system. These findings provide insights into the diverse roles of vtRNAs and their correlation with various cellular processes to improve our understanding of their biological functions.

A POLR3B-variant reveals a Pol III transcriptome response dependent on La protein/SSB

RNA polymerase III (Pol III, POLR3) synthesizes tRNAs and other small non-coding RNAs. Human POLR3 pathogenic variants cause a range of developmental disorders, recapitulated in part by mouse models, yet some aspects of POLR3 deficiency have not been explored. We characterized a human POLR3B:c.1625A>G;p.(Asn542Ser) disease variant that was found to cause mis-splicing of POLR3B. Genome-edited POLR3B1625A>G HEK293 cells acquired the mis-splicing with decreases in multiple POLR3 subunits and TFIIIB, although display auto-upregulation of the Pol III termination-reinitiation subunit POLR3E. La protein was increased relative to its abundant pre-tRNA ligands which bind via their U(n)U-3'-termini. Assays for cellular transcription revealed greater deficiencies for tRNA genes bearing terminators comprised of 4Ts than of ≥5Ts. La-knockdown decreased Pol III ncRNA expression unlinked to RNA stability. Consistent with these effects, small-RNAseq showed that POLR3B1625A>G and patient fibroblasts express more tRNA fragments (tRFs) derived from pre-tRNA 3'-trailers (tRF-1) than from mature-tRFs, and higher levels of multiple miRNAs, relative to control cells. The data indicate that decreased levels of Pol III transcripts can lead to functional excess of La protein which reshapes small ncRNA profiles revealing new depth in the Pol III system. Finally, patient cell RNA analysis uncovered a strategy for tRF-1/tRF-3 as POLR3-deficiency biomarkers.

New insights into the (epi)genetics of twinning

Spontaneous dizygotic (DZ) twins, i.e. twins conceived without the use of ARTs, run in families and their prevalence varies widely around the globe. In contrast, monozygotic (MZ) twins occur at a constant rate across time and geographical regions and, with some rare exceptions, do not cluster in families. The leading hypothesis for MZ twins, which arise when a zygote splits during preimplantation stages of development, is random occurrence. We have found the first series of genes underlying the liability of being the mother of DZ twins and have shown that being an MZ twin is strongly associated with a stable DNA methylation signature in child and adult somatic tissues. Because identical twins keep this molecular signature across the lifespan, this discovery opens up completely new possibilities for the retrospective diagnosis of whether a person is an MZ twin whose co-twin may have vanished in the early stages of pregnancy. Here, we summarize the gene finding results for mothers of DZ twins based on genetic association studies followed by meta-analysis, and further present the striking epigenetic results for MZ twins.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

Analysis of DNA methylation at birth and in childhood reveals changes associated with season of birth and latitude

BACKGROUND

Seasonal variations in environmental exposures at birth or during gestation are associated with numerous adult traits and health outcomes later in life. Whether DNA methylation (DNAm) plays a role in the molecular mechanisms underlying the associations between birth season and lifelong phenotypes remains unclear.

METHODS

We carried out epigenome-wide meta-analyses within the Pregnancy And Childhood Epigenetic Consortium to identify associations of DNAm with birth season, both at differentially methylated probes (DMPs) and regions (DMRs). Associations were examined at two time points: at birth (21 cohorts, N = 9358) and in children aged 1-11 years (12 cohorts, N = 3610). We conducted meta-analyses to assess the impact of latitude on birth season-specific associations at both time points.

RESULTS

We identified associations between birth season and DNAm (False Discovery Rate-adjusted p values < 0.05) at two CpGs at birth (winter-born) and four in the childhood (summer-born) analyses when compared to children born in autumn. Furthermore, we identified twenty-six differentially methylated regions (DMR) at birth (winter-born: 8, spring-born: 15, summer-born: 3) and thirty-two in childhood (winter-born: 12, spring and summer: 10 each) meta-analyses with few overlapping DMRs between the birth seasons or the two time points. The DMRs were associated with genes of known functions in tumorigenesis, psychiatric/neurological disorders, inflammation, or immunity, amongst others. Latitude-stratified meta-analyses [higher (≥ 50°N), lower (< 50°N, northern hemisphere only)] revealed differences in associations between birth season and DNAm by birth latitude. DMR analysis implicated genes with previously reported links to schizophrenia (LAX1), skin disorders (PSORS1C, LTB4R), and airway inflammation including asthma (LTB4R), present only at birth in the higher latitudes (≥ 50°N).

CONCLUSIONS

In this large epigenome-wide meta-analysis study, we provide evidence for (i) associations between DNAm and season of birth that are unique for the seasons of the year (temporal effect) and (ii) latitude-dependent variations in the seasonal associations (spatial effect). DNAm could play a role in the molecular mechanisms underlying the effect of birth season on adult health outcomes.

Child Survival and Annual Crop Yield Reductions in Rural Burkina Faso: Critical Windows of Vulnerability Around Early-Life Development

Populations that are reliant on subsistence farming are particularly vulnerable to climatic effects on crop yields. However, empirical evidence on the role of the timing of exposure to crop yield deficits in early-life development is limited. We examined the relationship between child survival and annual crop yield reductions at different stages of early-life development in a subsistence farming population in Burkina Faso. Using shared frailty Cox proportional hazards models adjusting for confounders, we analyzed 57,288 children under 5 years of age followed by the Nouna Health and Demographic Surveillance System (1994-2016) in relation to provincial food-crop yield levels experienced in 5 nonoverlapping time windows (12 months before conception, gestation, birth-age 5.9 months, ages 6.0 months-1.9 years, and ages 2.0-4.9 years) and their aggregates (birth-1.9 years, first 1,000 days from conception, and birth-4.9 years). Of the nonoverlapping windows, point estimates were largest for child survival related to food-crop yields for the time window of 6.0 months-1.9 years: The adjusted mortality hazard ratio was 1.10 (95% confidence interval: 1.03, 1.19) for a 90th-to-10th percentile yield reduction. These findings suggest that child survival in this setting is particularly vulnerable to cereal-crop yield reductions during the period of nonexclusive breastfeeding.

The influence of early environment and micronutrient availability on developmental epigenetic programming: lessons from the placenta

DNA methylation is the most commonly studied epigenetic mark in humans, as it is well recognised as a stable, heritable mark that can affect genome function and influence gene expression. Somatic DNA methylation patterns that can persist throughout life are established shortly after fertilisation when the majority of epigenetic marks, including DNA methylation, are erased from the pre-implantation embryo. Therefore, the period around conception is potentially critical for influencing DNA methylation, including methylation at imprinted alleles and metastable epialleles (MEs), loci where methylation varies between individuals but is correlated across tissues. Exposures before and during conception can affect pregnancy outcomes and health throughout life. Retrospective studies of the survivors of famines, such as those exposed to the Dutch Hunger Winter of 1944-45, have linked exposures around conception to later disease outcomes, some of which correlate with DNA methylation changes at certain genes. Animal models have shown more directly that DNA methylation can be affected by dietary supplements that act as cofactors in one-carbon metabolism, and in humans, methylation at birth has been associated with peri-conceptional micronutrient supplementation. However, directly showing a role of micronutrients in shaping the epigenome has proven difficult. Recently, the placenta, a tissue with a unique hypomethylated methylome, has been shown to possess great inter-individual variability, which we highlight as a promising target tissue for studying MEs and mixed environmental exposures. The placenta has a critical role shaping the health of the fetus. Placenta-associated pregnancy complications, such as preeclampsia and intrauterine growth restriction, are all associated with aberrant patterns of DNA methylation and expression which are only now being linked to disease risk later in life.

Small noncoding vault RNA2-1 disrupts gut epithelial barrier function via interaction with HuR

Vault RNAs (vtRNAs) are small noncoding RNAs and highly expressed in many eukaryotes. Here, we identified vtRNA2-1 as a novel regulator of the intestinal barrier via interaction with RNA-binding protein HuR. Intestinal mucosal tissues from patients with inflammatory bowel diseases and from mice with colitis or sepsis express increased levels of vtRNAs relative to controls. Ectopically expressed vtRNA2-1 decreases the levels of intercellular junction (IJ) proteins claudin 1, occludin, and E-cadherin and causes intestinal epithelial barrier dysfunction in vitro, whereas vtRNA2-1 silencing promotes barrier function. Increased vtRNA2-1 also decreases IJs in intestinal organoid, inhibits epithelial renewal, and causes Paneth cell defects ex vivo. Elevating the levels of tissue vtRNA2-1 in the intestinal mucosa increases the vulnerability of the gut barrier to septic stress in mice. vtRNA2-1 interacts with HuR and prevents HuR binding to claudin 1 and occludin mRNAs, thus decreasing their translation. These results indicate that vtRNA2-1 impairs intestinal barrier function by repressing HuR-facilitated translation of claudin 1 and occludin.

Epigenetics and the Exposome: DNA Methylation as a Proxy for Health Impacts of Prenatal Environmental Exposures

The accumulation of every day exposures can impact health across the life course, but our understanding of such exposures is impeded by our ability to delineate the relationship between an individual's early life exposome and later life health effects. Measuring the exposome is challenging. Exposure assessed at a given time point captures a snapshot of the exposome but does not represent the full spectrum of exposures across the life course. In addition, the assessment of early life exposures and their effects is often further challenged by lack of relevant samples and the time gap between exposures and related health outcomes in later life. Epigenetics, specifically DNA methylation, has the potential to overcome these barriers as environmental epigenetic perturbances can be retained through time. In this review, we describe how DNA methylation can be framed in the world of the exposome. We offer three compelling examples of common environmental exposures, including cigarette smoke, the endocrine active compound bisphenol A (BPA), and the metal lead (Pb), to illustrate the application of DNA methylation as a proxy to measure the exposome. We discuss areas for future explorations and current limitations of this approach. Epigenetic profiling is a promising and rapidly developing tool and field of study, offering us a unique and powerful way to assess the early life exposome and its effects across different life stages.

Systemic interindividual epigenetic variation in humans is associated with transposable elements and under strong genetic control

BACKGROUND

Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.

RESULTS

We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R²=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.

CONCLUSIONS

A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

Increased copy number of imprinted genes in the chromosomal region 20q11-q13.32 is associated with resistance to antitumor agents in cancer cell lines

BACKGROUND

Parent of origin-specific allelic expression of imprinted genes is epigenetically controlled. In cancer, imprinted genes undergo both genomic and epigenomic alterations, including frequent copy number changes. We investigated whether copy number loss or gain of imprinted genes in cancer cell lines is associated with response to chemotherapy treatment.

RESULTS

We analyzed 198 human imprinted genes including protein-coding genes and noncoding RNA genes using data from tumor cell lines from the Cancer Cell Line Encyclopedia and Genomics of Drug Sensitivity in Cancer datasets. We examined whether copy number of the imprinted genes in 35 different genome locations was associated with response to cancer drug treatment. We also analyzed associations of pretreatment expression and DNA methylation of imprinted genes with drug response. Higher copy number of BLCAP, GNAS, NNAT, GNAS-AS1, HM13, MIR296, MIR298, and PSIMCT-1 in the chromosomal region 20q11-q13.32 was associated with resistance to multiple antitumor agents. Increased expression of BLCAP and HM13 was also associated with drug resistance, whereas higher methylation of gene regions of BLCAP, NNAT, SGK2, and GNAS was associated with drug sensitivity. While expression and methylation of imprinted genes in several other chromosomal regions was also associated with drug response and many imprinted genes in different chromosomal locations showed a considerable copy number variation, only imprinted genes at 20q11-q13.32 had a consistent association of their copy number with drug response. Copy number values among the imprinted genes in the 20q11-q13.32 region were strongly correlated. They were also correlated with the copy number of cancer-related non-imprinted genes MYBL2, AURKA, and ZNF217 in that chromosomal region. Expression of genes at 20q11-q13.32 was associated with ex vivo drug response in primary tumor samples from the Beat AML 1.0 acute myeloid leukemia patient cohort. Association of the increased copy number of the 20q11-q13.32 region with drug resistance may be complex and could involve multiple genes.

CONCLUSIONS

Copy number of imprinted and non-imprinted genes in the chromosomal region 20q11-q13.32 was associated with cancer drug resistance. The genes in this chromosomal region may have a modulating effect on tumor response to chemotherapy.

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.

Tissue- and ethnicity-independent hypervariable DNA methylation states show evidence of establishment in the early human embryo

We analysed DNA methylation data from 30 datasets comprising 3474 individuals, 19 tissues and 8 ethnicities at CpGs covered by the Illumina450K array. We identified 4143 hypervariable CpGs ('hvCpGs') with methylation in the top 5% most variable sites across multiple tissues and ethnicities. hvCpG methylation was influenced but not determined by genetic variation, and was not linked to probe reliability, epigenetic drift, age, sex or cell heterogeneity effects. hvCpG methylation tended to covary across tissues derived from different germ-layers and hvCpGs were enriched for proximity to ERV1 and ERVK retrovirus elements. hvCpGs were also enriched for loci previously associated with periconceptional environment, parent-of-origin-specific methylation, and distinctive methylation signatures in monozygotic twins. Together, these properties position hvCpGs as strong candidates for studying how stochastic and/or environmentally influenced DNA methylation states which are established in the early embryo and maintained stably thereafter can influence life-long health and disease.

Roles of transposable elements in the regulation of mammalian transcription

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

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.

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.

Pregnancy exposure to synthetic phenols and placental DNA methylation - An epigenome-wide association study in male infants from the EDEN cohort

In utero exposure to environmental chemicals, such as synthetic phenols, may alter DNA methylation in different tissues, including placenta - a critical organ for fetal development. We studied associations between prenatal urinary biomarker concentrations of synthetic phenols and placental DNA methylation. Our study involved 202 mother-son pairs from the French EDEN cohort. Nine phenols were measured in spot urine samples collected between 22 and 29 gestational weeks. We performed DNA methylation analysis of the fetal side of placental tissues using the IlluminaHM450 BeadChips. We evaluated methylation changes of individual CpGs in an adjusted epigenome-wide association study (EWAS) and identified differentially methylated regions (DMRs). We performed mediation analysis to test whether placental tissue heterogeneity mediated the association between urinary phenol concentrations and DNA methylation. We identified 46 significant DMRs (≥5 CpGs) associated with triclosan (37 DMRs), 2,4-dichlorophenol (3), benzophenone-3 (3), methyl- (2) and propylparaben (1). All but 2 DMRs were positively associated with phenol concentrations. Out of the 46 identified DMRs, 7 (6 for triclosan) encompassed imprinted genes (APC, FOXG1, GNAS, GNASAS, MIR886, PEG10, SGCE), which represented a significant enrichment. Other identified DMRs encompassed genes encoding proteins responsible for cell signaling, transmembrane transport, cell adhesion, inflammatory, apoptotic and immunological response, genes encoding transcription factors, histones, tumor suppressors, genes involved in tumorigenesis and several cancer risk biomarkers. Mediation analysis suggested that placental cell heterogeneity may partly explain these associations. This is the first study describing the genome-wide modifications of placental DNA methylation associated with pregnancy exposure to synthetic phenols or their precursors. Our results suggest that cell heterogeneity might mediate the effects of triclosan exposure on placental DNA methylation. Additionally, the enrichment of imprinted genes within the DMRs suggests mechanisms by which certain exposures, mainly to triclosan, could affect fetal development.

DNA methylation at a nutritionally sensitive region of the PAX8 gene is associated with thyroid volume and function in Gambian children

PAX8 is a key thyroid transcription factor implicated in thyroid gland differentiation and function, and PAX8 gene methylation is reported to be sensitive to the periconceptional environment. Using a novel recall-by-epigenotype study in Gambian children, we found that PAX8 hypomethylation at age 2 years is associated with a 21% increase in thyroid volume and an increase in free thyroxine (T4) at 5 to 8 years, the latter equivalent to 8.4% of the normal range. Free T4 was associated with a decrease in DXA-derived body fat and bone mineral density. Furthermore, offspring PAX8 methylation was associated with periconceptional maternal nutrition, and methylation variability was influenced by genotype, suggesting that sensitivity to environmental exposures may be under partial genetic control. Together, our results demonstrate a possible link between early environment, PAX8 gene methylation and thyroid gland development and function, with potential implications for early embryonic programming of thyroid-related health and disease.

Environmental Exposures around Conception: Developmental Pathways Leading to Lifetime Disease Risk

Environment around conception can influence the developmental programme with lasting effects on gestational and postnatal phenotype and with consequences for adult health and disease risk. Peri-conception exposure comprises a crucial part of the 'Developmental Origins of Health and Disease' (DOHaD) concept. In this review, we consider the effects of maternal undernutrition experienced during the peri-conception period in select human models and in a mouse experimental model of protein restriction. Human datasets indicate that macronutrient deprivation around conception affect the epigenome, with enduring effects on cardiometabolic and neurological health. The mouse model, comprising maternal low protein diet exclusively during the peri-conception period, has revealed a stepwise progression in altered developmental programming following induction through maternal metabolite deficiency. This progression includes differential effects in extra-embryonic and embryonic cell lineages and tissues, leading to maladaptation in the growth trajectory and increased chronic disease comorbidities. The timeline embraces an array of mechanisms across nutrient sensing and signalling, cellular, metabolic, epigenetic and physiological processes with a coordinating role for mTORC1 signalling proposed. Early embryos appear active participants in environmental sensing to optimise the developmental programme for survival but with the trade-off of later disease. Similar adverse health outcomes may derive from other peri-conception environmental experiences, including maternal overnutrition, micronutrient availability, pollutant exposure and assisted reproductive treatments (ART) and support the need for preconception health before pregnancy.

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.

A machine learning case-control classifier for schizophrenia based on DNA methylation in blood

Epigenetic dysregulation is thought to contribute to the etiology of schizophrenia (SZ), but the cell type-specificity of DNA methylation makes population-based epigenetic studies of SZ challenging. To train an SZ case-control classifier based on DNA methylation in blood, therefore, we focused on human genomic regions of systemic interindividual epigenetic variation (CoRSIVs), a subset of which are represented on the Illumina Human Methylation 450K (HM450) array. HM450 DNA methylation data on whole blood of 414 SZ cases and 433 non-psychiatric controls were used as training data for a classification algorithm with built-in feature selection, sparse partial least squares discriminate analysis (SPLS-DA); application of SPLS-DA to HM450 data has not been previously reported. Using the first two SPLS-DA dimensions we calculated a "risk distance" to identify individuals with the highest probability of SZ. The model was then evaluated on an independent HM450 data set on 353 SZ cases and 322 non-psychiatric controls. Our CoRSIV-based model classified 303 individuals as cases with a positive predictive value (PPV) of 80%, far surpassing the performance of a model based on polygenic risk score (PRS). Importantly, risk distance (based on CoRSIV methylation) was not associated with medication use, arguing against reverse causality. Risk distance and PRS were positively correlated (Pearson r = 0.28, P = 1.28 × 10^-12), and mediational analysis suggested that genetic effects on SZ are partially mediated by altered methylation at CoRSIVs. Our results indicate two innate dimensions of SZ risk: one based on genetic, and the other on systemic epigenetic variants.

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.

Elevated methylation of the vault RNA2-1 promoter in maternal blood is associated with preterm birth

Background

Preterm birth, defined as parturition before 37 completed weeks of gestation, is associated with an increased risk of neonatal complications and death, as well as poor health and disease later in life. Epigenetics could contribute to the mechanism underlying preterm birth.

Results

Genome-wide DNA methylation analysis of whole blood cells from 10 women (5 term and 5 preterm deliveries) was performed using an Illumina Infinium HumanMethylation450 BeadChips array. We identified 1,581 differentially methylated CpG sites in promoter regions between term and preterm birth. Although the differences were not significant after correcting for multiple tests, seven CpGs on the genomically imprinted vault RNA2-1 (VTRNA2-1; also known as non-coding RNA, nc886 or miR-886) showed the largest differences (range: 26–39 %). Pyrosequencing verification was performed with blood samples from pregnant women recruited additionally (39 term and 43 preterm deliveries). In total, 28 (34.1 %) samples showed hypomethylation of the VTRNA2-1 promoter (< 13 % methylation), while 54 (65.9 %) samples showed elevated methylation levels between 30 and 60 %. Elevated methylation of VTRNA2-1 promoter was associated with an increased risk of preterm birth after adjusting for maternal age, season of delivery, parity and white blood cell count. The mRNA expression of VTRNA2-1 was 0.51-fold lower in women with preterm deliveries (n = 20) compared with women with term deliveries (n = 20).

Conclusions

VTRNA2-1 is a noncoding transcript to environmentally responsive epialleles. Our results suggest that elevated methylation of the VTRNA2-1 promoter may result in increased risk of PTB caused by the pro-inflammatory cytokines. Further studies are needed to confirm the association of VTRNA2-1 methylation with preterm birth in a large population, and to elucidate the underlying mechanism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07865-y.

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.

Maternal Microbiota, Early Life Colonization and Breast Milk Drive Immune Development in the Newborn

The innate immune system is the oldest protection strategy that is conserved across all organisms. Although having an unspecific action, it is the first and fastest defense mechanism against pathogens. Development of predominantly the adaptive immune system takes place after birth. However, some key components of the innate immune system evolve during the prenatal period of life, which endows the newborn with the ability to mount an immune response against pathogenic invaders directly after birth. Undoubtedly, the crosstalk between maternal immune cells, antibodies, dietary antigens, and microbial metabolites originating from the maternal microbiota are the key players in preparing the neonate’s immunity to the outer world. Birth represents the biggest substantial environmental change in life, where the newborn leaves the protective amniotic sac and is exposed for the first time to a countless variety of microbes. Colonization of all body surfaces commences, including skin, lung, and gastrointestinal tract, leading to the establishment of the commensal microbiota and the maturation of the newborn immune system, and hence lifelong health. Pregnancy, birth, and the consumption of breast milk shape the immune development in coordination with maternal and newborn microbiota. Discrepancies in these fine-tuned microbiota interactions during each developmental stage can have long-term effects on disease susceptibility, such as metabolic syndrome, childhood asthma, or autoimmune type 1 diabetes. In this review, we will give an overview of the recent studies by discussing the multifaceted emergence of the newborn innate immune development in line with the importance of maternal and early life microbiota exposure and breast milk intake.

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.

Genomic properties of variably methylated retrotransposons in mouse

BACKGROUND

Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function.

RESULTS

Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation.

CONCLUSION

Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

Nc886, a Novel Suppressor of the Type I Interferon Response Upon Pathogen Intrusion

Interferons (IFNs) are a crucial component in the innate immune response. Especially the IFN-β signaling operates in most cell types and plays a key role in the first line of defense upon pathogen intrusion. The induction of IFN-β should be tightly controlled, because its hyperactivation can lead to tissue damage or autoimmune diseases. Activation of the IFN-β promoter needs Interferon Regulatory Factor 3 (IRF3), together with Nuclear Factor kappa-light-chain-enhancer of activated B cells (NF-κB) and Activator Protein 1 (AP-1). Here we report that a human noncoding RNA, nc886, is a novel suppressor for the IFN-β signaling and inflammation. Upon treatment with several pathogen-associated molecular patterns and viruses, nc886 suppresses the activation of IRF3 and also inhibits NF-κB and AP-1 via inhibiting Protein Kinase R (PKR). These events lead to decreased expression of IFN-β and resultantly IFN-stimulated genes. nc886′s role might be to restrict the IFN-β signaling from hyperactivation. Since nc886 expression is regulated by epigenetic and environmental factors, nc886 might explain why innate immune responses to pathogens are variable depending on biological settings.

The genome-wide impact of trisomy 21 on DNA methylation and its implications for hematopoiesis

Down syndrome is associated with genome-wide perturbation of gene expression, which may be mediated by epigenetic changes. We perform an epigenome-wide association study on neonatal bloodspots comparing 196 newborns with Down syndrome and 439 newborns without Down syndrome, adjusting for cell-type heterogeneity, which identifies 652 epigenome-wide significant CpGs (P < 7.67 × 10-8) and 1,052 differentially methylated regions. Differential methylation at promoter/enhancer regions correlates with gene expression changes in Down syndrome versus non-Down syndrome fetal liver hematopoietic stem/progenitor cells (P < 0.0001). The top two differentially methylated regions overlap RUNX1 and FLI1, both important regulators of megakaryopoiesis and hematopoietic development, with significant hypermethylation at promoter regions of these two genes. Excluding Down syndrome newborns harboring preleukemic GATA1 mutations (N = 30), identified by targeted sequencing, has minimal impact on the epigenome-wide association study results. Down syndrome has profound, genome-wide effects on DNA methylation in hematopoietic cells in early life, which may contribute to the high frequency of hematological problems, including leukemia, in children with Down syndrome.

Identification of methylation changes associated with positive and negative growth deviance in Gambian infants using a targeted methyl sequencing approach of genomic DNA

Low birthweight and reduced height gain during infancy (stunting) may arise at least in part from adverse early life environments that trigger epigenetic reprogramming that may favor survival. We examined differential DNA methylation patterns using targeted methyl sequencing of regions regulating gene activity in groups of rural Gambian infants: (a) low and high birthweight (DNA from cord blood (n = 16 and n = 20, respectively), from placental trophoblast tissue (n = 21 and n = 20, respectively), and DNA from peripheral blood collected from infants at 12 months of age (n = 23 and n = 17, respectively)), and, (b) the top 10% showing rapid postnatal length gain (high, n = 20) and the bottom 10% showing slow postnatal length gain (low, n = 20) based on z score change between birth and 12 months of age (LAZ) (DNA from peripheral blood collected from infants at 12 months of age). Using BiSeq analysis to identify significant methylation marks, for birthweight, four differentially methylated regions (DMRs) were identified in trophoblast DNA, compared to 68 DMRs in cord blood DNA, and 54 DMRs in 12‐month peripheral blood DNA. Twenty‐five DMRs were observed to be associated with high and low length for age (LAZ) at 12 months. With the exception of five loci (associated with two different genes), there was no overlap between these groups of methylation marks. Of the 194 CpG methylation marks contained within DMRs, 106 were located to defined gene regulatory elements (promoters, CTCF‐binding sites, transcription factor‐binding sites, and enhancers), 58 to gene bodies (introns or exons), and 30 to intergenic DNA. Distinct methylation patterns associated with birthweight between comparison groups were observed in DNA collected at birth (at the end of intrauterine growth window) compared to those established by 12 months (near the infancy/childhood growth transition). The longitudinal differences in methylation patterns may arise from methylation adjustments, changes in cellular composition of blood or both that continue during the critical postnatal growth period, and in response to early nutritional and infectious environmental exposures with impacts on growth and longer‐term health outcomes.

DNA Methylation Profiling and Genomic Analysis in 20 Children with Short Stature Who Were Born Small for Gestational Age

PURPOSE

In a significant proportion of children born small for gestational age (SGA) with failure of catch-up growth, the etiology of short stature remains unclear after routine diagnostic workup. We wanted to investigate if extensive analysis of the (epi)genome can unravel the cause of growth failure in a significant portion of these children.

PATIENTS AND METHODS

Twenty SGA children treated with GH because of short stature were selected from the BELGROW database of the Belgian Society for Pediatric Endocrinology and Diabetology for exome sequencing, single-nucleotide polymorphism (SNP) array and genome-wide methylation analysis to identify the (epi)genetic cause. First-year response to GH was compared with the response of SGA patients in the KIGS database.

RESULTS

We identified (likely) pathogenic variants in 4 children (from 3 families) using exome sequencing and found pathogenic copy number variants in 2 probands using SNP array. In a child harboring a NSD1-containing microduplication, we identified a DNA methylation signature that is opposite to the genome-wide DNA methylation signature of Sotos syndrome. Moreover, we observed multilocus imprinting disturbances in 2 children in whom no other genomic alteration could be identified. Five of 6 children with a genetic diagnosis had an "above average" response to GH.

CONCLUSIONS

The study indicates that a more advanced approach with deep genotyping can unravel unexpected (epi)genomic alterations in SGA children with persistent growth failure. Most SGA children with a genetic diagnosis had a good response to GH treatment.

Effect of maternal preconceptional and pregnancy micronutrient interventions on children's DNA methylation: Findings from the EMPHASIS study

BACKGROUND

Maternal nutrition in pregnancy has been linked to offspring health in early and later life, with changes to DNA methylation (DNAm) proposed as a mediating mechanism.

OBJECTIVE

We investigated intervention-associated DNAm changes in children whose mothers participated in 2 randomized controlled trials of micronutrient supplementation before and during pregnancy, as part of the EMPHASIS (Epigenetic Mechanisms linking Preconceptional nutrition and Health Assessed in India and sub-Saharan Africa) study (ISRCTN14266771).

DESIGN

We conducted epigenome-wide association studies with blood samples from Indian (n = 698) and Gambian (n = 293) children using the Illumina EPIC array and a targeted study of selected loci not on the array. The Indian micronutrient intervention was food based, whereas the Gambian intervention was a micronutrient tablet.

RESULTS

We identified 6 differentially methylated CpGs in Gambians [2.5-5.0% reduction in intervention group, all false discovery rate (FDR) <5%], the majority mapping to ESM1, which also represented a strong signal in regional analysis. One CpG passed FDR <5% in the Indian cohort, but overall effect sizes were small (<1%) and did not have the characteristics of a robust signature. We also found strong evidence for enrichment of metastable epialleles among subthreshold signals in the Gambian analysis. This supports the notion that multiple methylation loci are influenced by micronutrient supplementation in the early embryo.

CONCLUSIONS

Maternal preconceptional and pregnancy micronutrient supplementation may alter DNAm in children measured at 7-9 y. Multiple factors, including differences between the nature of the intervention, participants, and settings, are likely to have contributed to the lack of replication in the Indian cohort. Potential links to phenotypic outcomes will be explored in the next stage of the EMPHASIS study.

Can Children Catch up from the Consequences of Undernourishment? Evidence from Child Linear Growth, Developmental Epigenetics, and Brain and Neurocognitive Development

Recovery from nutritionally induced height deficits continues to garner attention. The current literature on catch-up growth, however, has 2 important limitations: wide-ranging definitions of catch-up growth are used, and it remains unclear whether children can recover from the broader consequences of undernutrition. We addressed these shortcomings by reviewing the literature on the criteria for catch-up in linear growth and on the potential to recover from undernutrition early in life in 3 domains: linear growth, developmental epigenetics, and child brain and neurocognitive development. Four criteria must be met to demonstrate catch-up growth in height: after a period in which a growth-inhibiting condition (criterion 1) causes a reduction in linear growth velocity (criterion 2), alleviation of the inhibiting condition (criterion 3) leads to higher-than-normal velocity (criterion 4). Accordingly, studies that are observational, do not use absolute height, or have no alleviation of an inhibiting condition cannot be used to establish catch-up growth. Adoption and foster care, which provide dramatic improvements in children's living conditions not typically attained in nutrition interventions, led to some (but incomplete) recovery in linear growth and brain and neurocognitive development. Maternal nutrition around the time of conception was shown to have long-term (potentially permanent) effects on DNA methylation in the offspring. Undernourishment early in life may thus have profound irreversible effects. Scientific, program, and policy efforts should focus on preventing maternal and child undernutrition rather than on correcting its consequences or attempting to prove they can be corrected.

Periconceptional environment predicts leukocyte telomere length in a cross-sectional study of 7-9 year old rural Gambian children

Early life exposures are important predictors of adult disease risk. Although the underlying mechanisms are largely unknown, telomere maintenance may be involved. This study investigated the relationship between seasonal differences in parental exposures at time of conception and leukocyte telomere length (LTL) in their offspring. LTL was measured in two cohorts of children aged 2 yrs (N = 487) and 7–9 yrs (N = 218). The association between date of conception and LTL was examined using Fourier regression models, adjusted for age, sex, leukocyte cell composition, and other potential confounders. We observed an effect of season in the older children in all models [likelihood ratio test (LRT) χ²₂ = 7.1, p = 0.03; fully adjusted model]. LTL was greatest in children conceived in September (in the rainy season), and smallest in those conceived in March (in the dry season), with an effect size (LTL peak–nadir) of 0.60 z-scores. No effect of season was evident in the younger children (LRT χ²₂ = 0.87, p = 0.65). The different results obtained for the two cohorts may reflect a delayed effect of season of conception on postnatal telomere maintenance. Alternatively, they may be explained by unmeasured differences in early life exposures, or the increased telomere attrition rate during infancy.

DNA methylation at the crossroads of gene and environment interactions

DNA methylation is an epigenetic mark involved in regulating genome function and is critical for normal development in mammals. It has been observed that the developmental environment can lead to permanent changes in gene expression and DNA methylation, at least at 'metastable epialleles'. These are defined as regions of the genome that show a variable epigenetic state that is established early in development and maintained through subsequent cell divisions. However, the majority of the known genome does not behave in this manner. Here, we use the developmental origins of adult disease hypothesis to understand environmental epigenomics. Some challenges to studying how DNA methylation is influenced by the environment include identifying DNA methylation changes associated with an environmental exposure in tissues with a complex cellular composition and at genomic regions for which DNA methylation is dynamically regulated in a cell-type specific manner. We also offer a perspective of how emerging technologies may be useful for dissecting the functional contribution of exposure-associated epigenetic changes and highlight recent evidence that suggests that genomic regions that are absent from genome assemblies may be unappreciated hotspots for environmental modulation of the epigenetic state.

Noncoding RNA 886 alleviates tumor cellular immunological rejection in host C57BL/C mice

Non‐coding RNA 886 (nc886/VTRNA2‐1) is a Pol III transcript and an atypical imprinted gene. Its exact function as a negative regulator of protein kinase R establishes its connection with innate immunity. Studies have shown that nc886 silencing is closely associated with prostate cancer progression. Previous work has constructed a cell model of stable nc886 overexpression (“mimic” or “nc886⁺”) in PC‐3M‐1E8 cell lines (1E8), which are highly bone‐metastatic human prostate cancer cells with low expression of nc886, and cells expressing the mimic were validated to have lower invasive and metastatic abilities than cells expressing the scramble transcript in vitro and in vivo. In this study, we directly injected mimic or scramble cells into the left ventricle of C57BL/C mice, an immunocompetent animal model, to elucidate the immune mechanisms of tumor‐host interactions. Interestingly, we found that tumor cells induced the inflammation of many important organs due to xenogeneic antigen rejection; this inflammation was ultimately repaired by tissue fibrosis after 28 days, except for in the spleen. The reason is that mimic cells, as heterogeneous antigens, are mostly directly recognized by macrophages or T cells in blood, and few mimic cells enter the spleen compared with scramble cells. The induction of splenic macrophage polarization to M2 macrophages by scramble cells is a critical factor in maintaining chronic splenic inflammation. In addition, we recognize that nc886 broadly decreases the expression of some human leukocyte antigen molecules and antigen transporters. This evidence reveals the interesting role of nc886 in regulating tumor cell antigens.

Differential hypermethylation of the VTRNA2-1 promoter in hepatocellular carcinoma as a prognostic factor: Tumor marker prognostic study

BACKGROUND

Vault RNA 2-1 (VTRNA2-1, also called nc886) is a 108-nucleotide noncoding transcript that is epigenetically controlled via 18 CpG dinucleotide modifications of its promoter, and can exert either tumor suppressor or oncogenic functions depending on cell types of cancers. In hepatocellular carcinoma (HCC), the role of VTRNA2-1 in prognosis of patients remains unexplored. Here, we analysed the methylation status of the VTRNA2-1 promoter and its correlation with clinical parameters in patients with HCC.

PATIENTS AND METHODS

A total of 92 patients with HCC were enrolled, genomic DNA of tumor versus normal tissues were extracted and bisulfite modified. VTRNA2-1 promoter regions chr5: 135416381 (cg06536614), 135416388, 135416394 (cg26328633), and 135416398 (cg25340688) were PCR amplified and pyrosequenced. The methylation status of VTRNA2-1 in patients was further analysed with other clinical parameters via univariate and multivariate analysis.

RESULTS

The differential hypermethylation status (tumor- normal) of the VTRNA2-1 promoter in HCC correlated well with the presence of large tumor size (p = 0.001), pathological vascular invasion (p = 0.036), tumor recurrence (p = 0.007) and more advanced tumor stage (stage III AJCC) in patients (p = 0.03). In addition, the methylation of the VTRNA2-1 promoter increased in stage III HCC tumor compared with stage I & II tumor (64.7% versus 36.0%, p = 0.030). Furthermore, the differential hypermethylation status of the VTRNA2-1 promoter was an independent factor for patient outcome after partial hepatectomy using multivariate Cox regression analysis (p = 0.011, HR = 2.305). Using another public dataset (GSE89852), we found that the differential hypermethylation of the VTRNA2-1 promoter was also significantly associated with tumor recurrence.

CONCLUSIONS

Patients had unfavourable outcomes when the VTRNA2-1 promoter was differentially hypermethylated in tumor tissues compared to its adjacent normal tissues. These findings suggest that such patients should receive intensive follow-up care or possible adjuvant therapy after liver resection.

Making headway towards understanding how epigenetic mechanisms contribute to early-life effects

It has become clear that in addition to the DNA sequence there is another layer of information, termed epigenetic modifications, that can influence phenotypes and traits. In particular, environmental epigenomics, which addresses the effect of the environment on the epigenome and human health, is becoming an area of great interest for many researchers working in different scientific fields. In this review, we will consider the current evidence that early-life environmental signals can have long-term effects on the epigenome. We will further evaluate how recent technological advances may enable us to unravel the molecular mechanisms underlying these phenomena, which will be crucial for understanding heritability in health and disease. This article is part of the theme issue 'Developing differences: early-life effects and evolutionary medicine'.

Genetically Predicted Levels of DNA Methylation Biomarkers and Breast Cancer Risk: Data From 228 951 Women of European Descent

BACKGROUND

DNA methylation plays a critical role in breast cancer development. Previous studies have identified DNA methylation marks in white blood cells as promising biomarkers for breast cancer. However, these studies were limited by low statistical power and potential biases. Using a new methodology, we investigated DNA methylation marks for their associations with breast cancer risk.

METHODS

Statistical models were built to predict levels of DNA methylation marks using genetic data and DNA methylation data from HumanMethylation450 BeadChip from the Framingham Heart Study (n = 1595). The prediction models were validated using data from the Women's Health Initiative (n = 883). We applied these models to genomewide association study (GWAS) data of 122 977 breast cancer patients and 105 974 controls to evaluate if the genetically predicted DNA methylation levels at CpG sites (CpGs) are associated with breast cancer risk. All statistical tests were two-sided.

RESULTS

Of the 62 938 CpG sites CpGs investigated, statistically significant associations with breast cancer risk were observed for 450 CpGs at a Bonferroni-corrected threshold of P less than 7.94 × 10-7, including 45 CpGs residing in 18 genomic regions, that have not previously been associated with breast cancer risk. Of the remaining 405 CpGs located within 500 kilobase flaking regions of 70 GWAS-identified breast cancer risk variants, the associations for 11 CpGs were independent of GWAS-identified variants. Integrative analyses of genetic, DNA methylation, and gene expression data found that 38 CpGs may affect breast cancer risk through regulating expression of 21 genes.

CONCLUSION

Our new methodology can identify novel DNA methylation biomarkers for breast cancer risk and can be applied to other diseases.

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.

Epigenetic biomarkers and preterm birth

Preterm birth (PTB) is a major public health challenge, and novel, sensitive approaches to predict PTB are still evolving. Epigenomic markers are being explored as biomarkers of PTB because of their molecular stability compared to gene expression. This approach is also relatively new compared to gene-based diagnostics, which relies on mutations or single nucleotide polymorphisms. The fundamental principle of epigenome diagnostics is that epigenetic reprogramming in the target tissue (e.g. placental tissue) might be captured by more accessible surrogate tissue (e.g. blood) using biochemical epigenome assays on circulating DNA that incorporate methylation, histone modifications, nucleosome positioning, and/or chromatin accessibility. Epigenomic-based biomarkers may hold great potential for early identification of the majority of PTBs that are not associated with genetic variants or mutations. In this review, we discuss recent advances made in the development of epigenome assays focusing on its potential exploration for association and prediction of PTB. We also summarize population-level cohort studies conducted in the USA and globally that provide opportunities for genetic and epigenetic marker development for PTB. In addition, we summarize publicly available epigenome resources and published PTB studies. We particularly focus on ongoing genome-wide DNA methylation and epigenome-wide association studies. Finally, we review the limitations of current research, the importance of establishing a comprehensive biobank, and possible directions for future studies in identifying effective epigenome biomarkers to enhance health outcomes for pregnant women at risk of PTB and their infants.

Metastable epialleles and their contribution to epigenetic inheritance in mammals

Many epigenetic differences between individuals are driven by genetic variation. Mammalian metastable epialleles are unusual in that they show variable DNA methylation states between genetically identical individuals. The occurrence of such states across generations has resulted in their consideration by many as strong evidence for epigenetic inheritance in mammals, with the classic A^vy and Axin^Fu mouse models - each products of repeat element insertions - being the most widely accepted examples. Equally, there has been interest in exploring their use as epigenetic biosensors given their susceptibility to environmental compromise. Here we review the classic murine metastable epialleles as well as more recently identified candidates, with the aim of providing a more holistic understanding of their biology. We consider the extent to which epigenetic inheritance occurs at metastable epialleles and explore the limited mechanistic insights into the establishment of their variable epigenetic states. We discuss their environmental modulation and their potential relevance in genome regulation. In light of recent whole-genome screens for novel metastable epialleles, we point out the need to reassess their biological relevance in multi-generational studies and we highlight their value as a model to study repeat element silencing as well as the mechanisms and consequences of mammalian epigenetic stochasticity.

The role of environmental exposures and the epigenome in health and disease

The genetic material of every organism exists within the context of regulatory networks that govern gene expression, collectively called the epigenome. Epigenetics has taken center stage in the study of diseases such as cancer and diabetes, but its integration into the field of environmental health is still emerging. As the Environmental Mutagenesis and Genomics Society (EMGS) celebrates its 50th Anniversary this year, we have come together to review and summarize the seminal advances in the field of environmental epigenomics. Specifically, we focus on the role epigenetics may play in multigenerational and transgenerational transmission of environmentally induced health effects. We also summarize state of the art techniques for evaluating the epigenome, environmental epigenetic analysis, and the emerging field of epigenome editing. Finally, we evaluate transposon epigenetics as they relate to environmental exposures and explore the role of noncoding RNA as biomarkers of environmental exposures. Although the field has advanced over the past several decades, including being recognized by EMGS with its own Special Interest Group, recently renamed Epigenomics, we are excited about the opportunities for environmental epigenetic science in the next 50 years. Environ. Mol. Mutagen. 61:176-192, 2020. © 2019 Wiley Periodicals, Inc.