Heritable and non-genetic factors as variables of pharmacologic phenotypes in lymphoblastoid cell lines

Single-cell RNA-seq reveals transcriptomic heterogeneity mediated by host-pathogen dynamics in lymphoblastoid cell lines

Lymphoblastoid cell lines (LCLs) are generated by transforming primary B cells with Epstein–Barr virus (EBV) and are used extensively as model systems in viral oncology, immunology, and human genetics research. In this study, we characterized single-cell transcriptomic profiles of five LCLs and present a simple discrete-time simulation to explore the influence of stochasticity on LCL clonal evolution. Single-cell RNA sequencing (scRNA-seq) revealed substantial phenotypic heterogeneity within and across LCLs with respect to immunoglobulin isotype; virus-modulated host pathways involved in survival, activation, and differentiation; viral replication state; and oxidative stress. This heterogeneity is likely attributable to intrinsic variance in primary B cells and host–pathogen dynamics. Stochastic simulations demonstrate that initial primary cell heterogeneity, random sampling, time in culture, and even mild differences in phenotype-specific fitness can contribute substantially to dynamic diversity in populations of nominally clonal cells.

Extensive epigenetic and transcriptomic variability between genetically identical human B-lymphoblastoid cells with implications in pharmacogenomics research

Genotyped human B-lymphoblastoid cell lines (LCLs) are widely used models in mapping quantitative trait loci for chromatin features, gene expression, and drug response. The extent of genotype-independent functional genomic variability of the LCL model, although largely overlooked, may inform association study design. In this study, we use flow cytometry, chromatin immunoprecipitation sequencing and mRNA sequencing to study surface marker patterns, quantify genome-wide chromatin changes (H3K27ac) and transcriptome variability, respectively, among five isogenic LCLs derived from the same individual. Most of the studied LCLs were non-monoclonal and had mature B cell phenotypes. Strikingly, nearly one-fourth of active gene regulatory regions showed significantly variable H3K27ac levels, especially enhancers, among which several were classified as clustered enhancers. Large, contiguous genomic regions showed signs of coordinated activity change. Regulatory differences were mirrored by mRNA expression changes, preferentially affecting hundreds of genes involved in specialized cellular processes including immune and drug response pathways. Differential expression of DPYD, an enzyme involved in 5-fluorouracil (5-FU) catabolism, was associated with variable LCL growth inhibition mediated by 5-FU. The extent of genotype-independent functional genomic variability might highlight the need to revisit study design strategies for LCLs in pharmacogenomics.

The influence of Neanderthal alleles on cytotoxic response

Various studies have shown that people of Eurasian origin contain traces of DNA inherited from interbreeding with Neanderthals. Recent studies have demonstrated that these Neanderthal variants influence a range of clinically important traits and diseases. Thus, understanding the genetic factors responsible for the variability in individual response to drug or chemical exposure is a key goal of pharmacogenomics and toxicogenomics, as dose responses are clinically and epidemiologically important traits. It is well established that ethnic and racial differences are important in dose response traits, but to our knowledge the influence of Neanderthal ancestry on response to xenobiotics is unknown. Towards this aim, we examined if Neanderthal ancestry plays a role in cytotoxic response to anti-cancer drugs and toxic environmental chemicals. We identified common Neanderthal variants in lymphoblastoid cell lines (LCLs) derived from the globally diverse 1000 Genomes Project and Caucasian cell lines from the Children's Hospital of Oakland Research Institute. We analyzed the effects of these Neanderthal alleles on cytotoxic response to 29 anti-cancer drugs and 179 environmental chemicals at varying concentrations using genome-wide data. We identified and replicated single nucleotide polymorphisms (SNPs) from these association results, including a SNP in the SNORD-113 cluster. Our results also show that the Neanderthal alleles cumulatively lead to increased sensitivity to both the anti-cancer drugs and the environmental chemicals. Our results demonstrate the influence of Neanderthal ancestry-informative markers on cytotoxic response. These results could be important in identifying biomarkers for personalized medicine or in dissecting the underlying etiology of dose response traits.

Pharmacogenetics of Chemotherapy-Induced Cardiotoxicity

PURPOSE OF REVIEW

The goal of this review is to summarize current understanding of pharmacogenetics and pharmacogenomics in chemotherapy-induced cardiotoxicity.

RECENT FINDINGS

Most of the studies rely on in vitro cytotoxic assays. There have been several smaller scale candidate gene approaches and a handful of genome-wide studies linking genetic variation to susceptibility to chemotherapy-induced cardiotoxicity. Currently, pharmacogenomic testing of all childhood cancer patients with an indication for doxorubicin or daunorubicin therapy for RARG rs2229774, SLC28A3 rs7853758, and UGT1A6*4 rs17863783 variants is recommended. There is no recommendation regarding testing in adults. There is clear evidence pointing to the role of pharmacogenetics and pharmacogenomics in cardiotoxicity susceptibility to chemotherapeutic agents. Larger scale studies are needed to further identify susceptibility markers and to develop pharmacogenomics-based risk profiling to improve quality of life and life expectancy in cancer survivors.

Human iPSC-derived neurons and lymphoblastoid cells for personalized medicine research in neuropsychiatric disorders

The development and clinical implementation of personalized medicine crucially depends on the availability of high-quality human biosamples; animal models, although capable of modeling complex human diseases, cannot reflect the large variation in the human genome, epigenome, transcriptome, proteome, and metabolome. Although the biosamples available from public biobanks that store human tissues and cells may represent the large human diversity for most diseases, these samples are not always sufficient for developing biomarkers for patient-tailored therapies for neuropsychiatric disorders. Postmortem human tissues are available from many biobanks; nevertheless, collections of neuronal human cells from large patient cohorts representing the human diversity remain scarce. Two tools are gaining popularity for personalized medicine research on neuropsychiatric disorders: human induced pluripotent stem cell-derived neurons and human lymphoblastoid cell lines. This review examines and contrasts the advantages and limitations of each tool for personalized medicine research.

Comparable impediment of cognitive function in female and male rats subsequent to daily administration of haloperidol after traumatic brain injury

Antipsychotic drugs, such as haloperidol (HAL), are prescribed in the clinic to manage traumatic brain injury (TBI)-induced agitation. While preclinical studies have consistently shown that once-daily administration of HAL hinders functional recovery after TBI in male rats, its effects in females are unknown. Hence, the objective of this study was to directly compare neurobehavioral and histological outcomes in both sexes to determine whether the reported deleterious effects of HAL extend to females. Anesthetized adult female and male rats received either a controlled cortical impact (CCI) or sham injury and then were randomly assigned to a dosing regimen of HAL (0.5mg/kg, i.p.) or vehicle (VEH; 1mL/kg, i.p.) that was initiated 24h after injury and continued once daily for 19 consecutive days. Motor function was tested using established beam-balance/walk protocols on post-operative days 1-5 and acquisition of spatial learning was assessed with a well-validated Morris water maze task on days 14-19. Cortical lesion volume was quantified at 21days. No statistical differences were revealed between the HAL and VEH-treated sham groups and thus they were pooled for each sex. HAL only impaired motor recovery in males (p<0.05), but significantly diminished spatial learning in both sexes (p<0.05). Females, regardless of treatment, exhibited smaller cortical lesions vs VEH-treated males (p<0.05). Taken together, the data show that daily HAL does not prohibit motor recovery in females, but does negatively impact cognition. These task-dependent differential effects of HAL in female vs male rats may have clinical significance as they can direct therapy.

MicroRNA biogenesis and cellular proliferation

Given the fundamental roles of microRNAs (miRNAs) in physiological, developmental, and pathologic processes, we hypothesized that genes involved in miRNA biogenesis contribute to human complex traits. For 13 such genes, we evaluated the relationship between transcription and 2 classes of complex traits, namely cellular growth and sensitivity to various chemotherapeutic agents in a set of lymphoblastoid cell lines. We found a highly significant correlation between argonaute RNA-induced silencing complex catalytic component 2 (AGO2) expression and cellular growth rate (Bonferroni-adjusted P < 0.05), and report additional miRNA biogenesis genes with suggestive associations with either cellular growth rate or chemotherapeutic sensitivity. AGO2 expression was found to be correlated with multiple drug sensitivity phenotypes. Furthermore, small interfering RNA knockdown of AGO2 resulted in cellular growth inhibition in an ovarian cancer cell line (OVCAR-3), supporting the role of this miRNA biogenesis gene in cell proliferation in cancer cells. Expression quantitative trait loci mapping indicated that genetic variation (in the form of both single-nucleotide polymorphisms and copy number variations) that may regulate the expression of AGO2 can have downstream effects on cellular growth-dependent complex phenotypes.

Population-based in vitro hazard and concentration-response assessment of chemicals: the 1000 genomes high-throughput screening study

BACKGROUND

Understanding of human variation in toxicity to environmental chemicals remains limited, so human health risk assessments still largely rely on a generic 10-fold factor (10½ each for toxicokinetics and toxicodynamics) to account for sensitive individuals or subpopulations.

OBJECTIVES

We tested a hypothesis that population-wide in vitro cytotoxicity screening can rapidly inform both the magnitude of and molecular causes for interindividual toxicodynamic variability.

METHODS

We used 1,086 lymphoblastoid cell lines from the 1000 Genomes Project, representing nine populations from five continents, to assess variation in cytotoxic response to 179 chemicals. Analysis included assessments of population variation and heritability, and genome-wide association mapping, with attention to phenotypic relevance to human exposures.

RESULTS

For about half the tested compounds, cytotoxic response in the 1% most "sensitive" individual occurred at concentrations within a factor of 10½ (i.e., approximately 3) of that in the median individual; however, for some compounds, this factor was > 10. Genetic mapping suggested important roles for variation in membrane and transmembrane genes, with a number of chemicals showing association with SNP rs13120371 in the solute carrier SLC7A11, previously implicated in chemoresistance.

CONCLUSIONS

This experimental approach fills critical gaps unaddressed by recent large-scale toxicity testing programs, providing quantitative, experimentally based estimates of human toxicodynamic variability, and also testable hypotheses about mechanisms contributing to interindividual variation.

Reprogramming LCLs to iPSCs Results in Recovery of Donor-Specific Gene Expression Signature

Renewable in vitro cell cultures, such as lymphoblastoid cell lines (LCLs), have facilitated studies that contributed to our understanding of genetic influence on human traits. However, the degree to which cell lines faithfully maintain differences in donor-specific phenotypes is still debated. We have previously reported that standard cell line maintenance practice results in a loss of donor-specific gene expression signatures in LCLs. An alternative to the LCL model is the induced pluripotent stem cell (iPSC) system, which carries the potential to model tissue-specific physiology through the use of differentiation protocols. Still, existing LCL banks represent an important source of starting material for iPSC generation, and it is possible that the disruptions in gene regulation associated with long-term LCL maintenance could persist through the reprogramming process. To address this concern, we studied the effect of reprogramming mature LCL cultures from six unrelated donors to iPSCs on the ensuing gene expression patterns within and between individuals. We show that the reprogramming process results in a recovery of donor-specific gene regulatory signatures, increasing the number of genes with a detectable donor effect by an order of magnitude. The proportion of variation in gene expression statistically attributed to donor increases from 6.9% in LCLs to 24.5% in iPSCs (P < 10-15). Since environmental contributions are unlikely to be a source of individual variation in our system of highly passaged cultured cell lines, our observations suggest that the effect of genotype on gene regulation is more pronounced in iPSCs than in LCLs. Our findings indicate that iPSCs can be a powerful model system for studies of phenotypic variation across individuals in general, and the genetic association with variation in gene regulation in particular. We further conclude that LCLs are an appropriate starting material for iPSC generation.

Molecular insight into thiopurine resistance: transcriptomic signature in lymphoblastoid cell lines

Background

There has been considerable progress in the management of acute lymphoblastic leukemia (ALL) but further improvement is needed to increase long-term survival. The thiopurine agent 6-mercaptopurine (6-MP) used for ALL maintenance therapy has a key influence on clinical outcomes and relapse prevention. Genetic inheritance in thiopurine metabolism plays a major role in interindividual clinical response variability to thiopurines; however, most cases of thiopurine resistance remain unexplained.

Methods

We used lymphoblastoid cell lines (LCLs) from healthy donors, selected for their extreme thiopurine susceptibility. Thiopurine metabolism was characterized by the determination of TPMT and HPRT activity. We performed genome-wide expression profiling in resistant and sensitive cell lines with the goal of elucidating the mechanisms of thiopurine resistance.

Results

We determined a higher TPMT activity (+44%; P = 0.024) in resistant compared to sensitive cell lines, although there was no difference in HPRT activity. We identified a 32-gene transcriptomic signature that predicts thiopurine resistance. This signature includes the GTPBP4 gene coding for a GTP-binding protein that interacts with p53. A comprehensive pathway analysis of the genes differentially expressed between resistant and sensitive cell lines indicated a role for cell cycle and DNA mismatch repair system in thiopurine resistance. It also revealed overexpression of the ATM/p53/p21 pathway, which is activated in response to DNA damage and induces cell cycle arrest in thiopurine resistant LCLs. Furthermore, overexpression of the p53 target gene TNFRSF10D or the negative cell cycle regulator CCNG2 induces cell cycle arrest and may also contribute to thiopurine resistance. ARHGDIA under-expression in resistant cell lines may constitute a novel molecular mechanism contributing to thiopurine resistance based on Rac1 inhibition induced apoptosis and in relation with thiopurine pharmacodynamics.

Conclusion

Our study provides new insights into the molecular mechanisms underlying thiopurine resistance and suggests a potential research focus for developing tailored medicine.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0150-6) contains supplementary material, which is available to authorized users.

Genome-wide association and pharmacological profiling of 29 anticancer agents using lymphoblastoid cell lines

AIM

Association mapping with lymphoblastoid cell lines (LCLs) is a promising approach in pharmacogenomics research, and in the current study we utilized LCLs to perform association mapping for 29 chemotherapy drugs.

MATERIALS & METHODS

Currently, we use LCLs to perform genome-wide association mapping of the cytotoxic response of 520 European-Americans to 29 different anticancer drugs; the largest LCL study to date. A novel association approach using a multivariate analysis of covariance design was employed with the software program MAGWAS, testing for differences in the dose-response profiles between genotypes without making assumptions about the response curve or the biologic mode of association. Additionally, by classifying 25 of the 29 drugs into eight families according to structural and mechanistic relationships, MAGWAS was used to test for associations that were shared across each drug family. Finally, a unique algorithm using multivariate responses and multiple linear regressions across pairs of response curves was used for unsupervised clustering of drugs.

RESULTS

Among the single-drug studies, suggestive associations were obtained for 18 loci, 12 within/near genes. Three of these, MED12L, CHN2 and MGMT, have been previously implicated in cancer pharmacogenomics. The drug family associations resulted in four additional suggestive loci (three contained within/near genes). One of these genes, HDAC4, associated with the DNA alkylating agents, shows possible clinical interactions with temozolomide. For the drug clustering analysis, 18 of 25 drugs clustered into the appropriate family.

CONCLUSION

This study demonstrates the utility of LCLs in identifying genes that have clinical importance in drug response and for assigning unclassified agents to specific drug families, and proposes new candidate genes for follow-up in a large number of chemotherapy drugs.

Genomic architecture of pharmacological efficacy and adverse events

The pharmacokinetic and pharmacodynamic disciplines address pharmacological traits, including efficacy and adverse events. Pharmacogenomics studies have identified pervasive genetic effects on treatment outcomes, resulting in the development of genetic biomarkers for optimization of drug therapy. Pharmacogenomics-based tests are already being applied in clinical decision making. However, despite substantial progress in identifying the genetic etiology of pharmacological response, current biomarker panels still largely rely on single gene tests with a large portion of the genetic effects remaining to be discovered. Future research must account for the combined effects of multiple genetic variants, incorporate pathway-based approaches, explore gene-gene interactions and nonprotein coding functional genetic variants, extend studies across ancestral populations, and prioritize laboratory characterization of molecular mechanisms. Because genetic factors can play a key role in drug response, accurate biomarker tests capturing the main genetic factors determining treatment outcomes have substantial potential for improving individual clinical care.

Functional analysis and molecular characterization of spontaneously outgrown human lymphoblastoid cell lines

In vitro, the infection of human B-cells with the lymphotropic gammaherpesvirus Epstein-Barr virus (EBV) induces formation of permanently growing lymphoblastoid cell lines (LCL). In a spontaneously outgrown LCL (cell line CSIII), we detected nucleotide sequence variations of the EBV nuclear antigen 1 (EBNA1) RNA that was different from the reference sequence of EBNA1 in the prototypic EBV strain B95-8. In the present study, we molecularly and functionally characterized this virus isolate in comparison to LCL with the prototypic nucleotide sequence. Although we detected high functional similarity between CSIII and the other LCL, our data suggest that the lytic cycle might be ineffective in the CSIII LCL. DNA microarray analysis indicated that RNA binding motif, single stranded interacting protein 1 (RBMS1), which is typically expressed in latency III of EBV to prevent the lytic cycle, was the most overexpressed gene in CSIII LCL.

The impact of microRNA expression on cellular proliferation

As an important class of non-coding regulatory RNAs, microRNAs (miRNAs) play a key role in a range of biological processes. These molecules serve as post-transcriptional regulators of gene expression and their regulatory activity has been implicated in disease pathophysiology and pharmacological traits. We sought to investigate the impact of miRNAs on cellular proliferation to gain insight into the molecular basis of complex traits that depend on cellular growth, including, most prominently, cancer. We examined the relationship between miRNA expression and intrinsic cellular growth (iGrowth) in the HapMap lymphoblastoid cell lines derived from individuals of different ethnic backgrounds. We found a substantial enrichment for miRNAs (53 miRNAs, FDR < 0.05) correlated with cellular proliferation in pooled CEU (Caucasian of northern and western European descent) and YRI (individuals from Ibadan, Nigeria) samples. Specifically, 119 miRNAs (59 %) were significantly correlated with iGrowth in YRI; of these miRNAs, 18 were correlated with iGrowth in CEU. To gain further insight into the effect of miRNAs on cellular proliferation in cancer, we showed that over-expression of miR-22, one of the top iGrowth-associated miRNAs, leads to growth inhibition in an ovarian cancer cell line (SKOV3). Furthermore, over-expression of miR-22 down-regulates the expression of its target genes (MXI1 and SLC25A37) in this ovarian cancer cell line, highlighting an miRNA-mediated regulatory network potentially important for cellular proliferation. Importantly, our study identified miRNAs that can be used as molecular targets in cancer therapy.

Multiplexed DNA repair assays for multiple lesions and multiple doses via transcription inhibition and transcriptional mutagenesis

The capacity to repair different types of DNA damage varies among individuals, making them more or less susceptible to the detrimental health consequences of damage exposures. Current methods for measuring DNA repair capacity (DRC) are relatively labor intensive, often indirect, and usually limited to a single repair pathway. Here, we describe a fluorescence-based multiplex flow-cytometric host cell reactivation assay (FM-HCR) that measures the ability of human cells to repair plasmid reporters, each bearing a different type of DNA damage or different doses of the same type of DNA damage. FM-HCR simultaneously measures repair capacity in any four of the following pathways: nucleotide excision repair, mismatch repair, base excision repair, nonhomologous end joining, homologous recombination, and methylguanine methyltransferase. We show that FM-HCR can measure interindividual DRC differences in a panel of 24 cell lines derived from genetically diverse, apparently healthy individuals, and we show that FM-HCR may be used to identify inhibitors or enhancers of DRC. We further develop a next-generation sequencing-based HCR assay (HCR-Seq) that detects rare transcriptional mutagenesis events due to lesion bypass by RNA polymerase, providing an added dimension to DRC measurements. FM-HCR and HCR-Seq provide powerful tools for exploring relationships among global DRC, disease susceptibility, and optimal treatment.

Genetic variation is the major determinant of individual differences in leukocyte endothelial adhesion

Objective

To determine the genetic contribution to leukocyte endothelial adhesion.

Methods

Leukocyte endothelial adhesion was assessed through a novel cell-based assay using human lymphoblastoid cell lines. A high-throughput screening method was developed to evaluate the inter-individual variability in leukocyte endothelial adhesion using lymphoblastoid cell lines derived from different donors. To assess heritability, ninety-two lymphoblastoid cell lines derived from twenty-three monozygotic twin pairs and twenty-three sibling pairs were compared. These lymphoblastoid cell lines were plated with the endothelial cell line EA.hy926 and labeled with Calcein AM dye. Fluorescence was assessed to determine endothelial cell adhesion to each lymphoblastoid cell line. Intra-pair similarity was determined for monozygotic twins and siblings using Pearson pairwise correlation coefficients.

Results

A leukocyte endothelial adhesion assay for lymphoblastoid cell lines was developed and optimized (CV = 8.68, Z′-factor = 0.67, SNR = 18.41). A higher adhesion correlation was found between the twins than that between the siblings. Intra-pair similarity for leukocyte endothelial adhesion in monozygotic twins was 0.60 compared to 0.25 in the siblings. The extent to which these differences are attributable to underlying genetic factors was quantified and the heritability of leukocyte endothelial adhesion was calculated to be 69.66% (p-value<0.0001).

Conclusions

There is a heritable component to leukocyte endothelial adhesion. Underlying genetic predisposition plays a significant role in inter-individual variability of leukocyte endothelial adhesion.

Comprehensive genetic analysis of cytarabine sensitivity in a cell-based model identifies polymorphisms associated with outcome in AML patients

A whole-genome approach was used to investigate the genetic determinants of cytarabine-induced cytotoxicity. We performed a meta-analysis of genome-wide association studies involving 523 lymphoblastoid cell lines (LCLs) from individuals of European, African, Asian, and African American ancestry. Several of the highest-ranked single-nucleotide polymorphisms (SNPs) were within the mutated in colorectal cancers (MCC) gene. MCC expression was induced by cytarabine treatment from 1.7- to 26.6-fold in LCLs. A total of 33 SNPs ranked at the top of the meta-analysis (P < 10(-5)) were successfully tested in a clinical trial of patients randomized to receive low-dose or high-dose cytarabine plus daunorubicin and etoposide; of these, 18 showed association (P < .05) with either cytarabine 50% inhibitory concentration in leukemia cells or clinical response parameters (minimal residual disease, overall survival (OS), and treatment-related mortality). This count (n = 18) was significantly greater than expected by chance (P = .016). For rs1203633, LCLs with AA genotype were more sensitive to cytarabine-induced cytotoxicity (P = 1.31 × 10(-6)) and AA (vs GA or GG) genotype was associated with poorer OS (P = .015), likely as a result of greater treatment-related mortality (P = .0037) in patients with acute myeloid leukemia (AML). This multicenter AML02 study trial was registered at www.clinicaltrials.gov as #NCT00136084.

Multivariate methods and software for association mapping in dose-response genome-wide association studies

BACKGROUND

The large sample sizes, freedom of ethical restrictions and ease of repeated measurements make cytotoxicity assays of immortalized lymphoblastoid cell lines a powerful new in vitro method in pharmacogenomics research. However, previous studies may have over-simplified the complex differences in dose-response profiles between genotypes, resulting in a loss of power.

METHODS

The current study investigates four previously studied methods, plus one new method based on a multivariate analysis of variance (MANOVA) design. A simulation study was performed using differences in cancer drug response between genotypes for biologically meaningful loci. These loci also showed significance in separate genome-wide association studies. This manuscript builds upon a previous study, where differences in dose-response curves between genotypes were constructed using the hill slope equation.

CONCLUSION

Overall, MANOVA was found to be the most powerful method for detecting real signals, and was also the most robust method for detection using alternatives generated with the previous simulation study. This method is also attractive because test statistics follow their expected distributions under the null hypothesis for both simulated and real data. The success of this method inspired the creation of the software program MAGWAS. MAGWAS is a computationally efficient, user-friendly, open source software tool that works on most platforms and performs GWASs for individuals having multivariate responses using standard file formats.

Genome-wide transcriptomic variations of human lymphoblastoid cell lines: insights from pairwise gene-expression correlations

Aims: Human lymphoblastoid cell lines (LCLs) are a rich resource of information on human interindividual genomic, transcriptomic, proteomic and phenomic variations, and are therefore gaining popularity for pharmacogenomic studies. In the present study we demonstrate that genome-wide transcriptomic data from a small LCL panel from unrelated individuals is sufficient for detecting pairs of genes that exhibit highly correlated expression levels and may thus convey insights about coregulated genes. Materials & methods: RNA samples were prepared from LCLs representing 12 unrelated healthy adult female Caucasian donors. Transcript levels were determined with the Affymetrix Human Gene arrays. Expression-level correlations were searched using Partek® Genomics Suite™ and the R environment. Sequences of detected correlated gene pairs were compared for shared conserved 3´-UTR miRNA binding. Results: Most of the approximately 33,000 transcripts covered by the Affymetrix arrays showed closely similar expression levels in LCLs from unrelated donors. However, the expression levels of some transcripts showed large inter-individual variations. When comparing the expression levels of each of the top 1000 genes showing the largest interindividual expression variations against the others, two sets containing 156 and 4438 correlated gene pairs with false-discovery rates of 0.01 and 0.05 were detected, respectively. Similar analysis of another gene-expression data set from LCLs (GSE11582) indicated that 61 and 39% of identified pairs matched the pairs detected from our transcriptomic data, respectively. Shared conserved 3´-UTR miRNA binding sites were noted for 14–17% of the top 100 gene pairs, suggesting that regulation by miRNA may contribute to their coordinated expression. Conclusion: Probing genome-wide transcriptomic data sets of LCLs from unrelated individuals may detect coregulated genes, adding insights on cellular regulation by miRNAs.

Original submitted 11 July 2012; Revision submitted 4 September 2012

Integration of cell line and clinical trial genome-wide analyses supports a polygenic architecture of Paclitaxel-induced sensory peripheral neuropathy

PURPOSE

We sought to show the relevance of a lymphoblastoid cell line (LCL) model in the discovery of clinically relevant genetic variants affecting chemotherapeutic response by comparing LCL genome-wide association study (GWAS) results to clinical GWAS results.

EXPERIMENTAL DESIGN

A GWAS of paclitaxel-induced cytotoxicity was conducted in 247 LCLs from the HapMap Project and compared with a GWAS of sensory peripheral neuropathy in patients with breast cancer (n = 855) treated with paclitaxel in the Cancer and Leukemia Group B (CALGB) 40101 trial. Significant enrichment was assessed by permutation resampling analysis.

RESULTS

We observed an enrichment of LCL cytotoxicity-associated single-nucleotide polymorphisms (SNP) in the sensory peripheral neuropathy-associated SNPs from the clinical trial with concordant allelic directions of effect (empirical P = 0.007). Of the 24 SNPs that overlap between the clinical trial (P < 0.05) and the preclinical cytotoxicity study (P < 0.001), 19 of them are expression quantitative trait loci (eQTL), which is a significant enrichment of this functional class (empirical P = 0.0447). One of these eQTLs is located in RFX2, which encodes a member of the DNA-binding regulatory factor X family. Decreased expression of this gene by siRNA resulted in increased sensitivity of Neuroscreen-1(NS-1; rat pheochromocytoma) cells to paclitaxel as measured by reduced neurite outgrowth and increased cytotoxicity, functionally validating the involvement of RFX2 in nerve cell response to paclitaxel.

CONCLUSIONS

The enrichment results and functional example imply that cellular models of chemotherapeutic toxicity may capture components of the underlying polygenic architecture of related traits in patients.

Gemcitabine metabolic pathway genetic polymorphisms and response in patients with non-small cell lung cancer

BACKGROUND AND OBJECTIVE

Gemcitabine is widely used to treat non-small cell lung cancer (NSCLC). The aim of this study was to assess the pharmacogenomic effects of the entire gemcitabine metabolic pathway, we genotyped single nucleotide polymorphisms (SNPs) within the 17 pathway genes using DNA samples from patients with NSCLC treated with gemcitabine to determine the effect of genetic variants within gemcitabine pathway genes on overall survival (OS) of patients with NSCLC after treatment of gemcitabine.

METHODS

Eight of the 17 pathway genes were resequenced with DNA samples from Coriell lymphoblastoid cell lines (LCLs) using Sanger sequencing for all exons, exon-intron junctions, and 5'-, 3'-UTRs. A total of 107 tagging SNPs were selected on the basis of the resequencing data for the eight genes and on HapMap data for the remaining nine genes, followed by successful genotyping of 394 NSCLC patient DNA samples. Association of SNPs/haplotypes with OS was performed using the Cox regression model, followed by functional studies performed with LCLs and NSCLC cell lines.

RESULTS

Five SNPs in four genes (CDA, NT5C2, RRM1, and SLC29A1) showed associations with OS of those patients with NSCLC, as well as nine haplotypes in four genes (RRM1, RRM2, SLC28A3, and SLC29A1) with a P value of less than 0.05. Genotype imputation using the LCLs was performed for a region of 200 kb surrounding those SNPs, followed by association studies with gemcitabine cytotoxicity. Functional studies demonstrated that downregulation of SLC29A1, NT5C2, and RRM1 in NSCLC cell lines altered cell susceptibility to gemcitabine.

CONCLUSION

These studies help in identifying biomarkers to predict gemcitabine response in NSCLC, a step toward the individualized chemotherapy of lung cancer.

Sex differences in human lymphoblastoid cells sensitivities to antipsychotic drugs

Adverse drug reactions (ADRs) are a major concern in pharmacotherapy and are more common among women. Immortalized human lymphoblastoid cell lines (LCLs) are emerging as a novel tool for studying interindividual variability in drug response, including ADRs. In the present study, we compared sensitivities of LCLs from unrelated healthy male and female donors to growth inhibition by a panel of common drugs. We observed large interindividual drug sensitivity variations with similar mean sensitivities recorded for LCLs from male and female donors for most tested drugs. A notable exception was observed for the typical antipsychotic haloperidol and the atypical antipsychotic risperidone, which exhibited, on average, more robust in vitro growth inhibition in male as compared with female LCLs. An opposite finding was observed for the antidepressant paroxetine, which was more potent for inhibiting the growth of female as compared with male LCLs. These observations are discussed in the context of the higher incidence of dystonia reported for male schizophrenia patients treated with haloperidol and the higher efficacy of paroxetine in female major depression patients.

Integrating heterogeneous high-throughput data for meta-dimensional pharmacogenomics and disease-related studies

The current paradigm of human genetics research is to analyze variation of a single data type (i.e., DNA sequence or RNA levels) to detect genes and pathways that underlie complex traits such as disease state or drug response. While these studies have detected thousands of variations that associate with hundreds of complex phenotypes, much of the estimated heritability, or trait variability due to genetic factors, remain unexplained. We may be able to account for a portion of the missing heritability if we incorporate a systems biology approach into these analyses. Rapid technological advances will make it possible for scientists to explore this hypothesis via the generation of high-throughput omics data - transcriptomic, proteomic and methylomic to name a few. Analyzing this 'meta-dimensional' data will require clever statistical techniques that allow for the integration of qualitative and quantitative predictor variables. For this article, we examine two major categories of approaches for integrated data analysis, give examples of their use in experimental and in silico datasets, and assess the limitations of each method.

Lymphoblastoid cell lines in pharmacogenomic discovery and clinical translation

The ability to predict how an individual patient will respond to a particular treatment is the ambitious goal of personalized medicine. The genetic make up of an individual has been shown to play a role in drug response. For pharmacogenomic studies, human lymphoblastoid cell lines (LCLs) comprise a useful model system for identifying genetic variants associated with pharmacologic phenotypes. The availability of extensive genotype data for many panels of LCLs derived from individuals of diverse ancestry allows for the study of genetic variants contributing to interethnic and interindividual variation in susceptibility to drugs. Many genome-wide association studies for drug-induced phenotypes have been performed in LCLs, often incorporating gene-expression data. LCLs are also being used in follow-up studies to clinical findings to determine how an associated variant functions to affect phenotype. This review describes the most recent pharmacogenomic findings made in LCLs, including the translation of some findings to clinical cohorts.

Genetic architecture of microRNA expression: implications for the transcriptome and complex traits

We sought to comprehensively and systematically characterize the relationship between genetic variation, miRNA expression, and mRNA expression. Genome-wide expression profiling of samples of European and African ancestry identified in each population hundreds of miRNAs whose increased expression is correlated with correspondingly reduced expression of target mRNAs. We scanned 3' UTR SNPs with a potential functional effect on miRNA binding for cis-acting expression quantitative trait loci (eQTLs) for the corresponding proximal target genes. To extend sequence-based, localized analyses of SNP effect on miRNA binding, we proceeded to dissect the genetic basis of miRNA expression variation; we mapped miRNA expression levels-as quantitative traits-to loci in the genome as miRNA eQTLs, demonstrating that miRNA expression is under significant genetic control. We found that SNPs associated with miRNA expression are significantly enriched with those SNPs already shown to be associated with mRNA. Moreover, we discovered that many of the miRNA-associated genetic variations identified in our study are associated with a broad spectrum of human complex traits from the National Human Genome Research Institute catalog of published genome-wide association studies. Experimentally, we replicated miRNA-induced mRNA expression inhibition and the cis-eQTL relationship to the target gene for several identified relationships among SNPs, miRNAs, and mRNAs in an independent set of samples; furthermore, we conducted miRNA overexpression and inhibition experiments to functionally validate the miRNA-mRNA relationships. This study extends our understanding of the genetic regulation of the transcriptome and suggests that genetic variation might underlie observed relationships between miRNAs and mRNAs more commonly than has previously been appreciated.

Using germline genomics to individualize pediatric cancer treatments

The amazing successes in cure rates for children with cancer over the last century have come in large part from identifying clinical, genetic, and molecular variables associated with response to therapy in large cooperative clinical trials and stratifying therapies according to the predicted risk of relapse. There is an expanding interest in identifying germline genomic variants, as opposed to genetic variants within the tumor, that are associated with susceptibility to toxicity and for risk of relapse. This review highlights the most important germline pharmacogenetic and pharmacogenomic studies in pediatric oncology. Incorporating germline genomics into risk-adapted therapies will likely lead to safer and more effective treatments for children with cancer.

Mixed effects modeling of proliferation rates in cell-based models: consequence for pharmacogenomics and cancer

The International HapMap project has made publicly available extensive genotypic data on a number of lymphoblastoid cell lines (LCLs). Building on this resource, many research groups have generated a large amount of phenotypic data on these cell lines to facilitate genetic studies of disease risk or drug response. However, one problem that may reduce the usefulness of these resources is the biological noise inherent to cellular phenotypes. We developed a novel method, termed Mixed Effects Model Averaging (MEM), which pools data from multiple sources and generates an intrinsic cellular growth rate phenotype. This intrinsic growth rate was estimated for each of over 500 HapMap cell lines. We then examined the association of this intrinsic growth rate with gene expression levels and found that almost 30% (2,967 out of 10,748) of the genes tested were significant with FDR less than 10%. We probed further to demonstrate evidence of a genetic effect on intrinsic growth rate by determining a significant enrichment in growth-associated genes among genes targeted by top growth-associated SNPs (as eQTLs). The estimated intrinsic growth rate as well as the strength of the association with genetic variants and gene expression traits are made publicly available through a cell-based pharmacogenomics database, PACdb. This resource should enable researchers to explore the mediating effects of proliferation rate on other phenotypes.

Cell-based models provide a convenient system to conduct studies that would be impossible to apply to human subjects, but the phenotypes measured on these models can be marred with biological noise. We propose a method (MEM) to address this issue by statistically combining data from various sources, and we apply it to the proliferation rates of cell lines collected as part of the International HapMap project. We show that the proliferation rate computed using our method is a better measure of the true proliferation rate of the cells and produces a much stronger association with gene expression phenotypes on the same cell lines: more than 30% of the genes tested were significantly associated with proliferation rate. We also demonstrate that genetic variants have an effect on growth rate. Finally, we make these intrinsic proliferation rates and the strength of the association with gene expression phenotypes public, which should allow other researchers to explore the mediating effects of proliferation on other phenotypes.

The combination of a genome-wide association study of lymphocyte count and analysis of gene expression data reveals novel asthma candidate genes

Recent genome-wide association studies (GWAS) have identified a number of novel genetic associations with complex human diseases. In spite of these successes, results from GWAS generally explain only a small proportion of disease heritability, an observation termed the ‘missing heritability problem’. Several sources for the missing heritability have been proposed, including the contribution of many common variants with small individual effect sizes, which cannot be reliably found using the standard GWAS approach. The goal of our study was to explore a complimentary approach, which combines GWAS results with functional data in order to identify novel genetic associations with small effect sizes. To do so, we conducted a GWAS for lymphocyte count, a physiologic quantitative trait associated with asthma, in 462 Hutterites. In parallel, we performed a genome-wide gene expression study in lymphoblastoid cell lines from 96 Hutterites. We found significant support for genetic associations using the GWAS data when we considered variants near the 193 genes whose expression levels across individuals were most correlated with lymphocyte counts. Interestingly, these variants are also enriched with signatures of an association with asthma susceptibility, an observation we were able to replicate. The associated loci include genes previously implicated in asthma susceptibility as well as novel candidate genes enriched for functions related to T cell receptor signaling and adenosine triphosphate synthesis. Our results, therefore, establish a new set of asthma susceptibility candidate genes. More generally, our observations support the notion that many loci of small effects influence variation in lymphocyte count and asthma susceptibility.

Identification of novel germline polymorphisms governing capecitabine sensitivity

BACKGROUND

Capecitabine, an oral 5-fluorouracil (5-FU) prodrug, is widely used in the treatment of breast, colorectal, and gastric cancers. To guide the selection of patients with potentially the greatest benefit of experiencing antitumor efficacy, or, alternatively, of developing toxicities, identifying genomic predictors of capecitabine sensitivity could permit its more informed use.

METHODS

The objective of this study was to perform capecitabine sensitivity genome-wide association studies (GWAS) using 503 well genotyped human cell lines from individuals representing multiple different world populations. A meta-analysis that included all ethnic populations then enabled the identification of novel germline determinants (single nucleotide polymorphisms [SNPs]) of capecitabine susceptibility.

RESULTS

First, an intrapopulation GWAS of Caucasian individuals identified reference SNP 4702484 (rs4702484) (within adenylate cyclase 2 [ADCY2]) at a level reaching genome-wide significance (P = 5.2 × 10(-8) ). This SNP is located upstream of the 5 methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene, and it is known that the enzyme for MTRR is involved in the methionine-folate biosynthesis and metabolism pathway, which is the primary target of 5-FU-related compounds, although the authors were unable to identify a direct relation between rs4702484 and MTRR expression in a tested subset of cells. In the meta-analysis, 4 SNPs comprised the top hits, which, again, included rs4702484 and 3 additional SNPs (rs8101143, rs576523, and rs361433) that approached genome-wide significance (P values from 1.9 × 10(-7) to 8.8 × 10(-7) ). The meta-analysis also identified 1 missense variant (rs11722476; serine to asparagine) within switch/sucrose nonfermentable-related, matrix-associated, actin-dependent regulator of chromatin (SMARCAD1), a novel gene for association with capecitabine/5-FU susceptibility.

CONCLUSIONS

Toward the goal of individualizing cancer chemotherapy, the current study identified novel SNPs and genes associated with capecitabine sensitivity that are potentially informative and testable in any patient regardless of ethnicity.

A comparison of association methods for cytotoxicity mapping in pharmacogenomics

Cytotoxicity assays of immortalized lymphoblastoid cell lines (LCLs) represent a promising new in vitro approach in pharmacogenomics research. However, previous studies employing LCLs in gene mapping have used simple association methods, which may not adequately capture the true differences in non-linear response profiles between genotypes. Two common approaches summarize each dose-response curve with either the IC50 or the slope parameter estimates from a hill slope fit and treat these estimates as the response in a linear model. The current study investigates these two methods, as well as four novel methods, and compares their power to detect differences between the response profiles of genotypes under a variety of different alternatives. The four novel methods include two methods that summarize each dose-response by its area under the curve, one method based off of an analysis of variance (ANOVA) design, and one method that compares hill slope fits for all individuals of each genotype. The power of each method was found to depend not only on the choice of alternative, but also on the choice for the set of dosages used in cytotoxicity measurements. The ANOVA-based method was found to be the most robust across alternatives and dosage sets for power in detecting differences between genotypes.

Chemotherapeutic-induced apoptosis: a phenotype for pharmacogenomics studies

AIM

To determine whether cellular apoptosis is a suitable phenotypic trait for pharmacogenomics studies by evaluating caspase 3/7-mediated activity in lymphoblastoid cell lines after treatment with six chemotherapeutic agents: 5'-deoxyfluorouridine, pemetrexed, cytarabine, paclitaxel, carboplatin, and cisplatin.

MATERIALS AND METHODS

Using monozygotic twin pair and sibling pair lymphoblastoid cell lines, we identified conditions for measurement of caspase 3/7 activity in lymphoblastoid cell lines. Genome-wide association studies were performed with over 2 million single nucleotide polymorphisms (SNPs) and cisplatin-induced apoptosis in HapMap CEU cell lines (n=77).

RESULTS

Although treatment with 5'-deoxyfluorouridine and pemetrexed for up to 24 h resulted in low levels of apoptosis or interindividual variation in caspase-dependent cell death; paclitaxel, cisplatin, carboplatin, and cytarabine treatment for 24 h resulted in 9.4-fold, 9.1-fold, 7.0-fold, and 6.0-fold increases in apoptosis relative to control, respectively. There was a weak correlation between caspase activity and cytotoxicity (r(2)=0.03-0.29) demonstrating that cytotoxicity and apoptosis are two distinct phenotypes that may produce independent genetic associations. Estimated heritability (h(2)) for apoptosis was 0.57 and 0.29 for cytarabine (5 and 40 μmol/l, respectively), 0.22 for paclitaxel (12.5 nmol/l), and 0.34 for cisplatin (5 μmol/l). In the genome-wide association study using the HapMap CEU panel, we identified a significant enrichment of cisplatin-induced cytotoxicity SNPs within the significant cisplatin-induced apoptosis SNPs and an enrichment of expression quantitative trait loci (eQTL). Among these eQTLs, we identified several eQTLs with known function related to apoptosis and/or cytotoxicity.

CONCLUSION

Our study identifies apoptosis as a phenotype for pharmacogenomic studies in lymphoblastoid cell lines after treatment with paclitaxel, cisplatin, carboplatin, and cytarabine that may have utility for discovering biomarkers to predict response to certain chemotherapeutics.

Genomic profiling in CEPH cell lines distinguishes between the camptothecins and indenoisoquinolines

We have attempted to use a familial genetics strategy to study mechanisms of topoisomerase 1 (Top1) inhibition. Investigations have steadily been chipping away at the pathways involved in cellular response following Top1 inhibition for more than 20 years. Our system-wide approach, which phenotypes a collection of genotyped human cell lines for sensitivity to compounds and interrogates all genes and molecular pathways simultaneously. Previously, we characterized the in vitro sensitivity of 15 families of Centre d'Etude Polymorphisme Humain (CEPH) cell lines (n = 142) to 9 camptothecin analogues. Linkage analysis revealed a pattern of 7 quantitative trait loci (QTL) shared by all of the camptothecins. To identify which, if any, QTLs are related to the general mechanism of Top1 inhibition or should be considered camptothecin specific, we characterized the in vitro sensitivity of the same panel of CEPH cell lines to the indenisoquinolones, a structurally distinct class of Top1 inhibitors. Four QTLs on chromosomes 1, 5, 11, and 16 were shared by both the camptothecins and the indenoisoquinolines and are considered associated with the general mechanism of Top1 inhibition. The remaining 3 QTLs (chromosomes 6 and 20) are considered specific to camptothecin-induced cytotoxicity. Finally, 8 QTLs were identified, which were unique to the indenoisoquinolines.

Genome-wide local ancestry approach identifies genes and variants associated with chemotherapeutic susceptibility in African Americans

Chemotherapeutic agents are used in the treatment of many cancers, yet variable resistance and toxicities among individuals limit successful outcomes. Several studies have indicated outcome differences associated with ancestry among patients with various cancer types. Using both traditional SNP-based and newly developed gene-based genome-wide approaches, we investigated the genetics of chemotherapeutic susceptibility in lymphoblastoid cell lines derived from 83 African Americans, a population for which there is a disparity in the number of genome-wide studies performed. To account for population structure in this admixed population, we incorporated local ancestry information into our association model. We tested over 2 million SNPs and identified 325, 176, 240, and 190 SNPs that were suggestively associated with cytarabine-, 5'-deoxyfluorouridine (5'-DFUR)-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-4)). Importantly, some of these variants are found only in populations of African descent. We also show that cisplatin-susceptibility SNPs are enriched for carboplatin-susceptibility SNPs. Using a gene-based genome-wide association approach, we identified 26, 11, 20, and 41 suggestive candidate genes for association with cytarabine-, 5'-DFUR-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-3)). Fourteen of these genes showed evidence of association with their respective chemotherapeutic phenotypes in the Yoruba from Ibadan, Nigeria (p<0.05), including TP53I11, COPS5 and GAS8, which are known to be involved in tumorigenesis. Although our results require further study, we have identified variants and genes associated with chemotherapeutic susceptibility in African Americans by using an approach that incorporates local ancestry information.

Interactions between glucocorticoid treatment and cis-regulatory polymorphisms contribute to cellular response phenotypes

Glucocorticoids (GCs) mediate physiological responses to environmental stress and are commonly used as pharmaceuticals. GCs act primarily through the GC receptor (GR, a transcription factor). Despite their clear biomedical importance, little is known about the genetic architecture of variation in GC response. Here we provide an initial assessment of variability in the cellular response to GC treatment by profiling gene expression and protein secretion in 114 EBV-transformed B lymphocytes of African and European ancestry. We found that genetic variation affects the response of nearby genes and exhibits distinctive patterns of genotype-treatment interactions, with genotypic effects evident in either only GC-treated or only control-treated conditions. Using a novel statistical framework, we identified interactions that influence the expression of 26 genes known to play central roles in GC-related pathways (e.g. NQO1, AIRE, and SGK1) and that influence the secretion of IL6.

Discovery of regulatory elements in human ATP-binding cassette transporters through expression quantitative trait mapping

ATP-Binding Cassette (ABC) membrane transporters determine the disposition of many drugs, metabolites and endogenous compounds. Coding region variation in ABC transporters is the cause of many genetic disorders, but much less is known about the genetic basis and functional outcome of ABC transporter expression level variation. We used genotype and mRNA transcript level data from human lymphoblastoid cell lines to assess population and gender differences in ABC transporter expression, and to guide the discovery of genomic regions involved in transcriptional regulation. Nineteen of 49 ABC genes were differentially expressed between individuals of African, Asian and European descent suggesting an important influence of race on expression level of ABC transporters. Twenty-four significant associations were found between transporter transcript levels and proximally located genetic variants. Several of the associations were experimentally validated in reporter assays. Through influencing ABC expression levels, these SNPs may affect disease susceptibility and response to drugs.

Comprehensive evaluation of the contribution of X chromosome genes to platinum sensitivity

Using a genome-wide gene expression data set generated from Affymetrix GeneChip Human Exon 1.0ST array, we comprehensively surveyed the role of 322 X chromosome gene expression traits on cellular sensitivity to cisplatin and carboplatin. We identified 31 and 17 X chromosome genes whose expression levels are significantly correlated (after multiple testing correction) with sensitivity to carboplatin and cisplatin, respectively, in the combined HapMap CEU (Utah residents with ancestry from northern and western Europe) and YRI (Yoruba in Ibahan, Nigeria) populations (false discovery rate, FDR < 0.05). Of those, 14 overlap for both cisplatin and carboplatin. Using an independent gene expression quantification method, the Illumina Sentrix Human-6 Expression BeadChip, measured on the same HapMap cell lines, we found that 4 and 2 of these genes are significantly associated with carboplatin and cisplatin sensitivity, respectively, in both analyses. Two genes, CTPS2 and DLG3, were identified by both genome-wide gene expression analyses as correlated with cellular sensitivity to both platinating agents. The expression of DLG3 gene was also found to correlate with cellular sensitivity to platinating agents in NCI-60 cancer cell lines. In addition, we evaluated whether the expression of X chromosome genes contributed to the observed differences in sensitivity to the platinums between CEU and YRI-derived cell lines. Of the 34 distinct genes significantly correlated with either carboplatin or cisplatin sensitivity, 14 are differentially expressed (defined as P < 0.05) between CEU and YRI. Thus, sex chromosome genes play a role in cellular sensitivity to platinating agents and differences in the expression level of these genes are an important source of variation that should be included in comprehensive pharmacogenomic studies.

Population differences in the rate of proliferation of international HapMap cell lines

The International HapMap Project is a resource for researchers containing genotype, sequencing, and expression information for EBV-transformed lymphoblastoid cell lines derived from populations across the world. The expansion of the HapMap beyond the four initial populations of Phase 2, referred to as Phase 3, has increased the sample number and ethnic diversity available for investigation. However, differences in the rate of cellular proliferation between the populations can serve as confounders in phenotype-genotype studies using these cell lines. Within the Phase 2 populations, the JPT and CHB cell lines grow faster (p < 0.0001) than the CEU or YRI cell lines. Phase 3 YRI cell lines grow significantly slower than Phase 2 YRI lines (p < 0.0001), with no widespread genetic differences based on common SNPs. In addition, we found significant growth differences between the cell lines in the Phase 2 ASN populations and the Han Chinese from the Denver metropolitan area panel in Phase 3 (p < 0.0001). Therefore, studies that separate HapMap panels into discovery and replication sets must take this into consideration.

Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS

Although genome-wide association studies (GWAS) of complex traits have yielded more reproducible associations than had been discovered using any other approach, the loci characterized to date do not account for much of the heritability to such traits and, in general, have not led to improved understanding of the biology underlying complex phenotypes. Using a web site we developed to serve results of expression quantitative trait locus (eQTL) studies in lymphoblastoid cell lines from HapMap samples (http://www.scandb.org), we show that single nucleotide polymorphisms (SNPs) associated with complex traits (from http://www.genome.gov/gwastudies/) are significantly more likely to be eQTLs than minor-allele-frequency–matched SNPs chosen from high-throughput GWAS platforms. These findings are robust across a range of thresholds for establishing eQTLs (p-values from 10⁻⁴–10⁻⁸), and a broad spectrum of human complex traits. Analyses of GWAS data from the Wellcome Trust studies confirm that annotating SNPs with a score reflecting the strength of the evidence that the SNP is an eQTL can improve the ability to discover true associations and clarify the nature of the mechanism driving the associations. Our results showing that trait-associated SNPs are more likely to be eQTLs and that application of this information can enhance discovery of trait-associated SNPs for complex phenotypes raise the possibility that we can utilize this information both to increase the heritability explained by identifiable genetic factors and to gain a better understanding of the biology underlying complex traits.

We show here that single nucleotide polymorphisms (SNPs) associated with complex traits (as identified in the catalog of results from genome-wide association studies http://www.genome.gov/gwastudies/) are more likely than other SNPs chosen from high-throughput genotyping platforms to predict expression levels of genes. These observations confirm that genetic risk factors for complex traits will often affect phenotype by altering the amount or timing of protein production, rather than by changing the type of protein produced. This knowledge can be used to improve our ability to discover genetic risk factors for complex traits and to improve our understanding of their underlying biology.