Skip to main content

Tensor decomposition-based unsupervised feature extraction identifies candidate genes that induce post-traumatic stress disorder-mediated heart diseases



Although post-traumatic stress disorder (PTSD) is primarily a mental disorder, it can cause additional symptoms that do not seem to be directly related to the central nervous system, which PTSD is assumed to directly affect. PTSD-mediated heart diseases are some of such secondary disorders. In spite of the significant correlations between PTSD and heart diseases, spatial separation between the heart and brain (where PTSD is primarily active) prevents researchers from elucidating the mechanisms that bridge the two disorders. Our purpose was to identify genes linking PTSD and heart diseases.


In this study, gene expression profiles of various murine tissues observed under various types of stress or without stress were analyzed in an integrated manner using tensor decomposition (TD).


Based upon the obtained features,  400 genes were identified as candidate genes that may mediate heart diseases associated with PTSD. Various gene enrichment analyses supported biological reliability of the identified genes. Ten genes encoding protein-, DNA-, or mRNA-interacting proteins—ILF2, ILF3, ESR1, ESR2, RAD21, HTT, ATF2, NR3C1, TP53, and TP63—were found to be likely to regulate expression of most of these  400 genes and therefore are candidate primary genes that cause PTSD-mediated heart diseases. Approximately 400 genes in the heart were also found to be strongly affected by various drugs whose known adverse effects are related to heart diseases and/or fear memory conditioning; these data support the reliability of our findings.


TD-based unsupervised feature extraction turned out to be a useful method for gene selection and successfully identified possible genes causing PTSD-mediated heart diseases.


Post-traumatic stress disorder (PTSD) [1] is primarily a mental illness caused by stressors. Nevertheless, PTSD can cause additional symptoms apparently not directly related to the central nervous system. PTSD-mediated heart diseases are some of such examples [2]. Although it is believed that PTSD highly correlates with heart failure [3], the mechanisms by which PTSD mediates heart failure are still unclear. Because a study on twins revealed a strong correlation between PTSD and heart diseases [4], the genomic factors are believed to link the two disorders. Therefore, in this study, tensor decomposition (TD)-based unsupervised feature extraction (FE)—which is the extension of a recently proposed principal component analysis-based unsupervised FE that has been successfully applied to various bioinformatics problems [522]—was used for various gene expression profiles of murine tissues with the aim to find genes coexpressive or differentially expressed between stressful and unstressful conditions in the brain and heart. As shown in the text below, we identified approximately 400 genes using TD-derived features, and these genes are strongly related to neurodegenerative diseases as well as cardiac-muscle aberrations. Furthermore, the top 10 genes encoding protein-, DNA-, or mRNA-interacting proteins were identified as those governing expression of these  400 genes and are possible therapeutic targets in PTSD-mediated heart diseases according to other reports.


Gene expression

The gene expression profiles used in this study were downloaded from the Gene Expression Omnibus (GEO) database (GEO ID GSE68077). The file “GSE68077_series_matrix.txt” that is available as “Series Matrix File(s)” was downloaded. Probes whose names start with “EA” were removed. The gene expression profile was standardized (i.e., means and variances in each sample are 0 and 1, respectively).

Gene expression profiles were formatted as a tensor, \(x_{i,j_{1},j_{2},j_{3},j_{4}}\), of the ith probe, subjected to j 1th treatment (j 1=1: control, j 1=2: treated [stress-exposed] samples), in the j 2th tissue [ j 2=1: amygdala (AY), j 2=2: hippocampus (HC), j 2=3: medial prefrontal cortex (MPFC), j 2=4: septal nucleus (SE), j 2=5: striatum (ST), j 2=6: ventral striatum (VS), j 2=7: blood, j 2=8: heart, j 2=9: hemibrain, j 2=10: spleen], with the j 3th stress duration (j 3=1: 10 days, j 3=2: five days) and j 4th rest period after application of stress (j 4=1: 1.5 weeks, j 4=2: 24 hours, j 4=3: 6 weeks). Zero values were assigned to missing observations (e.g., measurements at 6 weeks after a 5-day period of stress are not available).

TD-based unsupervised FE

To perform gene selection using gene singular value vectors, \(x_{\ell _{1},i_{1}}\) for synthetic data and \(x_{\ell _{5},i}\) for real gene expression profiles, we have to decide which singular value vectors are used for the selection. If we denote this set of gene singular value vectors as Ω, then a P-value, \(P_{i_{1}}\), for synthetic data and P i for real gene expression profiles, is assigned to each gene by assuming that the singular value vectors of genes obey a χ 2 distribution,

$$P_{i_{1}} = P \left [ > \sum_{\ell_{5} \in \Omega} \left(\frac{x_{\ell_{1},i_{1}}}{\sigma_{\ell_{1}}}\right)^{2}\right] $$

for synthetic datasets and

$$P_{i} = P \left [ > \sum_{\ell_{5} \in \Omega} \left(\frac{x_{\ell_{5},i}}{\sigma_{\ell_{5}}}\right)^{2}\right] $$

for real gene expression profiles, where \(\sigma _{\ell _{1}}\) and \(\sigma _{\ell _{5}}\) are the standard deviation of the 1th and 5th singular value vectors (\(x_{\ell _{1},i_{1}}\) and \(x_{\ell _{5},i}\)), whereas P[>x] is the P-value for the hypothesis that the argument is greater than x, assuming a χ 2 distribution. After that, genes associated with adjusted P-values less than 0.01, 0.05, and 0.1 for the synthetic dataset and 0.01 for real gene expression profiles were selected, respectively.

Computation of adjusted P-values

The adjusted P-values were computed either by the p.adjust function in R [23] with the “BH” (Benjamini–Hochberg) option or the fdrtool in R [23] with the statistic=“pvalue” option. The AUC was computed by means of the colAUC function in the caTools package in R [23].

The synthetic dataset

This dataset is composed of a 30,000×10×10 tensor, \(x_{i_{1},i_{2},i_{3}},1 \le i_{1} \le 30,000, 1 \le i_{2},i_{3}\le 10\). The first (i 1), second (i 2), and third (i 3) modes represent genes, tissues, and treatments, respectively. For each of 10 treatments, 100 genes were expressed in four out of 10 tissues. The first gene through the 100th gene are expressive during the first treatment, the 101st to 200th genes are expressive during the second treatment, and so on. If a combination of a gene set and class falls into a blue filled square in Fig. 1 (e.g., the second gene set in the third class), then \(x_{i_{1},i_{2},i_{3}}\) obeys a Gaussian distribution of mean 4 and variance 1; otherwise, the mean is assumed to be zero. x ij is standardized within each sample as well.

Fig. 1

The gene expression pattern of 10 gene sets. Each of which includes 100 genes from synthetic data (thus, 1000 of the total of 30,000 genes are being considered). The remaining 29,000 genes do not have any class specificity. Blue squares represent classes where the genes in each gene set are expressive. Ten tissues are assumed to be treated in 10 distinct ways. Thus, in total, there are 100 classes

Enrichment analyses

These analyses were conducted by means of g:profiler r1622_e84_eg31 [24] and Enrichr [25]. All probe IDs were converted to gene symbols before uploading to the servers. For g:profiler, all the genes included in the microarray were uploaded as a background.

Enrichment analysis of MSigDB

A total of 457 gene symbols were uploaded to and login are needed). The option “C2: curated gene sets” was selected.

Clustering analysis

For synthetic data, two clustering methods were used. The first is hierarchical clustering (Ward method) using the Euclidean distance between the first gene’s and 10th gene’s singular value vectors, \(x_{\ell _{1},i_{1}},1 \le \ell _{1} \le 10\), as the distance. Then, the generated trees were partitioned into 11 clusters. The Ward method was implemented as the hclust function in R [23] with the method =“ward” option. Partitioning was performed using the cutree function in R using k =11 option (the number of clusters is 11). The second one is a Gaussian mixture, which was carried out by the Mclust function in the mclust [26] package in R [23] with the G =1:11 option (assuming 1 to 11 clusters).


TD applied to tensors of gene expression profiles

In this study, gene expression profiles were regarded as tensors. Gene expression profiles were analyzed in various tissues including the heart and brain, under various conditions (stressful or unstressful), with various periods of stress and rest time after application of a stressor. These datasets were naturally regarded as a tensor, \(x_{i,j_{1},\ldots,j_{m}}\), where i stands for genes and j k ,k=1,…,m denotes various tissues as well as experimental conditions. To reduce the number of degrees of freedom, tensors can be decomposed to smaller tensors, vectors, or matrices. Although there are multiple implementations, higher-order singular value decomposition [27] (HOSVD) was employed in the present study, and a tensor was decomposed as \(x_{i,j_{1},\ldots,j_{m}} = \sum _{\ell _{1},\ldots,\ell _{m+1}} G (\ell _{1}, \ldots, \ell _{m+1}) \cdot x_{\ell _{m+1},i} \prod _{k=1}^{m} x_{\ell _{k},j_{k}}\), where G is the core tensor and \(x_{\ell _{k},j_{k}}\) and \(x_{\ell _{m+1},i}\) are singular value matrices. In this implementation, singular value vectors, \(x_{\ell _{k},j_{k}}\) or \(x_{\ell _{m+1},j}\), associated with G with greater absolute values primarily and correlatively contribute to the original tensor, \(x_{i,j_{1},\ldots,j_{m}}\) (\(x_{\ell _{k},j_{k}}, k=1,\ldots,m, x_{\ell _{m+1},i}\) are supposed to be orthogonal matrices and thus have the same absolute values and equally contribute to the tensor, whereas only the amount of G counts for the contribution).

Synthetic data

Prior to application to gene expression profiles, TD was applied to a synthetic dataset to demonstrate usefulness of our strategy. In this synthetic dataset, 10 tissues were assumed to be treated with 10 experimental conditions (thus, there were assumed to be 10×10=100 samples). In each experimental condition, in four out of 10 tissues, 100 distinct genes were expressed (Fig. 1). Thus, in total, 1,000 genes were expressive in some of the 10 tissues under some of the conditions tested. The remaining majority of genes (as many as 29,000 because in total 30,000 genes are assumed to exist) were not expressed at all in any tissues under any conditions. The task was to identify separately 10 sets of 100 genes as being coexpressive in four tissues.

Applying TD to the synthetic dataset,

$$x_{i_{1},i_{2},i_{3}} = \sum_{\ell_{1},\ell_{2},\ell_{3}} G(\ell_{1},\ell_{2},\ell_{3}) x_{\ell_{1},i_{1}}x_{\ell_{2},i_{2}}x_{\ell_{3},i_{3}},$$

genes were embedded into the space spanned by the derived gene singular value vectors, \(x_{\ell _{1},i_{1}}\). We found that they are strictly clustered coincidently with the 10 presumed clusters (Fig. 2). Although genes identified as outliers by means of these gene singular value vectors, \(x_{\ell _{1},i_{1}}\), were extracted, it was obvious that TD-based unsupervised FE successfully identified some of the 1,000 genes with a relatively smaller number of false positives no matter which adjusted P-values were employed as threshold values (Fig. 3).

Fig. 2

Scatter plots involving the second gene’s through fifth gene’s singular value vectors. \(x_{\ell _{1},i_{1}}, 2 \le \ell _{1} \le 5\), of 1,000 genes (1≤i 1≤1000) that belong to one of the 10 gene sets. These 10 gene sets are represented by distinct combinations of colors and symbols. The 29,000 genes not included in any of the 10 gene sets are omitted for clarity

Fig. 3

Performance of synthetic data (averaged over 100 trials). BH: Benjamini-Hochberg, FDR: false discovery rate. Red curves: true positive rates, black curves: false positive rates, solid curves: P=0.01, dashed curves: P=0.05, dotted curves: P=0.1, blue dash-and-dot curves: area under the curve (AUC)

To test quantitatively whether the 10 clusters of the genes were identified correctly, hierarchical clustering as well as mixture Gaussian clustering were performed. Ten gene clusters were found to be identified correctly, and expression patterns were also correct (Table 1 and Additional file 1).

Table 1 Clustering of genes identified by TD-based unsupervised FE for synthetic data

In the next subsection, TD-based unsupervised FE is applied to real gene expression profiles.

Real gene expression profiles

In the previous subsection, the usefulness of our strategy was successfully demonstrated on synthetic data. In this subsection, TD-based unsupervised FE is applied to real data, i.e., gene expression profiles [28], which were formatted as a five-mode tensor that contains indices corresponding to genes (i) versus tissues (j 2) vs stress duration (j 3) vs rest period after application of stress (j 4) vs control or treatment (j 1) (Table 2). Replicates in each category were averaged within each category before TD was applied. It is definitely a challenging dataset because genes being sought must be differentially expressed between treated and control samples, not in all tissues but in some tissues. HOSVD was applied to the tensor and five singular value matrices, \(x_{\ell _{k},j_{k}}, 1\le k \le 4\), and \(x_{\ell _{5},i}\), were obtained. Figure 4a shows the second control-related or treatment-related singular value vectors, \(x_{\ell _{1}=2,j_{1}}, j_{1}=1,2\). The findings indicate that this expression represents a difference between control (j 1=1) and treated (j 1=2) samples. Next, tissue singular value vectors, \(x_{\ell _{2},j_{2}}, 1 \le \ell _{2}, j_{2} \le 10\), were studied (Fig. 4b and Additional file 2). Then, the fourth tissue singular value vector, \(x_{\ell _{2}=4,j_{2}}, 1 \le j_{2} \le 10\), was found to show coexpression among amygdala (AY, j 2=1), hippocampus (HC, j 2=1), and heart (j 2=8), which represents the phenotypes of interest. After that, we investigated which gene singular value vectors are associated with the fourth sample as well as the second control-related or treatment-related singular value vectors. Table 3 shows the top-ranked core tensor G( 1=2, 2=4, 3, 4, 5) with greater absolute values. Then, gene singular values vectors, \(x_{\ell _{5} \in (1,4,11),i} \), were identified as being associated with the fourth sample as well as the second control-related or treatment-related singular value vectors. After that, 801 probes (Additional file 3) associated with adjusted P-values less than 0.01 were selected as outliers using these three gene singular value vectors.

Fig. 4

Singular value vectors employed a The second control-related or treatment-related singular value vector, \(x_{\ell _{1}=2,j_{1}}\). Control: j 1=1, and treatment (stress): j 1=2. b The fourth tissue singular value vector, \(x_{\ell _{2}=4,j_{2}}\), AY: j 2=1, HC: j 2=2, heart: j 2=8, hemibrain: j 2=9, and spleen: j 2=10. Other tissue singular value vectors, 2≠4, can be found in Additional file 2

Table 2 Samples used in this study Numbers before/after comma are control/treated samples
Table 3 Top-ranked G( 1=2, 2=4, 3, 4, 5) with greater absolute values

To determine whether the 801 selected probes are selectively expressive in the AY, HC, and heart as expected, the t test was applied to all the 40 combinations of control and treated samples. Then, 13 combinations (Table 4) turned out to have the adjusted P-values less than 0.01. Because the AY, HC, and heart are abundantly represented in Table 4, it is obvious that our strategy, TD-based unsupervised FE, successfully identified probes selectively coexpressive in AY, HC, and heart between control and treated samples.

Table 4 Thirteen combinations of tissues and experimental conditions where the selected 801 probes are differentially expressed between Stress-exposed and control samples

Comparison with other methodologies

In contrast to the success of this strategy, which was applied to a synthetic dataset and real gene expression profiles, other tested methods failed to identify some of the 1000 genes correctly in synthetic data (Additional file 4) and failed to identify some of the significantly differentially expressed genes (Additional file 4).


Biological reliability was evaluated by means of 457 gene symbols (Additional file 5) that are associated with the 801 identified probes, uploaded to g:profiler. Various biological terms were enriched (Additional file 6): 198 Gene Ontology (GO) biological process (BP) terms, 79 GO cellular component (CC) terms, 49 GO molecular function (MF) terms, 7 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 38 REACTOME pathways. Among them, neurodegenerative-disease-related KEGG pathways were enriched (three categories): “Huntington’s disease,” “Parkinson’s disease,” and “Alzheimer’s disease” as well as one KEGG heart disease-related pathway, “Cardiac muscle contraction”, GO BP terms, “heart contraction” and “cardiac muscle contraction” as well as GO CC term “neuron part” were also identified. Thus, these findings suggest that the identified genes are potentially related to neuronal functions as well as heart anomalies.

To evaluate the relation between the selected genes and PTSD-mediated heart diseases, two biological terms, “nonsense-mediated decay” (NMD; REACTOME) and “SRP-dependent cotranslational protein targeted to membrane” (REACTOME, GO BP), were further analyzed because these two are reported to be specifically related to fear memory and heart failure (see below).

NMD was reported to regulate cardiac myosin-binding protein C mutant levels in cardiomyopathic mice [29]. Arc mRNA is targeted for NMD, and time-tependent expression of Arc and Zif268 after acquisition of fear conditioning is observed [30]. On the other hand, cardiac involvement is less common and survival is better among patients with anti-SRP [31]. Srp54 is also upregulated after contextual fear conditioning in the rat [32]. Furthermore, SRP is related to NMD [33]. SRP and NMD are associated with abnormal gene expression, e.g., miss-splicing, which is suggested to induce PTSD-mediated heart diseases.

Next, to identify what governed the processes overall, 457 gene symbols were uploaded to Enrichr [25] because 457 genes are too many to be considered primary factors of PTSD-mediated heart diseases; a smaller number of factors is preferable. As a result, many genes listed among “Transcription factors PPI” were found to be enriched (41 genes have the adjusted P-values less than 0.05; see Additional file 7 for the full list). As for the top 10 proteins (Table 5), many are related to heart diseases and fear memories as shown below. Mutation of the most significant gene, ILF3, is related to heart attacks [34]. ESR1, the second most significant gene, was reported to be related to fear conditioning [35] and heart diseases [36]. RAD21, the third most significant gene, is related to memory formation [37] (through genomic structure of BDNF and Arc) as well as to heart diseases [38]. HTT, the fourth most significant gene, is related to heart diseases [39] and its full name is Huntingtin, which is naturally related to the corresponding neurodegenerative disease. ATF2, the seventh most significant gene, was reported to be possibly involved in Alzheimer’s and heart diseases [40].

Table 5 Top 10 proteins interacting with the 457 genes identified by TD-based unsupervised FE, which were listed by Enrichr (“Transcription factors PPI”)

ATF2 and c-Jun are parts of AP-1, which is known to be related to fear extinction [41] as well as extinction of contextual fear memory [42]. NR3C1, the eighth most significant gene, has frequently been reported to be related to fear emmory [43], and its mutation is reported to be related to muscle strength [44]. GR, encoded by NR3C1, is required for fetal heart maturation [45]. NR3C1 was also reported to be one of the driver genes of PTSD [2] in a study involving a search for the genes causing PTSD-mediated heart diseases. The ninth most significant gene, TP53, is one of upstream regulators of pathways associated with the onset of memory deficits in mice [46], although nothing was reported to the tenth significant TP63.

Figure 5 shows the graphical representation of enriched biological terms as well as the 10 above-mentioned interacting proteins. The figure indicates that they are tightly inter-related. Thus, the 10 identified proteins are likely to regulate expression of genes enriched in these biological terms and PTSD-mediated heart diseases as well although additional experimental validation is needed.

Fig. 5

Graphical representation of relations between the identified biological terms and proteins. Biological temrs (orange) and various protein-, DNA-, or mRNA-binding proteins (cyan). A: “heart contraction” (GO BP), B: “cardiac muscle contraction” (GO BP), C: “protein targeted to ER” (GO BP), D: “SRP-dependent cotranslational protein targeted to membrane” (GO BP), E: “neuron part” (GO CC), F: “Huntington’s disease” (KEGG), G: “Parkinson’s disease” (KEGG), H: “Alzheimer’s disease” (KEGG), I: “Cardiac muscle contraction” (KEGG), J: “nonsense-mediated decay (NMD)” (REACTOME), K: “SRP-dependent cotranslational protein targeted to membrane” (REACTOME). Orange edges: genes shared by biological terms, blue edges: genes targeted by protein-, DNA-, or mRNA-binding proteins. Width of edges is proportional to the number of genes. Sizes of the orange circles representing biological terms are proportional to the number of genes enriched in each biological term among the 457 genes

To confirm correctness of identification of the enriched biological terms, the 457 gene symbols were also uploaded to MSigDB [47] (Additional file 8). As a result, “Nonsense-mediated decay enhanced by the exon junction complex,” “SRP-dependent cotranslational protein targeted to membrane” (REACTOME), “Parkinson’s disease,” “Alzheimer’s disease,” “Huntington’s disease,” and “Cardiac muscle contraction” (KEGG) were found to be significantly enriched. Therefore, the identified enrichment of these biological terms is trustworthy.

Finally, to confirm the reliability of the selected genes, DrugMatrix in Enricher was analyzed. Many compounds that affect gene expression in the rodent heart were identified; in total, 7,098 combinations of drugs with various dose densities and solvents were found to have the adjusted P-values less than 0.01 (Additional file 9). First of all, all the top 10 combinations (Table 6) were found to decrease expression of some genes in the heart, in agreement with the expectation that the identified genes should be related to the heart because they are supposed to contribute to heart failure. In addition, many adverse effects caused by these drugs, as shown in Table 6, are also associdated with PTSD-mediated heart diseases. For example, the most significant drug, low-dose prednisolone [48], increases long-term risk of ischemic cerebrovascular events. Long-term administration of the second most significant drug, ethosuximide, adversely affects fear memory [49]. The third most significant compound, caffeine, is known to be related to heart diseases [50] as well as acquisition and retention of Pavlovian conditioned freezing [51]. The fourth most significant drug, clomipramine, was once suggested to be used for the treatment of anxiety [52]. The fifth most significant drug, prednisolone, was reported to alleviate adverse cardiac effects [53]. Several cases of cardiac adverse reactions related to the seventh and ninth most significant drug, vinorelbine, have been reported in the literature [54]. The eighth most significant drug, atropine, affects heart rate [55]. Cardiac arrest was reported after administration of the tenth most significant drug, oxymetazoline (nasal spray) [56]. These relations between heart problems or fear memory and drugs downregulating expression of selected genes in the heart support the reliability of our findings, too.

Table 6 Top 10 significant drugs identified by DrugMatrix in Enricher

Before closing this subsection, I would like to comment on some points. First, comparisons with some related works. Since Vaccarino et al [4] has clearly denoted that there are limited number of mutated genes shared between PTSD and congenital heart defect (CHD), it might not look reasonable that I argued that genomic background was important. However, even if there are no shared mutated genes between PTSD and CHD, genomic background can induce association between PTSD and CHD. For examples, there are two genes A and B. A is a CHD causing genes and B is regulating A. Then, even if mutation of B takes place not in CHD but in PTSD, genomic background (i.e., mutation of gene B) still can induce CHD. This means that shared mutated genes is not only genomic background that can induce the association between PTSD and CHD. Cho et al [57] also investigated mRNA and miRNA expression of stressed mouse heart. Nevertheless, since we have extensively studied this study in our previous study [7], I did not discuss about it in the present paper. Finally, Pollard et al [2] identified 37 mutated genes shared between PTSD and cardiovascular disease (CVD). However, as in the case of PTSD and CHD, shared mutated genes are not only factors that can mediate PTSD mediated heart disease. Actually, there are no significant overlaps between these 37 genes and our 457 genes. Possibly, our identified gene expression alteration between PTSD and controls is not due to a direct effect of shared mutated genes but due to more complicated indirect effect. Second, I would like to briefly argue about how TD can figure out hidden relations among genes and diseases. From the mathematical point of views, TD is nothing but possible assumption. Thus, the validation of methodology can be done only based on the goodness of outcomes. Since our results are biologically reliable, our assumption that gene expression has hidden structure that can be figured out TD seems to be correct. More applications of this strategy will add more confidence to TD based unsupervised FE in the future.


In this paper, TD-based unsupervised FE was applied to murine tissue gene expression profiles with and without stress conditions. The resulting 457 genes associated with 801 probes identified as outliers using gene singular value vectors were subjected to various enrichment analyses. Ten proteins likely to regulate expression of these genes are proposed here as possible causal genes of PTSD-mediated heart diseases.







Biological process


Cellular component


Feature extraction


Gene expression omnibus


Gene Ontology




Higher-order singular value decomposition


Kyoto Encyclopedia of Genes and Genomes


Molecular function


Medial prefrontal cortex


Nonsense-mediated decay


Post-traumatic stress disorder


Septal nucleus


Signal recognition particle




Tensor decomposition


Ventral striatum


  1. 1

    Kirkpatrick HA, Heller GM. Post-traumatic stress disorder: theory and treatment update. Int J Psychiatry Med. 2014; 47(4):337–46.

    Article  PubMed  Google Scholar 

  2. 2

    Pollard HB, Shivakumar C, Starr J, Eidelman O, Jacobowitz DM, Dalgard CL, Srivastava M, Wilkerson MD, Stein MB, Ursano RJ. “Soldier’s Heart”: A Genetic Basis for Elevated Cardiovascular Disease Risk Associated with Post-traumatic Stress Disorder. Front Mol Neurosci. 2016; 9:87.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3

    Deng LX, Khan AM, Drajpuch D, Fuller S, Ludmir J, Mascio CE, Partington SL, Qadeer A, Tobin L, Kovacs AH, Kim YY. Prevalence and Correlates of Post-traumatic Stress Disorder in Adults With Congenital Heart Disease. Am J Cardiol. 2016; 117(5):853–7.

    Article  PubMed  Google Scholar 

  4. 4

    Vaccarino V, Goldberg J, Rooks C, Shah AJ, Veledar E, Faber TL, Votaw JR, Forsberg CW, Bremner JD. Post-traumatic stress disorder and incidence of coronary heart disease: a twin study. J Am Coll Cardiol. 2013; 62(11):970–8.

    Article  PubMed  Google Scholar 

  5. 5

    Taguchi YH. Principal component analysis based unsupervised feature extraction applied to publicly available gene expression profiles provides new insights into the mechanisms of action of histone deacetylase inhibitors. Neuroepigenetics. 2016; 8:1–18. doi:10.1016/j.nepig.2016.10.001.

    Article  Google Scholar 

  6. 6

    Taguchi YH. microrna-mrna interaction identification in wilms tumor using principal component analysis based unsupervised feature extraction. In: 2016 IEEE 16th International Conference on Bioinformatics and Bioengineering (BIBE): 2016. p. 71–8. doi:10.1109/BIBE.2016.14.

  7. 7

    Taguchi YH, Iwadate M, Umeyama H. Principal component analysis-based unsupervised feature extraction applied to in silico drug discovery for posttraumatic stress disorder-mediated heart disease. BMC Bioinformatics. 2015; 16:139.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8

    Taguchi YH, Okamoto A. Principal component analysis for bacterial proteomic analysis In: Shibuya T, Kashima H, Sese J, Ahmad S, editors. Pattern Recognition in Bioinformatics. LNCS. vol. 7632. Heidelberg: Springer: 2012. p. 141–52.

    Google Scholar 

  9. 9

    Ishida S, Umeyama H, Iwadate M, Taguchi YH. Bioinformatic Screening of Autoimmune Disease Genes and Protein Structure Prediction with FAMS for Drug Discovery. Protein Pept Lett. 2014; 21(8):828–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10

    Kinoshita R, Iwadate M, Umeyama H, Taguchi YH. Genes associated with genotype-specific DNA methylation in squamous cell carcinoma as candidate drug targets. BMC Syst Biol. 2014; 8 Suppl 1:4.

    Article  Google Scholar 

  11. 11

    Taguchi YH, Murakami Y. Principal component analysis based feature extraction approach to identify circulating microRNA biomarkers. PLoS ONE. 2013; 8(6):66714.

    Article  Google Scholar 

  12. 12

    Taguchi YH, Murakami Y. Universal disease biomarker: can a fixed set of blood microRNAs diagnose multiple diseases?. BMC Res Notes. 2014; 7:581.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13

    Murakami Y, Toyoda H, Tanahashi T, Tanaka J, Kumada T, Yoshioka Y, Kosaka N, Ochiya T, Taguchi YH. Comprehensive mimaNA expression analysis in peripheral blood can diagnose liver disease. PLoS ONE. 2012; 7(10):48366.

    Article  Google Scholar 

  14. 14

    Murakami Y, Tanahashi T, Okada R, Toyoda H, Kumada T, Enomoto M, Tamori A, Kawada N, Taguchi YH, Azuma T. Comparison of Hepatocellular Carcinoma miRNA Expression Profiling as Evaluated by Next Generation Sequencing and Microarray. PLoS ONE. 2014; 9(9):106314.

    Article  Google Scholar 

  15. 15

    Murakami Y, Kubo S, Tamori A, Itami S, Kawamura E, Iwaisako K, Ikeda K, Kawada N, Ochiya T, Taguchi YH. Comprehensive analysis of transcriptome and metabolome analysis in Intrahepatic Cholangiocarcinoma and Hepatocellular Carcinoma. Sci Rep. 2015; 5:16294.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16

    Umeyama H, Iwadate M, Taguchi YH. TINAGL1 and B3GALNT1 are potential therapy target genes to suppress metastasis in non-small cell lung cancer. BMC Genomics. 2014; 15 Suppl 9:2.

    Article  Google Scholar 

  17. 17

    Taguchi YH, Iwadate M, Umeyama H. Heuristic principal component analysis-based unsupervised feature extraction and its application to gene expression analysis of amyotrophic lateral sclerosis data sets. In: Computational Intelligence in Bioinformatics and Computational Biology (CIBCB), 2015 IEEE Conference On: 2015. p. 1–10. doi:10.1109/CIBCB.2015.7300274.

  18. 18

    Taguchi YH, Iwadate M, Umeyama H, Murakami Y, Okamoto A. Heuristic principal component analysis-aased unsupervised feature extraction and its application to bioinformatics In: Wang B, Li R, Perrizo W, editors. Big Data Analytics in Bioinformatics and Healthcare. Hershey: IGI Global: 2015. p. 138–62.

    Google Scholar 

  19. 19

    Taguchi YH. Identification of aberrant gene expression associated with aberrant promoter methylation in primordial germ cells between E13 and E16 rat F3 generation vinclozolin lineage. BMC Bioinformatics. 2015; 16 Suppl 18:16.

    Article  Google Scholar 

  20. 20

    Taguchi YH. Identification of More Feasible MicroRNA-mRNA Interactions within Multiple Cancers Using Principal Component Analysis Based Unsupervised Feature Extraction. Int J Mol Sci. 2016; 17(5):696.

    Article  PubMed Central  Google Scholar 

  21. 21

    Taguchi YH. Principal component analysis based unsupervised feature extraction applied to budding yeast temporally periodic gene expression. BioData Min. 2016; 9:22.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22

    Taguchi YH, Iwadate M, Umeyama H. SFRP1 is a possible candidate for epigenetic therapy in non-small cell lung cancer. BMC Med Genomics. 2016; 9 Suppl 1:28.

    Article  PubMed  Google Scholar 

  23. 23

    R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2015. R Foundation for Statistical Computing.

  24. 24

    Reimand J, Arak T, Adler P, Kolberg L, Reisberg S, Peterson H, Vilo J. g:Profiler-a web server for functional interpretation of gene lists (2016 update). Nucleic Acids Res. 2016; 44(W1):83–9.

    Article  Google Scholar 

  25. 25

    Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, Gundersen GW, Ma’ayan A. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016; 44(W1):90–7.

    Article  Google Scholar 

  26. 26

    Fraley C, Raftery AE. Model-based Clustering, Discriminant Analysis and Density Estimation. J Am Stat Assoc. 2002; 97:611–31.

    Article  Google Scholar 

  27. 27

    Lathauwer LD, Moor BD, Vandewalle J. A multilinear singular value decomposition. SIAM J Matrix Anal Appl. 2000; 21(4):1253–78. doi:10.1137/S0895479896305696

  28. 28

    Muhie S, Gautam A, Meyerhoff J, Chakraborty N, Hammamieh R, Jett M. Brain transcriptome profiles in mouse model simulating features of post-traumatic stress disorder. Mol Brain. 2015; 8:14.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29

    Vignier N, Schlossarek S, Fraysse B, Mearini G, Kramer E, Pointu H, Mougenot N, Guiard J, Reimer R, Hohenberg H, Schwartz K, Vernet M, Eschenhagen T, Carrier L. Nonsense-mediated mRNA decay and ubiquitin-proteasome system regulate cardiac myosin-binding protein C mutant levels in cardiomyopathic mice. Circ Res. 2009; 105(3):239–48.

    CAS  Article  PubMed  Google Scholar 

  30. 30

    Lonergan ME, Gafford GM, Jarome TJ, Helmstetter FJ. Time-dependent expression of Arc and zif268 after acquisition of fear conditioning. Neural Plast. 2010; 2010:139891.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31

    Kao AH, Lacomis D, Lucas M, Fertig N, Oddis CV. Anti-signal recognition particle autoantibody in patients with and patients without idiopathic inflammatory myopathy. Arthritis Rheum. 2004; 50(1):209–15.

    CAS  Article  PubMed  Google Scholar 

  32. 32

    Federighi G, Traina G, Macchi M, Ciampini C, Bernardi R, Baldi E, Bucherelli C, Brunelli M, Scuri R. Modulation of gene expression in contextual fear conditioning in the rat. PLoS ONE. 2013; 8(11):80037.

    Article  Google Scholar 

  33. 33

    Popp MW, Maquat LE. Organizing principles of mammalian nonsense-mediated mRNA decay. Annu Rev Genet. 2013; 47:139–65.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34

    Yoshida T, Kato K, Oguri M, Horibe H, Kawamiya T, Yokoi K, Fujimaki T, Watanabe S, Satoh K, Aoyagi Y, Tanaka M, Yoshida H, Shinkai S, Nozawa Y, Yamada Y. Association of polymorphisms of BTN2A1 and ILF3 with myocardial infarction in Japanese individuals with or without hypertension, diabetes mellitus or chronic kidney disease. Int J Mol Med. 2011; 27(5):745–52.

    CAS  PubMed  Google Scholar 

  35. 35

    Cho J, Yu NK, Choi JH, Sim SE, Kang SJ, Kwak C, Lee SW, Kim JI, Choi DI, Kim VN, Kaang BK. Multiple repressive mechanisms in the hippocampus during memory formation. Science. 2015; 350(6256):82–7.

    CAS  Article  PubMed  Google Scholar 

  36. 36

    Kunnas T, Silander K, Karvanen J, Valkeapaa M, Salomaa V, Nikkari S. ESR1 genetic variants, haplotypes and the risk of coronary heart disease and ischemic stroke in the Finnish population: a prospective follow-up study. Atherosclerosis. 2010; 211(1):200–2.

    CAS  Article  PubMed  Google Scholar 

  37. 37

    Sams DS, Nardone S, Getselter D, Raz D, Tal M, Rayi PR, Kaphzan H, Hakim O, Elliott E. Neuronal CTCF Is Necessary for Basal and Experience-Dependent Gene Regulation, Memory Formation, and Genomic Structure of BDNF and Arc. Cell Rep. 2016; 17(9):2418–30.

    CAS  Article  PubMed  Google Scholar 

  38. 38

    Schuster K, Leeke B, Meier M, Wang Y, Newman T, Burgess S, Horsfield JA. A neural crest origin for cohesinopathy heart defects. Hum Mol Genet. 2015; 24(24):7005–16.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39

    Comings DE, MacMurray JP, Gonzalez N, Ferry L, Peters WR. Association of the serotonin transporter gene with serum cholesterol levels and heart disease. Mol Genet Metab. 1999; 67(3):248–53.

    CAS  Article  PubMed  Google Scholar 

  40. 40

    Liu H, Deng X, Shyu YJ, Li JJ, Taparowsky EJ, Hu CD. Mutual regulation of c-Jun and ATF2 by transcriptional activation and subcellular localization. EMBO J. 2006; 25(5):1058–69.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41

    Guedea AL, Schrick C, Guzman YF, Leaderbrand K, Jovasevic V, Corcoran KA, Tronson NC, Radulovic J. ERK-associated changes of AP-1 proteins during fear extinction. Mol Cell Neurosci. 2011; 47(2):137–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42

    Radwanska K, Schenatto-Pereira G, Ziôkowska M, Łukasiewicz K, Giese KP. Mapping fear memory consolidation and extinction-specific expression of JunB. Neurobiol Learn Mem. 2015; 125:106–12.

    CAS  Article  PubMed  Google Scholar 

  43. 43

    Maddox SA, Schafe GE, Ressler KJ. Exploring epigenetic regulation of fear memory and biomarkers associated with post-traumatic stress disorder. Front Psychiatry. 2013; 4:62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44

    Ash GI, Kostek MA, Lee H, Angelopoulos TJ, Clarkson PM, Gordon PM, Moyna NM, Visich PS, Zoeller RF, Price TB, Devaney JM, Gordish-Dressman H, Thompson PD, Hoffman EP, Pescatello LS. Glucocorticoid Receptor (NR3C1) Variants Associate with the Muscle Strength and Size Response to Resistance Training. PLoS ONE. 2016; 11(1):0148112.

    Article  Google Scholar 

  45. 45

    Rog-Zielinska EA, Thomson A, Kenyon CJ, Brownstein DG, Moran CM, Szumska D, Michailidou Z, Richardson J, Owen E, Watt A, Morrison H, Forrester LM, Bhattacharya S, Holmes MC, Chapman KE. Glucocorticoid receptor is required for foetal heart maturation. Hum Mol Genet. 2013; 22(16):3269–82.

    CAS  Article  PubMed  Google Scholar 

  46. 46

    Neuner SM, Wilmott LA, Hoffmann BR, Mozhui K, Kaczorowski CC. Hippocampal proteomics defines pathways associated with memory decline and resilience in normal aging and Alzheimer’s disease mouse models. Behav Brain Res. 2017; 322(Pt B):288–298.

    CAS  Article  PubMed  Google Scholar 

  47. 47

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005; 102(43):15545–50. doi:10.1073/pnas.0506580102.

  48. 48

    Ajeganova S, Svensson B, Hafstrom I, Andersson M, Forslind K, Catharina K, Leden I, Lindell B, Petersson I, Schaufelberger C, Teleman A, Theander J. Low-dose prednisolone treatment of early rheumatoid arthritis and late cardiovascular outcome and survival: 10-year follow-up of a 2-year randomised trial. BMJ Open. 2014; 4(4):004259.

    Article  Google Scholar 

  49. 49

    Ponnusamy R, Pradhan N. The effects of chronic administration of ethosuximide on learning and memory: a behavioral and biochemical study on nonepileptic rats. Behav Pharmacol. 2006; 17(7):573–80.

    CAS  Article  PubMed  Google Scholar 

  50. 50

    Cornelis MC, El-Sohemy A. Coffee, caffeine, and coronary heart disease. Curr Opin Clin Nutr Metab Care. 2007; 10(6):745–51.

    CAS  Article  PubMed  Google Scholar 

  51. 51

    Dubroqua S, Low SR, Yee BK, Singer P. Caffeine impairs the acquisition and retention, but not the consolidation of Pavlovian conditioned freezing in mice. Psychopharmacol (Berl). 2015; 232(4):721–31.

    CAS  Article  Google Scholar 

  52. 52

    Ninan PT. The functional anatomy, neurochemistry, and pharmacology of anxiety. J Clin Psychiatry. 1999; 60 Suppl 22:12–17.

    CAS  PubMed  Google Scholar 

  53. 53

    Bauer R, Blain A, Greally E, Lochmuller H, Bushby K, MacGowan GA, Straub V. Attenuation of adverse cardiac effects in prednisolone-treated delta-sarcoglycan-deficient mice by mineralocorticoid-receptor-antagonism. Neuromuscul Disord. 2010; 20(1):21–8.

    Article  PubMed  Google Scholar 

  54. 54

    Lapeyre-Mestre M, Gregoire N, Bugat R, Montastruc JL. Vinorelbine-related cardiac events: a meta-analysis of randomized clinical trials. Fundam Clin Pharmacol. 2004; 18(1):97–105.

    CAS  Article  PubMed  Google Scholar 

  55. 55

    Greenberg S, Plummer C, Maisenbacher H, Friary J, Berg A. The effect of topical ophthalmic 1 rhythm in normal dogs. Vet Ophthalmol. 2015; 18(2):105–8.

    CAS  Article  PubMed  Google Scholar 

  56. 56

    Thrush DN. Cardiac arrest after oxymetazoline nasal spray. J Clin Anesth. 1995; 7(6):512–4.

    CAS  Article  PubMed  Google Scholar 

  57. 57

    Cho JH, Lee I, Hammamieh R, Wang K, Baxter D, Scherler K, Etheridge A, Kulchenko A, Gautam A, Muhie S, Chakraborty N, Galas DJ, Jett M, Hood L. Molecular evidence of stress-induced acute heart injury in a mouse model simulating posttraumatic stress disorder. Proc Natl Acad Sci U S A. 2014; 111(8):3188–93.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


Nothing needs to be acknowledged.


No funding supported this study. Publication charge was supported by Chuo university research project for specific subjects.

Availability of data and materials

All the datasets were retrieved from a public-domain database.

About this supplement

This article has been published as part of BMC Medical Genomics Volume 10 Supplement 4, 2017: 16th International Conference on Bioinformatics (InCoB 2017): Medical Genomics. The full contents of the supplement are available online at

Authors’ contributions

YHT planned the study, performed all the analyses, and wrote the paper.

Ethics approval and consent to participate

Nothing related to these is included in this study.

Consent for publication

Nothing related is included in this study.

Competing interests

The author declares that he has no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Y.-H. Taguchi.

Additional files

Additional file 1

Mean expression profiles of 10 clusters in synthetic data. Mean gene expression profiles of the 10 clusters from Table 1 for synthetic data when the Ward method was employed. (PDF 17 kb)

Additional file 2

Other tissue singular value vectors. \(x_{\ell 1_{2} eq 4,j_{2}}\) for gene expression profiles. (PDF 10 kb)

Additional file 3

801 probe IDs. Probe IDs identified by TD-based unsupervised FE. (CSV 11 kb)

Additional file 4

Comparison with other methods. A supporting document about the comparisons with other methods. (PDF 948 kb)

Additional file 5

457 gene symbols. Gene symbols associated with 801 probes. (CSV 3 kb)

Additional file 6

g:profiler. A full list of enrichment analyses of g:profiler. (XLSX 697 kb)

Additional file 7

Transcription factors PPI. A full list of enrichment analyses of transcription factors PPI provided by Enrichr. (CSV 70 kb)

Additional file 8

MSigDB. A full list of enrichment analyses of DrugMatrix provided by MSigDB. (CSV 193 kb)

Additional file 9

DrugMatrix. A full list of enrichment analyses of DrugMatrix provided by Enrichr. (CSV 4 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Taguchi, Y. Tensor decomposition-based unsupervised feature extraction identifies candidate genes that induce post-traumatic stress disorder-mediated heart diseases. BMC Med Genomics 10, 67 (2017).

Download citation


  • Tensor decomposition
  • Feature extraction
  • Post-traumatic stress disorder
  • Heart disease