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.

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.

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.

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.

Epigenetic Dysregulation of KCNK9 Imprinting and Triple-Negative Breast Cancer

Simple Summary

Genomic imprinting is an inherited form of parent-of-origin specific epigenetic gene regulation that is dysregulated by poor prenatal nutrition and environmental toxins. Here, we showed that KCNK9 is imprinted in breast tissue and identified the differentially methylated region (DMR) controlling its imprint status. Hypomethylation at the DMR, coupled with biallelic expression of KCNK9, occurred in 63% of triple-negative breast cancers (TNBC). The association between hypomethylation and TNBC status was highly significant in African-Americans (p = 0.006), but not in Caucasians (p = 0.70). The high frequency of KCNK9 DMR hypomethylation in TNBC and non-cancerous breast tissue from high-risk women provides evidence that hypomethylation of the KNCK9 DMR/TASK3 overexpression may provide a new target for prevention of TNBC.

Abstract

Genomic imprinting is an inherited form of parent-of-origin specific epigenetic gene regulation that is dysregulated by poor prenatal nutrition and environmental toxins. KCNK9 encodes for TASK3, a pH-regulated potassium channel membrane protein that is overexpressed in 40% of breast cancer. However, KCNK9 gene amplification accounts for increased expression in <10% of these breast cancers. Here, we showed that KCNK9 is imprinted in breast tissue and identified a differentially methylated region (DMR) controlling its imprint status. Hypomethylation at the DMR, coupled with biallelic expression of KCNK9, occurred in 63% of triple-negative breast cancers (TNBC). The association between hypomethylation and TNBC status was highly significant in African-Americans (p = 0.006), but not in Caucasians (p = 0.70). KCNK9 hypomethylation was also found in non-cancerous tissue from 77% of women at high-risk of developing breast cancer. Functional studies demonstrated that the KCNK9 gene product, TASK3, regulates mitochondrial membrane potential and apoptosis-sensitivity. In TNBC cells and non-cancerous mammary epithelial cells from high-risk women, hypomethylation of the KCNK9 DMR predicts for increased TASK3 expression and mitochondrial membrane potential (p < 0.001). This is the first identification of the KCNK9 DMR in mammary epithelial cells and demonstration that its hypomethylation in breast cancer is associated with increases in both mitochondrial membrane potential and apoptosis resistance. The high frequency of hypomethylation of the KCNK9 DMR in TNBC and non-cancerous breast tissue from high-risk women provides evidence that hypomethylation of the KNCK9 DMR/TASK3 overexpression may serve as a marker of risk and a target for prevention of TNBC, particularly in African American women.

Associations between maternal obesity, gestational cytokine levels and child obesity in the NEST cohort

BACKGROUND

Although maternal systemic inflammation is hypothesized to link maternal pre-pregnancy obesity to offspring metabolic dysfunction, patient empirical data are limited.

OBJECTIVES

In this study, we hypothesized that pre-pregnancy obesity alters systemic chemo/cytokines concentrations in pregnancy, and this alteration contributes to obesity in children.

METHODS

In a multi-ethnic cohort of 361 mother-child pairs, we measured prenatal concentrations of plasma TNF-α, IL-6, IL-8, IL-1β, IL-4, IFN-γ, IL-12 p70 subunit, and IL-17A using a multiplex ELISA and examined associations of pre-pregnancy obesity on maternal chemo/cytokine levels, and associations of these cytokine levels with offspring body mass index z score (BMI-z) at age 2-6 years using linear regression.

RESULTS

After adjusting for maternal smoking, ethnicity, age, and education, pre-pregnancy obesity was associated with increased concentrations of TNF-α (P = .026) and IFN-γ (P = .06). While we found no evidence for associations between TNF-α concentrations and offspring BMI-z, increased IFN-γ concentrations were associated with decreased BMI-z (P = .0002), primarily in Whites (P = .0011). In addition, increased maternal IL-17A concentrations were associated with increased BMI-z in offspring (P = .0005) with stronger associations in African Americans (P = .0042) than Whites (P = .24).

CONCLUSIONS

Data from this study are consistent with maternal obesity-related inflammation during pregnancy, increasing the risk of childhood obesity in an ethnic-specific manner.

Effects of Cadmium Exposure on DNA Methylation at Imprinting Control Regions and Genome-Wide in Mothers and Newborn Children

BACKGROUND

Imprinted genes are defined by their preferential expression from one of the two parental alleles. This unique mode of gene expression is dependent on allele-specific DNA methylation profiles established at regulatory sequences called imprinting control regions (ICRs). These loci have been used as biosensors to study how environmental exposures affect methylation and transcription. However, a critical unanswered question is whether they are more, less, or equally sensitive to environmental stressors as the rest of the genome.

OBJECTIVES

Using cadmium exposure in humans as a model, we aimed to determine the relative sensitivity of ICRs to perturbation of methylation compared to similar, nonimprinted loci in the genome.

METHODS

We assayed DNA methylation genome-wide using bisulfite sequencing of 19 newborn cord blood and 20 maternal blood samples selected on the basis of maternal blood cadmium levels. Differentially methylated regions (DMRs) associated with cadmium exposure were identified.

RESULTS

In newborn cord blood and maternal blood, 641 and 1,945 cadmium-associated DMRs were identified, respectively. DMRs were more common at the 15 maternally methylated ICRs than at similar nonimprinted loci in newborn cord blood (p=5.64×10^-8) and maternal blood (p=6.22×10^-14), suggesting a higher sensitivity for ICRs to cadmium. Genome-wide, Enrichr analysis indicated that the top three functional categories for genes that overlapped DMRs in maternal blood were body mass index (BMI) (p=2.0×10^-5), blood pressure (p=3.8×10^-5), and body weight (p=0.0014). In newborn cord blood, the top three functional categories were BMI, atrial fibrillation, and hypertension, although associations were not significant after correction for multiple testing (p=0.098). These findings suggest that epigenetic changes may contribute to the etiology of cadmium-associated diseases.

CONCLUSIONS

We analyzed cord blood and maternal blood DNA methylation profiles genome-wide at nucleotide resolution in individuals selected for high and low blood cadmium levels in the first trimester. Our findings suggest that ICRs may be hot spots for perturbation by cadmium, motivating further study of these loci to investigate potential mechanisms of cadmium action. https://doi.org/10.1289/EHP2085.

Epigenetics, obesity and early-life cadmium or lead exposure

Obesity is a complex and multifactorial disease, which likely comprises multiple subtypes. Emerging data have linked chemical exposures to obesity. As organismal response to environmental exposures includes altered gene expression, identifying the regulatory epigenetic changes involved would be key to understanding the path from exposure to phenotype and provide new tools for exposure detection and risk assessment. In this report, we summarize published data linking early-life exposure to the heavy metals, cadmium and lead, to obesity. We also discuss potential mechanisms, as well as the need for complete coverage in epigenetic screening to fully identify alterations. The keys to understanding how metal exposure contributes to obesity are improved assessment of exposure and comprehensive establishment of epigenetic profiles that may serve as markers for exposures.

Analysis of Imprinted Gene Regulation

Genetic studies have been well established for identifying sequence variants associated with phenotypes. With the expanding field of epigenetics, and the growing understanding of epigenetic regulation of gene expression, similar studies can be undertaken to also define associations between epigenetic variation and phenotypes. Of particular interest are imprinted genes, which have parent-of-origin specific regulation and expression, and are key regulators of early development. Herein, we describe methods for examining epigenetic regulation by the two major hallmarks of imprinted genes: differentially methylated regions (DMRs), regulatory DNA sequences with allele specific methylation; and monoallelic expression, the silencing and transcription of opposite alleles in a parent-of-origin specific manner.

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

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

Adaptive radiation-induced epigenetic alterations mitigated by antioxidants

Humans are exposed to low-dose ionizing radiation (LDIR) from a number of environmental and medical sources. In addition to inducing genetic mutations, there is concern that LDIR may also alter the epigenome. Such heritable effects early in life can either be positively adaptive or result in the enhanced formation of diseases, including cancer, diabetes, and obesity. Herein, we show that LDIR significantly increased DNA methylation at the viable yellow agouti (A(vy)) locus in a sex-specific manner (P=0.004). Average DNA methylation was significantly increased in male offspring exposed to doses between 0.7 and 7.6 cGy, with maximum effects at 1.4 and 3.0 cGy (P<0.01). Offspring coat color was concomitantly shifted toward pseudoagouti (P<0.01). Maternal dietary antioxidant supplementation mitigated both the DNA methylation changes and coat color shift in the irradiated offspring. Thus, LDIR exposure during gestation elicits epigenetic alterations that lead to positive adaptive phenotypic changes that are negated with antioxidants, indicating they are mediated in part by oxidative stress. These findings provide evidence that in the isogenic A(vy) mouse model, epigenetic alterations resulting from LDIR play a role in radiation hormesis, bringing into question the assumption that every dose of radiation is harmful.

DNA double-strand breaks relieve USF-mediated repression of Dβ2 germline transcription in developing thymocytes

Activation of germline promoters is central to V(D)J recombinational accessibility, driving chromatin remodeling, nucleosome repositioning, and transcriptional read-through of associated DNA. We have previously shown that of the two TCRβ locus (Tcrb) D segments, Dβ1 is flanked by an upstream promoter that directs its transcription and recombinational accessibility. In contrast, transcription within the DJβ2 segment cluster is initially restricted to the J segments and only redirected upstream of Dβ2 after D-to-J joining. The repression of upstream promoter activity prior to Tcrb assembly correlates with evidence that suggests DJβ2 recombination is less efficient than that of DJβ1. Because inefficient DJβ2 assembly offers the potential for V-to-DJβ2 recombination to rescue frameshifted V-to-DJβ1 joints, we wished to determine how Dβ2 promoter activity is modulated upon Tcrb recombination. In this study, we show that repression of the otherwise transcriptionally primed 5'Dβ2 promoter requires binding of upstream stimulatory factor (USF)-1 to a noncanonical E-box within the Dβ2 12-recombination signal sequence spacer prior to Tcrb recombination. USF binding is lost from both rearranged and germline Dβ2 sites in DNA-dependent protein kinase, catalytic subunit-competent thymocytes. Finally, genotoxic dsDNA breaks lead to rapid loss of USF binding and gain of transcriptionally primed 5'Dβ2 promoter activity in a DNA-dependent protein kinase, catalytic subunit-dependent manner. Together, these data suggest a mechanism by which V(D)J recombination may feed back to regulate local Dβ2 recombinational accessibility during thymocyte development.

Novel retrotransposed imprinted locus identified at human 6p25

Differentially methylated regions (DMRs) are stable epigenetic features within or in proximity to imprinted genes. We used this feature to identify candidate human imprinted loci by quantitative DNA methylation analysis. We discovered a unique DMR at the 5′-end of FAM50B at 6p25.2. We determined that sense transcripts originating from the FAM50B locus are expressed from the paternal allele in all human tissues investigated except for ovary, in which expression is biallelic. Furthermore, an antisense transcript, FAM50B-AS, was identified to be monoallelically expressed from the paternal allele in a variety of tissues. Comparative phylogenetic analysis showed that FAM50B orthologs are absent in chicken and platypus, but are present and biallelically expressed in opossum and mouse. These findings indicate that FAM50B originated in Therians after divergence from Prototherians via retrotransposition of a gene on the X chromosome. Moreover, our data are consistent with acquisition of imprinting during Eutherian evolution after divergence of Glires from the Euarchonta mammals. FAM50B expression is deregulated in testicular germ cell tumors, and loss of imprinting occurs frequently in testicular seminomas, suggesting an important role for FAM50B in spermatogenesis and tumorigenesis. These results also underscore the importance of accounting for parental origin in understanding the mechanism of 6p25-related diseases.

Identification of chromosome 7 inversion breakpoints in an autistic family narrows candidate region for autism susceptibility

Chromosomal breaks and rearrangements have been observed in conjunction with autism and autistic spectrum disorders. A chromosomal inversion has been previously reported in autistic siblings, spanning the region from approximately 7q22.1 to 7q31. This family is distinguished by having multiple individuals with autism and associated disabilities. The region containing the inversion has been strongly implicated in autism by multiple linkage studies, and has been particularly associated with language defects in autism as well as in other disorders with language components. Mapping of the inversion breakpoints by FISH has localized the inversion to the region spanning approximately 99-108.75 Mb of chromosome 7. The proximal breakpoint has the potential to disrupt either the coding sequence or regulatory regions of a number of cytochrome P450 genes while the distal region falls in a relative gene desert. Copy number variant analysis of the breakpoint regions detected no duplication or deletion that could clearly be associated with disease status. Association analysis in our autism data set using single nucleotide polymorphisms located near the breakpoints showed no significant association with proximal breakpoint markers, but has identified markers near the distal breakpoint ( approximately 108-110 Mb) with significant associations to autism. The chromosomal abnormality in this family strengthens the case for an autism susceptibility gene in the chromosome 7q22-31 region and targets a candidate region for further investigation.

Investigation of potential gene–gene interactions between apoe and reln contributing to autism risk

BACKGROUND

Several candidate gene studies support RELN as susceptibility gene for autism. Given the complex inheritance pattern of autism, it is expected that gene-gene interactions will exist. A logical starting point for examining potential gene-gene interactions is to evaluate the joint effects of genes involved in a common biological pathway. RELN shares a common biological pathway with APOE, and Persico et al. have observed transmission distortion of the APOE-2 allele in autism families.

OBJECTIVE

We evaluated RELN and APOE for joint effects in autism susceptibility.

METHODS

A total of 470 Caucasian autism families were analyzed (265 multiplex; 168 trios with no family history; 37 positive family history but only one sampled affected). These families were genotyped for 11 RELN polymorphisms, including the 5' untranslated region repeat previously associated with autism, as well as for the APOE functional allele. We evaluated single locus allelic and genotypic association with the pedigree disequilibrium test and geno-PDT, respectively. Multilocus effects were evaluated using the extended version of the multifactorial dimensionality reduction method.

RESULTS

For the single locus analyses, there was no evidence for an effect of APOE in our data set. Evidence for association with RELN (rs2,073,559; trio subset P=0.07 PDT; P=0.001 geno-PDT; overall geno-PDT P=0.05), however, was found. For multilocus geno-PDT analysis, the joint genotype of APOE and RELN rs2,073,559 was highly significant (trio subset, global P=0.0001), probably driven by the RELN single locus effect. Using the extended version of the multifactorial dimensionality reduction method to detect multilocus effects, there were no statistically significant associations for any of the n-locus combinations involving RELN or APOE in the overall or multiplex subset. In the trio subset, 1-locus and 2-locus models selected only markers in RELN as best models for predicting autism case status.

CONCLUSION

Thus, we conclude that there is no main effect of APOE in our autism data set, nor is there any evidence for a joint effect of APOE with RELN. RELN, however, remains a good candidate for autism susceptibility.

Analysis of the RELN gene as a genetic risk factor for autism

Several genome-wide screens have indicated the presence of an autism susceptibility locus within the distal long arm of chromosome 7 (7q). Mapping at 7q22 within this region is the candidate gene reelin (RELN). RELN encodes a signaling protein that plays a pivotal role in the migration of several neuronal cell types and in the development of neural connections. Given these neurodevelopmental functions, recent reports that RELN influences genetic risk for autism are of significant interest. The total data set consists of 218 Caucasian families collected by our group, 85 Caucasian families collected by AGRE, and 68 Caucasian families collected at Tufts University were tested for genetic association of RELN variants to autism. Markers included five single-nucleotide polymorphisms (SNPs) and a repeat in the 5'-untranslated region (5'-UTR). Tests for association in Duke and AGRE families were also performed on four additional SNPs in the genes PSMC2 and ORC5L, which flank RELN. Family-based association analyses (PDT, Geno-PDT, and FBAT) were used to test for association of single-locus markers and multilocus haplotypes with autism. The most significant association identified from this combined data set was for the 5'-UTR repeat (PDT P-value=0.002). These analyses show the potential of RELN as an important contributor to genetic risk in autism.

The RNA polymerase II CTD kinase CTDK-I affects pre-mRNA 3' cleavage/polyadenylation through the processing component Pti1p

There are several kinases in Saccharomyces cerevisiae that phosphorylate the CTD of RNA polymerase II, but specific and distinct functions of the phospho-CTDs generated by the different kinases are not well understood. A genetic screen for suppressors of loss of yeast CTD kinase I (CTDK-I) function (by deletion of the catalytic subunit gene CTK1) identified PTI1, a potential 3' cleavage/polyadenylation factor. Genetic and physical interactions connect Pti1p to components of CF IA and CF II/CPF, and mutations of PTI1 or CTK1 affect 3' cleavage site choice and transcript abundance of particular genes. Therefore, one important function of the CTDK-I-generated phospho-CTD appears to be the coupling of transcription to 3' processing of pre-mRNAs by a Pti1p-containing complex.