Dissecting the Genetic Architecture of Cystatin C in Diversity Outbred Mice

Total Page:16

File Type:pdf, Size:1020Kb

Dissecting the Genetic Architecture of Cystatin C in Diversity Outbred Mice MULTIPARENTAL POPULATIONS Dissecting the Genetic Architecture of Cystatin C in Diversity Outbred Mice M. Nazmul Huda,*,† Melissa VerHague,‡ Jody Albright,§ Tangi Smallwood,§ Timothy A. Bell,§ Excel Que,* Darla R. Miller,§ Baback Roshanravan,** Hooman Allayee,†† Fernando Pardo Manuel de Villena,§ and Brian J. Bennett*,†,1 *Obesity and Metabolism Research Unit, Western Human Nutrition Research Center, USDA, ARS, Davis, California 95616, †Department of Nutrition, University of California Davis, California 95616, ‡Nutrition Research Institute, University of North § Carolina Kannapolis, North Carolina 28081, Department of Genetics, University of North Carolina at Chapel Hill, North Carolina 27599, **Department of Medicine, Division of Nephrology, University of California, Davis, Davis, California 95616, and ††Departments of Preventive Medicine and Biochemistry & Molecular Medicine, Keck School of Medicine, University of Southern California, Los Angeles California 90033 ORCID IDs: 0000-0003-2775-4470 (M.N.H.); 0000-0002-5738-5795 (F.P.M.d.V.); 0000-0002-0766-3195 (B.J.B.) ABSTRACT Plasma concentration of Cystatin C (CysC) level is a biomarker of glomerular filtration rate in the KEYWORDS kidney. We use a Systems Genetics approach to investigate the genetic determinants of plasma CysC Quantitative trait concentration. To do so we perform Quantitative Trait Loci (QTL) and expression QTL (eQTL) analysis of loci 120 Diversity Outbred (DO) female mice, 56 weeks of age. We performed network analysis of kidney gene Multi parental expression to determine if the gene modules with common functions are associated with kidney biomarkers models of chronic kidney diseases. Our data demonstrates that plasma concentrations and kidney mRNA levels Cystatin C of CysC are associated with genetic variation and are transcriptionally coregulated by immune genes. Kidney Specifically, Type-I interferon signaling genes are coexpressed with Cst3 mRNA levels and associated biomarkers with CysC concentrations in plasma. Our findings demonstrate the complex control of CysC by genetic Type-I interferon polymorphisms and inflammatory pathways. signalling pathway Multiparent Advanced Generation Inter-Cross (MAGIC) multiparental populations MPP The kidney is a complex organ responsible for excretion (Finco 1997), sex, and race, limiting its accuracy to represent true GFR (Stevens secretion (Davis et al. 1961), reabsorption (Lemann et al. 1970), and et al. 2006). An alternative to creatinine is Cystatin C (CysC), often activating vitamin D (Fraser and Kodicek 1973). The gold standard used in research as the basis for estimating glomerular filtration rate. for assessing kidney function is the glomerular filtration rate (GFR), CysC is produced by all mammalian cells, secreted into the blood, but it is difficult to measure with precision. Therefore, GFR is often filtered through the glomerulus, and catabolized by tubular cells estimated from circulating creatinine (Ferguson et al. 2015). Creat- (Inker and Okparavero 2011). Plasma CysC had been found to be inine is an amino acid derivative, released by muscle, and freely unaltered by age, sex, race, and metabolic disorders and was proposed filtered by the kidney glomerulus (Narayanan and Appleton 1980). as a clinical biomarker for GFR (Filler et al. 2005; Newman et al. However, the level of plasma creatinine is influenced by a number of 1995). Several clinical trials reported CysC as a superior marker factors, including: diet, muscle mass, medication, chronic illness, age, compared to creatinine (Plebani et al. 1998; Kyhse-Andersen et al. 1994; Volume 10 | July 2020 | 2529 Harmoinen et al. 1999) while a few others could not show a the manufacturer’s instructions. In brief, plasma samples were diluted significant difference between CysC and creatinine (Donadio et al. 200-fold and incubated in a 96 well microplate pre-coated with CysC 2001; Oddoze et al. 2001). Although CysC is not influenced by specific antibody and CysC concentration was determined from the physiological factors, several recent studies (Köttgen et al. 2010; color intensity of oxidized Tetramethylbenzidine (TMB) measured at Köttgen et al. 2009; Hwang et al. 2007) have found associations 450 nm. Urinary total protein, creatinine and plasma blood urea between SNPs near the Cystatin C gene (Cst3) and plasma CysC levels nitrogen were measured by COBAS INTEGRA 400 plus analyzer or CysC-based estimated glomerular filtration rate (eGFRcys). How (Roche Diagnostics, Rotkreuz, Switzerland). To measure urinary Na+, these genetic variants relate to plasma CysC protein level or its mRNA samples were diluted 1500-fold with 1.0% v/v trace metal free nitric level remains to be determined. acid and analyzed by using a Varian VISTA AX CCD Simultaneous In the current study, we focus on the genetic determination of CysC Inductively Coupled Plasma Atomic Emission Spectroscopy (Varian, protein and mRNA levels using a “Systems Genetics” approach, which is CA, USA). Standards Na+ for Inductively Coupled Plasma Mass useful for the discovery of genes and pathways associated with a reduced Spectrometry (Spex CertiPrep, NJ, USA) was used to make standards kidney function (Keller et al. 2012). We also perform co-expression ranging from 0.05 ppm to 5.0 ppm. Certified urine controls analysis and expression QTL (eQTL) analysis to further investigate the (Seronorm, Stasjonsveien, Norway) and sample pool controls were control of CysC levels using the Diversity Outbred (DO) population, used to check the accuracy of the assay. which are derived from eight founder strains (A/J, C57BL/6J, 129S1/ SvImJ, NOD/ShiLtJ, NZO/HiLtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Kidney mRNA microarray Since the DO is a highly recombinant population with immense genetic Total RNA was extracted from about 25 mg of kidney tissue using and phenotypic variation (Svenson et al. 2012; Smallwood et al. 2014), automated instrumentation (Maxwell 16 Tissue LEV Total RNA DO mice have been successfully used in studies focused on the genetic Purification Kit, Promega). RNA concentration was measured by response to environmental toxin exposure and diseases (French et al. fluorometry (Picogreen Life Technologies), and RNA quality was 2015; Tyler et al. 2017).TheDOmicecapturethesamesetofallelic verified using a microfluidics platform (Bioanalyzer, Agilent Tech- variants as the eight founder strains, and their genetic backgrounds are nologies). 95 RNA samples were chosen for microarray analysis and well-studied, which make them an excellent model for researching the 1 sample was run in duplicate as a control. RNA was hybridized to genetic susceptibility to disease. In this study, we focus on the genetic Affymetrix Mouse Gene 2.1 ST 96-Array Plate using a GeneTitan variation associated with CysC gene expression and plasma CysC instrument from Affymetrix according to the manufacturer’s proto- concentration in DO mice. We then performed co-expression network cols. We used the robust multiarray average method (RMA) imple- analysis to identify gene modules associated with CysC and identified mented in the Affymetrix gene expression console with default two gene modules associated with CysC levels. settings (median polish and sketch-quantile normalization) to esti- mate the normalized expression levels of transcripts. All probes MATERIALS AND METHODS containing known single nucleotide polymorphisms (SNPs) from the eight founder inbred mouse strains of the DO mouse population Study design and sample collection were masked (165,204 probes) during normalization by downloading Female diversity outbred (DO) mice (n = 120; J: DO; G16; stock theSNPsfromtheSangersequencingwebsite(http://www.sanger.ac.uk/ number 009376) were obtained from the Jackson Laboratory (Bar science/data/mouse-genomes-project) and overlapping them with Harbor, ME) at 4 weeks of age. These mice were used to generate probe sequences. The data are available at NCBI Gene Expression progeny for another study (Didion et al. 2015) and then aged. At Omnibus (GEO) database under the accession ID GSE122061. 56 weeks of age, all mice were injected intraperitoneally with a volume of sterile isotonic saline equivalent to 10% of their body weight and Genotyping urine was collected in a chilled metabolic cage system (Hatteras Inc, DNA was extracted and purified from tail samples using Qiagen NC). On the following day, mice were fasted for 4 hr, plasma was DNeasy kit (Qiagen, MD, USA) according to the manufacturer’s collected from the retro-orbital sinus into EDTA-containing micro- instructions. Genotyping was performed using the Mega Mouse tubes and centrifuged. Mice were euthanized, dissected, and tissue Universal Genotyping Array (MegaMUGA) by GeneSeek (Neogen, samples were snap frozen in liquid nitrogen. Biological samples were Lansing, MI) (Welsh et al. 2012). The MegaMUGA array is built on ° stored at -80 until assayed. All procedures were approved by the the Illumina Infinium platform and contains 77,800 SNP markers IACUC at UNC, Chapel Hill (IACUC Protocol Number 13-103). that are distributed throughout the genome at an average spacing of 33 Kb. Genotypes of DO mice are accessible through UNC’s Plasma and urine biochemical assays Mutant Mouse Resource and Research Centers website (https:// Plasma CysC was measured by a commercially available quantitative www.med.unc.edu/mmrrc/genotypes/). We estimated heritability sandwich ELISA kit (R&D Systems, MN, USA) for mice according to (h2) from allele probability using a linear mixed model in R package QTL2 version 0.18 (Broman et al. 2019). Copyright
Recommended publications
  • Seq2pathway Vignette
    seq2pathway Vignette Bin Wang, Xinan Holly Yang, Arjun Kinstlick May 19, 2021 Contents 1 Abstract 1 2 Package Installation 2 3 runseq2pathway 2 4 Two main functions 3 4.1 seq2gene . .3 4.1.1 seq2gene flowchart . .3 4.1.2 runseq2gene inputs/parameters . .5 4.1.3 runseq2gene outputs . .8 4.2 gene2pathway . 10 4.2.1 gene2pathway flowchart . 11 4.2.2 gene2pathway test inputs/parameters . 11 4.2.3 gene2pathway test outputs . 12 5 Examples 13 5.1 ChIP-seq data analysis . 13 5.1.1 Map ChIP-seq enriched peaks to genes using runseq2gene .................... 13 5.1.2 Discover enriched GO terms using gene2pathway_test with gene scores . 15 5.1.3 Discover enriched GO terms using Fisher's Exact test without gene scores . 17 5.1.4 Add description for genes . 20 5.2 RNA-seq data analysis . 20 6 R environment session 23 1 Abstract Seq2pathway is a novel computational tool to analyze functional gene-sets (including signaling pathways) using variable next-generation sequencing data[1]. Integral to this tool are the \seq2gene" and \gene2pathway" components in series that infer a quantitative pathway-level profile for each sample. The seq2gene function assigns phenotype-associated significance of genomic regions to gene-level scores, where the significance could be p-values of SNPs or point mutations, protein-binding affinity, or transcriptional expression level. The seq2gene function has the feasibility to assign non-exon regions to a range of neighboring genes besides the nearest one, thus facilitating the study of functional non-coding elements[2]. Then the gene2pathway summarizes gene-level measurements to pathway-level scores, comparing the quantity of significance for gene members within a pathway with those outside a pathway.
    [Show full text]
  • Transcriptomic Signature and Metabolic Programming of Bovine Classical and Nonclassical Monocytes Indicate Distinct Functional Specializations
    bioRxiv preprint doi: https://doi.org/10.1101/2020.10.30.362731; this version posted November 1, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Transcriptomic signature and metabolic programming of bovine classical and nonclassical monocytes indicate distinct functional specializations Stephanie C. Talker1,2, G. Tuba Barut1,2, Reto Rufener3, Lilly von Münchow4, Artur Summerfield1,2 1Institute of Virology and Immunology, Bern and Mittelhäusern, Switzerland 2Department of Infectious Diseases and Pathobiology, Vetsuisse Faculty, University of Bern, Bern, Switzerland 3Institute of Parasitology, Vetsuisse Faculty, University of Bern, Bern, Switzerland 4 Bucher Biotec AG, Basel, Switzerland *Correspondence: Corresponding Author [email protected] Keywords: monocyte subsets, transcriptome, metabolism, cattle Abstract Similar to human monocytes, bovine monocytes can be split into CD14+CD16- classical and CD14-CD16+ nonclassical monocytes (cM and ncM, respectively). Here, we present an in-depth analysis of their steady-state transcriptomes, highlighting pronounced functional specializations. Gene transcription indicates that pro-inflammatory and antibacterial processes are associated with cM, while ncM appear to be specialized in regulatory/anti-inflammatory functions and tissue repair, as well as antiviral responses and T-cell immunomodulation. In support of these functional differences, we found that oxidative phosphorylation prevails in ncM, whereas cM are clearly biased towards aerobic glycolysis. Furthermore, bovine monocyte subsets differed in their responsiveness to TLR ligands, supporting an antiviral role of ncM. Taken together, these data clearly indicate a variety of subset-specific functions in cM and ncM that are likely to be transferable to monocyte subsets of other species, including humans.
    [Show full text]
  • In Silico Prediction and in Vitro Characterization of Multifunctional Human Rnase3
    Hindawi Publishing Corporation BioMed Research International Volume 2013, Article ID 170398, 12 pages http://dx.doi.org/10.1155/2013/170398 Research Article In Silico Prediction and In Vitro Characterization of Multifunctional Human RNase3 Pei-Chun Lien,1 Ping-Hsueh Kuo,1 Chien-Jung Chen,1 Hsiu-Hui Chang,1 Shun-lung Fang,1 Wei-Shuo Wu,2 Yiu-Kay Lai,2 Tun-Wen Pai,3 and Margaret Dah-Tsyr Chang1, 4 1 InstituteofMolecularandCellularBiology,NationalTsingHuaUniversity,No.101,Section2,KuangFuRoad, Hsinchu 30013, Taiwan 2 Institute of Biotechnology, National Tsing Hua University, No. 101, Section 2, Kuang Fu Road, Hsinchu 30013, Taiwan 3 Department of Computer Science and Engineering, National Taiwan Ocean University, 2 Pei Ning Road, Keelung 20224, Taiwan 4 Department of Medical Science, National Tsing Hua University, No. 101, Section 2, Kuang Fu Road, Hsinchu 30013, Taiwan Correspondence should be addressed to Margaret Dah-Tsyr Chang; [email protected] Received 31 October 2012; Accepted 2 December 2012 Academic Editor: Hao-Teng Chang Copyright © 2013 Pei-Chun Lien et al. is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Human ribonucleases A (hRNaseA) superfamily consists of thirteen members with high-structure similarities but exhibits divergent physiological functions other than RNase activity. Evolution of hRNaseA superfamily has gained novel functions which may be preserved in a unique region or domain to account for additional molecular interactions. hRNase3 has multiple functions including ribonucleolytic, heparan sulfate (HS) binding, cellular binding, endocytic, lipid destabilization, cytotoxic, and antimicrobial activities.
    [Show full text]
  • Unveiling the Multifaceted Antimicrobial Mechanism of Action of Human Host Defence Rnases
    ADVERTIMENT. Lʼaccés als continguts dʼaquesta tesi queda condicionat a lʼacceptació de les condicions dʼús establertes per la següent llicència Creative Commons: http://cat.creativecommons.org/?page_id=184 ADVERTENCIA. El acceso a los contenidos de esta tesis queda condicionado a la aceptación de las condiciones de uso establecidas por la siguiente licencia Creative Commons: http://es.creativecommons.org/blog/licencias/ WARNING. The access to the contents of this doctoral thesis it is limited to the acceptance of the use conditions set by the following Creative Commons license: https://creativecommons.org/licenses/?lang=en Departament de Bioquímica, Biología Molecular i Biomedicina Unveiling the multifaceted antimicrobial mechanism of action of human host defence RNases Tesis presentada por Javier Arranz Trullén para optar al grado de Doctor en Bioquímica, Biología Molecular y Biomedicina bajo la dirección de la Dra. Ester Boix Borràs y el Dr. David Pulido Gómez Unitat de Biociències del Departament de Bioquímica i Biología Molecular. Universitat Autònoma de Barcelona Dra. Ester Boix Borràs Dr. David Pulido Gómez Javier Arranz Trullén Campus de Bellaterra, Septiembre 2016 1 2 Biochemistry, Molecular Biology and Biomedicine Department Unveiling the multifaceted antimicrobial mechanism of action of human host defence RNases Thesis presented by Javier Arranz Trullén to obtain the PhD in Biochemistry, Molecular Biology and Biomedicine directed by Dr. Ester Boix Borràs and Dr. David Pulido Gómez Biosciences Unit of the Biochemistry and Molecular
    [Show full text]
  • Single-Cell Transcriptomes Reveal a Complex Cellular Landscape in the Middle Ear and Differential Capacities for Acute Response to Infection
    fgene-11-00358 April 9, 2020 Time: 15:55 # 1 ORIGINAL RESEARCH published: 15 April 2020 doi: 10.3389/fgene.2020.00358 Single-Cell Transcriptomes Reveal a Complex Cellular Landscape in the Middle Ear and Differential Capacities for Acute Response to Infection Allen F. Ryan1*, Chanond A. Nasamran2, Kwang Pak1, Clara Draf1, Kathleen M. Fisch2, Nicholas Webster3 and Arwa Kurabi1 1 Departments of Surgery/Otolaryngology, UC San Diego School of Medicine, VA Medical Center, La Jolla, CA, United States, 2 Medicine/Center for Computational Biology & Bioinformatics, UC San Diego School of Medicine, VA Medical Center, La Jolla, CA, United States, 3 Medicine/Endocrinology, UC San Diego School of Medicine, VA Medical Center, La Jolla, CA, United States Single-cell transcriptomics was used to profile cells of the normal murine middle ear. Clustering analysis of 6770 transcriptomes identified 17 cell clusters corresponding to distinct cell types: five epithelial, three stromal, three lymphocyte, two monocyte, Edited by: two endothelial, one pericyte and one melanocyte cluster. Within some clusters, Amélie Bonnefond, Institut National de la Santé et de la cell subtypes were identified. While many corresponded to those cell types known Recherche Médicale (INSERM), from prior studies, several novel types or subtypes were noted. The results indicate France unexpected cellular diversity within the resting middle ear mucosa. The resolution of Reviewed by: Fabien Delahaye, uncomplicated, acute, otitis media is too rapid for cognate immunity to play a major Institut Pasteur de Lille, France role. Thus innate immunity is likely responsible for normal recovery from middle ear Nelson L. S. Tang, infection. The need for rapid response to pathogens suggests that innate immune The Chinese University of Hong Kong, China genes may be constitutively expressed by middle ear cells.
    [Show full text]
  • Uniprot Acceprotiens 121 113 Ratio(113/12 114 Ratio
    Uniprot Acceprotiens 121 113 ratio(113/12 114 ratio(114/12 115 ratio(115/12 116 ratio(116/12 117 ratio(117/12 118 ratio(118/12 119 ratio(119/121) P02768 Serum albumin OS=Homo s666397.2 862466.6 1.29 593482.1 0.89 2220420.5 3.33 846469.3 1.27 634302.5 0.95 736961.1 1.11 842297.5 1.26 P02760 Protein AMBP OS=Homo s381627.7 294812.3 0.77 474165.8 1.24 203377.3 0.53 349197.6 0.92 346271.7 0.91 328356.1 0.86 411229.3 1.08 B4E1B2 cDNA FLJ53691, highly sim78511.8 107560.1 1.37 85218.8 1.09 199640.4 2.54 90022.3 1.15 73427.3 0.94 82722 1.05 102491.8 1.31 A0A0K0K1HEpididymis secretory sperm 3358.1 4584.8 1.37 4234.8 1.26 8496.1 2.53 4193.7 1.25 3507.1 1.04 3632.2 1.08 4873.3 1.45 D3DNU8 Kininogen 1, isoform CRA_302648.3 294936.6 0.97 257956.9 0.85 193831.3 0.64 290406.7 0.96 313453.3 1.04 279805.5 0.92 228883.9 0.76 B4E1C2 Kininogen 1, isoform CRA_167.9 229.7 1.37 263.2 1.57 278 1.66 326 1.94 265.5 1.58 290.3 1.73 341.5 2.03 O60494 Cubilin OS=Homo sapiens G40132.6 45037.5 1.12 38654.5 0.96 34055.8 0.85 39708.6 0.99 44702.9 1.11 45025.7 1.12 32701.3 0.81 P98164 Low-density lipoprotein rece40915.4 45344.8 1.11 35817.7 0.88 35721.8 0.87 42157.7 1.03 46693.4 1.14 48624 1.19 38847.7 0.95 A0A024RABHeparan sulfate proteoglyca46985.3 43536.1 0.93 49827.7 1.06 33964.3 0.72 44780.9 0.95 46858.6 1.00 47703.5 1.02 37785.7 0.80 P01133 Pro-epidermal growth factor 75270.8 73109.5 0.97 66336.1 0.88 56680.9 0.75 70877.8 0.94 76444.3 1.02 81110.3 1.08 65749.7 0.87 Q6N093 Putative uncharacterized pro47825.3 55632.5 1.16 48428.3 1.01 63601.5 1.33 65204.2 1.36 59384.5
    [Show full text]
  • Mechanisms of Sodium–Glucose Cotransporter 2 Inhibition
    Diabetes Care Volume 43, September 2020 2183 – Ele Ferrannini,1 Ashwin C. Murthy,2 Mechanisms of Sodium Glucose Yong-hoLee,3 ElzaMuscelli,1 SophieWeiss,4 Rachel M. Ostroff,4 Naveed Sattar,5 Cotransporter 2 Inhibition: Stephen A. Williams,4 and Peter Ganz6 Insights From Large-Scale Proteomics Diabetes Care 2020;43:2183–2189 | https://doi.org/10.2337/dc20-0456 OBJECTIVE To assess the effects of empagliflozin, a selective sodium–glucose cotransporter 2 (SGLT2) inhibitor, on broad biological systems through proteomics. RESEARCH DESIGN AND METHODS Aptamer-based proteomics was used to quantify 3,713 proteins in 144 paired plasma samples obtained from 72 participants across the spectrum of glucose tolerance before and after 4 weeks of empagliflozin 25 mg/day. The biology of the plasma proteins significantly changed by empagliflozin (at false discovery rate– corrected P < 0.05) was discerned through Ingenuity Pathway Analysis. RESULTS Empagliflozin significantly affected levels of 43 proteins, 6 related to cardiomyocyte function(fattyacid–bindingprotein3and4[FABPA],neurotrophic receptortyrosine kinase, renin, thrombospondin 4, and leptin receptor), 5 to iron handling (ferritin heavy chain 1, transferrin receptor protein 1, neogenin, growth differentiation CARDIOVASCULAR AND METABOLIC RISK factor 2 [GDF2], and b2-microglobulin), and 1 to sphingosine/ceramide metabolism 1CNR Institute of Clinical Physiology, Pisa, Italy (neutral ceramidase), a known pathway of cardiovascular disease. Among the 2 fi Cardiovascular Division, Department of Medi- protein changes achieving the strongest statistical signi cance, insulin-like binding cine, Hospital of the University of Pennsylvania, factor protein-1 (IGFBP-1), transgelin-2, FABPA, GDF15, and sulphydryl oxidase Philadelphia, PA 2 precursor were increased, while ferritin, thrombospondin 3, and Rearranged 3Department of Medicine, Severance Hospital, during Transfection (RET) were decreased by empagliflozin administration.
    [Show full text]
  • Differentially Expressed Genes in Aneurysm Tissue Compared With
    On-line Table: Differentially expressed genes in aneurysm tissue compared with those in control tissue Fold False Discovery Direction of Gene Entrez Gene Name Function Change P Value Rate (q Value) Expression AADAC Arylacetamide deacetylase Positive regulation of triglyceride 4.46 1.33E-05 2.60E-04 Up-regulated catabolic process ABCA6 ATP-binding cassette, subfamily A (ABC1), Integral component of membrane 3.79 9.15E-14 8.88E-12 Up-regulated member 6 ABCC3 ATP-binding cassette, subfamily C (CFTR/MRP), ATPase activity, coupled to 6.63 1.21E-10 7.33E-09 Up-regulated member 3 transmembrane movement of substances ABI3 ABI family, member 3 Peptidyl-tyrosine phosphorylation 6.47 2.47E-05 4.56E-04 Up-regulated ACKR1 Atypical chemokine receptor 1 (Duffy blood G-protein–coupled receptor signaling 3.80 7.95E-10 4.18E-08 Up-regulated group) pathway ACKR2 Atypical chemokine receptor 2 G-protein–coupled receptor signaling 0.42 3.29E-04 4.41E-03 Down-regulated pathway ACSM1 Acyl-CoA synthetase medium-chain family Energy derivation by oxidation of 9.87 1.70E-08 6.52E-07 Up-regulated member 1 organic compounds ACTC1 Actin, ␣, cardiac muscle 1 Negative regulation of apoptotic 0.30 7.96E-06 1.65E-04 Down-regulated process ACTG2 Actin, ␥2, smooth muscle, enteric Blood microparticle 0.29 1.61E-16 2.36E-14 Down-regulated ADAM33 ADAM domain 33 Integral component of membrane 0.23 9.74E-09 3.95E-07 Down-regulated ADAM8 ADAM domain 8 Positive regulation of tumor necrosis 4.69 2.93E-04 4.01E-03 Up-regulated factor (ligand) superfamily member 11 production ADAMTS18
    [Show full text]
  • Supplementary Table S2. List of Genes with Expression That Is Positively Correlated (Pearson Correlation Coefficient P>0.3)
    Supplementary Table S2. List of genes with expression that is positively correlated (Pearson correlation coefficient p>0.3) with HOXA9 expression in the study by Sun et al (1). Most HOXA genes (indicated by bold face) show highly significant positive correlation with HOXA9 expression, except for HOXA6 and HOXA13, similar to our findings in the UCSF and MDA tumor sets. Note that HOXA11, which did not demonstrate a statistically significant positive correlation with HOXA9 expression in both the UCSF and MDA tumor sets, demonstrates a substantially lower correlation coefficient relative to the other HOXA genes that positively correlate with HOXA9 expression in the study by Sun et al. These data were obtained from the online ONCOMINE database (www.oncomine.org) (2) searching for transcripts positively correlated with HOXA9 expression. Gene Name Gene Symbol Reporter ID Correlation (p) homeobox A9 HOXA9 214651_s_at .7682 homeobox A9 HOXA9 209905_at .7682 homeobox A10 HOXA10 213150_at .6058 homeobox A7 HOXA7 235753_at .6058 homeobox A10 HOXA10 213147_at .6058 homeobox A7 HOXA7 206847_s_at .6058 homeobox A4 HOXA4 206289_at .5741 homeobox A2 HOXA2 1557051_s_at .5379 homeobox A1 HOXA1 214639_s_at .5379 homeobox A3 HOXA3 235521_at .5379 homeobox B2 HOXB2 205453_at .5379 homeobox B3 HOXB3 228904_at .5379 EST EST 1555907_at .5379 homeobox A4 HOXA4 230080_at .5379 homeobox A2 HOXA2 228642_at .5379 homeobox A2 HOXA2 1557050_at .5379 homeobox A5 HOXA5 213844_at .5379 homeobox A2 HOXA2 214457_at .5113 homeobox B7 HOXB7 204778_x_at .5042 homeobox C6 HOXC6
    [Show full text]
  • Bioinformatics Analysis of Gene Expression Pro Les of Immune
    Bioinformatics Analysis of Gene Expression Proles of Immune-mediated Necrotizing Myopathy Jiaying Shi First Hospital of Shanxi Medical University Shule Wang First Hospital of Shanxi Medical University Jingfei Zhang First Hospital of Shanxi Medical University Xueli Chang First Hospital of Shanxi Medical University Juan Wang First Hospital of Shanxi Medical University Jing Zhang First Hospital of Shanxi Medical University Junhong Guo First Hospital of Shanxi Medical University Wei Zhang ( [email protected] ) Department of Neurology, The First Hospital, Shanxi Medical University, Taiyuan 030000, P. R. China Research Keywords: Immune-mediated necrotizing myopathy, bioinformatics analysis, cancer-related pathways, skeletal muscle construction, immune response Posted Date: January 8th, 2021 DOI: https://doi.org/10.21203/rs.3.rs-137686/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/21 Abstract Background: Immune-mediated necrotizing myopathy (IMNM) is a type of autoimmune myopathy with limited therapeutic measures. This study aims to elucidate the potential biomarkers and investigate the underlying mechanisms in IMNM. Materials and Methods: Microarray datasets in GSE128470 and GSE39454 were obtained from the Gene Expression Omnibus. Differentially expressed genes (DEGs) were ltrated by limma package in R statistical software. Functional enrichment analyses were performed using DAVID online tools. STRING database was used to construct protein‐protein interaction (PPI) networks. The module analysis and hub genes validation were performed using Cytoscape software. Results: Integrated analysis of two databases revealed 160 co-expressed DEGs in IMNM, including 80 downregulated genes and 80 upregulated genes. GO enrichment analysis revealed that sarcomere is the most signicantly enriched GO term within the DEGs.
    [Show full text]
  • Table S1. 103 Ferroptosis-Related Genes Retrieved from the Genecards
    Table S1. 103 ferroptosis-related genes retrieved from the GeneCards. Gene Symbol Description Category GPX4 Glutathione Peroxidase 4 Protein Coding AIFM2 Apoptosis Inducing Factor Mitochondria Associated 2 Protein Coding TP53 Tumor Protein P53 Protein Coding ACSL4 Acyl-CoA Synthetase Long Chain Family Member 4 Protein Coding SLC7A11 Solute Carrier Family 7 Member 11 Protein Coding VDAC2 Voltage Dependent Anion Channel 2 Protein Coding VDAC3 Voltage Dependent Anion Channel 3 Protein Coding ATG5 Autophagy Related 5 Protein Coding ATG7 Autophagy Related 7 Protein Coding NCOA4 Nuclear Receptor Coactivator 4 Protein Coding HMOX1 Heme Oxygenase 1 Protein Coding SLC3A2 Solute Carrier Family 3 Member 2 Protein Coding ALOX15 Arachidonate 15-Lipoxygenase Protein Coding BECN1 Beclin 1 Protein Coding PRKAA1 Protein Kinase AMP-Activated Catalytic Subunit Alpha 1 Protein Coding SAT1 Spermidine/Spermine N1-Acetyltransferase 1 Protein Coding NF2 Neurofibromin 2 Protein Coding YAP1 Yes1 Associated Transcriptional Regulator Protein Coding FTH1 Ferritin Heavy Chain 1 Protein Coding TF Transferrin Protein Coding TFRC Transferrin Receptor Protein Coding FTL Ferritin Light Chain Protein Coding CYBB Cytochrome B-245 Beta Chain Protein Coding GSS Glutathione Synthetase Protein Coding CP Ceruloplasmin Protein Coding PRNP Prion Protein Protein Coding SLC11A2 Solute Carrier Family 11 Member 2 Protein Coding SLC40A1 Solute Carrier Family 40 Member 1 Protein Coding STEAP3 STEAP3 Metalloreductase Protein Coding ACSL1 Acyl-CoA Synthetase Long Chain Family Member 1 Protein
    [Show full text]
  • Research Article in Silico Prediction and in Vitro Characterization of Multifunctional Human Rnase3
    Hindawi Publishing Corporation BioMed Research International Volume 2013, Article ID 170398, 12 pages http://dx.doi.org/10.1155/2013/170398 Research Article In Silico Prediction and In Vitro Characterization of Multifunctional Human RNase3 Pei-Chun Lien,1 Ping-Hsueh Kuo,1 Chien-Jung Chen,1 Hsiu-Hui Chang,1 Shun-lung Fang,1 Wei-Shuo Wu,2 Yiu-Kay Lai,2 Tun-Wen Pai,3 and Margaret Dah-Tsyr Chang1, 4 1 InstituteofMolecularandCellularBiology,NationalTsingHuaUniversity,No.101,Section2,KuangFuRoad, Hsinchu 30013, Taiwan 2 Institute of Biotechnology, National Tsing Hua University, No. 101, Section 2, Kuang Fu Road, Hsinchu 30013, Taiwan 3 Department of Computer Science and Engineering, National Taiwan Ocean University, 2 Pei Ning Road, Keelung 20224, Taiwan 4 Department of Medical Science, National Tsing Hua University, No. 101, Section 2, Kuang Fu Road, Hsinchu 30013, Taiwan Correspondence should be addressed to Margaret Dah-Tsyr Chang; [email protected] Received 31 October 2012; Accepted 2 December 2012 Academic Editor: Hao-Teng Chang Copyright © 2013 Pei-Chun Lien et al. is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Human ribonucleases A (hRNaseA) superfamily consists of thirteen members with high-structure similarities but exhibits divergent physiological functions other than RNase activity. Evolution of hRNaseA superfamily has gained novel functions which may be preserved in a unique region or domain to account for additional molecular interactions. hRNase3 has multiple functions including ribonucleolytic, heparan sulfate (HS) binding, cellular binding, endocytic, lipid destabilization, cytotoxic, and antimicrobial activities.
    [Show full text]