Mutant Nrf2E79Q enhances the promotion and progression of a subset of oncogenic Ras keratinocytes and skin tumors.

Squamous cell carcinomas (SCCs), including lung, head & neck, bladder, and skin SCCs often display constitutive activation of the KEAP1-NRF2 pathway. Constitutive activation is achieved through multiple mechanisms, including activating mutations in NFE2L2 (NRF2). To determine the functional consequences of Nrf2 activation on skin SCC development, we assessed the effects of mutant Nrf2E79Q expression, one of the most common activating mutations in human SCCs, on tumor promotion and progression in the mouse skin multistage carcinogenesis model using a DMBA-initiation/TPA-promotion protocol where the Hras A->T mutation (Q61L) is the canonical driver mutation. Nrf2E79Q expression was temporally and conditionally activated in the epidermis at two stages of tumor development: 1) after DMBA initiation in the epidermis but before cutaneous tumor development and 2) in pre-existing DMBA-initiated/TPA-promoted squamous papillomas. Expression of Nrf2E79Q in the epidermis after DMBA initiation but before tumor occurrence inhibited the development/promotion of 70% of squamous papillomas. However, the remaining papillomas often displayed non-canonical Hras and Kras mutations and enhanced progression to SCCs compared to control mice expressing wildtype Nrf2. Nrf2E79Q expression in pre-existing tumors caused rapid regression of 60% of papillomas. The remaining papillomas displayed the expected canonical Hras A->T mutation (Q61L) and enhanced progression to SCCs. These results demonstrate that mutant Nrf2E79Q enhances the promotion and progression of a subset of skin tumors and alters the frequency and diversity of oncogenic Ras mutations when expressed early after initiation.

DNA methylation of imprint control regions associated with Alzheimer's disease in non-Hispanic Blacks and non-Hispanic Whites.

Alzheimer's disease (AD) prevalence is twice as high in non-Hispanic Blacks (NHBs) as in non-Hispanic Whites (NHWs). The objective of this study was to determine whether aberrant methylation at imprint control regions (ICRs) is associated with AD. Differentially methylated regions (DMRs) were bioinformatically identified from whole-genome bisulfite sequenced DNA derived from brain tissue of 9 AD (5 NHBs and 4 NHWs) and 8 controls (4 NHBs and 4 NHWs). We identified DMRs located within 120 regions defined as candidate ICRs in the human imprintome ( https://genome.ucsc.edu/s/imprintome/hg38.AD.Brain_track ). Eighty-one ICRs were differentially methylated in NHB-AD, and 27 ICRs were differentially methylated in NHW-AD, with two regions common to both populations that are proximal to the inflammasome gene, NLRP1, and a known imprinted gene, MEST/MESTIT1. These findings indicate that early developmental alterations in DNA methylation of regions regulating genomic imprinting may contribute to AD risk and that this epigenetic risk differs between NHBs and NHWs.

Developmental regulation of cellular metabolism is required for intestinal elongation and rotation.

Malrotation of the intestine is a prevalent birth anomaly, the etiology of which remains poorly understood. Here, we show that late-stage exposure of Xenopus embryos to atrazine, a widely used herbicide that targets electron transport chain (ETC) reactions, elicits intestinal malrotation at high frequency. Interestingly, atrazine specifically inhibits the cellular morphogenetic events required for gut tube elongation, including cell rearrangement, differentiation and proliferation; insufficient gut lengthening consequently reorients the direction of intestine rotation. Transcriptome analyses of atrazine-exposed intestines reveal misexpression of genes associated with glycolysis and oxidative stress, and metabolomics shows that atrazine depletes key glycolytic and tricarboxylic acid cycle metabolites. Moreover, cellular bioenergetics assays indicate that atrazine blocks a crucial developmental transition from glycolytic ATP production toward oxidative phosphorylation. Atrazine-induced defects are phenocopied by rotenone, a known ETC Complex I inhibitor, accompanied by elevated reactive oxygen species, and rescued by antioxidant supplementation, suggesting that malrotation may be at least partly attributable to redox imbalance. These studies reveal roles for metabolism in gut morphogenesis and implicate defective gut tube elongation and/or metabolic perturbations in the etiology of intestinal malrotation.

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

Background: Differentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile these regions, but it is very expensive and it requires not only a high coverage but it is also computationally intensive to assess those regions. Findings: To address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 9,757 probes map to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites and obtained a mean R 2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array and obtained a mean R 2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799-0.945). Conclusions: Our custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.

Maternal educational attainment in pregnancy and epigenome-wide DNA methylation changes in the offspring from birth until adolescence.

Maternal educational attainment (MEA) shapes offspring health through multiple potential pathways. Differential DNA methylation may provide a mechanistic understanding of these long-term associations. We aimed to quantify the associations of MEA with offspring DNA methylation levels at birth, in childhood and in adolescence. Using 37 studies from high-income countries, we performed meta-analysis of epigenome-wide association studies (EWAS) to quantify the associations of completed years of MEA at the time of pregnancy with offspring DNA methylation levels at birth (n = 9 881), in childhood (n = 2 017), and adolescence (n = 2 740), adjusting for relevant covariates. MEA was found to be associated with DNA methylation at 473 cytosine-phosphate-guanine sites at birth, one in childhood, and four in adolescence. We observed enrichment for findings from previous EWAS on maternal folate, vitamin-B12 concentrations, maternal smoking, and pre-pregnancy BMI. The associations were directionally consistent with MEA being inversely associated with behaviours including smoking and BMI. Our findings form a bridge between socio-economic factors and biology and highlight potential pathways underlying effects of maternal education. The results broaden our understanding of bio-social associations linked to differential DNA methylation in multiple early stages of life. The data generated also offers an important resource to help a more precise understanding of the social determinants of health.

An epigenome-wide association study of child appetitive traits and DNA methylation.

The etiology of childhood appetitive traits is poorly understood. Early-life epigenetic processes may be involved in the developmental programming of appetite regulation in childhood. One such process is DNA methylation (DNAm), whereby a methyl group is added to a specific part of DNA, where a cytosine base is next to a guanine base, a CpG site. We meta-analyzed epigenome-wide association studies (EWASs) of cord blood DNAm and early-childhood appetitive traits. Data were from two independent cohorts: the Generation R Study (n = 1,086, Rotterdam, the Netherlands) and the Healthy Start study (n = 236, Colorado, USA). DNAm at autosomal methylation sites in cord blood was measured using the Illumina Infinium HumanMethylation450 BeadChip. Parents reported on their child's food responsiveness, emotional undereating, satiety responsiveness and food fussiness using the Children's Eating Behaviour Questionnaire at age 4-5 years. Multiple regression models were used to examine the association of DNAm (predictor) at the individual site- and regional-level (using DMRff) with each appetitive trait (outcome), adjusting for covariates. Bonferroni-correction was applied to adjust for multiple testing. There were no associations of DNAm and any appetitive trait when examining individual CpG-sites. However, when examining multiple CpGs jointly in so-called differentially methylated regions, we identified 45 associations of DNAm with food responsiveness, 7 associations of DNAm with emotional undereating, 13 associations of DNAm with satiety responsiveness, and 9 associations of DNAm with food fussiness. This study shows that DNAm in the newborn may partially explain variation in appetitive traits expressed in early childhood and provides preliminary support for early programming of child appetitive traits through DNAm. Investigating differential DNAm associated with appetitive traits could be an important first step in identifying biological pathways underlying the development of these behaviors.

AHRR Hypomethylation mediates the association between maternal smoking and metabolic profiles in children.

BACKGROUND: Tobacco smoking during pregnancy is associated with metabolic dysfunction in children, but mechanistic insights remain limited. Hypomethylation of cg05575921 in the aryl hydrocarbon receptor repressor (AHRR) gene is associated with in utero tobacco smoke exposure. In this study, we evaluated whether AHRR hypomethylation mediates the association between maternal smoking and metabolic dysfunction in children. METHODS: We assessed metabolic dysfunction using liver fat content (LFC), serum, and clinical data in children aged 7-12 years (n=78) followed since birth. Maternal smoking was self-reported at 12 weeks gestation. Methylation was measured by means of pyrosequencing at 3 sequential CpG sites, including cg05575921, at birth and at ages 7-12. Regression models were used to evaluate whether AHRR methylation mediated the association between maternal smoking and child metabolic dysfunction. RESULTS: Average AHRR methylation at birth was significantly higher among children of nonsmoking mothers compared with children of mothers who smoked (69.8% ± 4.4% vs. 63.5% ± 5.5, p=0.0006). AHRR hypomethylation at birth was associated with higher liver fat content (p=0.01), triglycerides (p=0.01), and alanine aminotransferase levels (p=0.03), and lower HDL cholesterol (p=0.01) in childhood. AHRR hypomethylation significantly mediated associations between maternal smoking and liver fat content (indirect effect=0.213, p=0.018), triglycerides (indirect effect=0.297, p=0.044), and HDL cholesterol (indirect effect = -0.413, p=0.007). AHRR methylation in childhood (n=78) was no longer significantly associated with prenatal smoke exposure or child metabolic parameters (p>0.05). CONCLUSIONS: AHRR hypomethylation significantly mediates the association between prenatal tobacco smoke exposure and features of childhood metabolic dysfunction, despite the lack of persistent hypomethylation of AHRR into childhood. Further studies are needed to replicate these findings and to explore their causal and long-term significance.

Machine learning reveals common transcriptomic signatures across rat brain and placenta following developmental organophosphate ester exposure.

Toxicogenomics is a critical area of inquiry for hazard identification and to identify both mechanisms of action and potential markers of exposure to toxic compounds. However, data generated by these experiments are highly dimensional and present challenges to standard statistical approaches, requiring strict correction for multiple comparisons. This stringency often fails to detect meaningful changes to low expression genes and/or eliminate genes with small but consistent changes particularly in tissues where slight changes in expression can have important functional differences, such as brain. Machine learning offers an alternative analytical approach for "omics" data that effectively sidesteps the challenges of analyzing highly dimensional data. Using 3 rat RNA transcriptome sets, we utilized an ensemble machine learning approach to predict developmental exposure to a mixture of organophosphate esters (OPEs) in brain (newborn cortex and day 10 hippocampus) and late gestation placenta of male and female rats, and identified genes that informed predictor performance. OPE exposure had sex specific effects on hippocampal transcriptome, and significantly impacted genes associated with mitochondrial transcriptional regulation and cation transport in females, including voltage-gated potassium and calcium channels and subunits. To establish if this holds for other tissues, RNAseq data from cortex and placenta, both previously published and analyzed via a more traditional pipeline, were reanalyzed with the ensemble machine learning methodology. Significant enrichment for pathways of oxidative phosphorylation and electron transport chain was found, suggesting a transcriptomic signature of OPE exposure impacting mitochondrial metabolism across tissue types and developmental epoch. Here we show how machine learning can complement more traditional analytical approaches to identify vulnerable "signature" pathways disrupted by chemical exposures and biomarkers of exposure.

An epigenome-wide association study of child appetitive traits and DNA methylation.

Childhood appetitive traits are consistently associated with obesity risk, and yet their etiology is poorly understood. Appetitive traits are complex phenotypes which are hypothesized to be influenced by both genetic and environmental factors, as well as their interactions. Early-life epigenetic processes, such as DNA methylation (DNAm), may be involved in the developmental programming of appetite regulation in childhood. In the current study, we meta-analyzed epigenome-wide association studies (EWASs) of cord blood DNAm and early-childhood appetitive traits. Data were from two independent cohorts: the Generation R Study (n=1,086, Rotterdam, the Netherlands) and the Healthy Start study (n=236, Colorado, USA). DNAm at autosomal methylation sites in cord blood was measured using the Illumina Infinium HumanMethylation450 BeadChip. Parents reported on their child's food responsiveness, emotional undereating, satiety responsiveness and food fussiness using the Children's Eating Behaviour Questionnaire at age 4-5 years. Multiple regression models were used to examine the association of DNAm (predictor) at the individual site- and regional-level (using DMRff) with each appetitive trait (outcome), adjusting for covariates. Bonferroni-correction was applied to adjust for multiple testing. There were no associations of DNAm and any appetitive trait at the individual site-level. However, at the regional level, we identified 45 associations of DNAm with food responsiveness, 7 associations of DNAm with emotional undereating, 13 associations of DNAm with satiety responsiveness, and 9 associations of DNAm with food fussiness. This study shows that DNAm in the newborn may partially explain variation in appetitive traits expressed in early childhood and provides preliminary support for early programming of child appetitive traits through DNAm. Investigating differential DNAm associated with appetitive traits could be an important first step in identifying biological pathways underlying the development of these behaviors.

LncRNA Tuna is activated in cadmium-induced placental insufficiency and drives the NRF2-mediated oxidative stress response.

Cadmium (Cd) is a toxic heavy metal found throughout the environment and one of the top ten toxicants of major public health concern identified by the World Health Organization. In utero Cd exposure causes fetal growth restriction, malformation, and spontaneous abortion; however, the mechanisms by which Cd impacts these outcomes are poorly understood. Cd accumulates in the placenta, suggesting that these negative outcomes may be a consequence of disrupted placental function and placental insufficiency. To understand the impact of Cd on gene expression within the placenta, we developed a mouse model of Cd-induced fetal growth restriction through maternal consumption of CdCl2 and performed RNA-seq on control and CdCl2 exposed placentae. The top differentially expressed transcript was the Tcl1 Upstream Neuron-Associated (Tuna) long non-coding RNA, which was upregulated over 25-fold in CdCl2 exposed placentae. Tuna has been shown to be critical for neural stem cell differentiation. However, within the placenta, there is no evidence that Tuna is normally expressed or functional at any developmental stage. To determine the spatial expression of Cd-activated Tuna within the placenta, we used in situ hybridization as well as placental layer-specific RNA isolation and analysis. Both methods confirmed the absence of Tuna expression in control samples and determined that Cd-induced Tuna expression is specific to the junctional zone. Since many lncRNAs regulate gene expression, we hypothesized that Tuna forms part of the mechanism of Cd-induced transcriptomic changes. To test this, we over-expressed Tuna in cultured choriocarcinoma cells and compared gene expression profiles to those of control and CdCl2 exposed cells. We demonstrate significant overlap between genes activated by Tuna overexpression and genes activated by CdCl2 exposure, with enrichment in the NRF2-mediated oxidative stress response. Herein we analyze the NRF2 pathway and show that Tuna increases NRF2/NRF2 both at the transcript and protein levels. Tuna drives increased NRF2 target gene expression, a result that is abrogated with the use of an NRF2 inhibitor, confirming that Tuna activates oxidative stress response genes through this pathway. This work identifies the lncRNA Tuna as a potential novel player in Cd-induced placental insufficiency.

C/EBPß deficiency enhances the keratinocyte innate immune response to direct activators of cytosolic pattern recognition receptors.

The skin is the first line of defense to cutaneous microbes and viruses, and epidermal keratinocytes play a critical role in preventing infection by viruses and pathogens through activation of the type I interferon (IFN) response. Using RNAseq analysis, here we report that the conditional deletion of C/EBPß transcription factor in mouse epidermis (CKOß mice) resulted in the upregulation of IFNß and numerous keratinocyte interferon-stimulated genes (ISGs). The expression of cytosolic pattern recognition receptors (cPRRs), that recognize viral RNA and DNA, were significantly increased, and enriched in the RNAseq data set. cPRRs stimulate a type I IFN response that can trigger cell death to eliminate infected cells. To determine if the observed increases in cPRRs had functional consequences, we transfected CKOß primary keratinocytes with the pathogen and viral mimics poly(I:C) (dsRNA) or poly(dA:dT) (synthetic B-DNA) that directly activate PRRs. Transfected CKOß primary keratinocytes displayed an amplified type I IFN response which was accompanied by increased activation of IRF3, enhanced ISG expression, enhanced activation of caspase-8, caspase-3 and increased apoptosis. Our results identify C/EBPß as a critical repressor of the keratinocyte type I IFN response, and demonstrates that the loss of C/EBPß primes keratinocytes to the activation of cytosolic PRRs by pathogen RNA and DNA to induce cell death mediated by caspase-8 and caspase-3.

Salinity-induced ionoregulatory changes in the gill proteome of the mayfly, Neocloeon triangulifer.

Ecologists have observed declines in the biodiversity of sensitive freshwater organisms in response to increasing concentrations of major ions (salinization). Yet, how changing salinities physiologically challenge aquatic organisms, such as mayflies, remains remarkably understudied. Moreover, it is not well understood the degree to which species respond and acclimate to salinity changes. Our lab is developing the Baetid mayfly, N. triangulifer, as a model organism for physiological research. We have previously described acclimatory changes in both ion flux rates and altered mRNA transcript levels in response to chronic exposures to elevated major ion concentrations at the whole animal level. In the present study, we use shotgun proteomics to identify the specific proteins associated with apical ion transport and how their abundance changes in response to chronic salinity exposures in gills. Gills were isolated from the penultimate nymphal stage of N. triangulifer reared under control culture conditions, elevated NaCl (157 mg L-1 Na), elevated CaCl2 (121 mg L-1 Ca), elevated Ca/MgSO4 (735 mg L-1 SO4). These conditions mirrored those from previously published physiological work. We also acutely exposed nymphs to dilute (50% dilution of culture water with deionized water) to explore proteomic changes in the gills in response to dilute conditions. We report 710 unique peptide sequences among treatment groups, including important apical ion transporters such as Ca-ATPase, Na/K ATPase, and V-ATPase. Treatment with elevated NaCl and Ca/MgSO4 appeared to cause more significant differential protein expression (452 and 345, respectively) compared to CaCl2 and dilute groups (134 and 17, respectively). Finally, we demonstrated the breadth of physiological functions in gills by exploring non-transport related pathways found in our dataset, including ATP synthesis, calcium signaling, and oxidative stress response. We discuss our results in the context of freshwater salinization and the challenges of working with non-model species without fully sequenced and annotated genomes.

TCDD alters essential transcriptional regulators of osteogenic differentiation in multipotent mesenchymal stem cells.

Differentiation of multipotent mesenchymal stem cells (MSCs) into bone-forming osteoblasts requires strict coordination of transcriptional pathways. Aryl hydrocarbon receptor ligands, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), have been shown to alter osteoblast differentiation in vitro and bone formation in multiple developmental in vivo models. The goal of the present study was to establish a global transcriptomic landscape during early, intermediate, and apical stages of osteogenic differentiation in vitro in response to TCDD exposure. Human bone-derived mesenchymal stem cells (hBMSCs) were cultured in growth media (GM), osteogenic differentiation media (ODM), or ODM containing 10 nM TCDD (ODM + TCDD), thus enabling a comparison of the transcriptomic profiles of undifferentiated, differentiated, and differentiated-TCDD-exposed hBMSCs, respectively. In this test system, exposure to TCDD attenuated the differentiation of hBMSCs into osteoblasts as evidenced by reduced alkaline phosphatase activity and mineralization. At various timepoints, we observed altered expression of genes that play a role in the Wnt, fibroblast growth factor, bone morphogenetic protein/transforming growth factor beta developmental pathways, as well as pathways related to extracellular matrix organization and deposition. Reconstruction of gene regulatory networks with the interactive dynamic regulatory event miner (iDREM) analysis revealed modulation of transcription factors (TFs) including POLR3G, NR4A1, RDBP, GTF2B, POU2F2, and ZEB1, which may putatively influence osteoblast differentiation and the requisite deposition and mineralization of bone extracellular matrix. We demonstrate that the combination of RNA-Seq data in conjunction with the iDREM regulatory model captures the transcriptional dynamics underlying MSC differentiation under different conditions in vitro. Model predictions are consistent with existing knowledge and provide a new tool to identify novel pathways and TFs that may facilitate a better understanding of the osteoblast differentiation process, perturbation by exogenous agents, and potential intervention strategies targeting those specific pathways.

The imprinted gene Zac1 regulates steatosis in developmental cadmium-induced nonalcoholic fatty liver disease.

Cadmium (Cd) exposure in adulthood is associated with nonalcoholic fatty liver disease (NAFLD), characterized by steatosis, inflammation, and fibrosis. The prevalence of NAFLD in children is increasing, suggesting a role for the developmental environment in programming susceptibility. However, the role of developmental Cd exposure in programming NAFLD and the underlying mechanisms remain unclear. We have proposed that imprinted genes are strong candidates for connecting the early life environment and later life disease. In support of this, we previously identified roles for the Imprinted Gene Network (IGN) and its regulator Zac1 in programming NAFLD in response to maternal metabolic dysfunction. Here, we test the hypothesis that developmental Cd exposure is sufficient to program NAFLD, and further, that this process is mediated by Zac1 and the IGN. Using mice, we show that developmental cadmium chloride (CdCl2) exposure leads to histological, biochemical, and molecular signatures of steatosis and fibrosis in juveniles. Transcriptomic analyses comparing livers of CdCl2-exposed and control mice show upregulation of Zac1 and the IGN coincident with disease presentation. Increased hepatic Zac1 expression is independent of promoter methylation and imprinting statuses. Finally, we show that over-expression of Zac1 in cultured hepatocytes is sufficient to induce lipid accumulation in a Pparγ-dependent manner and demonstrate direct binding of Zac1 to the Pparγ promoter. Our findings demonstrate that developmental Cd exposure is sufficient to program NAFLD in later life, and with our previous work, establish Zac1 and the IGN as key regulators of prosteatotic and profibrotic pathways, two of the major pathological hallmarks of NAFLD.

Genome-wide transcriptome analysis of the orphan crop tef (Eragrostis tef (Zucc.) Trotter) under long-term low calcium stress.

Calcium (Ca2+) is one of the essential mineral nutrients for plant growth and development. However, the effects of long-term Ca2+ deficiency in orphan crops such as Tef [(Eragrostis tef) (Zucc.) Trotter], which accumulate high levels of Ca in the grains, remained unknown. Tef is a staple crop for nearly 70 million people in East Africa, particularly in Ethiopia and Eritrea. It is one of the most nutrient-dense grains, and is also more resistant to marginal soils and climatic conditions than main cereals like corn, wheat, and rice. In this study, tef plants were grown in a hydroponic solution containing optimum (1 mM) or low (0.01 mM) Ca2+, and plant growth parameters and whole-genome transcriptome were analyzed. Ca+2-deficient plants exhibited leaf necrosis, leaf curling, and growth stunting symptoms. Ca2+ deficiency significantly decreased root and shoot Ca, potassium (K), and copper content in both root and shoots. At the same time, it greatly increased root iron (Fe) content, suggesting the role of Ca2+ in the uptake and/or translocation of these minerals. Transcriptomic analysis using RNA-seq revealed that members of Ca2+ channels, including the cyclic nucleotide-gated channels and glutamate receptor-like channels, Ca2+-transporters, Ca2+-binding proteins and Ca2+-dependent protein kinases were differentially regulated by Ca+2 treatment. Moreover, several Fe/metal transporters, including members of vacuolar Fe transporters, yellow stripe-like, natural resistance-associated macrophage protein, and oligo-peptide transporters, were differentially regulated between shoot and root in response to Ca2+ treatment. Taken together, our findings suggest that Ca2+ deficiency affects plant growth and mineral accumulation by regulating the transcriptomes of several transporters and signaling genes.

A resource for integrated genomic analysis of the human liver.

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Genomic map of candidate human imprint control regions: the imprintome.

Imprinted genes - critical for growth, metabolism, and neuronal function - are expressed from one parental allele. Parent-of-origin-dependent CpG methylation regulates this expression at imprint control regions (ICRs). Since ICRs are established before tissue specification, these methylation marks are similar across cell types. Thus, they are attractive for investigating the developmental origins of adult diseases using accessible tissues, but remain unknown. We determined genome-wide candidate ICRs in humans by performing whole-genome bisulphite sequencing (WGBS) of DNA derived from the three germ layers and from gametes. We identified 1,488 hemi-methylated candidate ICRs, including 19 of 25 previously characterized ICRs (https://humanicr.org/). Gamete methylation approached 0% or 100% in 332 ICRs (178 paternally and 154 maternally methylated), supporting parent-of-origin-specific methylation, and 65% were in well-described CTCF-binding or DNaseI hypersensitive regions. This draft of the human imprintome will allow for the systematic determination of the role of early-acquired imprinting dysregulation in the pathogenesis of human diseases and developmental and behavioural disorders.

Meta-analysis of epigenome-wide association studies in newborns and children show widespread sex differences in blood DNA methylation.

BACKGROUND: Among children, sex-specific differences in disease prevalence, age of onset, and susceptibility have been observed in health conditions including asthma, immune response, metabolic health, some pediatric and adult cancers, and psychiatric disorders. Epigenetic modifications such as DNA methylation may play a role in the sexual differences observed in diseases and other physiological traits. METHODS: We performed a meta-analysis of the association of sex and cord blood DNA methylation at over 450,000 CpG sites in 8438 newborns from 17 cohorts participating in the Pregnancy And Childhood Epigenetics (PACE) Consortium. We also examined associations of child sex with DNA methylation in older children ages 5.5-10 years from 8 cohorts (n = 4268). RESULTS: In newborn blood, sex was associated at Bonferroni level significance with differences in DNA methylation at 46,979 autosomal CpG sites (p < 1.3 × 10-7) after adjusting for white blood cell proportions and batch. Most of those sites had lower methylation levels in males than in females. Of the differentially methylated CpG sites identified in newborn blood, 68% (31,727) met look-up level significance (p < 1.1 × 10-6) in older children and had methylation differences in the same direction. CONCLUSIONS: This is a large-scale meta-analysis examining sex differences in DNA methylation in newborns and older children. Expanding upon previous studies, we replicated previous findings and identified additional autosomal sites with sex-specific differences in DNA methylation. Differentially methylated sites were enriched in genes involved in cancer, psychiatric disorders, and cardiovascular phenotypes.

Alterations in DNA methylation associate with fatty liver and metabolic abnormalities in a multi-ethnic cohort of pre-teenage children.

Non-Alcoholic fatty liver disease (NAFLD) is the leading cause of chronic liver disease in children. Epigenetic alterations, such as through DNA methylation (DNAm), may link adverse childhood exposures and fatty liver and provide non-invasive methods for identifying children at high risk for NAFLD and associated metabolic dysfunction. We investigated the association between differential DNAm and liver fat content (LFC) and liver injury in pre-adolescent children. Leveraging data from the Newborn Epigenetics Study (NEST), we enrolled 90    mother-child dyads and used linear regression to identify CpG sites and differentially methylated regions (DMRs) in peripheral blood associated with LFC and alanine aminotransferase (ALT) levels in 7-12yo children. DNAm was measured using Infinium HumanMethylationEPIC BeadChips (Illumina). LFC and fibrosis were quantified by magnetic resonance imaging proton density fat fraction and elastography. Median LFC was 1.4% (range, 0.3-13.4%) and MRE was 2.5 kPa (range, 1.5-3.6kPa). Three children had LFC ≥ 5%, while six (7.6%) met our definition of NAFLD (LFC ≥ 3.7%). All children with NAFLD were obese and five were Black. LFC was associated with 88 DMRs and 106 CpGs (FDR<5%). The top two CpGs, cg25474373 and cg07264203, mapped to or near RFTN2 and PRICKLE2 genes. These two CpG sites were also significantly associated with a NAFLD diagnosis. As higher LFC associates with an adverse cardiometabolic profile already in childhood, altered DNAm may identify these children early in disease course for targeted intervention. Larger, longitudinal studies are needed to validate these findings and determine mechanistic relevance.