Home LiteratureArticle Details
PMID: 21637794 Published · ppublish English Journal Article Research Support, N.I.H., Extramural Research Support, Non-U.S. Gov't

Identification, replication, and functional fine-mapping of expression quantitative trait loci in primary human liver tissue.

PLoS genetics ·Vol. 7 ·No. 5 ·2011-05-00 ·Pages e1002078

Innocenti F, Cooper GM, Stanaway IB, Gamazon ER, Smith JD, Mirkov S, Ramirez J, Liu W, Lin YS, Moloney C, Aldred SF, Trinklein ND, Schuetz E, Nickerson DA, Thummel KE, Rieder MJ, Rettie AE, Ratain MJ, Cox NJ, Brown CD

Abstract

The discovery of expression quantitative trait loci ("eQTLs") can help to unravel genetic contributions to complex traits. We identified genetic determinants of human liver gene expression variation using two independent collections of primary tissue profiled with Agilent (n = 206) and Illumina (n = 60) expression arrays and Illumina SNP genotyping (550K), and we also incorporated data from a published study (n = 266). We found that ∼30% of SNP-expression correlations in one study failed to replicate in either of the others, even at thresholds yielding high reproducibility in simulations, and we quantified numerous factors affecting reproducibility. Our data suggest that drug exposure, clinical descriptors, and unknown factors associated with tissue ascertainment and analysis have substantial effects on gene expression and that controlling for hidden confounding variables significantly increases replication rate. Furthermore, we found that reproducible eQTL SNPs were heavily enriched near gene starts and ends, and subsequently resequenced the promoters and 3'UTRs for 14 genes and tested the identified haplotypes using luciferase assays. For three genes, significant haplotype-specific in vitro functional differences correlated directly with expression levels, suggesting that many bona fide eQTLs result from functional variants that can be mechanistically isolated in a high-throughput fashion. Finally, given our study design, we were able to discover and validate hundreds of liver eQTLs. Many of these relate directly to complex traits for which liver-specific analyses are likely to be relevant, and we identified dozens of potential connections with disease-associated loci. These included previously characterized eQTL contributors to diabetes, drug response, and lipid levels, and they suggest novel candidates such as a role for NOD2 expression in leprosy risk and C2orf43 in prostate cancer. In general, the work presented here will be valuable for future efforts to precisely identify and functionally characterize genetic contributions to a variety of complex traits.

MeSH Terms
3' Untranslated Regions Age Factors Blacks Chromosome Mapping Cloning, Molecular Female Gene Expression Profiling Genetic Vectors Genome, Human Genome-Wide Association Study Genotype Hep G2 Cells Hispanic or Latino Humans Liver/metabolism Luminescent Measurements Male Polymorphism, Single Nucleotide Principal Component Analysis Promoter Regions, Genetic Quantitative Trait Loci Reproducibility of Results Sex Factors Transfection Whites
Chemicals
3' Untranslated Regions
Authors & Affiliations
20 authors, click to expand affiliations / ORCID
Innocenti Federico
Cancer Research Center, Committee on Clinical Pharmacology and Pharmacogenomics, Department of Medicine, The University of Chicago, Chicago, Illinois, United States of America.
Cooper Gregory M
Stanaway Ian B
Gamazon Eric R
Smith Joshua D
Mirkov Snezana
Ramirez Jacqueline
Liu Wanqing
Lin Yvonne S
Moloney Cliona
Aldred Shelly Force
Trinklein Nathan D
Schuetz Erin
Nickerson Deborah A
Thummel Ken E
Rieder Mark J
Rettie Allan E
Ratain Mark J
Cox Nancy J
Brown Christopher D
Conflict of Interest

The authors ND Trinklein and SF Aldred are employees of and hold stock in SwitchGear Genomics. SwitchGear Genomics sells reporter vectors commercially and is the owner by assignment of patents or patent applications related to the reporter platform. The author C Moloney is an employee of Merck Research Labs.

References (59)
59 references, click to expand
  1. Understanding mechanisms underlying human gene expression variation with RNA sequencing.
    Nature. 2010 Apr 1;464(7289):768-72 PMID: 20220758
  2. Gene-expression variation within and among human populations.
    Am J Hum Genet. 2007 Mar;80(3):502-9 PMID: 17273971
  3. Polymorphisms within the C-reactive protein (CRP) promoter region are associated with plasma CRP levels.
    Am J Hum Genet. 2005 Jul;77(1):64-77 PMID: 15897982
  4. Gene expression in skin and lymphoblastoid cells: Refined statistical method reveals extensive overlap in cis-eQTL signals.
    Am J Hum Genet. 2010 Dec 10;87(6):779-89 PMID: 21129726
  5. SCAN: SNP and copy number annotation.
    Bioinformatics. 2010 Jan 15;26(2):259-62 PMID: 19933162
  6. Transcriptome genetics using second generation sequencing in a Caucasian population.
    Nature. 2010 Apr 1;464(7289):773-7 PMID: 20220756
  7. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase.
    Am J Hum Genet. 2006 Apr;78(4):629-44 PMID: 16532393
  8. The pursuit of genome-wide association studies: where are we now?
    J Hum Genet. 2010 Apr;55(4):195-206 PMID: 20300123
  9. Genetic inheritance of gene expression in human cell lines.
    Am J Hum Genet. 2004 Dec;75(6):1094-105 PMID: 15514893
  10. Predicting in vivo drug interactions from in vitro drug discovery data.
    Nat Rev Drug Discov. 2005 Oct;4(10):825-33 PMID: 16224454
  11. Imputation-based analysis of association studies: candidate regions and quantitative traits.
    PLoS Genet. 2007 Jul;3(7):e114 PMID: 17676998
  12. Targeted capture and massively parallel sequencing of 12 human exomes.
    Nature. 2009 Sep 10;461(7261):272-6 PMID: 19684571
  13. Genetic dissection of transcriptional regulation in budding yeast.
    Science. 2002 Apr 26;296(5568):752-5 PMID: 11923494
  14. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.
    Nature. 2007 Jun 14;447(7146):799-816 PMID: 17571346
  15. High-resolution mapping of expression-QTLs yields insight into human gene regulation.
    PLoS Genet. 2008 Oct;4(10):e1000214 PMID: 18846210
  16. Specific transcription and RNA splicing defects in five cloned beta-thalassaemia genes.
    Nature. 1983 Apr 14;302(5909):591-6 PMID: 6188062
  17. Identification and functional analysis of human transcriptional promoters.
    Genome Res. 2003 Feb;13(2):308-12 PMID: 12566409
  18. The effect of single nucleotide polymorphisms from genome wide association studies in multiple sclerosis on gene expression.
    PLoS One. 2010 Apr 13;5(4):e10142 PMID: 20405052
  19. Genetic analysis of genome-wide variation in human gene expression.
    Nature. 2004 Aug 12;430(7001):743-7 PMID: 15269782
  20. The architecture of gene regulatory variation across multiple human tissues: the MuTHER study.
    PLoS Genet. 2011 Feb 03;7(2):e1002003 PMID: 21304890
  21. On the design and analysis of gene expression studies in human populations.
    Nat Genet. 2007 Jul;39(7):807-8; author reply 808-9 PMID: 17597765
  22. Genomewide association study of a rapid progression cohort identifies new susceptibility alleles for AIDS (ANRS Genomewide Association Study 03).
    J Infect Dis. 2009 Oct 15;200(8):1194-201 PMID: 19754311
  23. Genome-wide association study identifies five new susceptibility loci for prostate cancer in the Japanese population.
    Nat Genet. 2010 Sep;42(9):751-4 PMID: 20676098
  24. A Bayesian framework to account for complex non-genetic factors in gene expression levels greatly increases power in eQTL studies.
    PLoS Comput Biol. 2010 May 06;6(5):e1000770 PMID: 20463871
  25. Common genetic variants account for differences in gene expression among ethnic groups.
    Nat Genet. 2007 Feb;39(2):226-31 PMID: 17206142
  26. Expression quantitative trait loci: replication, tissue- and sex-specificity in mice.
    Genetics. 2010 Jul;185(3):1059-68 PMID: 20439777
  27. Mapping determinants of human gene expression by regional and genome-wide association.
    Nature. 2005 Oct 27;437(7063):1365-9 PMID: 16251966
  28. Genes mirror geography within Europe.
    Nature. 2008 Nov 6;456(7218):98-101 PMID: 18758442
  29. Practical issues in imputation-based association mapping.
    PLoS Genet. 2008 Dec;4(12):e1000279 PMID: 19057666
  30. A comparison of background correction methods for two-colour microarrays.
    Bioinformatics. 2007 Oct 15;23(20):2700-7 PMID: 17720982
  31. Population genomics of human gene expression.
    Nat Genet. 2007 Oct;39(10):1217-24 PMID: 17873874
  32. High-resolution analysis of DNA regulatory elements by synthetic saturation mutagenesis.
    Nat Biotechnol. 2009 Dec;27(12):1173-5 PMID: 19915551
  33. Pharmacokinetic and pharmacodynamic variability: a daunting challenge in drug therapy.
    Curr Drug Metab. 2007 Feb;8(2):109-36 PMID: 17305491
  34. Sequence polymorphisms cause many false cis eQTLs.
    PLoS One. 2007 Jul 18;2(7):e622 PMID: 17637838
  35. Effect of VKORC1 haplotypes on transcriptional regulation and warfarin dose.
    N Engl J Med. 2005 Jun 2;352(22):2285-93 PMID: 15930419
  36. Genomewide association study of leprosy.
    N Engl J Med. 2009 Dec 31;361(27):2609-18 PMID: 20018961
  37. A genome-wide gene expression signature of environmental geography in leukocytes of Moroccan Amazighs.
    PLoS Genet. 2008 Apr 11;4(4):e1000052 PMID: 18404217
  38. Candidate causal regulatory effects by integration of expression QTLs with complex trait genetic associations.
    PLoS Genet. 2010 Apr 01;6(4):e1000895 PMID: 20369022
  39. Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function.
    Nat Genet. 2005 Mar;37(3):233-42 PMID: 15711545
  40. PLINK: a tool set for whole-genome association and population-based linkage analyses.
    Am J Hum Genet. 2007 Sep;81(3):559-75 PMID: 17701901
  41. Large upward bias in estimation of locus-specific effects from genomewide scans.
    Am J Hum Genet. 2001 Dec;69(6):1357-69 PMID: 11593451
  42. Expression quantitative trait loci are highly sensitive to cellular differentiation state.
    PLoS Genet. 2009 Oct;5(10):e1000692 PMID: 19834560
  43. How replicable are mRNA expression QTL?
    Mamm Genome. 2006 Jun;17(6):643-56 PMID: 16783644
  44. Sample matching by inferred agonal stress in gene expression analyses of the brain.
    BMC Genomics. 2007 Sep 24;8:336 PMID: 17892578
  45. Principal components analysis corrects for stratification in genome-wide association studies.
    Nat Genet. 2006 Aug;38(8):904-9 PMID: 16862161
  46. Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans.
    Nat Genet. 2008 Feb;40(2):189-97 PMID: 18193044
  47. Overcoming the winner's curse: estimating penetrance parameters from case-control data.
    Am J Hum Genet. 2007 Apr;80(4):605-15 PMID: 17357068
  48. Capturing heterogeneity in gene expression studies by surrogate variable analysis.
    PLoS Genet. 2007 Sep;3(9):1724-35 PMID: 17907809
  49. Genetic analysis of human traits in vitro: drug response and gene expression in lymphoblastoid cell lines.
    PLoS Genet. 2008 Nov;4(11):e1000287 PMID: 19043577
  50. From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus.
    Nature. 2010 Aug 5;466(7307):714-9 PMID: 20686566
  51. Genetics of gene expression surveyed in maize, mouse and man.
    Nature. 2003 Mar 20;422(6929):297-302 PMID: 12646919
  52. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.
    Bioinformatics. 2003 Jan 22;19(2):185-93 PMID: 12538238
  53. Common regulatory variation impacts gene expression in a cell type-dependent manner.
    Science. 2009 Sep 4;325(5945):1246-50 PMID: 19644074
  54. A new multipoint method for genome-wide association studies by imputation of genotypes.
    Nat Genet. 2007 Jul;39(7):906-13 PMID: 17572673
  55. Mapping the genetic architecture of gene expression in human liver.
    PLoS Biol. 2008 May 6;6(5):e107 PMID: 18462017
  56. Importin 8 is a gene silencing factor that targets argonaute proteins to distinct mRNAs.
    Cell. 2009 Feb 6;136(3):496-507 PMID: 19167051
  57. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS.
    PLoS Genet. 2010 Apr 01;6(4):e1000888 PMID: 20369019
  58. Identification of common genetic variation that modulates alternative splicing.
    PLoS Genet. 2007 Jun;3(6):e99 PMID: 17571926
  59. Genome-wide analysis of transcript isoform variation in humans.
    Nat Genet. 2008 Feb;40(2):225-31 PMID: 18193047
Article Info
Journal
PLoS genetics
Abbr.
PLoS Genet
ISSN
1553-7404
Published
2011-05-00
Epub
2011-00-26
Pages
e1002078
Language
English
Region
United States
NLM ID
101239074
PMCID
PMC3102751
Subset
IM
Grants
NIDDK NIH HHS · N01-DK-7-0004/HHSN267200700004C · United States
NIDDK NIH HHS · R21 DK081157 · United States
NIGMS NIH HHS · P01 GM32165 · United States
NHLBI NIH HHS · U01 HL66682 · United States
NINDS NIH HHS · R01 NS053646 · United States
NINDS NIH HHS · NS053646 · United States
NHLBI NIH HHS · U01 HL066682 · United States
NCRR NIH HHS · 1KL2RR025015 · United States
NIGMS NIH HHS · P01 GM032165 · United States
NIGMS NIH HHS · U01 GM061393 · United States
NCI NIH HHS · P30 CA014599 · United States
NIGMS NIH HHS · R01GM094418 · United States
NCI NIH HHS · K07 CA140390 · United States
NCI NIH HHS · K07CA140390-01 · United States
NIDDK NIH HHS · HHSN267200700004C · United States
NLM NIH HHS · HHSN267200700004G · United States
NIDDK NIH HHS · R21DK081157-01A2 · United States
NIGMS NIH HHS · U01 GM61393 · United States
NCRR NIH HHS · KL2 RR025015 · United States
NIGMS NIH HHS · R01 GM094418 · United States
Analysis Services
Analysis Services

Contact

No. 2 Wenbo Road, Zhangqiu District, Jinan, Shandong

Qilu Normal University · Genelibs Bioinformatics Lab

750 Shunhua Rd, Jinan

2F, Bldg F, University Science Park

Tel: 0531-88819269

WeChat Official Account

Follow our WeChat subscription account for real-time updates and the latest in medical and biological research.


Business Email

E-mail: product@genelibs.com