Keywords
schizophrenia, miRNA, interleukin-2, cytokines, SNP, expression, pathway
schizophrenia, miRNA, interleukin-2, cytokines, SNP, expression, pathway
Schizophrenia (OMIM code: 181500) is a severe psychiatric disorder with a heterogeneous clinical phenotype1. While the etiology of this disorder remains largely unknown, it has become evident that immune-inflammatory processes play an important multilevel role in disease development and progression2–4. Both cellular and humoral components of the innate and adaptive immune system were shown to be altered at different stages of disease development, both locally within the central nervous system as well as at a systemic level. Cytokines and chemokines are essential signal mediators of the immune system that modulate and guide immune/inflammatory responses as well as perform a wide range of other functions related to cell survival, proliferation, differentiation and migration5. Considering their biological importance, cytokines, especially interleukins (IL), have received considerable attention in the context of schizophrenia. The association of interleukins and other cytokines and their receptors with schizophrenia has been previously reported both on the level of genetic variants as well as the levels of gene expression and protein abundance, indicating their essential role in disease predisposition, development, progression and treatment response6–9.
IL-2 controls a wide range of biological activities, largely depending on the biological context. It is essential for T lymphocyte proliferation and differentiation, but it is also implicated in the generation and maintenance of regulatory T (Treg) cells10,11. The role of the IL2 gene in schizophrenia is still unclear and conflicting. However, previous studies, including our own findings, suggest the involvement of IL-2 in the pathogenesis of schizophrenia. Decreased lymphocyte production of IL-2 and increased IL-2 receptors have been reported previously12–16. Interleukin-2 receptor gamma (IL2RG) is an important signaling component of receptors for many cytokines, including IL-2, -4, -7, -9, -15 and -2113. Moreover, the IL2RG gene is over-expressed in the blood of schizophrenia patients13. On the other hand, other groups reported increased IL-2 serum levels in schizophrenia17–19. Some IL2 genetic polymorphisms were also reported to be associated with schizophrenia20,21. Expression of the IL2 gene is controlled at multiple layers. The IL2 gene contains at least two cis elements for transcript stability regulation, located in both the 3’ and 5’ untranslated regions (UTRs)22,23. Single nucleotide polymorphisms (SNPs) in the promoter region of IL2 influence the expression levels of this cytokine. Finally, expression of IL2 is also regulated by micro-RNAs24.
Micro-RNAs (miRNAs) are class of small, non-coding RNAs (comprised of about 22 nucleotides). MiRNAs are found in animals, plants and some viruses. They function in the regulation of gene expression at posttranscriptional level and RNA silencing25–27. As miRNAs are involved in the normal functioning of eukaryotic cells, deregulation of miRNAs has become associated with disease. There is a manually curated “miR2Disease” database, which aims to provide a comprehensive resource of microRNA deregulation in various human diseases28. Also, a number of studies have reported altered miRNA expression in schizophrenia, bipolar disorder and major depression and anxiety disorders29–32.
The aim of our study was to explore the possible association of miR-31, miR-146a, miR-181c and miR-155 with schizophrenia pathogenesis, as well as their link to IL2 gene expression in disease. This study is the first to report genetic association between schizophrenia and mentioned above miRNAs; however, several studies have reported a role for the IL2 gene in schizophrenia. All of these four micro-RNAs play a major role in regulating expression of the cytokine network. Particularly, miR-31, miR-146a and miR-181c are regulators of IL2 gene expression33–35, while miR-155 expression is greatly enhanced following stimulation of macrophages and dendritic cells by Toll-like receptors36. We also studied the possible association of the IL2 rs2069778 SNP genotype with IL2 and miRNA levels. This SNP was chosen due to its high minor allele frequency, clinical significance in autoimmune diseases37, as well as its location near the regulatory elements of the IL2 gene.
Informed written consent was obtained from all study participants. The study has been approved by the Ethical Committee of the Institute of Molecular Biology of the National Academy of Sciences RA (IRB00004079, IORG0003427).
This case-control study was conducted from January 2016 to February 2017. A total of 225 patients with paranoid schizophrenia (SCZ) and 225 sex- and age-matched controls (CTRL) with no family history of schizophrenia were involved in this study (Table 1). This was the maximum available number of schizophrenia patients in Armenia who agreed to participate in this study. From these subjects, 61 patients and 60 controls were tested for micro-RNA expression and 66 patients and 99 controls for IL2 gene expression. There was no specific criteria for dividing subsets in this study; subjects were divided into groups according to availability of biological material (DNA and RNA). All subjects were genotyped for this study.
SCZ | CTRL | |
---|---|---|
Age (mean ± SD) | 46.34±13.78 | 46.16±15.01 |
Sex male/female | 132/93 | 137/88 |
Age of first onset of schizophrenia | 30.83±12.92 | NA |
Paranoid schizophrenia (OMIM code: 181500, ICD-10-CM code: F20.0, DSM-5 code: 295.90) was diagnosed by two independent psychiatrists. Schizophrenia patients were recruited from the clinics of the Psychiatric Medical Center of the Ministry of Health of the Republic of Armenia (MH RA). Healthy subjects with any psychiatric illness during their lifetime, any serious endocrine or neurological disorder, any treatment or medical condition known to affect the brain or meeting the DSM-5 criteria for intellectual disability38 were excluded from this study. Exclusion criteria for all study subjects included any treatment with immune-modulating drugs and serious medical disorder.
Healthy control subjects were recruited among the blood donors of the Erebouni Medical Center (MH RA) and were interviewed by psychiatrists.
A total of 10 ml of peripheral blood was collected in EDTA containing tubes (5ml for RNA and 5ml for DNA isolation) from each study participant.
Peripheral blood mononuclear cells were isolated from whole blood using the following protocol, as described in 13: 10 ml of Red Cell Lysis Buffer (RCLB) (containing 0.144M ammonium chloride, 1 mM sodium bicarbonate) was added to 5 ml of fresh blood. After 5 minutes, the mixture was centrifuged at 1000g for 10 minutes, the supernatant was discarded and the pellet was gently rinsed with RCLB. The pellet was re-suspended in 5 ml of the RCLB buffer and centrifuged at 1000g for 10 minutes. The final purified pellet was stored in RNAlater (Ambion, Austin, TX, USA) at -20°C for later use. Total RNA was extracted using High Pure miRNA Isolation Kit (Cat. No. 05080576001 Roche Applied Science, Penzberg, Germany) according to the manufacturer’s instructions. The quantity and quality of RNA and DNA samples were assessed using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s protocol and stored at -80C until further use. The average RNA yield was 25 mg per 5 ml of blood sample.
The method for genomic DNA isolation was a modification of Miller’s salting-out procedure39 where a chloroform extraction phase is added (2ml of chloroform added to each tube supernatant and centrifuged at 3000g for 10 minutes)40.
cDNA synthesis was performed by reverse transcription (0.5 μg total RNA, total volume of cDNA 20 μl) using Transcriptor First Strand cDNA Synthesis Kit (Cat. No. 04897030001, Roche Applied Science) with anchored dT primers (0.4 μg; ABgene, Waltham, MA, USA) at 47°C for 45 min. cDNA samples were stored at -20°C until further use.
mRNA levels of IL2 (target) and PSMB2 (housekeeping) genes and were measured by quantitative real-time polymerase chain reaction with the Rotor Gene 3000 instrument (Corbett Research, USA). cDNA (5 μl, corresponding to 20 ng calculated on input total RNA) was added to 20 μl PCR-mix. The final reaction mix contained 900 nM of each sense and antisense primer (Roche Applied Science), 100 nM LNA probe (Roche Applied Science), 3.5 mM MgCl2, 200 μM of each dNTP (ABgene), 0,2μl Thermo-Start TAQ polymerase with concentration 5 U/μl and 1x ThermoStart Buffer (ABgene).
qPCR was performed using following thermal cycling conditions: initial denaturation for 15 minutes at 95°C; 40 cycles of denaturation for 20 seconds at 95°C, annealing for 30 seconds at 60°C and extension for 60 seconds at 72°C; and a final extension for five minutes at 72°C.
The primers and fluorescently-labeled Locked Nucleic Acid (LNA) probe from the Universal Probe Library for the IL2 gene were selected using the Probe Finder web-based software as follows:
IL2 gene:
Left primer: 5'-AAG TTT TAC ATG CCC AAG AAG G-3'
Right primer: 5'AAG TGA AAG TTT TTG CTT TGA GCT A-3'
Probe: #65
PSMB2 gene:
Left primer 5'-GTG AGA GGG CAG TGG AAC TC-3'
Right primer 5'-GAA GGT TGG CAG ATT CAG GA-3'
Probe #50
Reverse transcription for selected micro-RNAs and measurement of micro-RNA expression by quantitative real-time polymerase chain reaction (RT-PCR) were performed using TaqMan Micro-RNA Assays (see Table 2) and TaqMan Universal PCR Master MIX II (no UNG) (Cat. No. 4440040, Thermo Fisher Scientific, Waltham, USA), according to the manufacturer’s instructions. qPCR was performed using Realist-DX IAB real-time PCR system (GeneTiCA, Czech Republic).
qPCR was performed using following thermal cycling conditions: polymerase activation for 10 minutes at 95°C; 40 cycles of denaturation for 15 seconds at 95°C and annealing/extension for 60 seconds at 60°C.
We used miR-16 as a housekeeping micro-RNA. Data were expressed as arbitrary units (miR-X/miR-16 ratio). miR-16 is currently one of generally accepted housekeeping micro-RNAs for normalizing micro-RNA expression in blood cells by qRT PCR41. In addition, there are no reports on miR-16 alterations in any psychiatric disease including schizophrenia.
DNA samples of all patients and controls were genotyped for IL2 rs2069778 SNP using a polymerase chain reaction with sequence-specific primers (PCR-SSP)42. The sequences of specific primers were designed based on relevant DNA sequences available in the NCBI GenBank database (RefSeq Accession: NG_016779.1). Primer sequences for the selected SNP were as follows:
Reverse standard: 5’-CAC CAC TAC AAA TTC TAC AAA TTC G-3'
Reverse mutant: 5’-CAC CAC TAC AAA TTC TAC AAA TTC A-3'
Forward constant: 5’-CTG GTG CCA GAA AGA GCT TG-3'
The presence/absence of allele-specific amplicons were visualized by electrophoresis using 2% agarose gel in 0.5x Tris-Borate-EDTA (TBE) buffer stained with ethidium bromide fluorescent dye. To check the reproducibility of results, randomly selected DNA samples of study subjects (10% of total) were genotyped twice.
Genotyping was carried out at Laboratory of Human Genomics and Immunomics, Institute of Molecular Biology NAS RA (Yerevan, Armenia).
The Shapiro–Wilk test for normality revealed non-parametric distribution of the obtained data. Therefore, the significance of difference in gene expression levels between each study group was analyzed by the Mann–Whitney U test. p-values less than 0.05 were considered as significant. Statistical analysis was performed using GraphPad Prism (version 5) software. Allele and genotype frequencies were checked for Hardy-Weinberg equilibrium and were in equilibrium. No investigation of potential sources of bias was undertaken.
For characterization of enriched functions and biological pathways of the studied miRNAs targets, we used miRsystem43, which performs enrichment analyses, accounting both for target genes as well as the levels of individual miRNA expression.
In total, 66 SCZ patients (male/female: 33/33, mean age±S.D.: 51±11.2 years) and 99 healthy controls (male/female: 45/44, mean age±S.D.: 50±13.9 years) participated in IL2 gene expression step44.
We studied mRNA expression levels of the IL2 gene in schizophrenia and its possible association with the IL2 rs2069778 C/T SNP. The median mRNA expression levels of IL2 in the patient group were significantly lower than in healthy control subjects (patients vs. controls, median [interquartile range]: 0.06889 [0.6499–0.007519] vs. 1.469 [3.858–0.000], p=0.0095) (Figure 1).
Furthermore, analysis revealed a significant difference in IL2 expression between those with and without the IL2 rs2069778 C/T SNP in both the control and schizophrenia groups. Particularly, in schizophrenic patients, IL2 mRNA expression levels were significantly higher in rs2069778*T minor allele carriers (CT+TT) than in CC homozygotes (CC vs. CT+TT, median [interquartile range]: 0.033711 [0.5433–000453] vs. 0.2178 [2.618–0.03726], p=0.0003) (Figure 2). The same difference was found in control groups (CC vs. CT+TT, median [interquartile range]: 1.083 [2.840–0.000] vs. 2.625 [4.966–0.000], p=0.0495) (Figure 2). It is worth noting that expression of IL2 in schizophrenic patients carrying the CC genotype was lower than in corresponding controls (p=0.0216), while the difference in expression between carriers of the T minor allele was not significant (p=0.22).
A total of 61 SCZ patients (male/female: 37/24, mean age±S.D.: 45.4±13.9 years) and 60 healthy controls (male/female: 37/23, mean age±S.D.: 44.5±13.6 years) were tested for micro-RNA expression.
Median expression levels of all studied miRNAs were significantly higher in schizophrenic patients as compared to healthy controls (Table 3 and Figure 3)45.
Further analysis indicated that in T allele carriers, miR-181c had significantly lower expression compared to CC homozygous variants, both in schizophrenia patients (TT+CT vs. CC, median [interquartile range]: 1.99 [2.92–1.39] vs. 3.46 [4.84–2.14], p=0.0045) and controls (TT+CT vs. CC, median [interquartile range]: 0.41 [2.41–0.10] vs. 1.90 [5.31–0.59], p=0.011) (Figure 4).
CC – homozygous for rs2069778 SNP, TT+CT – carrier for rs2069778 SNP.
Interestingly, in T allele carriers, miR-31 also had significantly lower expression compared to CC homozygous variants in controls (TT+CT vs. CC, median [interquartile range]: 0.09 [1.61–0.03] vs. 1.84 [3.35–0.48], p=0.0015) but in schizophrenia patients there was no significant difference (TT+CT vs. CC, median [interquartile range]: 3.93 [7.74–2.91] vs. 5.03 [8.12–3.16], p=0.36) (Figure 4). We have not found any association between IL2 rs2069778 variants and expressions of miR-155 and miR-146a.
Our analysis demonstrated significant up-regulation of IL2 expression-modulating miRNAs in schizophrenia. However, it is known that miRNAs can affect multiple targets. For characterization of enriched functions and biological pathways of the studied miRNAs targets, we used a freely available online integrated system called miRsystem43. The analysis resulted in 30 pathways from 4 databases (KEGG46, Biocarta, Reactome, Pathway Interaction Database) significantly enriched with miRNA targets. The majority of these were related to immune/inflammatory system pathways where IL2 is either an effector or a target gene (Table 4). Figure 5 is an example of miR target enriched pathway46.
Red boxes indicate target genes for studied miRNAs. This figure has been reproduced with permission from the KEGG PATHWAY database46.
In this study, we observed decreased levels of IL2 expression in peripheral blood mononuclear cells of schizophrenic patients, paralleled with increased expression of IL2-regulating miRNAs (miR-31, miR-146a, miR-155 and miR-181c). In addition, we demonstrated that carriage of the minor allele for IL2 rs2069778 is associated with increased IL2 expression levels which might suggest either a regulatory role for this SNP or a linkage with other SNPs that can modulate gene expression.
Conflicting results on the levels of IL2 have been reported for this disease both in medicated and undedicated patients. The data published by Singh et al.47 and Theodoropoulou et al.48 are in line with our findings. However, Ebrinc et al.49 and Zhang et al.50,51 have found elevated levels of IL2 in their subjects. Furthermore, there are studies reporting no difference in the IL2 levels of schizophrenic patients as compared to controls52. Because expression was studied in different populations, we can speculate that the observed discrepancies could be partially explained by different genetic backgrounds. Moreover, the validity of our results is supported by the detected increase of IL2 expression-modulating miRNAs measured using independent assay techniques. Finally, consistent with IL2’s role in maintenance of Treg cells, their low levels were also reported in schizophrenia53.
Though the studies of miRNA involvement in schizophrenia are a relatively new direction, there are already results implicating miRNA deregulation in the pathogenesis of schizophrenia. miR-137 is the micro-RNA best known for its role in schizophrenia pathogenesis54–56. This micro-RNA is also well known due to a genetic polymorphism (SNP variant) in its gene, which was described as a genetic risk factor for the development of schizophrenia in a European population57–59. There are other recent studies which confirmed the role of distinct micro-RNAs such as miR-195, miR-181b, miR-301a, miR-19, miR-206, miR-30a and miR-219 in the pathogenesis of schizophrenia60–64. In this study, we reported four miRNAs (miR-31, miR-146a, miR-155, miR-181c) that were up-regulated in schizophrenic patients. Besides targeting IL2 expression, these molecules have many other targets that are involved in immune/inflammatory pathways, confirming the essential role of immune system disturbances in disease development and progression.
The limitation of the present study is the inability to recruit medication-free patients for assessment of the effect of treatment on IL2 and miRNA expression. However, in many studies cited in this paper50,51, regardless the direction of difference in IL2 levels, no differences were observed between treated and untreated patient groups. It should also be noted that we measured IL2 gene and miRNA expression in two different patient groups with little overlap, which prevented us from performing direct correlation analysis between the levels of IL2 and miRNAs.
Overall, our findings further strengthen the role of immune system deregulation in the development and progression of schizophrenia and necessitate further research towards understanding the changes of the Th1/Th2/T-reg response in this disease and in response to antipsychotic treatment.
All studied miRNAs (miR-31, miR-146a, miR-155, miR-181c) were over-expressed in schizophrenic patients, suggesting a role for them in disease pathogenesis. IL2 gene expression was down-regulated in schizophrenic patients. The IL2 rs2069778 C/T SNP is not associated with schizophrenia but regulates expression of the IL2 gene. Abnormalities in studied miRNA expressions result in the deregulation of the T-Cell receptor signaling pathway in schizophrenia.
Figshare: F1000 1990 Ghazaryan et al..xlsx. https://doi.org/10.6084/m9.figshare.9816860.v144.
Figshare: F1000 1990 Ghazaryan et al. Raw Data miRNA.xlsx. https://doi.org/10.6084/m9.figshare.10012442.v345.
Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
We thank the administration and medical staff of the Psychiatric Medical Centers of the Ministry of Health of the Republic of Armenia for selection of patients and healthy subjects.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
References
1. Pandey G, Ren X, Rizavi H, Zhang H: Proinflammatory cytokines and their membrane-bound receptors are altered in the lymphocytes of schizophrenia patients. Schizophrenia Research. 2015; 164 (1-3): 193-198 Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: Neurology of mood disorders, schizophrenia and suicide.
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
References
1. Alural B, Genc S, Haggarty S: Diagnostic and therapeutic potential of microRNAs in neuropsychiatric disorders: Past, present, and future. Progress in Neuro-Psychopharmacology and Biological Psychiatry. 2017; 73: 87-103 Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: Pharmacogenomics
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 1 10 Dec 19 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)