Cis-regulatory Single Nucleotide Polymorphisms Research Articles

Transcription factor-binding k-mer analysis clarifies the cell type dependency of binding specificities and cis-regulatory SNPs in humans

BackgroundTranscription factors (TFs) exhibit heterogeneous DNA-binding specificities in individual cells and whole organisms under natural conditions, and de novo motif discovery usually provides multiple motifs, even from a single chromatin immunoprecipitation-sequencing (ChIP-seq) sample. Despite the accumulation of ChIP-seq data and ChIP-seq-derived motifs, the diversity of DNA-binding specificities across different TFs and cell types remains largely unexplored.ResultsHere, we applied MOCCS2, our k-mer-based motif discovery method, to a collection of human TF ChIP-seq samples across diverse TFs and cell types, and systematically computed profiles of TF-binding specificity scores for all k-mers. After quality control, we compiled a set of TF-binding specificity score profiles for 2,976 high-quality ChIP-seq samples, comprising 473 TFs and 398 cell types. Using these high-quality samples, we confirmed that the k-mer-based TF-binding specificity profiles reflected TF- or TF-family dependent DNA-binding specificities. We then compared the binding specificity scores of ChIP-seq samples with the same TFs but with different cell type classes and found that half of the analyzed TFs exhibited differences in DNA-binding specificities across cell type classes. Additionally, we devised a method to detect differentially bound k-mers between two ChIP-seq samples and detected k-mers exhibiting statistically significant differences in binding specificity scores. Moreover, we demonstrated that differences in the binding specificity scores between k-mers on the reference and alternative alleles could be used to predict the effect of variants on TF binding, as validated by in vitro and in vivo assay datasets. Finally, we demonstrated that binding specificity score differences can be used to interpret disease-associated non-coding single-nucleotide polymorphisms (SNPs) as TF-affecting SNPs and provide candidates responsible for TFs and cell types.ConclusionsOur study provides a basis for investigating the regulation of gene expression in a TF-, TF family-, or cell-type-dependent manner. Furthermore, our differential analysis of binding-specificity scores highlights noncoding disease-associated variants in humans.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle.

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = -3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Abstract 2905: Induction of hTERT and increased proliferative potential in conditionally reprogrammed normal bronchial epithelial cells

Abstract Background. Based on increasing evidence from this laboratory and genome wide association studies (GWAS), single nucleotide polymorphisms (SNPs) responsible for inter-individual variation in normal bronchial epithelial cell (NBEC) cis-regulation of antioxidant, DNA repair, and cell cycle control genes are key determinants of lung cancer risk. Thus, there is a need for NBEC culture methods that enable extended population doublings without genetic alteration to enable experimental investigation of putative cis-regulatory SNPs in NBEC. Toward this goal, we assessed the effect of previously reported conditional reprogrammed culture (CRC) conditions on regulation of human telomerase reverse transcriptase (hTERT) transcript abundance and proliferative potential in NBEC. Methods. NBEC were obtained by bronchoscopic brush from eight individuals after obtaining informed consent to an IRB-approved protocol. NBEC were incubated in three different culture conditions: bronchial epithelial cell growth media (BEGM) only, co-cultured with irradiated mouse embryonic fibroblasts (IRR-MEF) + Rho kinase inhibitor (ROCKi) in BEGM, and conditioned BEGM + ROCKi. Media were changed every three days and cells were passaged and sub-cultured after ten days. Human telomerase reverse transcriptase (hTERT) was measured in three individuals after each passage in triplicate via qPCR. The proliferative capacity of all eight individuals was assessed using cell count and morphology at passage >3. Results. Co-culturing NBEC with IRR-MEF in BEGM supplemented with ROCKi produced a highly proliferative cell population while maintaining lineage commitment evident after removal of CRC conditions. Transcript abundance of hTERT was elevated 6.5-fold in NBEC in co-cultured conditions and 4.3-fold in NBEC in conditioned media compared to BEGM alone. Cell count in CRC conditions were up to 22-fold higher compared to BEGM alone. Cells were passaged and sub-cultured up to passage 4, followed by being frozen down in cell culture freezing media for further assessment. Conclusion. NBEC hTERT transcript abundance was up-regulated and cell population proliferative potential was extended in CRC conditions. It is likely that hTERT functions to protect the ends of linear chromosomes in dividing cells, enabling increased cell divisions while maintaining normal genome. These cell populations will be used in future studies to assess the effect of putative cis-regulatory single nucleotide polymorphisms (SNPs) on gene expression in NBEC. Citation Format: Daniel J. Craig, Rose T. Zolondek, Xiaolu Zhang, Jiyoun Yeo, Erin L. Crawford, James C. Willey. Induction of hTERT and increased proliferative potential in conditionally reprogrammed normal bronchial epithelial cells. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2905.

Cell adhesion molecule pathway genes are regulated by cis-regulatory SNPs and show significantly altered expression in Alzheimer's disease brains

We previously identified the cell adhesion molecule (CAM) pathway as a consistent signal in 2 Alzheimer's disease (AD) genome-wide association studies (GWAS). However, the genetic mechanisms of the CAM pathway in AD are unclear. Here, we conducted pathway analysis using (1) Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathways; (2) 4 brain expression GWAS datasets; and (3) 2 whole-genome AD case-control expression datasets. Using the 4 brain expression GWAS datasets, we identified that genes regulated by cis-regulatory single-nucleotide polymorphisms (SNPs) were significantly enriched in the CAM pathway (p = 2.05E−06, p = 6.10E−07, p = 2.05E−06, and p = 1.47E−07 for each dataset). Interestingly, CAM is a significantly enriched pathway using down-regulated genes (raw p = 0.0235 and adjusted p = 0.0305) and all differentially expressed genes (raw p = 0.0105 and adjusted p = 0.0156) in dataset 5, and all differentially expressed genes (raw p = 0.0041 and adjusted p = 0.0062) in dataset 6. Collectively, our results show that CAM pathway genes are regulated by cis-regulatory SNPs and show significantly altered expression in AD. We believe that our results advance the understanding of AD mechanisms and will be useful for future genetic studies of AD.

Abstract 4156: Inter-individual variation in allele specific expression of catalase (CAT) in normal bronchial epithelial cells and association of putative cis-regulatory CAT SNP rs12807961 with lung cancer risk

Abstract Background: Catalase (CAT) is a key antioxidant gene expressed at high levels in most human tissues, including normal bronchial epithelial cells (NBEC). NBEC CAT expression is more disperse (more high and low extreme values) among subjects with cancer compared to controls. CAT shares this property with 14 other antioxidant and DNA repair genes comprised by the Lung Cancer Risk Test (LCRT) reported from this laboratory. We hypothesize that inter-individual variation in CAT regulation in NBEC is in part due to inter-individual variation at one or more cis-regulatory single nucleotide polymorphisms (SNPs). If so, this should manifest as inter-individual variation in allele-specific CAT expression in NBEC. Through funding in part from RC2 CA148572 and HL108016 we collected NBEC samples from over 500 subjects at risk for lung cancer. In this pilot study, we assessed allele-specific and total expression of multiple genes in NBEC samples from 85 subjects and assessed the genotype at putative cis-regulatory CAT SNP rs12807961 in gDNA from 95 subjects. Methods: RNA extracted from normal bronchial airway brush specimens of 85 subjects (26 cancer cases and 59 non-cancer controls) was reverse transcribed. Using next generation sequencing (NGS), allele-specific expression (ASE) was measured as allelic imbalance in each cDNA at three marker SNPs in the CAT coding region (rs1049982, rs769217, and rs704724) and one putative regulatory SNP (rs12807961) that was 4364 bases upstream of transcription start site, using gDNA as control. Specifically, each cDNA and matched peripheral blood cell gDNA sample was subjected to targeted competitive template multiplex PCR amplicon library generation followed by NGS (Blomquist et al, PLOS one, 2013) on Illumina Hiseq platform. The genotype at putative cis-regulatory SNP rs12807961 was assessed in gDNA from 95 subjects including those assessed for ASE (a total of 31 cancer cases and 64 non-cancer controls) using a TaqMan® SNP Genotyping Assay. Results: Among heterozygotes, there was significant inter-individual variation in cDNA allelic imbalance at rs1049982 (p<0.001, n=40) and rs769217 (p<0.001, n=28) as measured by F-test using matched gDNA controls. In this cohort there was insufficient number of heterozygotes at rs704724 (n=2) to assess ASE. Among all 95 subjects assessed for rs12807961 genotype, nine were homozygous minor allele at the rs12807961 CAT SNP. Of these, 7/31 cancer cases and 2/64 non-cancer controls were homozygous minor allele. This difference was significant by two-tailed Fisher exact test (P<0.05) following Bonferroni adjustment for multiple testing. Conclusions: These data support the hypothesis that cis-regulatory DNA variants contribute to inter-individual variation in CAT regulation in NBEC and that this is associated with lung cancer risk. Citation Format: Jiyoun Yeo, Xaiolu Zhang, Erin L. Crawford, James C. Willey. Inter-individual variation in allele specific expression of catalase (CAT) in normal bronchial epithelial cells and association of putative cis-regulatory CAT SNP rs12807961 with lung cancer risk. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4156. doi:10.1158/1538-7445.AM2014-4156

Abstract 4251: Identification of expression quantitative trait loci at lung cancer and COPD risk genes in normal bronchial epithelial cells

Abstract Background: Lung cancer and chronic obstructive pulmonary disease (COPD) and are leading causes of morbidity and mortality both in the United States and worldwide. Inhalation of cigarette smoke is the primary known and preventable cause but only 10-15% of heavy smokers develop these diseases. This suggests that cigarette smoke exposure interacts with inherited susceptibility factors to determine risk. While some susceptibility genes are known, they account for less than 5% of risk for either disease. There is urgent need to identify heritable susceptibility factors that explain the majority of lung cancer and COPD risk. In this study we focused on discovery of cis-regulatory expression quantitative trait loci (eQTL) responsible for inter-individual variation in the expression of genes that we and/or others have reported to display altered regulation in subjects with lung cancer or COPD. Through funding in part from RC2 CA148572 and HL108016 we collected normal bronchial epithelial cell (NBEC) samples from over 500 subjects with lung cancer or COPD or demographically at risk for lung cancer or COPD. In this pilot study, we assessed allele-specific expression (ASE) and total expression of multiple putative risk genes in NBEC samples from 85 subjects. Methods: A targeted competitive multiplex next generation sequencing (NGS) method (STAndardized RNA SEQuencing [STARSEQ]; Blomquist et al, PLOS one 2013) was used to quantify a) allele specific expression (ASE) and total expression at 140 single nucleotide polymorphism (SNP) sites among 41 target genes (including eleven of 14 genes comprised by the previously reported lung cancer risk test (Blomquist et al, Can. Res. 2009) in NBEC total RNA from 85 subjects (26 cancer cases and 59 non-cancer controls), and b) allelic representation and putative cis-regulatory single nucleotide polymorphisms (rSNPs) in matched peripheral blood leukocyte genomic (g)DNA from the same subjects. Heterozygosity was determined by gDNA analysis. Results: Significant (p<0.05) inter-individual variation in ASE measured as allelic imbalance relative to matched gDNA was observed for: CAT, CCND2, ERCC4, ERCC5, KEAP1, MUC5B and TP73. Among these genes, for some the same allele was typically expressed at a higher level among the subjects studied (e.g. ERCC5 rs17655) which suggests a closely proximal cis-regulatory SNP and for others the over expressed allele was random (e.g. CAT rs1049982) which suggests a more distal cis regulatory SNP. Conclusions: This small study of 85 subjects provides evidence for cis-regulatory eQTL in selected genes with high prior likelihood for contributing to lung cancer and COPD risk. The genes with ASE reported here will be among those included in additional studies in archived NBEC samples from 500 subjects aiming to further optimize current tests for lung cancer and COPD risk and better understand mechanisms of risk. Citation Format: Erin L. Crawford, Jiyoun Yeo, Xiaolu Zhang, Karan Padda, Taylor Arend, Thomas M. Blomquist, Albert M. Levin, Mei Lu, James C. Willey. Identification of expression quantitative trait loci at lung cancer and COPD risk genes in normal bronchial epithelial cells. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4251. doi:10.1158/1538-7445.AM2014-4251

Global survey on sequence characteristics of plant microRNA genes: Cis-regulatory SNPs in promoters and microRNA precursors

MicroRNAs (miRNAs) are tiny molecules with multiple essential roles in plant growth and development. So far, numerous miRNA families, both highly conserved and species-specific, have been identified and functionally characterized. However, we are still far from a thorough understanding of these molecules, including their sequence characteristics. In this bioinformatics study, cis-acting regulatory elements located within plant miRNA promoters were predicted, and a comprehensive list of single nucleotide polymorphisms (SNPs) residing within these cis-elements, designated as cis-regulatory SNPs, was identified. These SNPs potentially contribute to the distinct expression patterns of miRNA genes in various Arabidopsis (Arabidopsis thaliana) or rice (Oryza sativa) subspecies. Additionally, sequence conservation analysis was carried out for all the currently miRBase-registered pre-miRNAs (precursor miRNAs) belonging to numerous plant species. Together, these analyses provide a basis to gain deeper insights into the biological role of miRNA sequence characteristics.

Effects of BRCA2 cis-regulation in normal breast and cancer risk amongst BRCA2 mutation carriers

IntroductionCis-acting regulatory single nucleotide polymorphisms (SNPs) at specific loci may modulate penetrance of germline mutations at the same loci by introducing different levels of expression of the wild-type allele. We have previously reported that BRCA2 shows differential allelic expression and we hypothesize that the known variable penetrance of BRCA2 mutations might be associated with this mechanism.MethodsWe combined haplotype analysis and differential allelic expression of BRCA2 in breast tissue to identify expression haplotypes and candidate cis-regulatory variants. These candidate variants underwent selection based on in silico predictions for regulatory potential and disruption of transcription factor binding, and were functionally analyzed in vitro and in vivo in normal and breast cancer cell lines. SNPs tagging the expression haplotypes were correlated with the total expression of several genes in breast tissue measured by Taqman and microarray technologies. The effect of the expression haplotypes on breast cancer risk in BRCA2 mutation carriers was investigated in 2,754 carriers.ResultsWe identified common haplotypes associated with differences in the levels of BRCA2 expression in human breast cells. We characterized three cis-regulatory SNPs located at the promoter and two intronic regulatory elements which affect the binding of the transcription factors C/EBPα, HMGA1, D-binding protein (DBP) and ZF5. We showed that the expression haplotypes also correlated with changes in the expression of other genes in normal breast. Furthermore, there was suggestive evidence that the minor allele of SNP rs4942440, which is associated with higher BRCA2 expression, is also associated with a reduced risk of breast cancer (per-allele hazard ratio (HR) = 0.85, 95% confidence interval (CI) = 0.72 to 1.00, P-trend = 0.048).ConclusionsOur work provides further insights into the role of cis-regulatory variation in the penetrance of disease-causing mutations. We identified small-effect genetic variants associated with allelic expression differences in BRCA2 which could possibly affect the risk in mutation carriers through altering expression levels of the wild-type allele.

SNP discovery, expression and cis-regulatory variation in the UGT2B genes

UGT2B enzymes metabolize multiple endogenous and exogenous molecules, including steroid hormones and clinical drugs. However, little is known about the inter-individual variation in gene expression and its determinants. We re-sequenced candidate regulatory regions and the partial coding regions (41.1 kb) of UGT2B genes and identified 332 genetic variants. We measured gene expression in normal breast and liver samples and observed different patterns. The expression levels varied greatly across individuals in both tissues and were significantly correlated with each other in liver. Genotyping of tagging single-nucleotide polymorphisms (SNPs) in the same samples and association tests between genotype and transcript levels identified 62 variants that were associated with at least one UGT2B mRNA levels in either tissue. Most of these cis-regulatory SNPs were not shared between tissues, suggesting that this gene family is regulated in a tissue-specific manner. Our results provide insight into studying the role of UGT2B variation in hormone-dependent cancers and drug response.

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.

A Catalog of Neutral and Deleterious Polymorphism in Yeast

The abundance and identity of functional variation segregating in natural populations is paramount to dissecting the molecular basis of quantitative traits as well as human genetic diseases. Genome sequencing of multiple organisms of the same species provides an efficient means of cataloging rearrangements, insertion, or deletion polymorphisms (InDels) and single-nucleotide polymorphisms (SNPs). While inbreeding depression and heterosis imply that a substantial amount of polymorphism is deleterious, distinguishing deleterious from neutral polymorphism remains a significant challenge. To identify deleterious and neutral DNA sequence variation within Saccharomyces cerevisiae, we sequenced the genome of a vineyard and oak tree strain and compared them to a reference genome. Among these three strains, 6% of the genome is variable, mostly attributable to variation in genome content that results from large InDels. Out of the 88,000 polymorphisms identified, 93% are SNPs and a small but significant fraction can be attributed to recent interspecific introgression and ectopic gene conversion. In comparison to the reference genome, there is substantial evidence for functional variation in gene content and structure that results from large InDels, frame-shifts, and polymorphic start and stop codons. Comparison of polymorphism to divergence reveals scant evidence for positive selection but an abundance of evidence for deleterious SNPs. We estimate that 12% of coding and 7% of noncoding SNPs are deleterious. Based on divergence among 11 yeast species, we identified 1,666 nonsynonymous SNPs that disrupt conserved amino acids and 1,863 noncoding SNPs that disrupt conserved noncoding motifs. The deleterious coding SNPs include those known to affect quantitative traits, and a subset of the deleterious noncoding SNPs occurs in the promoters of genes that show allele-specific expression, implying that some cis-regulatory SNPs are deleterious. Our results show that the genome sequences of both closely and distantly related species provide a means of identifying deleterious polymorphisms that disrupt functionally conserved coding and noncoding sequences.

Allele-Specific Up-Regulation of FGFR2 Increases Susceptibility to Breast Cancer

The recent whole-genome scan for breast cancer has revealed the FGFR2 (fibroblast growth factor receptor 2) gene as a locus associated with a small, but highly significant, increase in the risk of developing breast cancer. Using fine-scale genetic mapping of the region, it has been possible to narrow the causative locus to a haplotype of eight strongly linked single nucleotide polymorphisms (SNPs) spanning a region of 7.5 kilobases (kb) in the second intron of the FGFR2 gene. Here we describe a functional analysis to define the causative SNP, and we propose a model for a disease mechanism. Using gene expression microarray data, we observed a trend of increased FGFR2 expression in the rare homozygotes. This trend was confirmed using real-time (RT) PCR, with the difference between the rare and the common homozygotes yielding a Wilcox p-value of 0.028. To elucidate which SNPs might be responsible for this difference, we examined protein–DNA interactions for the eight most strongly disease-associated SNPs in different breast cell lines. We identify two cis-regulatory SNPs that alter binding affinity for transcription factors Oct-1/Runx2 and C/EBPβ, and we demonstrate that both sites are occupied in vivo. In transient transfection experiments, the two SNPs can synergize giving rise to increased FGFR2 expression. We propose a model in which the Oct-1/Runx2 and C/EBPβ binding sites in the disease-associated allele are able to lead to an increase in FGFR2 gene expression, thereby increasing the propensity for tumour formation.

Allele-Specific KRT1 Expression Is a Complex Trait

The differential expression of alleles occurs commonly in humans and is likely an important genetic factor underlying heritable differences in phenotypic traits. Understanding the molecular basis of allelic expression differences is thus an important challenge. Although many genes have been shown to display differential allelic expression, this is the first study to examine in detail the cumulative effects of multiple cis-regulatory polymorphisms responsible for allele-specific expression differences. We have used a variety of experimental approaches to identify and characterize cis-regulatory polymorphisms responsible for the extreme allele-specific expression differences of keratin-1 (KRT1) in human white blood cells. The combined data from our analyses provide strong evidence that the KRT1 allelic expression differences result from the haplotypic combinations and interactions of five cis-regulatory single nucleotide polymorphisms (SNPs) whose alleles differ in their affinity to bind transcription factors and modulate KRT1 promoter activity. Two of these cis-regulatory SNPs bind transcriptional activators with the alleles on the high-expressing KRT1 haplotype pattern having a higher affinity than the alleles on the low-expressing haplotype pattern. In contrast, the other three cis-regulatory SNPs bind transcriptional inhibitors with the alleles on the low-expressing haplotype pattern having a higher affinity than the alleles on the high-expressing haplotype pattern. Our study provides important new insights into the degree of complexity that the cis-regulatory sequences responsible for allele-specific transcriptional regulation have. These data suggest that allelic expression differences result from the cumulative contribution of multiple DNA sequence polymorphisms, with each having a small effect, and that allele-specific expression can thus be viewed as a complex trait.