- 1Department of Occupational Therapy, La Trobe University, Bundoora, VIC, Australia
- 2Department of Psychology and Counselling, La Trobe University, Bundoora, VIC, Australia
- 3Neurorehabilitation and Recovery, Stroke, The Florey Institute of Neuroscience and Mental Health, Heidelberg, VIC, Australia
- 4Western Health, Department of Neurology, Sunshine, VIC, Australia
- 5Department of Biochemistry and Genetics, La Trobe University, Bundoora, VIC, Australia
- 6Medical Sciences Precinct, University of Tasmania, Hobart, TAS, Australia
- 7Faculty of Health and Medical Sciences, Internal Medicine, University of Western Australia, Perth, WA, Australia
- 8Clinical Research, Harry Perkins Institute of Medical Research, Perth, WA, Australia
- 9Department of Medicine, The University of Melbourne, Sunshine, VIC, Australia
- 10Monash Health, Neurology and Stroke, Clayton, VIC, Australia
- 11Department of Neurology, Royal Melbourne Hospital, Parkville, VIC, Australia
Currently the longitudinal proteomic profile of post-ischemic stroke recovery is relatively unknown with few well-accepted biomarkers or understanding of the biological systems that underpin recovery. We aimed to characterize plasma derived biological pathways associated with recovery during the first year post event using a discovery proteomics workflow coupled with a topological pathway systems biology approach. Blood samples (n = 180, ethylenediaminetetraacetic acid plasma) were collected from a subgroup of 60 first episode stroke survivors from the Australian START study at 3 timepoints: 3–7 days (T1), 3-months (T2) and 12-months (T3) post-stroke. Samples were analyzed by liquid chromatography mass spectrometry using label-free quantification (data available at ProteomeXchange with identifier PXD015006). Differential expression analysis revealed that 29 proteins between T1 and T2, and 33 proteins between T1 and T3 were significantly different, with 18 proteins commonly differentially expressed across the two time periods. Pathway analysis was conducted using Gene Graph Enrichment Analysis on both the Kyoto Encyclopedia of Genes and Genomes and Reactome databases. Pathway analysis revealed that the significantly differentiated proteins between T1 and T2 were consistently found to belong to the complement pathway. Further correlational analyses utilized to examine the changes in regulatory effects of proteins over time identified significant inhibitory regulation of clusterin on complement component 9. Longitudinal post-stroke blood proteomics profiles suggest that the alternative pathway of complement activation remains in a state of higher activation from 3-7 days to 3 months post-stroke, while simultaneously being regulated by clusterin and vitronectin. These findings also suggest that post-stroke induced sterile inflammation and immunosuppression could inhibit recovery within the 3-month window post-stroke.
Introduction
Ischemic stroke covers a variety of cerebrovascular events that affect up to 800,000 people in the United States every year, with 133,000 deaths reported in 2017 (16.74%) (1). Of the survivors, 30% are reported to experience prolonged cognitive impairment (2) and depressive symptoms at any point 5 years post-stroke (3). Currently there are few well-accepted biomarkers for recovery and comparatively little literature exploring the biological systems that drive recovery or even the most optimal times for monitoring biological and behavioral recovery. Evidence from stroke rehabilitation studies suggest the greatest efficacy for motor-based rehabilitation is within this 3 month time window (4), though recovery may continue at a slower rate over subsequent months and years. Although there has been increasing research examining the blood biomarkers of stroke recovery (5), the additional linking of biomarkers to biological systems remains speculative. Hence, we aimed to investigate the changes in the molecular profile of proteins in plasma samples via a mass spectrometry (MS) based discovery proteomics approach (6). Mass spectrometry and nuclear magnetic resonance (NMR) based techniques examining protein expression are among the most versatile techniques for protein identification and quantification, with the ability to address a wide range of biological samples, especially plasma, and serum (7, 8). Proteomics utilizes the advantage of systems biology techniques to quantify a large number of analytes in an exploratory fashion, with a computational bioinformatics approach to further categorize biomarkers into biosystems (9).
Proteomics have recently been used to pursue multiple clinical questions within stroke research, relating to differentiation of ischemic from hemorrhagic stroke (10–12) and investigations of potential biomarkers involved in post-stroke recovery (13, 14). Although traditional bioinformatics methods were originally developed to accommodate gene expression data, proteomics studies can utilize these methods to organize and visualize findings by adopting standardized change scores and using annotations that are common across proteomics and genomics (15). Indeed, our laboratory has previously used proteomic methods and Gene Set Enrichment Analysis (GSEA) to investigate the relationship between protein changes in plasma at 3-months post-stroke and affective (depression) outcomes (16). The results indicated that proteins involved in the complement but not the coagulation pathway of the immune system are likely to be associated with post-stroke depressive symptoms (14). The complement system is recognized as an innate immune pathway that contributes to primary host defence by encouraging phagocytosis of unwanted cells. This new study aims to expand upon our earlier single time point study by using discovery proteomics to identify longitudinal changes in blood plasma protein expression over the post-stroke timeline of recovery; specifically 3 timepoints post ischemic stroke: 3–7 days (T1), 3-months (T2) and 12-months (T3). We also aimed to improve upon our previous set-based functional annotation methods by utilizing Gene Graph Enrichment Analysis (GGEA). The GGEA approach differs from the GSEA approach by further addition to traditional set-based functional annotations by incorporation of networks of established biochemical relationships (gene to gene) using topological omics databases such as Reactome (17) and KEGG (18) for the analysis of structural representation of biological pathways in the analytical workflow. This more novel approach further addresses the regulatory mechanisms in gene and protein pathways by examining co-expression and co-regulation networks using correlation analyses (19). Examining the changes in correlational strength also allows for quantification of the changes in the regulatory effect of proteins between timepoints. Although discovery approaches are hypothesis-free by nature, based on our previous study (14) and others (20) that suggest inflammatory and immune homeostasis will be disrupted in the post-stroke recovery timeline, we hypothesize that the complement system will be dysregulated when comparing early 3–7 days post-stroke to later 3-month and 12–month post-stroke timepoints.
Materials and Methods
Subjects
Data from a subset of ischemic stroke patients was obtained from the longitudinal stroke cohort known as START, which comprised participants from START_PrePARE (STroke imAging pRevention and Treatment: Prediction and Prevention to Achieve optimal Recovery Endpoints; (21) (Neuroscience Trials Australia: NTA 0902) and START_EXTEND [STroke imAging pRevention and Treatment: EXtending the time for Thrombolysis in Emergency Neurological Deficits (22) (Neuroscience Trials Australia: NTA 0901, Clinicaltrials.gov number: NCT01580839)] studies. These prospective, integrated studies were longitudinal and provided long term plasma and serum samples at 3 time-points post-stroke. Selection of the subset ischemic stroke patients for the current study (n = 106) was initially based on the availability of ethylenediaminetetraacetic acid (EDTA) treated plasma samples at each of the 3 timepoints: 3–7 days (T1), 3-months (T2) and 12 months (T3). Laboratory batch processing limitations and time/cost considerations allowed for 180 EDTA plasma samples to be processed, i.e., samples from 60 stroke patients across 3 timepoints. We therefore conducted a stratified selection process involving patients with complete blood samples at the 3 repeated times and full clinical scores, to ensure a spread of scores across clinical outcomes of stroke severity, mood, and cognition. As expected, the current subset contained only participants who completed the 3 sessions in the 12 months post first acute stroke event. See Table 1 below for baseline characteristics for patients (43 male) and Table 2 for clinical characteristics at T1, T2, and T3. The total START_PrePARE and START_EXTEND cohort comprised 219 patients, with 21 deaths, 8 lost to follow-up and 2 withdrawals. See Supplementary Table 1 for a comparison between the selected subset and the total cohort. Healthy control data was not available as part of the original protocols for the START_PrePARE and START EXTEND studies. Ethics was approved by the Human Research Ethics Committee of Austin Hospital, Heidelberg (HREC code: H2010/03588), and relevant university and hospital sites.
Blood Collection and Storage
All samples were collected in Benton-Dickson (BD) EDTA-coated 4 ml vacutainers and were mixed and left to stand in ambient room temperature for 30 min. Average time to blood draw was 3.45 days (range = 2.03–7.58 days, SD = 1.174) after stroke onset for T1. The tubes were then centrifuged at 1100–1300 g at room temperature and the resulting plasma was aliquoted into cryotubes and immediately stored in a −80oC freezer. For transport from the central study freezer to the analysis site, temperature was kept at −70 to −80oC on dry ice before transfer into a −80oC freezer. This procedure was consistent for the 3-month and 12-month follow up periods.
Mass Spectrometry
Label-free quantitation (LFQ) proteomics was conducted on the Q Exactive HF Orbitrap instrument (Thermo-Fisher Scientific). Details of the sample preparation, instrument parameters and protein identification are available in the Supplementary Methods.
Bioinformatics and Statistical Analysis
The initial MaxQuant output consisted of 358 identified proteins across 180 samples. After removal of site-identified contaminants and filtering for at least 80% of data present across the protein, 163 proteins remained (Supplementary Data 1). This procedure regarding missing values is an expected step in addressing proteomics data as the protein identification search on MaxQuant (utilizing an FDR of 1%) will not output as an identified protein in the final data file unless there is a single unique peptide, regardless of peptide length. The Limma R package (version 3.34.9) with Benjamini–Hochberg (BH) adjustments for familywise error rate was used to compare which proteins were differentially expressed across T1, T2, and T3. To further understand the functional organization of significant protein sets, the EnrichmentBrowser R package (version 2.09.17) was employed with GGEA as the network-based enrichment method across the KEGG and Reactome databases using default settings (minimum set size = 3, maximum set size = 500, permutations = 1000). Correlations with BH corrections were used to further explore the changes in proteins that regulate identified pathways, with Fisher r to z transformations conducted to explore changes in regulatory effect across timepoints. Statistical significance levels were set to a =.05 after multiple comparisons adjustment.
Data Availability
Due to the requirements of ethics and the nature of ongoing clinical trials with the START cohorts, unidentified patient clinical data may only be made available upon request. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (23) partner repository with the dataset identifier PXD015006.
Results
The differential expression (DE) analysis revealed that 29 proteins significantly differed between T1 and T2, and that 33 proteins significantly differed between T1 and T3 (Table 3) with 8 proteins at FDR of 0.05 indicating that the effects are higher than would occur due to chance alone. The changes between T1 to T2 and T1 to T3 constitute 17.79 and 20.25%, respectively, in proportion to the total number of proteins identified. Eleven proteins were uniquely expressed between T1 and T2 and 15 proteins were uniquely expressed between T1 and T3, with 18 of the same proteins expressed both between T1 and T2 and T1 and T3. There were no significant differences in protein expression between T2 and T3, potentially suggesting that the currently identified proteome in post-stroke survivors does not change significantly between 3- and 12-month times. See Supplementary Data 1 for a full list of proteins and fold change values across all comparisons.
The list of DE proteins and the full expression matrix were submitted to EnrichmentBrowser using GGEA as the network enrichment method (1000 permutations and α = 0.05). This revealed significant functional annotation of our set of proteins in both KEGG and Reactome databases only between T1 and T2, after BH adjustment of p-values (Table 4) (see Supplementary Data 2 a full list of nominally significant sets from KEGG and Reactome). In the case of the Reactome database, the nested structures in the ontological organization of the pathways can lead to redundancy resulting in identification of multiple pathways, especially as the “Complement Cascade” is located as a subset of the “Innate Immune System” and the “Immune System.” Further interpretation is needed to account for the number and type of identified proteins between the reference pathways and experimental data and should occur at the level closest to the relevant biological processes. Both KEGG and Reactome databases concurred on the identification of the complement system as the significantly enriched network in these plasma samples.
Table 4. Significantly Enriched Gene Pathways between T1 (3–7 days) and T2 (3-month) timepoints in Stroke Survivors.
A visual illustration of this enriched pathway on the Reactome database reveals the changes and interactions of the detected molecules on the complement pathway (Supplementary Figure 1), displaying information on nodes (proteins) and edges (the correlation between nodes). Similarly, the results from the KEGG database includes elements of the complement pathway but also the coagulation cascades (Supplementary Figure 2). As the complexity of information is difficult to interpret without a significant degree of system specific understanding in these pathway diagrams, a composite network pathway was created to amalgamate the statistical relationships found in both diagrams pertaining to the complement system as shown in Figure 1.
Figure 1. A pathway diagram contrasting information from differential proteomic analysis of KEGG and Reactome GGEA complement system pathways at 3–7 days (T1) and 3-months (T2) post-stroke. This pathway is activated by the C1q complex antigen linked immunoglobulin through the classical pathway and pathogen surface linked mannose binding linked (MBL) proteins through the MASP complex. The alternative pathway of complement activation follows a constant and low state of activation by a feedback loop that is heavily regulated by complement factor I (CFI) and complement factor H (CFH). Simultaneous upregulation of complement factor D (CFD) and CFI suggest that the alternative pathway is undergoing activation but also being regulated to prevent autoimmune insult. Recently, the link between coagulation and complement has been established through thrombin, plasmin and factor XIIa, these factors were not identified or significantly related to complement in this model. All activation pathways lead to the promotion of cleavage of fluid phase complement component 3 (C3) in the bloodstream, with heavier β chains forming convertases downstream and lighter α chains such as C3a and complement component 5α (C5a) able to signal for potent local and systemic inflammation responses, with the most extreme being anaphylactic shock. The final stage of complement activation is the membrane attack complex (MAC), a large protein complex that is constructed on cells to disrupt the outer membrane and promotes active cytolysis. The profile of significant differences in MAC proteins such as C8G and C9 suggest plasma complement regulation by vitronectin (VTN) or clusterin (CLU). Image adapted from KEGG [map04610: Complement and coagulation cascades (homo sapiens)] with permission.
Some of the regulatory proteins in the complement pathway (CLU & VTN) were not included in the pathway diagrams for KEGG and Reactome databases used in this analysis. Therefore, a correlation matrix with BH corrected p-values was created to explore the direction and strength of the regulatory effects between regulator and related proteins (Table 5). Fisher's transformation of Pearson's r to normally distributed z scores was computed to compare changes in regulatory effects between T2 and T1. While not overall presenting significant results, examination of the rdiff(T2-T1) complement regulators on their target proteins show that the middle stage complement regulators shifted their effects forward, resulting in a reduced inhibitory effect. Only the inhibitory effect of CLU on C9 was significantly increased from T1 to T2, with VTN also showing a large effect size for changes in C9 regulation.
Table 5. Pearson's Correlation matrix with Fisher r to z transformations for complement proteins and their regulators between T1 and T2.
Discussion
Currently the interaction of the many biological mechanisms limiting post-stroke recovery remain poorly understood. To better understand the biochemical pathways impacting/affecting the post-stroke timeline of recovery, this study employed a discovery approach utilizing mass spectrometry to examine protein expression in patients at 3–7 days (T1), 3-months (T2), and 12-months (T3) post-stroke. Two sets of proteins were identified to be significantly different based on differential expression analysis from T1 to T2 and T1 to T3, but not between T2 and T3. Of these proteins, complement (C8A, C9, CFI), apolipoproteins (APOA4, APOA2, APOD) and membrane bound proteins (TF, ACTB) were highly overexpressed; consistent with high abundance typically found in human blood samples (24). The lists of significantly differently expressed proteins were analyzed using the GGEA bioinformatics algorithm based on the topological consistency of the observed data sets compared to database defined reference pathways. This revealed that the proteins identified in this longitudinal experiment significantly conformed to pathways central to the complement system in both KEGG and Reactome databases. To our knowledge, the current study is first published use of the GGEA approach in clinical proteomics of human blood.
There is currently limited knowledge of the changes in complement system in post-stroke recovery, especially relative to blood plasma concentrations in similar aged individuals (mean age 68 ± 14 years). Furthermore, most studies examining the complement system post-stroke have described its damage-exacerbating role in acute stroke human and animal models but have not examined changes in levels longitudinally over time post-stroke (25–29). Many of the studies exploring complement in stroke have examined the system in the context of danger associated molecular pattern (DAMP) signaling and subsequent neuronal repair in central nervous system (CNS) specific injury (30), with plasma and serum based studies examining levels in relation to clinical outcomes (31, 32). The endothelial junctions forming the blood brain barrier (BBB) have been traditionally viewed as preventative to the entry of large peripherally derived molecules such as complement into the cerebral parenchyma (33, 34). The CNS has also been recognized as able to endogenously biosynthesize complement proteins in glial cells (35). Although disruption to BBB permeability allows for the passage of complement proteins from periphery to CNS (34), to our knowledge, no study has simultaneously examined the levels and functional activity of the complement system in both domains. Therefore, the results of this study have been interpreted as referring to the circulating fluid phase proteins that enhance the decay of complement proteins by cleaving active proteins through allosteric binding (25).
The complement pathway has three traditional modes of activation, the classical, lectin and alternative pathways. The classical and lectin pathways are activated by antigen and pathogen binding complexes, respectively, whereas the alternative pathway is an internal biomechanical activation loop that regulates downstream pathway activity (36). In the classical pathway of complement activation, upregulation of SERPING1 or the C1 inhibitor protein is responsible for regulation of the C1 complex. Specifically, SERPING1 strongly binds to the C1r and C1s proteases and also to the activation units of the lectin pathway [mannose binding lectin, MBL, and associated proteases MASP1 and MASP2 (37)], effectively inhibiting the initial activation effects of the complement system in both antigen and pathogen related pathways. In addition to the effects on the complement system, high concentrations of this protease can further inhibit leukocyte-endothelial adhesion and vascular rolling via interference of cell adhesion molecules (CAMs) or selectins, the first physiological requirements for trans-endothelial leukocyte infiltration (38). Although a fold change of 0.008 or 0.8% for SERPING1 was detected between T1 and T2 in this study, previous reports have indicated that a fold change can be observed during acute inflammation and that higher levels may be therapeutic in preventing autoimmune injury by providing a stop mechanism to complement system activation (39, 40). These acute properties have been demonstrated in animal models of post-stroke middle cerebral artery occlusion (MCAO) and shown to attenuate ischemic reperfusion injury (41, 42). Limited evidence for longitudinal changes in levels of SERPING1 in human cardiovascular disease have previously been suggested as a link to low-grade levels of chronic inflammation (43). Furthermore, the incidence rates of post-stroke immunological trajectory of post-stroke recovery may also be explained by the immune inhibitory effects of SERPING1, following the initial phase of post-stroke immune flux and subsequent immunosuppression (44).
The cleavage of C3 is the central amplification step of the whole complement pathway, with active C3b forming the C5 convertase and C3a inducing further phagocytotic chemotaxis (29). Given that C3 is central in all complement activation pathways, especially with the alternative pathway requiring C3 hydrolysis, much of the theorized regulatory activity in this pathway is focused on inhibition of C3 convertase formation and decay rate (36). Of the complement regulatory proteins, CFH and CFI were overexpressed, with only plasma CFI levels shown to be statistically different. Additionally, the changes in regulatory effect for middle stage complement proteins show that even though CFH and CFI are upregulated between T1 and T2, their inhibitory effects were reduced. Although CFI and CFH are cofactors of C3 convertase (not identified in this study) inhibition and not C3 itself, the overall profile suggests a relative lack of effectiveness in the functional ability of these proteins to regulate alternative complement pathway activity in this older group. Furthermore, this profile may also be indicative of overexpression in levels of CD55 or decay accelerating factor (DAF), a protein that inhibits C3 and C5 convertase formation, albeit on cell surfaces (45). Complement C3a has previously been shown to be acutely elevated from 1 to 28 days post-stroke (46, 47). Indeed, ongoing anaphylatoxin signaling from components such as C3a and C5a are damaging to host systems, especially in cerebral ischemia (25). The quantification of properdin, a protein that acts as a positive regulator of C3 and C5 convertases (48), would further ameliorate understanding of alternative pathway activation in this sample.
Complement component C5 is the important functional unit of the pathway. Cleavage of soluble C5 protein into C5b is able to create the porous membrane attack complex (MAC) on targeted cells to induce cell death by cytolysis while C5a is able to produce local and system wide inflammatory cascades that are 20 times more potent than C3a (49). The profile here shows that C5 is upregulated at 3-months, despite reduced C3 levels. This suggests that the reduction in C3 levels may be indicative of active cleaving to produce downstream increased C3b and C3a levels to produce or maintain constant immune system homeostasis. The detection of increased C5 may also imply involvement of the novel complement pathway or the extrinsic complement pathway, linking coagulation, and complement systems. This pathway was originally thought to be activated based on thrombin functionally substituting for the C3 dependent C5 convertase (50). Recently, a study has also demonstrated that in thrombosis, interactions with plasmin in both surface bound and fluid phase complement is capable of upregulating C5a and C5b whether in venous or arterial thrombi (51). In our analysis however, plasminogen was not found to be significantly differentially expressed or related to the complement pathway (Supplementary Figure 2).
Our methodology identified and quantified all the molecular components in MAC formation pathway. Our findings demonstrate that C6 and C8A were significantly upregulated and C8G and C9 were significantly downregulated between T1 and T2. Regulators of MAC formation exhibit function by preventing polymerisation of the C8 or C9 complexes, thereby inhibiting cytolytic function on cell surfaces (52). Pore formation is heavily regulated by CD59 or MAC-inhibitory protein on host cell surfaces to prevent further immune insult resulting from early tissue induced hyperinflammation and gradual immune dysregulation (53). As this expression profile was obtained in plasma, the expression profile here is suggestive of fluid phase regulation of MAC formation, dependent on vitronectin (VTN or S-protein) and clusterin (CLU or apolipoprotein J). To create a successful MAC, from 12 to 16 molecules of C9 are needed as the final step to form a cyclical pore in addition to 1 of each previous MAC protein in the chain (54). Both VTN and CLU inhibit the formation of the MAC insertion C5b-7 complex and also inhibit the ability of C9 to bind to inserted MAC complexes to create functioning lysis pores (54). VTN and CLU were both found to be overexpressed, with the increases in VTN shown to be statistically significant (36). Overall, the profile of the complement system here appears to be dysregulated, with an increase in alternative pathway activation by CFB, and reduced inhibitory effectiveness of all middle stage complement proteins but especially CFI and CFH. The apparent contradiction in expected abundances of MAC proteins suggests that the pathway is also being inhibited by CLU and VTN.
Although few studies have presented this specific profile of complement expression, and currently none in a complex proteomics pathway model, there are some interpretations relevant to stroke survival and recovery that can be made here. From the upregulation of C5 and initial MAC proteins, it is possible that C5 is being continually cleaved to produce the pro-inflammatory chemokine C5a. This can be understood in theory of post-stroke immune profiles where both central and peripheral systems never fully recover to pre-stroke levels and establish new albeit dysfunctionally higher inflammation profiles (44).
The profile of complement activation described in this study is likely to indicate that
1) the internal activation loop that readies the immune system for pathological threats and tissue damage (the alternative pathway) is activated,
2) middle stage complement regulators that are normally expected to be activated to prevent indiscriminatory autoimmune damage are overexpressed but their intended inhibitory effects (CFI –> C3) are reduced and
3) although early chain MAC proteins (C6, C7, C8A, C8B) trend toward overexpression by time T2 are likely due to alternative pathway activation, CLU and VTN are also overexpressed, possibly as a compensatory mechanism for alternative pathway activation and thereby preventing autoimmune MAC formation by inhibiting C9 expression on host/friendly cells.
Our results did not find/demonstrate an association between C8G expression and CLU or VTN, although downregulation of C8G is theorized to be coincidental with immunosuppression (55). It is speculated that based on the positional role of C8G as the protein that provides a binding site for the first C9 molecule in the structural formation of the MAC pore [see (56) for an overview of MAC protein structure], the significant downregulation of C8G between T1 and T2 can likely be attributed to inhibitory CD55 or CD59 (that was not present in our proteomics analysis). Furthermore, there are additional molecules and mechanisms that may regulate MAC formation through C8G expression such as anti-apoptotic processes that promote cell survival (57). This may also explain the immunosuppressed state of post-stroke patients and infections (58), whereby the system is not able to react as efficiently to acquired pathogens such as viral pneumonia due to increased fluid phase MAC regulation (59).
Limitations and Future Research
These results are limited to interpretation based on the sample characteristics of a cohort of mild stroke patients over a 12-month period post-stroke. During this period, it is recognized that other potential concurrent processes such as recurrent stroke, infections, venous thromboembolism could impact on the results: however, we were not able to control or adjust for such factors in the current analyses. The lack of an age matched control sample represents a challenge to external validity, and it is uncertain if these changes could also occur in people without stroke, as a process of senescence and aging or other comorbidities. Additionally, our results were not controlled for by age and baseline stroke as the analyses using the EnrichmentBrowser package did not have a streamlined function to address confounds. The profile of complement expression in this sample may also have been influenced by the myriad of post-stroke medications that currently include aspirin, warfarin, and clopidergel. Although these pharmacological agents target aspects of the coagulation system, it is possible that they also exhibit regulatory effects on the complement system (60) via the intersection between hemostasis and innate immunity (61). These confounds remain difficult to statistically control in the context of bioinformatics models, and should form part of the conceptualization in the design phase in future studies (62). While the focus of the current study was to explore evidence of change in biological pathways over the first-year post-stroke, it is recommended that future studies investigate the relationship between changes in these pathways and changes with age as well as in neurological and functional recovery profiles.
The sample preparation methods used in this study did not include fractionation or depletion of the plasma samples, leading to a low number of proteins identified relative to the blood proteome. Blood samples represent a very complex sample with high dynamic range for proteomics analysis, with a proportionally high abundance of proteins such as albumin, haptoglobin, fibrinogen, and immunoglobulin G (63). Future studies could increase the range of identified proteins by depleting abundant plasma proteins with various methods before MS analysis (64, 65); although some popular depletion columns also remove proteins identified here as differentially expressed such as those within the complement system (66). In addition to this, limitations to the discovery approach imply that there is a degree of uncertainty in the ability to identify a complete group of proteins (67, 68), especially where set and network-based bioinformatics are the intended analyses. Therefore, after pathways are identified by robust techniques such as GGEA, the integrity of these pathways can be enhanced by examining proteins that are missing in the pathway and directly related to pathway activation and function by targeted proteomics techniques such as single or mass reaction monitoring (SRM or MRM) in MS with immunoaffinity enrichment (69).
In this study, the LC-MS methodology did not detect the specific formation of C3 or C5 chains or convertase protein complexes. Therefore, the pathway representations of relative abundance for C3 and C5 may not necessarily reflect the abundance of downstream protein chains of the active products. Specifically, C3a and C5a would be interesting proteins to target as they have potent effects on inflammation physiology and would be of concern especially if shown to be elevated at 3 months post-stroke. Fluid phase identification of inactive pre-MAC complexes with protein S such as sC5b-7, 8 and 9 in addition to VTN and CLU may also provide further information on the regulatory activity of the MAC in stroke survivors. Furthermore, the longitudinal profile of complement in this study and others (25, 47) suggests that complement dysregulation begins in the first week post-stroke and continues for at least 3 months, and may set a post-stroke level of immune homeostasis. This finding needs to be further validated at other timepoints within this window, especially in relation to infection based mortality (70). To this extent, it is possible that the post-stroke immune profile is settled at 3 weeks to 1 month post-stroke when adaptive immunity is engaged in phagocytizing dead cells (44) and resolution of cerebral microedema (71).
Conclusion
The biology of post-stroke recovery is not well-understood, with patients exhibiting varying profiles based on factors such as their age, location of lesion, and degree of stroke damage (72). This study aimed to characterize common peripheral biological systems involved in post-stroke recovery by examining the longitudinal proteomics profile of EDTA plasma in stroke survivors. Specifically, elements of the alternative pathway of complement system and MAC proteins, i.e., CFI, C6, C8A, C8G, C9 were found to be activated but also undergoing regulation at 3 months post-stroke. This increased turnover may lead to the upregulation of anaphylatoxins C3a and C5a that could explain the prolonged sterile inflammation profile of post-stroke survivors. These results suggest a biomechanism for post-stroke immunosuppression by complement system regulatory proteins. This knowledge may be useful to guide assessment in the clinical setting for post-stroke infections and immune recovery in the 3-month recovery window. For example, it may be used to inform future recommendations for in-hospital and post-stroke infections by suggesting clinical utility of complement panel blood tests. Future investigations should examine the clinical implications of this biological profile and determine the temporal trajectory of immune homeostasis post-stroke.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: http://www.proteomexchange.org/, PXD015006.
Ethics Statement
The studies involving human participants were reviewed and approved by Human Research Ethics Committee of Austin Health, all participating Hospitals and La Trobe University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
VN, NR, SC, LC, and GH contributed to project conception, data analysis, and draft editing. PF and HR contributed to data analysis and project conception. TW, HM, DH, SD, and GD contributed to project conception. All authors contributed to the article and approved the submitted version.
Funding
We acknowledge financial support for conduct of the START_EXTEND and START_PrePARE studies from the Commonwealth Scientific and Industrial Research Organization (CSIRO) of Australia, Flagship Collaboration Fund through the Preventative Health Flagship and support for analysis, write up and researchers from the James S. McDonnell Foundation 21st Century Science Initiative in Cognitive Rehabilitation -Collaborative Award (# 220020413); NHMRC Centre of Research Excellence in Stroke Rehabilitation and Brain Injury (#1077898); NHMRC program grant Saving brain and changing practice in stroke (#1113352); 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 awarded to VN. The funding sources had no role in the conduct of this study or the writing of this report.
Conflict of Interest
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.
Acknowledgments
The authors acknowledge the stroke patients and all research personnel associated with the START_EXTEND and START_PrePARE teams.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.00692/full#supplementary-material
References
1. Benjamin EJ, Blaha MJ, Chiuve SE, Cushman M, Das SR, Deo R, et al. Heart disease and stroke statistics−2017 update: a report from the American heart association. Circulation. (2017) 135:e146. doi: 10.1161/CIR.0000000000000491
2. Sun JH, Tan L, Yu JT. Post-stroke cognitive impairment: epidemiology, mechanisms and management. Ann Transl Med. (2014) 2:80. doi: 10.3978/j.issn.2305-5839.2014.08.05
3. Robinson RG, Jorge RE. Post-stroke depression: a review. Am J Psychiatry. (2016) 173:221–31. doi: 10.1176/appi.ajp.2015.15030363
4. Hatem SM, Saussez G, della Faille M, Prist V, Zhang X, Dispa D, et al. Rehabilitation of motor function after stroke: a multiple systematic review focused on techniques to stimulate upper extremity recovery. Front Hum Neurosci. (2016) 10:442. doi: 10.3389/fnhum.2016.00442
5. Bang OY. Advances in biomarker for stroke patients: from marker to regulator. Precis Future Med. (2017) 1:32–42. doi: 10.23838/pfm.2017.00052
6. Aslam B, Basit M, Nisar MA, Rasool MH, Khurshid M. Proteomics: technologies and their applications. J Chromatogr Sci. (2017) 55:182–96. doi: 10.1093/chromsci/bmw167
7. Xue J, Huang W, Chen X, Li Q, Cai Z, Yu T, et al. Neutrophil-to-lymphocyte ratio is a prognostic marker in acute ischemic stroke. J Stroke Cerebrovasc Dis. (2017) 26:650–7. doi: 10.1016/j.jstrokecerebrovasdis.2016.11.010
8. Gowda G, Djukovic D. Overview of mass spectrometry-based metabolomics: opportunities and challenges. Methods Mol Biol. (2014) 1198:3–12. doi: 10.1007/978-1-4939-1258-2_1
9. Chandramouli K, Qian P. Proteomics: challenges, techniques and possibilities to overcome biological sample complexity. Hum Genomics Proteomics. (2009) 2009:239204. doi: 10.4061/2009/239204
10. Penn AM, Bibok MB, Saly VK, Coutts SB, Lesperance ML, Balshaw RF, et al. Validation of a proteomic biomarker panel to diagnose minor-stroke and transient ischaemic attack: phase 2 of SpecTRA, a large scale translational study. Biomarkers. (2018) 23:793–803. doi: 10.1080/1354750X.2018.1499130
11. Penn AM, Saly V, Trivedi A, Lesperance ML, Votova K, Jackson AM, et al. Differential proteomics for distinguishing ischemic stroke from controls: a pilot study of the SpecTRA project. Transl Stroke Res. (2018) 9:590–9. doi: 10.1007/s12975-018-0609-z
12. Montaner J, Mendioroz M, Delgado P, García-Berrocoso T, Giralt D, Merino C, et al. Differentiating ischemic from hemorrhagic stroke using plasma biomarkers: the S100B/RAGE pathway. J Proteomics. (2012) 75:4758–65. doi: 10.1016/j.jprot.2012.01.033
13. Zhou F, Zhou L, Guo T, Wang N, Hao H, Zhou Y, et al. Plasma proteomics reveals coagulation, inflammation, and metabolic shifts in H-type hypertension patients with and without acute ischemic stroke. Oncotarget. (2017) 8:100384–95. doi: 10.18632/oncotarget.22233
14. Nguyen VA, Carey LM, Giummarra L, Faou P, Cooke I, Howells DW, et al. A pathway proteomic profile of ischemic stroke survivors reveals innate immune dysfunction in association with mild symptoms of depression – a pilot study. Front Neurol. (2016) 7:85. doi: 10.3389/fneur.2016.00085
15. Schmidt A, Forne I, Imhof A. Bioinformatic analysis of proteomics data. BMC Syst Biol. (2014) 8:S3. doi: 10.1186/1752-0509-8-S2-S3
16. 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 USA. (2005) 102:15545. doi: 10.1073/pnas.0506580102
17. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The reactome pathway knowledgebase. Nucleic Acids Res. (2018) 46:649–55. doi: 10.1093/nar/gkx1132
18. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
19. Geistlinger L, Csaba G, Küffner R, Mulder N, Zimmer R. From sets to graphs: towards a realistic enrichment analysis of transcriptomic systems. Bioinformatics. (2011) 27:366–73. doi: 10.1093/bioinformatics/btr228
20. Shichita T, Ito M, Yoshimura A. Post-ischemic inflammation regulates neural damage and protection. Front Cell Neurosci. (2014) 8:319. doi: 10.3389/fncel.2014.00319
21. 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:636–44. doi: 10.1111/ijs.12190
22. 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:74–80. doi: 10.1111/j.1747-4949.2011.00730.x
23. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. (2019) 47:D442–D50. doi: 10.1093/nar/gky1106
24. von Zychlinski A, Williams M, McCormick S, Kleffmann T. Absolute quantification of apolipoproteins and associated proteins on human plasma lipoproteins. J Proteomics. (2014) 106:181–90. doi: 10.1016/j.jprot.2014.04.030
25. Alawieh A, Elvington A, Tomlinson S. Complement in the homeostatic and ischemic brain. Front Immunol. (2015) 6:417. doi: 10.3389/fimmu.2015.00417
26. Alawieh A, Langley EF, Tomlinson S. Targeted complement inhibition salvages stressed neurons and inhibits neuroinflammation after stroke in mice. Sci Transl Med. (2018) 10:eaao6459. doi: 10.1126/scitranslmed.aao6459
27. Di Napoli M. Systemic complement activation in ischemic stroke. Stroke. (2001) 32:1443. doi: 10.1161/01.STR.32.6.1443-a
28. Stokowska A, Atkins AL, Morán J, Pekny T, Bulmer L, Pascoe MC, et al. Complement peptide C3a stimulates neural plasticity after experimental brain ischaemia. Brain. (2017) 140:353–69. doi: 10.1093/brain/aww314
29. Stokowska A, Pekna M. Complement C3a: shaping the plasticity of the post-stroke brain. In: Lapchak PA, Zhang JH, editors. Cellular and Molecular Approaches to Regeneration and Repair. Cham: Springer International Publishing (2018). p. 521–41. doi: 10.1007/978-3-319-66679-2_26
30. Alawieh A, Tomlinson S. Injury site-specific targeting of complement inhibitors for treating stroke. Immunol Rev. (2016) 274:270–80. doi: 10.1111/imr.12470
31. Stokowska A, Olsson S, Holmegaard L, Jood K, Blomstrand C, Jern C, et al. Plasma C3 and C3a levels in cryptogenic and large-vessel disease stroke: associations with outcome. Cerebrovasc Dis. (2011) 32:114–22. doi: 10.1159/000328238
32. Zhang B, Yang N, Gao C. Is plasma C3 and C4 levels useful in young cerebral ischemic stroke patients? Associations with prognosis at 3 months. J Thromb Thrombolysis. (2015) 39:209–14. doi: 10.1007/s11239-014-1100-7
33. Jacob A, Alexander JJ. Complement and blood–brain barrier integrity. Mol Immunol. (2014) 61:149–52. doi: 10.1016/j.molimm.2014.06.039
34. Veerhuis R, Nielsen HM, Tenner AJ. Complement in the brain. Mol Immunol. (2011) 48:1592–603. doi: 10.1016/j.molimm.2011.04.003
35. Gasque P, Fontaine M, Morgan BP. Complement expression in human brain. biosynthesis of terminal pathway components and regulators in human glial cells and cell lines. J Immunol. (1995) 154:4726–33.
36. Noris M, Remuzzi G. Overview of complement activation and regulation. Semin Nephrol. (2013) 33:479–92. doi: 10.1016/j.semnephrol.2013.08.001
37. Orsini F, Chrysanthou E, Dudler T, Cummings WJ, Takahashi M, Fujita T, et al. Mannan binding lectin-associated serine protease-2 (MASP-2) critically contributes to post-ischemic brain injury independent of MASP-1. J Neuroinflamm. (2016) 13:213. doi: 10.1186/s12974-016-0684-6
38. Cai S, Dole VS, Bergmeier W, Scafidi J, Feng H, Wagner DD, et al. A direct role for C1 inhibitor in regulation of leukocyte adhesion. J Immunol. (2005) 174:6462–6. doi: 10.4049/jimmunol.174.10.6462
39. Davis A. Biological effects of C1 inhibitor. Drug News Perspect. (2004) 17:439–46. doi: 10.1358/dnp.2004.17.7.863703
40. Cicardi M, Zingale L, Zanichelli A, Pappalardo E, Cicardi B. C1 inhibitor: molecular and clinical aspects. Springer Semin Immunopathol. (2005) 27:286–98. doi: 10.1007/s00281-005-0001-4
41. Heydenreich N, Nolte MW, Göb E, Langhauser F, Hofmeister M, Kraft P, et al. C1-inhibitor protects from brain ischemia-reperfusion injury by combined antiinflammatory and antithrombotic mechanisms. Stroke. (2012) 43:2457–67. doi: 10.1161/STROKEAHA.112.660340
42. Albert-Weissenberger C, Mencl S, Schuhmann MK, Salur I, Göb E, Langhauser F, et al. C1-Inhibitor protects from focal brain trauma in a cortical cryolesion mice model by reducing thrombo-inflammation. Front Cell Neurosci. (2014) 8:269. doi: 10.3389/fncel.2014.00269
43. Hertle E, Arts ICW, van der Kallen CJH, Feskens EJM, Schalkwijk CG, Stehouwer CDA, et al. Classical pathway of complement activation: longitudinal associations of c1q and C1-INH with cardiovascular outcomes. Arterioscler Thromb Vasc Biol. (2018) 38:1242–4. doi: 10.1161/ATVBAHA.118.310806
44. Kamel H, Iadecola C. Brain-immune interactions and ischemic stroke: clinical implications. Arch Neurol. (2012) 69:576–81. doi: 10.1001/archneurol.2011.3590
45. Carter AM. Complement activation: an emerging player in the pathogenesis of cardiovascular disease. Scientifica. (2012) 2012:14. doi: 10.6064/2012/402783
46. Cojocaru IM, Cojocaru M, Tanasescu R, Burcin C, Atanasiu AN, Petrescu AM, et al. Changes in plasma levels of complement in patients with acute ischemic stroke. Rom J Intern Med. (2008) 46:77–80.
47. Mocco J, Wilson DA, Komotar RJ, Sughrue ME, Coates K, Sacco RL, et al. Alterations in plasma complement levels following human ischemic stroke. Neurosurgery. (2006) 59:1–6. doi: 10.1227/01.NEU.0000219221.14280.65
48. Kouser L, Abdul-Aziz M, Nayak A, Stover CM, Sim RB, Kishore U. Properdin and factor H: opposing players on the alternative complement pathway “See-Saw”. Front Immunol. (2013) 4:93. doi: 10.3389/fimmu.2013.00093
49. Mak TW, Saunders ME. Complement. The Immune Response. Burlington: Academic Press (2006). p. 553–81. doi: 10.1016/B978-012088451-3.50021-1
50. Huber-Lang M, Sarma JV, Zetoune FS, Rittirsch D, Neff TA, McGuire SR, et al. Generation of C5a in the absence of C3: a new complement activation pathway. Nat Med. (2006) 12:682. doi: 10.1038/nm1419
51. Foley JH, Walton BL, Aleman MM, O'Byrne AM, Lei V, Harrasser M, et al. Complement activation in arterial and venous thrombosis is mediated by plasmin. EBioMedicine. (2016) 5:175–82. doi: 10.1016/j.ebiom.2016.02.011
52. Hadders Michael A, Bubeck D, Roversi P, Hakobyan S, Forneris F, Morgan BP, et al. Assembly and regulation of the membrane attack complex based on structures of C5b6 and sC5b9. Cell Rep. (2012) 1:200–7. doi: 10.1016/j.celrep.2012.02.003
53. Ricklin D, Reis ES, Lambris JD. Complement in disease: a defence system turning offensive. Nat Rev Nephrol. (2016) 12:383–401. doi: 10.1038/nrneph.2016.70
54. Parsons ES, Stanley GJ, Pyne ALB, Hodel AW, Nievergelt AP, Menny A, et al. Single-molecule kinetics of pore assembly by the membrane attack complex. Nat Commun. (2019) 10:2066. doi: 10.1038/s41467-019-10058-7
55. Morikis D, Lambris JD. Structural Biology of the Complement System. Boca Raton, FL: CRC Press (2005). doi: 10.1201/9780849350368
56. Serna M, Giles JL, Morgan BP, Bubeck D. Structural basis of complement membrane attack complex formation. Nat Commun. (2016) 7:10587. doi: 10.1038/ncomms10587
57. Portt L, Norman G, Clapp C, Greenwood M, Greenwood MT. Anti-apoptosis and cell survival: a review. Biochim Biophys Acta. (2011) 1813:238–59. doi: 10.1016/j.bbamcr.2010.10.010
58. Westendorp WF, Nederkoorn PJ, Vermeij JD, Dijkgraaf MG, de Beek Dv. Post-stroke infection: a systematic review and meta-analysis. BMC Neurol. (2011) 11:110. doi: 10.1186/1471-2377-11-110
59. Armstrong JR, Mosher BD. Aspiration pneumonia after stroke: intervention and prevention. Neurohospitalist. (2011) 1:85–93. doi: 10.1177/1941875210395775
60. Arachchillage DRJ, Mackie IJ, Efthymiou M, Chitolie A, Hunt BJ, Isenberg DA, et al. Rivaroxaban limits complement activation compared with warfarin in antiphospholipid syndrome patients with venous thromboembolism. J Thromb Haemost. (2016) 14:2177–86. doi: 10.1111/jth.13475
61. Keragala CB, Draxler DF, McQuilten ZK, Medcalf RL. Haemostasis and innate immunity – a complementary relationship. Br J Haematol. (2018) 180:782–98. doi: 10.1111/bjh.15062
62. Maes E, Cho WC, Baggerman G. Translating clinical proteomics: the importance of study design. Expert Rev Proteomics. (2015) 12:217–9. doi: 10.1586/14789450.2015.1041512
63. Dayon L, Kussmann M. Proteomics of human plasma: a critical comparison of analytical workflows in terms of effort, throughput and outcome. EuPA Open Proteom. (2013) 1:8–16. doi: 10.1016/j.euprot.2013.08.001
64. Cao Z, Yende S, Kellum JA, Robinson RA. Additions to the human plasma proteome via a tandem MARS depletion iTRAQ-based workflow. Int J Proteomics. (2013) 2013:8. doi: 10.1155/2013/654356
65. Tu C, Rudnick PA, Martinez MY, Cheek KL, Stein SE, Slebos RJC, et al. Depletion of abundant plasma proteins and limitations of plasma proteomics. J Proteome Res. (2010) 9:4982–91. doi: 10.1021/pr100646w
66. Smith MPW, Wood SL, Zougman A, Ho JTC, Peng J, Jackson D, et al. A systematic analysis of the effects of increasing degrees of serum immunodepletion in terms of depth of coverage and other key aspects in top-down and bottom-up proteomic analyses. Proteomics. (2011) 11:2222–35. doi: 10.1002/pmic.201100005
67. Faria SS, Morris CFM, Silva AR, Fonseca MP, Forget P, Castro MS, et al. A timely shift from shotgun to targeted proteomics and how it can be groundbreaking for cancer research. Front Oncol. (2017) 7:13. doi: 10.3389/fonc.2017.00013
69. Whiteaker JR, Zhao L, Schoenherr RM, Kennedy JJ, Ivey RG, Paulovich AG. Peptide immunoaffinity enrichment with targeted mass spectrometry: application to quantification of ATM kinase phospho-signaling. Methods Mol Biol. (2017) 1599:197–213. doi: 10.1007/978-1-4939-6955-5_15
70. Shi K, Wood K, Shi FD, Wang X, Liu Q. Stroke-induced immunosuppression and poststroke infection. Stroke Vasc Neurol. (2018) 3:34. doi: 10.1136/svn-2017-000123
71. Stokum JA, Gerzanich V, Simard JM. Molecular pathophysiology of cerebral edema. J Cereb Blood Flow Metab. (2016) 36:513–38. doi: 10.1177/0271678X15617172
Keywords: longitudinal, stroke, proteomics, immune system, complement system, bioinformatics, systems biology, pathway analysis
Citation: Nguyen VA, Riddell N, Crewther SG, Faou P, Rajapaksha H, Howells DW, Hankey GJ, Wijeratne T, Ma H, Davis S, Donnan GA and Carey LM (2020) Longitudinal Stroke Recovery Associated With Dysregulation of Complement System—A Proteomics Pathway Analysis. Front. Neurol. 11:692. doi: 10.3389/fneur.2020.00692
Received: 18 April 2020; Accepted: 09 June 2020;
Published: 28 July 2020.
Edited by:
Aristeidis H. Katsanos, McMaster University, CanadaReviewed by:
Marios K. Georgakis, LMU Munich University Hospital, GermanyMichael Chong, Population Health Research Institute (PHRI), Canada
Copyright © 2020 Nguyen, Riddell, Crewther, Faou, Rajapaksha, Howells, Hankey, Wijeratne, Ma, Davis, Donnan and Carey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vinh A. Nguyen, di5uZ3V5ZW4mI3gwMDA0MDtsYXRyb2JlLmVkdS5hdQ==