A prospective genome-wide study of prostate cancer metastases reveals association of wnt pathway activation and increased cell cycle proliferation with primary resistance to abiraterone acetate-prednisone.

Background: Genomic aberrations have been identified in metastatic castration-resistant prostate cancer (mCRPC), but molecular predictors of resistance to abiraterone acetate/prednisone (AA/P) treatment are not known. Patients and methods: In a prospective clinical trial, mCRPC patients underwent whole-exome sequencing (n = 82) and RNA sequencing (n = 75) of metastatic biopsies before initiating AA/P with the objective of identifying genomic alterations associated with resistance to AA/P. Primary resistance was determined at 12 weeks of treatment using criteria for progression that included serum prostate-specific antigen measurement, bone and computerized tomography imaging and symptom assessments. Acquired resistance was determined using the end point of time to treatment change (TTTC), defined as time from enrollment until change in treatment from progressive disease. Associations of genomic and transcriptomic alterations with primary resistance were determined using logistic regression, Fisher's exact test, single and multivariate analyses. Cox regression models were utilized for determining association of genomic and transcriptomic alterations with TTTC. Results: At 12 weeks, 32 patients in the cohort had progressed (nonresponders). Median study follow-up was 32.1 months by which time 58 patients had switched treatments due to progression. Median TTTC was 10.1 months (interquartile range: 4.4-24.1). Genes in the Wnt/ß-catenin pathway were more frequently mutated and negative regulators of Wnt/ß-catenin signaling were more frequently deleted or displayed reduced mRNA expression in nonresponders. Additionally, mRNA expression of cell cycle regulatory genes was increased in nonresponders. In multivariate models, increased cell cycle proliferation scores (≥ 50) were associated with shorter TTTC (hazard ratio = 2.11, 95% confidence interval: 1.17-3.80; P = 0.01). Conclusions: Wnt/ß-catenin pathway activation and increased cell cycle progression scores can serve as molecular markers for predicting resistance to AA/P therapy.

New insights into the pharmacogenomics of antidepressant response from the GENDEP and STAR*D studies: rare variant analysis and high-density imputation.

Genome-wide association studies have generally failed to identify polymorphisms associated with antidepressant response. Possible reasons include limited coverage of genetic variants that this study tried to address by exome genotyping and dense imputation. A meta-analysis of Genome-Based Therapeutic Drugs for Depression (GENDEP) and Sequenced Treatment Alternatives to Relieve Depression (STAR*D) studies was performed at the single-nucleotide polymorphism (SNP), gene and pathway levels. Coverage of genetic variants was increased compared with previous studies by adding exome genotypes to previously available genome-wide data and using the Haplotype Reference Consortium panel for imputation. Standard quality control was applied. Phenotypes were symptom improvement and remission after 12 weeks of antidepressant treatment. Significant findings were investigated in NEWMEDS consortium samples and Pharmacogenomic Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS) for replication. A total of 7062 950 SNPs were analyzed in GENDEP (n=738) and STAR*D (n=1409). rs116692768 (P=1.80e-08, ITGA9 (integrin α9)) and rs76191705 (P=2.59e-08, NRXN3 (neurexin 3)) were significantly associated with symptom improvement during citalopram/escitalopram treatment. At the gene level, no consistent effect was found. At the pathway level, the Gene Ontology (GO) terms GO: 0005694 (chromosome) and GO: 0044427 (chromosomal part) were associated with improvement (corrected P=0.007 and 0.045, respectively). The association between rs116692768 and symptom improvement was replicated in PGRN-AMPS (P=0.047), whereas rs76191705 was not. The two SNPs did not replicate in NEWMEDS. ITGA9 codes for a membrane receptor for neurotrophins and NRXN3 is a transmembrane neuronal adhesion receptor involved in synaptic differentiation. Despite their meaningful biological rationale for being involved in antidepressant effect, replication was partial. Further studies may help in clarifying their role.

TSPAN5, ERICH3 and selective serotonin reuptake inhibitors in major depressive disorder: pharmacometabolomics-informed pharmacogenomics.

Millions of patients suffer from major depressive disorder (MDD), but many do not respond to selective serotonin reuptake inhibitor (SSRI) therapy. We used a pharmacometabolomics-informed pharmacogenomics research strategy to identify genes associated with metabolites that were related to SSRI response. Specifically, 306 MDD patients were treated with citalopram or escitalopram and blood was drawn at baseline, 4 and 8 weeks for blood drug levels, genome-wide single nucleotide polymorphism (SNP) genotyping and metabolomic analyses. SSRI treatment decreased plasma serotonin concentrations (P<0.0001). Baseline and plasma serotonin concentration changes were associated with clinical outcomes (P<0.05). Therefore, baseline and serotonin concentration changes were used as phenotypes for genome-wide association studies (GWAS). GWAS for baseline plasma serotonin concentrations revealed a genome-wide significant (P=7.84E-09) SNP cluster on chromosome four 5' of TSPAN5 and a cluster across ERICH3 on chromosome one (P=9.28E-08) that were also observed during GWAS for change in serotonin at 4 (P=5.6E-08 and P=7.54E-07, respectively) and 8 weeks (P=1.25E-06 and P=3.99E-07, respectively). The SNPs on chromosome four were expression quantitative trait loci for TSPAN5. Knockdown (KD) and overexpression (OE) of TSPAN5 in a neuroblastoma cell line significantly altered the expression of serotonin pathway genes (TPH1, TPH2, DDC and MAOA). Chromosome one SNPs included two ERICH3 nonsynonymous SNPs that resulted in accelerated proteasome-mediated degradation. In addition, ERICH3 and TSPAN5 KD and OE altered media serotonin concentrations. Application of a pharmacometabolomics-informed pharmacogenomic research strategy, followed by functional validation, indicated that TSPAN5 and ERICH3 are associated with plasma serotonin concentrations and may have a role in SSRI treatment outcomes.

Pharmacogenomics of selective serotonin reuptake inhibitor treatment for major depressive disorder: genome-wide associations and functional genomics.

A genome-wide association (GWA) study of treatment outcomes (response and remission) of selective serotonin reuptake inhibitors (SSRIs) was conducted using 529 subjects with major depressive disorder. While no SNP associations reached the genome-wide level of significance, 14 SNPs of interest were identified for functional analysis. The rs11144870 SNP in the riboflavin kinase (RFK) gene on chromosome 9 was associated with 8-week treatment response (odds ratio (OR)=0.42, P=1.04 × 10⁻6). The rs915120 SNP in the G protein-coupled receptor kinase 5 (GRK5) gene on chromosome 10 was associated with 8-week remission (OR=0.50, P=1.15 × 10⁻5). Both SNPs were shown to influence transcription by a reporter gene assay and to alter nuclear protein binding using an electrophoretic mobility shift assay. This report represents an example of joining functional genomics with traditional GWA study results derived from a GWA analysis of SSRI treatment outcomes. The goal of this analytical strategy is to provide insights into the potential relevance of biologically plausible observed associations.

Multivariate models to detect genomic signatures for a class of drugs: application to thiopurines pharmacogenomics.

Often, analysis for pharmacogenomic studies involving multiple drugs from the same class is completed by analyzing each drug individually for association with genomic variation. However, by completing the analysis of each drug individually, we may be losing valuable information. When studying multiple drugs from the same drug class, one may wish to determine genomic variation that explains the difference in response between individuals for the drug class, as opposed to each individual drug. Therefore, we have developed a multivariate model to assess whether genomic variation impacts a class of drugs. In addition to determine genomic effects that are similar for the drugs, we will also be able to determine genomic effects that differ between the drugs (that is, interaction). We will illustrate the utility of this multivariate model for cytotoxicity and genomic data collected on the Coriell Human Variation Panel for the class of anti-purine metabolites (6-mercaptopurine and 6-thioguanine).

Catechol O-methyltransferase pharmacogenomics and selective serotonin reuptake inhibitor response.

We applied a systematic pharmacogenetic approach to investigate the role of genetic variation in the gene encoding catechol O-methyltransferase (COMT) in individual variation in selective serotonin reuptake inhibitor (SSRI) response among depressed patients. In all, 23 single-nucleotide polymorphisms (SNPs) in COMT were genotyped using DNA from the Sequenced Treatment Alternatives to Relieve Depression (STAR(*)D) study (N=1914). One SNP, rs13306278, located in the distal promoter region of COMT, showed significant association with remission in White non-Hispanic (WNH) subjects (P=0.038). Electromobility shift assay for rs13306278 showed alternation in the ability of the variant sequence to bind nuclear proteins. A replication study was performed using samples from the Mayo Clinic Pharmacogenetics Research Network Citalopram/Escitalopram Pharmacogenomic study (N=422) that demonstrated a similar trend for association. Our findings suggest that novel genetic markers in the COMT distal promoter may influence SSRI response phenotypes.

Catechol-O-methyltransferase: thermolabile enzyme in erythrocytes of subjects homozygous for allele for low activity.

Low catechol-O-methyltransferase (COMT) activity (less than 8 units per milliliter) in the human erythrocyte is inherited as an autosomal recessive trait (COMTL). The average half-life of COMT in erythrocyte lysates incubated at 48 degrees C was significantly shorter in lysates from three subjects with low enzyme activity than in lysates from three subjects with high enzyme activity (12.5 +/- 0.9 minutes compared with 21.2 +/- 1.4 minutes, P less than .01). When the ratios of COMT activities in lysates heated at 48 degrees C for 15 minutes to enzyme activities in unheated samples were used as a measure of enzyme thermostability in blood samples from 316 randomly selected subjects, the ratios were significantly less for subjects with low enzyme activity than for subjects with higher enzyme activity. The presense of thermolabile COMT in blood of individuals homozygous for COMTL raises the possibility that the locus COMT may represent the structural gene for the human enzyme.

Proportional release of norepinephrine and dopamine- -hydroxylase from sympathetic nerves.

Dopamine-beta-hydroxylase(DBH), the enzyme that catalyzes the conversion of dopamine to norepinephrine, is localized in the vesicles containing catecholamine in sympathetic nerves. This enzyme is released with norepinephrine when the nerves to the guinea pig vas deferens are stimulated in vitro, and the amount of enzyme discharged increases as the length of stimulation periods increases. The amount of DBH released is proportional to the amount of norepinephrine released, and the ratio of norepinephrine to DBH discharged into the incubation medium is similar to that in the soluble portion of the contents of the synaptic vesicles from the vas deferens. These data are compatible with the release of the neurotransmitter norepinephrine and DBH from symnpathetic nerves by a process of exocytosis.

Serum dopamine-beta-hydroxylase activity: sibling-sibling correlation.

Dopamine-beta-hydroxylase activity is released into the blood with catecholamines from the adrenal medulla and sympathetic nerves. This enzyme activity has been measured in the blood of 317 normal children and 227 normal adults. A significant sibling-sibling correlation of serum dopamine-beta-hydroxylase values was found in the 94 sibling pairs tested. Frequency distributions of serum enzyme values in both children and adults suggest the existence of two populations with regard to serum activity of this enzyme.

SLC6A4 variation and citalopram response.

The influence of genetic variations in SLC6A4 (serotonin transporter gene) on citalopram treatment of depression using the Sequenced Treatment to Relieve Depression (STAR*D) sample was assessed. Of primary interest were three previously studied polymorphisms: 1) the VNTR variation of the second intron, 2) the indel promoter polymorphism (5HTTLPR or SERT), and 3) a single nucleotide polymorphism (SNP) rs25531. Additionally, SLC6A4 was resequenced to identify new SNPs for exploratory analyses. DNA from 1914 subjects in the STAR*D study were genotyped for the intron 2 VNTR region, the indel promoter polymorphism, and rs25531. Associations of these variants with remission of depressive symptoms were evaluated following citalopram treatment. In white non-Hispanic subjects, variations in the intron 2 VNTR (point-wise P = 0.041) and the indel promoter polymorphism (point-wise P = 0.039) were associated with remission following treatment with citalopram. The haplotype composed of the three candidate loci was also associated with remission, with a global p-value of 0.040 and a maximum statistic simulation p-value of 0.0031 for the S-a-12 haplotype, under a dominant model. One SNP identified through re-sequencing the SLC6A4 gene, Intron7-83-TC, showed point-wise evidence of association, which did not remain significant after correction for the number of SNPs evaluated in this exploratory analysis. No associations between these SLC6A4 variations and remission were found in the white Hispanic or black subjects. These findings suggest that multiple variations in the SLC6A4 gene are associated with remission in white non-Hispanic depressed adults treated with citalopram. The mechanism of action of these variants remains to be determined.

Sulfotransferase gene copy number variation: pharmacogenetics and function.

Pharmacogenetics is the study of the role of inheritance in variation to drug response. Drug response phenotypes can vary from adverse drug reactions at one end of the spectrum to equally serious lack of the desired effect of drug therapy at the other. Many of the current important examples of pharmacogenetics involve inherited variation in drug metabolism. Sulfate conjugation catalyzed by cytosolic sulfotransferase (SULT) enzymes, particularly SULT1A1, is a major pathway for drug metabolism in humans. Pharmacogenetic studies of SULT1A1 began over a quarter of a century ago and have advanced from biochemical genetic experiments to include cDNA and gene cloning, gene resequencing, and functional studies of the effects of single nucleotide polymorphisms (SNPs). SNP genotyping, in turn, led to the discovery of functionally important copy number variations (CNVs) in the SULT1A1 gene. This review will briefly describe the evolution of our understanding of SULT1A1 pharmacogenetics and CNV, as well as challenges involved in utilizing both SNP and CNV data in an attempt to predict SULT1A1 function. SULT1A1 represents one example of the potential importance of CNV for the evolving disciplines of pharmacogenetics and pharmacogenomics.

The pharmacogenetics research network: from SNP discovery to clinical drug response.

The NIH Pharmacogenetics Research Network (PGRN) is a collaborative group of investigators with a wide range of research interests, but all attempting to correlate drug response with genetic variation. Several research groups concentrate on drugs used to treat specific medical disorders (asthma, depression, cardiovascular disease, addiction of nicotine, and cancer), whereas others are focused on specific groups of proteins that interact with drugs (membrane transporters and phase II drug-metabolizing enzymes). The diverse scientific information is stored and annotated in a publicly accessible knowledge base, the Pharmacogenetics and Pharmacogenomics Knowledge base (PharmGKB). This report highlights selected achievements and scientific approaches as well as hypotheses about future directions of each of the groups within the PGRN. Seven major topics are included: informatics (PharmGKB), cardiovascular, pulmonary, addiction, cancer, transport, and metabolism.

Inheritance of low serum immunoglobulin D.

Previous studies have shown that log IgD levels in normal individuals are distributed in a nonunimodal manner. Therefore, in this study we tested whether inheritance might play a role in determination of IgD levels. IgD levels were measured in serum or plasma from 301 randomly selected children aged 6-18 yr, 245 consecutive adult blood donors, and 134 first-degree relatives of subjects with low IgD levels. Comparison of serum and plasma from five individuals revealed no difference, so the two were used interchangeably. The distributions of log IgD levels in randomly selected populations of both adults and children were nonunimodal with nadirs at 2.15 IU/ml. In both of these randomly selected populations, 13-14% of the subjects had low IgD values (<2.15 IU/ml). In addition, there was a significant sibling-sibling correlation of log IgD values (r = 0.56, n = 72, P <0.01). Because of the nonunimodality of the frequency distribution histogram for IgD values and because of the familial aggregation of these values, the study was extended to include first-degree relatives of subjects with low plasma IgD. Blood samples from 92% of living first-degree relatives, 134 individuals, were analyzed for their level of IgD, and the results of segregation and pedigree analyses of these data were compatible with autosomal recessive inheritance of an allele for low plasma IgD levels. IgD values in plasma from siblings of probands for low IgD were also non-unimodal in distribution with a nadir at congruent with2.15 IU/ml. The results suggest that there is autosomal recessive inheritance of an allele for low plasma IgD.

Human thiopurine methyltransferase pharmacogenetics. Kindred with a terminal exon splice junction mutation that results in loss of activity.

Thiopurine methyltransferase (TPMT) catalyzes S-methylation of thiopurine drugs such as 6-mercaptopurine. Large variations in levels of TPMT activity in human tissue can result from a common genetic polymorphism with a series of alleles for low activity. This polymorphism is an important factor responsible for large individual variations in thiopurine toxicity and therapeutic efficacy. We now report a new variant allele, TPMT*4, that contains a G--> A transition that disrupts the intron/exon acceptor splice junction at the final 3' nucleotide of intron 9, the terminal intron of the TPMT gene. This new allele cosegregated within an extended kindred with reduced TPMT activity. We attempted to determine the mechanism(s) by which the presence of TPMT*4 might result in low enzyme activity. Although very few mature transcripts derived from allele TPMT*4 were detected, the mutation did lead to generation of at least two aberrant mRNA species. The first resulted from use of a novel splice site located one nucleotide 3' downstream from the original splice junction. That mRNA species contained a single nucleotide deletion and a frameshift within exon 10, the terminal exon of the gene. The second novel mRNA species resulted from activation of a cryptic splice site located within intron 9, leading to inclusion of 330 nucleotides of intron sequence. That sequence contained a premature translation termination codon. TPMT*4 is the first reported allele for low TPMT activity as a result of a mutation within an intron. These observations also provide insight into mechanisms of mRNA processing after disruption of a terminal exon splice junction.

Electronic Medical Record-Integrated Pharmacogenomics and Related Clinical Decision Support Concepts.

Advances in pharmacogenomics (PGx) have the potential to transform healthcare by allowing precision medicine to become a reality. However, PGx knowledge is new, complex, and evolving, and relying on the cognition of clinicians alone is insufficient for clinical implementation. Integrating clinical decision support (CDS) tools in the electronic health record (EHR) is critical for translating PGx into clinical practice. Herein, we review current strategies to implement PGx using EHR-CDS functionalities.

Genetic Breast Cancer Susceptibility Variants and Prognosis in the Prospectively Randomized SUCCESS A Study.

Large-scale genotyping studies have identified over 70 single nucleotide polymorphisms (SNPs) associated with breast cancer (BC) risk. However, knowledge regarding genetic risk factors associated with the prognosis is limited. The aim of this study was therefore to investigate the prognostic effect of nine known breast cancer risk SNPs. BC patients (n = 1687) randomly sampled in an adjuvant, randomized phase III trial (SUCCESS A study) were genotyped for nine BC risk SNPs: rs17468277 (CASP8) , rs2981582 (FGFR2) , rs13281615(8q24), rs3817198 (LSP1) , rs889312 (MAP3K1) , rs3803662 (TOX3) , rs13387042(2q35), rs4973768 (SLC4A7) , rs6504950 (COX11) . Cox proportional hazards models were used to test the SNPs' association with overall survival (OS) and progression-free survival (PFS). Additional analyses were carried out for molecular subgroups. rs3817198 in LSP1 (lymphocyte-specific protein 1) was the only SNP that significantly influenced OS (p = 0.01) and PFS (p < 0.01) in the likelihood ratio test comparing the genetic survival model with the clinical survival model. In the molecular subgroups, triple-negative patients with two minor alleles in rs3817198 had a much better prognosis relative to OS (adjusted HR 0.03; 95% CI 0.002 - 0.279) and PFS (HR 0.09; 95% CI 0.02 - 0.36) than patients with the common alleles. The same effect on PFS was shown for patients with luminal A tumors (HR 0.19; 95% CI 0.05 - 0.84), whereas patients with luminal B tumors had a poorer PFS with two minor alleles (HR 2.13; 95% CI 1.02 - 4.40). The variant in rs3817198 has a prognostic effect particularly in the subgroup of patients with triple-negative BC, suggesting a possible link with immunomodulation and BC.

The Pharmacogenomics Research Network Translational Pharmacogenetics Program: Outcomes and Metrics of Pharmacogenetic Implementations Across Diverse Healthcare Systems.

Numerous pharmacogenetic clinical guidelines and recommendations have been published, but barriers have hindered the clinical implementation of pharmacogenetics. The Translational Pharmacogenetics Program (TPP) of the National Institutes of Health (NIH) Pharmacogenomics Research Network was established in 2011 to catalog and contribute to the development of pharmacogenetic implementations at eight US healthcare systems, with the goal to disseminate real-world solutions for the barriers to clinical pharmacogenetic implementation. The TPP collected and normalized pharmacogenetic implementation metrics through June 2015, including gene-drug pairs implemented, interpretations of alleles and diplotypes, numbers of tests performed and actionable results, and workflow diagrams. TPP participant institutions developed diverse solutions to overcome many barriers, but the use of Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines provided some consistency among the institutions. The TPP also collected some pharmacogenetic implementation outcomes (scientific, educational, financial, and informatics), which may inform healthcare systems seeking to implement their own pharmacogenetic testing programs.

Genomewide Association Studies in Pharmacogenomics: Meeting Report of the NIH Pharmacogenomics Research Network-RIKEN (PGRN-RIKEN) Collaboration.

Genomewide association studies (GWAS) have resulted in the identification of many heritable genetic factors that underlie risk for human disease or variation in physiologic traits. In contrast, there are fewer GWAS of drug response phenotypes, despite extensive unexplained interindividual variability. To address this urgent need, the NIH Pharmacogenomics Research Network (PGRN) and the Center for Integrative Medical Sciences (IMS) at RIKEN support a collaboration, PGRN-RIKEN, with the goal of accelerating GWAS of drug response phenotypes.

Genetic polymorphism for human platelet thermostable phenol sulfotransferase (TS PST) activity.

Platelet TS PST basal activity and thermal stability were measured in blood samples from 237 individuals in 50 nuclear families. Significant correlations were found among first degree relatives, confirming the previously reported familial aggregation of TS PST basal activity and thermal stability. Commingling analysis of basal TS PST activity provided evidence for multiple component distributions, and after transformation to remove skewness, segregation analysis supported a major gene hypothesis. For TS PST thermal stability, commingling analysis also provided evidence for multiple component distributions. However, segregation analyses were equivocal with regard to the presence of a major gene for thermal stability, since support for a major gene model depended on skewness. Bivariate commingling analysis, which examined thermal stability by simultaneously considering basal activity and activity after heating, suggested that genotypes, as defined by the inferred component distributions for TS PST activity, differ in thermal stability. A three-allele model is proposed as one hypothesis that may account for the combined results of basal activity and thermal stability. The results of this study indicate that a major gene polymorphism in conjunction with polygenic inheritance plays an important role in the regulation of both level of activity and thermal stability of this important drug-metabolizing enzyme in humans.