Regulatory interaction between the ZPBP2-ORMDL3/Zpbp2-Ormdl3 region and the circadian clock.

Genome-wide association study (GWAS) loci for several immunity-mediated diseases (early onset asthma, inflammatory bowel disease (IBD), primary biliary cholangitis, and rheumatoid arthritis) map to chromosomal region 17q12-q21. The predominant view is that association between 17q12-q21 alleles and increased risk of developing asthma or IBD is due to regulatory variants. ORM sphingolipid biosynthesis regulator (ORMDL3) residing in this region is the most promising gene candidate for explaining association with disease. However, the relationship between 17q12-q21 alleles and disease is complex suggesting contributions from other factors, such as trans-acting genetic and environmental modifiers or circadian rhythms. Circadian rhythms regulate expression levels of thousands of genes and their dysregulation is implicated in the etiology of several common chronic inflammatory diseases. However, their role in the regulation of the 17q12-q21 genes has not been investigated. Moreover, the core clock gene nuclear receptor subfamily 1, group D, member 1 (NR1D1) resides about 200 kb distal to the GWAS region. We hypothesized that circadian rhythms influenced gene expression levels in 17q12-q21 region and conversely, regulatory elements in this region influenced transcription of the core clock gene NR1D1 in cis. To test these hypotheses, we examined the diurnal expression profiles of zona pellucida binding protein 2 (ZPBP2/Zpbp2), gasdermin B (GSDMB), and ORMDL3/Ormdl3 in human and mouse tissues and analyzed the impact of genetic variation in the ZPBP2/Zpbp2 region on NR1D1/Nr1d1 expression. We found that Ormdl3 and Zpbp2 were controlled by the circadian clock in a tissue-specific fashion. We also report that deletion of the Zpbp2 region altered the expression profile of Nr1d1 in lungs and ileum in a time-dependent manner. In liver, the deletion was associated with enhanced expression of Ormdl3. We provide the first evidence that disease-associated genes Zpbp2 and Ormdl3 are regulated by circadian rhythms and the Zpbp2 region influences expression of the core clock gene Nr1d1.

Loss of the zona pellucida-binding protein 2 (Zpbp2) gene in mice impacts airway hypersensitivity and lung lipid metabolism in a sex-dependent fashion.

The human chromosomal region 17q12-q21 is one of the best replicated genome-wide association study loci for childhood asthma. The associated SNPs span a large genomic interval that includes several protein-coding genes. Here, we tested the hypothesis that the zona pellucida-binding protein 2 (ZPBP2) gene residing in this region contributes to asthma pathogenesis using a mouse model. We tested the lung phenotypes of knock-out (KO) mice that carry a deletion of the Zpbp2 gene. The deletion attenuated airway hypersensitivity (AHR) in female, but not male, mice in the absence of allergic sensitization. Analysis of the lipid profiles of their lungs showed that female, but not male, KO mice had significantly lower levels of sphingosine-1-phosphate (S1P), very long-chain ceramides (VLCCs), and higher levels of long-chain ceramides compared to wild-type controls. Furthermore, in females, lung resistance following methacholine challenge correlated with lung S1P levels (Pearson correlation coefficient 0.57) suggesting a link between reduced AHR in KO females, Zpbp2 deletion, and S1P level regulation. In livers, spleens and blood plasma, however, VLCC, S1P, and sphingosine levels were reduced in both KO females and males. We also find that the Zpbp2 deletion was associated with gain of methylation in the adjacent DNA regions. Thus, we demonstrate that the mouse ortholog of ZPBP2 has a role in controlling AHR in female mice. Our data also suggest that Zpbp2 may act through regulation of ceramide metabolism. These findings highlight the importance of phospholipid metabolism for sexual dimorphism in AHR.

X chromosome dosage and presence of SRY shape sex-specific differences in DNA methylation at an autosomal region in human cells.

BACKGROUND: Sexual dimorphism in DNA methylation levels is a recurrent epigenetic feature in different human cell types and has been implicated in predisposition to disease, such as psychiatric and autoimmune disorders. To elucidate the genetic origins of sex-specific DNA methylation, we examined DNA methylation levels in fibroblast cell lines and blood cells from individuals with different combinations of sex chromosome complements and sex phenotypes focusing on a single autosomal region--the differentially methylated region (DMR) in the promoter of the zona pellucida binding protein 2 (ZPBP2) as a reporter. RESULTS: Our data show that the presence of the sex determining region Y (SRY) was associated with lower methylation levels, whereas higher X chromosome dosage in the absence of SRY led to an increase in DNA methylation levels at the ZPBP2 DMR. We mapped the X-linked modifier of DNA methylation to the long arm of chromosome X (Xq13-q21) and tested the impact of mutations in the ATRX and RLIM genes, located in this region, on methylation levels. Neither ATRX nor RLIM mutations influenced ZPBP2 methylation in female carriers. CONCLUSIONS: We conclude that sex-specific methylation differences at the autosomal locus result from interaction between a Y-linked factor SRY and at least one X-linked factor that acts in a dose-dependent manner.

Role of DNA methylation in expression control of the IKZF3-GSDMA region in human epithelial cells.

Chromosomal region 17q12-q21 is associated with asthma and harbors regulatory polymorphisms that influence expression levels of all five protein-coding genes in the region: IKAROS family zinc finger 3 (Aiolos) (IKZF3), zona pellucida binding protein 2 (ZPBP2), ORMDL sphingolipid biosynthesis regulator 3 (ORMDL3), and gasdermins A and B (GSDMA, GSDMB). Furthermore, DNA methylation in this region has been implicated as a potential modifier of the genetic risk of asthma development. To further characterize the effect of DNA methylation, we examined the impact of treatment with DNA methyltransferase inhibitor 5-aza-2'-deoxycytidine (5-aza-dC) that causes DNA demethylation, on expression and promoter methylation of the five 17q12-q21 genes in the human airway epithelium cell line NuLi-1, embryonic kidney epithelium cell line 293T and human adenocarcinoma cell line MCF-7. 5-aza-dC treatment led to upregulation of expression of GSDMA in all three cell lines. ZPBP2 was upregulated in NuLi-1, but remained repressed in 293T and MCF-7 cells, whereas ORMDL3 was upregulated in 293T and MCF-7 cells, but not NuLi-1. Upregulation of ZPBP2 and GSDMA was accompanied by a decrease in promoter methylation. Moreover, 5-aza-dC treatment modified allelic expression of ZPBP2 and ORMDL3 suggesting that different alleles may respond differently to treatment. We also identified a polymorphic CTCF-binding site in intron 1 of ORMDL3 carrying a CG SNP rs4065275 and determined its methylation level. The site's methylation was unaffected by 5-aza-dC treatment in NuLi-1 cells. We conclude that modest changes (8-13%) in promoter methylation levels of ZPBP2 and GSDMA may cause substantial changes in RNA levels and that allelic expression of ZPBP2 and ORMDL3 is mediated by DNA methylation.

Robertsonian translocations modify genomic distribution of γH2AFX and H3.3 in mouse germ cells.

Heterozygosity for Robertsonian translocations hampers pairing and synapsis between the translocated chromosome and its normal homologs during meiotic prophase I. This causes meiotic silencing of unsynapsed chromatin in pericentromeric regions. Several lines of evidence suggest that autosomal asynapsis leads to meiotic arrest in males and two underlying mechanisms have been proposed: (1) reactivation of the X and Y chromosomes due to competition for silencing factors and (2) meiotic silencing of genes that are located in the unsynapsed regions and are essential for meiotic progression. The latter mechanism requires that asynapsis and meiotic silencing spread beyond the p-arms of the normal homologs into gene-rich regions. We used chromatin immunoprecipitation assays to determine whether histones γH2AFX and H3.3, both marks of asynapsis and meiotic silencing, are enriched in gene-rich regions of the translocated chromosomes and their homologs in the spermatocytes of heterozygous carriers of Robertsonian translocations. We also asked if γH2AFX and H3.3 enrichment was reduced at the X chromosome and if γH2AFX and H3.3 enrichment was higher on the normal homolog. Our data show that γH2AFX enrichment extends as far as 9-15 Mb of the annotated genomic sequence of the q-arms of the translocated chromosomal trivalents and that both γH2AFX and H3.3 levels are reduced over the X chromosome. Our data are also suggestive of an asymmetry in γH2AFX and H3.3 enrichment with a bias toward the non-translocated homolog.

Sex- and age-dependent DNA methylation at the 17q12-q21 locus associated with childhood asthma.

Chromosomal region 17q12-q21 is one of the best-replicated genome-wide association study (GWAS) hits and associated with childhood-onset asthma. However, the mechanism by which the genetic association is restricted to childhood-onset disease is unclear. During childhood, more boys than girls develop asthma. Therefore, we tested the hypothesis that the 17q12-q21 genetic association was sex-specific. Indeed, a TDT test showed that in the Saguenay-Lac-Saint-Jean familial collection, the 17q12-q21 association was significant among male, but not among female asthmatic subjects. We next hypothesized that the bias in the genetic association resulted from sex-specific and/or age-dependent DNA methylation at regulatory regions and determined the methylation profiles of five 17q12-q21 gene promoters using the bisulfite sequencing methylation assay. We identified a single regulatory region within the zona pellucida binding protein 2 (ZPBP2) gene, which showed statistically significant differences between males and females with respect to DNA methylation. DNA methylation also varied with age and was higher in adult males compared to boys. We have recently identified two functionally important polymorphisms, both within the ZPBP2 gene that influence expression levels of neighboring genes. Combined with the results of the present work, these data converge pointing to the same 5 kb region within the ZPBP2 gene as a critical region for both gene expression regulation and predisposition to asthma. Our data show that sex- and age-dependent DNA methylation may act as a modifier of genetic effects and influence the results of genetic association studies.

DNA methyltransferase 1 (Dnmt1) mutation affects Snrpn imprinting in the mouse male germ line.

DNA methylation and DNA methyltransferases are essential for spermatogenesis. Mutations in the DNA methyltransferase Dnmt1 gene exert a paternal effect on epigenetic states and phenotypes of offspring, suggesting that DNMT1 is important for the epigenetic remodeling of the genome that takes place during spermatogenesis. However, the specific role of DNMT1 in spermatogenesis and the establishment of genomic imprints in the male germ line remains elusive. To further characterize the effect of DNMT1 deficiency on the resetting of methylation imprints during spermatogenesis, we analyzed the methylation profiles of imprinted regions in the spermatozoa of mice that were heterozygous for a Dnmt1 loss-of-function mutation. The mutation did not affect the H19 or IG differentially methylated regions (DMRs) that are usually highly methylated but led to a partial hypermethylation of the Snrpn DMR, a region that should normally be unmethylated in mature spermatozoa. This defect does not appear in mouse models with mutations in Dnmt3a and Mthfr genes and, therefore, it is specific for the Dnmt1 gene and is suggestive of a role of DNMT1 in imprint resetting or maintenance in the male germ line.

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.

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.

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).

Parental effect of DNA (Cytosine-5) methyltransferase 1 on grandparental-origin-dependent transmission ratio distortion in mouse crosses and human families.

Transmission ratio distortion (TRD) is a deviation from the expected Mendelian 1:1 ratio of alleles transmitted from parents to offspring and may arise by different mechanisms. Earlier we described a grandparental-origin-dependent sex-of-offspring-specific TRD of maternal chromosome 12 alleles closely linked to an imprinted region and hypothesized that it resulted from imprint resetting errors in the maternal germline. Here, we report that the genotype of the parents for loss-of-function mutations in the Dnmt1 gene influences the transmission of grandparental chromosome 12 alleles. More specifically, maternal Dnmt1 mutations restore Mendelian transmission ratios of chromosome 12 alleles. Transmission of maternal alleles depends upon the presence of the Dnmt1 mutation in the mother rather than upon the Dnmt1 genotype of the offspring. Paternal transmission mirrors the maternal one: live-born offspring of wild-type fathers display 1:1 transmission ratios, whereas offspring of heterozygous Dnmt1 mutant fathers tend to inherit grandpaternal alleles. Analysis of allelic transmission in the homologous region of human chromosome 14q32 detected preferential transmission of alleles from the paternal grandfather to grandsons. Thus, parental Dnmt1 is a modifier of transmission of alleles at an unlinked chromosomal region and perhaps has a role in the genesis of TRD.