title: Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing.
abstract: Many disease-associated variants identified by genome-wide association (GWA) studies are expected to regulate gene expression. Allele-specific expression (ASE) quantifies transcription from both haplotypes using individuals heterozygous at tested SNPs. We performed deep human transcriptome-wide resequencing (RNA-seq) for ASE analysis and expression quantitative trait locus discovery. We resequenced double poly(A)-selected RNA from primary CD4(+) T cells (n = 4 individuals, both activated and untreated conditions) and developed tools for paired-end RNA-seq alignment and ASE analysis. We generated an average of 20 million uniquely mapping 45 base reads per sample. We obtained sufficient read depth to test 1371 unique transcripts for ASE. Multiple biases inflate the false discovery rate which we estimate to be approximately 50% for random SNPs. However, after controlling for these biases and considering the subset of SNPs that pass HapMap QC, 4.6% of heterozygous SNP-sample pairs show evidence of imbalance (P < 0.001). We validated four findings by both bacterial cloning and Sanger sequencing assays. We also found convincing evidence for allelic imbalance at multiple reporter exonic SNPs in CD6 for two samples heterozygous at the multiple sclerosis-associated variant rs17824933, linking GWA findings with variation in gene expression. Finally, we show in CD4(+) T cells from a further individual that high-throughput sequencing of genomic DNA and RNA-seq following enrichment for targeted gene sequences by sequence capture methods offers an unbiased means to increase the read depth for transcripts of interest, and therefore a method to investigate the regulatory role of many disease-associated genetic variants.
authors.affiliation: Centre for Digestive Diseases, Blizard Institute of Cell and Molecular Science, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, London, UK.
authors.forename: Graham A
authors.initials: GA
authors.lastname: Heap
authors.forename: Jennie H M
authors.initials: JH
authors.lastname: Yang
authors.forename: Kate
authors.initials: K
authors.lastname: Downes
authors.forename: Barry C
authors.initials: BC
authors.lastname: Healy
authors.forename: Karen A
authors.initials: KA
authors.lastname: Hunt
authors.forename: Nicholas
authors.initials: N
authors.lastname: Bockett
authors.forename: Lude
authors.initials: L
authors.lastname: Franke
authors.forename: Patrick C
authors.initials: PC
authors.lastname: Dubois
authors.forename: Charles A
authors.initials: CA
authors.lastname: Mein
authors.forename: Richard J
authors.initials: RJ
authors.lastname: Dobson
authors.forename: Thomas J
authors.initials: TJ
authors.lastname: Albert
authors.forename: Matthew J
authors.initials: MJ
authors.lastname: Rodesch
authors.forename: David G
authors.initials: DG
authors.lastname: Clayton
authors.forename: John A
authors.initials: JA
authors.lastname: Todd
authors.forename: David A
authors.initials: DA
authors.lastname: van Heel
authors.forename: Vincent
authors.initials: V
authors.lastname: Plagnol
mesh_terms: D000483:Alleles; D022981:Allelic Imbalance / Q000235:genetics*; D020029:Base Pairing / Q000235:genetics; D015982:Bias; D002478:Cells, Cultured; D019295:Computational Biology; D004194:Disease / Q000235:genetics; D044127:Epigenesis, Genetic; D005189:False Positive Reactions; D020869:Gene Expression Profiling / Q000379:methods*; D056426:Genetic Loci / Q000235:genetics; D055106:Genome-Wide Association Study; D006579:Heterozygote; D057166:High-Throughput Screening Assays / Q000379:methods*; D006801:Humans; D020641:Polymorphism, Single Nucleotide / Q000235:genetics; D012333:RNA, Messenger / Q000235:genetics / Q000378:metabolism; D015203:Reproducibility of Results; D017422:Sequence Analysis, DNA / Q000379:methods*
chemical_list: D012333:RNA, Messenger