Predictive Biomarkers of Overall Survival in Patients with Metastatic Renal Cell Carcinoma Treated with IFNα ± Bevacizumab: Results from CALGB 90206 (Alliance).

PURPOSE: CALGB 90206 was a phase III trial of 732 patients with metastatic renal cell carcinoma (mRCC) comparing bevacizumab plus IFNα (BEV + IFN) with IFNα alone (IFN). No difference in overall survival (OS) was observed. Baseline samples were analyzed to identify predictive biomarkers for survival benefit. PATIENTS AND METHODS: A total of 32 biomarkers were assessed in 498 consenting patients randomly assigned into training (n = 279) and testing (n = 219) sets. The proportional hazards model was used to test for treatment arm and biomarker interactions of OS. The estimated coefficients from the training set were used to compute a risk score for each patient and to classify patients by risk in the testing set. The resulting model was assessed for predictive accuracy using the time-dependent area under the ROC curve (tAUROC). RESULTS: A statistically significant three-way interaction between IL6, hepatocyte growth factor (HGF), and bevacizumab treatment was observed in the training set and confirmed in the testing set (P < 0.0001). The model based on IL6, HGF, and bevacizumab treatment was predictive of OS (P < 0.001), with the high- and low-risk groups having a median OS of 10.2 [95% confidence interval (CI), 8.0-13.8] and 34.3 (95% CI, 28.5-40.5) months, respectively. The average tAUROC for the final model of OS based on 100 randomly split testing sets was 0.78 (first, third quartiles = 0.77, 0.79). CONCLUSIONS: IL6 and HGF are potential predictive biomarkers of OS benefit from BEV + IFN in patients with mRCC. The model based on key biological and clinical factors demonstrated predictive efficacy for OS. These markers warrant further validation in future anti-VEGF and immunotherapy in mRCC trials. See related commentaries by Mishkin and Kohn, p. 2722 and George and Bertagnolli, p. 2725.

Long-Term Recovery of the Adaptive Immune System in Rhesus Macaques After Total Body Irradiation.

PURPOSE: Ionizing radiation causes acute damage to hematopoietic and immune cells, but the long-term immunologic consequences of irradiation are poorly understood. We therefore performed a prospective study of the delayed immune effects of radiation using a rhesus macaque model. METHODS AND MATERIALS: Ten macaques received 4 Gy high-energy x-ray total body irradiation (TBI) and 6 control animals received sham irradiation. TBI caused transient lymphopenia that resolved over several weeks. Once white blood cell counts recovered, flow cytometry was used to immunophenotype the circulating adaptive immune cell populations 4, 9, and 21 months after TBI. Data were fit using a mixed-effects model to determine age-dependent, radiation-dependent, and interacting effects. T cell receptor (TCR) sequencing and quantification of TCR Excision Circles were used to determine relative contributions of thymopoiesis and peripheral expansion to T cell repopulation. Two years after TBI, the cohort was vaccinated with a 23-valent pneumococcal polysaccharide vaccine and a tetravalent influenza hemagglutinin vaccine. RESULTS: Aging, but not TBI, led to significant changes in the frequencies of dendritic cells, CD4 and CD8 T cells, and B cells. However, irradiated animals exhibited increased frequencies of central memory T cells and decreased frequencies of naïve T cells. These consequences of irradiation were time-dependent and more prolonged in the CD8 T cell population. Irradiation led to transient increases in CD8+ T cell TCR Excision Circles and had no significant effect on TCR sequence entropy, indicating T cell recovery was partially mediated by thymopoiesis. Animals that were irradiated and then vaccinated showed normal immunoglobulin G binding and influenza neutralization titers in response to the 4 protein antigens but weaker immunoglobulin G binding titers to 10 of the 23 polysaccharide antigens. CONCLUSIONS: These findings indicate that TBI causes subtle but long-lasting immune defects that are evident years after recovery from lymphopenia.

fastJT: An R package for robust and efficient feature selection for machine learning and genome-wide association studies.

BACKGROUND: Parametric feature selection methods for machine learning and association studies based on genetic data are not robust with respect to outliers or influential observations. While rank-based, distribution-free statistics offer a robust alternative to parametric methods, their practical utility can be limited, as they demand significant computational resources when analyzing high-dimensional data. For genetic studies that seek to identify variants, the hypothesis is constrained, since it is typically assumed that the effect of the genotype on the phenotype is monotone (e.g., an additive genetic effect). Similarly, predictors for machine learning applications may have natural ordering constraints. Cross-validation for feature selection in these high-dimensional contexts necessitates highly efficient computational algorithms for the robust evaluation of many features. RESULTS: We have developed an R extension package, fastJT, for conducting genome-wide association studies and feature selection for machine learning using the Jonckheere-Terpstra statistic for constrained hypotheses. The kernel of the package features an efficient algorithm for calculating the statistics, replacing the pairwise comparison and counting processes with a data sorting and searching procedure, reducing computational complexity from O(n2) to O(n log(n)). The computational efficiency is demonstrated through extensive benchmarking, and example applications to real data are presented. CONCLUSIONS: fastJT is an open-source R extension package, applying the Jonckheere-Terpstra statistic for robust feature selection for machine learning and association studies. The package implements an efficient algorithm which leverages internal information among the samples to avoid unnecessary computations, and incorporates shared-memory parallel programming to further boost performance on multi-core machines.

Identification of Novel Mast Cell Activators Using Cell-Based High-Throughput Screening.

Mast cells (MCs) are known to regulate innate and adaptive immunity. MC activators have recently been described as safe and effective vaccine adjuvants. Many currently known MC activators are inadequate for in vivo applications, however, and research on identifying novel MC activators is limited. In this study, we identified novel MC activators by using high-throughput screening (HTS) assays using approximately 55,000 small molecules. Data sets obtained by the primary HTS assays were statistically evaluated using quality control rules and the B-score calculation, and compounds with B-scores of >3.0 were chosen as mast cell activators (hits). These hits were re-evaluated with secondary and tertiary HTS assays, followed by further statistical analysis. From these hits, we selected 15 compounds that caused degranulation in murine and human MCs, with potential for flexible chemical modification for further study. Among these 15 compounds, ST101036, ST029248, and ST026567 exhibited higher degranulation potency than other hit compounds in both human and mouse MCs. In addition, the 15 compounds identified promote de novo synthesis of cytokines and induce the release of eicosanoids from human and mouse MCs. HTS enabled us to identify small-molecule MC activators with unique properties that may be useful as vaccine adjuvants.

Genetic variation determines VEGF-A plasma levels in cancer patients.

Angiogenesis is essential in tumor biology and is regulated by vascular endothelial growth factor (VEGF) ligands and receptors. Here we aimed to discover genetic variants associated with levels of circulating angiogenic proteins in cancer patients. Plasma was collected at baseline in 216 pancreatic and 114 colorectal cancer patients. Thirty-one angiogenic proteins were measured by ELISA. 484,523 Single Nucleotide Polymorphisms (SNP) were tested for association with plasma levels for each protein in pancreatic cancer patients. Three top-ranked hits were then genotyped in colorectal cancer patients, where associations with the same proteins were measured. The results demonstrated rs2284284 and MCP1 (P-value = 6.7e-08), rs7504372 and VEGF-C (P-value = 9.8e-09), and rs7767396 and VEGF-A (P-value = 5.8e-09) were SNP-protein pairs identified in pancreatic cancer patients. In colorectal cancer patients, only rs7767396 (A > G) and VEGF-A was validated (P-value = 5.18e-05). The AA genotype of rs7767396 exhibited 2.04-2.3 and 2.7-3.4-fold higher VEGF-A levels than those with AG and GG genotypes. The G allele of rs7767396 reduces binding of the NF-AT1 transcription factor. In conclusion, a common genetic variant predicts the plasma levels of VEGF-A in cancer patients through altered binding of NF-AT1.

Impact of early life exposure to ionizing radiation on influenza vaccine response in an elderly Japanese cohort.

The objective of this study was to evaluate effects of whole body radiation exposure early in life on influenza vaccination immune responses much later in life. A total of 292 volunteers recruited from the cohort members of ongoing Adult Health Study (AHS) of Japanese atomic bomb (A-bomb) survivors completed this observational study spanning two influenza seasons (2011-2012 and 2012-2013). Peripheral blood samples were collected prior to and three weeks after vaccination. Serum hemagglutination inhibition (HAI) antibody titers were measured as well as concentrations of 25 cytokines and chemokines in culture supernatant from peripheral blood mononuclear cells, with and without in vitro stimulation with influenza vaccine. We found that influenza vaccination modestly enhanced serum HAI titers in this unique cohort of elderly subjects, with seroprotection ranging from 18 to 48% for specific antigen/season combinations. Twelve percent of subjects were seroprotected against all three vaccine antigens post-vaccination. Males were generally more likely to be seroprotected for one or more antigens post-vaccination, with no differences in vaccine responses based on age at vaccination or radiation exposure in early life. These results show that early life exposure to ionizing radiation does not prevent responses of elderly A-bomb survivors to seasonal influenza vaccine.

Bayesian Multi-Plate High-Throughput Screening of Compounds.

High-throughput screening of compounds (chemicals) is an essential part of drug discovery, involving thousands to millions of compounds, with the purpose of identifying candidate hits. Most statistical tools, including the industry standard B-score method, work on individual compound plates and do not exploit cross-plate correlation or statistical strength among plates. We present a new statistical framework for high-throughput screening of compounds based on Bayesian nonparametric modeling. The proposed approach is able to identify candidate hits from multiple plates simultaneously, sharing statistical strength among plates and providing more robust estimates of compound activity. It can flexibly accommodate arbitrary distributions of compound activities and is applicable to any plate geometry. The algorithm provides a principled statistical approach for hit identification and false discovery rate control. Experiments demonstrate significant improvements in hit identification sensitivity and specificity over the B-score and R-score methods, which are highly sensitive to threshold choice. These improvements are maintained at low hit rates. The framework is implemented as an efficient R extension package BHTSpack and is suitable for large scale data sets.

Radiation- and Age-Associated Changes in Peripheral Blood Dendritic Cell Populations among Aging Atomic Bomb Survivors in Japan.

Previous immunological studies in atomic bomb survivors have suggested that radiation exposure leads to long-lasting changes, similar to immunological aging observed in T-cell-adaptive immunity. However, to our knowledge, late effects of radiation on dendritic cells (DCs), the key coordinators for activation and differentiation of T cells, have not yet been investigated in humans. In the current study, we hypothesized that numerical and functional decreases would be observed in relationship to radiation dose in circulating conventional DCs (cDCs) and plasmacytoid DCs (pDCs) among 229 Japanese A-bomb survivors. Overall, the evidence did not support this hypothesis, with no overall changes in DCs or functional changes observed with radiation dose. Multivariable regression analysis for radiation dose, age and gender effects revealed that total DC counts as well as subpopulation counts decreased in relationship to increasing age. Further analyses revealed that in women, absolute numbers of pDCs showed significant decreases with radiation dose. A hierarchical clustering analysis of gene expression profiles in DCs after Toll-like receptor stimulation in vitro identified two clusters of participants that differed in age-associated expression levels of genes involved in antigen presentation and cytokine/chemokine production in cDCs. These results suggest that DC counts decrease and expression levels of gene clusters change with age. More than 60 years after radiation exposure, we also observed changes in pDC counts associated with radiation, but only among women.

A robust microparticle platform for a STING-targeted adjuvant that enhances both humoral and cellular immunity during vaccination.

Most FDA-approved adjuvants for infectious agents boost humoral but not cellular immunity, and have poorly-understood mechanisms. Stimulator of interferon genes (STING, also known as MITA, MPYS, or ERIS) is an exciting adjuvant target due to its role in cyclic dinucleotide (CDN)-driven anti-viral immunity; however, a major hindrance is STING's cytosolic localization which requires intracellular delivery of its agonists. As a result, STING agonists administered in a soluble form have elicited suboptimal immune responses. Delivery of STING agonists via particle platforms has proven a more successful strategy, but the opportunity for improved formulations and bioactivity remains. In this study we evaluated the adjuvant activity of the potent STING agonist, CDN 3'3'-cGAMP (cGAMP), encapsulated in acid-sensitive acetalated dextran (Ace-DEX) polymeric microparticles (MPs) which passively target antigen-presenting cells for intracellular release. This formulation was superior to all particle delivery systems evaluated and maintained its bioactivity following a sterilizing dose of gamma irradiation. Compared to soluble cGAMP, the Ace-DEX cGAMP MPs enhanced type-I interferon responses nearly 1000-fold in vitro and 50-fold in vivo, caused up to a 104-fold boost in antibody titers, increased Th1-associated responses, and expanded germinal center B cells and memory T cells. Furthermore, the encapsulated cGAMP elicited no observable toxicity in animals and achieved protective immunity against a lethal influenza challenge seven months post-immunization when using CDN adjuvant doses up to 100-fold lower than previous reports. For these reasons, Ace-DEX MP-encapsulated cGAMP represents a potent vaccine adjuvant of humoral and cellular immunity.

Radiation- and Age-Associated Changes in Peripheral Blood Dendritic Cell Populations among Aging Atomic Bomb Survivors in Japan.

Previous immunological studies in atomic bomb survivors have suggested that radiation exposure leads to long-lasting changes, similar to immunological aging observed in T-cell-adaptive immunity. However, to our knowledge, late effects of radiation on dendritic cells (DCs), the key coordinators for activation and differentiation of T cells, have not yet been investigated in humans. In the current study, we hypothesized that numerical and functional decreases would be observed in relationship to radiation dose in circulating conventional DCs (cDCs) and plasmacytoid DCs (pDCs) among 229 Japanese A-bomb survivors. Overall, the evidence did not support this hypothesis, with no overall changes in DCs or functional changes observed with radiation dose. Multivariable regression analysis for radiation dose, age and gender effects revealed that total DC counts as well as subpopulation counts decreased in relationship to increasing age. Further analyses revealed that in women, absolute numbers of pDCs showed significant decreases with radiation dose. A hierarchical clustering analysis of gene expression profiles in DCs after Toll-like receptor stimulation in vitro identified two clusters of participants that differed in age-associated expression levels of genes involved in antigen presentation and cytokine/chemokine production in cDCs. These results suggest that DC counts decrease and expression levels of gene clusters change with age. More than 60 years after radiation exposure, we also observed changes in pDC counts associated with radiation, but only among women.

Sublethal Total Body Irradiation Causes Long-Term Deficits in Thymus Function by Reducing Lymphoid Progenitors.

Total body irradiation (TBI) damages hematopoietic cells in the bone marrow and thymus; however, the long-term effects of irradiation with aging remain unclear. In this study, we found that the impact of radiation on thymopoiesis in mice varied by sex and dose but, overall, thymopoiesis remained suppressed for ≥12 mo after a single exposure. Male and female mice showed a long-term dose-dependent reduction in thymic cKit+ lymphoid progenitors that was maintained throughout life. Damage to hematopoietic stem cells (HSCs) in the bone marrow was dose dependent, with as little as 0.5 Gy causing a significant long-term reduction. In addition, the potential for T lineage commitment was radiation sensitive with aging. Overall, the impact of irradiation on the hematopoietic lineage was more severe in females. In contrast, the rate of decline in thymic epithelial cell numbers with age was radiation-sensitive only in males, and other characteristics including Ccl25 transcription were unaffected. Taken together, these data suggest that long-term suppression of thymopoiesis after sublethal irradiation was primarily due to fewer progenitors in the BM combined with reduced potential for T lineage commitment. A single irradiation dose also caused synchronization of thymopoiesis, with a periodic thymocyte differentiation profile persisting for at least 12 mo postirradiation. This study suggests that the number and capability of HSCs for T cell production can be dramatically and permanently damaged after a single relatively low TBI dose, accelerating aging-associated thymic involution. Our findings may impact evaluation and therapeutic intervention of human TBI events.

Pharmacogenetic Discovery in CALGB (Alliance) 90401 and Mechanistic Validation of a VAC14 Polymorphism that Increases Risk of Docetaxel-Induced Neuropathy.

PURPOSE: Discovery of SNPs that predict a patient's risk of docetaxel-induced neuropathy would enable treatment individualization to maximize efficacy and avoid unnecessary toxicity. The objectives of this analysis were to discover SNPs associated with docetaxel-induced neuropathy and mechanistically validate these associations in preclinical models of drug-induced neuropathy. EXPERIMENTAL DESIGN: A genome-wide association study was conducted in metastatic castrate-resistant prostate cancer patients treated with docetaxel, prednisone and randomized to bevacizumab or placebo on CALGB 90401. SNPs were genotyped on the Illumina HumanHap610-Quad platform followed by rigorous quality control. The inference was conducted on the cumulative dose at occurrence of grade 3+ sensory neuropathy using a cause-specific hazard model that accounted for early treatment discontinuation. Genes with SNPs significantly associated with neuropathy were knocked down in cellular and mouse models of drug-induced neuropathy. RESULTS: A total of 498,081 SNPs were analyzed in 623 Caucasian patients, 50 (8%) of whom experienced grade 3+ neuropathy. The 1,000 SNPs most associated with neuropathy clustered in relevant pathways including neuropathic pain and axonal guidance. An SNP in VAC14 (rs875858) surpassed genome-wide significance (P = 2.12 × 10-8, adjusted P = 5.88 × 10-7). siRNA knockdown of VAC14 in stem cell-derived peripheral neuronal cells increased docetaxel sensitivity as measured by decreased neurite processes (P = 0.0015) and branches (P < 0.0001). Prior to docetaxel treatment, VAC14 heterozygous mice had greater nociceptive sensitivity than wild-type litter mate controls (P = 0.001). CONCLUSIONS: VAC14 should be prioritized for further validation of its potential role as a predictor of docetaxel-induced neuropathy and biomarker for treatment individualization. Clin Cancer Res; 22(19); 4890-900. ©2016 AACR.

Long term effects of radiation exposure on telomere lengths of leukocytes and its associated biomarkers among atomic-bomb survivors.

Ionizing radiation (IR) is a major source of cellular damage and the immediate cellular response to IR has been well characterized. But the long-term impact of IR on cell function and its relationship with aging are not known. Here, we examined the IR effects on telomere length and other biomarkers 50 to 68 years post-exposure (two time points per person) in survivors of the atomic bombing at Hiroshima during WWII. We found that telomere length of leukocytes was inversely correlated with the dose of IR (p=0.008), and this effect was primarily found in survivors who were exposed at younger ages; specifically those <12 years old (p=0.0004). Although a dose-related retardation of telomere shortening with age was observed in the cross-sectional data, longitudinal follow-up after 11 years did not show IR exposure-related alteration of the rate of telomere shortening with age. In addition, IR diminished the associations between telomere length and selected aging biomarkers that were observed in survivors with no dose. These included uric acid metabolism, cytokines, and blood T cell counts. These findings showed long-lasting detrimental effects of IR on telomere length of leukocytes in both dose- and age-at-exposure dependent manner, and on alterations of biomarkers with aging.

An interlaboratory comparison of dosimetry for a multi-institutional radiobiological research project: Observations, problems, solutions and lessons learned.

PURPOSE: An interlaboratory comparison of radiation dosimetry was conducted to determine the accuracy of doses being used experimentally for animal exposures within a large multi-institutional research project. The background and approach to this effort are described and discussed in terms of basic findings, problems and solutions. METHODS: Dosimetry tests were carried out utilizing optically stimulated luminescence (OSL) dosimeters embedded midline into mouse carcasses and thermal luminescence dosimeters (TLD) embedded midline into acrylic phantoms. RESULTS: The effort demonstrated that the majority (4/7) of the laboratories was able to deliver sufficiently accurate exposures having maximum dosing errors of ≤5%. Comparable rates of 'dosimetric compliance' were noted between OSL- and TLD-based tests. Data analysis showed a highly linear relationship between 'measured' and 'target' doses, with errors falling largely between 0 and 20%. Outliers were most notable for OSL-based tests, while multiple tests by 'non-compliant' laboratories using orthovoltage X-rays contributed heavily to the wide variation in dosing errors. CONCLUSIONS: For the dosimetrically non-compliant laboratories, the relatively high rates of dosing errors were problematic, potentially compromising the quality of ongoing radiobiological research. This dosimetry effort proved to be instructive in establishing rigorous reviews of basic dosimetry protocols ensuring that dosing errors were minimized.

Immunization with the Haemophilus ducreyi trimeric autotransporter adhesin DsrA with alum, CpG or imiquimod generates a persistent humoral immune response that recognizes the bacterial surface.

The Ducreyi serum resistance A (DsrA) protein of Haemophilus ducreyi belongs to a large family of multifunctional outer membrane proteins termed trimeric autotransporter adhesins responsible for resistance to the bactericidal activity of human complement (serum resistance), agglutination and adhesion. The ability of DsrA to confer serum resistance and bind extracellular matrix proteins lies in its N-terminal passenger domain. We have previously reported that immunization with a recombinant form of the passenger domain of DsrA, rNT-DsrA, in complete/incomplete Freund's adjuvant, protects against a homologous challenge in swine. We present herein the results of an immunogenicity study in mice aimed at investigating the persistence, type of immune response, and the effect of immunization route and adjuvants on surrogates of protection. Our results indicate that a 20 µg dose of rNT-DsrA administered with alum elicited antisera with comparable bacterial surface reactivity to that obtained with complete/incomplete Freund's adjuvant. At that dose, high titers and bacterial surface reactivity persisted for 211 days after the first immunization. Administration of rNT-DsrA with CpG or imiquimod as adjuvants elicited a humoral response with similar quantity and quality of antibodies (Abs) as seen with Freund's adjuvant. Furthermore, intramuscular administration of rNT-DsrA elicited high-titer Abs with significantly higher reactivity to the bacterial surface than those obtained with subcutaneous immunization. All rNT-DsrA/adjuvant combinations tested, save CpG, elicited a Th2-type response. Taken together, these findings show that a 20 µg dose of rNT-DsrA administered with the adjuvants alum, CpG or imiquimod elicits high-quality Abs with reactivity to the bacterial surface that could protect against an H. ducreyi infection.

Gene expression markers of efficacy and resistance to cetuximab treatment in metastatic colorectal cancer: results from CALGB 80203 (Alliance).

PURPOSE: Formalin-fixed, paraffin-embedded tumor samples from CALGB 80203 were analyzed for expression of EGFR axis-related genes to identify prognostic or predictive biomarkers for cetuximab treatment. PATIENTS AND METHODS: Patients (238 total) with first-line metastatic colorectal cancer (mCRC) were randomized to FOLFOX or FOLFIRI chemotherapy ± cetuximab. qRT-PCR analyses were conducted on tissues from 103 patients at baseline to measure gene expression levels of HER-related genes, including amphiregulin (AREG), betacellulin (BTC), NT5E (CD73), DUSP4, EGF, EGFR, epigen (EPGN), epiregulin (EREG), HBEGF, ERBB2 (HER2), ERBB3 (HER3), ERBB4 (HER4), PHLDA1, and TGFA. The interactions between expression levels and treatment with respect to progression-free survival (PFS) and overall survival (OS) were modeled using multiplicative Cox proportional hazards models. RESULTS: High tumor mRNA levels of HER2 [hazard ratio (HR), 0.64; P = 0.002] and EREG (HR, 0.89; P = 0.016) were prognostic markers associated with longer PFS across all patients. HER3 and CD73 expression levels were identified as potential predictive markers of benefit from cetuximab. In KRAS wild-type (WT) tumors, low HER3 expression was associated with longer OS from cetuximab treatment, whereas high HER3 expression was associated with shorter OS from cetuximab treatment (chemo + cetuximab: HR, 1.15; chemo-only: HR, 0.48; Pinteraction = 0.029). High CD73 expression was associated with longer PFS from cetuximab treatment in patients with KRAS-WT (chemo + cetuximab: HR, 0.91; chemo-only: HR, 1.57; Pinteraction = 0.026) and KRAS-mutant (Mut) tumors (chemo + cetuximab: HR, 0.80; chemo-only: HR, 1.29; P = 0.025). CONCLUSIONS: Gene expression of HER3 and CD73 was identified as a potential predictive marker for cetuximab. These data implicate HER axis signaling and immune modulation as potential mechanisms of cetuximab action and sensitivity.

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.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.