Guillermo E. Rivas

Effects of endocrine disruptors on imprinted gene expression in the mouse embryo

Environmental endocrine disruptors (EDs) are synthetic chemicals that resemble natural hormones and are known to cause epigenetic perturbations. EDs have profound effects on development and fertility. Imprinted genes had been identified as susceptible loci to environmental insults by EDs because they are functionally haploid, and because the imprints undergo epigenetic resetting between generations. To screen for possible epigenetic perturbations caused by EDs at imprinted loci, we treated pregnant mice daily between 8.5 and 12.5 days post coitum (dpc) with di-(2-ethylhexyl)-phthalate (DEHP), bisphenol A (BPA), vinclozolin (VZ), or control oil vehicle. After isolating RNA from the placenta, yolk sac, amnion, head, body, heart, liver, lung, stomach, and intestines of 13.5 dpc embryos we measured the allele-specific expression of 38 imprinted transcripts using multiplex single nucleotide primer extension (SNuPE) assays. In this representative data set we identified only a small number of transcripts that exhibited a substantial relaxation of imprinted expression with statistical significance: Slc22a18 with 10% relaxation in the embryo after BPA treatment; Rtl1as with 11 and 16% relaxation in the lung and placenta, respectively after BPA treatment; and Rtl1 with 12% relaxation in the yolk sac after DEHP treatment. Additionally, the standard deviation of allele-specificity increased in various organs after ED treatment for several transcripts including Igf2r, Rasgrf1, Usp29, Slc38a4, and Xist. Our data suggest that the maintenance of strongly biased monoallelic expression of imprinted genes is generally insensitive to EDs in the 13.5 dpc embryo and extra-embryonic organs, but is not immune to those effects.

Polymorphism in the Androgen Receptor and Mammographic Density in Women Taking and Not Taking Estrogen and Progestin Therapy

There is some evidence that women with a higher number of CAG repeat lengths on the androgen receptor (AR) gene have increased breast cancer risk. We evaluated the association between AR-CAG repeat length and mammographic density, a strong breast cancer risk factor, in 404 African-American and Caucasian breast cancer patients. In postmenopausal estrogen progestin therapy users, carriers of the less active AR-CAG had statistically significantly higher mean percentage of density (41.4%) than carriers of the more active AR-CAG (25.7%; P = 0.04). Our results raise the question of whether the number of AR-CAG repeats predicts breast cancer risk in estrogen progestin therapy users.

CTCF-Dependent Chromatin Bias Constitutes Transient Epigenetic Memory of the Mother at the H19-Igf2 Imprinting Control Region in Prospermatogonia

Genomic imprints—parental allele-specific DNA methylation marks at the differentially methylated regions (DMRs) of imprinted genes—are erased and reestablished in germ cells according to the individuals sex. Imprint establishment at paternally methylated germ line DMRs occurs in fetal male germ cells. In prospermatogonia, the two unmethylated alleles exhibit different rates of de novo methylation at the H19/Igf2 imprinting control region (ICR) depending on parental origin. We investigated the nature of this epigenetic memory using bisulfite sequencing and allele-specific ChIP–SNuPE assays. We found that the chromatin composition in fetal germ cells was biased at the ICR between the two alleles with the maternally inherited allele exhibiting more H3K4me3 and less H3K9me3 than the paternally inherited allele. We determined genetically that the chromatin bias, and also the delayed methylation establishment in the maternal allele, depended on functional CTCF insulator binding sites in the ICR. Our data suggest that, in primordial germ cells, maternally inherited allele-specific CTCF binding sets up allele-specific chromatin differences at the ICR. The erasure of these allele-specific chromatin marks is not complete before the process of de novo methylation imprint establishment begins. CTCF–dependent allele-specific chromatin composition imposes a maternal allele-specific delay on de novo methylation imprint establishment at the H19/Igf2 ICR in prospermatogonia.

Allele-Specific H3K79 Di- versus Trimethylation Distinguishes Opposite Parental Alleles at Imprinted Regions

ABSTRACT Imprinted gene expression corresponds to parental allele-specific DNA CpG methylation and chromatin composition. Histone tail covalent modifications have been extensively studied, but it is not known whether modifications in the histone globular domains can also discriminate between the parental alleles. Using multiplex chromatin immunoprecipitation-single nucleotide primer extension (ChIP-SNuPE) assays, we measured the allele-specific enrichment of H3K79 methylation and H4K91 acetylation along the H19/Igf2 imprinted domain. Whereas H3K79me1, H3K79me2, and H4K91ac displayed a paternal-specific enrichment at the paternally expressed Igf2 locus, H3K79me3 was paternally biased at the maternally expressed H19 locus, including the paternally methylated imprinting control region (ICR). We found that these allele-specific differences depended on CTCF binding in the maternal ICR allele. We analyzed an additional 11 differentially methylated regions (DMRs) and found that, in general, H3K79me3 was associated with the CpG-methylated alleles, whereas H3K79me1, H3K79me2, and H4K91ac enrichment was specific to the unmethylated alleles. Our data suggest that allele-specific differences in the globular histone domains may constitute a layer of the “histone code” at imprinted genes.

Coordinated allele-specific histone acetylation at the differentially methylated regions of imprinted genes

Genomic imprinting is an epigenetic inheritance system characterized by parental allele-specific gene expression. Allele-specific DNA methylation and chromatin composition are two epigenetic modification systems that control imprinted gene expression. To get a general assessment of histone lysine acetylation at imprinted genes we measured allele-specific acetylation of a wide range of lysine residues, H3K4, H3K18, H3K27, H3K36, H3K79, H3K64, H4K5, H4K8, H4K12, H2AK5, H2BK12, H2BK16 and H2BK46 at 11 differentially methylated regions (DMRs) in reciprocal mouse crosses using multiplex chromatin immunoprecipitation SNuPE assays. Histone acetylation marks generally distinguished the methylation-free alleles from methylated alleles at DMRs in mouse embryo fibroblasts and embryos. Acetylated lysines that are typically found at transcription start sites exhibited stronger allelic bias than acetylated histone residues in general. Maternally methylated DMRs, that usually overlap with promoters exhibited higher levels of acetylation and a 10% stronger allele-specific bias than paternally methylated DMRs that reside in intergenic regions. Along the H19/Igf2 imprinted domain, allele-specific acetylation at each lysine residue depended on functional CTCF binding sites in the imprinting control region. Our results suggest that many different histone acetyltransferase and histone deacetylase enzymes must act in concert in setting up and maintaining reciprocal parental allelic histone acetylation at DMRs.

Strong Signature of Natural Selection within an FHIT Intron Implicated in Prostate Cancer Risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, re-sequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearsons χ2 = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearsons χ2 = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (π = 0.0072, Tajimas D = 3.31, 14 SNPs) and the Japanese (π = 0.0049, Fay & Wus H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wus H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer.

Analysis of candidate colitis genes in the Gdac1 locus of mice deficient in glutathione peroxidase-1 and -2.

Background Mice that are deficient for glutathione peroxidases 1 and 2 (GPX) show large variations in the penetrance and severity of colitis in C57BL/6J and 129S1/SvImJ backgrounds. We mapped a locus contributing to this difference to distal chromosome 2 (∼119–133 mbp) and named it glutathione peroxidase-deficiency-associated colitis 1 (Gdac1). The aim of this study was to identify the best gene candidates within the Gdac1 locus contributing to the murine colitis phenotype. Method/Principal Findings We refined the boundaries of Gdac1 to 118–125 mbp (95% confidence interval) by increasing sample size and marker density across the interval. The narrowed region contains 128 well-annotated protein coding genes but it excludes Fermt1, a human inflammatory bowel disease candidate that was within the original boundaries of Gdac1. The locus we identified may be the Cdcs3 locus mapped by others studying IL10-knockout mice. Using in silico analysis of the 128 genes, based on published colon expression data, the relevance of pathways to colitis, gene mutations, presence of non-synonymous-single-nucleotide polymorphisms (nsSNPs) and whether the nsSNPs are predicted to have an impact on protein function or expression, we excluded 42 genes. Based on a similar analysis, twenty-five genes from the remaining 86 genes were analyzed for expression-quantitative-trait loci, and another 15 genes were excluded. Conclusion/Significance Among the remaining 10 genes, we identified Pla2g4f and Duox2 as the most likely colitis gene candidates, because GPX metabolizes PLA2G4F and DUOX2 products. Pla2g4f is a phospholipase A2 that has three potentially significant nsSNP variants and showed expression differences across mouse strains. PLA2G4F produces arachidonic acid, which is a substrate for lipoxygenases and, in turn, for GPXs. DUOX2 produces H2O2 and may control microbial populations. DUOX-1 and -2 control microbial populations in mammalian lung and in the gut of several insects and zebrafish. Dysbiosis is a phenotype that differentiates 129S1/SvImJ from C57BL/6J and may be due to strain differences in DUOX2 activity.

MIRA-SNuPE, a quantitative, multiplex method for measuring allele-specific DNA methylation

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss-of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.