Anna K. Naumova

Allele-Specific Chromatin Remodeling in the ZPBP2/GSDMB/ORMDL3 Locus Associated with the Risk of Asthma and Autoimmune Disease

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis.

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus–associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region.

Cell culture-induced aberrant methylation of the imprinted IG DMR in human lymphoblastoid cell lines

DNA methylation patterns are often poorly conserved through cell culturing. To determine the effect of cell immortalization and culture on DNA methylation profiles, we analyzed methylation in the differentially methylated regions (DMR) of five imprinted domains: the intergenic (IG) DMR on chromosome 14q32; potassium voltage-gated channel, KQT-like subfamily, member 1, (KCNQ1); small nuclear ribonucleoprotein polypeptide N (SNRPN), mesoderm specific transcript homolog (MEST); and H19 in lymphoblastoid cell lines (LCLs). In the IG DMR we found an aberrant methylation pattern that was consistent through all the cell lines tested, and significantly different from that of noncultured peripheral blood cells. Using a generalized linear mixed model to compare methylation profiles, we show that recently derived LCLs significantly differ from the CEPH LCLs. This implies a gradual cell-culture related deterioration of DNA methylation in the IG DMR with at least two steps that may be identified: loss of methylation at CG sites 1 and 8; and loss of allelic differences in DNA methylation. The IG DMR methylation profile also confirms the high level of clonality of the CEPH LCLs. We conclude that non-transformed primary cells may be less susceptible to epigenetic anomalies and therefore may provide a more accurate reflection of gene expression in vivo.

Novel imprinted transcripts from the Dlk1-Gtl2 intergenic region, Mico1 and Mico1os, show circadian oscillations

Most of the known imprinted genes are assembled into clusters that share common imprinting control regions (ICRs). Non-coding transcripts are often associated with ICRs and implicated in imprinting regulation. We undertook a systematic search for transcripts originating from the Dlk1-Gtl2 intergenic region that contains the ICR for the chromosome 12 imprinted cluster and identified two overlapping transcripts expressed from opposite strands exclusively from the maternal chromosome. These novel imprinted transcripts most likely represent non-coding RNAs and are located telomeric to the IG DMR, extending the proximal boundary of the region of maternal-specific transcription. Their expression is tissue-specific and shows diurnal and circadian oscillations. Therefore, we named these novel transcripts maternal intergenic circadian oscillating 1 (Mico1) and Mico1, opposite strand (Mico1os).

Genomic Imprinting Variations in the Mouse Type 3 Deiodinase Gene Between Tissues and Brain Regions

The Dio3 gene, which encodes for the type 3 deiodinase (D3), controls thyroid hormone (TH) availability. The lack of D3 in mice results in tissue overexposure to TH and a broad neuroendocrine phenotype. Dio3 is an imprinted gene, preferentially expressed from the paternally inherited allele in the mouse fetus. However, heterozygous mice with paternal inheritance of the inactivating Dio3 mutation exhibit an attenuated phenotype when compared with that of Dio3 null mice. To investigate this milder phenotype, the allelic expression of Dio3 was evaluated in different mouse tissues. Preferential allelic expression of Dio3 from the paternal allele was observed in fetal tissues and neonatal brain regions, whereas the biallelic Dio3 expression occurred in the developing eye, testes, and cerebellum and in the postnatal brain neocortex, which expresses a larger Dio3 mRNA transcript. The newborn hypothalamus manifests the highest degree of Dio3 expression from the paternal allele, compared with other brain regions, and preferential allelic expression of Dio3 in the brain relaxed in late neonatal life. A methylation analysis of two regulatory regions of the Dio3 imprinted domain revealed modest but significant differences between tissues, but these did not consistently correlate with the observed patterns of Dio3 allelic expression. Deletion of the Dio3 gene and promoter did not result in significant changes in the tissue-specific patterns of Dio3 allelic expression. These results suggest the existence of unidentified epigenetic determinants of tissue-specific Dio3 imprinting. The resulting variation in the Dio3 allelic expression between tissues likely explains the phenotypic variation that results from paternal Dio3 haploinsufficiency.

Defective imprint resetting in carriers of Robertsonian translocation Rb (8.12)

Meiotic silencing of unsynapsed chromatin (MSUC) occurs in the germ cells of translocation carriers and may cause meiotic arrest and infertility. We hypothesized that if bypassing meiotic checkpoints MSUC may cause epigenetic defects in sperm. We investigated the meiotic behavior of the Robertsonian translocation Rb (8.12) in mice. The unsynapsed 8 and 12 trivalent was associated with the XY body during early and mid-pachynema in heterozygous Rb (8.12) carriers, suggesting possible silencing of pericentromeric genes, such as the Dnmt3a gene. In wild-type mice, DNMT3A protein showed a dramatic accumulation in the nucleus during the mid-pachytene stage and distinct association with the XY body. In translocation carriers, DNMT3A was less abundant in a proportion of pachytene spermatocytes that also had unsynapsed pericentromeric regions of chromosomes 8 and 12. The same mice had incomplete methylation of the imprinted H19 differentially methylated region (DMR) in sperm. We propose that impaired H19 imprint establishment results from lack of synapsis in chromosomes 8 and 12 probably through transient silencing of a chromosome 8 or 12 gene during pachynema. Furthermore, our findings support the notion that imprint establishment at the H19 locus extends into pachynema.

Dynamics of response to asynapsis and meiotic silencing in spermatocytes from Robertsonian translocation carriers.

Failure of homologous synapsis during meiotic prophase triggers transcriptional repression. Asynapsis of the X and Y chromosomes and their consequent silencing is essential for spermatogenesis. However, asynapsis of portions of autosomes in heterozygous translocation carriers may be detrimental for meiotic progression. In fact, a wide range of phenotypic outcomes from meiotic arrest to normal spermatogenesis have been described and the causes of such a variation remain elusive. To better understand the consequences of asynapsis in male carriers of Robertsonian translocations, we focused on the dynamics of recruitment of markers of asynapsis and meiotic silencing at unsynapsed autosomal trivalents in the spermatocytes of Robertsonian translocation carrier mice. Here we report that the enrichment of breast cancer 1 (BRCA1) and histone γH2AX at unsynapsed trivalents declines during the pachytene stage of meiosis and differs from that observed in the sex body. Furthermore, histone variant H3.3S31, which associates with the sex chromosomes in metaphase I/anaphase I spermatocytes, localizes to autosomes in 12% and 31% of nuclei from carriers of one and three translocations, respectively. These data suggest that the proportion of spermatocytes with markers of meiotic silencing of unsynapsed chromatin (MSUC) at trivalents depends on both, the stage of meiosis and the number of translocations. This may explain some of the variability in phenotypic outcomes associated with Robertsonian translocations. In addition our data suggest that the dynamics of response to asynapsis in Robertsonian translocations differs from the response to sex chromosomal asynapsis in the male germ line.

Oocyte heterogeneity with respect to the meiotic silencing of unsynapsed X chromosomes in the XY female mouse

In the XY pachytene spermatocyte, the sex chromosomes do not synapse except for the pseudoautosomal region and become transcriptionally silenced. It has been suggested that the meiotic silencing of unsynapsed chromatin (MSUC) also occurs in oocytes. In the XY sex-reversed female mouse, the sex chromosomes fail to pair in the majority of oocytes and a greater number of oocytes are eliminated during the meiotic prophase compared to the XX female. Yet, many XY oocytes survive to reach the second meiotic metaphase. The goal of our current study was to determine whether the single X chromosome shows the characteristics of asynapsis and meiotic silencing in a proportion of XY oocytes, which can explain the survival of the remaining oocytes. We first examined the accumulation of markers associated with asynapsis or transcriptional silencing, i.e., BRCA1, γH2AX, H3K9me3, and H3K27me3, at the single X chromosome in the XY oocyte. We found that γH2AX and BRCA1 were enriched on the single X chromosome whereas H3K9me3 was not, and H3K27me3 was enriched at all chromosomes in the majority of XY oocytes. We next examined the meiotic silencing of the single X chromosome using enrichment of the X-encoded ATRX protein. On average, ATRX enrichment was lower in XY oocytes than in XX oocytes as expected from its half gene dosage. However, the intensity of ATRX staining in XY oocytes harboring γH2AX domains showed a remarkable heterogeneity. We conclude that MSUC occurs with varying consequences, resulting in a heterogeneous population of oocytes with respect to protein enrichment in the XY female mouse.

Local genotype influences DNA methylation at two asthma-associated regions, 5q31 and 17q21, in a founder effect population

Background Two asthma-associated regions 17q12–q21 and 5q31.1 harbour genes that show strong effect of genotype on expression levels. DNA methylation has an important role in gene regulation; therefore, we examined DNA methylation at promoters of 12 genes from 5q31 and 17q12–q21 regions. Our goal was to determine whether DNA methylation was associated with predisposition to asthma and whether such a relationship was independent from genetic association. Methods Using sodium bisulfite sequencing and pyrosequencing methylation assays, we examined the effect of genotype on DNA methylation in peripheral blood cells from individuals from the Saguenay-Lac-Saint-Jean asthma familial collection and lymphoblastoid cell lines. Results The local genotype influenced methylation levels of solute carrier family 22 (organic 3 cation/carnitine transporter) member 5 (SLC22A5), zona pellucida binding protein 2 (ZPBP2) and gasdermin A (GSDMA) promoter regions. The genotype had a dominant effect on ZPBP2 and GSDMA methylation with lower methylation levels in individuals that carry the asthma-predisposing alleles. Males also had lower methylation at the ZPBP2 promoter than females. We did not observe an effect of asthma status that would be independent of the genotype and the sex effects in the GSDMA, ZPBP2 and SLC22A5 regions; however, GSDMA and ZPBP2 data were suggestive of interaction between asthma and methylation levels in females and SLC22A5 in males. Conclusions The local genotype influences methylation levels at SLC22A5 and ZPBP2 promoters independently of the asthma status. Further studies are necessary to confirm the relationship between GSDMA-ZPBP2 and SLC22A5 methylation and asthma in females and males separately.