In vitro assays fail to predict in vivo effects of regulatory polymorphisms

Highly conserved non-coding elements on either side of SOX9 associated with Pierre Robin sequence

Pierre Robin sequence (PRS) is an important subgroup of cleft palate. We report several lines of evidence for the existence of a 17q24 locus underlying PRS, including linkage analysis results, a clustering of translocation breakpoints 1.06-1.23 Mb upstream of SOX9, and microdeletions both approximately 1.5 Mb centromeric and approximately 1.5 Mb telomeric of SOX9. We have also identified a heterozygous point mutation in an evolutionarily conserved region of DNA with in vitro and in vivo features of a developmental enhancer. This enhancer is centromeric to the breakpoint cluster and maps within one of the microdeletion regions. The mutation abrogates the in vitro enhancer function and alters binding of the transcription factor MSX1 as compared to the wild-type sequence. In the developing mouse mandible, the 3-Mb region bounded by the microdeletions shows a regionally specific chromatin decompaction in cells expressing Sox9. Some cases of PRS may thus result from developmental misexpression of SOX9 due to disruption of very-long-range cis-regulatory elements.

TNFSF15 transcripts from risk haplotype for Crohn's disease are overexpressed in stimulated T cells

TNFSF15 is a susceptibility gene for Crohn's disease (CD). It remains to be elucidated how the associated single nucleotide polymorphisms (SNPs) in TNFSF15 affect the susceptibility to CD. Because there are no non-synonymous SNPs in TNFSF15, we speculated that one or more of the SNPs associated with CD may act as cis-regulatory SNPs. To reveal the effects of the SNPs on the transcriptional activity of TNFSF15, we first examined the allelic expression imbalance of TNFSF15 in peripheral blood mononuclear cells (PBMCs). When PBMCs stimulated by phytohemagglutinin (PHA) were examined, the allelic ratio of mRNA transcribed from the risk haplotype to the non-risk haplotype increased, compared with the ratio without stimulation. When peripheral blood T cells and Jurkat cells stimulated by phorbol 12-myristate 13-acetate + ionomycin were examined, an allelic expression imbalance similar to that observed in PBMCs stimulated by PHA was confirmed. The promoter assay in stimulated Jurkat cells showed that the luciferase activity of the promoter region (-979 to +35) of the risk haplotype was significantly higher than that of the non-risk haplotype, and deletion and mutagenesis analysis demonstrated that this difference resulted from the -358T/C SNP. The promoter activity of -358C (risk allele) was higher than that of -358T (non-risk allele) in stimulated T cells. This effect of -358T/C on the transcriptional activity in stimulated T cells may confer susceptibility to CD.

PPP2R2B CAG repeat length in the Han Chinese in Taiwan: Association analyses in neurological and psychiatric disorders and potential functional implications

PPP2R2B, a protein widely expressed in neurons throughout the brain, regulates the protein phosphatase 2A (PP2A) activity for the microtubule-associated protein tau and other substrates. Altered PP2A activity has been implicated in spinocerebellar ataxia 12, Alzheimer's disease (AD), and other tauopathies. Through a case-control study and a reporter assay, we investigated the association of PPP2R2B CAG repeat polymorphism with Taiwanese AD, essential tremor (ET), Parkinson's disease (PD), and schizophrenia and clarified the functional implication of this polymorphism. The distribution of the alleles was not significantly different between patients and controls, with 68.6-76.1% alleles at lengths of 10, 13, and 16 triplets. No expanded alleles were detected in either group. However, the frequency of the individuals carrying the short 5-, 6-, and 7-triplet alleles was notably higher in patients with AD (5/180 [2.8%], Fisher's exact test, P = 0.003; including 2 homozygotes) and ET (4/132 [3.0%], Fisher's exact test, P < 0.001) than in the controls (1/625 [0.2%]). The PPP2R2B transcriptional activity was significantly lower in the luciferase reporter constructs containing the (CAG)(5-7) allele than in those containing the common 10-, 13-, and 16-triplet alleles in both neuroblastoma and embryonic kidney cells. Therefore, our preliminary results suggest that the PPP2R2B gene CAG repeat polymorphism may be functional and may, in part, play a role in conferring susceptibility to AD and ET in Taiwan.

Disruption of an AP-2alpha binding site in an IRF6 enhancer is associated with cleft lip

Previously we have shown that nonsyndromic cleft lip with or without cleft palate (NSCL/P) is strongly associated with SNPs in IRF6 (interferon regulatory factor 6). Here, we use multispecies sequence comparisons to identify a common SNP (rs642961, G>A) in a newly identified IRF6 enhancer. The A allele is significantly overtransmitted (P = 1 x 10(-11)) in families with NSCL/P, in particular those with cleft lip but not cleft palate. Further, there is a dosage effect of the A allele, with a relative risk for cleft lip of 1.68 for the AG genotype and 2.40 for the AA genotype. EMSA and ChIP assays demonstrate that the risk allele disrupts the binding site of transcription factor AP-2alpha and expression analysis in the mouse localizes the enhancer activity to craniofacial and limb structures. Our findings place IRF6 and AP-2alpha in the same developmental pathway and identify a high-frequency variant in a regulatory element contributing substantially to a common, complex disorder.

Polymorphisms affecting gene transcription and mRNA processing in pharmacogenetic candidate genes: detection through allelic expression imbalance in human target tissues

BACKGROUND

Genetic variation in mRNA expression plays a critical role in human phenotypic diversity, but it has proven difficult to detect regulatory polymorphisms - mostly single nucleotide polymorphisms (rSNPs). Additionally, variants in the transcribed region, termed here 'structural RNA SNPs' (srSNPs), can affect mRNA processing and turnover. Both rSNPs and srSNPs cause allelic mRNA expression imbalance (AEI) in heterozygous individuals. We have used AEI to discover and characterize regulatory polymorphisms in OPRM1, TPH2, MDR1, DRD2, and VKORC1. The objective of this study was to use AEI to determine the extent of cis-regulatory factors in pharmacogenetic genes.

METHODS

We applied a rapid and accurate AEI methodology for testing 42 genes implicated in cardiovascular and central nervous system diseases, and affecting drug metabolism and transport. Each gene was analyzed in physiologically relevant human autopsy tissues, including brain, heart, liver, intestines, and lymphocytes.

RESULTS

Substantial AEI was observed in approximately 55% of the surveyed genes. Focusing on cardiovascular candidate genes in human hearts, AEI analysis revealed frequent cis-acting regulatory factors in ACE and SOD2 mRNA expression, having potential clinical significance. SNP scanning to locate regulatory polymorphisms in a number of genes failed to support several previously proposed promoter SNPs discovered with use of reporter gene assays in heterologous tissues, while srSNPs appear more frequent than expected. Computational analysis of mRNA folding indicates that approximately 90% of srSNPs affect mRNA folding, and hence potentially function.

CONCLUSION

Our results indicate that both rSNPs and srSNPs represent a still largely untapped reservoir of variants that contribute to human phenotypic diversity.

Regulatory polymorphism in vitamin K epoxide reductase complex subunit 1 (VKORC1) affects gene expression and warfarin dose requirement

Warfarin dose requirements have been associated with 2 main haplotypes in VKORC1, but the responsible polymorphisms remain unknown. To search for regulatory polymorphisms, we measured allelic mRNA expression of VKORC1 in human liver, heart, and B lymphocytes. The observed 2-fold allelic mRNA expression imbalance narrowed possible candidate SNPs to -1639G>A and 1173C<T. This genotype effect was observed selectively in the liver but not in heart or lymphocytes. In vitro expression of VKORC1 gene constructs, including coding and promoter regions, failed to reveal any genotype effect on transcription and mRNA processing. Chromatin immunoprecipitation with antibodies against acetyl-histone3 and K4-trimethyl-histone3 revealed preferential association of the promoter -1639 G allele with active chromatin, consistent with enhanced mRNA expression. The minor -1639 A allele generates a suppressor E-box binding site, apparently regulating gene expression by a mechanism undetectable with reporter gene assays. A clinical association study demonstrated that promoter SNP -1639G>A, and the tightly linked intron1 SNP 1173C>T, predict warfarin dose more accurately than intron 2 SNP 1542G>C in blacks. Increased warfarin dose requirement in blacks was accounted for by lower frequency of the -1639 A allele. Therefore, -1639G>A is a suitable biomarker for warfarin dosing across ethnic populations.

Gene expression in the etiology of schizophrenia

Gene expression represents a fundamental interface between genes and environment in the development and ongoing plasticity of the human brain. Individual differences in gene expression are likely to underpin much of human diversity, including psychiatric illness. In the past decade, the development of microarray and proteomic technology has enabled global description of gene expression in schizophrenia. However, it is difficult on the basis of gene expression assays alone to distinguish between those changes that constitute primary etiology and those that reflect secondary pathology, compensatory mechanisms, or confounding influences. In this respect, tests of genetic association with schizophrenia will be instructive because changes in gene expression that result from gene variants that are associated with the disorder are likely to be of primary etiological significance. However, regulatory polymorphism is extremely difficult to recognize on the basis of sequence interrogation alone. Functional assays at the messenger RNA and/or protein level will be essential in elucidating the molecular mechanisms underlying genetic association with schizophrenia and are likely to become increasingly important in the identification of regulatory variants with which to test for association with the disorder and related traits. Once established, etiologically relevant changes in gene expression can be recapitulated in model systems in order to elucidate the molecular and physiological pathways that may ultimately give rise to the condition.

Bayesian meta-analysis of genetic association studies with different sets of markers

Robust assessment of genetic effects on quantitative traits or complex-disease risk requires synthesis of evidence from multiple studies. Frequently, studies have genotyped partially overlapping sets of SNPs within a gene or region of interest, hampering attempts to combine all the available data. By using the example of C-reactive protein (CRP) as a quantitative trait, we show how linkage disequilibrium in and around its gene facilitates use of Bayesian hierarchical models to integrate informative data from all available genetic association studies of this trait, irrespective of the SNP typed. A variable selection scheme, followed by contextualization of SNPs exhibiting independent associations within the haplotype structure of the gene, enhanced our ability to infer likely causal variants in this region with population-scale data. This strategy, based on data from a literature based systematic review and substantial new genotyping, facilitated the most comprehensive evaluation to date of the role of variants governing CRP levels, providing important information on the minimal subset of SNPs necessary for comprehensive evaluation of the likely causal relevance of elevated CRP levels for coronary-heart-disease risk by Mendelian randomization. The same method could be applied to evidence synthesis of other quantitative traits, whenever the typed SNPs vary among studies, and to assist fine mapping of causal variants.

Functional significance of lipoprotein lipase HindIII polymorphism associated with the risk of coronary artery disease

Lipoprotein lipase (LPL) plays a pivotal role in lipid metabolism by hydrolyzing triglyceride (TG)-rich lipoprotein particles. Abnormalities in normal LPL function are associated with the risk of coronary artery disease (CAD). A number of genetic variants have been identified in the LPL gene that affects different functions of the LPL protein. A common HindIII polymorphism in intron 8 (T/G) of the LPL gene has been found to be associated with altered plasma TG and HDL-cholesterol, and CAD risk in several studies, but its functional significance is unknown. It has been shown that certain intronic sequence contain regulatory elements that are important for transcription and translational regulation of a gene. In this study we tested the hypothesis that this polymorphism affects the binding site of a transcription factor that regulates the transcription of LPL gene. Electrophoretic mobility shift assays (EMSAs) revealed that the HindIII site binds to a transcription factor and that the mutant allele has lower binding affinity than the wild type allele. Transcription assays containing the entire intron 8 sequence along with full-length human LPL promoter were carried out in COS-1 and human vascular smooth muscle cells. The mutant allele was associated with significantly decreased luciferase expression level compared to the wild type allele in both the muscle (3.394+/-0.022 vs. 4.184+/-0.028; P=4.7 x 10(-6)) and COS-1 (11.603+/-0.409 vs. 14.373+/-1.096; P<0.0001) cells. In conclusion, this study demonstrates for the first time that the polymorphic HindIII site in the LPL gene is functional because it affects the binding of a transcription factor and it also has an impact on LPL expression.

CEBPG regulates ERCC5/XPG expression in human bronchial epithelial cells and this regulation is modified by E2F1/YY1 interactions

Marked inter-individual variation in lung cancer risk cannot be accounted for solely by cigarette smoke and other environmental exposures. Evidence suggests that variation in bronchial epithelial cell expression of key DNA repair genes plays a role. Variation in these genes correlates with variation in expression of CEBPG and E2F1 transcription factors. Here, we investigated the mechanistic basis for correlation of the DNA repair gene ERCC5 (previously known as XPG) with CEBPG and E2F1. CEBPG expression vector transfected into H23 or H460 cell lines up-regulated endogenous ERCC5 and also luciferase from a reporter construct containing 589 bp of ERCC5 5' regulatory region. A recognition site for CEBPG and a region containing sites for YY1 on the sense strand and E2F1 on the anti-sense strand participated in CEBPG up-regulation of ERCC5. CEBPG, E2F1 and YY1 binding to their respective sites were confirmed by electrophoretic mobility shift assay. Thus, we conclude that CEBPG regulates ERCC5 expression and this regulation is modified by E2F1/YY1 interactions. Several polymorphisms have been identified in these regions and, based on the data presented here, it is reasonable to hypothesize that they may contribute to risk for bronchogenic carcinoma.