Myocardial Recovery in Recent Onset Dilated Cardiomyopathy: Role of CDCP1 and Cardiac Fibrosis.

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.

Androgen receptor-mediated pharmacogenomic expression quantitative trait loci: implications for breast cancer response to AR-targeting therapy.

BACKGROUND: Endocrine therapy is the most important treatment modality of breast cancer patients whose tumors express the estrogen receptor α (ERα). The androgen receptor (AR) is also expressed in the vast majority (80-90%) of ERα-positive tumors. AR-targeting drugs are not used in clinical practice, but have been evaluated in multiple trials and preclinical studies. METHODS: We performed a genome-wide study to identify hormone/drug-induced single nucleotide polymorphism (SNP) genotype - dependent gene-expression, known as PGx-eQTL, mediated by either an AR agonist (dihydrotestosterone) or a partial antagonist (enzalutamide), utilizing a previously well characterized lymphoblastic cell line panel. The association of the identified SNPs-gene pairs with breast cancer phenotypes were then examined using three genome-wide association (GWAS) studies that we have published and other studies from the GWAS catalog. RESULTS: We identified 13 DHT-mediated PGx-eQTL loci and 23 Enz-mediated PGx-eQTL loci that were associated with breast cancer outcomes post ER antagonist or aromatase inhibitors (AI) treatment, or with pharmacodynamic (PD) effects of AIs. An additional 30 loci were found to be associated with cancer risk and sex-hormone binding globulin levels. The top loci involved the genes IDH2 and TMEM9, the expression of which were suppressed by DHT in a PGx-eQTL SNP genotype-dependent manner. Both of these genes were overexpressed in breast cancer and were associated with a poorer prognosis. Therefore, suppression of these genes by AR agonists may benefit patients with minor allele genotypes for these SNPs. CONCLUSIONS: We identified AR-related PGx-eQTL SNP-gene pairs that were associated with risks, outcomes and PD effects of endocrine therapy that may provide potential biomarkers for individualized treatment of breast cancer.

IL17RB genetic variants are associated with acamprosate treatment response in patients with alcohol use disorder: A proteomics-informed genomics study.

Acamprosate is a Food and Drug Administration (FDA) approved medication for the treatment of alcohol use disorder (AUD). However, only a subset of patients achieves optimal treatment outcomes. Currently, no biological measures are utilized to predict response to acamprosate treatment. We applied our established pharmaco-omics informed genomics strategy to identify potential biomarkers associated with acamprosate treatment response. Specifically, our previous open-label acamprosate clinical trial recruited 442 patients with AUD who were treated with acamprosate for three months. We first performed proteomics using baseline plasma samples to identify potential biomarkers associated with acamprosate treatment outcomes. Next, we applied our established "proteomics-informed genome-wide association study (GWAS)" research strategy, and identified 12 proteins, including interleukin-17 receptor B (IL17RB), associated with acamprosate treatment response.​ A GWAS for IL17RB concentrations identified several genome-wide significant signals. Specifically, the top hit single nucleotide polymorphism (SNP) rs6801605 with a minor allele frequency of 38% in the European American population mapped 4 kilobase (Kb) upstream of IL17RB, and intron 1 of the choline dehydrogenase (CHDH) gene on chromosome 3 (p: 4.8E-20). The variant genotype (AA) for the SNP rs6801605 was associated with lower IL17RB protein expression. In addition, we identified a series of genetic variants in IL17RB that were associated with acamprosate treatment outcomes. Furthermore, the variantgenotypes for all of those IL17RB SNPs were protective for alcohol relapse. Finally, we demonstrated that the basal level of mRNA expression of IL17RB was inversely correlated with those of nuclear factor-κB (NF-κB) subunits, and a significantly higher expression of NF-κB subunits was observed in AUD patients who relapsed to alcohol use. In summary, this study illustrates that IL17RB genetic variants might contribute to acamprosate treatment outcomes. This series of studies represents an important step toward generating functional hypotheses that could be tested to gain insight into mechanisms underlying acamprosate treatment response phenotypes. (The ClinicalTrials.gov Identifier: NCT00662571).

Glucocorticoids unmask silent non-coding genetic risk variants for common diseases.

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.

Dynamic assessment of serum chromogranin A and treatment response with abiraterone acetate in metastatic castration-resistant prostate cancer.

OBJECTIVE: Elevated 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. METHODS: Between 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. RESULTS: An 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). CONCLUSIONS: We 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.

Cytochrome P450 Transcriptional Regulation by Testis-Specific Y-Encoded-Like Protein: Identification of Novel Upstream Transcription Factors.

Cytochrome P450s (CYPs) display significant inter-individual variation in expression, much of which remains unexplained by known CYP single-nucleotide polymorphisms (SNPs). Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators for several drug-metabolizing CYPs including CYP3A4 However, transcription factors (TFs) that might influence CYP expression through an effect on TSPYL expression are unknown. Therefore, we studied regulators of TSPYL expression in hepatic cell lines and their possible SNP-dependent variation. Specifically, we identified candidate TFs that might influence TSPYL expression using the ENCODE ChIPseq database. Subsequently, the expression of TSPYL1/2/4 as well as that of selected CYP targets for TSPYL regulation were assayed in hepatic cell lines before and after knockdown of TFs that might influence CYP expression through TSPYL-dependent mechanisms. Those results were confirmed by studies of TF binding to TSPYL1/2/4 gene promoter regions. In hepatic cell lines, knockdown of the REST and ZBTB7A TFs resulted in decreased TSPYL1 and TSPYL4 expression and increased CYP3A4 expression, changes reversed by TSPYL1/4 overexpression. Potential binding sites for REST and ZBTB7A on the promoters of TSPYL1 and TSPYL4 were confirmed by chromatin immunoprecipitation. Finally, common SNP variants in upstream binding sites on the TSPYL1/4 promoters were identified and luciferase reporter constructs confirmed SNP-dependent modulation of TSPYL1/4 gene transcription. In summary, we identified REST and ZBTB7A as regulators of the expression of TSPYL genes which themselves can contribute to regulation of CYP expression and-potentially-of drug metabolism. SNP-dependent modulation of TSPYL transcription may contribute to individual variation in both CYP expression and-downstream-drug response phenotypes. SIGNIFICANCE STATEMENT: Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators of cytochrome P450 (CYP) gene expression. Here, we report that variation in TSPYL expression as a result of the effects of genetically regulated TSPYL transcription factors is an additional factor that could result in downstream variation in CYP expression and potentially, as a result, variation in drug biotransformation.

Single cell transcriptomics reveals distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with opioid use disorder.

The opioid epidemic represents a national crisis. Oxycodone is one of the most prescribed opioid medications in the United States, whereas buprenorphine is currently the most prescribed medication for opioid use disorder (OUD) pharmacotherapy. Given the extensive use of prescription opioids and the global opioid epidemic, it is essential to understand how opioids modulate brain cell type function at the single-cell level. We performed single nucleus RNA-seq (snRNA-seq) using iPSC-derived forebrain organoids from three male OUD subjects in response to oxycodone, buprenorphine, or vehicle for seven days. We utilized the snRNA-seq data to identify differentially expressed genes following drug treatment using the Seurat integrative analysis pipeline. We utilized iPSC-derived forebrain organoids and single-cell sequencing technology as an unbiased tool to study cell-type-specific and drug-specific transcriptional responses. After quality control filtering, we analyzed 25787 cells and identified sixteen clusters using unsupervised clustering analysis. Our results reveal distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with OUD. Specifically, buprenorphine displayed a significant influence on transcription regulation in glial cells. However, oxycodone induced type I interferon signaling in many cell types, including neural cells in brain organoids. Finally, we demonstrate that oxycodone, but not buprenorphine activated STAT1 and induced the type I interferon signaling in patients with OUD. These data suggest that elevation of STAT1 expression associated with OUD might play a role in transcriptional regulation in response to oxycodone. In summary, our results provide novel mechanistic insight into drug action at single-cell resolution.

Knowledge-guided analysis of "omics" data using the KnowEnG cloud platform.

We present Knowledge Engine for Genomics (KnowEnG), a free-to-use computational system for analysis of genomics data sets, designed to accelerate biomedical discovery. It includes tools for popular bioinformatics tasks such as gene prioritization, sample clustering, gene set analysis, and expression signature analysis. The system specializes in "knowledge-guided" data mining and machine learning algorithms, in which user-provided data are analyzed in light of prior information about genes, aggregated from numerous knowledge bases and encoded in a massive "Knowledge Network." KnowEnG adheres to "FAIR" principles (findable, accessible, interoperable, and reuseable): its tools are easily portable to diverse computing environments, run on the cloud for scalable and cost-effective execution, and are interoperable with other computing platforms. The analysis tools are made available through multiple access modes, including a web portal with specialized visualization modules. We demonstrate the KnowEnG system's potential value in democratization of advanced tools for the modern genomics era through several case studies that use its tools to recreate and expand upon the published analysis of cancer data sets.

TCF7L2 lncRNA: a link between bipolar disorder and body mass index through glucocorticoid signaling.

Bipolar disorder (BD) and obesity are highly comorbid. We previously performed a genome-wide association study (GWAS) for BD risk accounting for the effect of body mass index (BMI), which identified a genome-wide significant single-nucleotide polymorphism (SNP) in the gene encoding the transcription factor 7 like 2 (TCF7L2). However, the molecular function of TCF7L2 in the central nervous system (CNS) and its possible role in the BD and BMI interaction remained unclear. In the present study, we demonstrated by studying human induced pluripotent stem cell (hiPSC)-derived astrocytes, cells that highly express TCF7L2 in the CNS, that the BD-BMI GWAS risk SNP is associated with glucocorticoid-dependent repression of the expression of a previously uncharacterized TCF7L2 transcript variant. That transcript is a long non-coding RNA (lncRNA-TCF7L2) that is highly expressed in the CNS but not in peripheral tissues such as the liver and pancreas that are involved in metabolism. In astrocytes, knockdown of the lncRNA-TCF7L2 resulted in decreased expression of the parent gene, TCF7L2, as well as alterations in the expression of a series of genes involved in insulin signaling and diabetes. We also studied the function of TCF7L2 in hiPSC-derived astrocytes by integrating RNA sequencing data after TCF7L2 knockdown with TCF7L2 chromatin-immunoprecipitation sequencing (ChIP-seq) data. Those studies showed that TCF7L2 directly regulated a series of BD risk genes. In summary, these results support the existence of a CNS-based mechanism underlying BD-BMI genetic risk, a mechanism based on a glucocorticoid-dependent expression quantitative trait locus that regulates the expression of a novel TCF7L2 non-coding transcript.

ERICH3: vesicular association and antidepressant treatment response.

Selective serotonin reuptake inhibitors (SSRIs) are standard of care for major depressive disorder (MDD) pharmacotherapy, but only approximately half of these patients remit on SSRI therapy. Our previous genome-wide association study identified a single-nucleotide polymorphism (SNP) signal across the glutamate-rich 3 (ERICH3) gene that was nearly genome-wide significantly associated with plasma serotonin (5-HT) concentrations, which were themselves associated with SSRI response for MDD patients enrolled in the Mayo Clinic PGRN-AMPS SSRI trial. In this study, we performed a meta-analysis which demonstrated that those SNPs were significantly associated with SSRI treatment outcomes in four independent MDD trials. However, the function of ERICH3 and molecular mechanism(s) by which it might be associated with plasma 5-HT concentrations and SSRI clinical response remained unclear. Therefore, we characterized the human ERICH3 gene functionally and identified ERICH3 mRNA transcripts and protein isoforms that are highly expressed in central nervous system cells. Coimmunoprecipitation identified a series of ERICH3 interacting proteins including clathrin heavy chain which are known to play a role in vesicular function. Immunofluorescence showed ERICH3 colocalization with 5-HT in vesicle-like structures, and ERICH3 knock-out dramatically decreased 5-HT staining in SK-N-SH cells as well as 5-HT concentrations in the culture media and cell lysates without changing the expression of 5-HT synthesizing or metabolizing enzymes. Finally, immunofluorescence also showed ERICH3 colocalization with dopamine in human iPSC-derived neurons. These results suggest that ERICH3 may play a significant role in vesicular function in serotonergic and other neuronal cell types, which might help explain its association with antidepressant treatment response.

Evidence for machine learning guided early prediction of acute outcomes in the treatment of depressed children and adolescents with antidepressants.

BACKGROUND: The treatment of depression in children and adolescents is a substantial public health challenge. This study examined artificial intelligence tools for the prediction of early outcomes in depressed children and adolescents treated with fluoxetine, duloxetine, or placebo. METHODS: The study samples included training datasets (N = 271) from patients with major depressive disorder (MDD) treated with fluoxetine and testing datasets from patients with MDD treated with duloxetine (N = 255) or placebo (N = 265). Treatment trajectories were generated using probabilistic graphical models (PGMs). Unsupervised machine learning identified specific depressive symptom profiles and related thresholds of improvement during acute treatment. RESULTS: Variation in six depressive symptoms (difficulty having fun, social withdrawal, excessive fatigue, irritability, low self-esteem, and depressed feelings) assessed with the Children's Depression Rating Scale-Revised at 4-6 weeks predicted treatment outcomes with fluoxetine at 10-12 weeks with an average accuracy of 73% in the training dataset. The same six symptoms predicted 10-12 week outcomes at 4-6 weeks in (a) duloxetine testing datasets with an average accuracy of 76% and (b) placebo-treated patients with accuracies of 67%. In placebo-treated patients, the accuracies of predicting response and remission were similar to antidepressants. Accuracies for predicting nonresponse to placebo treatment were significantly lower than antidepressants. CONCLUSIONS: PGMs provided clinically meaningful predictions in samples of depressed children and adolescents treated with fluoxetine or duloxetine. Future work should augment PGMs with biological data for refined predictions to guide the selection of pharmacological and psychotherapeutic treatment in children and adolescents with depression.

ZNF423 modulates the AMP-activated protein kinase pathway and metformin response in a single nucleotide polymorphisms, estrogen and selective estrogen receptor modulator dependent fashion.

OBJECTIVES: We previously discovered that the single nucleotide polymorphisms (SNP) rs9940645 in the ZNF423 gene regulate ZNF423 expression and serve as a potential biomarker for response to selective estrogen receptor modulators (SERMs). Here we explored pathways involved in ZNF423-mediated SERMs response and drugs that potentially sensitize SERMs. METHODS: RNA sequencing and label-free quantitative proteomics were performed to identify genes and pathways that are regulated by ZNF423 and the ZNF423 SNP. Both cultured cells and mouse xenograft models with different ZNF423 SNP genotypes were used to study the cellular responses to metformin. RESULTS: We identified ribosome and AMP-activated protein kinase (AMPK) signaling as potential pathways regulated by ZNF423 or ZNF423 rs9940645 SNP. Moreover, using clustered regularly interspaced short palindromic repeats/Cas9-engineered ZR75-1 breast cancer cells with different ZNF423 SNP genotypes, striking differences in cellular responses to metformin, either alone or in the combination of tamoxifen, were observed in both cell culture and the mouse xenograft model. CONCLUSIONS: We found that AMPK signaling is modulated by the ZNF423 rs9940645 SNP in estrogen and SERM-dependent fashion. The ZNF423 rs9940645 SNP affects metformin response in breast cancer and could be a potential biomarker for tailoring the metformin treatment.

Single-nucleotide polymorphism biomarkers of adjuvant anastrozole-induced estrogen suppression in early breast cancer.

OBJECTIVES: Based on our previous findings that postmenopausal women with estrone (E1) and estradiol (E2) concentrations at or above 1.3 pg/ml and 0.5 pg/ml, respectively, after 6 months of adjuvant anastrozole therapy had a three-fold risk of recurrence, we aimed to identify a single-nucleotide polymorphism (SNP)-based model that would predict elevated E1 and E2 and then validate it in an independent dataset. PATIENTS AND METHODS: The test set consisted of 322 women from the M3 study and the validation set consisted of 152 patients from MA.27. All patients were treated with adjuvant anastrozole, had on-anastrozole E1 and E2 concentrations and genome-wide genotyping. RESULTS: SNPs were identified from the M3 genome-wide association study. The best model to predict the E1-E2 phenotype with high balanced accuracy was a support vector machine model using clinical factors plus 46 SNPs. We did not have an independent cohort that is similar to the M3 study with clinical, E1-E2 phenotypes and genotype data to test our model. Hence, we chose a nested matched case-control cohort (MA.27 study) for testing. Our E1-E2 model was not validated but we found the MA.27 validation cohort was both clinically and genomically different. CONCLUSIONS: We identified a SNP-based model that had excellent performance characteristics for predicting the phenotype of elevated E1 and E2 in women treated with anastrozole. This model was not validated in an independent dataset but that dataset was clinically and genomically substantially different. The model will need validation in a prospective study.

SLCO1B1: Application and Limitations of Deep Mutational Scanning for Genomic Missense Variant Function.

SLCO1B1 (solute carrier organic anion transporter family member 1B1) is an important transmembrane hepatic uptake transporter. Genetic variants in the SLCO1B1 gene have been associated with altered protein folding, resulting in protein degradation and decreased transporter activity. Next-generation sequencing (NGS) of pharmacogenes is being applied increasingly to associate variation in drug response with genetic sequence variants. However, it is difficult to link variants of unknown significance with functional phenotypes using "one-at-a-time" functional systems. Deep mutational scanning (DMS) using a "landing pad cell-based system" is a high-throughput technique designed to analyze hundreds of gene open reading frame (ORF) missense variants in a parallel and scalable fashion. We have applied DMS to analyze 137 missense variants in the SLCO1B1 ORF obtained from the Exome Aggregation Consortium project. ORFs containing these variants were fused to green fluorescent protein and were integrated into "landing pad" cells. Florescence-activated cell sorting was performed to separate the cells into four groups based on fluorescence readout indicating protein expression at the single cell level. NGS was then performed and SLCO1B1 variant frequencies were used to determine protein abundance. We found that six variants not previously characterized functionally displayed less than 25% and another 12 displayed approximately 50% of wild-type protein expression. These results were then functionally validated by transporter studies. Severely damaging variants identified by DMS may have clinical relevance for SLCO1B1-dependent drug transport, but we need to exercise caution since the relatively small number of severely damaging variants identified raise questions with regard to the application of DMS to intrinsic membrane proteins such as organic anion transporter protein 1B1. SIGNIFICANCE STATEMENT: The functional implications of a large numbers of open reading frame (ORF) "variants of unknown significance" (VUS) in transporter genes have not been characterized. This study applied deep mutational scanning to determine the functional effects of VUS that have been observed in the ORF of SLCO1B1(s olute carrier organic anion transporter family member 1B1). Several severely damaging variants were identified, studied, and validated. These observations have implications for both the application of deep mutational scanning to intrinsic membrane proteins and for the clinical effect of drugs and endogenous compounds transported by SLCO1B1.

Therapeutic potential of triterpenoid saponin anemoside B4 from Pulsatilla chinensis.

Pulsatilla Decoction (Bai-Tou-Weng-Tang) has been used medically in China for thousands of years for the treatment of diseases caused by bacteria. In recent decades, Pulsatilla Decoction is becoming a well-known formula prescription used for the treatment of ulcerative colitis in traditional Chinese medicine. Pulsatilla chinensis is the chief herbal source of Pulsatilla Decoction, and it is rich in triterpenoid saponins, such as anemoside B4, anemoside A3, and 23-hydroxybetulinic acid. Anemoside B4 is the most abundant of that group and has been used as a quality control marker for Pulsatilla chinensis. As the major active component of Pulsatilla chinensis, anemoside B4 has also received attention as a pure compound for its therapeutic potential. In this review, we systematically analyze the findings on triterpenoid saponins, especially anemoside B4, anemoside A3 and 23-hydroxybetulinic acid, included in Pulsatilla chinensis and Pulsatilla Decoction. We discuss the pharmacokinetics and tissue distribution of these triterpenoid saponins as well as their biological activities. We also summarize the pharmacological effects of anemoside B4 and its two possible metabolites, anemoside A3 and 23-hydroxybetulinic acid, as pure compounds. In summary, this review sketches a profile of the state of existing knowledge with regard to the pharmacological effects of anemoside B4, especially its anti-inflammatory and immunomodulatory effects. These findings point to the possibility that anemoside B4 has potential to be studied further as a natural compound-originated immunomodulatory agent for the treatment of inflammatory diseases such as ulcerative colitis and thus, may represent one of the most important active components of Pulsatilla Decoction responsible for its anti-ulcerative colitis efficacy.

The lncRNA MIR2052HG regulates ERα levels and aromatase inhibitor resistance through LMTK3 by recruiting EGR1.

BACKGROUND: Our previous genome-wide association study using the MA.27 aromatase inhibitors adjuvant trial identified SNPs in the long noncoding RNA MIR2052HG associated with breast cancer-free interval. MIR2052HG maintained ERα both by promoting AKT/FOXO3-mediated ESR1 transcription and by limiting ubiquitin-mediated ERα degradation. Our goal was to further elucidate MIR2052HG's mechanism of action. METHODS: RNA-binding protein immunoprecipitation assays were performed to demonstrate that the transcription factor, early growth response protein 1 (EGR1), worked together with MIR2052HG to regulate that lemur tyrosine kinase-3 (LMTK3) transcription in MCF7/AC1 and CAMA-1 cells. The location of EGR1 on the LMTK3 gene locus was mapped using chromatin immunoprecipitation assays. The co-localization of MIR2052HG RNA and the LMTK3 gene locus was determined using RNA-DNA dual fluorescent in situ hybridization. Single-nucleotide polymorphisms (SNP) effects were evaluated using a panel of human lymphoblastoid cell lines. RESULTS: MIR2052HG depletion in breast cancer cells results in a decrease in LMTK3 expression and cell growth. Mechanistically, MIR2052HG interacts with EGR1 and facilitates its recruitment to the LMTK3 promoter. LMTK3 sustains ERα levels by reducing protein kinase C (PKC) activity, resulting in increased ESR1 transcription mediated through AKT/FOXO3 and reduced ERα degradation mediated by the PKC/MEK/ERK/RSK1 pathway. MIR2052HG regulated LMTK3 in a SNP- and aromatase inhibitor-dependent fashion: the variant SNP increased EGR1 binding to LMTK3 promoter in response to androstenedione, relative to wild-type genotype, a pattern that can be reversed by aromatase inhibitor treatment. Finally, LMTK3 overexpression abolished the effect of MIR2052HG on PKC activity and ERα levels. CONCLUSIONS: Our findings support a model in which the MIR2052HG regulates LMTK3 via EGR1, and LMTK3 regulates ERα stability via the PKC/MEK/ERK/RSK1 axis. These results reveal a direct role of MIR2052HG in LMTK3 regulation and raise the possibilities of targeting MIR2052HG or LMTK3 in ERα-positive breast cancer.

Deep sequencing across germline genome-wide association study signals relating to breast cancer events in women receiving aromatase inhibitors for adjuvant therapy of early breast cancer.

OBJECTIVE: To identify additional genetic variants beyond those observed in a previous genome-wide association study (GWAS) in women treated on the MA.27 clinical trial in which women were randomized to 5 years of adjuvant therapy with anastrozole or exemestane. PATIENTS AND METHODS: We performed a matched case-control study in 234 women who had a recurrence of breast cancer (cases) and 649 women who had not (controls). The analysis was restricted to White women with an estrogen receptor-positive breast cancer. Multiplex PCR-based targeted deep sequencing was performed of the MIR2052HG region on chromosome 8 between positions 75.4 and 75.7, a span of 300 kb, in an attempt to identify additional functional single nucleotide polymorphisms (SNPs). RESULTS: A total of 4677 unique variants were identified that had not been identified in the previous GWAS. Clinical Annotation of Variants analysis revealed 10 variants, including eight SNPs and two insertion-deletion mutations with moderate or high impact. However, none of the common and variant regions was significant after adjustment for the most significant SNP (rs13260300) identified in our previous GWAS. We performed haplotype analysis that revealed two regions in which the haplotypes lost significance when adjusted for this prior GWAS SNP and one region with two significant haplotypes (P = 0.046 and 0.031) after adjusting for the GWAS SNP. CONCLUSION: We were unable to identify common or rare variant regions that added value to the findings from our previous GWAS. We did find two haplotypes that were significant after adjusting for our top GWAS SNP but these were considered to be of marginal value.

4-Hydroxytamoxifen enhances sensitivity of estrogen receptor α-positive breast cancer to docetaxel in an estrogen and ZNF423 SNP-dependent fashion.

PURPOSE: In early stage, ERα-positive breast cancer, concurrent use of endocrine therapy and chemotherapy has not been shown to be superior to sequential use. We hypothesized that genetic biomarkers can aid in selecting patients who would benefit from chemo-endocrine therapy. Our previous studies revealed that ZNF423 is a transcription factor for BRCA1 and an intronic single nucleotide polymorphism (SNP) in ZNF423, rs9940645, determines tamoxifen response. Here, we identified mitosis-related genes that are regulated by ZNF423 which led us to investigate taxane response in a rs9940645 SNP- and tamoxifen-dependent fashion. METHODS: The Cancer Genome Atlas (TCGA) breast cancer dataset was used to identify genes correlated with ZNF423. Quantitative reverse transcription PCR, chromatin immunoprecipitation, and luciferase reporter assays were used to validate the gene regulation. We used CRISPR/Cas9 to engineer paired ZR-75-1 cells which differ only in ZNF423 rs9940645 SNP genotype to test SNP-dependent phenotypes including cell cycle and cell viability. We validated our findings in an additional two breast cancer cell lines, Hs578T-ERα and HCC1500. RESULTS: Mitosis-related genes VRK1 and PBK, which encode histone H3 kinases, were experimentally validated to be regulated by ZNF423. ZNF423 knockdown decreased VRK1 and PBK expression and activity. Additionally, ZNF423 knockdown enhanced docetaxel-induced G2/M arrest and cytotoxicity through VRK1 or PBK regulation. Lastly, cells carrying the rs9940645 variant genotype had increased G2/M arrest and decreased cell viability when treated with docetaxel in combination with estradiol and 4-OH-TAM. CONCLUSIONS: We identified ZNF423 regulated genes involved in the G2/M phase of the cell cycle. 4-OH-TAM sensitized ERα-positive breast cancer cells to docetaxel in a ZNF423 SNP-dependent manner. Our findings suggest that patients with rs9940645 variant genotype may benefit from concurrent tamoxifen and docetaxel. This would impact a substantial proportion of patients because this SNP has a minor allele frequency of 0.47.

Single Nucleotide Polymorphisms at a Distance from Aryl Hydrocarbon Receptor (AHR) Binding Sites Influence AHR Ligand-Dependent Gene Expression.

Greater than 90% of significant genome-wide association study (GWAS) single-nucleotide polymorphisms (SNPs) are in noncoding regions of the genome, but only 25.6% are known expression quantitative trait loci (eQTLs). Therefore, the function of many significant GWAS SNPs remains unclear. We have identified a novel type of eQTL for which SNPs distant from ligand-activated transcription factor (TF) binding sites can alter target gene expression in a SNP genotype-by-ligand-dependent fashion that we refer to as pharmacogenomic eQTLs (PGx-eQTLs)-loci that may have important pharmacotherapeutic implications. In the present study, we integrated chromatin immunoprecipitation-seq with RNA-seq and SNP genotype data for a panel of lymphoblastoid cell lines to identify 10 novel cis PGx-eQTLs dependent on the ligand-activated TF aryl hydrocarbon receptor (AHR)-a critical environmental sensor for xenobiotic (drug) and immune response. Those 10 cis PGx-eQTLs were eQTLs only after AHR ligand treatment, even though the SNPs did not create/destroy an AHR response element-the DNA sequence motif recognized and bound by AHR. Additional functional studies in multiple cell lines demonstrated that some cis PGx-eQTLs are functional in multiple cell types, whereas others displayed SNP-by-ligand-dependent effects in just one cell type. Furthermore, four of those cis PGx-eQTLs had previously been associated with clinical phenotypes, indicating that those loci might have the potential to inform clinical decisions. Therefore, SNPs across the genome that are distant from TF binding sites for ligand-activated TFs might function as PGx-eQTLs and, as a result, might have important clinical implications for interindividual variation in drug response. SIGNIFICANCE STATEMENT: More than 90% of single-nucleotide polymorphisms (SNPs) that are associated with clinical phenotypes are located in noncoding regions of the genome. However, the mechanisms of action of many of those SNPs have not been elucidated, and drugs may unmask functional expression quantitative trail loci (eQTLs). In the current study, we used drugs that bind to the ligand-activated transcription factor aryl hydrocarbon receptor (AHR) and identified SNPs that were associated with interindividual variation in gene expression following drug exposure-termed pharmacogenomic (PGx)-eQTLs. Possibly of greater significance, those PGx-eQTL SNPs were outside of AHR binding sites, indicating that they do not interrupt AHR DNA recognition. PGx-eQTLs such as those described in this work may have crucial implications for interindividual variation in drug.

Pharmacogenomic Next-Generation DNA Sequencing: Lessons from the Identification and Functional Characterization of Variants of Unknown Significance in CYP2C9 and CYP2C19.

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P < 0.05) and CYP2C19 (P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.