Qingqin S. Li
David Galbraith3†
Madhukar H. Trivedi- 1Neuroscience Therapeutic Area, Janssen Research and Development, LLC, Titusville, NJ, United States
- 2JRD Data Science, Janssen Research and Development, LLC, Titusville, NJ, United States
- 3Rancho Biosciences, San Diego, CA, United States
- 4Department of Psychiatry, Peter O’Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, United States
- 5Neuroscience Therapeutic Area, Janssen Research and Development, LLC, San Diego, CA, United States
Major depressive disorder (MDD) is an episodic condition with relapsing and remitting disease course. Elucidating biomarkers that can predict future relapse in individuals responding to an antidepressant treatment holds the potential to identify those patients who are prone to illness recurrence. The current study explored relationships between relapse risk in recurrent MDD and circulating microRNAs (miRNAs) that participate in RNA silencing and post-transcriptional regulation of gene expression. Serum samples were acquired from individuals with a history of recurrent MDD who were followed longitudinally in the observational study, OBSERVEMDD0001 (ClinicalTrials.gov Identifier: NCT02489305). Circulating miRNA data were obtained in 63 participants who relapsed (“relapsers”) and 154 participants who did not relapse (“non-relapsers”) during follow-up. The miRNA was quantified using the ID3EAL™ miRNA Discovery Platform from MiRXES measuring 575 circulating miRNAs using a patented qPCR technology and normalized with a standard curve from spike-in controls in each plate. The association between miRNAs and subsequent relapse was tested using a linear model, adjusting for age, gender, and plate. Four miRNAs were nominally associated with relapse status during the observational follow-up phase with a false discover rate adjusted p-value < 0.1. Enrichment analysis of experimentally validated targets revealed 112 significantly enriched pathways, including neurogenesis, response to cytokine, neurotrophin signaling, vascular endothelial growth factor signaling, relaxin signaling, and cellular senescence pathways. These data suggest these miRNAs putatively associated with relapse status may have the potential to regulate genes involved in multiple signaling pathways that have previously been associated with MDD. If shown to be significant in a larger, independent sample, these data may hold potential for developing a miRNA signature to identify patients likely to relapse, allowing for earlier intervention.
Introduction
Major depressive disorder (MDD) is an episodic condition with a relapsing and remitting disease course that has been linked to altered function and genetic variation in neural circuits and molecular pathways which regulate neuroimmune, neuroendocrine, and behavioral responses to stressors and opportunities for reward (1, 2). It has been reported that clinical features such as rumination, subthreshold inter-episodic depressive symptoms, family history, early MDD onset, and family history of mental illness predict relapse/recurrence (3, 4). The relationships between biological pathways implicated in MDD and episodes of illness relapse or recurrence remain inconclusive. A recent meta-analysis examining the role of cortisol in predicting future depressive episodes with a small effect size, however, sensitivity analysis revealed that the results became insignificant once applying outlier removal and/or exclusion of low-quality studies (5). Elucidating the neurobiological mechanisms which govern transitions in and out of illness episodes in MDD holds importance both for optimizing the clinical management of MDD and for guiding the development of novel treatments that more effectively prevent relapse and recurrence (6). Moreover, the identification of novel peripheral biomarker signatures that predict or correlate with relapse risk in MDD holds the potential to complement clinical assessments in precision medicine approaches to prognostic predictions.
Genetics plays an important role in the disease etiology of MDD and this has been characterized by large-scale consortium and cohort studies (7, 8). Gene-environment interactions appear to play a particularly important role in the pathogenesis of MDD, with early-life stress (ELS) constituting an important risk factor for many psychiatric disorders including MDD (9). ELS exerts its effect via epigenetic regulations and results in sustained hypothalamic-pituitary-adrenal (HPA)-axis dysregulation (10–12). Epigenetic regulations including deoxyribonucleic acid (DNA) methylation, histone modification, and non-coding RNA regulation play an important role in MDD (13–18). microRNA (miRNA) is the most common form of non-coding RNA and binds to the 3’-untranslated region (UTR) of its target mRNA and modulates the mRNA transcription (19). miRNA exists intracellularly but also in bodily fluids such as plasma, serum, and cerebrospinal fluid (20, 21). Extracellular exosome vesicles have been described in bodily fluid and are signaling messages mediating cell-cell communication (22, 23). Circulating RNAs including both the miRNA in extracellular vesicles including exosome vesicles and vesicle-free miRNA can be quantified easily without explicitly isolating the exosome fraction and emerge as a new class of biomarkers (24).
The roles of circulating miRNAs have been reported in other disease conditions such as coronary artery disease (CAD), diabetes, cancer, and autism spectrum disorder (25–28). In MDD, the roles of circulating miRNAs are also being characterized. Circulating hsa-let-7e-5p, hsa-miR-125a-5p, miR-34b-5p, miR-34c-5p, miR-451a, miR-17-5p, and miR-223-3p were up-regulated in plasma samples from patients with MDD and/or bipolar disorders compared to healthy controls (29–31), while miR-320a, miR-134 and miR-144-5p were down-regulated (29, 32, 33). The downregulation of miR-134 was also observed in the chronic unpredictable mild stress (CUMS) rat model of depression (33). The plasma miR-144-5p expression level was inversely related to the Montgomery-Asberg Depression Rating Scale (MADRS) depression score (32). Lower plasma levels of hsa-let-7b-5p are associated with a higher future risk of MDD (within 5 years) (34). let-7b-5p, let-7c-5p, miR-374b, and miR-10a were also downregulated in blood, and miR-508-3p and miR-152 were downregulated in the prefrontal cortex in patients with MDD compared to controls (35–37). Among these miRNAs, miR-34b-5p and miR-34c-5p expression levels were negatively correlated with NOTCH1 (30), while GRIN2A and DISC1 are the predicted targets for miR-320a and SLC17A7 is the predicted target for miR-451a, miR-17-5p, and miR-223-3p (29). There are quite some miRNAs implicated, and consistent players are yet to be revealed as more studies are conducted allowing for meta-analysis.
In this study, we used an MDD sample to study the relationship between circulating miRNA in serums and the future risk of relapse. Understanding biomarkers predicting future relapse episodes are useful in identifying the subjects with a worse disease course and has prognostic potential.
Materials and methods
Study cohort
A total of 217 blood serum samples from patients with MDD were collected from the observational clinical study, OBSERVEMDD0001 (ClinicalTrials.gov Identifier: NCT02489305), in which patients were prospectively followed for a maximum of 2.8 years for relapse detection. The entrance criteria included requirements that participants met the Diagnostic and Statistical Manual of Mental Disorders (DSM)-V criteria for non-psychotic, recurrent MDD, and had experienced the onset of the most recent major depressive episode within the 24 months before screening. In addition, patients must have been receiving an oral antidepressant pharmacotherapy (at an approved therapeutic dose) to which they had recently responded or remitted (within the past 3 months), as documented by a MADRS total score of ≤ 14 at both the screening and baseline visits. The serum samples used for the miRNA assays were drawn from the participants at either a baseline visit or follow-up visit. The mean sample storage duration prior to the assay is presented in Table 1. For both subjects who relapsed and those who did not relapse during the observational follow up phase, the miRNA assay was obtained for only a single time point; for the participants who ultimately relapsed, the biosample used for the assay was drawn before the relapse event (the mean time interval between sampling and relapse declaration appears in Table 1).
Table 1. Demographic and clinical characteristics of the samples.
A relapse event was defined by the occurrence of any of the following: (a) MADRS total score greater than or equal to 22 on at least two consecutive visits, with an interval of 1-2 weeks; (b) Hospitalization for worsening of depression; (c) Suicidal ideation with intent or suicidal behavior, and (d) as determined by the clinical investigator. In addition, a single occurrence of MADRS total score ≥ 22 not verified by a follow-up visit within two weeks was also considered a relapse if it was followed by a verification visit (within 23 days) and (i) CGI-S change-from-baseline = 2 or (ii) Medication change suggestive of worsening symptoms (increased dose or augmentation) within 14 days of MADRS = 22 (14-day window preceding or following MADRS = 22). The relapse events were confirmed by adjudication from clinicians.
Blood collection and serum processing
Patients were instructed to adhere to a low-fat diet on the day of sample collection. A volume of 8.5 ml blood was collected via venipuncture using a gold top serum separation tube (SST). The clinical sites were instructed to thoroughly mix the blood with the clotting activation agent by inverting the tube no less than five times and to allow the blood to clot for at least 30 min. The blood tube then was centrifuged at room temperature [minimum of 1,500 g for 15–20 min] until clot and serum were separated by a well-formed polymer barrier. The serum was aliquoted to four cryovials to minimize the freeze and thaw cycle and frozen.
Serum microRNA quantitative polymerase chain reaction assay
The absolute expression (copy numbers) of 575 candidate miRNAs were quantified in each patient and control biospecimen using ID3EAL miRNA Discovery Platform using a miRNA-specific RT-qPCR assay (MiRXES, Singapore) via a controlled workflow described in detail previously (38). Total RNA from 200 μl of patient serum specimen was isolated using miRNeasy serum/plasma miRNA isolation kit (Qiagen, Germany). Two sets of synthetic spike-in miRNA controls (three each at high, medium, and low concentration) were added to samples before RNA isolation and before RT-qPCR to monitor and normalize technical variations throughout the entire experiment. The absolute expression of each miRNA was normalized using calibrated spike-in controls and determined for each patient’s serum sample (expressed as log 2 copy number/ml serum). The samples were assayed across two plates in this study.
Data preprocessing and statistical analysis
For samples below the lowest level of detection (and therefore missing), it was set to a plate-level minimum where feasible. miRNA with greater than 90% missingness were excluded from the statistical analysis.
Although this is an absolute miRNA quantification assay, we used R package variancePartition (v1.20.0) (39) to understand the source of variation in the qPCR experiment. Factors such as age, gender, race, plate, sample age, and relapse status were considered (Supplementary Figure 1). The relapse status was analyzed both in a linear regression framework and in a cox regression survival analysis framework. For the binary endpoint, a linear model was fitted using R package limma (v3.46.0) (40) to identify miRNA predicting relapse status in the observational follow-up phase, while correcting for age, gender, and plate. miRNA associated with relapse status with a false discovery rate of less than 0.1 was reported. For the top miRNA identified, a Cox proportional hazards regression model using R package survival (v3.2-7) (41, 42) was also fitted to leverage the time to event information, and the miRNA level was also categorized as low, medium, and high expressed level for Kaplan–Meier curve plotting using R package survminer (v0.4.9) (43).
Over-representation analysis of putative targets of differentially expressed microRNAs
Ingenuity Knowledge Base (QIAGEN, Redwood City, CA, United States) is a database with curated and integrated miRNA targets from various sources and literature [miRecords (44), TarBase (45), TargetScan (46), and Ingenuity Expert Findings] and classified the targets into three groups: experimentally observed, predicted with high confidence [cumulative weighted context score (CWCS) less than -0.4 for TargetScan v7.2], and predicted with moderate confidence (CWCS between -0.2 and -0.4 for TargetScan v7.2). Over-representation analysis (ORA) (47) of experimentally observed targets from Ingenuity Knowledge Base and a few manually curated targets were performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (48) and R package clusterProfiler v3.18.1 (49). ORA was performed using additionally performed using the Gene Set Enrichment Analysis (GSEA) resource1 which assumes a broader background gene set and test overrepresentation at higher levels of the ontology hierarchy. Gene ontology databases used included c2.cp and c5 subsets of Molecular signatures database (MSigDB) (50).
Results
The demographic and clinical characteristics of the participant sample are described in Table 1. Of 217 participants, 63 (29.0%) relapsed. The average MADRS score at the time of blood draw for subjects who relapsed during the follow-up phase is 8.25 ± 5.17, and the relapse event occurred 16.9 ± 11.6 weeks later on average. The average MADRS score at the time of blood draw for subjects who were censored (i.e., did not relapse during the follow-up phase) is 5.47 ± 5.37, and the censor occurred 8.6 ± 9.1 weeks later on average.
Of the 575 miRNAs profiled with this technology, 287 miRNAs passing the 90% missingness filtering were included in the statistical analysis. Four miRNAs were nominally associated with relapse status during the follow-up phase (false discovery rate (FDR) adjusted p-value < 0.1, Figure 1 [volcano plot] and Table 2). Specifically, lower levels of hsa-miR-199b-5p (Figure 2A and Supplementary Figure S3A), hsa-miR-215-5p (Figures 2B, 3B), hsa-miR-200a-3p (Figures 2C, 3A), hsa-miR-143-3p miRNA (Figure 2D and Supplementary Figure S3B) were associated with risk of future relapse (Table 2). Overall, more miRNAs trended toward being down-regulated than up-regulation in the subgroup who relapsed versus the subgroup who did not relapse. The experiment used a balanced block design (Supplementary Figure 1 and the design table for the study) to minimize technical artifacts. Except a few miRNAs, none of the factors examined explained a significant proportion of the total variability of miRNA expression levels. For the four miRNAs identified, the variance fractions explained by plate or other factors did not dominate the total variance explained (Supplementary Figure 2). For example, for hsa-miR-199b-5p, 5% of the variance was explained by the risk of future relapse, while 0.074, 2.73, 3.32, 0.18, and 0.038% of the variances were explained by sex, age, race, sample age, and plate, respectively. Using the Cox proportional hazard regression framework, all four miRNAs were significant with p < 0.05 (Supplementary Table 1).
Figure 1. Volcano plot. logFC for relapser vs. non-relapser was plotted against -logP, where P is the association p-value.
Table 2. microRNA (miRNA) associated with relapse status with FDR <0.1.
Figure 2. Violin plots for the differentially expressed miRNA (FDR < 0.1) between patients who relapsed during the follow-up phase (1) vs. those who did not (0).
Figure 3. Kaplan–Meier curves for (A) hsa-miR-200a-3p, (B) hsa-miR-215-5p.
Ingenuity Knowledge Base revealed 44 unique experimentally validated targets of the four miRNAs identified and two of them, ZEB1 and ZEB2, were targets for both hsa-miR-215-5p and hsa-miR-200a-3p (Supplementary Table 2). Enrichment analysis of these experimentally validated targets augmented by additional manual curation of the literature revealed 112 significantly enriched pathways from KEGG database using clusterProfiler (Supplementary Table 3 and Figures 4A–C), including the MAPK signaling (p = 6.18 × 10–6, q-value = 3.84 × 10–5), neurotrophin signaling (p = 1.66 × 10–5, q-value = 7.26 × 10–5), vascular endothelial growth factor (VEGF) signaling (p = 1.55 × 10–4, q-value = 4.22 × 10–4), relaxin signaling (p = 3.11 × 10–4, q-value = 6.80 × 10–4), and cellular senescence pathways (p = 6.52 × 10–6, q-value = 3.84 × 10–5). A full list of both experimentally validated targets and predicted targets from Ingenuity Knowledge Base is provided in Supplementary Table 4. ORA using GSEA additionally revealed gene sets such as neurogenesis (p = 3.33 × 10–13, FDR q-value = 1.91 × 10–10) and response to cytokine (p = 8.30 × 10–10, FDR q-value = 1.55 × 10–7; see Supplementary Table 5). Additional enriched pathways from canonical databases are provided in Supplementary Table 6 and these included the MAPK pathway from BIOCARTA,2 the MAPK signaling pathway from WIKI Pathways,3 and the MAPK signaling pathway from KEGG (48, 51) that was also identified by the clusterProfiler analysis using the same KEGG database.
Figure 4. KEGG gene set enrichment analysis results. Targets of miRNA were plotted against the enriched pathway (A). Selected enriched gene sets were also plotted with p-value (B) and gene ratio (C) on x-axis.
Discussion
In this study, we evaluated the relationship between circulating miRNA and the future relapse status and identified four miRNAs nominally associated with the relapse status with FDR less than 0.1. Enrichment analysis of experimentally validated targets of these four miRNAs revealed 112 statistically significant pathways enriched among the experimentally validated targets of these four miRNAs, including neurogenesis, response to cytokine, neurotrophin signaling, VEGF signaling, relaxin signaling, and cellular senescence pathways, which were implicated in MDD previously. Higher expression level of four miRNAs (hsa-miR-199b-5p, hsa-miR-143-3p, hsa-miR-200a-3p, and hsa-miR-215-5p) identified as being nominally associated with lower risk of future MDD relapse in this study. Hsa-miR-199b-5p was implicated in neurogenesis, MAPK signaling, and BDNF signaling pathways, all of which were statistically enriched among the experimentally validated targets of the four miRNAs identified in the current study. The hsa-miR-199b-5p expression level increases during early brain development (52). Compared to healthy control, hsa-miR-199b-5p was up-regulated in neural progenitor cells and neurons differentiated from iPSC of patients with Rett syndrome (RTT) carrying mutations in the methyl-CpG-binding protein 2 (MECP2) gene (52). MECP2 has been implicated in RTT, a neurodevelopment disorder, and other neuropsychiatric conditions including depression and cocaine abuse (53). iPSC derived from MECP2 mutants was defective in neurogenesis, and overexpression of hsa-miR-199b-5p resulted in a similar phenotype (52). In contrast, inhibition of hsa-miR-199b-5p could partially rescue the MECP2 deficiency (52). Thus hsa-miR-199b-5p was placed in a neurogenesis pathway downstream of MECP2 and its function was mediated via extracellular signal-regulated kinase (ERK/MAPK) (52). As such MAPK signaling pathway was enriched among the experimentally validated miRNA targets including mitogen-activated protein kinase 1 (MAPK1), mitogen-activated protein kinase 7 (MAPK7), mitogen-activated protein kinase 12 (MAPK12), mitogen-activated protein kinase kinase 4 (MAP2K4), protein kinase cAMP-activated catalytic subunit beta (PRKACB), KRAS proto-oncogene, GTPase (KRAS), transforming growth factor beta 2 (TGFB2), and Rac family small GTPase 1 (RAC1). MAP2K4 is a target of hsa-miR-200a-3p, while MAPK7 and MAPK12 are targets of hsa-miR-143-3p according to the curations from Ingenuity Knowledge Base.
MeCP2 protein contributes to ELS-dependent epigenetic programming as phosphorylation of MeCP2 leads to hypomethylation of corticotropin releasing hormone (CRH), arginine vasopressin (AVP), and proopiomelanocortin (POMC, encoding a prohormone for adrenocorticotropic hormone (ACTH), a key mediator of the adrenocortical response to stress), and upregulation of these transcripts and ultimately enhancement of HPA-axis activity (54–57).
MECP2 regulates the expression of the brain-derived neurotrophic factor (BDNF) gene and genes modulating neuronal physiology such as calcium/calmodulin-dependent kinase CAMK2D and the voltage-gated potassium channel KCNH7, and genes involved in axon guidance and synapse formation such as EPHA7, SDK1 and CNTN4 (58–60). Mice expressing a truncated version of MeCP2 display anxiety-like phenotype and has abnormal stress response and elevated serum corticosterone levels, as the truncated form of MeCP2 fails to bind to the CpG rich promoter region of the CRH gene (54). In addition, mice carrying a genetic knock-in mutation that eliminates the phosphorylation site of MeCP2 also exhibit depressive-like behaviors and do not respond to chronic imipramine treatment (61). Decreased levels of MeCP2 and BDNF were detected in the hippocampus in a preclinical CUMS-induced rat model for depression and in the blood of depressed patients (62). Escitalopram improves the expression of MeCP2 in the CUMS depression model (63). Among the targets of differentially expressed miRNA, the BDNF signaling pathway was also enriched (Supplementary Table 6).
Hsa-miR-199b-5p is likely implicated in inflammation as it was reported to be correlated with absolute neutrophil count (64) and upregulated in tuberculosis patients (65). The immune pathways (Th1 and Th2 cell differentiation, IL-17 signaling pathway, tuberculosis, TNF signaling pathway, and T cell receptor signaling pathway in Supplementary Table 2) were also enriched in our study.
Hsa-miR-215-5p is another miRNA differentially expressed in the current study. It was up-regulated in the synaptosome in MDD in a human postmortem brain study where synaptic and total tissue fractions were obtained from the dorsolateral prefrontal cortex (dlPFC) of MDD and healthy patients (66). In our study lower expression levels of hsa-miR-199b-5p and hsa-miR-215-5p were associated with worse disease course (i.e., relapse in the follow-up phase).
MiR-200a-3p, another miRNA differentially expressed in the current study, was proposed to be a biomarker of central sensitization in chronic pain and depression. In unpredictable chronic mild stress (UCMS) and spared nerve injury (SNI) to illicit depressive-like and chronic pain behavior, miR-200a-3p was reduced 4 weeks after chronic stress in the prefrontal cortex (67). Pathways such as the glucocorticoid receptor (GR) signaling pathway and neurogenesis were positively regulated by miR-200a-3p after 4 weeks of stress, reflecting the activation of the HPA axis to overcome the stress (67). In our study, miR-200a-3p was downregulated in patients prone to relapse, consistent with the lower expression level when the rats exhibited depressive-like phenotype.
In rheumatoid arthritis (RA), increased expression of miR−143−3p is accompanied by promoting cell proliferation, inflammatory cytokine secretion, and inhibiting apoptosis (68), which is in contrast to a study using an Alzheimer’s disease (AD) cell model where miR-143-3p inhibition promotes neuronal survival (69). In the current study, a lower expression level of miR−143−3p was associated with future relapse. miR−143−3p inhibition also suppressed the Ras/p38 mitogen activated protein kinase (MAPK) signaling pathway, whereas TNF−α treatment stimulated it (68). TNF signaling pathway was also enriched in this study (Supplementary Table 2).
The study has several limitations which merit comment. Firstly, even though we report four miRNAs associated with future relapse status with a false discovery rate of less than 10%, we did not observe any miRNAs with FDR less than 5%. Secondly, the field has half a dozen studies on circulating miRNA for comparing MDD vs. healthy control, however, there is still no consensus on the definitive miRNAs associated with MDD, most likely owing to the effect size and that individual study did not have the power to unequivocally nominate the miRNA that can be replicated. Our study sample size is modest and will certainly benefit from future replication and/or meta-analysis studies examining the same endpoint. Third, the differences in miRNAs between the group who relapsed versus those that did not were based on a single time point, and the specific time course of the findings in relation to the course of depression was not established. Fourth, the differences associated with subsequent relapse have multiple possible interpretations. For example, they may reflect pathophysiological correlates which contribute to the process of impending relapse or compensatory changes in response to other factors which drove the relapse propensity.
In summary, we have identified four candidate miRNAs associated with risk of future risk and identified statistically significant pathways enriched among genes targeted by these miRNAs. Future experiments are needed to further shed light on the role and dysregulation of these candidate miRNAs in depression.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the respective local or central Institutional Review Boards (IRBs) overseeing the clinical sites participating in the study, these included the University of Pennsylvania Office of Regulatory Affairs IRB, University of Iowa IRB, Baylor College Of Medicine IRB, University of Michigan IRB, University of Cincinnati IRB, Sharp HealthCare IRB, Springfield Committee for Research Involving Human Subjects (SCRIHS), Western IRB, University of Kansas School of Medicine – Wichita Human Subjects Committee, Rush University Medical Center IRB, Hartford Hospital IRB, University of Massachusetts Medical School IRB, Baylor College Of Medicine IRB, Stanford University Administrative Panel on Human Subjects in Medical Research, Duke University Health System IRB, Johns Hopkins Medicine IRB, IRB-III – Medical University of South Carolina, Weill Cornell Medical College IRB, and Sterling IRB. The patients/participants provided their written informed consent to participate in this study.
Author contributions
RM, WD, and MT were involved in the conception and design of the OBSERVEMDD0001 clinical study. RM was involved in the OBSERVEMDD0001 study execution. QL was involved in the conception, design, and execution of this miRNA study and wrote the first draft of the manuscript. QL and DG were involved in the data analysis. All authors contributed to manuscript editing and revision.
Funding
This study received funding from Janssen Research & Development, LLC. MT received funding from the Agency for Healthcare Research and Quality (AHRQ), Cyberonics, Inc., the National Alliance for Research in Schizophrenia and Depression, and the National Institute of Mental Health and National Institute on Drug Abuse. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Acknowledgments
We thank the study participants without whom this study would not have been accomplished, and the investigational site staff, and Janssen cross-function study team (Global Clinical Operation, Data Management, Statistical Programming, and Statistics) for their contribution to this study. We also thank the staff at MiRXES for generating the assay data.
Conflict of Interest
QL and WD are employees of Janssen Research & Development, LLC, of Johnson & Johnson, and may hold equity in Johnson & Johnson. RM is an employee of Janssen Research & Development, LLC when the study was conducted; he is now retired and a consultant to Janssen Research & Development, LLC. MT is or has been an advisor/consultant and received fee from (lifetime disclosure): Abbott Laboratories, Inc., Abdi Ibrahim, Akzo (Organon Pharmaceuticals Inc.), Alkermes, AstraZeneca, Axon Advisors, Bristol-Myers Squibb Company, Cephalon, Inc., Cerecor, CME Institute of Physicians, Concert Pharmaceuticals, Inc., Eli Lilly and Company, Evotec, Fabre Kramer Pharmaceuticals, Inc., Forest Pharmaceuticals, GlaxoSmithKline, Janssen Global Services, LLC, Janssen Pharmaceutica Products, LP, Johnson & Johnson PRD, Libby, Lundbeck, Meade Johnson, MedAvante, Medtronic, Merck, Mitsubishi Tanabe Pharma Development America, Inc., Naurex, Neuronetics, Otsuka Pharmaceuticals, Pamlab, Parke-Davis Pharmaceuticals, Inc., Pfizer Inc., PgxHealth, Phoenix Marketing Solutions, Rexahn Pharmaceuticals, Ridge Diagnostics, Roche Products Ltd., Sepracor, SHIRE Development, Sierra, SK Life and Science, Sunovion, Takeda, Tal Medical/Puretech Venture, Targacept, Transcept, VantagePoint, Vivus, and Wyeth-Ayerst Laboratories.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.937360/full#supplementary-material
Footnotes
- ^ https://www.gsea-msigdb.org/gsea/msigdb/compute_overlaps.jsp
- ^ http://www.gsea-msigdb.org/gsea/msigdb/genesets.jsp?collection=CP:BIOCARTA
- ^ http://www.gsea-msigdb.org/gsea/msigdb/genesets.jsp?collection=CP:WIKIPATHWAYS
References
1. Price JL, Drevets WC. Neural circuits underlying the pathophysiology of mood disorders. Trends Cogn Sci. (2012) 16:61–71.
2. Drevets WC, Wittenberg GM, Bullmore ET, Manji HK. Immune targets for therapeutic development in depression: towards precision medicine. Nat Rev Drug Discov. (2022) 21:224–44. doi: 10.1038/s41573-021-00368-1
3. Kang SG, Cho SE. Neuroimaging biomarkers for predicting treatment response and recurrence of major depressive disorder. Int J Mol Sci. (2020) 21:2148.
4. Ronold EH, Joormann J, Hammar A. Facing recovery: emotional bias in working memory, rumination, relapse, and recurrence of major depression; an experimental paradigm conducted five years after first episode of major depression. Appl Neuropsychol Adult. (2020) 27:299–310. doi: 10.1080/23279095.2018.1550406
5. Kennis M, Gerritsen L, Van Dalen M, Williams A, Cuijpers P, Bockting C. Prospective biomarkers of major depressive disorder: a systematic review and meta-analysis. Mol Psychiatry. (2020) 25:321–38.
7. Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. (2018) 50:668–81. doi: 10.1038/s41588-018-0090-3
8. Howard DM, Adams MJ, Clarke TK, Hafferty JD, Gibson J, Shirali M, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci. (2019) 22:343–52. doi: 10.1038/s41593-018-0326-7
9. Zhao MZ, Song XS, Ma JS. Gene x environment interaction in major depressive disorder. World J Clin Cases. (2021) 9:9368–75. doi: 10.12998/wjcc.v9.i31.9368
10. Bockmuhl Y, Patchev AV, Madejska A, Hoffmann A, Sousa JC, Sousa N, et al. Methylation at the CpG island shore region upregulates Nr3c1 promoter activity after early-life stress. Epigenetics. (2015) 10:247–57. doi: 10.1080/15592294.2015.1017199
11. Van Bodegom M, Homberg JR, Henckens M. Modulation of the hypothalamic-pituitary-adrenal axis by early life stress exposure. Front Cell Neurosci. (2017) 11:87. doi: 10.3389/fncel.2017.00087
12. Silva RC, Maffioletti E, Gennarelli M, Baune BT, Minelli A. Biological correlates of early life stressful events in major depressive disorder. Psychoneuroendocrinology. (2021) 125:105103.
13. Dwivedi Y. Emerging role of microRNAs in major depressive disorder: diagnosis and therapeutic implications. Dialogues Clin Neurosci. (2014) 16:43–61.
14. Kichukova TM, Popov NT, Ivanov HY, Vachev TI. Circulating microRNAs as a novel class of potential diagnostic biomarkers in neuropsychiatric disorders. Folia Med. (2015) 57:159–72. doi: 10.1515/folmed-2015-0035
15. Dwivedi Y. Pathogenetic and therapeutic applications of microRNAs in major depressive disorder. Prog Neuropsychopharmacol Biol Psychiatry. (2016) 64:341–8. doi: 10.1016/j.pnpbp.2015.02.003
16. Yuan H, Mischoulon D, Fava M, Otto MW. Circulating microRNAs as biomarkers for depression: many candidates, few finalists. J Affect Disord. (2018) 233:68–78. doi: 10.1016/j.jad.2017.06.058
17. Roy B, Yoshino Y, Allen L, Prall K, Schell G, Dwivedi Y. Exploiting Circulating MicroRNAs as Biomarkers in Psychiatric Disorders. Mol Diagn Ther. (2020) 24:279–98.
18. Miao C, Chang J. The important roles of microRNAs in depression: new research progress and future prospects. J Mol Med. (2021) 99:619–36. doi: 10.1007/s00109-021-02052-8
20. Gui Y, Liu H, Zhang L, Lv W, Hu X. Altered microRNA profiles in cerebrospinal fluid exosome in Parkinson disease and Alzheimer disease. Oncotarget. (2015) 6:37043–53. doi: 10.18632/oncotarget.6158
22. Kourembanas S. Exosomes: vehicles of intercellular signaling, biomarkers, and vectors of cell therapy. Annu Rev Physiol. (2015) 77:13–27. doi: 10.1146/annurev-physiol-021014-071641
23. Simeone P, Bologna G, Lanuti P, Pierdomenico L, Guagnano MT, Pieragostino D, et al. Extracellular vesicles as signaling mediators and disease biomarkers across biological barriers. Int J Mol Sci. (2020) 21:2514.
24. Szilagyi M, Pos O, Marton E, Buglyo G, Soltesz B, Keseru J, et al. Circulating Cell-Free Nucleic Acids: main Characteristics and Clinical Application. Int J Mol Sci. (2020) 21:6827.
25. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. (2008) 105:10513–8.
26. Wang X, Sundquist J, Zoller B, Memon AA, Palmer K, Sundquist K, et al. Determination of 14 circulating microRNAs in Swedes and Iraqis with and without diabetes mellitus type 2. PLoS One. (2014) 9:e86792. doi: 10.1371/journal.pone.0086792
27. Bye A, Rosjo H, Nauman J, Silva GJ, Follestad T, Omland T, et al. Circulating microRNAs predict future fatal myocardial infarction in healthy individuals – The HUNT study. J Mol Cell Cardiol. (2016) 97:162–8. doi: 10.1016/j.yjmcc.2016.05.009
28. Salloum-Asfar S, Satheesh NJ, Abdulla SA. Circulating miRNAs, small but promising biomarkers for autism spectrum disorder. Front Mol Neurosci. (2019) 12:253. doi: 10.3389/fnmol.2019.00253
29. Camkurt MA, Acar S, Coskun S, Gunes M, Gunes S, Yilmaz MF, et al. Comparison of plasma MicroRNA levels in drug naive, first episode depressed patients and healthy controls. J Psychiatr Res. (2015) 69:67–71.
30. Sun N, Lei L, Wang Y, Yang C, Liu Z, Li X, et al. Preliminary comparison of plasma notch-associated microRNA-34b and -34c levels in drug naive, first episode depressed patients and healthy controls. J Affect Disord. (2016) 194:109–14. doi: 10.1016/j.jad.2016.01.017
31. Gecys D, Dambrauskiene K, Simonyte S, Patamsyte V, Vilkeviciute A, Musneckis A, et al. Circulating hsa-let-7e-5p and hsa-miR-125a-5p as possible biomarkers in the diagnosis of major depression and bipolar disorders. Dis Markers. (2022) 2022:3004338. doi: 10.1155/2022/3004338
32. Wang X, Sundquist K, Hedelius A, Palmer K, Memon AA, Sundquist J. Circulating microRNA-144-5p is associated with depressive disorders. Clin Epigenet. (2015) 7:69. doi: 10.1186/s13148-015-0099-8
33. Zhang HP, Liu XL, Chen JJ, Cheng K, Bai SJ, Zheng P, et al. Circulating microRNA 134 sheds light on the diagnosis of major depressive disorder. Transl Psychiatry. (2020) 10:95. doi: 10.1038/s41398-020-0773-2
34. Roumans S, Sundquist K, Memon AA, Hedelius A, Sundquist J, Wang X. Association of circulating let-7b-5p with major depressive disorder: a nested case-control study. BMC Psychiatry. (2021) 21:616. doi: 10.1186/s12888-021-03621-4
35. Smalheiser NR, Lugli G, Zhang H, Rizavi H, Cook EH, Dwivedi Y. Expression of microRNAs and other small RNAs in prefrontal cortex in schizophrenia, bipolar disorder and depressed subjects. PLoS One. (2014) 9:e86469. doi: 10.1371/journal.pone.0086469
36. Gururajan A, Naughton ME, Scott KA, O’connor RM, Moloney G, Clarke G, et al. MicroRNAs as biomarkers for major depression: a role for let-7b and let-7c. Transl Psychiatry. (2016) 6:e862. doi: 10.1038/tp.2016.131
37. Liu W, Zhang F, Zheng Y, He S, Zhang T, Guo Q, et al. The role of circulating blood microRNA-374 and microRNA-10 levels in the pathogenesis and therapeutic mechanisms of major depressive disorder. Neurosci Lett. (2021) 763:136184. doi: 10.1016/j.neulet.2021.136184
38. So JBY, Kapoor R, Zhu F, Koh C, Zhou L, Zou R, et al. Development and validation of a serum microRNA biomarker panel for detecting gastric cancer in a high-risk population. Gut. (2021) 70:829–37. doi: 10.1136/gutjnl-2020-322065
39. Hoffman GE, Schadt EE. variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformat. (2016) 17:483. doi: 10.1186/s12859-016-1323-z
40. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
41. Therneau TM, Grambsch PM. Modeling Survival Data: Extending the Cox Model. New York, NY: Springer (2000).
42. Therneau T. A Package for Survival Analysis in R. R package version 3.2-7. Vienna: R Core Team (2020).
43. Kassambara A, Kosinski M, Biecek P. survminer: Drawing Survival Curves using “ggplot2”. R package version 0.4.9. Vienna: R Core Team (2021).
44. Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. (2009) 37:D105–10. doi: 10.1093/nar/gkn851
45. Karagkouni D, Paraskevopoulou MD, Chatzopoulos S, Vlachos IS, Tastsoglou S, Kanellos I. DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. (2018) 46:D239–45. doi: 10.1093/nar/gkx1141
46. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. (2015) 4:e05005.
47. Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. (2004) 20:3710–5. doi: 10.1093/bioinformatics/bth456
48. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30.
49. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
50. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. (2011) 27:1739–40. doi: 10.1093/bioinformatics/btr260
51. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. (2017) 45:D353–61.
52. Mellios N, Feldman DA, Sheridan SD, Ip JPK, Kwok S, Amoah SK, et al. MeCP2-regulated miRNAs control early human neurogenesis through differential effects on ERK and AKT signaling. Mol Psychiatry. (2018) 23:1051–65. doi: 10.1038/mp.2017.86
53. Ausio J. MeCP2 and the enigmatic organization of brain chromatin. Implications for depression and cocaine addiction. Clin Epigenet. (2016) 8:58. doi: 10.1186/s13148-016-0214-5
54. Mcgill BE, Bundle SF, Yaylaoglu MB, Carson JP, Thaller C, Zoghbi HY. Enhanced anxiety and stress-induced corticosterone release are associated with increased Crh expression in a mouse model of Rett syndrome. Proc Natl Acad Sci USA. (2006) 103:18267–72. doi: 10.1073/pnas.0608702103
55. Murgatroyd C, Patchev AV, Wu Y, Micale V, Bockmuhl Y, Fischer D, et al. Dynamic DNA methylation programs persistent adverse effects of early-life stress. Nat Neurosci. (2009) 12:1559–66.
56. Wu Y, Patchev AV, Daniel G, Almeida OF, Spengler D. Early-life stress reduces DNA methylation of the Pomc gene in male mice. Endocrinology. (2014) 155:1751–62. doi: 10.1210/en.2013-1868
57. Zimmermann CA, Hoffmann A, Raabe F, Spengler D. Role of mecp2 in experience-dependent epigenetic programming. Genes. (2015) 6:60–86.
58. Chen WG, Chang Q, Lin Y, Meissner A, West AE, Griffith EC, et al. Derepression of BDNF transcription involves calcium-dependent phosphorylation of MeCP2. Science. (2003) 302:885–9. doi: 10.1126/science.1086446
59. Zhou Z, Hong EJ, Cohen S, Zhao WN, Ho HY, Schmidt L, et al. Brain-specific phosphorylation of MeCP2 regulates activity-dependent Bdnf transcription, dendritic growth, and spine maturation. Neuron. (2006) 52:255–69. doi: 10.1016/j.neuron.2006.09.037
60. Gabel HW, Kinde B, Stroud H, Gilbert CS, Harmin DA, Kastan NR, et al. Disruption of DNA-methylation-dependent long gene repression in Rett syndrome. Nature. (2015) 522:89–93. doi: 10.1038/nature14319
61. Hutchinson AN, Deng JV, Cohen S, West AE. Phosphorylation of MeCP2 at Ser421 contributes to chronic antidepressant action. J Neurosci. (2012) 32:14355–63. doi: 10.1523/JNEUROSCI.2156-12.2012
62. Su M, Hong J, Zhao Y, Liu S, Xue X. MeCP2 controls hippocampal brain-derived neurotrophic factor expression via homeostatic interactions with microRNA132 in rats with depression. Mol Med Rep. (2015) 12:5399–406. doi: 10.3892/mmr.2015.4104
63. Dionisie V, Ciobanu AM, Toma VA, Manea MC, Baldea I, Olteanu D, et al. Escitalopram targets oxidative stress, caspase-3, BDNF and MeCP2 in the hippocampus and frontal cortex of a rat model of depression induced by chronic unpredictable mild stress. Int J Mol Sci. (2021) 22:7483. doi: 10.3390/ijms22147483
64. Van Der Sijde F, Vietsch EE, Mustafa DAM, Li Y, Van Eijck CHJ. Serum miR-338-3p and miR-199b-5p are associated with the absolute neutrophil count in patients with resectable pancreatic cancer. Clin Chim Acta. (2020) 505:183–9. doi: 10.1016/j.cca.2020.03.005
65. Zhang Y, Zhang X, Zhao Z, Zheng Y, Xiao Z, Li F. Integrated bioinformatics analysis and validation revealed potential immune-regulatory miR-892b, miR-199b-5p and miR-582-5p as diagnostic biomarkers in active tuberculosis. Microb Pathog. (2019) 134:103563. doi: 10.1016/j.micpath.2019.103563
66. Yoshino Y, Roy B, Dwivedi Y. Differential and unique patterns of synaptic miRNA expression in dorsolateral prefrontal cortex of depressed subjects. Neuropsychopharmacology. (2021) 46:900–10. doi: 10.1038/s41386-020-00861-y
67. Satyanarayanan SK, Shih YH, Wen YR, Palani M, Lin YW, Su H, et al. miR-200a-3p modulates gene expression in comorbid pain and depression: molecular implication for central sensitization. Brain Behav Immun. (2019) 82:230–8. doi: 10.1016/j.bbi.2019.08.190
68. Yang Z, Wang J, Pan Z, Zhang Y. miR-143-3p regulates cell proliferation and apoptosis by targeting IGF1R and IGFBP5 and regulating the Ras/p38 MAPK signaling pathway in rheumatoid arthritis. Exp Ther Med. (2018) 15:3781–90. doi: 10.3892/etm.2018.5907
Keywords: circulating miRNA, depression, relapse, hsa-miR-199b-5p, hsa-miR-215-5p
Citation: Li QS, Galbraith D, Morrison RL, Trivedi MH and Drevets WC (2022) Circulating microRNA associated with future relapse status in major depressive disorder. Front. Psychiatry 13:937360. doi: 10.3389/fpsyt.2022.937360
Received: 10 May 2022; Accepted: 25 July 2022;
Published: 17 August 2022.
Edited by:
Mira Jakovcevski, Max Planck Institute of Psychiatry (MPI), GermanyReviewed by:
Daniele Caprioli, Sapienza University of Rome, ItalySajad Karampoor, Iran University of Medical Sciences, Iran
Copyright © 2022 Li, Galbraith, Morrison, Trivedi and Drevets. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingqin S. Li, qli2@its.jnj.com
†Present address: David Galbraith, Janssen Research and Development, LLC, Spring House, PA, United States