Investigating cis-acting regulatory variation using assays of relative allelic expression

Identification of rs11615992 as a novel regulatory SNP for human P2RX7 by allele-specific expression

P2RX7 (purinergic receptor P2X 7) is an important membrane ion channel and involved in multiple physiological processes. One non-synonymous SNP on P2RX7, rs3751143, had been proven to reduce ion channel function and further associated with multiple diseases. However, it was still unclear whether there were other cis-regulatory elements for P2RX7, which might further contribute to related diseases. Allele-specific expression (ASE) is a robust and sensitive approach to identify the potential functional region in human genome. In the current study, we measured ASE on rs3751143 in lung tissues and observed a consistent excess of A allele over C (P = 0.001), which indicated that SNP(s) in linkage disequilibrium (LD) could regulate P2RX7 expression. By analyzing the 1000 genomes project data for Chinese, one SNP locating ~ 5 kb away and downstream of P2RX7, rs11615992, was disclosed to be in strong LD with rs3751143. The dual-luciferase assay confirmed that rs11615992 could alter target gene expression in lung cell line. Through chromosome conformation capture, it was verified that the region surrounding rs11615992 could interact with P2RX7 promoter and effect as an enhancer. By chromatin immunoprecipitation, the related transcription factor POU2F1 (POU class 2 homeobox 1) was recognized to bind the region spanning rs11615992. Our work identified a novel long-distance cis-regulatory SNP for P2RX7, which might contribute to multiple diseases.

Allele-specific expression identified rs2509956 as a novel long-distance cis-regulatory SNP for SCGB1A1, an important gene for multiple pulmonary diseases

SCGB1A1 (secretoglobin family 1A member 1) is an important protein for multiple pulmonary diseases, especially asthma, chronic obstructive pulmonary disease, and lung cancer. One single-nucleotide polymorphism (SNP) at 5'-untranslated region of SCGB1A1, rs3741240, has been suggested to be associated with reduced protein expression and further asthma susceptibility. However, it was still unclear whether there were other cis-regulatory elements for SCGB1A1 that might further contribute to pulmonary diseases. Allele-specific expression (ASE) is a novel approach to identify the functional region in human genome. In the present study, we measured ASE on rs3741240 in lung tissues and observed a consistent excess of G allele over A (P < 10^-6), which indicated that this SNP or the one(s) in linkage disequilibrium (LD) could regulate SCGB1A1 expression. By analyzing 1000 Genomes Project data for Chinese, one SNP locating ~10.2 kb away and downstream of SCGB1A1, rs2509956, was identified to be in strong LD with rs3741240. Reporter gene assay confirmed that both SNPs could regulate gene expression in the lung cell. By chromosome conformation capture, it was verified that the region surrounding rs2509956 could interact with SCGB1A1 promoter region and act as an enhancer. Through chromatin immunoprecipitation and overexpression assay, the related transcription factor RELA (RELA proto-oncogene, NF-kB subunit) was recognized to bind the region spanning rs2509956. Our work identified a novel long-distance cis-regulatory SNP for SCGB1A1, which might contribute to multiple pulmonary diseases.

Evidence of allelic imbalance in the schizophrenia susceptibility gene ZNF804A in human dorsolateral prefrontal cortex

The rs1344706, an intronic SNP within the zinc-finger protein 804A gene (ZNF804A), was identified as one of the most compelling risk SNPs for schizophrenia (SZ) and bipolar disorder (BD). It is however not clear by which molecular mechanisms ZNF804A increases disease risk. We evaluated the role of ZNF804A in SZ and BD by genotyping the originally associated rs1344706 SNP and an exonic SNP (rs12476147) located in exon four of ZNF804A in a sample of 422 SZ, 382 BD, and 507 controls from the isolated population of the Costa Rica Central Valley. We also investigated the rs1344706 SNP for allelic specific expression (ASE) imbalance in the dorsolateral prefrontal cortex (DLPFC) of 46 heterozygous postmortem brains. While no significant association between rs1344706 and SZ or BD was observed in the Costa Rica sample, we observed an increased risk of SZ for the minor allele (A) of the exonic rs12476147 SNP (p=0.026). Our ASE assay detected a significant over-expression of the rs12476147 A allele in DLPFC of rs1344706 heterozygous subjects. Interestingly, cDNA allele ratios were significantly different according to the intronic rs1344706 genotypes (p-value=0.03), with the rs1344706 A allele associated with increased ZNF804A rs12476147 A allele expression (average 1.06, p-value=0.02, for heterozygous subjects vs. genomic DNA). In conclusion, we have demonstrated a significant association of rs12476147 with SZ, and using a powerful within-subject design, an allelic expression imbalance of ZNF804A exonic SNP rs12476147 in the DLPFC. Although this data does not preclude the possibility of other functional variants in ZNF804A, it provides evidence that the rs1344706 SZ risk allele is the cis-regulatory variant directly responsible for this allelic expression imbalance in adult cortex.

Germline variants in MRE11/RAD50/NBN complex genes in childhood leukemia

BACKGROUND

The MRE11, RAD50, and NBN genes encode proteins of the MRE11-RAD50-NBN (MRN) complex involved in cellular response to DNA damage and the maintenance of genome stability. In our previous study we showed that the germline p.I171V mutation in NBN may be considered as a risk factor in the development of childhood acute lymphoblastic leukemia (ALL) and some specific haplotypes of that gene may be associated with childhood leukemia. These findings raise important questions about the role of mutations in others genes of the MRN complex in childhood leukemia. The aim of this study was to answer the question whether MRE11 and RAD50 alterations may be associated with childhood ALL or AML.

METHODS

We estimated the frequency of constitutional mutations and polymorphisms in selected regions of MRE11, RAD50, and NBN in the group of 220 children diagnosed with childhood leukemias and controls (n=504/2200). The analysis was performed by specific amplification of region of interest by PCR and followed by multi-temperature single-strand conformation polymorphism (PCR-MSSCP) technique. We performed two molecular tests to examine any potential function of the detected the c.551+19G>A SNP in RAD50 gene. To our knowledge, this is the first analysis of the MRE11, RAD50 and NBN genes in childhood leukemia.

RESULTS

The frequency of either the AA genotype or A allele of RAD50_rs17166050 were significantly different in controls compared to leukemia group (ALL+AML) (p<0.0019 and p<0.0019, respectively). The cDNA analysis of AA or GA genotypes carriers has not revealed evidence of splicing abnormality of RAD50 pre-mRNA. We measured the allelic-specific expression of G and A alleles at c.551+19G>A and the statistically significant overexpression of the G allele has been observed. Additionally we confirmed the higher incidence of the p.I171V mutation in the leukemia group (7/220) than among controls (12/2400) (p<0.0001).

CONCLUSION

The formerly reported sequence variants in the RAD50 and MRE11 gene may not constitute a risk factor of childhood ALL in Polish population. The RAD50_rs17166050 variant allele is linked to decreased ALL risk (p<0.0009, OR=0.6358 (95%CI: 0.4854-0.8327)). Despite the fact that there is no splicing abnormality in carriers of the variant allele but an excess of the G over the A allele was consistently observed. This data demonstrate that some specific alternations of the RAD50 gene may be associated with childhood ALL.

Next-generation DNA sequencing-based assay for measuring allelic expression imbalance (AEI) of candidate neuropsychiatric disorder genes in human brain

BACKGROUND

Common genetic variants that regulate gene expression are widely suspected to contribute to the etiology and phenotypic variability of complex diseases. Although high-throughput, microarray-based assays have been developed to measure differences in mRNA expression among independent samples, these assays often lack the sensitivity to detect rare mRNAs and the reproducibility to quantify small changes in mRNA expression. By contrast, PCR-based allelic expression imbalance (AEI) assays, which use a "marker" single nucleotide polymorphism (mSNP) in the mRNA to distinguish expression from pairs of genetic alleles in individual samples, have high sensitivity and accuracy, allowing differences in mRNA expression greater than 1.2-fold to be quantified with high reproducibility. In this paper, we describe the use of an efficient PCR/next-generation DNA sequencing-based assay to analyze allele-specific differences in mRNA expression for candidate neuropsychiatric disorder genes in human brain.

RESULTS

Using our assay, we successfully analyzed AEI for 70 candidate neuropsychiatric disorder genes in 52 independent human brain samples. Among these genes, 62/70 (89%) showed AEI ratios greater than 1 ± 0.2 in at least one sample and 8/70 (11%) showed no AEI. Arranging log2AEI ratios in increasing order from negative-to-positive values revealed highly reproducible distributions of log2AEI ratios that are distinct for each gene/marker SNP combination. Mathematical modeling suggests that these log2AEI distributions can provide important clues concerning the number, location and contributions of cis-acting regulatory variants to mRNA expression.

CONCLUSIONS

We have developed a highly sensitive and reproducible method for quantifying AEI of mRNA expressed in human brain. Importantly, this assay allowed quantification of differential mRNA expression for many candidate disease genes entirely missed in previously published microarray-based studies of mRNA expression in human brain. Given the ability of next-generation sequencing technology to generate large numbers of independent sequencing reads, our method should be suitable for analyzing from 100- to 200-candidate genes in 100 samples in a single experiment. We believe that this is the appropriate scale for investigating variation in mRNA expression for defined sets candidate disorder genes, allowing, for example, comprehensive coverage of genes that function within biological pathways implicated in specific disorders. The combination of AEI measurements and mathematical modeling described in this study can assist in identifying SNPs that correlate with mRNA expression. Alleles of these SNPs (individually or as sets) that accurately predict high- or low-mRNA expression should be useful as markers in genetic association studies aimed at linking candidate genes to specific neuropsychiatric disorders.

Differential allelic expression of c.1568C > A at UGT2B15 is due to variation in a novel cis-regulatory element in the 3'UTR

Differential allelic expression (DAE) is a powerful tool to identify cis-regulatory elements for gene expression. The UDP-glucuronosyltransferase 2 family, polypeptide B15 (UGT2B15), is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In the present study, we measured the relative expression of two alleles at SNP c.1568C>A (rs4148269) in this gene, which causes an amino acid substitution (T523K). An excess of the C over the A allele was consistently observed in both liver (P=0.0021) and breast (P=0.012) samples, suggesting that SNP(s) in strong linkage disequilibrium (LD) with c.1568C>A can regulate UGT2B15 expression in both tissues. By resequencing, one such SNP, c.1761T>C (rs3100) in 3' untranslated region (UTR), was identified. Reporter gene assays showed that the 1761T allele results in a significantly higher gene expression level than the 1761C allele in HepG2, MCF-7, LNCaP, and Caco-2 cell lines (all P<0.001), thus indicating that this variation can regulate UGT2B15 gene expression in liver, breast, colon, and prostate tissues. Considering its location, we postulated that this SNP is within an unknown microRNA binding site and can influence microRNA targeting. Considering the importance of UGT2B15 in metabolism, we proposed that this SNP might contribute to multiple cancer risk and variability in drug response.

Study of transcriptional effects in Cis at the IFIH1 locus

BACKGROUND

The Thr allele at the non-synonymous single-nucleotide polymorphism (nsSNP) Thr946Ala in the IFIH1 gene confers risk for Type 1 diabetes (T1D). The SNP is embedded in a 236 kb linkage disequilibrium (LD) block that includes four genes: IFIH1, GCA, FAP and KCNH7. The absence of common nsSNPs in the other genes makes the IFIH1 SNP the strongest functional candidate, but it could be merely a marker of association, due to LD with a variant regulating expression levels of IFIH1 or neighboring genes.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the effect of the T1D-associated variation on mRNA transcript expression of these genes. Heterozygous mRNA from lymphoblastoid cell lines (LCLs), pancreas and thymus was examined by allelic expression imbalance, to detect effects in cis on mRNA expression. Using single-nucleotide primer extension, we found no difference between mRNA transcripts in 9 LCLs, 6 pancreas and 13 thymus samples, suggesting that GCA and FAP are not involved. On the other hand, KCNH7 was not expressed at a detectable level in all tissues examined. Moreover, the association of the Thr946Ala SNP with T1D is not due to modulation of IFIH1 expression in organs involved in the disease, pointing to the IFIH1 nsSNP as the causal variant.

CONCLUSIONS/SIGNIFICANCE

The mechanism of the association of the nsSNP with T1D remains to be determined, but does not involve mRNA modulation. It becomes necessary to study differential function of the IFIH1 protein alleles at Thr946Ala to confirm that it is responsible for the disease association.

Allelic imbalance (AI) identifies novel tissue-specific cis-regulatory variation for human UGT2B15

Allelic imbalance (AI) is a powerful tool to identify cis-regulatory variation for gene expression. UGT2B15 is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In this study, we measured the relative expression of two alleles at this gene by using SNP rs1902023:G>T. An excess of the G over the T allele was consistently observed in liver (P<0.001), but not in breast (P=0.06) samples, suggesting that SNPs in strong linkage disequilibrium with G253T regulate UGT2B15 expression in liver. Seven such SNPs were identified by resequencing the promoter and exon 1, which define two distinct haplotypes. Reporter gene assays confirmed that one haplotype displayed approximately 20% higher promoter activity compared to the other major haplotype in liver HepG2 (P<0.001), but not in breast MCF-7 (P=0.540) cells. Reporter gene assays with additional constructs pointed to rs34010522:G>T and rs35513228:C>T as the cis-regulatory variants; both SNPs were also evaluated in LNCaP and Caco-2 cells. By ChIP, we showed that the transcription factor Nrf2 binds to the region spanning rs34010522:G>T in all four cell lines. Our results provide a good example for how AI can be used to identify cis-regulatory variation and gain insights into the tissue specific regulation of gene expression.

A functional haplotype implicated in vulnerability to develop cocaine dependence is associated with reduced PDYN expression in human brain

Dynorphin peptides and the kappa-opioid receptor are important in the rewarding properties of cocaine, heroin, and alcohol. We tested polymorphisms of the prodynorphin gene (PDYN) for association with cocaine dependence and cocaine/alcohol codependence. We genotyped six single nucleotide polymorphisms (SNPs), located in the promoter region, exon 4 coding, and 3' untranslated region, in 106 Caucasians and 204 African Americans who were cocaine dependent, cocaine/alcohol codependent, or controls. In Caucasians, we found point-wise significant associations of 3'UTR SNPs (rs910080, rs910079, and rs2235749) with cocaine dependence and cocaine/alcohol codependence. These SNPs are in high linkage disequilibrium, comprising a haplotype block. The haplotype CCT was significantly experiment-wise associated with cocaine dependence and with combined cocaine dependence and cocaine/alcohol codependence (false discovery rate, q=0.04 and 0.03, respectively). We investigated allele-specific gene expression of PDYN, using SNP rs910079 as a reporter, in postmortem human brains from eight heterozygous subjects, using SNaPshot assay. There was significantly lower expression for C allele (rs910079), with ratios ranging from 0.48 to 0.78, indicating lower expression of the CCT haplotype of PDYN in both the caudate and nucleus accumbens. Analysis of total PDYN expression in 43 postmortem brains also showed significantly lower levels of preprodynorphin mRNA in subjects having the risk CCT haplotype. This study provides evidence that a 3'UTR PDYN haplotype, implicated in vulnerability to develop cocaine addiction and/or cocaine/alcohol codependence, is related to lower mRNA expression of the PDYN gene in human dorsal and ventral striatum.

Cis-regulatory variation is typically polyallelic in Drosophila

Gene expression levels vary heritably, with approximately 25-35% of the loci affecting expression acting in cis. We characterized standing cis-regulatory variation among 16 wild-derived strains of Drosophila melanogaster. Our experiment's robust biological and technical replication enabled precise estimates of variation in allelic expression on a high-throughput SNP genotyping platform. We observed concordant, significant differential allelic expression (DAE) in 7/10 genes queried with multiple SNPs, and every member of a set of eight additional, one-assay genes suggest significant DAE. Four of the high-confidence, multiple-assay genes harbor three or more statistically distinguishable allelic classes, often at intermediate frequency. Numerous intermediate-frequency, detectable regulatory polymorphisms cast doubt on a model in which cis-acting variation is a product of deleterious mutations of large effect. Comparing our data to predictions of population genetics theory using coalescent simulations, we estimate that a typical gene harbors 7-15 cis-regulatory sites (nucleotides) at which a selectively neutral mutation would elicit an observable expression phenotype. If standing cis-regulatory variation is actually slightly deleterious, the true mutational target size is larger.

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.

Expression of human PTPN22 alleles

Considering the female predominance in most of the autoimmune disorders that associate with the PTPN22 Trp620 variant and the complexity by which this variant influences immunologic tolerance, the objective of this study was to ascertain if the allele-specific expression of the disease-associated Arg620Trp polymorphism is affected by cis-acting or sex-specific trans-acting factor/s (e.g. sex-hormones). The use of the allele-specific transcript quantification of the Arg620Trp encoding 1858T polymorphism revealed no difference in the expression of the 1858C- and T-alleles in non-stimulated peripheral blood mononuclear cells (PBMCs) from non-pregnant female subjects, male subjects or pregnant female subjects in first or third trimester (P=0.70), respectively. While the transcription of PTPN22 in anti-CD3/anti-CD28 stimulated PBMCs increased fourfold (P<0.0001) and 13-fold (P<0.0001) after 48 and 72 h of activation, respectively, the expression of PTPN22 1858C- and T-alleles increased to the same extent (P=0.64). The present result essentially excludes such phenomena as a partial explanation for the female predominance in most of the autoimmune disorders that associate with the PTPN22 Trp620 variant.