Original Research ARTICLE
A Pathway Proteomic Profile of Ischemic Stroke Survivors Reveals Innate Immune Dysfunction in Association with Mild Symptoms of Depression – A Pilot Study
- 1Occupational Therapy, College of Science Health and Engineering, School of Allied Health, La Trobe University, Melbourne, VIC, Australia
- 2Neurorehabilitation and Recovery, Stroke, The Florey Institute of Neuroscience and Mental Health, Melbourne, VIC, Australia
- 3School of Psychology and Public Health, La Trobe University, Melbourne, VIC, Australia
- 4School of Molecular Sciences, La Trobe University, Melbourne, VIC, Australia
- 5School of Medicine, University of Tasmania, Hobart, TAS, Australia
- 6Commonwealth Science and Industrial Research Organisation (CSIRO), Melbourne, VIC, Australia
- 7The Florey Institute of Neuroscience and Mental Health, Parkville, VIC, Australia
- 8Monash University, Clayton, VIC, Australia
- 9The University of Melbourne, Parkville, VIC, Australia
- 10Department of Medicine, Melbourne Brain Centre, Royal Melbourne Hospital, Melbourne, VIC, Australia
Depression after stroke is a common occurrence, raising questions as to whether depression could be a long-term biological and immunological sequela of stroke. Early explanations for post-stroke depression (PSD) focused on the neuropsychological/psychosocial effects of stroke on mobility and quality of life. However, recent investigations have revealed imbalances of inflammatory cytokine levels in association with PSD, though to date, there is only one published proteomic pathway analysis testing this hypothesis. Thus, we examined the serum proteome of stroke patients (n = 44, mean age = 63.62 years) and correlated these with the Montgomery–Åsberg Depression Rating Scale (MADRS) scores at 3 months post-stroke. Overall, the patients presented with mild depression symptoms on the MADRS, M = 6.40 (SD = 7.42). A discovery approach utilizing label-free relative quantification was employed utilizing an LC-ESI–MS/MS coupled to a LTQ-Orbitrap Elite (Thermo-Scientific). Identified peptides were analyzed using the gene set enrichment approach on several different genomic databases that all indicated significant downregulation of the complement and coagulation systems with increasing MADRS scores. Complement and coagulation systems are traditionally thought to play a key role in the innate immune system and are established precursors to the adaptive immune system through pro-inflammatory cytokine signaling. Both systems are known to be globally affected after ischemic or hemorrhagic stroke. Thus, our results suggest that lowered complement expression in the periphery in conjunction with depressive symptoms post-stroke may be a biomarker for incomplete recovery of brain metabolic needs, homeostasis, and inflammation following ischemic stroke damage. Further proteomic investigations are now required to construct the temporal profile, leading from acute lesion damage to manifestation of depressive symptoms. Overall, the findings provide support for the involvement of inflammatory and immune mechanisms in PSD symptoms and further demonstrate the value and feasibility of the proteomic approach in stroke research.
Over 15 million people worldwide experience a stroke each year; 5 million of those events are fatal, and another 5 million people are left with a permanent disability (1). Previous epidemiological reviews have concluded that ~30% of stroke survivors are likely to experience post-stroke depression (PSD) (2, 3). Prevalence of depression is stroke survivors is reported to peak at 3 months post-stroke based on testing using the Diagnostic and Statistical Manual V [DSM-V (4)] (5). Studies have also characterized depression symptoms on the same criteria as early as 15 days and as late as 12 months post-stroke (5). Stroke patients with PSD show poorer functional and recovery outcomes compared to patients not suffering depression (6, 7). PSD is also found to contribute to poorer quality of life and increased mortality rates (1, 6), thus highlighting the need for patient care that extends beyond that of physical and cognitive rehabilitation. Investigation of biomarkers linked with underlying mechanisms has potential to guide the targeting of therapy, both prevention and treatment.
Currently, there is no consensus regarding the etiology of PSD with much debate as to the extent to which PSD stems from a purely biological origin and/or more likely incorporates elements of psychosocial response (7, 8). Following the biological argument for PSD progression, it was first proposed that the location of the stroke lesion could predict PSD presentation (9). This has been long debated without resolution (10). More recently, basal ganglia or frontal lobe lesions (7), white matter hyperintensities (11), and interruption of connecting pathways (12) have been linked to PSD. Furthermore, a recent review suggests that right hemisphere stroke predicts PSD incidence in the sub-acute, 1–6 months post-stroke period (13). Other explanations for correlations between PSD symptoms and stroke lesions have suggested that the overall balance of monoamines, such as serotonin, dopamine, and norepinephrine, are disrupted following cerebrovascular damage (14). These same amines have previously been associated with depression, consistent with a link between stroke and depressive symptoms (15).
Other reviews and meta-analyses have failed to find clear evidence for lesion location as a risk factor for PSD (3, 5). For instance, similar incidence rates of comorbid depression can be seen in other cerebrovascular diseases such as vascular dementia in the elderly where brain lesions are common (16, 17). Surprisingly, prevalence rates of depression are also similar in patients who have been subject to transient ischemic attacks (TIAs) (18–20) or carotid artery stenosis patients (21), where there is little to no presentation of lesions visible on computerized axial tomography or magnetic resonance imaging (MRI). In TIA especially, prevalence of PSD is comparable to stroke even at 12 months post-stroke (19). Thus, while currently detectable brain lesions may contribute to the etiology of PSD, they are unlikely to be the primary cause with evidence suggesting other pathogenesis involving immune disruption and circulating cytokines long after the stroke (22).
It is possible that PSD could be viewed as depression resulting first from damage caused by the initial ischemia and reperfusion injury (23) and exacerbated by psychosocial issues such as anxiety, loss of confidence, and apathy associated with even mild loss of mobility or functioning (24). In line with this view, the cytokine hypothesis proposes a role for pro-inflammatory cytokines such as tumor necrosis factor α, interleukin-1 (IL-1), IL-6, IL-8, and anti-inflammatory IL-10. This is consistent with evidence that pro-inflammatory cytokines have been found to be increased in both depressive mood (25, 26) and cerebrovascular events (8, 16, 27, 28). Pro-inflammatory cytokine levels are elevated during and long after the therapeutic window of typical tissue plasminogen activator treatment (tPA; <4.5 h) (16, 29, 30). Recent immunoassay research has established that cytokine levels often remained significantly elevated at 1 year post-stroke, leading the authors to suggest that the homeostatic balance of pro- and anti-inflammatory cytokines may be disrupted in the long term (31). Elevated pro-inflammatory cytokine expression is also likely to influence glucocorticoid resistance and sensitivity, which is reported to lead to overactivation of the hypothalamus–pituitary–adrenal (HPA) axis, thereby inducing depressive mood (32–34). Furthermore, persistently high concentrations of pro-inflammatory cytokine levels may provide functional markers for long-term impaired immunological responses as the body attempts to resolve the lesion-induced neurodegeneration.
Given the multifaceted nature of PSD, it is difficult to ascertain the underlying molecular basis of stroke damage associated with depressive symptoms at the 3-month peak prevalence period in human patients. In recent years, proteomics has provided a powerful platform for pathology research. Downstream from transcriptomic processes, the proteome describes the totality of proteins that can be produced by the organism’s genome. The complexity of the proteome accounts for the biological functioning of an organism and as such includes not only the primary function of the protein itself but also how it fits into the biological environment that is further determined by protein–protein interactions, posttranslational modifications (35), and protein degradation. With the addition of alternative splicing to these processes, the 25,000 genes in the human genome can theoretically encode up to 1,000,000 different proteins (36). To this end, proteomics approaches coupled with bioinformatics have proven useful in identifying potential therapeutic biomarkers in pathologies including cancer (37) and cardiovascular disease (36). In the case of proteomic disease research, biomarkers refer to a subset of proteins that are upregulated, downregulated, activated, or deactivated and are detectable using a particular methodology in a disease phenotype.
The proteomic composition of a biological organism can reflect not only rapid changes in response to external challenges (38, 39) but also can be studied to understand disease progression and biological changes in the long term (40, 41). In stroke, the damage caused by the event would be expected to show both acute and later chronic consequences. A discovery proteomics approach utilizing a mass spectrometer with high dynamic range in order to assay as much of the human proteome as possible is therefore warranted in the context of proteome complexity and the multifactorial nature of PSD etiology. Previous studies adopting this approach have been successful in providing new insights and support for the interpretation of molecular biological mechanisms following a stroke (42). The majority of these studies focused on revealing biomarkers that might be useful aids to rapid diagnosis, for example, in the case of differentiation between hemorrhagic stroke versus ischemic stroke (43) and ischemic stroke versus healthy controls (44). Research has also been successful in characterizing the treatment effects of tPA (45) and electroacupuncture (46) on the human proteome, demonstrating that proteomics approaches can objectively characterize the molecular changes associated with certain treatments.
A review of the literature up to December 2015 revealed only a single published study examining the proteomic profile of PSD. Isobaric tagging for relative and absolute quantification (iTRAQ) was employed by Zhan and colleagues (47) to compare the ethylenediaminetetraacetic acid (EDTA) anticoagulated blood of patients with stroke and PSD, stroke without PSD, and healthy controls. The iTRAQ approach involves isobaric tagging of tryptic digested samples based on phenotypes of interest for comparison through commercially available iTRAQ reagents and thus facilitates proteome to proteome comparisons of disease versus control groups (48). This is conducted by introducing a pooled mixture of all the samples to the mass spectrometer, whereby relative quantitation can be assessed by comparing differences in the peak intensity between labeled peptides (49), an approach well suited for investigations of disease phenotypes. The Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics tool for gene enrichment revealed that of the peptides that statistically differentiated stroke patients with and without PSD, only the complement and coagulation cascade pathway accounted for the clustering of identified proteins. Given that this research was conducted on patients at 1 month post-stroke, it was suggested that these findings indicate a homeostatic imbalance of pro-inflammatory and anti-inflammatory processes in PSD (47). These findings are in line with the principles of the cytokine hypothesis of PSD, albeit complement pathways are more reliably identifiable in blood as biomarkers of immune disruption (50, 51). Zhan and colleagues (47) also further examined the top protein candidates by protein immunoblotting of apolipoproteins A IV (ApoA-4) and C-II (ApoC-2), C-reactive protein (CRP), gelsolin, haptoglobin, and leucine-rich glycoprotein (LRG). All of these proteins were significantly altered in PSD patients compared to stroke without depression. As there are limited data on the biological processes directly involved in PSD, Zhan et al.’s (47) study represents a successful first step in applying proteomics approaches to PSD, a disorder that has been previously difficult to characterize in molecular terms.
As the molecular basis of depression after stroke is not well understood, we aimed to take a discovery approach to identifying molecular pathways impacted by depression in stroke survivors. We sought to use a proteomics technique that is capable of identifying and quantifying a large number of biological entities. Our research aim was to identify and quantify peptides to be used to understand the biological mechanisms associated with PSD symptoms. Thus, a data-driven discovery proteomics methodology was used to investigate the biological mechanisms associated with depression etiology 3 months post-stroke. Given the recent finding of Zhan et al. (47), we also sought to investigate whether this finding could be replicated using a different proteomics and analytic approach to that previously reported, in order to strengthen the validity of the proteomic profile of PSD. Our analysis was conducted at 3 months post-stroke.
First, a label-free quantitation (LFQ) technique was employed to examine the blood samples of stroke patients. The LFQ technically differs from the iTRAQ approach in that the protein samples are not pooled or tagged for analysis, and each sample generates its own proteomic profile (52). The relationship between protein expression and various clinical measurements can then be explored instead of comparisons being planned a priori, thus conforming to a discovery approach better suiting our aim. Furthermore, it is debatable that iTRAQ provides better quantitation than LFQ approaches, with LFQ shown to perform similarly (53) or more accurately (54) than iTRAQ approaches.
Second, blood serum was chosen compared to EDTA anticoagulated plasma. To date, there has been little or no research published comparing the differences between serum and plasma blood samples in stroke proteomics. EDTA has been shown to be the least proteolytically active of the plasma samples and has previously been shown to reveal low abundance proteins in healthy patients (55). However, the proteomic and secretomic profile of serum may be more suitable for profiling of coagulation and complement systems (56). Furthermore, in a small batch optimization analysis of these bloods, we found overall increased expression of coagulation and complement in serum compared to EDTA, as expected (57).
Third, biological pathways and mechanisms were ascertained using the gene set enrichment analysis (GSEA) statistical and bioinformatics tool. This tool compares the current protein expression data from our study to those from the Molecular Signatures Database (MSIGDB) (58). The approach of DAVID considers gene annotation clustering after statistical procedures have identified significant differences in protein expression (59). However, the GSEA approach develops an enriched gene set from the original expression data and compares them to sets from the databases on MSIGDB (58). Using pathway analysis affords greater explanatory power in highlighting relationships between gene sets and phenotypes that might go unnoticed in a comparison of individual proteins. While this is inherently a technique for genes and genomic profiling, the statistical approach of database list ranking in GSEA has been demonstrated to be applicable to mass spectrometry (MS)-based proteomics (60–62). The databases cannot yet interpret the full complexity of the molecular proteome as they do not consider protein/protein interactions and post-translational modifications. However, despite these limitations, GSEA represents one of the best tools currently available for understanding cellular pathways (63). Ultimately, we aim to show that the proteomics approach is a reliable and viable method informing our understanding of PSD and stroke pathophysiology.
Subjects and Methods
A subset of 44 stroke patients (30 females and 14 males) were obtained consecutively and prospectively from the START_PrePARE (STroke imAging pRevention and Treatment_Prediction and Prevention to Achieve optimal Recovery Endpoints) cohort, a longitudinal stroke cohort study with advanced clinical and neuroimaging data conducted in Australia (64). The participants were recruited following a first ischemic stroke and were over 18 years of age. A diagnosis of ischemic stroke was determined by an experienced neurologist using clinical assessments supplemented with computerized tomography (CT) or MRI. Further information on inclusion and exclusion criteria for this study can be found in protocol papers for START_PrePARE (64) and START_EXTEND (65). Patients included were selected consecutively from the prospectively collected START_PrePARE cohort as patients were being recruited. Patients were not excluded if they had a prior history of depression, as our focus was to identify biological factors associated with the presence of depression in stroke survivors at a particular point in time, irrespective of whether depression was present prior to or post-stroke (64). Inclusion of patients with prior history of depression was also important to improve ecological validity and permit generalization of findings to clinical populations, given prior history of depression in stroke patients.
The data presented in this manuscript were obtained primarily at 3 months post-stroke (±7 days), when clinical data and blood were collected. Ethics was approved by the Human Research Ethics Committee of Austin Hospital, Heidelberg (HREC code: H2010/03588), and relevant university and hospital sites. Further details on study design are provided in related protocol papers for START_PrePARE (64) and START_EXTEND (65).
The National Institute of Health Stroke Scale (NIHSS) was administered by a trained neurologist or health-care professional. This measure is designed to be a test of patient neurological status and correlates highly with stroke severity (66). Prior history of depression was obtained using the two-item Patient Health Questionnaire (PHQ-2) (67) at 3–7 (±1) days post-stroke and again at 12 months (±7 days) post-stroke. The PHQ-2 is a brief screening tool that employs the first two questions of the nine-item Patient Health Questionnaire (PHQ-9) (68) and shares high correlation, interrater reliability, and internal consistency with the PHQ-9 (69).
Depressive symptoms were assessed at 3 months post-stroke using the Montgomery–Åsberg Depression Rating Scale (MADRS). The MADRS is a validated measure in clinical depression research (70, 71), and was delivered using the structured interview guide (SIGMA) (72). This is a standardized interview format for the MADRS, providing the clinician with more versatility to probe the circumstances surrounding depression symptomatology compared to a self-report method. The SIGMA format is reported to have higher reliability than the self-report MADRS (72). The MADRS_SIGMA was selected as it is a standardized and validated observer-rating of depression at a point in time (72). It was administered by a health-care professional (stroke nurse, occupational therapist, or doctor) specifically trained in the administration of this assessment. Assessors used SIGMAs, detailed protocol manuals, and training videos to enhance standardization in the administration of the tool. Non-verbal supports were available to patients who were aphasic. Higher scores on the MADRS indicate more depressive symptoms.
Global cognitive impairment was screened using the Montreal Cognitive Assessment (MoCA), a validated screening measure of cognitive impairment in dementia and post-stroke recovery (73). While positive scores on the MoCA indicate optimal cognitive functioning, an adjusted cutoff score of <23 in stroke populations suggests cognitive impairment (74). The modified Rankin Scale (mRS) was used to measure functional disability. It is an interview based measure, with lower scores indicating lesser levels of observable functional impairment in daily life (75). All assessments were administered by a health-care professional, specifically trained in the administration of these assessments. Background details on age, gender, subtype of acute ischemic stroke, and thrombolysis were also obtained.
Blood Collection and Serum Separation
Blood samples were obtained by venipuncture at the 3-month follow-up assessment. All samples were collected in plastic serum-separating tube (SST) vacutainers and were allowed to clot at ambient temperate for 30 min. The tubes were then centrifuged at 1100–1300 g at room temperature, and the resulting serum was aliquoted into Eppendorf 4 × 2.0 mL tubes and immediately stored at −80°C. Upon moving blood samples from hospital sites to the central laboratory, temperature was kept at −20°C prior to transfer into a −80°C freezer.
Sample Preprocessing and Trypsination
Ten microliters of serum from each patient was first stabilized in 100 μL of 8M urea pH = 8.3 and stored at −80°C until used. For proteomic analysis, the stabilized samples were processed as follows: 20 μL protein solution was added to 90 μL of 8M urea pH = 8.3 and reduced for 5 h with 1 μL of 200 mM tris(2-carboxyethyl)phosphine (TCEP). After this, samples were alkylated for 1 h at 25°C in the dark with 4 μL of 1M iodoacetamide (IAA). In sample, digests were performed overnight (37°C) by addition of 1 μg of trypsin (Promega, Madison, WI, USA) and 900 μL of 50 mM Tris pH = 8.3, followed by a second digestion step with 1 μg trypsin and an additional incubation of 4 h at 37°C. Two hundred microliters of the digested solution were collected and dried by SpeedVac centrifugation. The digested proteins were resuspended in 100 μL of 1% (v/v) formic acid and centrifuged at 14,000 rpm for 2 min. The solid-phase extraction was performed with Empore reversed-phase extraction disks (SDB-XC reversed-phase material, 3M) according to Ishihama et al. (76) with the following modifications: the membrane was conditioned with 50 μL of 80% (v/v) acetonitrile, 0.1% (w/v) trifluroacetic acid, and then washed with 50 μL of 0.1% trifluroacetic acid before the tryptic peptides were bound to the membrane. The bound peptides were eluted by 50 μL 80% (v/v) acetonitrile, 0.1% (w/v) trifluroacetic acid, and dried in a SpeedVac centrifuge.
Tryptic peptides reconstituted in 0.1% formic acid and 2% acetonitrile (buffer A) were analyzed by LC-ESI–MS/MS on a LTQ-Orbitrap Elite (Thermo-Fisher Scientific). Peptides were loaded onto a trap column (C18 PepMap 100 μm i.d. × 2 cm trapping column, Thermo-Fisher Scientific) at 5 μL/min for 6 min before switching the precolumn in line with the analytical column (Easy-Spray 75 μm i.d. × 50 cm, Thermo-Fisher Scientific). The separation of peptides was performed at 250 nL/min using a linear acetonitrile gradient of buffer A and buffer B (0.1% formic acid and 80% acetonitrile), starting from 5% buffer B to 60% over 300 min. This final separation step is equivalent to fractionation and was conducted in order to avoid potential biases and increased sample variability due to depletion. Although it is impossible to eliminate dynamic range issues in serum samples, this technique greatly increases the dynamic range detectable in our experiment.
Data were collected in data-dependent acquisition mode using m/z 300–1500 as MS scan range; CID MS/MS spectra were collected for the 20 most intense ions. Dynamic exclusion parameters were set as follows: repeat count 1, duration 90 s, and the exclusion list size was set at 500 with early expiration disabled. Other instrument parameters for the Orbitrap were the following: MS scan at 120,000 resolution, maximum injection time 150 ms, AGC target 1 × 106, and CID at 35% energy for a maximum injection time of 150 ms with AGT target of 5000. The Orbitrap Elite was operated in dual analyzer mode with the Orbitrap analyzer being used for MS and the linear trap being used for MS/MS. This procedure was performed on two technical replicates. The samples were then analyzed using the in-house Mascot server for protein identification. MaxQuant (Max-Planck Institute for Biochemistry, Martinsried, Germany) was used to obtain the relative quantification of identified proteins in the samples. Relative intensity or quantification is a measurement of peak height in a single sample that is compared to the same measurement in other samples. The absolute concentrations of the proteins in the sample are not known and require other methodologies to obtain, and thus, relative protein expression data cannot be generalized to other protein assays. Absolute quantification is possible on a mass spectrometer, but requires prior knowledge of target proteins and extensive methodologies.
Protein Identification and Label-Free Quantitation
Identification and LFQ of obtained spectra across all 44 samples was performed using MaxQuant version 18.104.22.168 to obtain identified proteins (77). All raw data and complete details of MaxQuant parameters and result files have been deposited in ProteomeXChange and are available with accession number PDX003494. Identification of peptides and proteins was performed internally by MaxQuant using the Andromeda (77) search engine to search against all reviewed and unreviewed human proteins in the Uniprot database (August 2013; 133,798 entries in total). Common contaminants and decoys (reversed sequences) were included automatically by Andromeda. Prior to searching, MS/MS spectra were filtered according to MaxQuant default settings for ion trap MS/MS spectra by retaining only the top eight peaks per 100 Da. Main search, precursor mass tolerance was set to 4.5 ppm, and MS/MS tolerance to 0.5 Da. Carbamidomethylation of cysteines was set as a fixed modification, and N-term acetylation and oxidation of methionine were included as variable modifications. Up to two missed cleavages were allowed, and peptides were required to be at least seven amino acids in length. False discovery rate (FDR) cutoffs for both peptides and proteins in the database search were set to 1%. Both unique and razor peptides were used for quantitation with a minimum of two peptides including at least one unique peptide required to calculate a protein quantitative value. The “match between runs” setting in MaxQuant was used to transfer peptide identifications from one run to another on the basis of matching retention time and mass-to-charge ratio.
Data Preprocessing and Analysis
The initial output from MaxQuant consisted of 515 protein groups. After removal of contaminants, the list was shortened to 475 protein groups. The LFQ signal intensity was log2 normalized to account for naturally skewed intensity values (78) and averaged over technical replicates. The discovery approach here employs relative quantification instead of absolute, and as such, it was better suited to examine a collection of genes/proteins as ontologies of biological structures or pathway processes to extract biological significance from the proteomics data and understand the systems involved in depression symptoms post-stroke (63). As relative protein expression is limited to comparison within the study or similar MS-based studies, pathway analysis provides greater explanatory power than traditional statistical biomarker approaches. Thus, the resulting protein group expression data were reduced to gene names then paired with the continuous MADRS scores to prepare for GSEA (Broad Institute, MIT). GSEA allows for a robust comparison of continuous phenotypes to gene expression with the selection of the “Pearson” metric. All analyses in GSEA were conducted with the GENE_SYMBOL chip and default number of permutations (1000) on full gene sets from Hallmark, Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), Biocarta, and Reactome Positional and Immunologic Signatures, acquired from MSIGDB v5.00. The Hallmark database is recommended as an entry to GSEA analysis, as it collects gene sets that represent well-defined biological states and processes with expression scores computed from many existing gene sets to reduce noise and redundancy, acting as the searchable “meta-analysis” of gene sets. GO is the earliest but most up-to-date functional gene annotation database and encompasses the largest variety of annotations under three headings: biological processes, molecular function, and cellular component (79). KEGG (80), Biocarta (81), and Reactome (82) gene sets are curated from external databases, with each database representing different approaches to compiling, such as genome sequencing, microarray profiling, and computational methods, to build complex networks. Positional gene sets represent the locations of genes on chromosomes and cytogenic bands (83), while the Immunologic Signatures database is a collection of immune responses curated from separate microarray studies. Included gene sets were no smaller than 15 and no larger than 500. As per recommendation by the GSEA program for discovery experiments, a FDR of 25% was deemed acceptable for statistical significance of enriched gene sets, where there is a 75% chance of rejecting a false positive (58). A nominal p (nom-p) value based on the statistical significance of each individual database is also included for reference. This statistic is not adjusted for multiple testing, whereas the FDR value is.
An enrichment score is a Kolmogorov–Smirnov-like ranking statistic that reflects the degree to which a set is overrepresented at the top (positive enrichment) or bottom (negative enrichment) of the list when compared to another list (58). For example, when adjusted for correlation with increasing score on a clinical scale, such as the MADRS, positive enrichment suggests biological upregulation of a given gene list, while negative enrichment suggests downregulation of a gene list with increasing presentation of depression symptoms. The normalized enrichment score (NES) is the primary statistic used for evaluating and comparing gene sets and is understood as (58):
This analysis yields the “leading edge subset” that refers to the cluster of proteins that contribute most to the enrichment score and can be interpreted as the genes that are most likely to affect change in complex pathway function.
To minimize the possible confound of prior history of depression on the analysis, the appropriate T-test or non-parametric (Mann–Whitney) comparisons were conducted and compared with group comparisons with Monte Carlo simulations of the MADRS scores in subgroups with and without prior history of depression.
Patient Demographics and Clinical Status
The average age of ischemic stroke patients in the sample was 63.62 years, range 34–87 (SD = 13.52), within the range of previous epidemiological findings, although younger than the mean age of 74.40 years (84). Thirty were males and 14 females. Based on the TOAST classification (85), 11 patients had large artery atherosclerosis, 7 had cardioembolism, and 11 had small vessel occlusion. The remaining 15 patients were unclassified. Nine patients were administered thrombolytic tPA within the 4.5-h time window and two were unknown. Patient demographics and NIHSS score at time of recruitment (i.e., baseline M = 1.39, SD = 1.26 days post-stroke), and clinical characteristics at 3–7 days and at 3 months post-stroke are presented in Table 1.
The average MADRS score for the sample was 6.40 (SD = 7.42, 95% CI = [4.17, 8.63], range 0–26). Thirteen (29.55%) patients had a score of >6 but <18, indicating that this sample mostly had mild depressive symptoms. Only 4 (9.09%) patients had scores within the severe symptoms range. However, a person scoring 12 in major depressive disorder or 8 post-stroke on the MADRS can be considered for treatment (86). The sample presented with mild stroke severity at 3 months, with a median NIHSS score of 0 and range 0–5. Similarly, scores on the MoCA and mRS suggest that this sample had relatively few cognitive problems and had recovered in daily activities at 3 months from the stroke. A correlation analysis of the MADRS with these clinical measures and patient age revealed no significant results at 3 months post-stroke.
Screening on the PHQ-2 revealed that 11 of the 44 patients had a prior history of depression. As it was possible that patients with a pre-stroke history of depression may be predisposed to more severe PSD symptoms, a Mann–Whitney comparison was conducted to test heterogeneity of MADRS scores between those with and without prior history of depression. Preliminary analysis suggested that the distributions of both groups were similar, and thus, the combined median metric is used (87). After running a Monte Carlo simulation for 1,000,000 samples at 99% CI to account for low sample sizes, it was found that patients with a previous history of depression (Mdn = 5.00) did not significantly differ from patients without any history (Mdn = 2.50), U = 165.00, z = −0.591, Monte Carlo p = 0.565 [99% CI = (0.563, 0.566)]. Following this validation, it was possible to continue with the available range of MADRS scores without having to categorize based on prior history of depression.
Gene Set Enrichment Analysis
Peptides (n = 475) from the proteomic analysis were analyzed to reveal underlying molecular pathways and gene ontologies associated with depressive symptoms. Five different data sets (Hallmark, GO, KEGG, Biocarta, and Reactome) were interrogated to identify enriched data sets. Multiple databases were used in our discovery approach to not only explore a comprehensive range of possible pathways (given the fact that each database is constructed differently) but also to identify commonalities across databases, thus strengthening the robustness and generalizability of our findings.
Of the gene set databases that were entered in GSEA (Table 2), Positional and Immunologic Signatures did not return enriched sets. With the exception of GO, all databases showed significantly enriched gene sets pertaining to complement cascade activation or general immune upregulation/downregulation. All sets were negatively enriched when compared with increasing MADRS score, suggesting significant protein downregulation with increasing level of depressive symptoms. There was also no clear enrichment pattern that would associate individual genes with MADRS scores. However, there are slight variations in set size and top gene contributors that reflect the differences in the compiling of the databases itself. Thus, although the statistics and bioinformatics analysis show that these pathways are associated with depression, biological interpretation is dependent upon further information from the individual databases. Positional gene sets collect data about chromosomal cytogenic band positioning of the genes involved, while Immunological Signatures collect published examples of specific immune activity against an immune challenge such as dendritic cell activity in human immunodeficiency virus.
Table 2. Size of gene set databases searched in GSEA and corresponding number of enriched sets adjusted for FDR.
As the Hallmark gene sets are computed from a collection of similar biological processes, this set provides strong support for the negative enrichment of this set and associated gene expressions (Table 3). A negatively skewed NES for both coagulation and complement Hallmark gene sets correlated against increasing MADRS scores indicates that there is an overexpression of downregulated genes in these sets. These genes can be seen in the leading edge subset (highlighted in bold), while the other genes comprise the structure of the set and are important in affirming the construct validity of the obtained sets, which do not significantly contribute to the overall magnitude and direction of the enrichment score.
Table 3. Significantly enriched gene sets and individual gene contributors from the hallmark database.
Kyoto Encyclopedia of Genes
Coagulation and complement cascades and systemic lupus erythematous (SLE) were implicated in the KEGG database (Table 4). In this database, the SLE set is defined by the antigen-activated complement pathway, demonstrating the characteristic molecular cascades of immune dysfunction. Even though SLE is a significant risk factor of ischemic stroke (88), as seen from the leading edge (Table 4), it is more likely that this set was merely expressing statistical enrichment of a gene set that is largely comprised of complement. Central to all complement pathways, complement component 3 (C3) is one of the main activators and mediators of this pathway and was found to contribute highly to the negative enrichment score.
The set size of Biocarta complement was only 17, with the majority of genes contributing to the various complement pathways and formation of the membrane attack complex (MAC) that is involved in attacking target cells. This gene list suggests that both mannan-binding lectin 2 (MBL2) and complement component 1 (C1Q), precursors of both the lectin and classical pathways, are downregulated in association with depressive symptoms (Table 5). C3 is central to the entire cascade as its activation is required for both lectin and classical pathways and attraction of the adaptive immune system.
Table 5. Significantly enriched gene sets and individual gene contributors from the Biocarta database.
The Reactome database implicated three gene sets with significant negative enrichment in association with MADRS scores, all with similar individual gene contributions (Table 6). This database consists of a large pathway map of which complement cascades are categorized under the innate immune system. The results for the complement pathways were similar to those obtained previously; both classical and lectin pathways are downregulated, and MAC activation genes, such as C9 and C6, are underrepresented. These genes are essential to the overall function of the complement system (Figure 1).
Table 6. Significantly enriched gene sets and individual gene contributors from the Reactome database.
Figure 1. The reactome pathway diagram, demonstrating the primary signaling and functional molecular pathways of the complement system on a target cell. In this diagram, multiple activation pathways are shown, eventually effecting cell lysis and breakdown of the target cell. This also shows secondary signaling pathways for C3a signaling for pro-inflammation and cytokine activity with C5a-promoting chemotaxic recruitment of inflammatory cells. Retrieved from: http://www.reactome.org/figures/Pathway_Illustrations/Complement_Cascade_72_oicr.png; reused under the Creative Commons license.
Depression after stroke is a common occurrence, negatively impacting functional outcome, response to rehabilitation, and quality of life. Utilizing a discovery proteomics approach, this study was aimed to find biological associates of depressive symptoms (MADRS scores) post-stroke in human serum. This was achieved by employing a label-free proteomics workflow, utilizing an LTQ Orbitrap with high dynamic range to maximize the range of protein identification. Statistical and bioinformatics analysis was completed on the GSEA platform on a wide range of curated gene set databases. This approach revealed consensus among well-developed databases for decreased gene expression that was associated with the complement pathways and partial support for precursors of the coagulation pathway.
The most important findings of this study are downregulation of coagulation and complement cascades in serum bloods from stroke patients with increasing level of MADRS-defined depressive symptoms, as depicted by Hallmark, KEGG, Reactome, and Biocarta databases. Most PSD studies examine a severe depression phenotype. However, we have shown here that even patients who have mild depressive symptoms and have recovered well can still exhibit biological changes that can be characterized by a proteomics approach. Furthermore, even though the analytical approach was different, our findings provide additional support for those obtained by Zhan et al. (47). Our findings also support the dysregulation of the coagulation and complement cascade pathway that was identified in Zhan et al. (47), in serum as opposed to EDTA blood samples. The contributing differences in individual gene expression are likely due to a number of factors, such as stroke timeline, study design, and statistical approaches, but most importantly, the use of serum contrasted with EDTA for blood preparation. Considering these methodological differences between the two studies, finding similar results only provides further confirmation for the proteomic profile of stroke and strengthens biological understanding of PSD.
From a molecular standpoint, the primary function of the complement system can be summarized as formation of the MAC on target cells, while secondary functions include signaling for toll-like receptor 4 (TLR4)-mediated inflammatory response (89), leading to increased pro-inflammatory cytokine levels (90). Complement has been traditionally viewed as a versatile system, conforming to three distinct activation pathways, leading to similar outcomes (91). Of the leading edge subsets of the databases searched, it is clear that a majority of the molecular determinants of the complement system are downregulated in relation to mild depressive symptoms. Of note, complement component 1q (C1q) and complement component 3 (C3) are both integral to the activation of this system, while downstream products that contribute to inflammatory signaling and formation on the MAC (Figure 1) are all downregulated. This may indicate perturbation in brain homeostasis following stroke, characterized by an ongoing state of lowered resistance to oxidative stress (92) and immune-related changes (93) in maintaining depression symptoms post-stroke. This is further implied by the Reactome database analysis, demonstrating negative enrichment of calcium-binding protein A12 (S100A12) and lipopolysaccharide-binding protein (LBP), attributed to innate immune functioning. S100A12 is an established cytoskeletal protein involved in the signaling of neutrophil response and a candidate biomarker of inflammation (94).
Coagulation is a key process in vascular diseases, especially in ischemic stroke where coagulation factors have been identified as pre-stroke risk factors (95). This process is heavily involved in the stroke itself and treatment with tPA (96, 97). Coagulation is an ongoing process in serum and deficiencies of several proteins, such as protein Z (PROZ) (98), fibronectin 1 (FN1) (97), and gelsolin (GSN) (99), have been identified to be involved in cerebrovascular and cardiovascular diseases. The overlap of genes between coagulation and complement sets is expected as the coagulation processes have been identified to cleave into the central components of the complement system (100). The inclusion of both complement and coagulation cascades into one gene set in KEGG compared to the presentation of two separate sets in Hallmark illustrates a previous debate as to whether the pathways should be considered separately (101, 102). However, further studies have shown that while complement deficiencies alone do not increase bleeding frequency and coagulation deficiencies alone do not impair immune responses, the functions of the two pathways are linked and indeed in the case of innate immunity (103, 104). It has been suggested that coagulation and complement are involved in a feedback loop, where complement activation increases platelet activation area on the target cell that in turn augments complement activity (105). Coagulation has independently been shown to be involved in the process of cell death and is associated with activation of the kinin–kallikrein pathway and facilitates defensive inflammatory responses by enhancing leukocyte activity (106). Finally, binding platelets also have the ability to enhance neutrophil activity via TLR4 signaling, a process that is vital to phagocytosis (105).
Although a robust association between relative presence of depressive symptoms and downregulation of complement and coagulation pathway was found in our necessarily small sample of mild stroke survivors living in the community, our findings need to be interpreted with care, especially as some patients had a prior history of depression. We did use a group comparison utilizing Monte Carlo simulations that indicated that prior depression did not significantly impact our findings. Due to the discovery approach and the patient’s characteristics of this sample, the generalizability of the findings is limited to survivors with mild stroke severity.
Complement Post-Stroke: From Acute Cell Death to Immunodepression and Depression
Cell death in the brain is the inevitable consequence of any stroke damage. In stroke, the first wave of cell death occurs as a result of hypoxia. This releases damage-associated molecular patterns (DAMPs) that can begin and perpetuate apoptotic and necrotic cell death cascades in neighboring neurons. Such molecules can include, but are not limited to, intracellular adenosine triphosphate and uridine triphosphate that has been leaked into the extracellular space, nitric oxide, heat shock proteins, S100 proteins, extracellular calcium ion levels, and cytokines (107, 108). The presence of DAMPs on the central side of the blood–brain barrier (BBB) are also associated with endothelial damage, which further contributes to increased permeability and infiltration of immune cells from the periphery (109).
Transcriptomic research has not yet identified pathways involved in the timeline of PSD pathogenesis in bloods; however, an overview of studies examining the whole blood ribonucleic acid (RNA) profiles of different stroke subtypes, including TIA in <1 week, has yielded interesting results (110). Functional analysis of the genes in these studies has shown that immune and homeostasis pathway expression can differentiate between cardioembolic and atherosclerosis stroke (111), and Gene Ontology clustering in cardioembolic ischemic stroke suggests gene expression indicative of cell death, lipid metabolism, and metal ion transport. Thus, here we propose that the proteomically detectable state of the peripheral immune system may be indicative of unresolved ischemic damage in the central parenchyma, resulting in disruptions to inflammatory, metabolic, and homeostatic balance that is related to the transcriptomic profile in the early stages of stroke damage (110). While there is no neurobiological evidence presented in our study, current knowledge of the etiology of PSD suggests that persistent neuroinflammation, driven by an increase in pro-inflammatory but also a decrease in anti-inflammatory cytokine signaling, is responsible for depressive symptoms post-stroke (22, 112). Indeed, it is possible that peripheral immunodepression may be caused by central and upstream molecular cascades that include cross BBB cytokine signaling (113) or bioavailability of immunoglobulin antibodies (114) and leukocyte immune cells (115). Thus, PSD can be considered as a natural sequelae of incomplete recovery from stroke that can be further exacerbated by anxiety from psychosocial issues (8).
We have shown that two innate immune pathways in peripheral bloods, complement and coagulation, trend toward downregulation at the 3-month phase in correlation with mild symptoms of depression. The complement system is largely involved in an array of normal and immunoregulatory functions, although its role in immunodepression post-stroke has not been explored (Table 7). As a whole however, immunodepression after stroke is a well-documented and a natural consequence of ischemia, but is poorly understood in functional terms (116). Initially, it appears that immunodepression is counterproductive after stroke as it may increase the chances of commonly reported secondary infections such as urinary tract infection and pneumonia (117, 118). Additionally, immune suppression is thought to be an adaptive response to central inflammation as an autoimmune response against the brain would be detrimental to recovery outcomes and possibly exacerbate damage (116, 119). There are many factors that may maintain poor immune responses and recovery in animal models such as induced psychosocial stress (120), abnormal BBB permeability (121), and possible ongoing antigen-related responses that have yet to be fully characterized (122, 123).
Table 7. A brief overview of the complement system in both central and peripheral functioning related to stroke risk, the stroke, and post-stroke outcomes, as understood from both human and animal studies.
There has been little research exploring the role of complement and indeed immunity post-stroke in relation to depression. The results here suggest that there is ongoing immunosuppression in the periphery with mild depressive symptoms, even at 3 months post-stroke. Previous research into antigen-presenting cells, such as macrophages and dendritic cells, has established that their expression is elevated centrally, but downregulated peripherally post-stroke (138, 139). A trend toward immune recovery would stipulate that a balance has been reached in these cell levels for both central and periphery. Indeed, it has also been theorized that this discrepancy accounts for the recruitment of peripheral dendritic cells into the brain to maintain a central inflammatory and immune response (140). Thus, given these lines of evidence, it is possible that our current findings reflect an ongoing sub-acute state of pro-inflammation that has not yet transitioned to anti-inflammation (109), manifesting behaviorally as depression symptoms and molecularly as peripheral immunodepression of the complement and coagulation systems.
Limitations and Future Studies
The patients in this study were not severely depressed and had recovered well after their stroke, with some patients recording a history of pre-stroke depression. Furthermore, there was no age-matched control group present in this study. The current study was a pilot to investigate the feasibility of the proteomics approach in stroke. While we found a significant association even in a sample with mild depressive symptoms and mild stroke severity, future studies that employ larger sample sizes with depressive symptoms and neurological severity that represent the range commonly experienced are recommended. Future studies may consider basic or laboratory assessment of the immunological condition of patients post-stroke as well as the relationship between clinical and/or psychosocial factors on PSD. We used a discovery approach where the findings can stand alone for the stroke cohort, although it is recognized that this approach is not as robust as with age-matched controls. Therefore, we recommend future comparison with age-matched non-stroke controls, with and without depression, to enable a more comprehensive interpretation of findings. In addition, future studies may compare the proteomes of other biofluids, such as CSF and urine, of the same patients to develop a better understanding of the compartmentalization or relationships of BBB and kidney physiology post-stroke. Complement and coagulation functioning can also be assayed in blood in the traditional hematological laboratory setting. From a technical perspective, further studies could optimize the comparison of ionic chelation properties of anticoagulants as this is not well understood in whole blood samples. It is also feasible to employ immunodepletion or different fractionation techniques on the mass spectrometer in the sample preparation stages to resolve dynamic range and peak detection issues for low abundance compounds.
This study examined the serum proteomic profile of stroke survivors at 3 months post-stroke using a label-free approach. The findings here and in Zhan et al. (47) are complementary and provide a basis for further research into blood proteins recently identified to be involved in the pathophysiology of PSD and possibly in other cerebrovascular diseases with comorbid anxiety. Analysis by GSEA on various databases has revealed enriched gene sets that are identifiable as complement and innate immune processes. As all of these gene sets were negatively enriched when correlated with increasing depressive symptoms of ischemic stroke survivors, this was interpreted as peripheral immunodepression, indicative of unresolved ongoing inflammatory processes in the brain. This interpretation is supported by the substantial body of literature that has linked PSD etiology to overactive immunologic processes, leading to increased inflammatory processes in both peripheral and central compartments. While these findings add to the growing body of evidence for differentially expressed proteins in PSD, more research is needed to characterize their molecular processes and how their expressions may change as result of the stroke and in development of PSD.
Patient, family member, or legally responsible person, depending on local ethics requirements, have given informed consent for participation in the START_PrePARE study.
Main author: VN. Senior authors: SC and LC. Cohort Project (START and PrePARE) conception and design: LC, GD, SD, DH, SC, LM, and HM. Data analysis application and interpretation: VN, SC, LG, TT, and DH. Proteomics methodology and data acquisition: SC, LC, DH, PF, and IC. Patient testing, project co-ordination: LC, HM and TT. All of the authors listed have read and approved the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors acknowledge and thank the participants, lead investigators and study co-ordinators from the START recruiting sites. In particular we acknowledge and thank those from the START-PrePARE recruiting sites: Austin Hospital (Helen Dewey), Epworth Hospital (Richard Gerraty), Monash Medical Centre (Henry Ma), Royal Melbourne Hospital (Bruce C. V. Campbell), and Western Hospital (Tissa Wijeratne). We also thank the Commonwealth Scientific and Industrial Research Organization (CSIRO) management and research teams involved in this research, in particular the CSIRO data analytics team (Simon McBride, Karen Harrap, Chris Stanbridge). Finally, we acknowledge and thank the management and research team of the START program of research (www.START.csiro.au).
We acknowledge financial support for conduct of the research from the Commonwealth Scientific and Industrial Research Organization (CSIRO) of Australia, Flagship Collaboration Fund through the Preventative Health Flagship; data collection and project planning by the Stroke Imaging, Prevention and Treatment research team (www.START.csiro.au); support for analysis and write up and/or researchers from the James S. McDonnell Foundation 21st Century Science Initiative in Cognitive Rehabilitation – Collaborative Award (# 220020413); Victorian Government’s Operational Infrastructure Support Program; an Australian Research Council Future Fellowship awarded to LC (#FT0992299); and a La Trobe University Post Graduate Scholarship supported by the Understanding Diseases Research Focus Group at La Trobe University.
ApoA4, ApoC2, APOE, apolipoproteins A4, C2, E; BBB, blood–brain barrier; C6, C3, C3a, C5a, complement component 6, 3, 3 a, 5a; EDTA, ethylenediaminetetraacetic acid; ESI, electrospray ionization; FDR, false discovery rate; GO, Gene Ontology; GSEA, gene set enrichment analysis; GSN, gelsolin; ITRAQ, isobaric tagging for relative and absolute quantification; KEGG, Kyoto Encyclopedia of Genes and Genomes; LC, liquid chromatography; LFQ, label-free quantification; MADRS, Montgomery–Åsberg Depression Rating Scale; MAC, membrane attack complex; MS, MS/MS, mass spectrometry, tandem mass spectrometry; MBL2, mannan-binding lectin 2; MoCA, Montreal Cognitive Assessment; mRS, modified Rankin Scale; NES, normalized enrichment score; NIHSS, National Institute of Health Stroke Scale; PHQ-9, PHQ-2, Patient Health Questionnaire 9, 2; PrePARE, Prediction and Prevention to Achieve optimal Recovery Endpoints; PSD, post-stroke depression; RNA, ribonucleic acid; S100A12, S100 calcium-binding protein A12; SIGMA, structured interview guide for the Montgomery–Åsberg Depression Rating Scale; SLE, systemic lupus erythematous; START, STroke imAging, pRevention and Treatment longitudinal stroke cohort study; TF, transferring; TIA, transient ischemic attack; TLR4, toll-like receptor 4; TOAST, trial of ORG 10172 in acute treatment; tPA, tissue plasminogen activator.
6. Ahn D-H, Lee Y-J, Jeong J-H, Kim Y-R, Park J-B. The effect of post-stroke depression on rehabilitation outcome and the impact of caregiver type as a factor of post-stroke depression. Ann Rehabil Med (2015) 39(1):74–80. doi:10.5535/arm.2015.39.1.74
7. Schmid AA, Kroenke K, Hendrie HC, Bakas T, Sutherland JM, Williams LS. Poststroke depression and treatment effects on functional outcomes. Neurology (2011) 76(11):1000–5. doi:10.1212/WNL.0b013e318210435e
8. Pascoe MC, Crewther SG, Carey LM, Crewther DP. Inflammation and depression: why poststroke depression may be the norm and not the exception. Int J Stroke (2011) 6(2):128–35. doi:10.1111/j.1747-4949.2010.00565.x
10. Vataja R, Leppävuori A, Pohjasvaara T, Mäntylä R, Aronen HJ, Salonen O, et al. Poststroke depression and lesion location revisited. J Neuropsychiatry Clin Neurosci (2004) 16(2):156–62. doi:10.1176/jnp.16.2.156
11. Tang W, Chen Y, Lu J, Chu WC, Mok V, Ungvari GS, et al. White matter hyperintensities in post-stroke depression: a case control study. J Neurol Neurosurg Psychiatry (2010) 81(12):1312–5. doi:10.1136/jnnp.2009.203141
12. Mitra J, Shen K-K, Ghose S, Bourgeat P, Fripp J, Salvado O, et al. Predicting poststroke depression from brain connectivity. In: O’Donnell L, Nedjati-Gilani G, Rathi Y, Reisert M, Schneider T, editors. Computational Diffusion MRI. Boston: Springer International Publishing (2014). p. 89–99.
14. Loubinoux I, Kronenberg G, Endres M, Schumann-Bard P, Freret T, Filipkowski RK, et al. Post-stroke depression: mechanisms, translation and therapy. J Cell Mol Med (2012) 16(9):1961–9. doi:10.1111/j.1582-4934.2012.01555.x
15. Robinson RG, Shoemaker WJ, Schlumpf M, Valk T, Bloom FE. Effect of experimental cerebral infarction in rat brain on catecholamine and behaviour. Nature (1975) 255(5506):332–4. doi:10.1038/255332a0
16. Perini F, Morra M, Alecci M, Galloni E, Marchi M, Toso V. Temporal profile of serum anti-inflammatory and pro-inflammatory interleukins in acute ischemic stroke patients. Neurol Sci (2001) 22(4):289–96. doi:10.1007/s10072-001-8170-y
17. Vu NQ, Aizenstein HJ. Depression in the elderly: brain correlates, neuropsychological findings, and role of vascular lesion load. Curr Opin Neurol (2013) 26(6):656–61. doi:10.1097/wco.0000000000000028
18. Broomfield NM, Quinn TJ, Abdul-Rahim AH, Walters MR, Evans JJ. Depression and anxiety symptoms post-stroke/TIA: prevalence and associations in cross-sectional data from a regional stroke registry. BMC Neurol (2014) 14:198. doi:10.1186/s12883-014-0198-8
19. El Husseini N, Goldstein LB, Peterson ED, Zhao X, Pan W, Olson DM, et al. Depression and antidepressant use after stroke and transient ischemic attack. Stroke (2012) 43(6):1609–16. doi:10.1161/strokeaha.111.643130
20. Luijendijk HJ, Stricker BHC, Wieberdink RG, Koudstaal PJ, Hofman A, Breteler MM, et al. Transient ischemic attack and incident depression. Stroke (2011) 42(7):1857–61. doi:10.1161/strokeaha.110.604405
21. Rao R, Jackson S, Howard R. Depression in older people with mild stroke, carotid stenosis and peripheral vascular disease: a comparison with healthy controls. Int J Geriatr Psychiatry (2001) 16(2):175–83. doi:10.1002/1099-1166(200102)16:2<175:AID-GPS298>3.0.CO;2-0
22. Spalletta G, Bossu P, Ciaramella A, Bria P, Caltagirone C, Robinson RG. The etiology of poststroke depression: a review of the literature and a new hypothesis involving inflammatory cytokines. Mol Psychiatry (2006) 11(11):984–91. doi:10.1038/sj.mp.4001879
23. Vogelgesang A, May VEL, Grunwald U, Bakkeboe M, Langner S, Wallaschofski H, et al. Functional status of peripheral blood T-cells in ischemic stroke patients. PLoS One (2010) 5(1):e8718. doi:10.1371/journal.pone.0008718
24. Hama S, Yamashita H, Yamawaki S, Kurisu K. Post-stroke depression and apathy: interactions between functional recovery, lesion location, and emotional response. Psychogeriatrics (2011) 11(1):68–76. doi:10.1111/j.1479-8301.2011.00358.x
27. Kim J-M, Stewart R, Kim S-W, Shin I-S, Kim J-T, Park M-S, et al. Associations of cytokine gene polymorphisms with post-stroke depression. World J Biol Psychiatry (2011) 13(8):579–87. doi:10.3109/15622975.2011.588247
28. Noonan K, Carey LM, Crewther SG. Meta-analyses indicate associations between neuroendocrine activation, deactivation in neurotrophic and neuroimaging markers in depression after stroke. J Stroke Cerebrovasc Dis (2013) 22(7):e124–35. doi:10.1016/j.jstrokecerebrovasdis.2012.09.008
30. Worthmann H, Tryc AB, Goldbecker A, Ma YT, Tountopoulou A, Hahn A, et al. The temporal profile of inflammatory markers and mediators in blood after acute ischemic stroke differs depending on stroke outcome. Cerebrovasc Dis (2010) 30(1):85–92. doi:10.1159/000314624
32. Leonard B, Maes M. Mechanistic explanations how cell-mediated immune activation, inflammation and oxidative and nitrosative stress pathways and their sequels and concomitants play a role in the pathophysiology of unipolar depression. Neurosci Biobehav Rev (2012) 36(2):764–85. doi:10.1016/j.neubiorev.2011.12.005
33. Maes M. Depression is an inflammatory disease, but cell-mediated immune activation is the key component of depression. Prog Neuropsychopharmacol Biol Psychiatry (2011) 35(3):664–75. doi:10.1016/j.pnpbp.2010.06.014
34. Zunszain PA, Anacker C, Cattaneo A, Carvalho LA, Pariante CM. Glucocorticoids, cytokines and brain abnormalities in depression. Prog Neuropsychopharmacol Biol Psychiatry (2011) 35(3):722–9. doi:10.1016/j.pnpbp.2010.04.011
35. Wetie AN, Woods A, Darie C. Mass spectrometric analysis of post-translational modifications (PTMs) and protein-protein interactions (PPIs). In: Woods AG, Darie CC, editors. Advancements of Mass Spectrometry in Biomedical Research (Vol. 806), Cham: Springer International Publishing (2014). p. 205–35.
37. Koomen JM, Haura EB, Bepler G, Sutphen R, Remily-Wood ER, Benson K, et al. Proteomic contributions to personalized cancer care. Mol Cell Proteomics (2008) 7(10):1780–94. doi:10.1074/mcp.R800002-MCP200
40. Kumar Y, Liang C, Bo Z, Rajapakse JC, Ooi EE, Tannenbaum SR. Serum proteome and cytokine analysis in a longitudinal cohort of adults with primary dengue infection reveals predictive markers of DHF. PLoS Negl Trop Dis (2012) 6(11):e1887. doi:10.1371/journal.pntd.0001887
41. Wildsmith KR, Schauer SP, Smith AM, Arnott D, Zhu Y, Haznedar J, et al. Identification of longitudinally dynamic biomarkers in Alzheimer’s disease cerebrospinal fluid by targeted proteomics. Mol Neurodegener (2014) 9(1):1–14. doi:10.1186/1750-1326-9-22
43. Lopez MF, Sarracino DA, Prakash A, Athanas M, Krastins B, Rezai T, et al. Discrimination of ischemic and hemorrhagic strokes using a multiplexed, mass spectrometry-based assay for serum apolipoproteins coupled to multi-marker ROC algorithm. Proteomics Clin Appl (2012) 6(3–4):190–200. doi:10.1002/prca.201100041
44. Sharma R, Gowda H, Chavan S, Advani J, Kelkar D, Kumar GSS, et al. Proteomic signature of endothelial dysfunction identified in the serum of acute ischemic stroke patients by the iTRAQ-based LC-MS approach. J Proteome Res (2015) 14(6):2466–79. doi:10.1021/pr501324n
45. Ning M, Sarracino DA, Buonanno FS, Krastins B, Chou S, McMullin D, et al. Proteomic protease substrate profiling of tPA treatment in acute ischemic stroke patients: a step toward individualizing thrombolytic therapy at the bedside. Transl Stroke Res (2010) 1(4):268–75. doi:10.1007/s12975-010-0047-z
46. Pan S, Zhan X, Su X, Guo L, Lv L, Su B. Proteomic analysis of serum proteins in acute ischemic stroke patients treated with acupuncture. Exp Biol Med (2011) 236(3):325–33. doi:10.1258/ebm.2011.010041
47. Zhan Y, Yang YT, You HM, Cao D, Liu CY, Zhou CJ, et al. Plasma-based proteomics reveals lipid metabolic and immunoregulatory dysregulation in post-stroke depression. Eur Psychiatry (2014) 29(5):307–15. doi:10.1016/j.eurpsy.2014.03.004
48. Pottiez G, Wiederin J, Fox HS, Ciborowski P. Comparison of 4-plex to 8-plex iTRAQ quantitative measurements of proteins in human plasma samples. J Proteome Res (2012) 11(7):3774–81. doi:10.1021/pr300414z
49. Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, et al. Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol Cell Proteomics (2004) 3(12):1154–69. doi:10.1074/mcp.M400129-MCP200
52. Neilson KA, Ali NA, Muralidharan S, Mirzaei M, Mariani M, Assadourian G, et al. Less label, more free: approaches in label-free quantitative mass spectrometry. Proteomics (2011) 11(4):535–53. doi:10.1002/pmic.201000553
53. Wang H, Alvarez S, Hicks LM. Comprehensive comparison of iTRAQ and label-free LC-based quantitative proteomics approaches using two Chlamydomonas reinhardtii strains of interest for biofuels engineering. J Proteome Res (2012) 11(1):487–501. doi:10.1021/pr2008225
54. Trinh HV, Grossmann J, Gehrig P, Roschitzki B, Schlapbach R, Greber UF, et al. iTRAQ-based and label-free proteomics approaches for studies of human adenovirus infections. Int J Proteomics (2013) 2013:16. doi:10.1155/2013/581862
57. Yi J, Kim C, Gelfand CA. Inhibition of intrinsic proteolytic activities moderates preanalytical variability and instability of human plasma. J Proteome Res (2007) 6(5):1768–81. doi:10.1021/pr060550h
58. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A (2005) 102(43):15545–50. doi:10.1073/pnas.0506580102
62. Sridhar S, Schembri F, Zeskind J, Shah V, Gustafson AM, Steiling K, et al. Smoking-induced gene expression changes in the bronchial airway are reflected in nasal and buccal epithelium. BMC Genomics (2008) 9(1):259. doi:10.1186/1471-2164-9-259
64. Carey LM, Crewther S, Salvado O, Lindén T, Connelly A, Wilson W, et al. STroke imAging pRevention and treatment (START): a longitudinal stroke cohort study: clinical trials protocol. Int J Stroke (2015) 10(4):636–44. doi:10.1111/ijs.12190
65. Ma H, Parsons MW, Christensen S, Campbell BCV, Churilov L, Connelly A, et al. A multicentre, randomized, double-blinded, placebo-controlled phase III study to investigate extending the time for thrombolysis in emergency neurological deficits (EXTEND). Int J Stroke (2012) 7(1):74–80. doi:10.1111/j.1747-4949.2011.00730.x
66. Brott T, Adams HP, Olinger CP, Marler JR, Barsan WG, Biller J, et al. Measurements of acute cerebral infarction: a clinical examination scale. Stroke (1989) 20(7):864–70. doi:10.1161/01.str.20.7.864
69. Janneke M, Gooskens F, Schepers VP, Schuurmans MJ, Lindeman E, Hafsteinsdóttir TB. Screening for poststroke depression using the patient health questionnaire. Nurs Res (2012) 61(5):333–41. doi:10.1097/NNR.0b013e31825d9e9e
72. Williams JBW, Kobak KA. Development and reliability of a structured interview guide for the Montgomery–Åsberg Depression Rating Scale (SIGMA). Br J Psychiatry (2008) 192:52–8. doi:10.1192/bjp.bp.106.032532
73. Dong Y, Sharma VK, Chan BP-L, Venketasubramanian N, Teoh HL, Seet RCS, et al. The Montreal Cognitive Assessment (MoCA) is superior to the mini-mental state examination (MMSE) for the detection of vascular cognitive impairment after acute stroke. J Neurol Sci (2010) 299(1):15–8. doi:10.1016/j.jns.2010.08.051
74. Cumming TB, Churilov L, Linden T, Bernhardt J. Montreal cognitive assessment and mini-mental state examination are both valid cognitive tools in stroke. Acta Neurol Scand (2013) 128(2):122–9. doi:10.1111/ane.12084
76. Ishihama Y, Rappsilber J, Mann M. Modular stop and go extraction tips with stacked disks for parallel and multidimensional peptide fractionation in proteomics. J Proteome Res (2006) 5(4):988–94. doi:10.1021/pr050385q
77. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotech (2008) 26:1367–72. doi:10.1038/nbt.1511
78. Cho H, Smalley DM, Theodorescu D, Ley K, Lee JK. Statistical identification of differentially labeled peptides from liquid chromatography tandem mass spectrometry. Proteomics (2007) 7(20):3681–92. doi:10.1002/pmic.200601034
84. Lee S, Shafe AC, Cowie MR. UK stroke incidence, mortality and cardiovascular risk management 1999-2008: time-trend analysis from the general practice research database. BMJ Open (2011) 1(2):e000269. doi:10.1136/bmjopen-2011-000269
85. Adams HP, Bendixen BH, Kappelle LJ, Biller J, Love BB, Gordon DL, et al. Classification of subtype of acute ischemic stroke. Definitions for use in a multicenter clinical trial. TOAST. Trial of Org 10172 in acute stroke treatment. Stroke (1993) 24(1):35–41. doi:10.1161/01.str.24.1.35
86. Sagen U, Vik TG, Moum T, Mørland T, Finset A, Dammen T. Screening for anxiety and depression after stroke: comparison of the hospital anxiety and depression scale and the Montgomery and Åsberg Depression Rating Scale. J Psychosom Res (2009) 67(4):325–32. doi:10.1016/j.jpsychores.2009.03.007
88. Chiu C-C, Huang C-C, Chan W-L, Chung C-M, Huang P-H, Lin S-J, et al. Increased risk of ischemic stroke in patients with systemic lupus erythematosus: a nationwide population-based study. Intern Med (2012) 51(1):17–21. doi:10.2169/internalmedicine.51.6154
92. Collard CD, Väkevä A, Morrissey MA, Agah A, Rollins SA, Reenstra WR, et al. Complement activation after oxidative stress: role of the lectin complement pathway. Am J Pathol (2000) 156(5):1549–56. doi:10.1016/S0002-9440(10)65026-2
93. Roos A, Bouwman LH, van Gijlswijk-Janssen DJ, Faber-Krol MC, Stahl GL, Daha MR. Human IgA activates the complement system via the mannan-binding lectin pathway. J Immunol (2001) 167(5):2861–8. doi:10.4049/jimmunol.167.5.2861
95. Stankovic S, Majkic-Singh N. Genetic aspects of ischemic stroke: coagulation, homocysteine, and lipoprotein metabolism as potential risk factors. Crit Rev Clin Lab Sci (2010) 47(2):72–123. doi:10.3109/10408361003791520
96. Fassbender K, Dempfle C-E, Mielke O, Schwartz A, Daffertshofer M, Eschenfelder C, et al. Changes in coagulation and fibrinolysis markers in acute ischemic stroke treated with recombinant tissue plasminogen activator. Stroke (1999) 30(10):2101–4. doi:10.1161/01.str.30.10.2101
97. Uchiyama S, Yamazaki M, Hara Y, Iwata M. Alterations of platelet, coagulation, and fibrinolysis markers in patients with acute ischemic stroke. Semin Thromb Hemost (1997) 23(6):535–41. doi:10.1055/s-2007-996132
98. Sofi F, Cesari F, Abbate R, Gensini GF, Broze G, Fedi S. A meta-analysis of potential risks of low levels of protein Z for diseases related to vascular thrombosis. Thromb Haemost (2010) 103(4):749–56. doi:10.1160/TH09-09-0645
99. Chou SH-Y, Lee P-S, Konigsberg RG, Gallacci D, Chiou T, Arai K, et al. Plasma-type gelsolin is decreased in human blood and cerebrospinal fluid after subarachnoid hemorrhage. Stroke (2011) 42(12):3624–7. doi:10.1161/strokeaha.111.631135
103. Agarwal V, Talens S, Grandits AM, Blom AM. A novel interaction between complement inhibitor C4b-binding protein and plasminogen that enhances plasminogen activation. J Biol Chem (2015) 290(30):18333–42. doi:10.1074/jbc.M114.619494
104. Amara U, Flierl MA, Rittirsch D, Klos A, Chen H, Acker B, et al. Molecular intercommunication between the complement and coagulation systems. J Immunol (2010) 185(9):5628–36. doi:10.4049/jimmunol.0903678
110. Sharp FR, Jickling GC, Stamova B, Tian Y, Zhan X, Liu D, et al. Molecular markers and mechanisms of stroke: RNA studies of blood in animals and humans. J Cereb Blood Flow Metab (2011) 31(7):1513–31. doi:10.1038/jcbfm.2011.45
111. Xu H, Tang Y, Liu D-Z, Ran R, Ander BP, Apperson M, et al. Gene expression in peripheral blood differs after cardioembolic compared with large-vessel atherosclerotic stroke: biomarkers for the etiology of ischemic stroke. J Cereb Blood Flow Metab (2008) 28(7):1320–8. doi:10.1038/jcbfm.2008.22
114. St-Amour I, Paré I, Alata W, Coulombe K, Ringuette-Goulet C, Drouin-Ouellet J, et al. Brain bioavailability of human intravenous immunoglobulin and its transport through the murine blood-brain barrier. J Cereb Blood Flow Metab (2013) 33(12):1983–92. doi:10.1038/jcbfm.2013.160
116. Macrez R, Ali C, Toutirais O, Le Mauff B, Defer G, Dirnagl U, et al. Stroke and the immune system: from pathophysiology to new therapeutic strategies. Lancet Neurol (2011) 10(5):471–80. doi:10.1016/S1474-4422(11)70066-7
117. Dirnagl U, Klehmet J, Braun JS, Harms H, Meisel C, Ziemssen T, et al. Stroke-induced immunodepression: experimental evidence and clinical relevance. Stroke (2007) 38(2):770–3. doi:10.1161/01.STR.0000251441.89665.bc
118. Offner H, Vandenbark AA, Hurn PD. Effect of experimental stroke on peripheral immunity: CNS ischemia induces profound immunosuppression. Neuroscience (2009) 158(3):1098–111. doi:10.1016/j.neuroscience.2008.05.033
120. Kubera M, Obuchowicz E, Goehler L, Brzeszcz J, Maes M. In animal models, psychosocial stress-induced (neuro)inflammation, apoptosis and reduced neurogenesis are associated to the onset of depression. Prog Neuropsychopharmacol Biol Psychiatry (2011) 35(3):744–59. doi:10.1016/j.pnpbp.2010.08.026
124. Shinjyo N, Ståhlberg A, Dragunow M, Pekny M, Pekna M. Complement-derived anaphylatoxin C3a regulates in vitro differentiation and migration of neural progenitor cells. Stem Cells (2009) 27(11):2824–32. doi:10.1002/stem.225
125. Kim GH, Mocco J, Hahn DK, Kellner CP, Komotar RJ, Ducruet AF, et al. Protective effect of C5a receptor inhibition after murine reperfused stroke. Neurosurgery (2008) 63(1):122–6. doi:10.1227/01.NEU.0000335079.70222.8D
126. Arumugam TV, Woodruff TM, Lathia JD, Selvaraj PK, Mattson MP, Taylor SM. Neuroprotection in stroke by complement inhibition and immunoglobulin therapy. Neuroscience (2009) 158(3):1074–89. doi:10.1016/j.neuroscience.2008.07.015
127. Van Beek J, Bernaudin M, Petit E, Gasque P, Nouvelot A, MacKenzie ET, et al. Expression of receptors for complement anaphylatoxins C3a and C5a following permanent focal cerebral ischemia in the mouse. Exp Neurol (2000) 161(1):373–82. doi:10.1006/exnr.1999.7273
128. Stephan AH, Barres BA, Stevens B. The complement system: an unexpected role in synaptic pruning during development and disease. Annu Rev Neurosci (2012) 35(1):369–89. doi:10.1146/annurev-neuro-061010-113810
129. Atkinson C, Zhu H, Qiao F, Varela JC, Yu J, Song H, et al. Complement-dependent P-selectin expression and injury following ischemic stroke. J Immunol (2006) 177(10):7266–74. doi:10.4049/jimmunol.177.10.7266
130. Rahpeymai Y, Hietala MA, Wilhelmsson U, Fotheringham A, Davies I, Nilsson A-K, et al. Complement: a novel factor in basal and ischemia-induced neurogenesis. EMBO J (2006) 25(6):1364–74. doi:10.1038/sj.emboj.7601004
131. Alawieh A, Elvington A, Zhu H, Yu J, Kindy MS, Atkinson C, et al. Modulation of post-stroke degenerative and regenerative processes and subacute protection by site-targeted inhibition of the alternative pathway of complement. J Neuroinflammation (2015) 12(1):247. doi:10.1186/s12974-015-0464-8
132. Elvington A, Atkinson C, Zhu H, Yu J, Takahashi K, Stahl GL, et al. The alternative complement pathway propagates inflammation and injury in murine ischemic stroke. J Immunol (2012) 189(9):4640–7. doi:10.4049/jimmunol.1201904
133. Mack WJ, Ducruet AF, Hickman ZL, Garrett MC, Albert EJ, Kellner CP, et al. Early plasma complement C3a levels correlate with functional outcome after aneurysmal subarachnoid hemorrhage. Neurosurgery (2007) 61(2):255–61. doi:10.1227/01.NEU.0000255518.96837.8E
134. Cavusoglu E, Eng C, Chopra V, Ruwende C, Yanamadala S, Clark LT, et al. Usefulness of the serum complement component C4 as a predictor of stroke in patients with known or suspected coronary artery disease referred for coronary angiography. Am J Cardiol (2007) 100(2):164–8. doi:10.1016/j.amjcard.2007.02.075
137. Muscari A, Bozzoli C, Puddu GM, Sangiorgi Z, Dormi A, Rovinetti C, et al. Association of serum C3 levels with the risk of myocardial infarction. Am J Med (1995) 98(4):357–64. doi:10.1016/S0002-9343(99)80314-3
138. Gelderblom M, Leypoldt F, Steinbach K, Behrens D, Choe C-U, Siler DA, et al. Temporal and spatial dynamics of cerebral immune cell accumulation in stroke. Stroke (2009) 40(5):1849–57. doi:10.1161/strokeaha.108.534503
Keywords: ischemic stroke, proteomics, post-stroke depression, complement and coagulation, immunity and inflammation, stroke neurological recovery, blood biomarkers
Citation: Nguyen VA, Carey LM, Giummarra L, Faou P, Cooke I, Howells DW, Tse T, Macaulay SL, Ma H, Davis SM, Donnan GA and Crewther SG (2016) A Pathway Proteomic Profile of Ischemic Stroke Survivors Reveals Innate Immune Dysfunction in Association with Mild Symptoms of Depression – A Pilot Study. Front. Neurol. 7:85. doi: 10.3389/fneur.2016.00085
Received: 23 February 2016; Accepted: 23 May 2016;
Published: 14 June 2016
Edited by:Jean-Claude Baron, University of Cambridge, UK
Reviewed by:Craig James Smith, Salford Royal NHS Foundation Trust, UK
Carine Ali, University of Caen Normandy, France
Teresa García-Berrocoso, Vall d’Hebron Research Institute (VHIR), Spain
Copyright: © 2016 Nguyen, Carey, Giummarra, Faou, Cooke, Howells, Tse, Macaulay, Ma, Davis, Donnan and Crewther. 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) or licensor 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: Leeanne M. Carey, firstname.lastname@example.org
†Joint senior authors.