Weinshilboum, Richard M.

BACKGROUND: Dilated cardiomyopathy (DCM) is a major cause of heart failure and carries a high mortality rate. Myocardial recovery in DCM-related heart failure patients is highly variable, with some patients having little or no response to standard drug therapy. A genome-wide association study may agnostically identify biomarkers and provide novel insight into the biology of myocardial recovery in DCM. METHODS: A genome-wide association study for change in left ventricular ejection fraction was performed in 686 White subjects with recent-onset DCM who received standard pharmacotherapy. Genome-wide association study signals were subsequently functionally validated and studied in relevant cellular models to understand molecular mechanisms that may have contributed to the change in left ventricular ejection fraction. RESULTS: The genome-wide association study identified a highly suggestive locus that mapped to the 5′-flanking region of the CDCP1 (CUB [complement C1r/C1s, Uegf, and Bmp1] domain containing protein 1) gene (rs6773435; P =7.12×10 −7 ). The variant allele was associated with improved cardiac function and decreased CDCP1 transcription. CDCP1 expression was significantly upregulated in human cardiac fibroblasts (HCFs) in response to the PDGF (platelet-derived growth factor) signaling, and knockdown of CDCP1 significantly repressed HCF proliferation and decreased AKT (protein kinase B) phosphorylation. Transcriptomic profiling after CDCP1 knockdown in HCFs supported the conclusion that CDCP1 regulates HCF proliferation and mitosis. In addition, CDCP1 knockdown in HCFs resulted in significantly decreased expression of soluble ST2 (suppression of tumorigenicity-2), a prognostic biomarker for heart failure and inductor of cardiac fibrosis. CONCLUSIONS: CDCP1 may play an important role in myocardial recovery in recent-onset DCM and mediates its effect primarily by attenuating cardiac fibrosis.

AbstractObjectiveElevated serum chromogranin A (CGA) is associated with intrinsic or treatment‐related neuroendocrine differentiation (NED) in men with metastatic castration‐resistant prostate cancer (mCRPC). Fluctuations in serum CGA during treatment of mCRPC have had conflicting results. We analyzed the impact of (i) rising serum CGA and (ii) baseline CGA/PSA ratio during treatment to identify associations with abiraterone acetate (AA) therapy.MethodsBetween June 2013 and August 2015, 92 men with mCRPC were enrolled in a prospective trial with uniform serum CGA processing performed before initiating abiraterone acetate/prednisone (AA/P) and serially after 12 weeks of AA/P treatments. Serum CGA was measured using a homogenous automated immunofluorescent assay. Patients receiving proton pump inhibitors or with abnormal renal function were excluded due to possible false elevations of serum CGA (n = 21 excluded), therefore 71 patients were analyzed. All patients underwent a composite response assessment at 12‐weeks. Kaplan–Meier estimates and Cox Regression models were used to calculate the association with time‐to‐treatment failure analyses and overall survival.ResultsAn increase in chromogranin was associated with a lower risk of treatment failure (hazard ratio [HR]: 0.52, p = 0.0181). The median CGA/PSA ratio was 7.8 (2.6–16.0) and an elevated pretreatment CGA/PSA ratio above the median was associated with a lower risk of treatment failure (HR: 0.54 p value = 0.0185). An increase in CGA was not found to be associated with OS (HR: 0.71, 95% CI: 0.42–1.21, p = 0.207). An elevated baseline CGA/PSA ratio was not associated with OS (HR: 0.62, 95% CI: 0.37–1.03, p = 0.062). An increase in PSA after 12 weeks of treatment was associated with an increased risk of treatment failure (HR: 4.14, CI: 2.21–7.73, p = < 0.0001) and worse OS (HR: 2.93, CI: 1.57–4.45, p = < 0.0001).ConclusionsWe show that an increasing chromogranin on AA/P and an elevated baseline CGA/PSA in patients with mCRPC were associated with a favorable response to AA/P with no changes in survival. There may be limited clinical utility in serum CGA testing to evaluate for lethal NED as AA/P did not induce lethal NED in this cohort. This highlights that not all patients with an increasing CGA have a worse OS.

201 Background: Previous studies have evaluated peripheral whole blood mRNA to identify prognostic biomarkers in metastatic prostate cancer, with no consistent genes of prognostic significance between studies and sample collection done at various points of therapy. We aim to identify the prognostic value of 20 previously identified genes in patients mCRPC from a prospective clinical trial. Methods: Between June 2013 and August 2015, whole blood was prospectively collected in PAXgene RNA tubes from 92 men with mCRPC (NCT01953640). Gene expression from samples with RIN > 4.0 were quantified using a customized NanoString CodeSet of 20 candidate genes (C1QA, MCM2, STOM, GABARAPL2, R10K3, HMBS, TFDP1, CDKN1A, ITGAL, TMCC2, TERF2IP, MMP9, SNCA, TIMP1, SLC4A1, SELENBP1, TMEM66, SEMA4D, PROS1, CD22). The raw data from the nCounter Analysis system were normalized to the positive controls and log2 transformed. Univariable Cox proportional hazards models were performed on each target gene sufficiently above negative control (5 times the maximum value of the negative control) to evaluate associations with overall survival (OS). For each gene we considered a linear, quadratic, or categorical model with quantile cutpoints. Categorical models with cutpoints at the 25th percentile for the genes CD22, GABARAPL2, R10K3, and STOM were used as there was a nonlinear relationship between gene expression and death hazard. The remaining genes were modeled with a linear model. Findings were considered statistically significant based on p-values adjusted using the Holm-Bonferroni correction. Results: 18/20 genes were identified and passed quality control (QC) measures, which excluded C1QA and MCM2. 72/92 men had an evaluable pretreatment specimen that passed QC. The median age was 72 years. The study cohort was evenly balanced in terms of Gleason score ( < = 7 v. > 7) and 60% had high volume disease. The median follow-up time was 5.73 years (IQR: 5.39, 6.49) and 58 patients had died. Of the 18 genes evaluated, STOM (Stomatin) [hazard ratio (HR) 3.00, adjusted p = 0.007] and GABARAPL2 (GABA Type A Receptor Associated Protein Like 2) [HR 2.62, adjusted p = 0.033] were significantly associated with overall survival in univariate analysis. In multivariate analysis, after adjusting for volume of metastatic disease and Gleason score, STOM [HR 2.89, adjusted p = 0.014] was significantly associated with overall survival. Conclusions: Elevated whole blood mRNA expression of STOM is adversely associated with OS in mCRPC. Further studies investigating elucidating the molecular significance of STOM in mCRPC are planned. Clinical trial information: NCT01953640 .

Abstract 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. Significance: The computational method TraRe built on Bayesian machine learning models for investigating transcriptional network structures shows that disruption of ELK3, MXD1, and MYB signaling cascades impacts abiraterone resistance in prostate cancer.

Abstract Androgen receptor (AR) is expressed in 80% to 90% of estrogen receptor α–positive (ER+) breast cancers. Accumulated evidence has shown that AR is a tumor suppressor and that its expression is associated with improved prognosis in ER+ breast cancer. However, both a selective AR agonist (RAD140) and an AR inhibitor (enzalutamide, ENZ) have shown a therapeutic effect on ER+ breast cancer, so the potential for clinical application of AR-targeting therapy for ER+ breast cancer is still in dispute. In this study, we evaluated the efficacy of ENZ and RAD140 in vivo and in vitro in AR+/ER+ breast cancer models, characterizing the relationship of AR and ER levels to response to AR-targeting drugs and investigating the alterations of global gene expression and chromatin binding of AR and ERα after ENZ treatment. In the AR-low setting, ENZ directly functioned as an ERα antagonist. Cell growth inhibition by ENZ in breast cancer with low AR expression was independent of AR and instead dependent on ER. In AR-high breast cancer models, AR repressed ERα signaling and ENZ promoted ERα signaling by antagonizing AR. In contrast, RAD140 activated AR signaling and suppressed AR-high tumor growth by deregulating ERα expression and blocking ERα function. Overall, analysis of the dynamic efficacies and outcomes of AR agonist, and antagonist in the presence of different AR and ERα levels reveals regulators of response and supports the clinical investigation of ENZ in selected ER+ tumors with a low AR/ER ratio and AR agonists in tumors with a high AR/ER ratio. Significance: The ratio of androgen receptor to estrogen receptor in breast cancer dictates the response to AR-targeted therapies, providing guidelines for developing AR-directed treatment strategies for patients with breast cancer.

Abstract Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were ‘unmasked’ only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

Abstract We identified resistance mechanisms to abiraterone acetate/prednisone (AA/P) in patients with metastatic castration-resistant prostate cancer (mCRPC) in the Prostate Cancer Medically Optimized Genome-Enhanced Therapy (PROMOTE) study. We analyzed whole-exome sequencing (WES) and RNA-sequencing data from 83 patients with metastatic biopsies before (V1) and after 12 weeks of AA/P treatment (V2). Resistance was determined by time to treatment change (TTTC). At V2, 18 and 11 of 58 patients had either short-term (median 3.6 months; range 1.4–4.5) or long-term (median 29 months; range 23.5–41.7) responses, respectively. Nonresponders had low expression of TGFBR3 and increased activation of the Wnt pathway, cell cycle, upregulation of AR variants, both pre- and posttreatment, with further deletion of AR inhibitor CDK11B posttreatment. Deletion of androgen processing genes, HSD17B11, CYP19A1 were observed in nonresponders posttreatment. Genes involved in cell cycle, DNA repair, Wnt-signaling, and Aurora kinase pathways were differentially expressed between the responder and non-responder at V2. Activation of Wnt signaling in nonresponder and deactivation of MYC or its target genes in responders was detected via SCN loss, somatic mutations, and transcriptomics. Upregulation of genes in the AURKA pathway are consistent with the activation of MYC regulated genes in nonresponders. Several genes in the AKT1 axis had increased mutation rate in nonresponders. We also found evidence of resistance via PDCD1 overexpression in responders. Implications: Finally, we identified candidates drugs to reverse AA/P resistance: topoisomerase inhibitors and drugs targeting the cell cycle via the MYC/AURKA/AURKB/TOP2A and/or PI3K_AKT_MTOR pathways.

AbstractAlternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3′ ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome‐wide the impact of glucocorticoids on APA in 30 human B‐lymphoblastoid cell lines. We found that glucocorticoids could regulate APA for a subset of genes, possibly by changing the expression of 142 RNA‐binding proteins, some with known APA‐regulating properties. Interestingly, genes with glucocorticoid‐mediated APA were enriched in viral translation‐related pathways, while genes with glucocorticoid‐mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid‐mediated APA might result in functional consequences distinct from gene expression. For example, glucocorticoids, a pharmacotherapy for severe COVID‐19, were found to change the APA but not the expression of LY6E, an important antiviral inhibitor in coronavirus diseases. Glucocorticoid‐mediated APA was also cell‐type‐specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype‐dependent glucocorticoid‐mediated APA (referred to as pharmacogenomic‐alterative polyadenylation quantitative trait loci), providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of observations regarding the impact of glucocorticoids on APA, raising the possibility that this mechanism might have implications for both disease pathophysiology and drug therapy.

AbstractBackground:Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype. Patients with TNBC are primarily treated with neoadjuvant chemotherapy (NAC). Response to NAC is prognostic, with reductions in low overall survival and disease-free survival rates in those patients who do not achieve pathological complete response (pCR). Based on this premise, we focused on patients with residual disease and hypothesized that paired analysis of primary and residual TNBC tumors following NAC could identify unique biomarkers associated with post-NAC recurrence.Methods and Results: We investigated 24 TNBC tumors from 12 patients with pre- and post-NAC data, including four patients who recurred shortly after surgery(<24 months) and eight who remained recurrence-free (> 48 months). These tumors were collected from a prospective NAC breast cancer study (BEAUTY) conducted at the Mayo Clinic. Differential expression analysis of pre-NAC biopsies showed minimal gene expression differences between early recurrent and non-recurrent TNBC tumors; however, post-NAC samples demonstrated a significant alteration in expression patterns in response to intervention. Topological-level differences associated with early recurrence were implicated in 251 gene sets, and an independent assessment of microarray gene expression data from the paired NAC I-SPY1 trial confirmed 56 genesets. Underlying these 56 gene sets, 113 genes were observed to be differentially expressed in the I-SPY1 and BEAUTY post-NAC studies. An independent (n=392) breast cancer data set with relapse-free survival (RFS) data was used to refine our gene list to a 17 gene signature. A three-fold cross-validation analysis of the gene signature with the combined BEAUTY and I-SPY1 data yielded an average AUC of 0.88 for six machine learning models. Due to the limited studies with pre-and post-NAC TNBC tumor data, further validation of the signature is needed.Conclusion:Analysis of multi-omics data from post-NAC TNBC chemoresistant tumors showed downregulation of mismatch repair, purin, and tubulin pathways. Additionally, weidentified a 17 gene signature in TNBC associated with post-NAC recurrence enriched with down-regulation of immune genes.