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.

Population genomics in a disease targeted primary cell model.

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

The proportion of mutations predicted to have a deleterious effect differs between gain and loss of function genes in neurodegenerative disease.

As more studies are turning to bioinformatic prediction programs to assess the potential impact of amino acid substitutions, it is relevant to evaluate the prediction results these programs give in genes that have been well-characterized for Mendelian diseases. Eight genes responsible for neurodegenerative disease with many identified mutations were sub-grouped into those that either have a gain or loss of function disease mechanism. Three prediction programs, PolyPhen, Panther and SIFT, were queried for the reported missense mutations. The mean percent of benign mutations was significantly higher in gain of function genes using the PolyPhen program (38% versus 21%, p=0.007). The probability that a gain of function mutation was predicted to have a damaging role was also significantly less using the Panther program (p=4.86x10(-12)). In contrast, there was no difference between gain and loss of function gene when the SIFT program was used. However, the most accurate distinction between gain and loss of function genes could be obtained when considering the mutations for which all three programs predicted the same result. Further, stratification of SOD1 mutations indicated that only the PolyPhen program could distinguish mutations that impaired enzymatic activity of SOD1 from those with near wildtype activity. The profile of benign and damaging changes from these genes will aid in the interpretation of bioinformatic prediction program results from missense mutations identified in novel genes.

A mutation in the HSN2 gene causes sensory neuropathy type II in a Lebanese family.

Hereditary sensory and autonomic neuropathy (HSAN) type II is an autosomal recessive disorder clinically characterized by distal and proximal sensory loss that is caused by the reduction or absence of peripheral sensory nerves. Recently, a novel gene called HSN2 has been found to be the cause of HSAN type II in five families from Newfoundland and Quebec. Screening of this gene in an HSAN type II Lebanese family showed a 1bp deletion mutation found in a homozygous state in all affected individuals. This novel mutation supports the hypothesis that HSN2 is the causative gene for HSAN type II.

Linkage to the CCM2 locus and genetic heterogeneity in familial cerebral cavernous malformation.

BACKGROUND: Cerebral cavernous malformation (CCM) is a form of intracranial vascular disease that may arise sporadically or be dominantly inherited. Linkage studies have revealed genetic heterogeneity among the dominantly inherited forms suggesting the existence of at least three loci called CCM1, CCM2 and CCM3. METHODS: In the present study, we screened five families with dominantly inherited CCM for CCM1 gene mutations with denaturing high performance liquid chromatography (DHPLC). Then, we performed linkage analysis and haplotyping on these five families using highly polymorphic markers at the candidate CCM loci. RESULTS: None of the five families tested with DHPLC were found to have mutations in the CCM1 gene. Based on haplotyping, we identified three families segregating alleles for CCM2, while two families segregated alleles for CCM3. Using linkage analysis, we could confirm that one family (IFCAS-1) had a positive Lod score of 2.03 (p<0.0001) at the CCM2 locus using marker D7S678. CONCLUSIONS: The present study is the first one to replicate linkage at the CCM2 locus and provides a fifth family identified as such. It also supports the concept of genetic heterogeneity in CCM, identifying four other families that showed no mutations in the CCM1 gene.

Krit1 missense mutations lead to splicing errors in cerebral cavernous malformation.

At least 40% of families affected with cerebral cavernous malformation have a mutation in Krit1. We previously identified two point mutations in Krit1 leading to changes in amino acids (D137G and Q210E) in two different families. Further RNA analysis reveals that both point mutations actually activate cryptic splice-donor sites, causing aberrant splicing and leading to a frameshift and protein truncation. To date, no simple missense mutations have been detected in Krit1.