Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer 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.     

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.     

Imprinting analysis of the mouse chromosome 7C region in DNMT1-null embryos

The mouse chromosome 7C, orthologous to the human 15q11–q13 has an imprinted domain, where most of the genes are expressed only from the paternal allele. The imprinted domain contains paternally expressed genes, Snurf/Snrpn, Ndn, Magel2, Mkrn3, and Frat3, C/D-box small nucleolar RNAs (snoRNAs), and the maternally expressed gene, Ube3a. Imprinted expression in this large (approximately 3–4 Mb) domain is coordinated by a bipartite cis-acting imprinting center (IC), located upstream of the Snurf/Snrpn gene. The molecular mechanism how IC regulates gene expression of the whole domain remains partially understood. Here we analyzed the relationship between imprinted gene expression and DNA methylation in the mouse chromosome 7C using DNA methyltransferase 1 (DNMT1)-null mutant embryos carrying Dnmt1^ps alleles, which show global loss of DNA methylation and embryonic lethality. In the DNMT1-null embryos at embryonic day 9.5, the paternally expressed genes were biallelically expressed. Bisulfite DNA methylation analysis revealed loss of methylation on the maternal allele in the promoter regions of the genes. These results demonstrate that DNMT1 is necessary for monoallelic expression of the imprinted genes in the chromosome 7C domain, suggesting that DNA methylation in the secondary differentially methylated regions (DMRs), which are acquired during development serves primarily to control the imprinted expression from the maternal allele in the mouse chromosome 7C.     

A catalogue of validated single nucleotide polymorphisms in bovine orthologs of mammalian imprinted genes and associations with beef production traits

Genetic (or ‘genomic’) imprinting, a feature of approximately 100 mammalian genes, results in monoallelic expression from one of the two parentally inherited chromosomes. To date, most studies have been directed on imprinted genes in murine or human models; however, there is burgeoning interest in the effects of imprinted genes in domestic livestock species. In particular, attention has focused on imprinted genes that influence foetal growth and development and that are associated with several economically important production traits in cattle, sheep and pigs. We have re-sequenced regions in 20 candidate bovine imprinted genes in order to validate single nucleotide polymorphisms (SNPs) that may influence important production traits in cattle. Putative SNPs detected via re-sequencing were subsequently re-formatted for high-throughput SNP genotyping in 185 cattle samples comprising 138 performance-tested European Bos taurus (all Limousin bulls), 29 African B. taurus and 18 Indian B. indicus samples. Analysis of the resulting genotypic data identified 117 validated SNPs. Preliminary genotype–phenotype association analyses using 83 SNPs that were polymorphic in the Limousin samples with minor allele frequencies ⩾0.05 revealed significant associations between two candidate bovine imprinted genes and a range of important beef production traits: average daily gain, average feed intake, live weight, feed conversion ratio, residual feed intake and residual gain. These genes were the Ras protein-specific guanine nucleotide releasing factor gene (RASGRF1) and the zinc finger, imprinted 2 gene (ZIM2). Despite the relatively small sample size used in these analyses, the observed associations with production traits are supported by the purported biological function of the RASGRF1 and ZIM2 gene products. These results support the hypothesis that imprinted genes contribute significantly to important complex production traits in cattle. Furthermore, these SNPs may be usefully incorporated into future marker-assisted and genomic selection breeding schemes.     

Sequence variation in 3′UTR region of crtRB1 gene and its effect on β-carotene accumulation in maize kernel

β-carotene fortification of maize has emerged as a potential, long-term and sustainable approach to alleviate vitamin A deficiency in humans. Among the several genes involved in the carotenoid biosynthetic pathway, the 543 bp allele at crtRB1 3′TE (Transposable Element) gene (allele 1, without insertion) is associated with higher β-carotene accumulation. Estimation of β-carotene through high performance liquid chromatography showed that the CIMMYT genotypes with allele 1 had high kernel β-carotene content whereas the Indian inbreds with the same allele had low β-carotene content. To know the reason for this variation, allele 1 of crtRB1 3′TE gene was sequenced from a set of 11 diverse maize inbreds collected from CIMMYT and Indian germplasm. The sequence data of the allele 1 revealed the presence of 13 single nucleotide polymorphisms (SNPs) and 7 insertions and deletions (InDels). Exonic region had two SNPs, intronic region had one SNP and one InDel, whereas 3′-untranslated region (UTR) region of the gene showed 10 SNPs and 6 InDels. Among the several SNPs and InDels, SNP4, SNP13, InDel6 and InDel7 identified in the 3′-UTR region clearly differentiated the high and the low β-carotene genotypes. These 3′-UTR polymorphisms in allele 1 of the crtRB1 3′TE gene could be associated with the variation in kernel β-carotene accumulation by regulating the translation and stability of the mRNA. The SNPs and the InDels associated with higher level of β-carotene will be used as a gene-based marker(s) in selection of genotypes and to develop biofortified maize hybrids to alleviate vitamin A deficiency in humans.     

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.     

Identification of functional DNA variants in the constitutive promoter region of MDM2

Although mutations in the oncoprotein murine double minute 2 (MDM2) are rare, MDM2 gene overexpression has been observed in several human tumors. Given that even modest changes in MDM2 levels might influence the p53 tumor suppressor signaling pathway, we postulated that sequence variation in the promoter region of MDM2 could lead to disregulated expression and variation in gene dosage. Two promoters have been reported for MDM2; an internal promoter (P2), which is located near the end of intron 1 and is p53-responsive, and an upstream constitutive promoter (P1), which is p53-independent. Both promoter regions contain DNA variants that could influence the expression levels of MDM2, including the well-studied single nucleotide polymorphism (SNP) SNP309, which is located in the promoter P2; i.e., upstream of exon 2. In this report, we screened the promoter P1 for DNA variants and assessed the functional impact of the corresponding SNPs. Using the dbSNP database and genotyping validation in individuals of European descent, we identified three common SNPs (?1494?G?>?A; indel 40?bp; and ?182?C?>?G). Three major promoter haplotypes were inferred by using these three promoter SNPs together with rs2279744 (SNP309). Following subcloning into a gene reporter system, we found that two of the haplotypes significantly influenced MDM2 promoter activity in a haplotype-specific manner. Site-directed mutagenesis experiments indicated that the 40?bp insertion/deletion variation is causing the observed allelic promoter activity. This study suggests that part of the variability in the MDM2 expression levels could be explained by allelic p53-independent P1 promoter activity.     

Comparison of multimarker logistic regression models, with application to a genomewide scan of schizophrenia

Background

Imprinted genes show expression from one parental allele only and are important for development and behaviour. This extreme mode of allelic imbalance has been described for approximately 56 human genes. Imprinting status is often disrupted in cancer and dysmorphic syndromes. More subtle variation of gene expression, that is not parent-of-origin specific, termed 'allele-specific gene expression' (ASE) is more common and may give rise to milder phenotypic differences. Using two allele-specific high-throughput technologies alongside bioinformatics predictions, normal term human placenta was screened to find new imprinted genes and to ascertain the extent of ASE in this tissue.

Results

Twenty-three family trios of placental cDNA, placental genomic DNA (gDNA) and gDNA from both parents were tested for 130 candidate genes with the Sequenom MassArray system. Six genes were found differentially expressed but none imprinted. The Illumina ASE BeadArray platform was then used to test 1536 SNPs in 932 genes. The array was enriched for the human orthologues of 124 mouse candidate genes from bioinformatics predictions and 10 human candidate imprinted genes from EST database mining. After quality control pruning, a total of 261 informative SNPs (214 genes) remained for analysis. Imprinting with maternal expression was demonstrated for the lymphocyte imprinted gene ZNF331 in human placenta. Two potential differentially methylated regions (DMRs) were found in the vicinity of ZNF331. None of the bioinformatically predicted candidates tested showed imprinting except for a skewed allelic expression in a parent-specific manner observed for PHACTR2, a neighbour of the imprinted PLAGL1 gene. ASE was detected for two or more individuals in 39 candidate genes (18%).

Conclusions

Both Sequenom and Illumina assays were sensitive enough to study imprinting and strong allelic bias. Previous bioinformatics approaches were not predictive of new imprinted genes in the human term placenta. ZNF331 is imprinted in human term placenta and might be a new ubiquitously imprinted gene, part of a primate-specific locus. Demonstration of partial imprinting of PHACTR2 calls for re-evaluation of the allelic pattern of expression for the PHACTR2-PLAGL1 locus. ASE was common in human term placenta.     

Impact of cytokine and cytokine receptor gene polymorphisms on cellular immunity after smallpox vaccination

We explored associations between SNPs in cytokine/cytokine receptor genes and cellular immunity in subjects following primary smallpox vaccination. We also analyzed the genotype–phenotype associations discovered in the Caucasian subjects among a cohort of African-Americans. In Caucasians we found 277 associations (p < 0.05) between gene SNPs and inter-individual variations in IFN-α, IL-12p40, IL-1β, IL-2, and TNF-α secretion levels. A collection of SNPs in the IL1RN, IL2RB, IL4R, IL6, IL10RB, IL12A, and IL12RB2 genes had consistent associations among both Caucasians and African-Americans. A regulatory SNP (rs452204) in the IL1RN gene was significantly associated with higher levels of IL-2 secretion in an allele dose-dependent manner in both race groups (p = 0.05 for Caucasians and p = 0.002 for African-Americans). IL12RB2 polymorphism rs3790567 was associated with a dose-related decrease in IL-1β secretion (p = 0.009 for Caucasians and p = 0.01 for African-Americans). Our results demonstrate that variations in smallpox vaccine-induced cytokine responses are modulated by genetic polymorphisms in cytokine and cytokine receptor genes.     

Use of Allele-Specific FAIRE to Determine Functional Regulatory Polymorphism Using Large-Scale Genotyping Arrays

Following the widespread use of genome-wide association studies (GWAS), focus is turning towards identification of causal variants rather than simply genetic markers of diseases and traits. As a step towards a high-throughput method to identify genome-wide, non-coding, functional regulatory variants, we describe the technique of allele-specific FAIRE, utilising large-scale genotyping technology (FAIRE-gen) to determine allelic effects on chromatin accessibility and regulatory potential. FAIRE-gen was explored using lymphoblastoid cells and the 50,000 SNP Illumina CVD BeadChip. The technique identified an allele-specific regulatory polymorphism within NR1H3 (coding for LXR-α), rs7120118, coinciding with a previously GWAS-identified SNP for HDL-C levels. This finding was confirmed using FAIRE-gen with the 200,000 SNP Illumina Metabochip and verified with the established method of TaqMan allelic discrimination. Examination of this SNP in two prospective Caucasian cohorts comprising 15,000 individuals confirmed the association with HDL-C levels (combined beta = 0.016; p = 0.0006), and analysis of gene expression identified an allelic association with LXR-α expression in heart tissue. Using increasingly comprehensive genotyping chips and distinct tissues for examination, FAIRE-gen has the potential to aid the identification of many causal SNPs associated with disease from GWAS.     

Harnessing expression data to identify novel candidate genes in polycystic ovary syndrome

Novel pathways in polycystic ovary syndrome (PCOS) are being identified in gene expression studies in PCOS tissues; such pathways may contain key genes in disease etiology. Previous expression studies identified both dickkopf homolog 1 (DKK1) and DnaJ (Hsp40) homolog, subfamily B, member 1 (DNAJB1) as differentially expressed in PCOS tissue, implicating them as candidates for PCOS susceptibility. To test this, we genotyped a discovery cohort of 335 PCOS cases and 198 healthy controls for three DKK1 single nucleotide polymorphisms (SNPs) and four DNAJB1 SNPs and a replication cohort of 396 PCOS cases and 306 healthy controls for 1 DKK1 SNP and 1 DNAJB1 SNP. SNPs and haplotypes were determined and tested for association with PCOS and component phenotypes. We found that no single nucleotide polymorphisms were associated with PCOS risk; however, the major allele of rs1569198 from DKK1 was associated with increased total testosterone (discovery cohort P = 0.0035) and dehydroepiandrosterone sulfate (replication cohort P = 0.05). Minor allele carriers at rs3962158 from DNAJB1 had increased fasting insulin (discovery cohort P = 0.003), increased HOMA-IR (discovery cohort P = 0.006; replication cohort P = 0.036), and increased HOMA-%B (discovery cohort P = 0.004). Carriers of haplotype 2 at DNAJB1 also had increased fasting insulin, HOMA-IR, and HOMA-%B. These findings suggest that genetic variation in DKK1 and DNAJB1 may have a role in the hyperandrogenic and metabolic dysfunction of PCOS, respectively. Our results also demonstrate the utility of gene expression data as a source of novel candidate genes in PCOS, a complex and still incompletely defined disease, for which alternative methods of gene identification are needed.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.