An eQTL-based method identifies CTTN and ZMAT3 as pemetrexed susceptibility markers.

Pemetrexed, approved for the treatment of non-small cell lung cancer and malignant mesothelioma, has adverse effects including neutropenia, leucopenia, thrombocytopenia, anemia, fatigue and nausea. The results we report here represent the first genome-wide study aimed at identifying genetic predictors of pemetrexed response. We utilized expression quantitative trait loci (eQTLs) mapping combined with drug-induced cytotoxicity data to gain mechanistic insights into the observed genetic associations with pemetrexed susceptibility. We found that CTTN and ZMAT3 expression signature explained >30% of the pemetrexed susceptibility phenotype variation for pemetrexed in the discovery population. Replication using PCR and a semi-high-throughput, scalable assay system confirmed the initial discovery results in an independent set of samples derived from the same ancestry. Furthermore, functional validation in both germline and tumor cells demonstrates a decrease in cell survival following knockdown of CTTN or ZMAT3. In addition to our particular findings on genetic and gene expression predictors of susceptibility phenotype for pemetrexed, the work presented here will be valuable to the robust discovery and validation of genetic determinants and gene expression signatures of various chemotherapeutic susceptibilities.

Population differences in microRNA expression and biological implications.

Population differences observed for complex traits may be attributed to the combined effect of socioeconomic, environmental, genetic and epigenetic factors. To better understand population differences in complex traits, genome-wide genetic and gene expression differences among ethnic populations have been studied. Here we set out to evaluate population differences in small non-coding RNAs through an evaluation of microRNA (miRNA) baseline expression in HapMap lymphoblastoid cell lines (LCLs) derived from 53 CEU (Utah residents with northern and western European ancestry) and 54 YRI (African from Ibadan, Nigeria). Using the Exiqon miRCURYTM LNA arrays, we found that 16% of all miRNAs evaluated in our study differ significantly between these 2 ethnic groups (pBonferroni corrected< 0.05). Furthermore, we explored the potential biological function of these observed differentially expressed miRNAs by comprehensively examining their effect on the transcriptome and their relationship with cellular sensitivity drug phenotypes. After multiple testing adjustment (false discovery rate (FDR)< 0.1), we found that 55% and 88% of the differentially expressed miRNAs were significantly and inversely correlated with an mRNA expression phenotype in the CEU and YRI samples, respectively. Interestingly, a substantial proportion (64%) of these miRNAs correlated with cellular sensitivity to chemotherapeutic agents (FDR< 0.05). Lastly, upon performing a genome-wide association study between SNPs and miRNA expression, we identified a large number of SNPs exhibiting different allele frequencies that affect the expression of these differentially expressed miRNAs, suggesting the role of genetic variants in mediating the observed population differences.

Expression and alternative splicing of folate pathway genes in HapMap lymphoblastoid cell lines.

AIM: Folate is vital for cell growth and development through its important role in one-carbon metabolism - an essential process in the synthesis of amino acids and nucleic acids. Folate pathway genes have been considered as therapeutic targets of drugs for the treatment of cancer and other diseases. Racial and ethnic disparities of folate metabolism and outcome of antifolate therapies have been reported. In this study, we evaluate the genetic regulation for expression and alternative splicing of folate related genes in HapMap lymphoblastoid cell lines (LCLs) of individuals of European and African descent. MATERIALS & METHODS: Gene and exon level expression and alternative splicing of folate pathway genes were compared in LCLs derived from the Centre d'Etude du Polymorphisme Humain (CEPH) from Utah (CEU) and the Yoruba from Ibadan (YRI) using a permutation-based test. A genome-wide association study was performed to search for SNPs associated with folate pathway gene expressions and alternative splicing in the combined population samples. RESULTS: A total of 52 folate pathway genes were evaluated in the analysis of which 46 were expressed in the LCLs. There were 12 genes (26%) with differential gene-level expression and 23 genes (50%) with differential alternative splicing for exons or UTRs between the CEU and the YRI (permutation p

The Werner's syndrome 4330T>C (Cys1367Arg) gene variant does not affect the in vitro cytotoxicity of topoisomerase inhibitors and platinum compounds.

PURPOSE: Werner's syndrome (WS) is a recessive disorder of premature onset of processes associated with aging. Defective DNA repair has been reported after exposure of cells isolated from WS patients to DNA-damaging agents. The germline 4330T>C (Cys1367Arg) variant in the WS gene (WRN) has been associated with protection from age-related diseases, suggesting it has a functional role. We studied whether the 4330T>C variant confers altered drug sensitivity in vitro. METHODS: 4330T>C was genotyped in 372 human lymphoblastoid cell lines (LCLs) from unrelated healthy Caucasian individuals using a TaqMan-based method. The study was powered to detect the effect of the 4330T>C genotypes after exposure to camptothecin (based upon preliminary data). The effect of the 4330T>C variant on the cytotoxicity of etoposide, carboplatin, cisplatin and daunorubicin was also tested. WRN expression in 57 LCLs was measured by microarray. RESULTS: No significant difference between the IC50 of the cells was observed among genotypes (P = 0.46) after exposure to camptothecin. No association was also observed for etoposide, carboplatin, cisplatin, and daunorubicin (ANOVA, P > 0.05). WRN expression also did not vary across genotypes (ANOVA, P = 0.37). CONCLUSION: These results suggest that this nonsynonymous variant has relatively normal function at the cellular level.

Identification of common genetic variants that account for transcript isoform variation between human populations.

In addition to the differences between populations in transcriptional and translational regulation of genes, alternative pre-mRNA splicing (AS) is also likely to play an important role in regulating gene expression and generating variation in mRNA and protein isoforms. Recently, the genetic contribution to transcript isoform variation has been reported in individuals of recent European descent. We report here results of an investigation of the differences in AS patterns between human populations. AS patterns in 176 HapMap lymphoblastoid cell lines derived from individuals of European and African ancestry were evaluated using the Affymetrix GeneChip Human Exon 1.0 ST Array. A variety of biological processes such as response to stimulus and transcription were found to be enriched among the differentially spliced genes. The differentially spliced genes also include some involved in human diseases that have different prevalence or susceptibility between populations. The genetic contribution to the population differences in transcript isoform variation was then evaluated by a genome-wide association using the HapMap genotypic data on single nucleotide polymorphisms (SNPs). The results suggest that local and distant genetic variants account for a substantial fraction of the observed transcript isoform variation between human populations. Our findings provide new insights into the complexity of the human genome as well as the health disparities between the two populations.

Identification of genomic regions contributing to etoposide-induced cytotoxicity.

Etoposide is routinely used in combination-based chemotherapy for testicular cancer and small-cell lung cancer; however, myelosuppression, therapy-related leukemia and neurotoxicity limit its utility. To determine the genetic contribution to cellular sensitivity to etoposide, we evaluated cell growth inhibition in Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines from 24 multi-generational pedigrees (321 samples) following treatment with 0.02-2.5 microM etoposide for 72 h. Heritability analysis showed that genetic variation contributes significantly to the cytotoxic phenotypes (h (2) = 0.17-0.25, P = 4.9 x 10(-5)-7.3 x 10(-3)). Whole genome linkage scans uncovered 8 regions with peak LOD scores ranging from 1.57 to 2.55, with the most significant signals being found on chromosome 5 (LOD = 2.55) and chromosome 6 (LOD = 2.52). Linkage-directed association was performed on a subset of HapMap samples within the pedigrees to find 22 SNPs significantly associated with etoposide cytotoxicity at one or more treatment concentrations. UVRAG, a DNA repair gene, SEMA5A, SLC7A6 and PRMT7 are implicated from these unbiased studies. Our findings suggest that susceptibility to etoposide-induced cytotoxicity is heritable and using an integrated genomics approach we identified both genomic regions and SNPs associated with the cytotoxic phenotypes.

Genetic variants contributing to daunorubicin-induced cytotoxicity.

Identifying heritable genetic variants responsible for chemotherapeutic toxicities has been challenging due in part to its multigenic nature. To date, there is a paucity of data on genetic variants associated with patients experiencing severe myelosuppression or cardiac toxicity following treatment with daunorubicin. We present a genome-wide model using International HapMap cell lines that integrate genotype and gene expression to identify genetic variants that contribute to daunorubicin-induced cytotoxicity. A cell growth inhibition assay was used to measure variations in the cytotoxicity of daunorubicin. Gene expression was determined using the Affymetrix GeneChip Human Exon 1.0ST Array. Using sequential analysis, we evaluated the associations between genotype and cytotoxicity, those significant genotypes with gene expression and correlated gene expression of the identified candidates with cytotoxicity. A total of 26, 9, and 18 genetic variants were identified to contribute to daunorubicin-induced cytotoxicity through their effect on 16, 9, and 36 gene expressions in the combined, Centre d' Etude du Polymorphisme Humain (CEPH), and Yoruban populations, respectively. Using 50 non-HapMap CEPH cell lines, single nucleotide polymorphisms generated through our model predicted 29% of the overall variation in daunorubicin sensitivity and the expression of CYP1B1 was significantly correlated with sensitivity to daunorubicin. In the CEPH validation set, rs120525235 and rs3750518 were significant predictors of transformed daunorubicin IC(50) (P = 0.005 and P = 0.0008, respectively), and rs1551315 trends toward significance (P = 0.089). This unbiased method can be used to elucidate genetic variants contributing to a wide range of cellular phenotypes.

Identification of genetic variants and gene expression relationships associated with pharmacogenes in humans.

OBJECTIVES: The very important pharmacogenes (VIPs) were selected by Pharmacogenetic Research Network (National Institutes of Health-PGRN) owing to their significant effects on drug treatment both at the pharmacokinetic and pharmacodynamic levels. Our objective was to identify single nucleotide polymorphisms (SNPs) that potentially affected the expression of these genes or potential SNP-gene interactions involved to improve our understanding of genetic effects on drug therapy. BASIC METHODS: Gene expression was evaluated in 176 International HapMap lymphoblastoid cell lines derived from CEU (CEPH, Utah residents with ancestry from northern and western Europe; n=87) and YRI (Yoruba in Ibadan, Nigeria; n=89) using Affymetrix GeneChip Human Exon 1.0 ST arrays (Affymetrix Laboratory, Affymetrix Inc., Santa Clara, California, USA) with interrogation of greater than 17,000 human genes. Genome-wide association was performed between over two million publicly available HapMap SNPs and gene expression. MAIN RESULTS: The expression of two PGRN-VIPs (GSTT1 and GSTM1) are significantly associated with SNPs within 2.5 Mb of the genes; whereas the expression of three and ten PGRN-VIPs are significantly associated with distant-acting SNPs in CEU and YRI, respectively. In addition, three and four PGRN-VIPs harbor SNPs that are distantly associated with other gene expressions in CEU and YRI, respectively. PRINCIPAL CONCLUSION: Using this information, one may identify genetic variants that are significantly associated with the expression of any set of genes of interest; or evaluate potential gene-gene interaction through SNP expression relationships.

Genetic architecture of transcript-level variation in humans.

We report here the results of testing the pairwise association of 12,747 transcriptional gene-expression values with more than two million single-nucleotide polymorphisms (SNPs) in samples of European (CEPH from Utah; CEU) and African (Yoruba from Ibadan; YRI) ancestry. We found 4,677 and 5,125 significant associations between expression quantitative nucleotides (eQTNs) and transcript clusters in the CEU and the YRI samples, respectively. The physical distance between an eQTN and its associated transcript cluster was referred to as the intrapair distance. An association with 4 Mb or less intrapair distance was defined as local; otherwise, it was defined as distant. The enrichment analysis of functional categories shows that genes harboring the local eQTNs are enriched in the categories related to nucleosome and chromatin assembly; the genes harboring the distant eQTNs are enriched in the categories related to transmembrane signal transduction, suggesting that these biological pathways are likely to play a significant role in regulation of gene expression. We highlight in the EPHX1 gene a deleterious nonsynonymous SNP that is distantly associated with gene expression of ORMDL3, a susceptibility gene for asthma.

Evaluation of genetic variation contributing to differences in gene expression between populations.

Gene expression is a complex quantitative trait partially regulated by genetic variation in DNA sequence. Population differences in gene expression could contribute to some of the observed differences in susceptibility to common diseases and response to drug treatments. We characterized gene expression in the full set of HapMap lymphoblastoid cell lines derived from individuals of European and African ancestry for 9156 transcript clusters (gene-level) evaluated with the Affymetrix GeneChip Human Exon 1.0 ST Array. Gene expression was found to differ significantly between these samples for 383 transcript clusters. Biological processes including ribosome biogenesis and antimicrobial humoral response were found to be enriched in these differential genes, suggesting their possible roles in contributing to the population differences at a higher level than that of mRNA expression and in response to environmental information. Genome-wide association studies for local or distant genetic variants that correlate with the differentially expressed genes enabled identification of significant associations with one or more single-nucleotide polymorphisms (SNPs), consistent with the hypothesis that genetic factors and not simply population identity or other characteristics (age of cell lines, length of culture, etc.) contribute to differences in gene expression in these samples. Our results provide a comprehensive view of the genes differentially expressed between populations and the enriched biological processes involved in these genes. We also provide an evaluation of the contributions of genetic variation and nongenetic factors to the population differences in gene expression.

A genome-wide approach to identify genetic variants that contribute to etoposide-induced cytotoxicity.

Large interindividual variance has been observed in sensitivity to drugs. To comprehensively decipher the genetic contribution to these variations in drug susceptibility, we present a genome-wide model using human lymphoblastoid cell lines from the International HapMap consortium, of which extensive genotypic information is available, to identify genetic variants that contribute to chemotherapeutic agent-induced cytotoxicity. Our model integrated genotype, gene expression, and sensitivity of HapMap cell lines to drugs. Cell lines derived from 30 trios of European descent (Center d'Etude du Polymorphisme Humain population) and 30 trios of African descent (Yoruban population) were used. Cell growth inhibition at increasing concentrations of etoposide for 72 h was determined by using alamarBlue assay. Gene expression on 176 HapMap cell lines (87 Center d'Etude du Polymorphisme Humain population and 89 Yoruban population) was determined by using the Affymetrix GeneChip Human Exon 1.0ST Array. We evaluated associations between genotype and cytotoxicity, genotype and gene expression and correlated gene expression of the identified candidates with cytotoxicity. The analysis identified 63 genetic variants that contribute to etoposide-induced toxicity through their effect on gene expression. These include genes that may play a role in cancer (AGPAT2, IL1B, and WNT5B) and genes not yet known to be associated with sensitivity to etoposide. This unbiased method can be used to elucidate genetic variants contributing to a wide range of cellular phenotypes induced by chemotherapeutic agents.

Mapping genes that contribute to daunorubicin-induced cytotoxicity.

Daunorubicin is an anthracycline antibiotic agent used in the treatment of hematopoietic malignancies. Toxicities associated with this agent include myelosuppression and cardiotoxicity; however, the genes or genetic determinants that contribute to these toxicities are unknown. We present an unbiased genome-wide approach that incorporates heritability, whole-genome linkage analysis, and linkage-directed association to uncover genetic variants contributing to the sensitivity to daunorubicin-induced cytotoxicity. Cell growth inhibition in 324 Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines (24 pedigrees) was evaluated following treatment with daunorubicin for 72 h. Heritability analysis showed a significant genetic component contributing to the cytotoxic phenotypes (h2 = 0.18-0.63 at 0.0125, 0.025, 0.05, 0.1, 0.2, and 1.0 mumol/L daunorubicin and at the IC50, the dose required to inhibit 50% cell growth). Whole-genome linkage scans at all drug concentrations and IC50 uncovered 11 regions with moderate peak LOD scores (> 1.5), including 4q28.2 to 4q32.3 with a maximum LOD score of 3.18. The quantitative transmission disequilibrium tests were done using 31,312 high-frequency single-nucleotide polymorphisms (SNP) located in the 1 LOD confidence interval of these 11 regions. Thirty genes were identified as significantly associated with daunorubicin-induced cytotoxicity (P < or = 2.0 x 10(-4), false discovery rate < or = 0.1). Pathway and functional gene ontology analysis showed that these genes were overrepresented in the phosphatidylinositol signaling system, axon guidance pathway, and GPI-anchored proteins family. Our findings suggest that a proportion of susceptibility to daunorubicin-induced cytotoxicity may be controlled by genetic determinants and that analysis using linkage-directed association studies with dense SNP markers can be used to identify the genetic variants contributing to cytotoxicity.

Gender-specific differences in expression in human lymphoblastoid cell lines.

Women and men have different risks for certain diseases and they often respond differently to treatment. These differences could be due to the sex-specific differences in the expression of genes related to primary disease susceptibility or pharmacodynamic targets. To evaluate the sex-specific pattern of gene expression, we compared gene expression levels using a publicly available microarray dataset of 233 (115 women and 118 men) lymphoblastoid cell lines. From the 4799 probes meeting a specified minimal level of expression, 10 genes (P<0.005, permutation adjusted false discovery rate less than 50%) located on autosomal chromosomes were identified using a permutation-based approach. These genes were found to be over-represented in certain gene ontology terms of biological process (cell adhesion, apoptosis, transcription and signal transduction), and molecular function (structural molecule activity, zinc ion binding, transcription factor activity and protein binding). A Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that two known pathways are over-represented: adherens junction and cytokine-cytokine receptor interaction.

Effect of population and gender on chemotherapeutic agent-induced cytotoxicity.

Large interindividual variance is observed in both response and toxicity associated with chemotherapy. Our goal is to identify factors that contribute to chemotherapy-induced toxicity. To this end, we used EBV-transformed B-lymphoblastoid HapMap cell lines derived from 30 Yoruban trios (African descent) and 30 Centre d' Etude du Polymorphisme Humain (CEPH) trios (European descent) to evaluate population- and gender-specific differences in cytotoxicity of carboplatin, cisplatin, daunorubicin, and etoposide using a high-throughput, short-term cytotoxicity assay. The IC(50) was compared for population- and gender-specific differences for the four drugs. We observed large interindividual variance in IC(50) values for carboplatin, cisplatin, daunorubicin, and etoposide for both Yoruban and CEPH populations (range from 8- to 433-fold). Statistically significant differences in carboplatin and daunorubicin IC(50) were shown when comparing Yoruban cell lines (n = 89) to CEPH cell lines (n = 87; P = 0.002 and P = 0.029, respectively). This population difference in treatment induced cytotoxicity was not seen for either cisplatin or etoposide. In the Yoruban population, cell lines derived from females were less sensitive to platinating agents than males [median carboplatin IC(50), 29.1 versus 24.6 micromol/L (P = 0.012); median cisplatin IC(50), 7.0 versus 6.0 micromol/L (P = 0.020) in female and male, respectively]. This difference was not observed in the CEPH population. These results show that population and gender may affect risk for toxicities associated with certain chemotherapeutic agents.

Role of GADD34 in modulation of cisplatin cytotoxicity.

Cisplatin and carboplatin are widely used clinical chemotherapeutic agents, especially against testicular, ovarian, and head and neck cancers. O(6)-Benzylguanine (BG) has been shown to result in enhanced cytotoxicity, apoptosis, and DNA platination when used in conjunction with cisplatin and carboplatin in head and neck cancer cell lines. Microarray expression data indicated overexpression of 19 genes and underexpression of 22 genes specific to treatment with the combination of BG+/-cisplatin compared to cisplatin alone treatment in SQ20b head and neck cancer cells (p<0.05) using the Affymetrix HG-U133A GeneChip((R)). Among the overexpressed probe sets were genes involved in DNA damage and apoptosis, including GADD34, DDIT4, ATF4, and PTHLH. A similarly structured analog, 9-CH(3)-BG, does not enhance cisplatin-induced cytotoxicity or apoptosis nor is there enhanced expression of GADD34 in cisplatin or 9-CH(3)-BG+/-cisplatin-treated cells compared to control cells. Analysis of cells exposed to 9-CH(3)-BG+/-cisplatin allowed us to focus our array list on 32 probe sets specific to BG+cisplatin versus cisplatin, ruling out differentially expressed probe sets common to 9-CH(3)-BG+cisplatin versus cisplatin. Similarly, 14 probe sets were specific to BG+/-cisplatin versus BG, ruling out differentially expressed probe sets common to 9-CH(3)-BG+/-cisplatin versus 9-CH(3)-BG. Quantitative real-time PCR demonstrated a dose dependent increase in GADD34 expression in cells exposed to BG+/-cisplatin with levels approximately >2-fold for cells exposed to BG+cisplatin compared to cisplatin alone. Levels of GADD34 transcripts were determined with both cisplatin and BG+cisplatin at several different time points concomitant with and following drug treatment. At all timepoints, GADD34 transcript levels are approximately two-fold elevated in cells treated with BG+cisplatin compared to cisplatin alone. Furthermore, significant changes in GADD34 expression levels in SQ20b, SCC35, and SCC61 cells, with approximately three-fold, two-fold, and 3.5-fold increases in expression, respectively, upon treatment with BG+/-cisplatin compared with control. Elucidation of these molecular pathways will aid in our goal of synthesizing more powerful modulators to increase efficacy of platinum agents.