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.

A probabilistic approach for SNP discovery in high-throughput human resequencing data.

New high-throughput sequencing technologies are generating large amounts of sequence data, allowing the development of targeted large-scale resequencing studies. For these studies, accurate identification of polymorphic sites is crucial. Heterozygous sites are particularly difficult to identify, especially in regions of low coverage. We present a new strategy for identifying heterozygous sites in a single individual by using a machine learning approach that generates a heterozygosity score for each chromosomal position. Our approach also facilitates the identification of regions with unequal representation of two alleles and other poorly sequenced regions. The availability of confidence scores allows for a principled combination of sequencing results from multiple samples. We evaluate our method on a gold standard data genotype set from HapMap. We are able to classify sites in this data set as heterozygous or homozygous with 98.5% accuracy. In de novo data our probabilistic heterozygote detection ("ProbHD") is able to identify 93% of heterozygous sites at a <5% false call rate (FCR) as estimated based on independent genotyping results. In direct comparison of ProbHD with high-coverage 1000 Genomes sequencing available for a subset of our data, we observe >99.9% overall agreement for genotype calls and close to 90% agreement for heterozygote calls. Overall, our data indicate that high-throughput resequencing of human genomic regions requires careful attention to systematic biases in sample preparation as well as sequence contexts, and that their impact can be alleviated by machine learning-based sequence analyses allowing more accurate extraction of true DNA variants.

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.

A cis-acting regulatory variant in the IL2RA locus.

The mechanism for the association of type 1 diabetes (T1D) with IL2RA remains to be clarified. Neither of the two distinct, transmission-disequilibrium confirmed loci mapping to this gene can be explained by a coding variant. An effect on the levels of the soluble protein product sIL-2RA has been reported but its cause and relationship to disease risk is not clear. To look for an allelic effect on IL2RA transcription in cis, we examined RNA from 48 heterozygous lymphocyte samples for differential allele expression. Of the 48 samples, 32 showed statistically significant allelic imbalance. No known single nucleotide polymorphism (SNP) had perfect correlation with this transcriptional effect but the one that showed the most significant (p = 1.6 x 10(-5)) linkage disequilibrium with it was the SNP rs3118470. We had previously shown rs3118470 to confer T1D susceptibility in a Canadian dataset, independently of rs41295061 as the major reported locus (p = 5 x 10(-3), after accounting for rs41295061 by conditional regression). Lower IL2RA levels consistently originated from the T1D predisposing allele. We conclude that an as yet unidentified variant or haplotype, best marked by rs3118470, is responsible for this independent effect and increases T1D risk through diminished expression of the IL-2R, likely by interfering with the proper development of regulatory T cells.

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 novel locus for idiopathic generalized epilepsy in French-Canadian families maps to 10p11.

Idiopathic generalized epilepsy (IGE) has evidence of a strong genetic etiology. We conducted genomewide linkage analysis for genes responsible for familial IGE in French-Canadian pedigrees. Twenty families segregating autosomal dominant epilepsy were collected. Four larger IGE families sufficiently powerful for independent linkage analysis were genome-scanned and follow-up fine mapping was performed over regions with LOD scores >3.0. The genotyping of 16 smaller families was carried out at significantly linked loci for supportive linkage analysis and haplotype comparisons. One of the four families provided a significant linkage result at marker D10S1426 on chromosome 10 (two-point LOD score = 3.05, theta = 0, multipoint LOD score = 3.18). Fine mapping revealed a segregating haplotype and key recombination breakpoints, suggesting a candidate gene interval of 6.5 Mb. Multipoint linkage analyses using the additional 16 families yielded a maximum LOD score under heterogeneity of 4.23 (alpha = 0.34) at this locus. Evaluation of recombination breakpoints in these families narrowed the candidate region to 1.7 Mb. Sequencing of the two known genes in this region, NRP1 and PARD3, was negative for mutation. Replication of linkage to this locus in other cohorts of IGE families is essential to characterize the underlying genetic mechanism for the disease.

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.

Genome-wide assessment of imprinted expression in human cells.

BACKGROUND: Parent-of-origin-dependent expression of alleles, imprinting, has been suggested to impact a substantial proportion of mammalian genes. Its discovery requires allele-specific detection of expressed transcripts, but in some cases detected allelic expression bias has been interpreted as imprinting without demonstrating compatible transmission patterns and excluding heritable variation. Therefore, we utilized a genome-wide tool exploiting high density genotyping arrays in parallel measurements of genotypes in RNA and DNA to determine allelic expression across the transcriptome in lymphoblastoid cell lines (LCLs) and skin fibroblasts derived from families. RESULTS: We were able to validate 43% of imprinted genes with previous demonstration of compatible transmission patterns in LCLs and fibroblasts. In contrast, we only validated 8% of genes suggested to be imprinted in the literature, but without clear evidence of parent-of-origin-determined expression. We also detected five novel imprinted genes and delineated regions of imprinted expression surrounding annotated imprinted genes. More subtle parent-of-origin-dependent expression, or partial imprinting, could be verified in four genes. Despite higher prevalence of monoallelic expression, immortalized LCLs showed consistent imprinting in fewer loci than primary cells. Random monoallelic expression has previously been observed in LCLs and we show that random monoallelic expression in LCLs can be partly explained by aberrant methylation in the genome. CONCLUSIONS: Our results indicate that widespread parent-of-origin-dependent expression observed recently in rodents is unlikely to be captured by assessment of human cells derived from adult tissues where genome-wide assessment of both primary and immortalized cells yields few new imprinted loci.

Targeted screening of cis-regulatory variation in human haplotypes.

Regulatory cis-acting variants account for a large proportion of gene expression variability in populations. Cis-acting differences can be specifically measured by comparing relative levels of allelic transcripts within a sample. Allelic expression (AE) mapping for cis-regulatory variant discovery has been hindered by the requirements of having informative or heterozygous single nucleotide polymorphisms (SNPs) within genes in order to assign the allelic origin of each transcript. In this study we have developed an approach to systematically screen for heritable cis-variants in common human haplotypes across >1,000 genes. In order to achieve the highest level of information per haplotype studied, we carried out allelic expression measurements by using both intronic and exonic SNPs in primary transcripts. We used a novel RNA pooling strategy in immortalized lymphoblastoid cell lines (LCLs) and primary human osteoblast cell lines (HObs) to allow for high-throughput AE. Screening hits from RNA pools were further validated by performing allelic expression mapping in individual samples. Our results indicate that >10% of expressed genes in human LCLs show genotype-linked AE. In addition, we have validated cis-acting variants in over 20 genes linked with common disease susceptibility in recent genome-wide studies. More generally, our results indicate that RNA pooling coupled with AE read-out by second generation sequencing or by other methods provides a high-throughput tool for cataloging the impact of common noncoding variants in the human genome.

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.

A novel duplication confirms the involvement of 5q23.2 in autosomal dominant leukodystrophy.

OBJECTIVE: To identify the underlying locus and disease-causing mutation for adult-onset autosomal dominant leukodystrophy (ADLD). DESIGN: Previously, an adult-onset ADLD locus on chromosome 5q23 was mapped between markers D5S1495 and CTT/CCT15. This region contains 13 known and putative candidate genes. A 2-point linkage analysis confirmed linkage of a large multigenerational French Canadian family to chromosome 5q23. In addition, screening of the 13 genes within the candidate interval as well as 5 neighboring genes was completed, followed by comparative genomic hybridization. SUBJECTS: A multigenerational French Canadian family with ADLD mimicking progressive multiple sclerosis was identified and studied. Eight affected family members were available for the study and presented with autonomic dysfunction as well as upper motorneuron signs affecting gait. RESULTS: The thorough candidate gene approach did not identify any mutation. Consequently, a whole-chromosome comparative genomic hybridization for chromosome 5 identified a 280-kilobase duplication within the chromosomal band 5q23.2 in 2 affected individuals. This duplication contains 3 genes: LMNB1, FLJ36242, and MARCH3. CONCLUSION: We have identified a novel duplication on chromosomal band 5q23.2 in a French Canadian family with ADLD that supports the implication of duplicated LMNB1 as the disease-causing mutation. However, additional functional studies of lamin B1 overexpression are necessary to elucidate the involvement of lamin B1 in myelination and in degenerative disorders such as ADLD and multiple sclerosis.

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.

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.

The 14q restless legs syndrome locus in the French Canadian population.

A new restless legs syndrome locus on chromosome 14 recently has been reported in one family of Italian origin. Our study aimed to replicate this finding and determine the importance of this locus in the French Canadian population. Markers spanning the region were genotyped in 14 large families and linkage assessed using two-point and multipoint logarithm of odds scores. Possible linkage to this locus was found in one of our kindreds providing support for the existence of this locus and indicating that this locus may be responsible for a small fraction of French Canadian restless legs syndrome.