Molecular Profiling of Aromatase Inhibitor–Treated Postmenopausal Breast Tumors Identifies Immune-Related Correlates of Resist

Total Page:16

File Type:pdf, Size:1020Kb

Molecular Profiling of Aromatase Inhibitor–Treated Postmenopausal Breast Tumors Identifies Immune-Related Correlates of Resist Published OnlineFirst March 14, 2013; DOI: 10.1158/1078-0432.CCR-12-1000 Clinical Cancer Predictive Biomarkers and Personalized Medicine Research Molecular Profiling of Aromatase Inhibitor–Treated Postmenopausal Breast Tumors Identifies Immune-Related Correlates of Resistance Anita K. Dunbier1,2,5, Zara Ghazoui1,2, Helen Anderson1,2, Janine Salter1, Ashutosh Nerurkar1, Peter Osin1, Roger A'hern3, William R. Miller4, Ian E. Smith1, and Mitch Dowsett1,2 Abstract Purpose: Estrogen withdrawal by treatment with aromatase inhibitors is the most effective form of endocrine therapy for postmenopausal estrogen receptor–positive (ERþ) breast cancer. However, response to therapy varies markedly and understanding of the precise molecular effects of aromatase inhibitors and causes of resistance is limited. We aimed to identify in clinical breast cancer those genes and pathways most associated with resistance to aromatase inhibitors by examining the global transcriptional effects of AI treatment. Experimental Design: Baseline and 2-week posttreatment biopsies were obtained from 112 postmen- opausal women with ERþ breast cancer receiving neoadjuvant anastrozole. Gene expression data were obtained from 81 baseline and 2-week paired samples. Pathway analysis identified (i) the most prevalent changes in expression and (ii) the pretreatment genes/pathways most related to poor antiproliferative response. Results: A total of 1,327 genes were differentially expressed after 2-week treatment (false discovery rate < 0.01). Proliferation-associated genes and classical estrogen-dependent genes were strongly downregulated whereas collagens and chemokines were upregulated. Pretreatment expression of an inflammatory signature correlated with antiproliferative response to anastrozole and this observation was validated in an inde- pendent study. Higher expression of immune-related genes such as SLAMF8 and TNF as well as lymphocytic infiltration were associated with poorer response (P < 0.001) and validated in an independent cohort. Conclusions: The molecular response to aromatase inhibitor treatment varies greatly between patients consistent with the variable clinical benefit from aromatase inhibitor treatment. Higher baseline expression of an inflammatory signature is associated with poor antiproliferative response and should be assessed further as a novel biomarker and potential target for aromatase inhibitor-treated patients. Clin Cancer Res; 19(10); 2775–86. Ó2013 AACR. Introduction inhibitors is the most effective form of endocrine therapy, Approximately 80% of human breast carcinomas present but response to them varies markedly (1–3). Mechanisms of as estrogen receptor a–positive (ERþ). In postmenopausal resistance appear to be multiple but are poorly character- In vitro women, estrogen withdrawal by treatment with aromatase ized. studies of acquired resistance to estrogen deprivation have identified several putative mechanisms which largely involve growth factor–related signal trans- Authors' Affiliations: 1Royal Marsden Hospital; 2Breakthrough Breast duction pathways (4–7) but there is limited clinical evi- 3 Cancer Research Centre, Institute of Cancer Research; Cancer Research dence to support these. It is also notable that these models UK Clinical Trials and Statistics Unit, Section of Clinical Trials, Institute of Cancer Research, London; 4Breast Research Group, University of Edin- do not involve an assessment of the contribution of human burgh, Edinburgh, United Kingdom; and 5Department of Biochemistry, stromal elements. University of Otago, Dunedin, New Zealand The presurgical/neoadjuvant setting provides an excep- Note: Supplementary data for this article are available at Clinical Cancer tionally valuable scenario for linking biology to clinical Research Online (http://clincancerres.aacrjournals.org/). response and to study mechanisms of resistance; within this A.K. Dunbier, Z. Ghazoui, and H. Anderson contributed equally to this work. setting, the proliferation marker Ki67 is a validated phar- Corresponding Author: Anita K. Dunbier, Department of Biochemistry, macodynamic indicator of response to endocrine therapy University of Otago, PO Box 56, Dunedin 9016, New Zealand. Phone: (8–10). Treatment-dependent changes in sequential mea- 643479-9258; Fax: 643479-7738; E-mail: [email protected] surements of Ki67 in neoadjuvant trials of endocrine ther- doi: 10.1158/1078-0432.CCR-12-1000 apy in postmenopausal women (IMPACT, p024, and Ó2013 American Association for Cancer Research. Z1031; refs. 9, 11, 12) have all revealed differences or lack www.aacrjournals.org 2775 Downloaded from clincancerres.aacrjournals.org on September 28, 2021. © 2013 American Association for Cancer Research. Published OnlineFirst March 14, 2013; DOI: 10.1158/1078-0432.CCR-12-1000 Dunbier et al. Materials and Methods Translational Relevance Patient samples Most postmenopausal women with estrogen receptor Fourteen-gauge core-cut tumor biopsies were obtained (ER)–positive breast cancer receive an aromatase inhib- from 112 postmenopausal women with stage I to IIIB ERþ itor at some stage during their treatment. However, early breast cancer before and after 2-week anastrozole responsiveness to aromatase inhibitors varies greatly. monotherapy in a neoadjuvant trial (17, 18). Patient demo- Previous studies of resistance mechanisms have been graphics are shown in Supplementary Table S1. Tissue was conducted in cell lines and have highlighted the impor- stored in RNAlater at À20C. Two 4-mm sections from the tance of signal transduction pathways, but these studies core were stained with hematoxylin and eosin (H&E) and omit the stromal influences that are increasingly recog- examined by a pathologist (A. Nerurkar and P. Osin) to nized to have a major influence on tumor biology. In confirm the presence of cancerous tissue, assess histopa- contrast, we report genome-wide expression profiling of thology, and the presence or absence of lymphocytic infil- 81 ER-positive breast carcinomas including stromal tis- tration. Tumors were deemed to be positive for lymphocytic sues before and during treatment with an aromatase infiltration if intraepithelial mononuclear cells could be inhibitor in the neoadjuvant setting. We identify an seen within tumor cell nests or in direct contact with tumor inflammatory signature as the strongest correlate of poor cells. Biopsies in which lymphocytic infiltrate could be seen antiproliferative response and confirm this finding in an without direct contact to tumor cells were not considered to independent set of tumors as well as using gross lym- have lymphocytic infiltration. Samples were assessed blind- phocytic infiltration as an alternative measure of inflam- ly by both pathologists and the concordance rate was 81%. matory activity. This work provides a potential new Discordant samples were then reassessed by A. Nerurkar avenue for drug development and for identifying who was blinded to the result of the initial assessment and patients with a reduced likelihood of response to aro- the more frequent of the 3 assessments was used in the matase inhibition. analysis. Total RNA was extracted using RNeasy (Qiagen). RNA quality was checked using an Agilent Bioanalyser: samples with RNA integrity values of less than 5 were excluded from further analysis. ER (H-score) and Ki67 (% cells positive) values by centralized immunohistochem- of differences that were parallel to recurrence-free survival istry were already available (17). (RFS) differences in the equivalent adjuvant trials (ATAC, An additional 71 stored sections from paraffin-embed- BIG 1-98, and MA27; refs. 1, 13, 14), respectively. In ded, formalin-fixed core-cut biopsies were obtained from a addition, Ki67 of patients after 2 weeks’ treatment predicts subset of patients from the IMPACT trial who received RFS more closely than pretreatment values (8). neoadjuvant anastrozole or tamoxifen (15). Sections were Nearly all patients on an aromatase inhibitor show H&E-stained, and the presence or absence of detectable reduction in Ki67 expression, suggesting that some benefit lymphocytic infiltration was assessed by a pathologist (A. is derived, although this may be modest (9). This contin- Nerurkar). Ki67 data were obtained from previous analysis uous marker of response is well-suited to assessment of of these patients (8, 9). mechanisms of resistance that are also likely to have a nonbinary effect. This may explain why changes in Ki67 have proven to be better predictors of benefit from endo- Gene expression analysis and data preprocessing crine therapy than clinical response during neoadjuvant RNA amplification, labeling, and hybridization on endocrine therapy (15, 16). HumanWG-6 v2 Expression BeadChips (Illumina) were Gene expression profiling of tumor biopsies before and conducted according to the manufacturer’s instructions at during treatment has the potential to enable the identifica- a single Illumina BeadStation facility. Tumor RNA of suf- tion of the most important genes/pathways involved in the ficient quality and quantity was available to generate expres- response to estrogen deprivation therapy and the pretreat- sion data from 104 pretreatment biopsies and 85 two-week ment determinants of response and resistance. Availability biopsies (Supplementary Fig. S1). Data were extracted using of comprehensive expression datasets, as provided here, will BeadStudio software and normalized with variance-stabi- allow the evaluation of candidate genes from experimental lizing transformation (VST) and Robust Spline Normaliza- research
Recommended publications
  • Murine Osteoclasts Secrete Serine Protease Htra1 Capable of Degrading Osteoprotegerin in the Bone Microenvironment
    ARTICLE https://doi.org/10.1038/s42003-019-0334-5 OPEN Murine osteoclasts secrete serine protease HtrA1 capable of degrading osteoprotegerin in the bone microenvironment Nagahiro Ochiai1,2, Yutaka Nakachi1, Tomotaka Yokoo1, Takahiro Ichihara2, Tore Eriksson 3, Yuki Yonemoto2, 1234567890():,; Takehiko Kato2, Hitoshi Ogata2, Natsuko Fujimoto2, Yasuhiro Kobayashi4, Nobuyuki Udagawa5, Shinsuke Kaku2, Tomokazu Ueki2, Yasushi Okazaki 1,6, Naoyuki Takahashi4 & Tatsuo Suda1 Osteoclasts are multinucleated cells responsible for bone resorption. The differentiation of osteoclasts from bone marrow macrophages (BMMs) is induced by receptor activator of NF- κB ligand (RANKL). Osteoprotegerin (OPG), a decoy receptor of RANKL, inhibits osteo- clastogenesis by blocking RANKL signaling. Here we investigated the degradation of OPG in vitro. Osteoclasts, but not BMMs, secreted OPG-degrading enzymes. Using mass spec- trometry and RNA-sequencing analysis, we identified high-temperature requirement A serine peptidase 1 (HtrA1) as an OPG-degrading enzyme. HtrA1 did not degrade OPG pre-reduced by dithiothreitol, suggesting that HtrA1 recognizes the three-dimensional structure of OPG. HtrA1 initially cleaved the amide bond between leucine 90 and glutamine 91 of OPG, then degraded OPG into small fragments. Inhibitory activity of OPG on RANKL-induced osteo- clastogenesis was suppressed by adding HtrA1 in RAW 264.7 cell cultures. These results suggest that osteoclasts potentially prepare a microenvironment suitable for osteoclasto- genesis. HtrA1 may be a novel drug target for osteoporosis. 1 Research Center for Genomic Medicine, Saitama Medical University, Saitama 350-1298, Japan. 2 Pharmacology Laboratories, Taisho Pharmaceutical Co., Ltd, Saitama 331-9530, Japan. 3 Chemistry Laboratories, Taisho Pharmaceutical Co., Ltd, Saitama 331-9530, Japan. 4 Institutes for Oral Science, Matsumoto Dental University, Nagano 399-0781, Japan.
    [Show full text]
  • An Animal Model with a Cardiomyocyte-Specific Deletion of Estrogen Receptor Alpha: Functional, Metabolic, and Differential Netwo
    Washington University School of Medicine Digital Commons@Becker Open Access Publications 2014 An animal model with a cardiomyocyte-specific deletion of estrogen receptor alpha: Functional, metabolic, and differential network analysis Sriram Devanathan Washington University School of Medicine in St. Louis Timothy Whitehead Washington University School of Medicine in St. Louis George G. Schweitzer Washington University School of Medicine in St. Louis Nicole Fettig Washington University School of Medicine in St. Louis Attila Kovacs Washington University School of Medicine in St. Louis See next page for additional authors Follow this and additional works at: https://digitalcommons.wustl.edu/open_access_pubs Recommended Citation Devanathan, Sriram; Whitehead, Timothy; Schweitzer, George G.; Fettig, Nicole; Kovacs, Attila; Korach, Kenneth S.; Finck, Brian N.; and Shoghi, Kooresh I., ,"An animal model with a cardiomyocyte-specific deletion of estrogen receptor alpha: Functional, metabolic, and differential network analysis." PLoS One.9,7. e101900. (2014). https://digitalcommons.wustl.edu/open_access_pubs/3326 This Open Access Publication is brought to you for free and open access by Digital Commons@Becker. It has been accepted for inclusion in Open Access Publications by an authorized administrator of Digital Commons@Becker. For more information, please contact [email protected]. Authors Sriram Devanathan, Timothy Whitehead, George G. Schweitzer, Nicole Fettig, Attila Kovacs, Kenneth S. Korach, Brian N. Finck, and Kooresh I. Shoghi This open access publication is available at Digital Commons@Becker: https://digitalcommons.wustl.edu/open_access_pubs/3326 An Animal Model with a Cardiomyocyte-Specific Deletion of Estrogen Receptor Alpha: Functional, Metabolic, and Differential Network Analysis Sriram Devanathan1, Timothy Whitehead1, George G. Schweitzer2, Nicole Fettig1, Attila Kovacs3, Kenneth S.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Properties of Human Genes Guided by Their Enrichment in Rare and Common Variants
    Properties of human genes guided by their enrichment in rare and common variants Authors: Eman Alhuzimi, Luis G. Leal, Michael J.E. Sternberg, Alessia David Affiliation: Structural Bioinformatics Group, Department of Life Sciences, Imperial College London, London, SW7 2AZ, UK SUPPLEMENTARY MATERIAL Construction of the dataset Genetic variants occurring in protein coding genes were extracted from ExAC (version 0.3, Release: 13-Jan-2015), UniProt (humsavar.txt, release: 04-Feb-2015) and ClinVar (release:7-Jan-2015). Variants were classified as ‘disease-causing’ if a disease association was reported in humsavar.txt or ClinVar. For variants reported in ClinVar, we defined the variant as disease-causing only if it was annotated as “pathogenic”. In order to avoid a potential bias, variants annotated as “likely pathogenic” were not included in the analysis. Variants were classified as ‘neutral’ when no association with disease was present (variants reported as “polymorphisms” in humsavar.txt and variants from ExAC, not reported as disease-causing in other databases). Non-disease variants were divided according to their global minor allele frequencies (MAF) into: ‘rare variants’ (MAF < 0.01) and ‘common variants’ (MAF ≥ 0.01). Global MAF data were extracted from Ensembl using the BioMart data-mining tool. We used the global MAF calculated in the ExAC project. For variants not reported in ExAC database we used the global MAF reported in dbSNP (which is calculated from the 1000Genomes project), when available. Variants with no MAF information or reported as “unclassified” in humsavar.txt, were not included in the analysis. When the gene enrichment analysis (described below) was performed, one gene overlapped between the disease- and rare- EVsets and three genes between the disease- and common- EVsets.
    [Show full text]
  • 4-6 Weeks Old Female C57BL/6 Mice Obtained from Jackson Labs Were Used for Cell Isolation
    Methods Mice: 4-6 weeks old female C57BL/6 mice obtained from Jackson labs were used for cell isolation. Female Foxp3-IRES-GFP reporter mice (1), backcrossed to B6/C57 background for 10 generations, were used for the isolation of naïve CD4 and naïve CD8 cells for the RNAseq experiments. The mice were housed in pathogen-free animal facility in the La Jolla Institute for Allergy and Immunology and were used according to protocols approved by the Institutional Animal Care and use Committee. Preparation of cells: Subsets of thymocytes were isolated by cell sorting as previously described (2), after cell surface staining using CD4 (GK1.5), CD8 (53-6.7), CD3ε (145- 2C11), CD24 (M1/69) (all from Biolegend). DP cells: CD4+CD8 int/hi; CD4 SP cells: CD4CD3 hi, CD24 int/lo; CD8 SP cells: CD8 int/hi CD4 CD3 hi, CD24 int/lo (Fig S2). Peripheral subsets were isolated after pooling spleen and lymph nodes. T cells were enriched by negative isolation using Dynabeads (Dynabeads untouched mouse T cells, 11413D, Invitrogen). After surface staining for CD4 (GK1.5), CD8 (53-6.7), CD62L (MEL-14), CD25 (PC61) and CD44 (IM7), naïve CD4+CD62L hiCD25-CD44lo and naïve CD8+CD62L hiCD25-CD44lo were obtained by sorting (BD FACS Aria). Additionally, for the RNAseq experiments, CD4 and CD8 naïve cells were isolated by sorting T cells from the Foxp3- IRES-GFP mice: CD4+CD62LhiCD25–CD44lo GFP(FOXP3)– and CD8+CD62LhiCD25– CD44lo GFP(FOXP3)– (antibodies were from Biolegend). In some cases, naïve CD4 cells were cultured in vitro under Th1 or Th2 polarizing conditions (3, 4).
    [Show full text]
  • Download Download
    Supplementary Figure S1. Results of flow cytometry analysis, performed to estimate CD34 positivity, after immunomagnetic separation in two different experiments. As monoclonal antibody for labeling the sample, the fluorescein isothiocyanate (FITC)- conjugated mouse anti-human CD34 MoAb (Mylteni) was used. Briefly, cell samples were incubated in the presence of the indicated MoAbs, at the proper dilution, in PBS containing 5% FCS and 1% Fc receptor (FcR) blocking reagent (Miltenyi) for 30 min at 4 C. Cells were then washed twice, resuspended with PBS and analyzed by a Coulter Epics XL (Coulter Electronics Inc., Hialeah, FL, USA) flow cytometer. only use Non-commercial 1 Supplementary Table S1. Complete list of the datasets used in this study and their sources. GEO Total samples Geo selected GEO accession of used Platform Reference series in series samples samples GSM142565 GSM142566 GSM142567 GSM142568 GSE6146 HG-U133A 14 8 - GSM142569 GSM142571 GSM142572 GSM142574 GSM51391 GSM51392 GSE2666 HG-U133A 36 4 1 GSM51393 GSM51394 only GSM321583 GSE12803 HG-U133A 20 3 GSM321584 2 GSM321585 use Promyelocytes_1 Promyelocytes_2 Promyelocytes_3 Promyelocytes_4 HG-U133A 8 8 3 GSE64282 Promyelocytes_5 Promyelocytes_6 Promyelocytes_7 Promyelocytes_8 Non-commercial 2 Supplementary Table S2. Chromosomal regions up-regulated in CD34+ samples as identified by the LAP procedure with the two-class statistics coded in the PREDA R package and an FDR threshold of 0.5. Functional enrichment analysis has been performed using DAVID (http://david.abcc.ncifcrf.gov/)
    [Show full text]
  • Supp Material.Pdf
    Simon et al. Supplementary information: Table of contents p.1 Supplementary material and methods p.2-4 • PoIy(I)-poly(C) Treatment • Flow Cytometry and Immunohistochemistry • Western Blotting • Quantitative RT-PCR • Fluorescence In Situ Hybridization • RNA-Seq • Exome capture • Sequencing Supplementary Figures and Tables Suppl. items Description pages Figure 1 Inactivation of Ezh2 affects normal thymocyte development 5 Figure 2 Ezh2 mouse leukemias express cell surface T cell receptor 6 Figure 3 Expression of EZH2 and Hox genes in T-ALL 7 Figure 4 Additional mutation et deletion of chromatin modifiers in T-ALL 8 Figure 5 PRC2 expression and activity in human lymphoproliferative disease 9 Figure 6 PRC2 regulatory network (String analysis) 10 Table 1 Primers and probes for detection of PRC2 genes 11 Table 2 Patient and T-ALL characteristics 12 Table 3 Statistics of RNA and DNA sequencing 13 Table 4 Mutations found in human T-ALLs (see Fig. 3D and Suppl. Fig. 4) 14 Table 5 SNP populations in analyzed human T-ALL samples 15 Table 6 List of altered genes in T-ALL for DAVID analysis 20 Table 7 List of David functional clusters 31 Table 8 List of acquired SNP tested in normal non leukemic DNA 32 1 Simon et al. Supplementary Material and Methods PoIy(I)-poly(C) Treatment. pIpC (GE Healthcare Lifesciences) was dissolved in endotoxin-free D-PBS (Gibco) at a concentration of 2 mg/ml. Mice received four consecutive injections of 150 μg pIpC every other day. The day of the last pIpC injection was designated as day 0 of experiment.
    [Show full text]
  • NATURAL KILLER CELLS, HYPOXIA, and EPIGENETIC REGULATION of HEMOCHORIAL PLACENTATION by Damayanti Chakraborty Submitted to the G
    NATURAL KILLER CELLS, HYPOXIA, AND EPIGENETIC REGULATION OF HEMOCHORIAL PLACENTATION BY Damayanti Chakraborty Submitted to the graduate degree program in Pathology and Laboratory Medicine and the Graduate Faculty of the University of Kansas in partial fulfillment ofthe requirements for the degree of Doctor of Philosophy. ________________________________ Chair: Michael J. Soares, Ph.D. ________________________________ Jay Vivian, Ph.D. ________________________________ Patrick Fields, Ph.D. ________________________________ Soumen Paul, Ph.D. ________________________________ Michael Wolfe, Ph.D. ________________________________ Adam J. Krieg, Ph.D. Date Defended: 04/01/2013 The Dissertation Committee for Damayanti Chakraborty certifies that this is the approved version of the following dissertation: NATURAL KILLER CELLS, HYPOXIA, AND EPIGENETIC REGULATION OF HEMOCHORIAL PLACENTATION ________________________________ Chair: Michael J. Soares, Ph.D. Date approved: 04/01/2013 ii ABSTRACT During the establishment of pregnancy, uterine stromal cells differentiate into decidual cells and recruit natural killer (NK) cells. These NK cells are characterized by low cytotoxicity and distinct cytokine production. In rodent as well as in human pregnancy, the uterine NK cells peak in number around mid-gestation after which they decline. NK cells associate with uterine spiral arteries and are implicated in pregnancy associated vascular remodeling processes and potentially in modulating trophoblast invasion. Failure of trophoblast invasion and vascular remodeling has been shown to be associated with pathological conditions like preeclampsia syndrome, hypertension in mother and/or fetal growth restriction. We hypothesize that NK cells fundamentally contribute to the organization of the placentation site. In order to study the in vivo role of NK cells during pregnancy, gestation stage- specific NK cell depletion was performed in rats using anti asialo GM1 antibodies.
    [Show full text]
  • 140503 IPF Signatures Supplement Withfigs Thorax
    Supplementary material for Heterogeneous gene expression signatures correspond to distinct lung pathologies and biomarkers of disease severity in idiopathic pulmonary fibrosis Daryle J. DePianto1*, Sanjay Chandriani1⌘*, Alexander R. Abbas1, Guiquan Jia1, Elsa N. N’Diaye1, Patrick Caplazi1, Steven E. Kauder1, Sabyasachi Biswas1, Satyajit K. Karnik1#, Connie Ha1, Zora Modrusan1, Michael A. Matthay2, Jasleen Kukreja3, Harold R. Collard2, Jackson G. Egen1, Paul J. Wolters2§, and Joseph R. Arron1§ 1Genentech Research and Early Development, South San Francisco, CA 2Department of Medicine, University of California, San Francisco, CA 3Department of Surgery, University of California, San Francisco, CA ⌘Current address: Novartis Institutes for Biomedical Research, Emeryville, CA. #Current address: Gilead Sciences, Foster City, CA. *DJD and SC contributed equally to this manuscript §PJW and JRA co-directed this project Address correspondence to Paul J. Wolters, MD University of California, San Francisco Department of Medicine Box 0111 San Francisco, CA 94143-0111 [email protected] or Joseph R. Arron, MD, PhD Genentech, Inc. MS 231C 1 DNA Way South San Francisco, CA 94080 [email protected] 1 METHODS Human lung tissue samples Tissues were obtained at UCSF from clinical samples from IPF patients at the time of biopsy or lung transplantation. All patients were seen at UCSF and the diagnosis of IPF was established through multidisciplinary review of clinical, radiological, and pathological data according to criteria established by the consensus classification of the American Thoracic Society (ATS) and European Respiratory Society (ERS), Japanese Respiratory Society (JRS), and the Latin American Thoracic Association (ALAT) (ref. 5 in main text). Non-diseased normal lung tissues were procured from lungs not used by the Northern California Transplant Donor Network.
    [Show full text]
  • Identification of Key Genes and Crucial Modules Associated with Coronary Artery Disease by Bioinformatics Analysis
    INTERNATIONAL JOURNAL OF MOLECULAR MEDICINE 34: 863-869, 2014 Identification of key genes and crucial modules associated with coronary artery disease by bioinformatics analysis XUEMEI ZHANG1, XIAOSHU CHENG2, HUIFENG LIU3, CHUNHUA ZHENG1, KUNRUI RAO1, YI FANG1, HAIRONG ZHOU4 and SHENGHE XIONG1 1Department of Cardiology, Third Affiliated Hospital of Nanchang University; 2Department of Cardiology, Second Affiliated Hospital of Nanchang University, Nanchang, Jiangxi 330008; 3Department of Cardiology, Xiaolan People's Hospital, Zhongshan, Guangdong 528415; 4Department of Cardiology, Second People's Hospital, Mudanjiang, Heilongjiang 157000, P.R. China Received December 19, 2013; Accepted May 28, 2014 DOI: 10.3892/ijmm.2014.1817 Abstract. The aim of this study was to identify key genes encoding GRB2 and the heat shock 70 kDa protein 8 (HSPA8). associated with coronary artery disease (CAD) and to explore Moreover, the chemokine and focal adhesion signaling path- the related signaling pathways. Gene expression profiles ways were selected based on their relative abundance in CAD. of 110 CAD and 112 non-CAD, healthy patients [CAD The GRB2 and HSPA8 proteins, as well as the chemokine and index (CADi) >23 and =0, respectively] were downloaded from focal adhension signaling pathways, might therefore be critical the Gene Expression Omnibus (GEO) database (accession: for the development of CAD. GSE12288). The differentially expressed genes (DEGs) in CAD were identified usingt -tests, and protein-protein interac- Introduction tion (PPI) networks for these DEGs were constructed using the Search Tool for the Retrieval of InteractiNg Genes (STRING) Coronary artery disease (CAD), also called coronary arte- database. The Database for Annotation, Visualization and riosclerosis, is the most common type of heart disease (1).
    [Show full text]
  • Supplementary Table 1
    Supplementary Table 1. Large-scale quantitative phosphoproteomic profiling was performed on paired vehicle- and hormone-treated mTAL-enriched suspensions (n=3). A total of 654 unique phosphopeptides corresponding to 374 unique phosphoproteins were identified. The peptide sequence, phosphorylation site(s), and the corresponding protein name, gene symbol, and RefSeq Accession number are reported for each phosphopeptide identified in any one of three experimental pairs. For those 414 phosphopeptides that could be quantified in all three experimental pairs, the mean Hormone:Vehicle abundance ratio and corresponding standard error are also reported. Peptide Sequence column: * = phosphorylated residue Site(s) column: ^ = ambiguously assigned phosphorylation site Log2(H/V) Mean and SE columns: H = hormone-treated, V = vehicle-treated, n/a = peptide not observable in all 3 experimental pairs Sig. column: * = significantly changed Log 2(H/V), p<0.05 Log (H/V) Log (H/V) # Gene Symbol Protein Name Refseq Accession Peptide Sequence Site(s) 2 2 Sig. Mean SE 1 Aak1 AP2-associated protein kinase 1 NP_001166921 VGSLT*PPSS*PK T622^, S626^ 0.24 0.95 PREDICTED: ATP-binding cassette, sub-family A 2 Abca12 (ABC1), member 12 XP_237242 GLVQVLS*FFSQVQQQR S251^ 1.24 2.13 3 Abcc10 multidrug resistance-associated protein 7 NP_001101671 LMT*ELLS*GIRVLK T464, S468 -2.68 2.48 4 Abcf1 ATP-binding cassette sub-family F member 1 NP_001103353 QLSVPAS*DEEDEVPVPVPR S109 n/a n/a 5 Ablim1 actin-binding LIM protein 1 NP_001037859 PGSSIPGS*PGHTIYAK S51 -3.55 1.81 6 Ablim1 actin-binding
    [Show full text]
  • Gene Ontology Functional Annotations and Pleiotropy
    Network based analysis of genetic disease associations Sarah Gilman Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy under the Executive Committee of the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2014 © 2013 Sarah Gilman All Rights Reserved ABSTRACT Network based analysis of genetic disease associations Sarah Gilman Despite extensive efforts and many promising early findings, genome-wide association studies have explained only a small fraction of the genetic factors contributing to common human diseases. There are many theories about where this “missing heritability” might lie, but increasingly the prevailing view is that common variants, the target of GWAS, are not solely responsible for susceptibility to common diseases and a substantial portion of human disease risk will be found among rare variants. Relatively new, such variants have not been subject to purifying selection, and therefore may be particularly pertinent for neuropsychiatric disorders and other diseases with greatly reduced fecundity. Recently, several researchers have made great progress towards uncovering the genetics behind autism and schizophrenia. By sequencing families, they have found hundreds of de novo variants occurring only in affected individuals, both large structural copy number variants and single nucleotide variants. Despite studying large cohorts there has been little recurrence among the genes implicated suggesting that many hundreds of genes may underlie these complex phenotypes. The question
    [Show full text]