Efavirenz (EFV) based antiretroviral therapy is expanding worldwide. However discontinuation of EFV containing regimens is common in some patients, particularly black patients, due most often to neuropsychiatric side effects. These adverse drug effects often result in premature drug discontinuation, as well as considerable morbidity.

To better address the problem of drug resistance during cancer chemotherapy and explore the possibility of manipulating drug response phenotypes, we developed a network-based phenotype mapping approach (P-Map) to identify gene candidates that upon perturbed can alter sensitivity to drugs. We used basal transcriptomics data from a panel of human lymphoblastoid cell lines (LCL) to infer drug response networks (DRNs) that are responsible for conferring response phenotypes for anthracycline and taxane, two common anticancer agents use in clinics. We further tested selected gene candidates that interact with phenotypic differentially expressed genes (PDEGs), which are up-regulated genes in LCL for a given class of drug response phenotype in triple-negative breast cancer (TNBC) cells. Our results indicate that it is possible to manipulate a drug response phenotype, from resistant to sensitive or vice versa, by perturbing gene candidates in DRNs and suggest plausible mechanisms regulating directionality of drug response sensitivity. More important, the current work highlights a new way to formulate systems-based therapeutic design: supplementing therapeutics that aim to target disease culprits with phenotypic modulators capable of altering DRN properties with the goal to re-sensitize resistant phenotypes.

Patient-derived xenografts (PDXs) are increasingly used in cancer research as a tool to inform cancer biology and drug response. Most available breast cancer PDXs have been generated in the metastatic setting. However, in the setting of operable breast cancer, PDX models both sensitive and resistant to chemotherapy are needed for drug development and prospective data are lacking regarding the clinical and molecular characteristics associated with PDX take rate in this setting.

Breast cancer patients with residual disease after neoadjuvant chemotherapy (NAC) have increased recurrence risk. Molecular characterization, knowledge of NAC response, and simultaneous generation of patient-derived xenografts (PDXs) may accelerate drug development. However, the feasibility of this approach is unknown.

Cytochrome P450 enzymes including CYP2C19 and CYP2D6 are important for antidepressant metabolism and polymorphisms of these genes have been determined to predict metabolite levels. Nonetheless, more evidence is needed to understand the impact of genetic variations on antidepressant response. In this study, individual clinical and genetic data from 13 studies of European and East Asian ancestry populations were collected. The antidepressant response was clinically assessed as remission and percentage improvement. Imputed genotype was used to translate genetic polymorphisms to metabolic phenotypes (poor, intermediate, normal, and rapid+ultrarapid) of CYP2C19 and CYP2D6. The association of CYP2C19 and CYP2D6 metabolic phenotypes with treatment response was examined using normal metabolizers as the reference. Among 5843 depression patients, a higher remission rate was found in CYP2C19 poor metabolizers compared to normal metabolizers at nominal significance but did not survive after multiple testing correction (OR=1.46, 95% CI [1.03, 2.06], p=0.033, heterogeneity I2=0%, subgroup difference p=0.72). No metabolic phenotype was associated with percentage improvement from baseline. After stratifying by antidepressants primarily metabolized by CYP2C19 and CYP2D6, no association was found between metabolic phenotypes and antidepressant response. Metabolic phenotypes showed differences in frequency, but not effect, between European- and East Asian-ancestry studies. In conclusion, metabolic phenotypes imputed from genetic variants using genotype were not associated with antidepressant response. CYP2C19 poor metabolizers could potentially contribute to antidepressant efficacy with more evidence needed. CYP2D6 structural variants cannot be imputed from genotype data, limiting inference of pharmacogenetic effects. Sequencing and targeted pharmacogenetic testing, alongside information on side effects, antidepressant dosage, depression measures, and diverse ancestry studies, would more fully capture the influence of metabolic phenotypes.

Developing patient derived models from individual tumors that capture the biological heterogeneity and mutation landscape in advanced prostate cancer is challenging, but essential for understanding tumor progression and delivery of personalized therapy in metastatic castrate resistant prostate cancer stage. To demonstrate the feasibility of developing patient derived xenograft models in this stage, we present a case study wherein xenografts were derived from cancer metastases in a patient progressing on androgen deprivation therapy and prior to initiating pre-chemotherapy enzalutamide treatment. Tissue biopsies from a metastatic rib lesion were obtained for sequencing before and after initiating enzalutamide treatment over a twelve-week period and also implanted subcutaneously as well as under the renal capsule in immuno-deficient mice. The genome and transcriptome landscapes of xenografts and the original patient tumor tissues were compared by performing whole exome and transcriptome sequencing of the metastatic tumor tissues and the xenografts at both time points. After comparing the somatic mutations, copy number variations, gene fusions and gene expression we found that the patient's genomic and transcriptomic alterations were preserved in the patient derived xenografts with high fidelity. These xenograft models provide an opportunity for predicting efficacy of existing and potentially novel drugs that is based on individual metastatic tumor expression signature and molecular pharmacology for delivery of precision medicine.

Oral anticoagulant (OAC) therapy has been the main treatment approach for stroke prevention for decades. Warfarin is the most widely prescribed OAC in the United States, but is difficult to manage due to variability in dose requirements across individuals. Pharmacogenomics may mitigate risk concerns related to warfarin use by fostering the opportunity to facilitate individualized medicine approaches to warfarin treatment (e.g., genome-guided dosing). While various economic evaluations exist examining the cost-effectiveness of pharmacogenomics testing for warfarin, few observational studies exist to support these studies, with even fewer using genotype as the main exposure of interest. We examined a cohort of individuals initiating warfarin therapy between 2004 and 2017 and examined bleeding and cost outcomes for the year following initiation using Mayo Clinic's billing and administrative data, as well the Mayo Clinic Rochester Cost Data Warehouse. Analyses included descriptive summaries, comparison of characteristics across exposure groups, reporting of crude outcomes, and multivariate analyses. We included N = 1,143 patients for analyses. Just over a third of our study population (34.9%) carried a warfarin-sensitive phenotype. Sensitive individuals differed in their baseline characteristics by being of older age and having a higher number of comorbid conditions; myocardial infarction, diabetes, and cancer in particular. The occurrence of bleeding events was not significantly different across exposure groups. No significant differences across exposure groups existed in either the likelihood of incurring all-cause healthcare costs or in the magnitude of those costs. Warfarin-sensitive individuals were no more likely to utilize cardiovascular-related healthcare services; however, they had lower total and inpatient cardiovascular-related costs compared to warfarin-insensitive patients. No significant differences existed in any other categories of costs. We found limited evidence that warfarin-sensitive individuals have different healthcare spending than warfarin-insensitive individuals. Additional real-world studies are needed to support the traditional economic evaluations currently existing in the literature.

Alcohol consumption can increase circulating levels of fibroblast growth factor 21 (FGF21). The effects of FGF21 in the central nervous system are associated with the regulation of catecholamines, neurotransmitters that play a crucial role in reward pathways. This study aims to identify genetic variants associated with FGF21 levels and evaluate their functional role in alcohol use disorder (AUD).

We have utilized a two-stage study design to search for SNPs associated with the survival of breast cancer patients treated with adjuvant chemotherapy. Our initial GWS data set consisted of 805 Finnish breast cancer cases (360 treated with adjuvant chemotherapy). The top 39 SNPs from this stage were analyzed in three independent data sets: iCOGS (n=6720 chemotherapy-treated cases), SUCCESS-A (n=3596), and POSH (n=518). Two SNPs were successfully validated: rs6500843 (any chemotherapy; per-allele HR 1.16, 95% C.I. 1.08-1.26, p=0.0001, p(adjusted)=0.0091), and rs11155012 (anthracycline therapy; per-allele HR 1.21, 95% C.I. 1.08-1.35, p=0.0010, p(adjusted)=0.0270). The SNP rs6500843 was found to specifically interact with adjuvant chemotherapy, independently of standard prognostic markers (p(interaction)=0.0009), with the rs6500843-GG genotype corresponding to the highest hazard among chemotherapy-treated cases (HR 1.47, 95% C.I. 1.20-1.80). Upon trans-eQTL analysis of public microarray data, the rs6500843 locus was found to associate with the expression of a group of genes involved in cell cycle control, notably AURKA, the expression of which also exhibited differential prognostic value between chemotherapy-treated and untreated cases in our analysis of microarray data. Based on previously published information, we propose that the eQTL genes may be connected to the rs6500843 locus via a RBFOX1-FOXM1 -mediated regulatory pathway.

When sequencing blood and tumor samples to identify targetable somatic variants for cancer therapy, clinically relevant germline variants may be uncovered. We evaluated the prevalence of deleterious germline variants in cancer susceptibility genes in women with breast cancer referred for neoadjuvant chemotherapy and returned clinically actionable results to patients. Exome sequencing was performed on blood samples from women with invasive breast cancer referred for neoadjuvant chemotherapy. Germline variants within 142 hereditary cancer susceptibility genes were filtered and reviewed for pathogenicity. Return of results was offered to patients with deleterious variants in actionable genes if they were not aware of their result through clinical testing. 124 patients were enrolled (median age 51) with the following subtypes: triple negative (n = 43, 34.7%), HER2+ (n = 37, 29.8%), luminal B (n = 31, 25%), and luminal A (n = 13, 10.5%). Twenty-eight deleterious variants were identified in 26/124 (21.0%) patients in the following genes: ATM (n = 3), BLM (n = 1), BRCA1 (n = 4), BRCA2 (n = 8), CHEK2 (n = 2), FANCA (n = 1), FANCI (n = 1), FANCL (n = 1), FANCM (n = 1), FH (n = 1), MLH3 (n = 1), MUTYH (n = 2), PALB2 (n = 1), and WRN (n = 1). 121/124 (97.6%) patients consented to return of research results. Thirteen (10.5%) had actionable variants, including four that were returned to patients and led to changes in medical management. Deleterious variants in cancer susceptibility genes are highly prevalent in patients with invasive breast cancer referred for neoadjuvant chemotherapy undergoing exome sequencing. Detection of these variants impacts medical management.

We demonstrate that model-based unsupervised learning can uniquely discriminate single-cell subpopulations by their gene expression distributions, which in turn allow us to identify specific genes for focused functional studies. This method was applied to MDA-MB-231 breast cancer cells treated with the antidiabetic drug metformin, which is being repurposed for treatment of triple-negative breast cancer. Unsupervised learning identified a cluster of metformin-treated cells characterized by a significant suppression of 230 genes (p-value < 2E-16). This analysis corroborates known studies of metformin action: a) pathway analysis indicated known mechanisms related to metformin action, including the citric acid (TCA) cycle, oxidative phosphorylation, and mitochondrial dysfunction (p-value < 1E-9); b) 70% of these 230 genes were functionally implicated in metformin response; c) among remaining lesser functionally-studied genes for metformin-response was CDC42, down-regulated in breast cancer treated with metformin. However, CDC42's mechanisms in metformin response remained unclear. Our functional studies showed that CDC42 was involved in metformin-induced inhibition of cell proliferation and cell migration mediated through an AMPK-independent mechanism. Our results points to 230 genes that might serve as metformin response signatures, which needs to be tested in patients treated with metformin and, further investigation of CDC42 and AMPK-independence's role in metformin's anticancer mechanisms.

Breast cancer is the most common invasive cancer among women. Currently, there are only a few models used for therapy selection, and they are often poor predictors of therapeutic response or take months to set up and assay. In this report, we introduce a microfluidic OrganoPlate® platform for extracellular matrix (ECM) embedded tumor culture under perfusion as an initial study designed to investigate the feasibility of adapting this technology for therapy selection.

Biomedical research is undergoing rapid change, with the development of a series of analytical omics techniques that are capable of generating Biomedical Big Data. These developments provide an unprecedented opportunity to gain novel insight into disease pathophysiology and mechanisms of drug action and response-but they also present significant challenges. Pharmacogenomics is a discipline within Clinical Pharmacology that has been at the forefront in defining, taking advantage of, and dealing with the opportunities and challenges of this aspect of the Post-Genome Project world. This overview will describe the evolution of germline pharmacogenomic research strategies as we have moved from an era of candidate genes to agnostic genome-wide association studies (GWAS) coupled with the functional and mechanistic pursuit of GWAS signals. Germline pharmacogenomic studies of breast cancer endocrine therapy will be used to illustrate research strategies that are being applied broadly to omics studies of drug response phenotypes.

There are currently no reliable approaches for correctly identifying which patients with major depressive disorder (MDD) will respond well to antidepressant therapy. However, recent genetic advances suggest that Polygenic Risk Scores (PRS) could allow MDD patients to be stratified for antidepressant response. We used PRS for MDD and PRS for neuroticism as putative predictors of antidepressant response within three treatment cohorts: The Genome-based Therapeutic Drugs for Depression (GENDEP) cohort, and 2 sub-cohorts from the Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomics Study PRGN-AMPS (total patient number = 760). Results across cohorts were combined via meta-analysis within a random effects model. Overall, PRS for MDD and neuroticism did not significantly predict antidepressant response but there was a consistent direction of effect, whereby greater genetic loading for both MDD (best MDD result, p < 5*10-5 MDD-PRS at 4 weeks, β = -0.019, S.E = 0.008, p = 0.01) and neuroticism (best neuroticism result, p < 0.1 neuroticism-PRS at 8 weeks, β = -0.017, S.E = 0.008, p = 0.03) were associated with less favourable response. We conclude that the PRS approach may offer some promise for treatment stratification in MDD and should now be assessed within larger clinical cohorts.

Programmed death ligand-1 (PD-L1 or B7-H1) is well known for its role in immune checkpoint regulation, but its function inside the tumor cells has rarely been explored. Here we report that nuclear PD-L1 is important for cancer cell sister chromatid cohesion. We found that depletion of PD-L1 suppresses cancer cell proliferation, colony formation in vitro, and tumor growth in vivo in immune-deficient NSG mice independent of its role in immune checkpoint. Specifically, PD-L1 functions as a subunit of the cohesin complex, and its deficiency leads to formation of multinucleated cells and causes a defect in sister chromatid cohesion. Mechanistically, PD-L1 compensates for the loss of Sororin, whose expression is suppressed in cancer cells overexpressing PD-L1. PD-L1 competes with Wing Apart-Like (WAPL) for binding to PDS5B, and secures proper sister chromatid cohesion and segregation. Our findings suggest an important role for nuclear PD-L1 in cancer cells independent of its function in immune checkpoint.

Survival rates of patients with metastatic castration-resistant prostate cancer (mCRPC) are low due to lack of response or acquired resistance to available therapies, such as abiraterone (Abi). A better understanding of the underlying molecular mechanisms is needed to identify effective targets to overcome resistance. Given the complexity of the transcriptional dynamics in cells, differential gene expression analysis of bulk transcriptomics data cannot provide sufficient detailed insights into resistance mechanisms. Incorporating network structures could overcome this limitation to provide a global and functional perspective of Abi resistance in mCRPC. Here, we developed TraRe, a computational method using sparse Bayesian models to examine phenotypically driven transcriptional mechanistic differences at three distinct levels: transcriptional networks, specific regulons, and individual transcription factors (TF). TraRe was applied to transcriptomic data from 46 patients with mCRPC with Abi-response clinical data and uncovered abrogated immune response transcriptional modules that showed strong differential regulation in Abi-responsive compared with Abi-resistant patients. These modules were replicated in an independent mCRPC study. Furthermore, key rewiring predictions and their associated TFs were experimentally validated in two prostate cancer cell lines with different Abi-resistance features. Among them, ELK3, MXD1, and MYB played a differential role in cell survival in Abi-sensitive and Abi-resistant cells. Moreover, ELK3 regulated cell migration capacity, which could have a direct impact on mCRPC. Collectively, these findings shed light on the underlying transcriptional mechanisms driving Abi response, demonstrating that TraRe is a promising tool for generating novel hypotheses based on identified transcriptional network disruptions.

The immune microenvironment of tumor plays a critical role in therapeutic responses to chemotherapy. Cancer tissues are composed of a complex network between antitumor and pro-tumor immune cells and molecules; therefore a comprehensive analysis of the tumor immune condition is imperative for better understanding of the roles of the immune microenvironment in anticancer treatment response. In this study, we performed T cell receptor (TCR) repertoire analysis of tumor infiltrating T cells (TILs) in cancer tissues of pre- and post-neoadjuvant chemotherapy (NAC) from 19 breast cancer patients; five cases showed CR (complete response), ten showed PR (partial response), and four showed SD/PD (stable disease/progressive disease) to the treatment. From the TCR sequencing results, we calculated the diversity index of the TCRβ chain and found that clonal expansion of TILs could be detected in patients who showed CR or PR to NAC. Noteworthy, the diversity of TCR was further reduced in the post-NAC tumors of CR patients. Our quantitative RT-PCR also showed that expression ratio of CD8/Foxp3 was significantly elevated in the post-NAC tumors of CR cases (p=0.0032), indicating that antitumor T cells were activated and enriched in these tumors. Collectively, our findings suggest that the clonal expansion of antitumor T cells may be a critical factor associated with response to chemotherapy and that their TCR sequences might be applicable for the development of TCR-engineered T cells treatment for individual breast cancer patients when their tumors relapse.

To determine the frequency of pathogenic inherited mutations in 157 select genes from patients with metastatic castrate-resistant prostate cancer (mCRPC).

The role of tumor protein D54 in breast cancer has not been studied and its function in breast cancer remains unclear. In our previous pharmacogenomic studies using lymphoblastoid cell line (LCL), this protein has been identified to affect metformin response. Although metformin has been widely studied as a prophylactic and chemotherapeutic drug, there is still a lack of biomarkers predicting the response to metformin in breast cancer. In this study, we revealed the novel function of TPD54 in breast cancer through understanding how TPD54 altered the cancer cell sensitivity to metformin.

Tamoxifen is one of the most commonly employed endocrine therapies for patients with estrogen receptor α (ERα)-positive breast cancer. Unfortunately the clinical benefit is limited due to intrinsic and acquired drug resistance. We previously reported a genome-wide association study that identified common SNPs near the CTSO gene and in ZNF423 associated with development of breast cancer during tamoxifen therapy in the NSABP P-1 and P-2 breast cancer prevention trials. Here, we have investigated their roles in ERα-positive breast cancer growth and tamoxifen response, focusing on the mechanism of CTSO. We performed in vitro studies including luciferase assays, cell proliferation, and mass spectrometry-based assays using ERα-positive breast cancer cells and a panel of genomic data-rich lymphoblastoid cell lines. We report that CTSO reduces the protein levels of BRCA1 and ZNF423 through cysteine proteinase-mediated degradation. We also have identified a series of transcription factors of BRCA1 that are regulated by CTSO at the protein level. Importantly, the variant CTSO SNP genotypes are associated with increased CTSO and decreased BRCA1 protein levels that confer resistance to tamoxifen. Characterization of the effect of both CTSO SNPs and ZNF423 SNPs on tamoxifen response revealed that cells with different combinations of CTSO and ZNF423 genotypes respond differently to Tamoxifen, PARP inhibitors or the combination of the two drugs due to SNP dependent differential regulation of BRCA1 levels. Therefore, these genotypes might be biomarkers for selection of individual drug to achieve the best efficacy.