SNP-based large-scale identification of allele-specific gene expression in human B cells

Polymorphism and variations in gene expression provide the genetic basis for human variation. Allelic variation of gene expression, in particular, may play a crucial role in phenotypic variation and disease susceptibility. To identify genes with allelic expression in human cells, we genotyped genomic DNA and cDNA isolated from 31 immortalized B cell lines from three Centre d'Etude du Polymorphisme Humain (CEPH) families using high-density single-nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs. We identified seven SNPs in five genes with monoallelic expression, 146 SNPs in 125 genes with allelic imbalance in expression with preferentially higher expression of one allele in a heterozygous individual. The monoallelically expressed genes (ERAP2, MDGA1, LOC644422, SDCCAG3P1 and CLTCL1) were regulated by cis-acting, non-imprinted differential allelic control. In addition, all monoallelic gene expression patterns and allelic imbalances in gene expression in B cells were transmitted from parents to offspring in the pedigree, indicating genetic transmission of allelic gene expression. Furthermore, frequent allele substitution, probably due to RNA editing, was also observed in 21 genes in 23 SNPs as well as in 48 SNPs located in regions containing no known genes. In this study, we demonstrated that allelic gene expression is frequently observed in human B cells, and SNP chips are very useful tools for detecting allelic gene expression. Overall, our data provide a valuable framework for better understanding allelic gene expression in human B cells.     

Genome-wide detection of allelic gene expression in hepatocellular carcinoma cells using a human exome SNP chip

Allelic variations in gene expression influence many biological responses and cause phenotypic variations in humans. In this study, Illumina Human Exome BeadChips containing more than 240,000 single nucleotide polymorphisms (SNPs) were used to identify changes in allelic gene expression in hepatocellular carcinoma cells following lipopolysaccharide (LPS) stimulation. We found 17 monoallelically expressed genes, 58 allelic imbalanced genes, and 7 genes showing allele substitution. In addition, we also detected 33 differentially expressed genes following LPS treatment in vitro using these human exome SNP chips. However, alterations in allelic gene expression following LPS treatment were detected in only three genes (MLXIPL, TNC, and MX2), which were observed in one cell line sample only, indicating that changes in allelic gene expression following LPS stimulation of liver cells are rare events. Among a total of 75 genes showing allelic expression in hepatocellular carcinoma cells, either monoallelic or imbalanced, 43 genes (57.33%) had expression quantitative trait loci (eQTL) data, indicating that high-density exome SNP chips are useful and reliable for studying allelic gene expression. Furthermore, most genes showing allelic expression were regulated by cis-acting mechanisms and were also significantly associated with several human diseases. Overall, our study provides a better understanding of allele-specific gene expression in hepatocellular carcinoma cells with and without LPS stimulation and potential clues for the cause of human disease due to alterations in allelic gene expression.     

The impact of cis-acting polymorphisms on the human phenotype

Cis-acting polymorphisms that affect gene expression are now known to be frequent, although the extent and mechanisms by which such variation affects the human phenotype are, as yet, only poorly understood. Key signatures of cis-acting variation are differences in gene expression that are tightly associated with regulatory SNPs or expression Quantitative Trait Loci (eQTL) and an imbalance of allelic expression (AEI) in heterozygous samples. Such cis-acting sequence differences appear often to have been under selection within and between populations and are also thought to be important in speciation. Here we describe the example of lactase persistence. In medical research, variants that affect regulation in cis have been implicated in both monogenic and polygenic disorders, and in the metabolism of drugs. In this review we suggest that by further understanding common regulatory variations and how they interact with other genetic and environmental variables it will be possible to gain insight into important mechanisms behind complex disease, with the potential to lead to new methods of diagnosis and treatments.     

Genetic Variants Contribute to Gene Expression Variability in Humans

Expression quantitative trait loci (eQTL) studies have established convincing relationships between genetic variants and gene expression. Most of these studies focused on the mean of gene expression level, but not the variance of gene expression level (i.e., gene expression variability). In the present study, we systematically explore genome-wide association between genetic variants and gene expression variability in humans. We adapt the double generalized linear model (dglm) to simultaneously fit the means and the variances of gene expression among the three possible genotypes of a biallelic SNP. The genomic loci showing significant association between the variances of gene expression and the genotypes are termed expression variability QTL (evQTL). Using a data set of gene expression in lymphoblastoid cell lines (LCLs) derived from 210 HapMap individuals, we identify cis-acting evQTL involving 218 distinct genes, among which 8 genes, ADCY1, CTNNA2, DAAM2, FERMT2, IL6, PLOD2, SNX7, and TNFRSF11B, are cross-validated using an extra expression data set of the same LCLs. We also identify ∼300 trans-acting evQTL between >13,000 common SNPs and 500 randomly selected representative genes. We employ two distinct scenarios, emphasizing single-SNP and multiple-SNP effects on expression variability, to explain the formation of evQTL. We argue that detecting evQTL may represent a novel method for effectively screening for genetic interactions, especially when the multiple-SNP influence on expression variability is implied. The implication of our results for revealing genetic mechanisms of gene expression variability is discussed.     

cis-Expression QTL analysis of established colorectal cancer risk variants in colon tumors and adjacent normal tissue

Genome-wide association studies (GWAS) have identified 19 risk variants associated with colorectal cancer. As most of these risk variants reside outside the coding regions of genes, we conducted cis-expression quantitative trait loci (cis-eQTL) analyses to investigate possible regulatory functions on the expression of neighboring genes. Forty microsatellite stable and CpG island methylator phenotype-negative colorectal tumors and paired adjacent normal colon tissues were used for genome-wide SNP and gene expression profiling. We found that three risk variants (rs10795668, rs4444235 and rs9929218, using near perfect proxies rs706771, rs11623717 and rs2059252, respectively) were significantly associated (FDR q-value ≤0.05) with expression levels of nearby genes (<2 Mb up- or down-stream). We observed an association between the low colorectal cancer risk allele (A) for rs10795668 at 10p14 and increased expression of ATP5C1 (q = 0.024) and between the colorectal cancer high risk allele (C) for rs4444235 at 14q22.2 and increased expression of DLGAP5 (q = 0.041), both in tumor samples. The colorectal cancer low risk allele (A) for rs9929218 at 16q22.1 was associated with a significant decrease in expression of both NOL3 (q = 0.017) and DDX28 (q = 0.046) in the adjacent normal colon tissue samples. Of the four genes, DLGAP5 and NOL3 have been previously reported to play a role in colon carcinogenesis and ATP5C1 and DDX28 are mitochondrial proteins involved in cellular metabolism and division, respectively. The combination of GWAS findings, prior functional studies, and the cis-eQTL analyses described here suggest putative functional activities for three of the colorectal cancer GWAS identified risk loci as regulating the expression of neighboring genes.     

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.     

Powerful Identification of Cis-regulatory SNPs in Human Primary Monocytes Using Allele-Specific Gene Expression

A large number of genome-wide association studies have been performed during the past five years to identify associations between SNPs and human complex diseases and traits. The assignment of a functional role for the identified disease-associated SNP is not straight-forward. Genome-wide expression quantitative trait locus (eQTL) analysis is frequently used as the initial step to define a function while allele-specific gene expression (ASE) analysis has not yet gained a wide-spread use in disease mapping studies. We compared the power to identify cis-acting regulatory SNPs (cis-rSNPs) by genome-wide allele-specific gene expression (ASE) analysis with that of traditional expression quantitative trait locus (eQTL) mapping. Our study included 395 healthy blood donors for whom global gene expression profiles in circulating monocytes were determined by Illumina BeadArrays. ASE was assessed in a subset of these monocytes from 188 donors by quantitative genotyping of mRNA using a genome-wide panel of SNP markers. The performance of the two methods for detecting cis-rSNPs was evaluated by comparing associations between SNP genotypes and gene expression levels in sample sets of varying size. We found that up to 8-fold more samples are required for eQTL mapping to reach the same statistical power as that obtained by ASE analysis for the same rSNPs. The performance of ASE is insensitive to SNPs with low minor allele frequencies and detects a larger number of significantly associated rSNPs using the same sample size as eQTL mapping. An unequivocal conclusion from our comparison is that ASE analysis is more sensitive for detecting cis-rSNPs than standard eQTL mapping. Our study shows the potential of ASE mapping in tissue samples and primary cells which are difficult to obtain in large numbers.     

Validating Discovered Cis-Acting Regulatory Genetic Variants: Application of an Allele Specific Expression Approach to HapMap Populations

Background

Localising regulatory variants that control gene expression is a challenge for genome research. Several studies have recently identified non-coding polymorphisms associated with inter-individual differences in gene expression. These approaches rely on the identification of signals of association against a background of variation due to other genetic and environmental factors. A complementary approach is to use an Allele-Specific Expression (ASE) assay, which is more robust to the effects of environmental variation and trans-acting genetic factors.

Methodology/Principal Findings

Here we apply an ASE method which utilises heterozygosity within an individual to compare expression of the two alleles of a gene in a single cell. We used individuals from three HapMap population groups and analysed the allelic expression of genes with cis-regulatory regions previously identified using total gene expression studies. We were able to replicate the results in five of the six genes tested, and refined the cis- associated regions to a small number of variants. We also showed that by using multi-populations it is possible to refine the associated cis-effect DNA regions.

Conclusions/Significance

We discuss the efficacy and drawbacks of both total gene expression and ASE approaches in the discovery of cis-acting variants. We show that the ASE approach has significant advantages as it is a cleaner representation of cis-acting effects. We also discuss the implication of using different populations to map cis-acting regions and the importance of finding regulatory variants which contribute to human phenotypic variation.     

Global spread and genetic variants of the two CYP9M10 haplotype forms associated with insecticide resistance in Culex quinquefasciatus Say

Insecticide resistance develops as a genetic factor (allele) conferring lower susceptibility to insecticides proliferates within a target insect population under strong positive selection. Intriguingly, a resistance allele pre-existing in a population often bears a series of further adaptive allelic variants through new mutations. This phenomenon occasionally results in replacement of the predominating resistance allele by fitter new derivatives, and consequently, development of greater resistance at the population level. The overexpression of the cytochrome P450 gene CYP9M10 is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. Previously, we have found two genealogically related overexpressing CYP9M10 haplotypes, which differ in gene copy number (duplicated and non-duplicated). The duplicated haplotype was derived from the non-duplicated overproducer probably recently. In the present study, we investigated allelic series of CYP9M10 involved in three C. quinquefasciatus laboratory colonies recently collected from three different localities. Duplicated and non-duplicated overproducing haplotypes coexisted in African and Asian colonies indicating a global distribution of both haplotype lineages. The duplicated haplotypes both in the Asian and African colonies were associated with higher expression levels and stronger resistance than non-duplicated overproducing haplotypes. There were slight variation in expression level among the non-duplicated overproducing haplotypes. The nucleotide sequences in coding and upstream regions among members of this group also showed a little diversity. Non-duplicated overproducing haplotypes with relatively higher expression were genealogically closer to the duplicated haplotypes than the other non-duplicated overproducing haplotypes, suggesting multiple cis-acting mutations before duplication.     

Asynchronous Replication,Mono-Allelic Expression,and Long Range Cis-Effects of ASAR6

Mammalian chromosomes initiate DNA replication at multiple sites along their length during each S phase following a temporal replication program. The majority of genes on homologous chromosomes replicate synchronously. However, mono-allelically expressed genes such as imprinted genes, allelically excluded genes, and genes on female X chromosomes replicate asynchronously. We have identified a cis-acting locus on human chromosome 6 that controls this replication-timing program. This locus encodes a large intergenic non-coding RNA gene named Asynchronous replication and Autosomal RNA on chromosome 6, or ASAR6. Disruption of ASAR6 results in delayed replication, delayed mitotic chromosome condensation, and activation of the previously silent alleles of mono-allelic genes on chromosome 6. The ASAR6 gene resides within an ∼1.2 megabase domain of asynchronously replicating DNA that is coordinated with other random asynchronously replicating loci along chromosome 6. In contrast to other nearby mono-allelic genes, ASAR6 RNA is expressed from the later-replicating allele. ASAR6 RNA is synthesized by RNA Polymerase II, is not polyadenlyated, is restricted to the nucleus, and is subject to random mono-allelic expression. Disruption of ASAR6 leads to the formation of bridged chromosomes, micronuclei, and structural instability of chromosome 6. Finally, ectopic integration of cloned genomic DNA containing ASAR6 causes delayed replication of entire mouse chromosomes.     

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.     

Evaluation of Parkinson Disease Risk Variants as Expression-QTLs

The recent Parkinson Disease GWAS Consortium meta-analysis and replication study reports association at several previously confirmed risk loci SNCA, MAPT, GAK/DGKQ, and HLA and identified a novel risk locus at RIT2. To further explore functional consequences of these associations, we investigated modification of gene expression in prefrontal cortex brain samples of pathologically confirmed PD cases (N = 26) and controls (N = 24) by 67 associated SNPs in these 5 loci. Association between the eSNPs and expression was evaluated using a 2-degrees of freedom test of both association and difference in association between cases and controls, adjusted for relevant covariates. SNPs at each of the 5 loci were tested for cis-acting effects on all probes within 250 kb of each locus. Trans-effects of the SNPs on the 39,122 probes passing all QC on the microarray were also examined. From the analysis of cis-acting SNP effects, several SNPs in the MAPT region show significant association to multiple nearby probes, including two strongly correlated probes targeting the gene LOC644246 and the duplicated genes LRRC37A and LRRC37A2, and a third uncorrelated probe targeting the gene DCAKD. Significant cis-associations were also observed between SNPs and two probes targeting genes in the HLA region on chromosome 6. Expanding the association study to examine trans effects revealed an additional 23 SNP-probe associations reaching statistical significance (p<2.8×10⁻⁸) including SNPs from the SNCA, MAPT and RIT2 regions. These findings provide additional context for the interpretation of PD associated SNPs identified in recent GWAS as well as potential insight into the mechanisms underlying the observed SNP associations.