Introduction

Molecular Medicine Reports

1791-2997 1791-3004

D.A. Spandidos

10.3892/mmr.2020.11513

mmr-22-05-4340

Articles

Investigating aberrantly expressed microRNAs in peripheral blood mononuclear cells from patients with treatment-resistant schizophrenia using miRNA sequencing and integrated bioinformatics

You

1 * Zhang

Yunqiao

1 * Long

Qing

1 * Liu

Zijun

1 Ma

Xiao

1 Lu

Zixiang

2 Yang

Wei

3 Feng

Ziqiao

1 Zhang

Wengyu

1 Teng

Zhaowei

1 Zeng

Yong

1Research Management Department, The Sixth Affiliated Hospital, Kunming Medical University, Yuxi, Yunnan 653100, P.R. China 2Psychiatric Ward, Honghe Second People's Hospital, Honghe, Yunnan 654399, P.R. China 3Psychiatric Ward, Yuxi Second People's Hospital, Yuxi, Yunnan 653100, P.R. China

Correspondence to: Dr Yong Zeng or Dr Zhaowei Teng, Research Management Department, The Sixth Affiliated Hospital, Kunming Medical University, 21 Nieer Road, Yuxi, Yunnan 653100, P.R. China, E-mail: zengyong1388851037@163.com, E-mail: tengzhaowei2002@163.com *

Contributed equally

112020

15092020

22 5 4340 4350 31032020 31072020

2020

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Treatment-resistant schizophrenia (TRS) is a common phenotype of schizophrenia that places a considerable burden on patients as well as on society. TRS is known for its tendency to relapse and uncontrollable nature, with a poor response to antipsychotics other than clozapine. Therefore, it is urgent to identify objective biological markers, so as to guide its treatment and associated clinical work. In the present study, the peripheral blood mononuclear cells (PBMCs) of patients with TRS and a healthy control group, which were gender-, age- and ethnicity-matched, were subjected to microRNA (miRNA/miR) sequencing to screen out the top three miRNAs with the highest fold change values. These were then validated in the TRS (n=34) and healthy control (n=31) groups by reverse transcription-quantitative PCR. For two of the top three miRNAs, the PCR results were in accordance with the sequencing result (P<0.01), while the third miRNA exhibited the opposite trend (P<0.01). To elucidate the functions of these two miRNAs, Homo sapiens (hsa)-miR-218-5p and hsa-miR-1262 and their regulatory network, target gene prediction was first performed using online TargetScan and Diana-micro T software. Bioinformatics analysis was then performed using functional enrichment analysis to determine the Gene Ontology terms in the category biological process and the Kyoto Encyclopedia of Genes and Genomes pathways. It was revealed that these target genes were markedly associated with the nervous system and brain function, and it was obvious that the differentially expressed miRNAs most likely participated in the pathogenesis of TRS. A receiver operating characteristic curve was generated to confirm the distinct diagnostic value of these two miRNAs. It was concluded that aberrantly expressed miRNAs in PMBCs may be implicated in the pathogenesis of TRS and may serve as specific peripheral blood-based biomarkers for the early diagnosis of TRS.

treatment-resistant schizophrenia microRNA biomarker diagnosis

Introduction

Schizophrenia (SCZ) is a complex and severe psychiatric disorder that affects ~1% of the general population worldwide, placing a heavy burden on the patients, the patients' family and society (1). Although administration of antipsychotics remains the major therapeutic strategy for SCZ, ~ one-third of patients are insensitive to them (2). Patients who only have a minor or no response to at least two non-clozapine drugs at adequate doses and treatment cycles are diagnosed with treatment-resistant SCZ (TRS) according to the guidelines of the American Psychiatric Association (3). At present, frequent rehospitalization and repeated episodes of symptoms in patients with TRS increase the burden on patients and their families (4). It is currently thought that a combination of environmental and genetic factors is most likely to be the major cause of SCZ. However, the precise diagnosis of SCZ still only depends on clinical symptoms, without the help of sensitive and objective molecular markers. It is therefore urgent and necessary to identify sensitive biomarkers for the early and accurate diagnosis of SCZ, so that earlier and effective therapeutic intervention may be provided (5–7). Nowadays, strong evidence from psychiatric clinical studies has indicated that environmental factors may result in aberrant gene expression in the central nervous system (CNS), which may be the major cause of symptom onset (8,9). A study has indicated that ≥50% of changes in gene expression in the CNS are involved in the changes in the expression of the same genes in peripheral tissues (10). Since brain tissues, several of which have deeply located and complex structures with difficult access, are challenging to directly investigate, the present study attempted to examine and elucidate potential biomarkers for SCZ via a number of blood-derived gene expression analyses (10–14). By using this non-invasive method, effective peripheral blood biomarkers may be discovered that may be considerably helpful in promoting earlier TRS diagnosis, as well as enabling the concise assessment of treatment response and even precisely predicting prognosis.

Recently, increasing attention has been paid to the regulatory mechanisms of microRNAs (miRNAs/miRs) (15–17), which are considered to be in charge of gene expression monitoring through the inhibition of mRNA translation (18). Previous studies have indicated that changes in the expression of certain mRNAs reflect parallel changes in the CNS (19) and peripheral blood (20,21). The notion that certain miRNAs may be sensitive biomarkers for SCZ has gained much research interest (22–27), since alterations in their expression in peripheral blood mononuclear cells (PBMCs) may provide valuable clues for the earlier and more precise diagnosis and more effective treatment of SCZ (28). It is worth noting that abnormally expressed miRNAs have been detected in the PBMCs and plasma of patients with SCZ. In a recent study, Lai et al (28) indicated that 7 miRNAs, namely miR-31, miR-431, miR-433, miR-107, miR-134, miR-99b and miR-487b, were downregulated in the PBMCs derived from 112 patients with mental illness, as compared with those of a control group of 76 non-psychotic patients. Of note, the imprinted delta-like homologue 1 and the type III iodothyronine deiodinase (DLK1-DIO3) region of 14q32 from PBMCs was considered an SCZ-associated miRNA signature (28). In that study, a total of 7 differentially expressed miRNAs [Homo sapiens (hsa)-miR-34a, miR-449a, miR-564, miR-432, miR-548d, miR-572 and miR-652] were identified to be closely associated with emotional retardation, restricted language and other neurocognitive dysfunctions in 30 patients with SCZ, as compared with 30 controls (28). In another noteworthy study, specific miRNA30e in the plasma was identified as a biomarker of SCZ diagnosis from the expressed profiles of 9 SCZ-associated miRNAs in PBMCs by means of reverse transcription-quantitative PCR (RT-qPCR) (29). Subsequently, miR-130b and miR-193a-3p were identified as biomarkers for the diagnosis of SCZ by high-throughput sequencing and RT-qPCR (30). The early growth response protein-1 (EGR1)/miR-30a-5p/neurogenic differentiation 1 (NEUROD1) axis was proposed as a promising biomarker for the diagnosis and treatment monitoring of patients with SCZ with acute mental states based on RT-qPCR technology (31).

miRNA analysis may therefore be used to effectively elucidate the etiology of SCZ and identify potential differences among mental illnesses with similar but unclear clinical symptoms (32). The present study aimed to identify novel and specific biomarkers to earlier and more precisely diagnose patients with TRS and implement effective therapeutic strategies. Second-generation high-throughput sequencing technology and bioinformatics were used to screen out and analyze the top three abnormally expressed miRNAs in PBMCs from 6 patients with TRS and 6 healthy controls. Next, RT-qPCR was used to validate the significantly differentially expressed miRNAs in a larger sample of cases and healthy controls. The newly identified target genes of the two confirmed differentially expressed miRNAs were predicted with the online software TargetScan and Diana-micro T, and bioinformatics methods further indicated that these target genes were closely associated with functions of the nervous system. Finally, these two miRNAs were proposed as diagnostic biomarkers for TRS, thereby providing a novel, simple and convenient diagnostic tool for this condition.

Materials and methods <sec> <title>Screening and collection of cases

The subjects of the present study were divided into two groups, the case group and the control group. For the case group, a total of 34 patients with TRS (17 males and 17 females) were recruited from the Psychiatric inpatient ward of Honghe Second People's Hospital and Yuxi Second People's Hospital (Yunnan, China). The inclusion criteria were as follows: i) All patients met the mental illness criteria outlined in the Diagnostic and Statistical Manual Fifth Edition Criteria for SCZ (DSM-5; American Psychiatric Association; 2013) (33); ii) Complete clinical data and iii) Voluntary participation in the study. These 34 patients were screened by 2 independent experienced deputy chief physicians, according to the inclusion criteria. The exclusion criteria were conditions such as physical or cerebral organic disease, other types of mental illness, drug abuse prior to hospitalization, patients who received modified electroconvulsive therapy within 3 months from enrollment, and patients with a history of blood transfusion within the last month.

The control group included 31 healthy individuals with no family history of major mental illness recruited at the Physical Examination Center of the Sixth Affiliated Hospital of Kunming Medical University (Kunming, China). Of the 31 control individuals, 16 were males and 15 females, all of whom had a normal physical examination and psychological examination to exclude the presence of SCZ, none of whom had suffered any recent adverse life events. The age and sex ratio of patients and controls was 2:1.

All included subjects were older than 18 years. According to the requirements of the ethics committee, written informed consent for the case group was required not only from the patient but also one family member/guardian, while for the healthy control group, informed consent was provided by each subject themselves. The present study was approved by the ethics committees of Honghe Second People's Hospital (Honghe, Yunnan; approval no. 20190013) and The Sixth Affiliated Hospital of Kunming Medical University (Yuxi, Yunnan; approval no. 2020kmykdx6h58).

Extraction of PBMCs

Approximately 5 ml peripheral venous blood was collected by quick blood sampling from 65 cases and processed within 1 h. This procedure was carried out according to the manufacturer's protocol (Dakwei Biological Engineering). In brief, fresh whole blood was collected in EDTA-K2 anti-coagulation tubes and then diluted with an equal volume of isotonic solution (normal saline). Subsequently, 15 ml separation liquid was added to the centrifuge tube. The volume ratio of the separation liquid, anticoagulated undiluted whole blood and isotonic solution (physiological saline) was 1:1:1. Next, the sample was centrifuged at 700–800 × g for 20–30 min at room temperature. Between the plasma and separation layers was a thin and dense white film. The white film layer was carefully pipetted into another centrifuge tube. The horizontal rotor was then centrifuged at 250 × g at room temperature for 10 min and the supernatant was discarded. Washing was repeated 1–2 times. Finally, 1 ml TRIzol^® (Invitrogen; Thermo Fisher Scientific, Inc.) was added and the sample was mixed. A total of 1.5 ml was added to the labeled cryotube. The cryotube was tightly sealed with parafilm and stored at −80°C.

RNA extraction and miRNA sequencing

RNA samples were harvested from 6 patients with TRS (3 males, 3 females) and 6 healthy controls (3 males, 3 females). The patients in the experimental groups were age-matched. Samples were delivered to Guangzhou Gidio Biotechnology Co., Ltd. for the extraction of small RNA and the construction of small RNA libraries. In brief, RNA was extracted using the TRIzol method, followed by detection of RNA purity with a NanoDrop 2000 (Thermo Fisher Scientific, Inc.) and RNA integrity with an Agilent 2100 Bioanalyzer (Agilent Technologies GmbH). Small RNAs were obtained by gel electrophoresis from total RNA, and 5′ and 3′ ends were joined. RT-qPCR was performed to construct small RNA libraries. Subsequently, the constructed libraries were tested for quality and yield using an Agilent 2100 and ABI StepOnePlus Real-Time PCR System (Thermo Fisher Scientific, Inc.). An Illumina HiSeq2000 (Illumina, Inc.) was used for sequencing. All procedures were in strict accordance with the manufacturer's protocols (Gidio Biotechnology Co., Ltd.)

RT-qPCR

RT-qPCR was used to validate the top three aberrantly expressed miRNAs in 34 patients with TRS and 31 healthy controls. In brief, after total RNA was isolated from PBMCs, complementary (c)DNA was synthetized using the miRcute Plus miRNA First-Strand cDNA kit and miRcute Plus miRNA qPCR kit (SYBR Green) according to the manufacturer's protocol. RT-qPCR was performed in the Quantitative PCR Reaction Amplifier (CFX96 Touch; Bio-Rad Laboratories, Inc.), with the reaction system including 2X miRcute Plus miRNA Premix (SYBR and ROX), PCR forward primer and PCR reverse primer. The primers for three miRNAs are summarized in Table SI. Detailed experimental instruments and reagents are included in Tables SII and SIII. h-actin for miRNA was applied to normalize the results. The expression levels of miRNAs were calculated using the 2^−ΔΔCq method (34). All RT-qPCR validation procedures were performed by Wuhan Genecreate Biological Engineering Co., Ltd.

miRNA target gene prediction and bioinformatics analysis

TargetScan release 7.2 (http://www.targetscan.org/) (35) and Diana-microT v5.0 (http://diana.imis.athena-innovation.gr/DianaTools/index.php?r=microT_CDS/index) (36), which are two major and widely accepted open databases, were used to predict the target genes of the significantly differentially expressed miRNAs. Subsequently, the target genes obtained from the two databases were intersected and underwent further GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. GO annotation and KEGG enrichment analysis of the target genes were performed using the OmicShare tools (37), a free online platform for data analysis (http://www.omicshare.com/tools). The GO project covers three classical domains, namely biological process (BP), cellular component and molecular function. The KEGG is a database that integrates genome, chemistry, and system function information, and aims to reveal the genetic and chemical blueprint of biological phenomena. P<0.05 was considered to indicate a statistically significant result. The data were analyzed and obtained using the OmicShare cloud platform (https://www.omicshare.com/).

Statistical analysis

SPSS 22 (IBM Corp.) and GraphPad Prism 7.0 (GraphPad Software, Inc.) were used for data analysis. The χ² test for the qualitative variables and the t-test for the quantitative variables were respectively used for the demographic analysis of the case and control groups. The Mann-Whitney U test was used to compare the differential expression levels of miRNAs. Receiver operating characteristic (ROC) curve analysis was used to assess the area under the curve (AUC). P<0.05 was considered to indicate a statistically significant difference.

Results <sec> <title>Demographic data of the subjects

All subjects were Han Chinese. The mean age (mean ± standard deviation) of TRS cases and controls was 41.82±11.414 and 37.06±11.363 years, respectively. There were no significant differences in terms of age or sex between the cases and healthy controls (Table I).

miRNA expression profile analysis

Analysis of the miRNA sequencing data revealed that a total of 34 miRNAs were significantly differentially expressed in the 6 cases of TRS as compared with the 6 controls (|log2 (fold change)|>1; P<0.05). Among the six miRNAs with significantly differential expression, two were markedly downregulated and four were upregulated (Table II). A heat map was constructed to display the differential hierarchical clustering distribution between TRS and control samples (Fig. 1).

RT-qPCR validation

The three miRNAs (hsa-miR-1262, hsa-miR-1304-3p and hsa-miR-218-5p) with the highest log₂FC values (|log₂FC|>4) were verified in 34 TRS cases and 31 controls by RT-qPCR. Of these three miRNAs, hsa-miR-218-5p and hsa-miR-1262 were markedly upregulated according to the next-generation sequencing results, while hsa-miR-1304-3p was markedly downregulated. The RT-qPCR results indicated that the expression of hsa-miR-218-5p and hsa-miR-1262 in the TRS cases was markedly higher than that in the controls, which was consistent with the sequencing results (P<0.01). However, the expression of hsa-miR-1304-3p in the TRS cases was also markedly higher than that in the controls, which was inconsistent with the sequencing results (P<0.01; Fig. 2).

Prediction of miRNA target genes

By using two publicly available and widely accepted miRNA target gene prediction databases, TargetScan and Diana-micro T, the target genes of hsa-miR-218-5p and hsa-miR-1262 were predicted (Fig. 3A and B, respectively). For hsa-miR-218-5p, 187 target genes were predicted, while for hsa-miR-1262 had only 51 target genes were obtained. Next, a miRNA-gene co-expression network was generated to summarize the interaction between miRNAs and target genes (Fig. 4).

GO and pathway enrichment analysis

Target genes of these two miRNAs were enriched by GO annotation and KEGG analysis. The results for the top 20 GO terms enriched by the target genes of hsa-miR-218-5p are presented in Fig. 5A. In the GO category BP, the target genes of hsa-miR-218-5p were prominently involved in various physiological and pathophysiological processes, including nervous system development, generation of neurons, cell morphogenesis regarding neuron differentiation, synapse organization, neuron projection morphogenesis and neuron differentiation. BP terms associated with brain development and function are presented in Table III. KEGG pathway analysis for the target genes of hsa-miR-218-5p indicated a significant enrichment in several important pathways including the PI3K/Akt, dopaminergic synapse and AMPK signaling pathways, as well as the Hippo signaling pathway-multiple species (Fig. 5B).

In the category BP, the GO analysis for the target genes of hsa-miR-1262 revealed their close association with BP regulation, biological regulation, cellular process regulation and cellular protein modification process (Fig. 5C). BP terms closely associated with brain development and function are presented in Table IV. KEGG pathway analysis for the target genes of hsa-miR-1262 indicated an obvious enrichment in crucial pathways, such as nonhomologous end-joining dopaminergic synapse oocyte meiosis and adrenergic signaling in cardiomyocytes (Fig. 5D).

Diagnostic evaluation of hsa-miR-218-5p and hsa-miR-1262

ROC curve analysis was performed to demonstrate the potential of the expression levels of the 2 miRNAs confirmed by RT-qPCR as biomarkers for distinguishing patients with TRS from healthy controls. As displayed in Fig. 6, the AUC for hsa-miR-218-5p was 0.80 (95% CI, 0.7301–0.8604; P<0.01; Fig. 6A), the best cutoff value was 0.62 (specificity, 0.826; sensitivity, 0.794). The AUC for hsa-miR-1262 was 0.69 (95% CI, 0.6129–0.7669; P<0.01; Fig. 6B), the best cutoff value was 0.47 (specificity, 0.828; sensitivity, 0.637). The AUC for the combination of the 2 miRNAs was able to clearly distinguish patients with TRS from healthy controls with an AUC of 0.7501 (95% CI, 0.7005–0.7998; P<0.01; Fig. 6C), the best cutoff value was 0.51 (specificity, 0.801; sensitivity, 0.711).

Discussion

TRS accounts for 30% of all cases of SCZ (38). Due to its refractory and complex nature, TRS is a unique and severe phenotype of SCZ and places a burden on the patients' family and society (39). Therefore, precision medicine for this particular type of SCZ requires to be urgently implemented in the clinic. Although precision medicine is a novel healthcare discipline (40) that gives priority to the unique characteristics of specific patients (41), its use in psychiatry is still lagging compared with other medical fields, since it relies on neuroimaging and/or biological technology to provide a benefit (42,43). Several studies have made outstanding contributions to the development of precise psychiatric medical treatments. The ENIGMA-Diffusion Tensor Imaging working group confirmed that the normal variation of human white matter is genetically associated with mental, neurological and developmental disorders (44,45). The ENIGMA consortium performed a meta-analysis of 760 patients with SCZ and 957 healthy participants and revealed that the cognitive function of patients with SCZ is associated with the connectivity of the overall brain structure (46). Another study indicated that antipsychotic drugs resulted in an acute extrapyramidal response in the body by acting on two single nucleotide polymorphisms (SNPs; rs938112 and rs2987902) of two genes (LSMAP and ABL1) and destroying the binding sites of key transcription factors (47), LSAMP is an adhesion molecule of neuronal cells that participates in Neurogenesis, neurite outgrowth and plasticity (48), Abelson non-receptor tyrosine kinase (ABL1) is related to neurodegenerative diseases such as Parkinson's disease (48); the schizophrenia genomics platform is able to further elucidate the pathogenesis of SCZ (49) and provide therapeutic prospects for personalized therapy (50).

It is well known that the therapeutic effect of antipsychotic drugs is mainly exerted through the blocking of dopamine D2 receptor and/or the serotonin type 2A receptor to reduce dopaminergic neurotransmission in the limbic midbrain. However, even though the molecular targets of psychotropic drugs have remained to be fully elucidated, several psychotropic drugs used in the clinic still comprise the cornerstone of psychosis treatment. Furthermore, clozapine is the only evidence-based treatment for TRS (51,52), with 60–70% of patients reporting relief of their symptoms following this treatment (53); however, at present, there are no evidence-based pharmacotherapies for the remaining 30% of TRS patients who are not responsive to clozapine (54). Therefore, it is important to maximize the effectiveness of the drug through personalized treatment (55). The aim of pharmacogenomics (PGx) is to study how genes control the pharmacokinetics and pharmacodynamics of drugs (56). The concept may be used interchangeably with pharmacogenetics: A specific gene or group of genes may affect the patient's response to drugs. Ongoing randomized controlled trials have evaluated the clinical utility of PGx testing on antipsychotics in patients with SCZ to probe the association between genes and drugs on SCZ (trial registration nos. NCT02573168 and NCT02566057). Furthermore, the feasibility of genome-wide association studies, determination of polygenic risk scores and epigenetics of DNA methylation may render PGx more suitable for clinical research (57).

Based on the above concepts of precision medicine and PGx, even the application of epigenetics in SCZ research (58), combined with the fact that TRS is a highly heritable mental disorder and its heritability is markedly stronger than that of drug-responsive SCZ (59) the aim of the present study was to explore more sensitive treatment targets for TRS. Multiple studies have indicated that changes in cells and molecules in the brain may be relatively reflected by changes in the transcription in PBMCs (60–62), which may possess conjunct miRNA expression patterns with the brain (63). A previous study suggested that overexpressed miR-223 directly acts on mRNAs in plasma by the method of RT-qPCR and luciferase assay, revealing that it may be involved in the pathogenesis of SCZ (64). miR-30e, miR-181b, miR-34a, miR-346 and miR-7 in plasma may be useful non-invasive biomarkers for the diagnosis of SCZ, miR-132, miR-181b, miR-30e and miR-432 in plasma may be the improvement of symptoms in patients with SCZ, revealing potential indicators of their prognosis (65). Therefore, we can speculate that the abnormal expression levels of certain miRNAs in plasma may directly reflect the gene expression in the brain. Whether the expression of miRNA in PBMC is indirectly related to the expression of genes in the brain and the permeability of the blood-brain barrier of miRNAs need to be further studied. Fortunately, certain studies have confirmed the value of abnormally expressed miRNAs in the peripheral blood as biomarkers for SCZ (66,67). Therefore, this novel, simple and simultaneously much cheaper strategy was explored by miRNA sequencing data analysis in the MCs derived from the peripheral blood of patients with TRS. With the assistance of bioinformatics analysis and RT-qPCR for validation, two biomarkers, hsa-miR-1262 and hsa-miR-218-5p, were first identified for the diagnosis of TRS, which exhibited significant alterations in patients with TRS, as compared with healthy controls. As there are no reports of these 2 miRNAs being associated with SCZ or other psychiatric diseases, hsa-miR-218-5p and hsa-miR-1262 underwent target gene predictions to provide further insight into the function of these two miRNAs using not only the TargetScan and Diana-microT online databases, but also GO and KEGG pathway analysis along with bioinformatics methods. The functions of these two miRNAs were indicated to be closely associated with brain development and function, including processes such as nervous system development, neuronal generation, neurogenesis and synapse function, and they were also mainly involved in dopaminergic synapses, Hippo signaling pathway-multiple species and GABAergic synapses, as demonstrated by KEGG analysis. Of note, certain studies on TRS candidate genes have focused on the dopaminergic and serotoninergic systems (68–72). SCZ is a complex disease that may be of neurodevelopmental origin and is characterized by abnormal synaptic plasticity at the cortical and subcortical levels, abnormal regulation of neurotransmitters (including dopamine, glutamate and serotonin) and abnormal signaling (73–75). Therefore, the ROC curve analysis reflected improved ability to differentiate patients with TRS from healthy controls (AUC value, 0.7501; 95% CI, 0.7005–0.7998; specificity, 0.801; sensitivity, 0.711) in the present study, hsa-miR-218-5p and hsa-miR-1262, which have a critical role in the pathophysiological changes of TRS, were indicated to be promising sensitive biomarkers for TRS diagnosis, therapy and prediction.

The miRNA-gene network analysis of the present study indicated that there were three genes, namely, RNF165, CBX5, and CACUL1 targeted by the 2 miRNAs. Among these three genes, the identification of ring finger protein RNF165 was in accordance with previous studies, in which it was identified as a good candidate for SCZ diagnosis. In the study by Pickard et al (76), five key genes, including TTMA, TCF4, DLGAP1, ARKL1 and ARKL2, were confirmed as biomarkers for SCZ diagnosis. Along with the development of SCZ gene diagnosis, more biomarkers were gradually identified, including CBX5 (also termed HP1), which serves as a platform for other repressive proteins, such as DNA methyltransferases, and has been indicated to have a role in the functioning of dopamine neurons and the construction of memory as a member of heterochromatin in the brain. In addition, recently identified CACUL1, also known as CAC1 and C10orf46, is a novel gene of the cullin family that was reported to promote cell proliferation in a gastric cancer cell line (77), serves as a potential antitumor target in prostate cancer (78), promotes the proliferation of human lung cancer (79) and has an important role in protecting cells from Aβ toxicity and oxidative stress in the hippocampus of patients with Alzheimer's disease (80). The Pharmacogenomics Knowledgebase (PharmGKB) may provide candidate antipsychotic drug target genes for research (81). A drug-gene interaction study indicated that the rs6313 position of the serotonin receptor gene HTR2A has a regulatory mechanism involving risperidone and olanzapine in SCZ (82). Another study indicated that gene variation of GABRB1 may be significantly associated with the antipsychotic drug dose (83). A strong connection has been reported between genetic polymorphisms (in the HSPG2, CNR1, DPP6 and SLC18A2 genes) and antipsychotic drugs, and even the side effects of these drugs (84). Another study identified 12 SNPs from 143 genes of 79 Han Chinese patients with SCZ via HaloPlex technology, with 2/12 SNPs determined to be significantly associated with the response to olanzapine, using the MassARRAY platform in two independent cohorts (85). Other studies confirmed that the effectiveness of antipsychotic drugs is closely associated with the variation of drug target genes in patients with TRS (86–89). Although these three genes (RNF165, CBX5 and CACUL1) have not been reported in the PharmGKB to interact with medicines, the above information may be useful for future drug research.

In conclusion, in the present study, a simple, fast and effective strategy, mainly based on the detection of the expression of the top three miRNAs in peripheral blood MCs, was applied. hsa-miR-218-5p and hsa-miR-1262 were proposed as two major biomarkers for earlier and more precise diagnosis of TRS, shedding new light on the development of both gene-based diagnosis and treatment. However, the present study had certain limitations. First, the sample size of the study was small. Furthermore, limitations regarding the diagnostic accuracy may lead to misdiagnosis of certain samples. For example, the AUC value in ROC analysis was not very high, and there may be false positives and false negatives. Another limitation of this study is that the inclusion criteria were not designed to tell apart patients with normal schizophrenia that may respond to treatment from those with TRS. Finally, although the two target genes were confirmed in a relatively large cohort, there was a lack of functional verification, such as gene knockout experiments. In the future, studies with larger samples and a multi-center design are required to further verify the present results.

Supplementary Material

Supporting Data

Acknowledgements

Not applicable.

Funding

The current study was supported by The National Natural Science Fund (grant no. 81760253).

Availability of data and materials

The datasets used and/or analyzed in this study are available from the corresponding author on reasonable request.

Authors' contributions

YZ and ZWT designed and supervised the study. XY, YQZ and QL were responsible for analyzing the data, writing the manuscript and revising it, and designing the figures. ZJL, XM, ZXL, WY, ZF and WYZ contributed to the statistical analysis of the data and text correction. All authors read and approved the final version.

Ethics approval and consent to participate

Not applicable.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References 1

Kahn

Sommer

Murray

Meyer-Lindenberg

Weinberger

Cannon

O'Donovan

Correll

Kane

van Os

Insel

Schizophrenia

Nat Rev Dis Primers1150672015

10.1038/nrdp.2015.67

27189524

Gillespie

Samanaite

Mill

Egerton

MacCabe

Is treatment-resistant schizophrenia categorically distinct from treatment-responsive schizophrenia? A systematic review

BMC Psychiatry17122017

10.1186/s12888-016-1177-y

28086761

Lehman

Lieberman

Dixon

McGlashan

Miller

Perkins

Kreyenbuhl

American Psychiatric Association; Steering Committee on Practice Guidelines

Practice guideline for the treatment of patients with schizophrenia, second edition

Am J Psychiatry1611562004

15000267

Samara

Dold

Gianatsi

Nikolakopoulou

Helfer

Salanti

Leucht

Efficacy, acceptability, and tolerability of antipsychotics in treatment-resistant schizophrenia: A network meta-analysis

JAMA Psychiatry731992102016

10.1001/jamapsychiatry.2015.2955

26842482

Lakhan

Kramer

Schizophrenia genomics and proteomics: Are we any closer to biomarker discovery?

Behav Brain Funct522009

10.1186/1744-9081-5-2

19128481

Schwarz

Bahn

Biomarker discovery in psychiatric disorders

Electrophoresis29288428902008

18512679

Schwarz

Bahn

The utility of biomarker discovery approaches for the detection of disease mechanisms in psychiatric disorders

Br J Pharmacol153S133S1362008

10.1038/sj.bjp.0707658

18193070

Qureshi

Mehler

Emerging roles of non-coding RNAs in brain evolution, development, plasticity and disease

Nat Rev Neurosci135285412012

10.1038/nrn3234

22814587

Tsankova

Renthal

Kumar

Nestler

Epigenetic regulation in psychiatric disorders

Nat Rev Neurosci83553672007

10.1038/nrn2132

17453016

Sullivan

Fan

Perou

Evaluating the comparability of gene expression in blood and brain

Am J Med Genet B Neuropsychiatr Genet1412612682006

10.1002/ajmg.b.30272

Vawter

Ferran

Galke

Cooper

Bunney

Byerley

Microarray screening of lymphocyte gene expression differences in a multiplex schizophrenia pedigree

Schizophr Res6741522004

10.1016/S0920-9964(03)00151-8

14741323

Glatt

Everall

Kremen

Corbeil

Sásik

Khanlou

Han

Liew

Tsuang

Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia

Proc Natl Acad Sci USA10215533155382005

10.1073/pnas.0507666102

16223876

Kuzman

Medved

Terzic

Krainc

Genome-wide expression analysis of peripheral blood identifies candidate biomarkers for schizophrenia

J Psychiatr Res43107310772009

10.1016/j.jpsychires.2009.03.005

19358997

Maffioletti

Tardito

Gennarelli

Bocchio-Chiavetto

Micro spies from the brain to the periphery: New clues from studies on microRNAs in neuropsychiatric disorders

Front Cell Neurosci8752014

10.3389/fncel.2014.00075

24653674

Yan

Chen

Wang

Zhao

Zeng

Wei

Yang

Qian

Zhou

The effects and the mechanisms of autophagy on the cancer-associated fibroblasts in cancer

J Exp Clin Cancer Res381712019

10.1186/s13046-019-1172-5

31014370

Zhou

Chang

Wang

Liu

Dai

Zhang

Cai

Liu

Interactome Analyses implicated CAMK2A in the genetic predisposition and pharmacological mechanism of bipolar disorder

J Psychiatr Res1151651752019

10.1016/j.jpsychires.2019.05.024

31150948

Yan

Gong

Commentary: Lico A causes ER stress and apoptosis via up-regulating miR-144-3p in human lung cancer cell line H292

Biomed J413913922018

10.1016/j.bj.2018.11.002

30709582

Ambros

The functions of animal microRNAs

Nature4313503552004

10.1038/nature02871

15372042

Yan

Zeng

Dai

Chen

Wei

Gong

Circular RNAs: Clinical relevance in cancer

Oncotarget9144414602018

10.18632/oncotarget.22846

29416705

Perkins

Jeffries

Sullivan

Expanding the ‘central dogma’: The regulatory role of nonprotein coding genes and implications for the genetic liability to schizophrenia

Mol Psychiatry1069782005

10.1038/sj.mp.4001577

15381925

Cheng

Papp

Varlamova

Dziema

Russell

Curfman

Nakazawa

Shimizu

Okamura

Impey

Obrietan

MicroRNA modulation of circadian-clock period and entrainment

Neuron548138292007

10.1016/j.neuron.2007.05.017

17553428

Kocerha

Faghihi

Lopez-Toledano

Huang

Ramsey

Caron

Sales

Willoughby

Elmen

Hansen

MicroRNA-219 modulates NMDA receptor-mediated neurobehavioral dysfunction

Proc Natl Acad Sci USA106350735122009

10.1073/pnas.0805854106

19196972

Coyle

MicroRNAs suggest a new mechanism for altered brain gene expression in schizophrenia

Proc Natl Acad Sci USA106297529762009

10.1073/pnas.0813321106

19251669

Miller

Wahlestedt

MicroRNA dysregulation in psychiatric disease

Brain Res133889992010

10.1016/j.brainres.2010.03.035

20303342

Gladkevich

Kauffman

Korf

Lymphocytes as a neural probe: Potential for studying psychiatric disorders

Prog Neuropsychopharmacol Biol Psychiatry285595762004

10.1016/j.pnpbp.2004.01.009

15093964

Marques-Deak

Cizza

Sternberg

Brain-Immune interactions and disease susceptibility

Mol Psychiatry102392502005

10.1038/sj.mp.4001643

15685252

Marazziti

Catena Dell'osso

Baroni

Masala

Dell'Osso

Consoli

Giannaccini

Betti

Lucacchini

Alterations of the dopamine transporter in resting lymphocytes of patients with different psychotic disorders

Psychiatry Res17554572010

10.1016/j.psychres.2009.03.009

19910055

Lai

Hsieh

Chen

Wen

Huang

Hsiao

Liu

MicroRNA expression aberration as potential peripheral blood biomarkers for schizophrenia

PLoS One6e216352011

10.1371/journal.pone.0021635

21738743

Sun

Zhang

Song

Zhao

Fan

Zhong

Niu

Guo

Dai

Aberrant microRNA expression in peripheral plasma and mononuclear cells as specific blood-based biomarkers in schizophrenia patients

J Clin Neurosci225705742015

10.1016/j.jocn.2014.08.018

25487174

Wei

Yuan

Liu

Wang

Yang

Wang

Deng

Zhao

Shen

Detection of circulating miRNA levels in schizophrenia

Am J Psychiatry172114111472015

10.1176/appi.ajp.2015.14030273

26183697

Liu

Zhang

Shugart

Yang

Liu

Sun

Yang

Guo

Shi

The early growth response protein 1-miR-30a-5p-neurogenic differentiation factor 1 axis as a novel biomarker for schizophrenia diagnosis and treatment monitoring

Transl Psychiatry7e9982017

10.1038/tp.2016.268

28072411

Kalia

Silva

Biomarkers of psychiatric diseases: Current status and future prospects

Metabolism64S11S152015

10.1016/j.metabol.2014.10.026

25467847

Tandon

Gaebel

Barch

Bustillo

Gur

Heckers

Malaspina

Owen

Schultz

Tsuang

Definition and description of schizophrenia in the DSM-5

Schizophr Res1503102013

10.1016/j.schres.2013.05.028

23800613

Rao

Huang

Zhou

Lin

An improvement of the 2^(-delta delta CT) method for quantitative real-time polymerase chain reaction data analysis

Biostat Bioinforma Biomath371852013

25558171

Agarwal

Bell

Nam

Bartel

Predicting effective microRNA target sites in mammalian mRNAs

ELife4e050052015

10.7554/eLife.05005

Paraskevopoulou

Georgakilas

Kostoulas

Vlachos

Vergoulis

Reczko

Filippidis

Dalamagas

Hatzigeorgiou

DIANA-microT web server v5.0: Service integration into miRNA functional analysis workflows

Nucleic Acids Res41W169W1732013

10.1093/nar/gkt393

23680784

Wang

Guo

Liu

Ruan

Wang

Competing endogenous RNA (ceRNA) hypothetic model based on comprehensive analysis of long non-coding RNA expression in lung adenocarcinoma

Peer J7e80242019

10.7717/peerj.8024

31720124

Meltzer

Treatment-resistant schizophrenia-the role of clozapine

Curr Med Res Opin141201997

10.1185/03007999709113338

9524789

Wimberley

Gasse

Meier

Agerbo

MacCabe

Horsdal

Polygenic risk score for schizophrenia and treatment-resistant schizophrenia

Schizophr Bull43106410692017

10.1093/schbul/sbx007

28184875

Ashley

Towards precision medicine

Nat Rev Genet175075222016

10.1038/nrg.2016.86

27528417

Tonelli

Shirts

Knowledge for precision medicine: Mechanistic reasoning and methodological pluralism

JAMA318164916502017

10.1001/jama.2017.11914

29052713

Alda

Personalized psychiatry: Many questions, fewer answers

J Psychiatry Neurosci383633652013

10.1503/jpn.130221

24148844

Manchia

Pisanu

Squassina

Carpiniello

Challenges and future prospects of precision medicine in psychiatry

Pharmgenomics Pers Med131271402020

32425581

Kochunov

Hong

Dennis

Morey

Tate

Wilde

Logue

Kelly

Donohoe

Favre

ENIGMA-DTI: Translating reproducible white matter deficits into personalized vulnerability metrics in cross-diagnostic psychiatric research

Human Brain Mapp16249982020

10.1002/hbm.24998

Thompson

Jahanshad

Ching

Salminen

Thomopoulos

Bright

Baune

Bertolín

Bralten

Bruin

ENIGMA and global neuroscience: A decade of large-scale studies of the brain in health and disease across more than 40 countries

Transl Psychiatry101002020

10.1038/s41398-020-0705-1

32198361

Holleran

Kelly

Alloza

Agartz

Andreassen

Arango

Banaj

Calhoun

Cannon

Carr

The relationship between white matter microstructure and general cognitive ability in patients with schizophrenia and healthy participants in the ENIGMA consortium

Am J Psychiatry1775375472020

10.1176/appi.ajp.2019.19030225

32212855

Boloc

Rodríguez

Torres

García-Cerro

Parellada

Saiz-Ruiz

Cuesta

Bernardo

Gassó

Lafuente

Identifying key transcription factors for pharmacogenetic studies of antipsychotics induced extrapyramidal symptoms

Psychopharmacology (Berl)237215121592020

10.1007/s00213-020-05526-8

32382784

Heinla

Leidmaa

Kongi

Pennert

Innos

Nurk

Tekko

Singh

Vanaveski

Reimets

Gene expression patterns and environmental enrichment-induced effects in the hippocampi of mice suggest importance of lsamp in plasticity

Front Neurosci92052015

10.3389/fnins.2015.00205

26136648

Rees

Owen

Translating insights from neuropsychiatric genetics and genomics for precision psychiatry

Genome Med12432020

10.1186/s13073-020-00734-5

32349784

Birnbaum

Weinberger

Special article: Translational science update. Pharmacological implications of emerging schizophrenia genetics: Can the bridge from ‘Genomics’ to ‘Therapeutics’ be defined and traversed?

J Clin Psychopharmacol403233292020

10.1097/JCP.0000000000001215

32433256

Chakos

Lieberman

Hoffman

Bradford

Sheitman

Effectiveness of second-generation antipsychotics in patients with treatment-resistant schizophrenia: A review and meta-analysis of randomized trials

Am J Psychiatry1585185262001

10.1176/appi.ajp.158.4.518

11282684

Wahlbeck

Cheine

Essali

Adams

Evidence of clozapine's effectiveness in schizophrenia: A systematic review and meta-analysis of randomized trials

Am J Psychiatry1569909991999

10401441

Meltzer

Treatment of the neuroleptic-nonresponsive schizophrenic patient

Schizophr Bull185155421992

10.1093/schbul/18.3.515

1357741

Siskind

Lee

Ravindran

Zhang

Motamarri

Kisely

Augmentation strategies for clozapine refractory schizophrenia: A systematic review and meta-analysis

Aust N Z J Psychiatry527517672018

10.1177/0004867418772351

29732913

van Westrhenen

Aitchison

Ingelman-Sundberg

Jukić

Pharmacogenomics of antidepressant and antipsychotic treatment: How far have we got and where are we going?

Front Psychiatry11942020

10.3389/fpsyt.2020.00094

32226396

Relling

Evans

Pharmacogenomics in the clinic

Nature5263433502015

10.1038/nature15817

26469045

Zhang

Malhotra

Recent progress in pharmacogenomics of antipsychotic drug response

Curr Psychiatry Rep20242018

10.1007/s11920-018-0886-y

29589131

Feinberg

Fallin

Epigenetics at the crossroads of genes and the environment

JAMA314112911302015

10.1001/jama.2015.10414

26372577

Nucifora

Woznica

Lee

Cascella

Sawa

Treatment resistant schizophrenia: Clinical, biological, and therapeutic perspectives

Neurobiol Dis1311042572019

10.1016/j.nbd.2018.08.016

30170114

Fisar

Raboch

Depression, antidepressants, and peripheral blood components

Neuro Endocrinol Lett2917282008

18283265

Cattaneo

Sesta

Calabrese

Nielsen

Riva

Gennarelli

The expression of VGF is reduced in leukocytes of depressed patients and it is restored by effective antidepressant treatment

Neuropsychopharmacology35142314282010

10.1038/npp.2010.11

20164831

Chen

Sun

Niu

Kong

Fan

Zhong

Zhang

A preliminary analysis of microRNA-21 expression alteration after antipsychotic treatment in patients with schizophrenia

Psychiatry Res2443243322016

10.1016/j.psychres.2016.04.087

27512922

Liang

Ridzon

Wong

Chen

Characterization of microRNA expression profiles in normal human tissues

BMC Genomics81662007

10.1186/1471-2164-8-166

17565689

Zhao

Jinde

Koike

Tada

Satomura

Yoshikawa

Nishimura

Takizawa

Kinoshita

Sakakibara

Altered expression of microRNA-223 in the plasma of patients with first-episode schizophrenia and its possible relation to neuronal migration-related genes

Transl Psychiatry92892019

10.1038/s41398-019-0609-0

31712567

Sun

Zhang

Niu

Guo

Song

Fan

Zhao

Zhong

Dai

A preliminary analysis of microRNA as potential clinical biomarker for schizophrenia

Am J Med Genet B Neuropsychiatr Genet1681701782015

10.1002/ajmg.b.32292

Guo

Shi

miRNAs of peripheral blood as the biomarker of schizophrenia

Hereditas15592018

10.1186/s41065-017-0044-2

28860957

Liu

Zhang

Wang

Shugart

Zhao

Liu

Sun

Yang

Zhang

Diagnostic value of blood-derived microRNAs for schizophrenia: Results of a meta-analysis and validation

Sci Rep7153282017

10.1038/s41598-017-15751-5

29127368

Inada

Nakamura

Iijima

Relationship between catechol-O-methyltransferase polymorphism and treatment-resistant schizophrenia

Am J Med Genet B Neuropsychiatr Genet12035392003

10.1002/ajmg.b.20023

Takahashi

Saito

Ishihara

Maeno

Inada

Ozaki

Relationship between three serotonin receptor subtypes (HTR3A, HTR2A and HTR4) and treatment-resistant schizophrenia in the Japanese population

Neurosci Lett43595982008

10.1016/j.neulet.2008.01.083

18359159

Kohlrausch

Gama

Lobato

Belmonte-de-Abreu

Callegari-Jacques

Gesteira

Barros

Carracedo

Hutz

Naturalistic pharmacogenetic study of treatment resistance to typical neuroleptics in European-Brazilian schizophrenics

Pharmacogenet Genomics185996092008

10.1097/FPC.0b013e328301a763

18551040

Bilic

Jukic

Vilibic

Savic

Bozina

Treatment- Resistant schizophrenia and DAT and SERT polymorphisms

Gene5431251322014

10.1016/j.gene.2014.03.050

24680725

Terzić

Kastelic

Dolžan

Plesničar

Influence of 5-HT1A and 5-HTTLPR genetic variants on the schizophrenia symptoms and occurrence of treatment-resistant schizophrenia

Neuropsychiatr Dis Treat114534592015

25759587

Catts

Fung

Long

Joshi

Vercammen

Allen

Fillman

Rothmond

Sinclair

Tiwari

Rethinking schizophrenia in the context of normal neurodevelopment

Front Cell Neurosci7602013

10.3389/fncel.2013.00060

23720610

Alexander-Bloch

Reiss

Rapoport

McAdams

Giedd

Bullmore

Gogtay

Abnormal cortical growth in schizophrenia targets normative modules of synchronized development

Biol Psychiatry764384462014

10.1016/j.biopsych.2014.02.010

24690112

Horváth

Mirnics

Schizophrenia as a disorder of molecular pathways

Biol Psychiatry7722282015

10.1016/j.biopsych.2014.01.001

24507510

Pickard

Malloy

Clark

Lehellard

Ewald

Mors

Porteous

Blackwood

Muir

Candidate psychiatric illness genes identified in patients with pericentric inversions of chromosome 18

Psychiatr Genet1537442005

10.1097/00041444-200503000-00007

15722956

Zheng

Zhao

Kong

Nan

Yao

Liao

CDK-Associated cullin 1 can promote cell proliferation and inhibit cisplatin-induced apoptosis in the AGS gastric cancer cell line

World J Surg Oncol1152013

10.1186/1477-7819-11-5

23311997

Choi

Lee

Kim

CACUL1 functions as a negative regulator of androgen receptor in prostate cancer cells

Cancer Lett3763603662016

10.1016/j.canlet.2016.04.019

27085459

Chen

Gao

Yang

Thakur

Ren

Zhang

Wang

Chen

CDK-Associated cullin 1 promotes cell proliferation with activation of ERK1/2 in human lung cancer A549 cells

Biochem Biophys Res Commun4371081132013

10.1016/j.bbrc.2013.06.048

23806693

Kong

Bai

Sun

Nan

Expression of the newly identified gene CAC1 in the hippocampus of alzheimer's disease patients

J Mol Neurosci472072182012

10.1007/s12031-012-9717-5

22415352

Drögemöller

Wright

Niehaus

Emsley

Warnich

Next-generation sequencing of pharmacogenes: A critical analysis focusing on schizophrenia treatment

Pharmacogenet Genomics236666742013

10.1097/FPC.0000000000000006

24141736

Maffioletti

Valsecchi

Minelli

Magri

Bonvicini

Barlati

Sacchetti

Vita

Gennarelli

Association study between HTR2A rs6313 polymorphism and early response to risperidone and olanzapine in schizophrenia patients

Drug Dev Res27doi:102020

Hettige

de Morae

Kennedy

De Luca

Candidate gene analysis of pharmacodynamic targets for antipsychotic dosage

Pharmacogenomics171992082016

10.2217/pgs.15.171

26780210

Zai

Tiwari

Zai

Maes

Kennedy

New findings in pharmacogenetics of schizophrenia

Curr Opin Psychiatry312002122018

10.1097/YCO.0000000000000417

29528898

Zhou

Sheng

Chen

Shen

Huai

Cui

Qin

Genetic association of olanzapine treatment response in Han Chinese schizophrenia patients

Front Pharmacol101772019

10.3389/fphar.2019.00177

30886581

Ruderfer

Charney

Readhead

Kidd

Kähler

Kenny

Keiser

Moran

Hultman

Scott

Polygenic overlap between schizophrenia risk and antipsychotic response: A genomic medicine approach

Lancet Psychiatry33503572016

10.1016/S2215-0366(15)00553-2

26915512

Gressier

Porcelli

Calati

Serretti

Pharmacogenetics of clozapine response and induced weight gain: A comprehensive review and meta-analysis

Eur Neuropsychopharmacol261631852016

10.1016/j.euroneuro.2015.12.035

26792444

Zhang

Malhotra

Pharmacogenetics and antipsychotics: Therapeutic efficacy and side effects prediction

Exp Opin Drug Metab Toxicol79372011

10.1517/17425255.2011.532787

Sriretnakumar

Huang

Müller

Pharmacogenetics of clozapine treatment response and side-effects in schizophrenia: An update

Exp Opin Drug Metab Toxicol11170917312015

10.1517/17425255.2015.1075003

Figure 1.

Heat map of 34 differentially expressed miRNAs between 6 cases of TRS and 6 controls. The columns display miRNA species and the rows individual samples. The color labels of each cell represent the expression level of miRNAs. Red represents upregulation and blue downregulation. Con, control; TRS, treatment-resistant schizophrenia; miRNA/miR, microRNA; hsa, Homo sapiens.

Figure 2.

Reverse transcription-quantitative PCR was performed to verify the expression of miRNAs from peripheral blood mononuclear cells in a large sample (TRS, n=34; control, n=31). miRNA levels are normalized to H-actin. GraphPad Prism 7 was used to plot the data and the Mann-Whitney U-test was performed and analyze significant differences. Data are presented as the mean ± SD. TRS, treatment-resistant schizophrenia; miRNA/miR, microRNA; hsa, Homo sapiens.

Figure 3.

Venn diagram of genes from three miRNA target prediction datasets. Target genes of (A) hsa-miR-218-5p and (B) hsa-miR-1262. miRNA/miR, microRNA; hsa, Homo sapiens.

Figure 4.

miRNA-gene co-expression network. Red represents the miRNA and blue represents the targets. miRNA/miR, microRNA; hsa, Homo sapiens.

Figure 5.

GO annotation and KEGG analysis of target genes of hsa-miR-218-5p and hsa-miR-1262. (A) Top 20 GO terms in the categories BB, CC and MF of the target genes of hsa-miR-218-5p (P<0.05). (B) A total of 14 statistically significantly enriched KEGG pathways were identified (P<0.05). (C) Top 20 GO terms (BB, CC and MF) of the target genes of hsa-miR-1262 (P<0.05). (D) A total of four statistically significantly enriched KEGG pathways are presented (P<0.05). GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; BP, Biological Process; CC, Cellular Component; MF, Molecular Function; miRNA/miR, microRNA; hsa, Homo sapiens.

Figure 6.

Diagnostic evaluation of 2 miRNAs. ROC curves of (A) hsa-miR-218-5P and (B) hsa-miR-1262, and (C) their combined ROC curve. ROC, receiver operating characteristic; AUC, area under the curve; miRNA/miR, microRNA; hsa, Homo sapiens.

Table I.

Demographic data of the cases and healthy controls.

Item	Cases (n=34)	Controls (n=31)	P-value^a
Sex			0.897
Male	17 (50.0)	16 (51.6)
Female	17 (50.0)	15 (48.4)
Age (years)	41.82±11.414	37.06±11.363	0.956

Values are expressed as n (%) or the mean ± standard deviation.

A t-test was performed for quantitative variables and a χ² test was performed for categorical variables.

Table II.

Differentially expressed miRNAs in peripheral blood mononuclear cells from cases vs. healthy controls.

miRNA	Log2FC	P-value	Direction of change
hsa-miR-1185-1-3p	1.38018786	0.014888805	Up
hsa-miR-1185-5p	2.055813817	0.01030376	Up
hsa-miR-1197	1.669823494	0.003827819	Up
hsa-miR-1262	4.449671416	0.01862468	Up
hsa-miR-127-5p	1.663957827	0.010180048	Up
hsa-miR-1304-3p	−7.019942788	0.004047489	Down
hsa-miR-136-3p	1.536408035	0.010754284	Up
hsa-miR-154-3p	2.635532322	0.001485503	Up
hsa-miR-181b-3p	2.399752284	0.037210273	Up
hsa-miR-218-5p	4.289711531	0.046989353	Up
hsa-miR-3653-3p	−3.106775941	0.031498569	Down
hsa-miR-369-3p	1.219176104	0.035538469	Up
hsa-miR-370-3p	1.383232956	0.034880938	Up
hsa-miR-376a-3p	1.265825554	0.041422801	Up
hsa-miR-376c-3p	1.281059616	0.019343549	Up
hsa-miR-377-3p	1.68135408	0.004843865	Up
hsa-miR-379-3p	1.33354957	0.02582764	Up
hsa-miR-382-5p	1.164588341	0.042177088	Up
hsa-miR-409-3p	1.215682697	0.031732848	Up
hsa-miR-409-5p	1.33800193	0.02587475	Up
hsa-miR-410-3p	1.144646395	0.04797714	Up
hsa-miR-411-3p	1.188991694	0.049323611	Up
hsa-miR-431-3p	1.433665135	0.037350442	Up
hsa-miR-431-5p	1.703040789	0.003189292	Up
hsa-miR-432-5p	1.265157654	0.035781474	Up
hsa-miR-485-3p	1.118394154	0.046032191	Up
hsa-miR-490-3p	1.792169025	0.040858315	Up
hsa-miR-493-5p	1.499534768	0.008400669	Up
hsa-miR-539-3p	1.558249279	0.010618161	Up
hsa-miR-543	1.112920365	0.044806101	Up
hsa-miR-654-5p	1.433656533	0.014340331	Up
hsa-miR-655-3p	1.943235022	0.013931082	Up
hsa-miR-6810-3p	2.183414515	0.042813132	Up
hsa-miR-758-3p	1.339956277	0.01416857	Up

FC, fold change; miRNA/miR, microRNA; hsa, Homo sapiens.

Table III.

Top 20 enriched GO Biological Process terms associated with brain development and functions of hsa-miR-218-5p target genes.

GO ID	GO term	Number of enriched genes	P-value
GO:0007399	Nervous system development	46	3.79×10⁻⁷
GO:0048699	Generation of neurons	34	7.75×10⁻⁷
GO:0048667	Cell morphogenesis involved in neuron differentiation	19	1.15×10⁻⁶
GO:0050808	Synapse organization	16	1.34×10⁻⁶
GO:0048812	Neuron projection morphogenesis	20	1.40×10⁻⁶
GO:0030182	Neuron differentiation	31	1.73×10⁻⁶
GO:0031175	Neuron projection development	25	2.26×10⁻⁶
GO:0048666	Neuron development	27	2.52×10⁻⁶
GO:0099173	Postsynapse organization	10	2.78×10⁻⁶
GO:0022008	Neurogenesis	34	3.36×10⁻⁶
GO:0050767	Regulation of neurogenesis	22	5.84×10⁻⁶
GO:0045664	Regulation of neuron differentiation	19	7.09×10⁻⁶
GO:0050773	Regulation of dendrite development	9	8.09×10⁻⁶
GO:0048523	Negative regulation of cellular process	72	1.06×10⁻⁵
GO:0016358	Dendrite development	11	1.08×10⁻⁵
GO:0060996	Dendritic spine development	7	2.48×10⁻⁵
GO:0050807	Regulation of synapse organization	10	4.98×10⁻⁵
GO:0060998	Regulation of dendritic spine development	6	5.14×10⁻⁵
GO:0099504	Synaptic vesicle cycle	9	6.40×10⁻⁵
GO:0050803	Regulation of synapse structure or activity	10	6.76×10⁻⁵

GO, gene ontology; miR, microRNA; hsa, Homo sapiens.

Table IV.

Top 20 enriched GO Biological Process terms associated with brain development and functions of hsa-miR-1262 target genes.

GO ID	GO term	Number of enriched genes	P-value
GO:0045196	Establishment or maintenance of neuroblast polarity	1	0.006
GO:0045200	Establishment of neuroblast polarity	1	0.007
GO:0055059	Asymmetric neuroblast division	1	0.014
GO:0022008	Neurogenesis	9	0.014
GO:0048812	Neuron projection morphogenesis	5	0.018
GO:0031630	Regulation of synaptic vesicle fusion to presynaptic active zone membrane	1	0.018
GO:2000465	Regulation of glycogen (starch) synthase activity	1	0.021
GO:0007399	Nervous system development	11	0.024
GO:0007611	Learning or memory	3	0.025
GO:2001223	Negative regulation of neuron migration	1	0.027
GO:0048699	Generation of neurons	8	0.028
GO:0031175	Neuron projection development	6	0.029
GO:1990403	Embryonic brain development	1	0.030
GO:0036445	Neuronal stem cell division	1	0.032
GO:0055057	Neuroblast division	1	0.032
GO:0007613	Memory	2	0.033
GO:0007272	Ensheathment of neurons	2	0.038
GO:0050890	Cognition	3	0.039
GO:0016082	Synaptic vesicle priming	1	0.043
GO:0030182	Neuron differentiation	7	0.043

GO, gene ontology; miR, microRNA; hsa, Homo sapiens.