A Review on the Expression Pattern of Non-coding RNAs in Patients With Schizophrenia: With a Special Focus on Peripheral Blood as a Source of Expression Analysis

Schizophrenia is a destructive neuropsychiatric disease with a median prevalence of 4.0 per 1,000 during the whole life. Genome-wide association studies have shown the role of copy number variants (generally deletions) and certain alleles of common single nucleotide polymorphisms in the pathogenesis of schizophrenia. This disorder predominantly follows the polygenic inheritance model. Schizophrenia has also been linked with various alterations in the transcript and protein content of the brain tissue. Recent studies indicate that alterations in non-coding RNAs (ncRNAs) signature underlie a proportion of this dysregulation. High throughput microarray investigations have demonstrated momentous alterations in the expression of long non-coding RNAs (lncRNA) and microRNAs (miRNAs) in the circulation or post-mortem brain tissues of patients with schizophrenia compared with control samples. While Gomafu, PINT, GAS5, TCONS_l2_00021339, IFNG-AS1, FAS-AS1, PVT1, and TUG1 are among down-regulated lncRNAs in schizophrenia, MEG3, THRIL, HOXA-AS2, Linc-ROR, SPRY4-IT1, UCA1, and MALAT1 have been up-regulated in these patients. Moreover, several miRNAs, such as miR-30e, miR-130b, hsa-miR-130b, miR-193a-3p, hsa-miR-193a-3p, hsa-miR-181b, hsa-miR-34a, hsa-miR-346, and hsa-miR-7 have been shown to be dysregulated in blood or brain samples of patients with schizophrenia. Dysregulation of these transcripts in schizophrenia not only provides insight into the pathogenic processes of this disorder, it also suggests these transcripts could serve as diagnostic markers for schizophrenia. In the present paper, we explore the changes in the expression of miRNAs and lncRNAs in patients with schizophrenia.

Schizophrenia is a destructive neuropsychiatric disease with a median prevalence of 4.0 per 1,000 during the whole life. Genome-wide association studies have shown the role of copy number variants (generally deletions) and certain alleles of common single nucleotide polymorphisms in the pathogenesis of schizophrenia. This disorder predominantly follows the polygenic inheritance model. Schizophrenia has also been linked with various alterations in the transcript and protein content of the brain tissue. Recent studies indicate that alterations in non-coding RNAs (ncRNAs) signature underlie a proportion of this dysregulation. High throughput microarray investigations have demonstrated momentous alterations in the expression of long non-coding RNAs (lncRNA) and microRNAs (miRNAs) in the circulation or post-mortem brain tissues of patients with schizophrenia compared with control samples. While Gomafu, PINT, GAS5, TCONS_l2_00021339, IFNG-AS1, FAS-AS1, PVT1, and TUG1 are among down-regulated lncRNAs in schizophrenia, MEG3, THRIL, HOXA-AS2, Linc-ROR, SPRY4-IT1, UCA1, and MALAT1 have been up-regulated in these patients. Moreover, several miRNAs, such as miR-30e, miR-130b, hsa-miR-130b, miR-193a-3p, hsa-miR-193a-3p, hsa-miR-181b, hsa-miR-34a, hsa-miR-346, and hsa-miR-7 have been shown to be dysregulated in blood or brain samples of patients with schizophrenia. Dysregulation of these transcripts in schizophrenia not only provides insight into the pathogenic processes of this disorder, it also suggests these transcripts could serve as diagnostic markers for schizophrenia. In the present paper, we explore the changes in the expression of miRNAs and lncRNAs in patients with schizophrenia.
Keywords: lncRNA, miRNA, review, schizophrenia spectrum disorder, development of schizophrenia INTRODUCTION Schizophrenia is a destructive neuropsychiatric disease with a median prevalence of 4.0 per 1,000 during the whole life and a lifetime morbid risk of 7.2 per 1,000 (1). The onset of disorder is usually in youth or early adulthood (2) and rarely after the forties or in the childhood (3,4). In spite of comparable prevalence in both sexes (1), the disorder is usually begins earlier with a more severe course in male subjects (2,5). The disorder is associated with an increased mortality rate (1) due to suicide and cardiovascular comorbidities. Schizophrenia is diagnosed by the detection of symptoms in three major domains: positive symptoms, such as hallucination, delusion, and disorganized thought and movement, negative symptoms, i.e., lack of interest and enthusiasm, poverty of content of speech, lack of motivation or aptitude to do tasks and affective annihilation, and cognitive symptoms including deficiencies in executive function and attention (6). Genome-wide association studies have shown the role of copy number variants (generally deletions) and certain alleles of common single nucleotide polymorphisms (SNPs) in the pathogenesis of schizophrenia. This disorder predominantly follows the polygenic inheritance model with a significant level of overlap in genetic factors with other psychiatric disorders namely autism and bipolar disorder (7).
Schizophrenia has also been associated with dysregulation of several transcripts and proteins in the brain tissue. These alterations are caused by an intricate dysregulation of gene expression and protein synthesis. Both spatial and temporal elements contribute in determination of this dysregulation. Several studies indicate that alterations in non-coding RNAs (ncRNAs) signature underlie the mechanisms of this dysregulation (8). Changes in expression profile of ncRNAs in the brain tissues of patients with schizophrenia and the observed association between this disorder and certain SNPs in genomic regions coding these transcripts further highlight the role of ncRNAs in the pathogenesis of schizophrenia (8). These transcripts have been shown to govern the complex pattern of gene expression, thus being considered as one of the principal epigenetic mechanisms of gene expression (8). It is worth mentioning epigenetic alterations in some brain areas and neural tracks characterize an important route through which environmental parameters interplay with personal genetic composition to influence susceptibility to psychiatric disorders during the lifespan (9). Moreover, two classes of ncRNAs namely microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are extremely abundant in the human brain, implying their crucial role in the appropriate function of this tissue (10,11). While miRNAs mainly repress gene expression at post-transcriptional stage (12), lncRNAs can either enhance or suppress expression of gene through acting at chromatin, transcriptional an post-transcriptional levels (13). The important roles of these two classes of ncRNAs in the regulation of expression of genes and their developmental and tissue-specific signature suggest that they might underlie the observed aberrations in the transcriptome and proteome of brain tissue in patients with schizophrenia. Besides, alterations in the expression of lncRNAs might explain the plasticity of the organization of evolving neurons and their role in development of brain (14). In the current review, we explore the changes in the expression of miRNAs and lncRNAs in patients with schizophrenia. For the purpose of preparing this narrative review, we searched PubMed and google scholar with the key words "microRNA" or "miRNA" or "long non-coding RNA" or "lncRNA" AND "schizophrenia." Subsequently, we assessed the relevance of the obtained material through reading the full texts of the articles. Finally, we tabulated the retrieved data in distinct tables.

LncRNAs AND SCHIZOPHRENIA
LncRNAs are a huge and dissimilar group of ncRNA which have more than 200 nucleotides. LncRNAs signify the bulk of the noncoding transcriptome and utilize numerous mechanisms to exert their regulatory functions on gene expression among them being suppression or recruitment of transcription factors, modulation of chromatin structure and regulation of the stability of transcripts (15,16). High throughput microarray investigations have demonstrated significant changes in lncRNA signature in the peripheral blood or post-mortem brain tissues of patients with schizophrenia compared with control samples (17)(18)(19)(20).

LncRNA PROFILE IN CENTRAL TISSUES
Expression of the lncRNA Gomafu has been regulated by neuronal activation. This lncRNA has been shown to directly interact with two splicing factors namely QKI and SRSF1. Therefore, aberrant expression of Gomafu changes the splicing patterns of DISC1 and ERBB4 genes to a pattern which is similar to what is reported in schizophrenia. Expression of Gomafu is substantially decreased in post-mortem cortical gray matter of patients with schizophrenia (20). Level of Gomafu activity has been correlated with neuronal structural plasticity which is strong in the course of development and subsequently is decreased in the adulthood (14). Therefore, abnormal activity of Gomafu in the brain tissues of subjects with schizophrenia might reflect abnormal neurodevelopment. Moreover, abnormal levels of this lncRNA might also affect human behavior, since it has been associated with susceptibility to substance abuse (21).
Hu et al. have assessed RNA profile of post-mortem brain samples in patients with schizophrenia and those with bipolar disorder and control subjects. They reported differential expression of 20 long intergenic non-coding RNAs (lincRNAs) in orbitofrontal cortex of bipolar patients and aberrant expression of 34 and 1 lincRNAs in anterior cingulate cortex and dorsolateral prefrontal cortex of patients with schizophrenia, respectively. Thus, they reported brain area-specific profiles for lincRNAs. Differentially expressed lincRNAs were enriched in pathways namely immune system development and oligodendrocyte differentiation. Moreover, they reported altered DNA methylation as a possible mechanism for dysregulation of lincRNAs (18).

LncRNA PROFILE IN PERIPHERAL BLOOD
Chen et al. have examined expression profile of lncRNAs in the peripheral blood mononuclear cells (PBMCs) of patients with schizophrenia compared with healthy subjects. They reported differential expression of 125 lncRNAs between these two subgroups. Notably, expression levels of ENST00000394742, TCONS_l2_00025502, ENST00000563823, ENST00000521622, and TCONS_l2_00021339 were suggestively down-regulated in patients (17). LncRNA profiling has also revealed up-regulation of three lncRNAs in schizophrenia, down-regulation of six lncRNAs in major depressive disorder, and up-regulation of three lncRNAs in generalized anxiety disorder (GAD). Notably, lncRNAs observed to be up-regulated in schizophrenia were significantly decreased in patients with GAD. Furthermore, down-regulated lncRNAs in major depressive disorder were up-regulated in patients with schizophrenia. Finally, there were significant differences in the expression levels of lncRNAs between patients with schizophrenia and GAD. Therefore, a number of lncRNAs are putative biomarkers for differentiation of schizophrenia from major depressive disorder and generalized anxiety disorder (19). Sudhalkar et al. have examined expression levels of MEG3, PINT, and GAS5 in the PBMCs of patients with psychosis compared with healthy controls. They reported diagnostic differences with MEG3, PINT, and GAS5, and symptom acuity effect with MEG3 and GAS5. Moreover, there was significant difference in the expression of MEG3 between drug naïve patients and patients received risperidone (22). IFNG-AS1 expression has been shown to be down-regulated in patients with schizophrenia compared with healthy subjects in correlation with IFNG expression indicating a putative role for inflammation in this disorder (23). We have recently demonstrated downregulation of FAS-AS1, PVT1, and TUG1 in patients with schizophrenia compared with controls. Yet, expressions of GAS5, NEAT1, and OIP5-AS1 were similar between patients and controls (24).

Blood
The association between FAS-AS1 expression and schizophrenia was remarkable in a subgroup of men.
PVT1 PVT1 and TUG1 were appropriate biomarkers in male patients.

TUG1
Frontiers in Psychiatry | www.frontiersin.org schizophrenia compared with including 62 over-expressed and 63 under-expressed lncRNAs in patients. Antipsychotic treatment has resulted in reduction in NONHSAT089447 and NONHSAT041499 levels, parallel with decrease in the post-treatment Positive And Negative Syndrome Scale (PANSS) scores. Moreover, reduction in NONHSAT041499 levels have been associated with improvement of several clinical manifestations and better response to therapies (17). Chen et al. have reported the effects of olanzapine treatment in suppression of expression of the NONHSAT089447 lncRNA. Small interfering RNA-mediated NONHSAT089447 silencing has reduced expression of dopamine receptors DRD3 and DRD5. In addition, Western blot studies verified the role of this lncRNA in regulation of DRD signaling (25). We have demonstrated up-regulation of HOXA-AS2, Linc-ROR, MEG3, SPRY4-IT1, and UCA1 in patients with schizophrenia compared with healthy subjects. Yet, when assessing their expressions in sex-based subclasses, the differences in their expressions were significant just among females. Moreover, we reported correlations between expressions of Linc-ROR and SPRY4-IT1 and age of patients (26). Ni et al. have profiled peripheral blood transcriptome of monozygotic twins discordant for schizophrenia. Using this approach, authors have demonstrated up-regulation of AC006129.1 lncRNA in patients. This lncRNA regulates inflammatory reactions through promoting expression of SOCS3 and CASP1. Further experiments showed that AC006129.1 interacts with the promoter region of the transcriptional repressor Capicua (CIC) to enhance the interactions of DNA methyltransferases with its promoter and decrease CIC expression, thus reversing CIC-associated SOCS3 and CASP1 suppression. Activation of SOCS3 increases the anti-inflammatory reactions by obstructing JAK/STAT pathway (27). Table 2 exhibits the list of up-regulated lncRNAs in schizophrenia.

Association Between LncRNA(s) SNP(s) and Risk of Schizophrenia
Few studies have appraised association between lncRNAs SNPs and risk of schizophrenia. Rao et al. have performed a twophase association study on 8 tag SNPs that encompass the entire MIAT region in two independent cohorts from Han Chinese population. They demonstrated significant association between paranoid schizophrenia and the rs1894720. Moreover, there was a weak association between rs4274 and this condition. No specific haplotype was detected that modulate risk of paranoid schizophrenia in the assessed population (29).

Diagnostic Value of LncRNAs in Schizophrenia
Expression levels of lncRNAs could be used as diagnostic markers for schizophrenia. We have assessed this possibility in a cohort of Iranian patients with schizophrenia (24,26). Based on the obtained results in a limited number of patients, GAS5 and OIP5-AS1 have been proposed as appropriate biomarkers in female subjects (24). Table 3 summarizes the results of these studies.

miRNAs AND SCHIZOPHRENIA
These small transcripts have intricate temporospatial signature in the brain tissue and can modulate expressions of myriad of genes by working as the specificity elements for genesilencing apparatus in the cells. Based on the observed association between miRNA dysregulation and substantial alterations in the network organization in the course of neurodevelopment, miRNAs are considered as important regulators of several neurological processes (30). Numerous studies have indicated aberrant expression of miRNAs in the brain and peripheral blood of patients with schizophrenia. Sun et al. have demonstrated over-expression of miR-132, miR-195, miR-30e, and miR-7 in plasma samples of patients with schizophrenia, and up-regulation of miR-212, miR-34a, and miR-30e in their PBMCs {#164}. They suggested that miRNA signature is more distinctive in plasma samples compared with PBMCs (34). A high throughput miRNA profiling revealed up-regulation of eight miRNAs in plasma samples obtained from schizophrenia patients, among them were miR-130b and miR-193a-3p which were up-regulated in schizophrenia but not in non-schizophrenia psychotic disorders. These results indicated these miRNAs as state-independent biomarkers for schizophrenia (35).

miRNA LEVELS IN PERIPHERAL BLOOD
Lai et al. have subsequently assessed the impact of hospitalization on the expression levels of these miRNAs. Notably, expression of none of these miRNAs did not change after 2 months hospitalization of patients even when clinical symptoms were remarkably ameliorated. Thus, these miRNAs have been suggested as trait biomarkers instead of statedependent biomarkers. Assessment of expression profile of CT value of NONHSAT041499 was significantly higher in patients after the treatment, representing the substantial down-regulation of this lncRNA expression by the treatment. The symptomatology score and total score were meaningfully reduced following treatment.

NONHSAT098126
These transcripts have been suggested as markers for the diagnosis and prognostic evaluations.
hsa-miR-34a and hsa-miR-548d in post-mortem brain samples showed no difference between patients and controls (36).

miRNA PROFILE IN CENTRAL TISSUES/CELL LINES
Expression profile of miRNAs has also been assessed in the olfactory epithelium as one of the limited available neural tissues that have neurons and neural stem cells. Mor et al. have detected over-expression of miR-382-5p in cultured olfactory cells obtained from schizophrenia patients compared with controls. Up-regulation of this miRNA was also verified in microdissected olfactory epithelium neuronal tissues of these patients. However, miR-382 was not expressed in lymphoblastoid cell lines originated from either patients with schizophrenia or control subjects. This miRNA has been shown to target FGFR1 and SPRY4, two genes whose expressions were decreased in the olfactory cells obtained from patients with schizophrenia (37). Other experiments in HEK293 and SH-SY5Y cell lines have validated CALN1 as a target of miR-137 (38). A number of other miRNAs can also regulate development of dendrites, thus being implicated in the functionality of synapses and neuronal interactions. For instance, miR-214 has been shown to increase dendrite dimension, complexity and morphogenesis. Its involvement in the pathogenesis of schizophrenia has been suggested through the observed interaction between this miRNA and quaking (Qki), a candidate gene in schizophrenia. The miR-214-Qki axis is a fundamental axis in the modulation of dendritic development in the neurons (39). Experiments in a mouse neuroblastoma cell line (Neuro2A) have shown the regulatory role of miR-137 on the expression of several important proteins in the PI3K-Akt-mTOR axis which have functions downstream of neuregulin/ErbB and BDNF. Therefore, this miRNA controls neuronal reactions to these factors and dendritic development, thus contributing in the risk of schizophrenia (40). Figure 1 shows the mechanism of involvement of miR-214 in the pathogenesis of schizophrenia. Tables 4, 5 show the list of down-/up-regulated miRNAs in the schizophrenia patients, respectively.

miRNAS VARIATIONS IN SCHIZOPHRENIA
Some studies have verified associations between SNPs within miRNA-coding genes and risk of schizophrenia. For instance, in a genome-wide association study, Ripke et al. have shown significant association between the rs1625579 of the miR-137 and risk of schizophrenia (57). Notably, others have confirmed associations between miR-137 target genes CUB, CSMD1, C10orf26, CACNA1C, TCF4, and ZNF804A, and risk of schizophrenia (58)(59)(60)(61). Notably, the mentioned SNP in miR-137 has been shown to influence cognitive activity. The T/T genotype of this SNP is correlated with working memory defects in patients with schizophrenia as reflected by their reduced scores in the brief assessment of cognition in schizophrenia instrument (62). The risk genotype has also been associated with dorsolateral prefrontal cortex hyperactivation (DLPFC) in patients with schizophrenia (63). In addition, the risk allele of this SNP has been associated with down-regulation of miR-137 in schizophrenia and is potentially involved in the modulation of expression of the schizophrenia risk locus TCF4, more emphasizing on the participation of miR-137 and its downstream molecules in this disorder (64). These studies further support the role of this miRNA in the pathogenesis of schizophrenia.

DIAGNOSTIC VALUE OF miRNAs IN SCHIZOPHRENIA
Diagnostic value of miRNAs has been assessed in different biological sources of patients with schizophrenia using statistical methods, such as plotting the receiver operating characteristic curve and calculation of the area under the curve (AUC) values. Such method has shown the diagnostic values of 0.767 and 0.756 for plasma and PBMC expression patterns of miR-30e, respectively. Plasma levels of this miRNA had sensitivity and specificity of 90.90 and 60.00%, respectively. These values for its PBMC levels were 81.80 and 68.00%, FIGURE 1 | miR-214 is overexpressed in patients with schizophrenia. This miRNA binds with the 3 ′ UTR of Qki transcript to reduce its expression. Qki is an RNA-binding protein which modulates stability of several mRNAs among them are AIP1, Wnt2a, and Rab5. These transcripts are involved in the dendritic development (39). Qki can also bind with FEZ1 mRNA to enhance its stability (41). FEZ1 binds with DISC1 and increases neurite outgrowth in response to NGF (42). Qki can bind with intronic regions and participate in the biogenesis of circRNAs (43). Abnormal expression of Qki might be associated with the observed dysregulation of cricRNAs in schizophrenia (44).
respectively. Moreover, logistic regression analysis showed higher sensitivity of plasma levels of this miRNA in differentiating schizophrenia patients from healthy subjects compared with its PBMC levels (34). The diagnostic value of a panel of over-expressed miRNAs including miR-30e, miR-181b, miR-34a, miR-346, and miR-7 has been reported to be 0.713. Notably, suitable pharmacotherapy resulted in down-regulation of miR-132, miR-181b, miR-432, and miR-30e expressions. Besides, there was a remarkable correlation between amelioration of clinical symptoms and alterations in the expression levels of miR-132, miR-181b, miR-212, and miR-30e (54). A certain molecular axis including the transcription factor the early growth response protein 1 (EGR1), miR-30a-5p, and its target gene NEUROD1 has been shown to differentiate schizophrenia patients from healthy subjects with diagnostic accuracy of 0.962 which was far higher than the diagnostic power of miR-30a-5p alone. This axis has also been proved useful for monitoring of patients in acute psychotic phase (65). Table 6 shows the results of studies which assessed diagnostic role of miRNAs in the schizophrenia.

DISCUSSION
Schizophrenia is complex disorder caused by interaction between several genomic loci. The speculation that this disorder is caused by one or a few common principal gene effects has been experimentally examined in genome-wide linkage studies yet results generally showed no genome-wide significance (66). In addition to genomic variants that contribute in the pathogenesis of this disorder, some other mechanisms might modify the risk of development of schizophrenia. Expressions of ncRNAs are modulated by neuronal activation (20), suggesting a role for these transcripts in the pathophysiology of neuropsychiatric disorders. In the current study, we reviewed the literature about the role of ncRNAs in the pathophysiology of schizophrenia. Based on the above-mentioned evidence, several lncRNAs and miRNAs have been dysregulated in blood or brain samples of patients with schizophrenia. Not discounting the role of ncRNAs as biomarkers for schizophrenia, peripheral expression profile of these transcripts does not necessarily reflect their expression in the brain tissues (36). Aberrant expression of a number of Prefrontal cortex these transcripts has been associated with PANSS score (17). Moreover, clinical diagnoses of psychosis and symptom severity have been shown to alter expression of a number of lncRNAs (22) indicating a substantial role for these transcripts in the pathogenesis of schizophrenia. Although there is a considerable level of overlap between psychiatric disorders in the terms of contributing genetic factors, expression levels of some ncRNAs could be used to differentiate a number of these conditions (19). Moreover, expression profile of certain ncRNAs can distinguish patients with schizophrenia from healthy subjects. At the present, schizophrenia is principally diagnosed based on the clinical symptoms and signs instead of on the pathophysiological biomarkers (67). Identification of such biomarkers would facilitate detection of malingering, thus has practical significance in the forensic medicine in the assessment of cases pretending psychological disorders for a particular gain. Moreover, molecular biomarkers have implications in the establishment of targeted therapies. Notably,    expression profile of ncRNAs in the circulation has the potential to improve the diagnostic and prognostic assessment of patients with schizophrenia (67). A former study has demonstrated up-regulation of NONHSAT089447, NONHSAT021545, and NONHSAT041499 lncRNAs in patients with schizophrenia, while down-regulation of these transcripts in patients with GAD. Notably, authors have reported six lncRNAs with opposite expression patterns in schizophrenia and major depressive disorder. Moreover, they have identified three GADrelated lncRNAs whose expressions were significantly differences between patients with schizophrenia and GAD patients (28). However, most of these studies do not test if expression profile of ncRNAs can distinguish people with schizophrenia from people with other disorders, such as major depression, bipolar disorder, or autism. Based on the similarities in many of clinical symptoms between subjects with schizophrenia and other disorders, including bipolar disorder, autism, and major depression, identification of specific markers for each disease has practical significance. Therefore, this field should be explored in future investigations. Particularly, miRNA signature has been correlated with clinical course, patients' response to pharmacologic interventions and prognosis of patients with schizophrenia (54). Perhaps, the most intensively assessed miRNA in this regard is miR-30 family. Consistent with this observation, expression of EGR1 which regulates expressions of this miRNA family has been decreased in PBMCs of patients with schizophrenia. On the other hand, expression of NEUROD1 as a target gene of miR-30a-5p has been increased in these patients (65). Based on the results of this study, identification of transcription factor/ miRNA/ target gene axes would be a practical method for recognition of molecular pathways in the pathogenesis of schizophrenia and development of diagnostic/ prognostic panels for this disorder.
Several schizophrenia-associated ncRNAs have been shown to modulate immune responses and inflammatory pathways. Recent studies have highlighted the presence of intricate interplay between the immune system, systemic inflammatory responses, and the central nervous system, which can result alterations in mood, cognitive functions, and behavior. All of these aspects contribute in the pathogenesis of schizophrenia (68). Immune responses can influence activity of neurotransmitters as well as neurodegenerative and neurodevelopmental processes all of which are related with this disorder (68). Assessment of lncRNA signature in Amygdala samples from schizophrenia patients has further endorsed dysregulation of immune-associated lncRNAs in these samples (69). Thus, the functional axes between ncRNAs and immune-related genes provide an explanation for involvement of these transcripts in the pathophysiology of schizophrenia.
Besides, high throughput studies in schizophrenia subjects have shown a tendency toward global up-regulation of miRNAs normally abundant in infants, whereas down-regulation of those normally abundant in prepuberty. Therefore, dysregulation of miRNAs dynamic changes might be involved in the pathogenesis of schizophrenia (70). A number of ncRNAs might explain the observed dissimilarity in the brain activation modes between patients with schizophrenia and controls. For instance, Gomafu has been shown to be intensely regulated in reaction to neuronal activation (71). This lncRNA is also implicated in schizophreniarelated alternative splicing (71). Therefore, aberrant expression of this lncRNA might reflect the abnormal brain activity in these patients.
Finally, it is worth mentioning that miRNAs have distinctive roles as trait-dependent markers or state-dependent markers. This speculation is supported by amelioration of dysregulated expression of a number of miRNAs after pharmacotherapy (54), while no change in the expression of other miRNAs following suitable treatments (36). Therefore, it is necessary to define the role of each miRNA as trait-or state-dependent marker to design distinctive diagnostic/ prognostic miRNA panels for schizophrenia.
Taken together, lncRNAs and miRNAs are potential transcripts that can explain the difference in expression profile of protein coding genes in brain and blood tissues of patients with schizophrenia and healthy subjects. Moreover, as a number of above-mentioned ncRNAs are functionally related with dopamine neurotransmission, these ncRNAs might alter response of patients to some types of antipsychotic drugs. In addition, several ncRNAs can be used as disease markers in schizophrenia.
Several questions should be addressed about the role of ncRNAs in the development of schizophrenia. Studies reviewed in this article have mostly assessed expression of ncRNAs in adult patients. Although differentially expressed ncRNAs among patients and controls have been enriched in neurodevelopment, neurotransmission, and synaptic plasticity, the functional impact of these ncRNAs in the development of neurons should be assessed through knock-out/-in studies in animal models and cell lines.
Finally, most of studies reviewed here have not appraised the effects of antipsychotic drugs on expression of ncRNAs. The data regarding the therapeutic regimens and types of antipsychotic drugs have not been presented in the main articles. This is possibly because patients have been under treatment with different antipsychotic drugs. We mention this point as a limitation of these studies.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.