Back to Journals » International Journal of General Medicine » Volume 15
Integrated Analysis of miRNA and mRNA Expression Profiles Reveals the Molecular Mechanism of Posttraumatic Stress Disorder and Therapeutic Drugs
Authors Gan C, Jin Z, Hu G, Li Z, Yan M
Received 11 October 2021
Accepted for publication 28 January 2022
Published 8 March 2022 Volume 2022:15 Pages 2669—2680
DOI https://doi.org/10.2147/IJGM.S334877
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 2
Editor who approved publication: Dr Scott Fraser
Chunchun Gan,1 Zhan Jin,1 Gaobo Hu,1 Zheming Li,2 Minli Yan3
1Quzhou College of Technology, Quzhou, Zhejiang, 324000, People’s Republic of China; 2College of Pharmacy, Hangzhou Medical College, Hangzhou, Zhejiang, People’s Republic of China; 3Department of Neurology, The First Affiliated Hospital of Zhejiang Chinese Medical University (Zhejiang Provincial Hospital of Traditional Chinese Medicine), Hangzhou, 310000, Zhejiang, People’s Republic of China
Correspondence: Minli Yan, Department of Neurology, The First Affiliated Hospital of Zhejiang Chinese Medical University (Zhejiang Provincial Hospital of Traditional Chinese Medicine), Hangzhou, 310000, Zhejiang, People’s Republic of China, Tel +86-571-87077785, Email [email protected] Zheming Li, College of Pharmacy, Hangzhou Medical College, Hanzhou, People’s Republic of China, Email [email protected]
Purpose: Post-traumatic stress disorder (PTSD) is a result of trauma exposure and is related to psychological suffering as a long-lasting health issue. Further analysis of the networks and genes involved in PTSD are critical to the molecular mechanisms of PTSD.
Methods: In this study, we aimed to identify key genes and molecular interaction networks involved in the pathogenesis of PTSD by integrating mRNA and miRNA data.
Results: By integrating three high-throughput datasets, 5606 differentially expressed genes (DEGs) were detected, including five differentially expressed miRNAs (DEmiRNAs) and 5525 differentially expressed mRNAs (DEmRNAs). Nineteen upregulated and 46 downregulated DEmRNAs were identified in both GSE64813 and GSE89866 datasets, while five upregulated DEmiRNAs were found in the GSE87768 dataset. Functional annotations of these DEmRNAs indicated that they were mainly enriched in blood coagulation, cell adhesion, platelet activation, and extracellular matrix (ECM)-receptor interaction. Integrated protein-protein and miRNA-protein interaction networks among the DEGs were established with the help of 65 nodes and 121 interactions. Finally, 286 small molecules were obtained based on the Drug-Gene Interaction database (DGIdb). Three genes, prostaglandin-endoperoxide synthase 1 (PTGS1), beta-tubulin gene (TUBB1), and cyclin-dependent kinase inhibitor 1A (CDKN1A), were the most promising targets for PTSD therapy. Additionally, the present study also provided a higher performance diagnostic model for PTSD based on 17 DEmRNAs, which was validated in two independent datasets, GSE109409 and GSE63878.
Conclusion: Our data provides a new molecular aspect that ECM-receptor interaction and the platelet activation process could be the potential molecular mechanism of PTSD, and the genes involved in this process may be promising therapeutic targets. A higher-performance diagnostic model for PTSD has also been identified.
Keywords: differentially expressed genes, post-traumatic stress disorder, TCGA, prognosis, drug-gene interaction database
Introduction
Posttraumatic stress disorder (PTSD) often occurs after exposure to trauma. The PTSD patients often re-experience traumatic events in the form of nightmares, flashbacks, or intrusive memories. Patients also develop persistent avoidance of talking or thinking about the traumatic event, show negative cognitive and mood changes, as well as emotional and physical reactions.1 Around 20% of military members suffer from PTSD, which is also common among the general population, affecting approximately 7% of adult Americans.2 In China, the prevalence rate of PTSD was reported to be as high as 62.8% (257 of 430) one month after the 2008 Wenchuan earthquake, and the symptoms of 3 years of PTSD were also reported, ranging from 8.8% to 22.7%.3 In addition to the exacerbating psychological symptoms and negative effects on the quality of life, multiple studies have shown that patients with PTSD have an increased risk of hypertension and cardiovascular disease.4–6 Genetic and molecular studies on PTSD have discovered some potential genes, which have also been confirmed by other related studies. Studies using candidate-based approaches have focused more on biologically indicative genes. However, high-throughput technology has greatly facilitated the screening of potential genome-wide genes, which have already been widely used in PTSD research.7 The molecular mechanism for the development of PTSD is still unknown, but studies from multi-omics data have revealed that changes are mainly involved in the disruption of glucocorticoid signaling and inflammatory systems, and these alterations occur at the level of gene expression.8
Genome-wide transcriptomic studies provide enormous evidence of peripheral blood transcriptional markers for PTSD; however, many limitations still exist. Previous studies based on transcriptome data have shown the important role of genes such as ATP6APIL and DSCAM,9 FKBP5 and STAT5B10–12 and the biological pathways involved, as well as the subsequent cellular processes and immune-related pathway13 were significantly enriched in differentially expressed genes. Investigations at the system-level regarding the behaviors and properties of these genes revealed important roles of not only the modules of innate immune and interferon signaling transduction but also the modules of hemostasis and wound responsiveness.14
As the amount of publicly available high-throughput data in global databases continues to grow, an open question arises: How do we appropriately use this large-scale data to achieve a comprehensive understanding of disease at the molecular level? The Random Forest algorithm is a common machine learning algorithm often used as a routine bioinformatics protocol for disease and biology research and plays a key role in various areas of human life, especially it provides methods for the diagnosis and prognosis of human diseases.15–17 Several studies have used a variety of biomarkers to establish diagnostic or prognostic models for clinical patients.18,19 However, there are few diagnosis and prognostic analyses of PTSD.
In the current study, we integrated transcriptomic data of PTSD from public databases to identify differentially expressed genes at the transcriptome level. A miRNA-mRNA regulatory network was established based on the miRNA-mRNA interaction database. Drugs interacting with PTSD DEGs were further identified by searching the drug-gene interaction database, which may provide helpful information for the treatment of PTSD. Based on PTSD DEGs, we constructed a diagnostic model that showed good performance for the classification of PTSD and control samples. The design of this study is illustrated in Figure 1.
Materials and Methods
Data Preparation
The post-traumatic brain injury gene expression profiles and the corresponding clinical datasets were obtained from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). Dataset with next-generation high-throughput sequencing platform were only included to avoid the effect of different platforms on data analysis, resulting in a total of 182 PTSD patients and 151 control samples from five datasets. Raw read counts of the datasets GSE64813,14 GSE89866, GSE8776820 and GSE10940921 were downloaded and preprocessed before differential expression analysis. For GSE63878,14 the raw CEL data were downloaded and normalized by the robust multichip average (RMA) method. The criteria for selecting these datasets are as follows: 1) all blood samples; 2) all case-control data sets; 3) for differential expression analysis, the three datasets GSE64813, GSE89866, and GSE87768 are second-generation sequencing data and are the same platform.
Information on the platforms and number of samples of each dataset is provided in Table 1.
Table 1 Dataset Used in This Study |
Differential Analysis of Gene Expression
Datasets GSE64813 and GSE89866 performed on the same platform of Illumina high-throughput were used to screen differentially expressed mRNAs (DEmRNAs). Dataset GSE87768 was also performed on the same platform, including six case – control samples matched in age and gender, and were used to identify differentially expressed miRNAs (DEmiRNAs). DESeq2 package22 (v1.18.1) was applied to screen the DEGs (DEmRNAs and DEmiRNAs) using the raw read count data, according to the R code shown in the manual. miRNAs and mRNAs with a mean of reads < 1 were excluded. DEGs referred to miRNAs and mRNAs with an adjusted P<0.05, using Benjamini-Hochberg false discovery rate.
miRNA-mRNA Interaction Analysis
The miRNA targets and interactions between selected DEmiRNAs and DEmRNAs were obtained using the R package multMiR23 (v1.2.0) and the reference manual. We considered the genes detected as miRNA targets only in the validated database available in the package. STRING database24 was used to predict interactive proteins using DEmRNAs as a query. The PPI network with a combined score > 0.4 was considered to be a significant result. The biological network of the protein interactions and miRNA-mRNA regulatory network was visualized using the open source Cytoscape25 (v3.7.1).
Gene Ontology and Pathway Enrichment Analysis
Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway enrichment analyses were carried out using GeneCodis,26 a web-based tool that allows for the simultaneous evaluation of annotations from various sources to depict the potential biological functions of the DEmRNAs. For this purpose, prior to computing p-values, the a-priori algorithm was used to extract annotations frequently co-occurring in the analyzed genes.
Construction of PTSD Diagnostic Risk Model
The GSE64813 and GSE89866 datasets were merged because of their small number of samples. Batch effect was removed by using Bioconductor “combat”27 and no difference is seen in both PTSD and control samples between these 2 datasets after batch correction. The least absolute shrinkage and selection operator (LASSO) regression model, and random forest algorithm were used to further identify the important DEmRNAs with the ability to differentiate control samples and PTSD.28 For the random forest algorithm, we first set mtry = 1–50 and ntree = 500, to search for the best mtry parameter with the lowest error rate. The best ntree parameter is then selected by fixing the mtry parameter. Finally, all DEmRNs were ordered according to their importance, and only DEmRNs with cumulative importance > 95% were kept. For the LASSO regression model, 10-fold cross-validation was performed to determine the optimal lambda parameter, which corresponds to the lowest error rate. The overlapping DEmRNAs screened by the above two methods were used to establish a logistic diagnostic risk model for PTSD. The merged dataset was randomly grouped into a validation and training set with equal proportions. Additionally, in order to avoid the bias of a single random grouping, we repeated this process 100 times and calculated the mean AUC, specificity, and sensitivity. This research was conducted strictly following the statement of transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD).29
Drug Target Analysis
Based on searches of the Drug-Gene Interaction database (DGIdb), potential therapeutic targets of DEmRNAs were identified. The DGIdb30 database contains the interaction information of more than 10,000 drugs and 40,000 specific genes. Currently, the latest version of the DGIdb consists of more than 56,309 interaction claims collected from 30 sources. Such storage of information greatly expands on existing anti-neoplastic drug–gene interactions and the catalogue of druggable genes in the DGIdb.
Statistical Analysis
Statistical analyses in this study were conducted using the R software31 (version 3.6.1). Group comparisons were performed using the Mann-Whitney U-test for variables showing an abnormal distribution. For continuous variables, an independent t-test was conducted on normally distributed variables. The model specificity and sensitivity were analyzed with a receiver operating characteristic (ROC) curve and quantified by the area under the ROC curve32 (AUC) 27. All statistical analyses were two-sided. Statistical significance was set at p <0.05. We used “*” to stand for p < 1e-5, “**” for p < 0.01, and “***” for p < 0.05 in this study.
Results
Identification of DEmRNAs and miRNAs
The GSE64813 and GSE89866 datasets identified 669 and 4856 differentially expressed mRNAs in the PTSD case and control samples, respectively (Figure 2A and B). A total of 65 DEmRNAs were obtained in both datasets, which accounted for 9.72% and 1.17% of the total DEmRNAs, respectively (Figure 2C). Further analysis of the common DEmRNAs revealed that 19 genes showed a significant upregulation, and 46 genes were sharply downregulated in PTSD samples (Figure 2D and E), indicating that the activation and inhibition of some biological processes regulated by these DEmRNAs may be important molecular mechanisms for the formation of PTSD. The Genotype-Tissue Expression (GTEx) project allows users to access the transcriptomic data of healthy tissues collected from autopsies.33 The utility of GTEx databases revealed a different expression pattern of these overlapped DEmRNAs between healthy brain tissues and other healthy tissues (Figure 2F), which implied a high degree of similarity between gene expression in peripheral blood and that in brain tissue. The expression of five miRNAs expressions were significantly upregulated in the PTSD group compared to the control group (Figure 3A). By searching the multiMiR database, a total of 8348 miRNA-targets connections were obtained, including four DEmiRNAs and 10 DEmRNAs screened (Figure 3B).
GO Term and KEGG Pathway Enrichment Analyses on DEmRNAs
KEGG pathway and gene ontology analysis results indicated that DEmRNAs from GSE64813 were significantly enriched in the immune system process, innate immune response, cell adhesion, and defense response to virus (Figure 4A). It was found from the results of functional enrichment analysis that DEmRNAs from GSE89866 were particularly enriched in herpes simplex virus 1 infection, shigellosis, ribosome, endocytosis, and Salmonella infection (Figure 4B). Although the existence of difference between the function enrichments of two datasets DEmRNAs, interestingly, we found the GO term and KEGG pathway of the overlapping DEmRNAs were significantly enriched in Platelet activation, cell adhesion, and ECM-receptor interaction (Figure 4C). Also, we added functional enrichment analysis of 19 upregulated genes and 46 downregulated genes respectively. The analysis results (Tables S1–S4) showed that KEGG of downregulated genes also enriched in Platelet activation, cell adhesion, and ECM-receptor interaction pathways.
miRNA Interactome of PTSD
The integrated approach outlined in Figure 1 was applied using DEmRNAs and DEmiRNAs to generate a PTSD-specific miRNA interactome, with potential miRNA–mRNA target pairs. Eighteen prioritized miRNA–mRNA pairs, including four DEmiRNAs and ten DEmRNAs validated from the multiMiR database, were determined. The finding that 10 DE target genes interacted with four DEmiRNAs indicated that in a complex structure, miRNAs bound to target many genes. The network was composed of a large cluster of inter-connected miRNA–mRNA pairs in which the targets were downregulated, and miRNAs were upregulated, which contains 81 nodes and 175 edges (Figure 5A). Within the upregulated miRNA cluster, previously unknown PTSD-related CDKN1A showed two connected miRNAs. A total of 286 potential drugs interacting with 20 DEmRNAs were identified when drug-gene interactions were explored using DGIdb (Table S5). Most potential drugs might interact with prostaglandin-endoperoxide synthase 1 (PTGS1) (129/286), beta-tubulin gene (TUBB1) (50/286), and cyclin-dependent kinase inhibitor 1A (CDKN1A) (37/286) either in an unknown manner or as an inhibitor (Figure 5B, 310/337). Additionally, the drugs that interact with the CLIP2 gene were all antagonists and three up-regulated genes shared the same drugs that were all inhibitors (Figure 5B), indicating that these genes may be important potential therapeutic targets for PTSD.
DEmRNAs for Diagnostic Prediction of PTSD
After merging the two datasets (GSE64813 and GSE89866 datasets), the control and PTSD samples are randomly (seed=12345) divided into the same proportion of training set (control/PTSD: 38/38) and validation set (control/PTSD: 38/38). Batch correction analysis showed that batch effects were eliminated after batch correction when merging the GSE64813 and GSE89866 datasets (Figure S1A and B). The random forest (Figure 6A and B) and LASSO (Figure 6C) methods yielded 48 and 17 DEmRNAs, of which all the DEmRNAs were common to both methods. Using the logistic regression method, we further used these 17 DEmRNAs to construct a risk score model for PTSD diagnosis. The Violin diagram (Figure 6D) shows that the risk scores of the PTSD samples of the two datasets were significantly higher than those of the control samples. In addition, we also evaluated the ability of the diagnostic model’s risk score to classify PTSD and control samples, and found that on the training and validation sets, the diagnostic model’s risk score classified the two samples well, and the AUC values reached 0.867 and 0.845 (Figure 6E). Two independent datasets, GSE109409 and GSE63878, were then used to assess the model performance and showed an AUC of 0.767 and 0.795, respectively, indicating a good discrimination between PTSD and controls (Figure 6F).
Discussion
Differential expression analysis of two independent datasets was performed to identify pathways and genes involved in PTSD etiology. This study found a small number of differentially expressed genes in peripheral blood collected from individuals with and without PTSD in the two datasets. To understand the role of DEmRNAs in PTSD, we detected an over-representation of genes associated with ECM-receptor interaction and platelet activation pathways. These results are consistent with other studies,34–36 and show a robust relationship between the ECM-receptor interaction pathway and PTSD.
PTSD is a disorder that originates in the brain, despite the existence of a blood-brain barrier between the periphery and the brain, and may affect the gene expression patterns between them. However, many studies have confirmed that there is a good correspondence between gene expression levels in peripheral blood and brain tissue, as demonstrated by similar gene expression profiles observed in the postmortem brain.37 The gene expression profiles in peripheral blood can be used as effective biomarkers for many psychiatric diseases.38 Furthermore, given the persistence of psychosomatic symptoms of PTSD in many patients, gene expression in the peripheral blood is also associated with brain tissue expression.39 Finally, we observed similar expression characteristics of these differentially expressed genes identified in peripheral blood samples from normal brain tissues (Figure 2F). These results suggest that peripheral blood gene expression profiles could also reflect changes in gene expression in the brains of patients with PTSD.
From the current miRNA-mRNA interaction network analyses, one of the novel findings was the identification of 4 DEmiRNAs associated with 12 DEmRNAs. Interestingly, gene CDKN1A network analysis also identified modules correlated with wound response and hemostasis with significant overlap.40 We also identified several novel differentially expressed miRNAs in PTSD subjects and normal subjects, including hsa-miR-374b-5p, hsa-miR-629-5p, hsa-miR-101-3p, hsa-miR-20b-5p, and hsa-miR-485-3p. Three of the five DEmiRNAs, hsa-miR-20b-5p, hsa-miR-629-5p, and hsa-miR-101-3p, were also reported in previous studies.20,41 The top hit was hsa-miR-20b-5p, which has 14 validated targets (Figure 3B). This indicates that hsa-miR-20b-5p is altered and is likely to be a candidate marker for PTSD, a stress-related disorder. It has been shown that hsa-miR-20b-5p is increased in the hippocampus of Appswe/PSΔE 9 mice, and inhibition of miR-20b-5p attenuates apoptosis induced by Aβ25-35 in Alzheimer’s disease by targeting RhoC.42 The expression of hsa-miR-485-3p was found to be upregulated in mice models that showed susceptibility to chronic unpredictable mild stress (CUMS) compared to CUMS- resilience.43 In Cohen’s study, miR-485, in conjunction with synaptic vesicle protein SV2A, has been found to regulate dendritic spine number in an activity-dependent manner.44 The expression of hsa-miR-374b-5p, which was newly identified in this study, was found to be a poor prognostic factor in human glioma.45 These results suggested the reliability of our study.
In the current study, a diagnostic prediction model for PTSD was established using the screened DEmRNAs. The model risk score increased significantly from a control to PTSD, and a high AUC indicated good performance of the model in identifying patients with PTSD from individuals and healthy controls. This result may facilitate the development of a new diagnostic strategy from a molecular perspective. Nonetheless, this study also has the following limitations. First, more studies are warranted to confirm the robust performance of the model on a larger cohort, and the consistency between mRNAs and miRNAs in the circulation also needs to be established to determine whether the detection of DEmiRNAs from peripheral blood could serve as a novel tool for PTSD screening. Second, more efforts based on experimental data are required.
Conclusions
In summary, our data provide a new molecular aspect that ECM-receptor interaction and platelet activation biological process may be the potential molecular mechanism of PTSD, and the genes involved in this process may be promising therapeutic targets. The present study also provided a high-performance diagnostic model for PTSD. Specifically, DEmiRNAs and DEmRNAs associated with cell proliferation appear to be involved in PTSD. The data showed that modules related to wound healing and platelet activation may contribute to resilience against PTSD development. Future studies could confirm the role of platelets in the stress response. Due to limitations including small sample size and lack of epigenetic testing, in future studies, participants with and without PTSD will be studied more closely in relation to important potential confounding variables and multivariate data sets by performing a meta-analysis.
Ethical Statement
According to the guidelines of the Ethics Committee of The First Affiliated Hospital of Zhejiang Chinese Medical University (Zhejiang Provincial Hospital of Traditional Chinese Medicine), any research involving human body (Declaration of Helsinki) and animal experiments shall be subject to ethical review.
TCGA and GEO belong to public database. The patients involved in database have obtained ethical approval. This work did not include any experiments on humans or animals. Users can download relevant data for free for research and publish relevant articles, so there are no ethical issues and other conflicts of interest. The waived ethics approval was approved by Ethics Committee of The First Affiliated Hospital of Zhejiang Chinese Medical University (Zhejiang Provincial Hospital of Traditional Chinese Medicine).
Funding
This study was supported by the Medical Health Science and Technology Planning Project of Zhejiang Provincial Health Commission (2020KY665), Quzhou Science and Technology Project (2018K23) and Science and Technology Planning Project of Traditional Chinese Medicine of Zhejiang Provincial (2021ZB101).
Disclosure
The authors report no conflicts of interest in this work.
References
1. Sorrell JM. Diagnostic and statistical manual of mental disorders-5: implications for older adults and their families. J Psychosoc Nurs Ment Health Serv. 2013;51(3):19–22. doi:10.3928/02793695-20130207-01
2. Kessler RC, Berglund P, Demler O, Jin R, Merikangas KR, Walters EE. Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the national comorbidity survey replication. Arch Gen Psychiatry. 2005;62(6):593–602. doi:10.1001/archpsyc.62.6.593
3. Hong C, Efferth T. Systematic review on post-traumatic stress disorder among survivors of the Wenchuan earthquake. Trauma Violence Abuse. 2016;17(5):542–561. doi:10.1177/1524838015585313
4. Edmondson D, Cohen BE. Posttraumatic stress disorder and cardiovascular disease. Prog Cardiovasc Dis. 2013;55(6):548–556.
5. Bedi US, Arora R. Cardiovascular manifestations of posttraumatic stress disorder. J Natl Med Assoc. 2007;99(6):642–649.
6. Shah AJ, Lampert R, Goldberg J, Veledar E, Bremner JD, Vaccarino V. Posttraumatic stress disorder and impaired autonomic modulation in male twins. Biol Psychiatry. 2013;73(11):1103–1110.
7. Sumner JA, Duncan LE, Wolf EJ, et al. Letter to the Editor: posttraumatic stress disorder has genetic overlap with cardiometabolic traits. Psychol Med. 2017;47(11):2036–2039.
8. Girgenti MJ, Hare BD, Ghosal S, Duman RS. Molecular and cellular effects of traumatic stress: implications for PTSD. Curr Psychiatry Rep. 2017;19(11):85.
9. Logue MW, Smith AK, Baldwin C, et al. An analysis of gene expression in PTSD implicates genes involved in the glucocorticoid receptor pathway and neural responses to stress. Psychoneuroendocrinology. 2015;57:1–13.
10. Sarapas C, Cai G, Bierer LM, et al. Genetic markers for PTSD risk and resilience among survivors of the World Trade Center attacks. Dis Markers. 2011;30(2–3):101–110. doi:10.1155/2011/328054
11. Mehta D, Gonik M, Klengel T, et al. Using polymorphisms in FKBP5 to define biologically distinct subtypes of posttraumatic stress disorder: evidence from endocrine and gene expression studies. Arch Gen Psychiatry. 2011;68(9):901–910. doi:10.1001/archgenpsychiatry.2011.50
12. Yehuda R, Cai G, Golier JA, et al. Gene expression patterns associated with posttraumatic stress disorder following exposure to the World Trade Center attacks. Biol Psychiatry. 2009;66(7):708–711. doi:10.1016/j.biopsych.2009.02.034
13. Mehta D, Voisey J, Bruenig D, et al. Transcriptome analysis reveals novel genes and immune networks dysregulated in veterans with PTSD. Brain Behav Immun. 2018;74:133–142. doi:10.1016/j.bbi.2018.08.014
14. Breen MS, Maihofer AX, Glatt SJ, et al. Gene networks specific for innate immunity define post-traumatic stress disorder. Mol Psychiatry. 2015;20(12):1538–1545. doi:10.1038/mp.2015.9
15. Yang Y, Wang C, Wei N, et al. Identification of prognostic chromatin-remodeling genes in clear cell renal cell carcinoma. Aging. 2020;12(24):25614–25642. doi:10.18632/aging.104170
16. Zhao T, Wang D, Hu Y, Zhang N, Zang T, Wang Y. Identifying Alzheimer’s disease-related mirna based on semi-clustering. Curr Gene Ther. 2019;19(4):216–223. doi:10.2174/1566523219666190924113737
17. Huang GH, Cao XY, Li YY, et al. Gene expression profile of the hippocampus of rats subjected to traumatic brain injury. J Cell Biochem. 2019;120(9):15776–15789. doi:10.1002/jcb.28848
18. Shi M, Chen MS, Sekar K, Tan CK, Ooi LL, Hui KM. A blood-based three-gene signature for the non-invasive detection of early human hepatocellular carcinoma. Eur J Cancer. 2014;50(5):928–936. doi:10.1016/j.ejca.2013.11.026
19. Kim DJ, Cho EJ, Yu KS, et al. Comprehensive metabolomic search for biomarkers to differentiate early stage hepatocellular carcinoma from cirrhosis. Cancers. 2019;11(10):1497. doi:10.3390/cancers11101497
20. Martin CG, Kim H, Yun S, et al. Circulating miRNA associated with posttraumatic stress disorder in a cohort of military combat veterans. Psychiatry Res. 2017;251:261–265. doi:10.1016/j.psychres.2017.01.081
21. Boscarino C, Nalpathamkalam T, Pellecchia G, Li W, Thiruvahindrapuram B, Merico D. Using next-generation sequencing transcriptomics to determine markers of post-traumatic symptoms: preliminary findings from a post-deployment cohort of soldiers. G3. 2019;9(2):463–471. doi:10.1534/g3.118.200516
22. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8
23. Ru Y, Kechris KJ, Tabakoff B, et al. The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res. 2014;42(17):e133. doi:10.1093/nar/gku631
24. Franceschini A, Szklarczyk D, Frankild S, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(Databaseissue):D808–815. doi:10.1093/nar/gks1094
25. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–2504. doi:10.1101/gr.1239303
26. Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40(WebServer issue):W478–483. doi:10.1093/nar/gks402
27. Leek JT, Storey JD, Gibson G. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007;3(9):1724–1735. doi:10.1371/journal.pgen.0030161
28. Xu RH, Wei W, Krawczyk M, et al. Circulating tumour DNA methylation markers for diagnosis and prognosis of hepatocellular carcinoma. Nat Mater. 2017;16(11):1155–1161. doi:10.1038/nmat4997
29. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMC Med. 2015;13:1. doi:10.1186/s12916-014-0241-z
30. Cotto KC, Wagner AH, Feng YY, et al. DGIdb 3.0: a redesign and expansion of the drug-gene interaction database. Nucleic Acids Res. 2018;46(D1):D1068–d1073.
31. Chan BKC. Data analysis using R programming. Adv Exp Med Biol. 2018;1082:47–122.
32. Robin X, Turck N, Hainard A, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011;12(1):77. doi:10.1186/1471-2105-12-77
33. Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, Hasz R, Walters G, Garcia F, Young N, Foster B. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45(6):580–585.
34. Kim GS, Uddin M. Sex-specific and shared expression profiles of vulnerability and resilience to trauma in brain and blood. Biol Sex Differ. 2020;11(1):13.
35. Bian Y, Yang L, Zhao M, et al. Identification of key genes and pathways in post-traumatic stress disorder using microarray analysis. Front Psychol. 2019;10:302.
36. Vidović A, Grubišić-Ilić M, Kozarić-Kovačić D, et al. Exaggerated platelet reactivity to physiological agonists in war veterans with posttraumatic stress disorder. Psychoneuroendocrinology. 2011;36(2):161–172.
37. Tylee DS, Kawaguchi DM, Glatt SJ. On the outside, looking in: a review and evaluation of the comparability of blood and brain “-omes”. Am J Med Genet B Neuropsychiatr Genet. 2013;162b(7):595–603.
38. Tylee DS, Chandler SD, Nievergelt CM, et al. Blood-based gene-expression biomarkers of post-traumatic stress disorder among deployed marines: a pilot study. Psychoneuroendocrinology. 2015;51:472–494.
39. Zieker J, Zieker D, Jatzko A, et al. Differential gene expression in peripheral blood of patients suffering from post-traumatic stress disorder. Mol Psychiatry. 2007;12(2):116–118.
40. van Loo KM, Martens GJ. Genetic and environmental factors in complex neurodevelopmental disorders. Curr Genomics. 2007;8(7):429–444.
41. Bam M, Yang X, Zumbrun EE, et al. Dysregulated immune system networks in war veterans with PTSD is an outcome of altered miRNA expression and DNA methylation. Sci Rep. 2016;6:31209.
42. Tian Z, Dong Q, Wu T, Guo J. MicroRNA-20b-5p aggravates neuronal apoptosis induced by β-amyloid via down-regulation of ras homolog family member C in Alzheimer’s disease. Neurosci Lett. 2021;742:135542.
43. Sun X, Song Z, Si Y, Wang JH. microRNA and mRNA profiles in ventral tegmental area relevant to stress-induced depression and resilience. Prog Neuropsychopharmacol Biol Psychiatry. 2018;86:150–165.
44. Cohen JE, Lee PR, Chen S, Li W, Fields RD. MicroRNA regulation of homeostatic synaptic plasticity. Proc Natl Acad Sci U S A. 2011;108(28):11650–11655.
45. Dong Q, Yuan G, Liu M, et al. Downregulation of microRNA-374a predicts poor prognosis in human glioma. Exp Ther Med. 2019;17(3):2077–2084.
© 2022 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.