This was a prospective, observational study, part of two separate studies conducted at the Karolinska University Hospital, and Karolinska Institutet, Stockholm, Sweden. The first study included TBI patients between 2007 and 2015. Oral informed consent was granted by next-of-kin. The second study included healthy volunteers, used as control subjects here, between 2014 and 2015. All control subjects provided written, informed consent. All research activities were in accordance with Swedish law and the Declaration of Helsinki. Ethical approvals (#2005/1526-31/2; #2014/1201-31/1) were granted by the Swedish Ethical Review Authority.
Study participant inclusion and exclusion criteria
Inclusion criteria for TBI patients were: (i) severe TBI (as per Glasgow Coma Scale [GCS] 3-8 upon hospital admission or else a higher GCS score but with a significant risk for deterioration) in need of NCCU treatment and invasive intracranial monitoring, and (ii) age 18-75 years. Exclusion criteria comprised: (i) desolate prognosis precluding NCCU treatment, (ii) penetrating TBI, (iii) unconsciousness due to etiology other than TBI, (iv) underlying chronical condition precluding follow-up, or (v) other reason precluding follow-up. Inclusion criteria for control subjects were: (i) previously healthy, (ii) age 18-50 years, (iii) sufficient linguistic knowledge to participate in self-evaluation forms. Exclusion criteria were: (i) ongoing, or history of, psychiatric illness, (ii) family history of serious psychiatric comorbidity, (iii) somatic illness precluding physical activity, (iv) current pharmacological treatment interacting with the study intervention, (v) substance abuse (smoking or narcotic substances), or (vi) pregnancy. Sample size calculation was based on expected protein level difference between TBI patients and control subjects and was exerted as a two-sample t-test. We utilized Cohen’s d (30,31) as effect size metric and set it to 0.8 (large effect) (30,31) in a power calculation utilizing the R package pwr (32). In order to obtain 80% power at the 0.05 significance level with n = 15 control patients, we needed to recruit n = 77 TBI patients. As this was not based on empirical data, we included patients continuously throughout the study period.
Clinical management, data, and sample acquisition
NCCU management of severe TBI at Karolinska University Hospital has been described elsewhere (33). In brief, Karolinska University Hospital employs an intracranial pressure (ICP-) driven approach, in accordance with the Brain Trauma Foundation Guidelines (34). ICP is monitored either through a closed external ventricular drain (EVD) (Medtronic, USA), or an intraparenchymal pressure monitor (Codman & Shurtleff Inc. Raynham, MA, USA or Rehau AG + CO, Rehay, Germany). While EVDs may be used to drain CSF in order to decrease ICP, the choice between monitoring device is multifactorial and not exclusively reliant on injury severity. At the NCCU, multi-modal monitoring data is automatically collected. Through the Karolinska University Hospital TBI Database, additional data is collected prospectively and comprise neurological variables, injury severity score variables, radiological variables, and outcome data, described in detail elsewhere (11). Functional outcome data (Glasgow Outcome Score, GOS) was collected at 6-12 months following hospital discharge, through structured questionnaires, or follow-up assessments in the outpatient clinic at the Neurosurgical Department. We collected CSF and serum, used for APOE genotyping, proteomic, and albumin analysis. The latter was assessed as QA, i.e. the CSF/serum albumin quotient (10), with the reference intervals (35): 15-29 years < 0.006; 30-49 years < 0.007; and ³ 50 years < 0.009. Sampling time points were not identical for albuminCSF, albuminserum and the proteomic samples from CSF and serum. Time discrepancies were in median (interquartile range [IQR]): 4.3 (0-11.8) hours for albuminCSF and albuminserum samples; 0.88 (-2.27-9.15) hours for albuminCSF and the proteomic sample; and -2.83 (-3.82 - -2.08) hours for albuminserum and the proteomic sample.
Sample acquisition
Control subjects were recruited to a study on effects of a physical exercise intervention (36), of which only baseline samples were used. Participants were instructed to abstain from physical exercise seven days before sampling, performed by lumbar puncture and venipuncture, between 7.30 and 9 AM while fasting since midnight after a full night of bed rest. For TBI patients, blood was sampled through an arterial line and CSF through an EVD. TBI sample acquisition occurred in median at 60.8 hours (IQR 36.6-109.1) following trauma for CSF samples and 53.3 hours (30.5-91.1) for serum samples (Figure S1A). Samples were stored locally in 4°C in median 1 day (0-1) for both CSF and serum (Figure S1B), until delivery to a local biobank, where samples were vertically incubated for 30 min before centrifugation for 15 min at 2000g, aliquoting, and storage at -80°C until further analysis (37). No protein content alteration was seen per sample (Figure S2A) or analyte (Figure S2B, representative example) due to delayed biobank delivery.
Genotyping
Whole blood was collected together with serum in ethylenediaminetetraacetic acid (EDTA) tubes, and was frozen in the biobank until DNA extraction. Genotyping was performed with the SNP markers rs429358 (ApoE112) and rs7412 (ApoE158) using single base primer extension (SBE) with detection of the incorporated allele by ¨Fluorescent Polarization Template Dye Incorporation¨ (FP-TDI) (38). Signal intensities were read using a Tecan Genios Pro fluorescence absorbance reader. Raw data from the fluorescence polarization was converted to genotype data using the software AlleleCaller 4.0.0.1 and alleles ε2, ε3 or ε4 were identified.
Proteomic analysis
In total, 177 protein depicted through 220 antibodies were examined (Table S1, where the full protein name is provided). For 43 proteins, two antibodies targeted different regions of the same protein, i.e. sibling antibodies (39). The protein panel was chosen based on CNS-enrichment (40), previous clinical/experimental/mass-spectrometry TBI studies, or previous neuroinflammation studies (20,24,26,41–45). Antibodies were selected from the Human Protein Atlas (HPA) (www.proteinatlas.org) (46).
Antibodies were immobilized onto color-coded magnetic beads (MagPlex, Luminex Corporation) as previously described (28). Briefly, the beads surface was activated by using 0.1 M sodium hydrogen phosphate (Sigma), 0.5mg of N-hydroxysulfosuccinimide (sulfo-NHS) (Nordic Biolabs) and 0.5mg 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochloride (EDC) (ProteoChem). Beads were then incubated with antibodies (16 μg/ml in 2-(N-morpholino)ethanesulfonic acid [MES] buffer, Sigma) for 2h at room temperature. Each antibody type was immobilized on a different bead identity (bead type with specific color-code). After incubation, the beads were washed with phosphate-buffered saline (PBS, Fisher Scientific) 0.05% Tween-20 (Fisher Scientific) (PBS-T) to eliminate the antibody excess, stored overnight in blocking buffer (Roche blocking reagent for ELISA, Roche), and combined into a suspension bead array.
Samples were processed as previously described, with minor adjustments (27,47). Serum and CSF samples were separately randomized into 96-well microtiter plates. CSF samples were diluted 0.6:1 in PBS (Fisher Scientific) with 0.5% bovine serum albumin (BSA, Sigma), 0.1% rabbit IgG (Nordic Biosite), and labeled with biotin (Fisher Scientific). The samples were then further diluted 1:8 in assay buffer (0.1% casein [Fisher Scientific], 0.5% polyvinyl alcohol [Sigma], 0.8% polyvinylpyrrolidone [Sigma] in PBS-T (0.05% Tween-20 [Fisher Scientific]), supplemented with 0.5 mg/ml rabbit IgG [Nordic Biosite]), heat treated (56◦C for 30min), and incubated with the bead array overnight at room temperature. Serum samples were diluted 1:10 in PBS (Fisher Scientific) prior to labeling with biotin (Fisher Scientific), and further diluted 1:50 in assay buffer (0.1% casein [Fisher Scientific], 0.5% polyvinyl alcohol [Sigma], 0.8% polyvinylpyrrolidone [Sigma] in PBS-T (0.05% Tween-20 [Fisher Scientific]), supplemented with 0.5 mg/ml rabbit IgG [Nordic Biosite]) after labeling, heat treated (56◦C for 30min), and incubated with the bead array for 2 hours at room temperature.
The captured proteins were cross-linked to the antibodies for 10 min at room temperature using 0.4% paraformaldehyde (Thermo Scientific). The antibody-protein immunocomplexes were detected by using a streptavidin-conjugated phycoerythrine (Fisher Scientific) and a FlexMap3D instrument (Luminex Corporation). The relative protein abundance was reported as median fluorescence intensity (MFI) for each bead identity and sample. Quality control assessments are described in Supplementary Methods. Briefly, bead counts were evaluated per sample and analyte (Figure S3A-S3B). Due to a small systematic increase in MFICSF samples (Figure S4A), background subtraction was conducted (Figure S4B). MFI values varied across analytes (Figure S4C), of which one was excluded due to borderline non-detected signal (Figure S4C, inset). Antigen profiles were assessed per sample and analyte (Figure S5-S6, Table S2), resulting in the exclusion of a few sibling antibodies (Supplementary Results).
Statistical analysis
For inferential analysis, matched CSF-serum patient samples were compared. Validation analysis was exerted in the non-matched TBI cohort with serum-samples only. We used R (version 4.0.2) (48), through RStudio® (version 1.3.1056) and the tidyverse (49), RColorBrewer (50), cowplot (51), and gridExtra (52) packages. Continuous data were presented as median (IQR). Categorical data were presented as count (%). For multiple testing correction, we used the Bonferroni, Holm (53) or the false-discovery rate (FDR) (54) method. A p-value < 0.05 was considered significant, unless otherwise stated.
A few variables (pre-hospital hypotension, QA, APOE allele status) had a substantial number of missing values (Table 1, Figure S7). When applicable, we conducted multiple imputation using n = 200 imputations in the mice package (55). Reported p-values were calculated as the unadjusted median p-value from all imputations.
Protein Characterization
Analytes were characterized using the HPA (46,56) version 19.1 (release date 2019/12/19, Ensembl version 92.38), using the protein tissue data, RNA tissue data (Consensus data set), and Brain Atlas (57) RNA data (Supplementary Methods, Supplementary Results).
Parallel assessments in CSF, serum, and relationship with BBB disruption
T-distributed stochastic neighbor embedding (t-SNE) (58,59) was employed to examine if proteins pertained to compartment (CSF or blood) and disease characteristics among study subjects (Supplementary Methods). We assessed protein levels in CSF and serum under control conditions and following TBI using the Wilcoxon rank sum test (FDR, padjusted < 0.05) and the Wilcoxon signed rank test (FDR, padjusted < 0.01).
Cluster analysis within CSF and serum was conducted for proteins that had a CSF/serum ratio significantly correlated (Kendall correlation, Holm method, padjusted < 0.05) with QA (Supplementary Methods). Clusters were visualized using the ComplexHeatmap package (60). Proteins significant upon linear regression (FDR, padjusted ≤ 0.01) compared with the reference cluster (containing the majority of control patients) were deemed significantly altered. For CSF (n = 3 clusters), proteins needed to be concurrently significant in all clusters compared with the reference cluster. Protein levels between TBI patients with disrupted/intact BBB were compared using the Wilcoxon Rank Sum Test (FDR, padjusted < 0.05). Linear regression models were used to examine if APOE4 carriership was important for QA, or protein levels (FDR, p £ 0.05). Age, gender and injury scores were used as covariates in addition to APOE variant.
Pathway and outcome analysis
Pathway analysis through the pathfindR package (61) and pipeline (62), was conducted for proteins altered following TBI or that pertained to a BBB integrity related cluster. For protein input, p-value thresholds were set to 0.05. For enrichment analyses, the Biocarta gene set and the Bonferroni method (padjusted ≤ 0.05) for multiple correction were used.
Proteins of interest for outcome analysis were: i) protein intersects between CSF cluster analysis and TBI-induced altered proteins in CSF, ii) protein intersects between CSF cluster analysis and TBI-induced altered proteins in serum, and iii) significantly elevated/decreased proteins following BBB disruption. Protein intersects were visualized using the VennDiagram package (63) in R. We used GOS as dependent variable and protein levels of an individual protein (or other variable of interest such as QA) as independent variable in a proportional odds regression analysis, using the rms package (64). Only TBI patients were included, as healthy control subjects by definition had no GOS data. We conducted univariable analysis, and if significant (FDR, padjusted £ 0.05 or ≤ 0.01 if multiple testing, the latter for dichotomized GOS/short-term mortality), multivariable analysis (FDR, padjusted < 0.05 if multiple testing or pimputed ≤ 0.05 if imputed). We used age, GCS motor score, pupillary reactions, hypoxia, hypotension and Stockholm computerized tomography (CT) score as covariates in accordance with the International Mission for Prognosis and Clinical Trial (IMPACT) database studies (65). We used the Stockholm instead of the Marshall CT score, as the former has been shown to be superior (66,67). When applicable, we conducted step-down modelling to see how the proteins performed jointly in the regression models.