Abstract
Background Antipsychotics suppress expression of inflammatory cytokines and inducible inflammatory enzymes. Elopiprazole (a phenylpiperazine antipsychotic drug in phase 1) has been characterized as a therapeutic drug to treat SARS-CoV-2 infection in a repurposing study. We aim to investigate the potential effects of aripiprazole (an FDA approved phenylpiperazine) on COVID19-related immunological parameters.
Methods Differential gene expression profiles of non-COVID versus COVID RNA-Seq samples (CRA002390 project in GSA database) and drug-naïve patients with psychosis at baseline and after three months of aripiprazole treatment was identified. An integrative analysis between COVID and aripiprazole immunomodulatory antagonist effects was performed.
Findings 82 out the 377 genes (21.7%) with expression significantly altered by aripiprazole have also their expression altered in COVID-19 patients and in 93.9% of these genes their expression is reverted by aripiprazole. The number of common genes with expression altered in both analyses is significantly higher than expected (Fisher’s Exact Test, two tail; P value=3.2e-11). 11 KEGG pathways were significantly enriched with genes with altered expression both in COVID-19 patients and aripiprazole medicated schizophrenia patients (P adj<0.05). The most significant pathways were associated to the immune system such as the “inflammatory bowel disease (IBD)” (the most significant pathway with a P adj of 0.00021), “Th1 and Th2 cell differentiation” and “B cell receptor signaling pathway”, all three related to the defense against infections.
Interpretation This exploratory investigation may provide further support to the notion that protective effect is exerted by phenylpiperazine by modulating the immunological dysregulation associated to COVID-19. Along with many ongoing studies and clinical trials, repurposing available medications could be of use in countering SARS-CoV-2 infection, but require further studies and trials.
1. Introduction
The SARS-CoV-2 epidemic has become the greatest challenge facing medicine today. Infected patients present with a wide range of clinical severity varying from asymptomatic to fatal condition (Wu et al., 2020). Advanced age, gender (male) and suffering comorbidities (diabetes, cardiovascular or chronic respiratory diseases) are risk factors for higher clinical severity, hospitalization rate and death from COVID-19 (Rubino et al., 2020; Zhou et al., 2020).
The presence of these comorbidities may decrease resilience and lower the ability to tolerate additional cytokine storm (Mangalmurti & Hunter, 2020). The high prevalence of these comorbidities among individuals with severe mental illnesses should increase the risk of COVID-19 complications and therefore might prone to a poorer outcome. However, somewhat unexpectedly low incidence and lack of severity of SARS-CoV-2 infection among people with schizophrenia have been reported (de Lusignan et al., 2020; van der Meer et al., 2020).
COVID-19 individuals who become critically and fatally ill seem to experience an indiscriminate and runaway immune response with an unchecked systemic overproduction of cytokines and immunological disbalance(Bhaskar et al., 2020; Manjili et al., 2020; Zeng et al., 2020).
Although, COVID-19 related immunopathogenesis is not understood just as an emergent cytokine storm but also as an impairment of protective T cell immunity (Chen et al., 2020).
Anti-inflammatory drugs such as dexamethasone have been shown to reduce deaths (Vabret et al., 2020). The crucial role of NLRP3 inflammasome activation in the pathogenesis of diseases caused by SARS-CoVs draws also attention toward potential role of its inhibitors in the treatment of COVID-19 (van den Berg & Te Velde, 2020).
The last decade has drastically increased the evidence linking psychosis spectrum disorders and altered immune function (Radhakrishnan et al., 2017). An imbalance in the host immune response seems to be associated with the activation of microglia in the pathophysiology of schizophrenia (Barichello et al., 2020). Antipsychotics were largely reported to suppress expression of inflammatory cytokines and inducible inflammatory enzymes (i.e., cyclooxygenase) and microglia activation (Dinesh et al., 2020). These anti-inflammatory effects are elicited via the reduction of proinflammatory cytokines production, modulating monocytes response through TLR and the inhibition of the microglial activation by reducing the levels of inducible nitric oxide synthase (iNOS), IL-1β, IL-6, and TNF-α16-18. The immunomodulatory effect of risperidone (pyridopyrimidines) and aripiprazole (marketed phenylpiperazine) has been demonstrated (Juncal-Ruiz et al., 2018), although aripiprazole had greater anti-inflammatory effect on TNF-α, IL-13, IL-17α and fractalkine. Wide range of different drug classes, such as cancer therapeutics, antipsychotics, and antimalarials, seem to have a beneficial effect against MERS and SARS coronaviruses (Dyall et al., 2017). A recent paper (Weston et al., 2020) described that, although infection cannot be prevented, chlorpromazine protects mice from severe clinical disease and SARS-CoV. Clozapine (atypical antipsychotic) has revealed to be effective in suppressing the proinflammatory cytokine expression by limiting the NLRP3 inflammasome activation in an in vitro model of schizophrenia (Giridharan et al., 2020). Interestingly, a research(L. Riva et al., 2020) investigating approximately 12,000 drugs in clinical-stage or Food and Drug Administration (FDA)-approved small molecules to identify candidate drugs to treat COVID-19, reported that elopiprazole (a never marketed phenylpiperazine antipsychotic drug) was listed among the 21 most potent compounds to inhibit SARS-CoV infection. Thus, the protective effect of phenylpiperazine marketed antipsychotics (aripiprazole) against a pernicious cytokine storm is a hypothesis that warrants further investigation.
The aim of the present study was to examine the potential immunomodulatory effects of antipsychotics by analyzing the profile of gene expression of drug-naïve patients with psychosis at baseline and after three months of treatment with aripiprazole (a phenylpiperazine antipsychotic drug).
2 Materials and Methods
2.1 Setting and Sample study
The cohort analyzed to study the effect of aripiprazole was obtained at the University Hospital Marques de Valdecilla (Cantabria, Spain). Conforming to international standards for research ethics, this study was approved by the Cantabria Ethics Institutional Review Board (IRB). Patients meeting inclusion criteria for a first episode of psychosis (drug-naïve) and their families provided written informed consent to be included in the study. After informed consent was signed, patients were included in a prospective, randomized, flexible-dose, open-label study (Crespo-Facorro et al., 2017; Mayoral Van-Son et al, 2020).
2.2 Laboratory assessments
Blood samples were obtained from 57 fasting schizophrenia subjects from 8:00 to 10:00 a.m. by the same staff and in the same setting. A detailed description of methodology followed to assess biochemical variables is available under request. None of the patients had a chronic inflammation or infection, or were taking medication that could influence the results of blood tests.
2.3 RNA extraction
Total RNA was extracted from blood using the Tempus™ Blood RNA Tube and the Tempus™ Spin RNA Isolation Kit (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s protocols. To select only high-quality RNA, the RNA integrity number (RIN) was characterized with a Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and samples with a RIN of at least 7.6 were used.
2.4 RNA next-generation sequencing
Total RNA was extracted from peripheral blood of each individual. The messenger RNA (mRNA) obtained from blood was sequenced at the Centro Nacional de Análisis Genómico (CNAG) using Illumina HiSeq instruments (San Diego, CA, USA). The mRNA was isolated from the total RNA and was fragmented once transformed into complementary DNA (cDNA). Fragments of 300bp on average were selected to construct the cDNA libraries for sequencing. Pair-end sequences of 70 nucleotides for each end were produced. The mRNA from blood samples of 57 drug-naïve psychotic patients at baseline and after 3 months of continuous treatment with aripiprazole was sequenced.
Sequence files were aligned to the GRCh38 human reference genome (Gencode release 25) using the STAR aligner(Dobin et al., 2013; Harrow et al., 2012). Reads count were normalized with the Voom algorithm using the cyclic loess method (Law et al., 2014) and significant gene expression changes between treated and naïve patients were identified with limma(Smyth et al., 2005). We performed a paired analysis using a Wald test and a parametric fit type with an adjusted p-value cutoff of 0.01.
2.5 Wuhan COVID19 dataset
COVID data was downloaded from the GSA server(Wang et al., 2017) with the CRA002390 identifier (Xiong et al., 2020). This research collected RNA-Seq samples of peripheral blood mononuclear cells (PBMC) from three COVID-19 patients and three healthy donors that were included in this study. These data were analyzed with RaNA-Seq (Prieto & Barrios, 2019) cloud platform and differential expression genes were detected by means of DESeq2 (Love et al., 2014) cloud using a Wald test, a parametric fit type and setting an adjusted p-value cutoff of 0.01.
2.6 Statistical and bioinformatic Analysis
Statistical significance of differential expression co-occurrence between COVID and Aripiprazole studies was calculated with a Fisher Exact test. Functional enrichment analyses, to identify biological pathways in the Biosystem database (downloaded on January 2020)(Geer et al., 2010) with a significant presence of differential expressed genes, were carried out with a Fisher Exact test. Visualization and final representation of pathways were performed with the pathview R package(Luo & Brouwer, 2013).
3 Results
Integrative transcriptomic analyses of aripiprazole effects on differentially expressed genes in COVID-19
3.1 COVID-19 patients (Wuhan dataset) versus controls
We found 2,137 genes with significant differential expression between COVID-19 patients and controls (P adj value <0.05). The most significant gene was the Charcot-Leyden crystal galectin gene (CLC) with a P adj value of 7.8e-23.
3.2 Drug-naïve schizophrenia patients at entry and after 3 months of aripiprazole treatment
We found 377 genes with significant differential expression before and after medication (P adj value <0.05). The two most significant genes were the LIM domain only 4 (LMO4) and the ATP binding cassette subfamily A member 9 (ABCA9) with a P adj value of 0.0039.
3.3 KEGG pathways significantly enriched with common genes with altered expression in COVID-19 patients and schizophrenia patients medicated with aripiprazole
We found that 11 KEGG pathways were significantly enriched with genes with altered expression both in COVID-19 patients and aripiprazole medicated schizophrenia patients (P adj<0.05). The most significant pathway was the Inflammatory bowel disease (IBD) with a P adj of 0.00021.
4 Discussion
Counterintuitively, the prevalence of psychosis among severely ill SARS-CoV-2 infected individuals is reduced. Severe COVID-19 infection requiring hospitalization does not appear to be associated to schizophrenia spectrum disorders (de Lusignan et al., 2020; van der Meer et al., 2020). Diverse speculative explanations such as higher degree of social isolation among psychosis individuals have been proposed to justify this lack of cases of COVID-19 among this population. The fact that antipsychotics have a demonstrated effect on immunological pathway may lead to the speculation about the beneficial effects of this drugs on controlling the acute hyper-inflammatory response that may be responsible for critical illness(Dinesh et al., 2020).
The present comparison of the genes with altered expression in COVID-19 patients versus controls and schizophrenia patients before and after medication with aripiprazole reveals that 82 out the 377 genes (21.7%) altered by aripiprazole medication are significantly altered also in the cohort of COVID-19 patients. The number of common genes to both analyses is significantly higher than expected by chance (Fisher’s Exact Test, two tail; P value=3.2e-11).
Interestingly, out of the 82 genes with expression altered in both analyzed cohorts 55 genes have decreased expression after aripiprazole medication and increased expression in COVID-19 patients; also, out of the 82 genes common to both cohorts 22 have increased expression after aripiprazole medication and decreased expression in COVID-19 patients. In total 77 genes out of 82 (93.9%) have altered expression in different direction when we compared the effects of COVID-19 and aripiprazole medication, suggesting that aripiprazole may have an effect that reverse many of the changes in expression caused by the coronavirus infection.
Analysis of pathways for enrichment of the 77 common genes with altered expression in COVID-19 and aripiprazole medicated patients shows 11 pathways significantly enriched (P value Fisher <0.05). Several of those pathways are related to the immune system such as the “inflammatory bowel disease (IBD)” (the most significant pathway), “Th1 and Th2 cell differentiation” and “B cell receptor signaling pathway”, both related to the defense against infections. Based on actual evidence showing hyperinflammation as well as T cell deficiencies and coagulation abnormalities, associated with life-threatening organ dysfunction, rather than with a mere hyper-inflammatory disease (Elhusseiny et al., 2020; G. Riva et al., 2020).
COVID-19 individuals who become critically and fatally ill experience an unchecked systemic overproduction of cytokines (Cytokines storm) and immuno-pathology has been suggested as a primary driver of morbidity and mortality with COVID-19 (Vabret et al., 2020). In addition to anti-inflammatory drugs (e.g., dexamethasone) other strategies including repurposing available medications aimed at quelling the inflammatory storm in critically ill patients (Guy et al., 2020).
A recent study (L. Riva et al., 2020) profiled a library of approximately 12,000 drugs in clinical-stage or FDA-approved small molecules to identify candidate drugs to treat COVID-19. In the list of the 21 most potent compounds to inhibit infection validated in dose response across multiple cell lines there are two antipsychotic drugs, elopiprazole and 8-(3-Chlorostyryl) caffeine which are in phase I and preclinical stages respectively. Elopiprazole and aripiprazole belong to the class of organic compounds known as phenylpiperazines. Phenylpiperazines are compounds containing a phenylpiperazine skeleton, which consists of a piperazine bound to a phenyl group.
Conclusions
The global effort to find therapies and develop vaccines capable of stopping the spread, and end the COVID-19 pandemic is a priority. This combined observational and exploratory investigations may provide further support to the notion that protective effect exerted by phenylpiperazines by modulating the immunological dysregulation associated to COVID-19. Along with many ongoing studies and clinical trials, repurposing available medications could be of use in countering SARS-CoV-2 infection, but require further studies and trials.
Data Availability
The cohort to study the effect of aripiprazole on gene expression was obtained from Cantabria Biobank. The patterns of gene expression of non-COVID and COVID were obtained from CRA002390 project in GSA database
Author’s contributions
Study concept and design: BC-F, MR-V, CP and JS; acquisition, analysis or interpretation of data: all authors; drafting of the manuscript: all authors; critical revision of the manuscript for important intellectual content: all authors; obtained funding: BC-F and JS; study supervision: BC-F, MR-V, CP and JS.
Role of Funding Sources
The present study was part of a larger prospective longitudinal study, the “First Episode Psychosis Clinical Program 10” (PAFIP10) study. ClinicalTrials.gov Identifiers: NCT02200588, NCT03481465, and NCT03476473.
No pharmaceutical industry or institutional sponsors participated in the study conception and design, data collection, analysis and interpretation of the results, or drafting of the manuscript.
Disclosure of potential conflict of interest
Prof. Crespo-Facorro has received unrestricted research funding from Instituto de Salud Carlos III, MINECO, Gobierno de Cantabria, Spanish Network for Research in Mental Health (CIBERSAM), from the 7th European Union Framework Program and Lundbeck. He has also received honoraria for his participation as a consultant and/or as a speaker at educational events from Janssen Johnson & Johnson, Mylan, Lundbeck, and Otsuka Pharmaceuticals.
Dr. Ruiz-Veguilla has received unrestricted research funding from Instituto de Salud Carlos III. He has also received honoraria for his participation as a consultant and/or as a speaker at educational events from Janssen, Lundbeck, and Otsuka Pharmaceuticals.
Dr. Vázquez-Bourgon has received unrestricted research funding from Instituto de Investigación Marqués de Valdecilla (IDIVAL). He has also received honoraria for his participation as a consultant and/or as a speaker at educational events from Janssen-Cilag and Lundbeck.
Carlos Prieto, Nathalia Garrido-Torres, Ana C. Sánchez-Hidalgo, and Jesus Sainz declare no conflicts of interest.
Dr. José Miguel Cisneros has received honoraria as a speaker from Novartis, Astellas Pharma, Pfizer, MSD, Janssen Pharmaceuticals, and AstraZeneca, outside the submitted work. He has also received report grants from Instituto de Salud Carlos III, Spanish Government, co-financed by the European Development Regional Fund “A way to achieve Europe”, during the conduct of the study.
Acknowledgements
We are highly indebted to the participants and their families for their cooperation in this study. We also thank IDIVAL biobank (Inés Santiuste and Jana Arozamena) for clinical samples and data as well as the PAFIP members (Marga Corredera) for the data collection. This work was supported by: SAF2016-76046-R and SAF2013-46292-R (MINECO and FEDER) to B.C.F.
We kindly thank all clinical staff at the Hospital Universitario Virgen del Rocio for support to collect clinical records and provide clinical care to COVID-19 patients.
We also kindly thank Dra. Marisa Barrigon for helpful discussions regarding clinical data analysis, and Idalino Rocha for manuscript editing and formatting.
Footnotes
↵† These authors might share senior authorship
benedicto.crespo.sspa{at}juntadeandalucia.es
npdunit{at}gmail.com
javier.vazquez{at}scsalud.es
anasanhid{at}gmail.com
nagatorr{at}gmail.com
josem.cisneros.sspa{at}juntadeandalucia.es
cprietos{at}usal.es
sainzjv{at}gmail.com
Section on Authors contributions updated to clarify role of author. Author affiliations updated to add ORCID numbers References revised and corrected