Mapping of PARK2 and PACRG Overlapping Regulatory Region Reveals LD Structure and Functional Variants in Association with Leprosy in Unrelated Indian Population Groups

Leprosy is a chronic infectious disease caused by Mycobacterium Leprae, where the host genetic background plays an important role toward the disease pathogenesis. Various studies have identified a number of human genes in association with leprosy or its clinical forms. However, non-replication of results has hinted at the heterogeneity among associations between different population groups, which could be due to differently evolved LD structures and differential frequencies of SNPs within the studied regions of the genome. A need for systematic and saturated mapping of the associated regions with the disease is warranted to unravel the observed heterogeneity in different populations. Mapping of the PARK2 and PACRG gene regulatory region with 96 SNPs, with a resolution of 1 SNP per 1 Kb for PARK2 gene regulatory region in a North Indian population, showed an involvement of 11 SNPs in determining the susceptibility towards leprosy. The association was replicated in a geographically distinct and unrelated population from Orissa in eastern India. In vitro reporter assays revealed that the two significantly associated SNPs, located 63.8 kb upstream of PARK2 gene and represented in a single BIN of 8 SNPs, influenced the gene expression. A comparison of BINs between Indian and Vietnamese populations revealed differences in the BIN structures, explaining the heterogeneity and also the reason for non-replication of the associated genomic region in different populations.     

Genetic meta-analysis of 15,901 African Americans identifies variation in EXOC3L1 is associated with HDL concentration

Meta-analyses of European populations has successfully identified genetic variants in over 150 loci associated with lipid levels, but results from additional ethnicities remain limited. Previously, we reported two novel lipid loci identified in a sample of 7,657 African Americans using a gene-centric array including 50,000 SNPs in 2,100 candidate genes. Initial discovery and follow-up of signals with P < 10⁻⁵ in additional African American samples confirmed CD36 and ICAM1. Using an additional 8,244 African American female samples from the Women’s Health Initiative SNP Health Association Resource genome-wide association study dataset, we further examined the previous meta-analyses results by attempting to replicate 20 additional putative lipid signals with P < 10⁻⁴. Replication confirmed rs868213, located in a splice donor region of exocyst complex component 3-like 1 (EXOC3L1) as a novel signal for HDL (additive allelic effect β = 0.02; P = 1.4 × 10⁻⁸; meta-analyses of discovery and replication). EXOC3L1 is strongly expressed in vascular endothelium and forms part of the exocyst complex, a key facilitator of the trafficking of lipid receptors. Increasing sample sizes for genetic studies in nonEuropean populations will continue to improve our understanding of lipid metabolism.     

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.     

A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10^?5) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.     

Association between the g.296596G > A genetic variant of RELN gene and susceptibility to autism in a Chinese Han population

Autism is a childhood neuro-developmental disorder, and Reelin (RELN) is an important candidate gene for influencing autism. This study aimed at investigating the influence of genetic variants of the RELN gene on autism susceptibility. In this study, 205 autism patients and 210 healthy controls were recruited and the genetic variants of the RELN gene were genotyped by the created restriction site-polymerase chain reaction (CRS-PCR) method. The influence of genetic variants on autism susceptibility was analyzed by association analysis, and the g.296596G > A genetic variant in exon10 of the RELN gene was detected. The frequencies of allele/genotype in autistic patients were significantly different from those in healthy controls, and a statistically significant association was detected between this genetic variant and autism susceptibility. Our data lead to the inference that the g.296596G > A genetic variant in the RELN gene has a potential influence on autism susceptibility in the Chinese Han population.     

Combination Testing Using a Single MSH5 Variant alongside HLA Haplotypes Improves the Sensitivity of Predicting Coeliac Disease Risk in the Polish Population

Assessment of non-HLA variants alongside standard HLA testing was previously shown to improve the identification of potential coeliac disease (CD) patients. We intended to identify new genetic variants associated with CD in the Polish population that would improve CD risk prediction when used alongside HLA haplotype analysis. DNA samples of 336 CD and 264 unrelated healthy controls were used to create DNA pools for a genome wide association study (GWAS). GWAS findings were validated with individual HLA tag single nucleotide polymorphism (SNP) typing of 473 patients and 714 healthy controls. Association analysis using four HLA-tagging SNPs showed that, as was found in other populations, positive predicting genotypes (HLA-DQ2.5/DQ2.5, HLA-DQ2.5/DQ2.2, and HLA-DQ2.5/DQ8) were found at higher frequencies in CD patients than in healthy control individuals in the Polish population. Both CD-associated SNPs discovered by GWAS were found in the CD susceptibility region, confirming the previously-determined association of the major histocompatibility (MHC) region with CD pathogenesis. The two most significant SNPs from the GWAS were rs9272346 (HLA-dependent; localized within 1 Kb of DQA1) and rs3130484 (HLA-independent; mapped to MSH5). Specificity of CD prediction using the four HLA-tagging SNPs achieved 92.9%, but sensitivity was only 45.5%. However, when a testing combination of the HLA-tagging SNPs and the MSH5 SNP was used, specificity decreased to 80%, and sensitivity increased to 74%. This study confirmed that improvement of CD risk prediction sensitivity could be achieved by including non-HLA SNPs alongside HLA SNPs in genetic testing.     

Heterozygous Reelin Mutations Cause Autosomal-Dominant Lateral Temporal Epilepsy

Autosomal-dominant lateral temporal epilepsy (ADLTE) is a genetic epilepsy syndrome clinically characterized by focal seizures with prominent auditory symptoms. ADLTE is genetically heterogeneous, and mutations in LGI1 account for fewer than 50% of affected families. Here, we report the identification of causal mutations in reelin (RELN) in seven ADLTE-affected families without LGI1 mutations. We initially investigated 13 ADLTE-affected families by performing SNP-array linkage analysis and whole-exome sequencing and identified three heterozygous missense mutations co-segregating with the syndrome. Subsequent analysis of 15 small ADLTE-affected families revealed four additional missense mutations. 3D modeling predicted that all mutations have structural effects on protein-domain folding. Overall, RELN mutations occurred in 7/40 (17.5%) ADLTE-affected families. RELN encodes a secreted protein, Reelin, which has important functions in both the developing and adult brain and is also found in the blood serum. We show that ADLTE-related mutations significantly decrease serum levels of Reelin, suggesting an inhibitory effect of mutations on protein secretion. We also show that Reelin and LGI1 co-localize in a subset of rat brain neurons, supporting an involvement of both proteins in a common molecular pathway underlying ADLTE. Homozygous RELN mutations are known to cause lissencephaly with cerebellar hypoplasia. Our findings extend the spectrum of neurological disorders associated with RELN mutations and establish a link between RELN and LGI1, which play key regulatory roles in both the developing and adult brain.     

Finding Missing Heritability in Less Significant Loci and Allelic Heterogeneity: Genetic Variation in Human Height

Genome-wide association studies (GWAS) have identified many common variants associated with complex traits in human populations. Thus far, most reported variants have relatively small effects and explain only a small proportion of phenotypic variance, leading to the issues of ‘missing’ heritability and its explanation. Using height as an example, we examined two possible sources of missing heritability: first, variants with smaller effects whose associations with height failed to reach genome-wide significance and second, allelic heterogeneity due to the effects of multiple variants at a single locus. Using a novel analytical approach we examined allelic heterogeneity of height-associated loci selected from SNPs of different significance levels based on the summary data of the GIANT (stage 1) studies. In a sample of 1,304 individuals collected from an island population of the Adriatic coast of Croatia, we assessed the extent of height variance explained by incorporating the effects of less significant height loci and multiple effective SNPs at the same loci. Our results indicate that approximately half of the 118 loci that achieved stringent genome-wide significance (p-value<5×10⁻⁸) showed evidence of allelic heterogeneity. Additionally, including less significant loci (i.e., p-value<5×10⁻⁴) and accounting for effects of allelic heterogeneity substantially improved the variance explained in height.     

Identification of a 31-bp Deletion in the RELN Gene Causing Lissencephaly with Cerebellar Hypoplasia in Sheep

Lissencephaly is an inherited developmental disorder in which neuronal migration is impaired. A type of lissencephaly associated with cerebellar hypoplasia (LCH) was diagnosed in a commercial flock of Spanish Churra sheep. The genotyping of 7 affected animals and 33 controls with the OvineSNP50 BeadChip enabled the localization of the causative mutation for ovine LCH to a 4.8-Mb interval on sheep chromosome 4 using genome-wide association and homozygosity mapping. The RELN gene, which is located within this interval, was considered a strong positional and functional candidate because it plays critical roles in neuronal migration and layer formation. By performing a sequencing analysis of this gene’s specific mRNA in a control lamb, we obtained the complete CDS of the ovine RELN gene. The cDNA sequence from an LCH-affected lamb revealed a deletion of 31 bp (c.5410_5440del) in predicted exon 36 of RELN, resulting in a premature termination codon. A functional analysis of this mutation revealed decreased levels of RELN mRNA and a lack of reelin protein in the brain cortex and blood of affected lambs. This mutation showed a complete concordance with the Mendelian recessive pattern of inheritance observed for the disease. The identification of the causal mutation of LCH in Churra sheep will facilitate the implementation of gene-assisted selection to detect heterozygous mutants, which will help breeders avoid at-risk matings in their flocks. Moreover, the identification of this naturally occurring RELN mutation provides an opportunity to use Churra sheep as a genetically characterized large animal model for the study of reelin functions in the developing and mature brain.     

Evaluation of European Schizophrenia GWAS Loci in Asian Populations via Comprehensive Meta-Analyses

Schizophrenia is a severe and highly heritable neuropsychiatric disorder. Recent genetic analyses including genome-wide association studies (GWAS) have implicated multiple genome-wide significant variants for schizophrenia among European populations. However, many of these risk variants were not largely validated in other populations of different ancestry such as Asians. To validate whether these European GWAS significant loci are associated with schizophrenia in Asian populations, we conducted a systematic literature search and meta-analyses on 19 single nucleotide polymorphisms (SNPs) in Asian populations by combining all available case-control and family-based samples, including up to 30,000 individuals. We employed classical fixed (or random) effects inverse variance weighted methods to calculate summary odds ratios (ORs) and 95 % confidence intervals (CIs). Among the 19 GWAS loci, we replicated the risk associations of nine markers (e.g., SNPs at VRK2, ITIH3/4, NDST3, NOTCH4) surpassing significance level (two-tailed P?<?0.05), and three additional SNPs in MIR137 and ZNF804A also showed trend associations (one-tailed P?<?0.05). These risk associations are in the same directions of allelic effects between Asian replication samples and initial European GWAS findings, and the successful replications of these GWAS loci in a different ethnic group provide stronger evidence for their clinical associations with schizophrenia. Further studies, focusing on the molecular mechanisms of these GWAS significant loci, will become increasingly important for understanding of the pathogenesis to schizophrenia.     

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus—a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST. The number of SNPs typically obtained from next‐generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.     

A Genome-Wide Association Study Suggests Novel Loci Associated with a Schizophrenia-Related Brain-Based Phenotype

Patients with schizophrenia and their siblings typically show subtle changes of brain structures, such as a reduction of hippocampal volume. Hippocampal volume is heritable, may explain a variety of cognitive symptoms of schizophrenia and is thus considered an intermediate phenotype for this mental illness. The aim of our analyses was to identify single-nucleotide polymorphisms (SNP) related to hippocampal volume without making prior assumptions about possible candidate genes. In this study, we combined genetics, imaging and neuropsychological data obtained from the Mind Clinical Imaging Consortium study of schizophrenia (n = 328). A total of 743,591 SNPs were tested for association with hippocampal volume in a genome-wide association study. Gene expression profiles of human hippocampal tissue were investigated for gene regions of significantly associated SNPs. None of the genetic markers reached genome-wide significance. However, six highly correlated SNPs (rs4808611, rs35686037, rs12982178, rs1042178, rs10406920, rs8170) on chromosome 19p13.11, located within or in close proximity to the genes NR2F6, USHBP1, and BABAM1, as well as four SNPs in three other genomic regions (chromosome 1, 2 and 10) had p-values between 6.75×10⁻⁶ and 8.3×10⁻⁷. Using existing data of a very recently published GWAS of hippocampal volume and additional data of a multicentre study in a large cohort of adolescents of European ancestry, we found supporting evidence for our results. Furthermore, allelic differences in rs4808611 and rs8170 were highly associated with differential mRNA expression in the cis-acting region. Associations with memory functioning indicate a possible functional importance of the identified risk variants. Our findings provide new insights into the genetic architecture of a brain structure closely linked to schizophrenia. In silico replication, mRNA expression and cognitive data provide additional support for the relevance of our findings. Identification of causal variants and their functional effects may unveil yet unknown players in the neurodevelopment and the pathogenesis of neuropsychiatric disorders.     

Association of PDE4B Polymorphisms with Susceptibility to Schizophrenia: A Meta-Analysis of Case-Control Studies

Background

The PDE4B single nucleotide polymorphisms (SNPs) have been reported to be associated with schizophrenia risk. However, current findings are ambiguous or even conflicting. To better facilitate the understanding the genetic role played by PDE4B in susceptibility to schizophrenia, we collected currently available data and conducted this meta-analysis.

Methods

A comprehensive electronic literature searching of PubMed, Embase, Web of Science and Cochrane Library was performed. The association between PDE4B SNPs and schizophrenia was evaluated by odds ratios (ORs) and 95% confidence intervals (CIs) under allelic, dominant and recessive genetic models. The random effects model was utilized when high between-study heterogeneity (I² > 50%) existed, otherwise the fixed effects model was used.

Results

Five studies comprising 2376 schizophrenia patients and 3093 controls were finally included for meta-analysis. The rs1040716 was statistically significantly associated with schizophrenia risk in Asian and Caucasian populations under dominant model (OR = 0.87, 95% CI: 0.76–0.99, P = 0.04). The rs2180335 was significantly related with schizophrenia risk in Asian populations under allelic (OR = 0.82, 95% CI: 0.72–0.93, P = 0.003) and dominant (OR = 0.75, 95% CI: 0.64–0.88, P < 0.001) models. A significant association was also observed between rs4320761 and schizophrenia in Asian populations under allelic model (OR = 0.87, 95% CI: 0.75–1.00, P = 0.048). In addition, a strong association tendency was found between rs6588190 and schizophrenia in Asian populations under allelic model (OR = 0.87, 95% CI: 0.76–1.00, P = 0.055).

Conclusion

This meta-analysis suggests that PDE4B SNPs are genetically associated with susceptibility to schizophrenia. However, due to limited sample size, more large-scale, multi-racial association studies are needed to further clarify the genetic association between various PDE4B variants and schizophrenia.     

Association of common variants in KIF21B and ankylosing spondylitis in a Chinese Han population: a replication study

KIF21B polymorphisms were found associated with susceptibility to multiple sclerosis and ankylosing spondylitis (AS) in populations of white European ancestry. We aimed to replicate the association of polymorphisms around KIF21B and AS in a Chinese Han population. This case–control study included 665 patients with AS and 1,042 healthy controls genotyped for seven single nucleotide polymorphisms (SNPs) of KIF21B—rs12118246, rs4915464, rs502658, rs10494829, rs12089839, rs6687260, and rs957957—by TaqMan genotyping assay; statistical analyses involved the use of PLINK. We also estimated the linkage disequilibrium and haplotypes of these SNPs. Two SNPs—rs502658 (allelic p?=?0.0002, odds ratio [OR] 0.60, 95 % confidence interval [95 % CI] 0.47–0.76) and rs10494829 (allelic p?=?0.003, OR 1.30, 95 % CI 1.12–1.52)—were significantly associated with AS in the Chinese Han population. In addition, a linear regression test showed that they have independent contribution to disease susceptibility. For both SNPs, haplotype AT was strongly associated with AS and increased the risk of the disease (p?=?0.045, OR 1.183, 95 % CI 1.004–1.395), and the genotype GC reduced the risk (p?=?0.011, OR 0.715, 95 % CI 0.55–0.928). This work identified a significant association of two SNPs in KIF21B and AS in the Chinese Han population. KIF21B may play an important role in the pathogenesis of AS in the Chinese population and might be a new therapeutic target for AS.     

The First Report on Single Nucleotide Polymorphisms of the <Emphasis Type="Italic">HSPA13</Emphasis> Gene in Koreans

Heat shock proteins (HSPs) are known as molecular chaperones, and they function in response to cell stress. HSPA13, also called STCH, is a member of the HSP70 family. In general, HSP70 family may play a protective role in prion diseases. In a recent study, the overexpression of HSPA13 was shown to shorten the incubation time of prion diseases. Although the exact role of HSPA13 in the pathogenesis of prion diseases remains unknown, the expression level of HSPA13 is significantly associated with the latent phase of prion diseases. It has been known that single nucleotide polymorphisms (SNPs) in promoter and open reading frame (ORF) region of genes can affect either gene expression or gene function. The purpose of this study was to investigate genotype and allele frequencies of SNPs found in the promoter and ORF of HSPA13 in healthy Korean population to obtain the information for subsequent population genetics and prion diseases studies. We observed four SNPs in the promoter region of HSPA13, of which two have previous identified (c.-608C>G; rs2242662 and c.-381G>A; rs2242661) and two are novel (c.-321C>T and c.-300A>G). Interestingly, we did not observe any polymorphisms in the ORF of this gene. To our knowledge, this is the first study of polymorphisms in the human HSPA13 gene.     

Association of single-nucleotide polymorphisms in RHOB and TXNDC3 with knee osteoarthritis susceptibility: two case-control studies in East Asian populations and a meta-analysis

Introduction

Conflicting findings on the association of single nucleotide polymorphisms (SNPs) in RHOB and TXNDC3 with susceptibility to knee osteoarthritis (OA) have been reported in European Caucasians. To examine the associations of these SNPs with OA in East Asian populations and to evaluate their global significance, we conducted two case-control studies in 955 Chinese and 750 Japanese patients.

Methods

We genotyped the previously implicated SNPs rs585017 (in RHOB) and rs4720262 (in TXNDC3) in patients with primary symptomatic knee OA with radiographic confirmation and in matched control individuals, and analyzed their associations. We further conducted a meta-analysis of the study findings together with those of previously reported European studies using the DerSimonian-Laird procedure.

Results

A significant association of RHOB with knee OA was observed in male Chinese patients (P = 0.02). No significant associations were found for RHOB in any other comparisons in the East Asian populations. The association of TXNDC3 was replicated in Chinese female (P = 0.04) and Japanese (P = 0.03) patients, although none of these associations persisted after Bonferroni correction. Significant association (P = 0.02 for the allelic frequency) with nonsignificant heterogeneity was found in the East Asian replication study. No significant association was found in any comparison in the meta-analysis for all studies.

Conclusion

Our study replicates the association, previously reported in European Caucasians, of TXNDC3 with knee OA susceptibility in an East Asian population.     

Role of Lung Function Genes in the Development of Asthma

Although our previous GWAS failed to identify SNPs associated with pulmonary function at the level of genomewide significance, it did show that the heritability for FEV₁/FVC was 41.6% in a Japanese population, suggesting that the heritability of pulmonary function traits can be explained by the additive effects of multiple common SNPs. In addition, our previous study indicated that pulmonary function genes identified in previous GWASs in non-Japanese populations accounted for 4.3% to 12.0% of the entire estimated heritability of FEV₁/FVC in a Japanese population. Therefore, given that many loci with individual weak effects may contribute to asthma risk, in this study, we created a quantitative score of genetic load based on 16 SNPs implicated in lower lung function in both Japanese and non-Japanese populations. This genetic risk score (GRS) for lower FEV₁/FVC was consistently associated with the onset of asthma (P = 9.6 × 10⁻⁴) in 2 independent Japanese populations as well as with the onset of COPD (P = 0.042). Clustering of asthma patients based on GRS levels indicated that an increased GRS may be responsible for the development of a particular phenotype of asthma characterized by early onset, atopy, and severer airflow obstruction.     

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.