Participants
Participants were recruited from the Transitional Living Center at Casa Colina Hospital and Centers for Healthcare, Pomona, California. To be eligible, the patients had to 1) be between 18 and 65 years old, 2) have suffered a moderate to severe TBI (GCS 9-12 and GCS 3-8, respectively) and be at least one year after injury post-injury at the time of enrollment (i.e., more than a year post-injury but less than two years post-injury), 3) have grossly intact motor use of dominant hand and no gross receptive language impairment/aphasia. Participants were excluded if 1) they were neither speakers of English or Spanish, 2) they had a diagnosis other than TBI (e.g., spinal cord injury, other primary neurologic condition), a developmental disorder (i.e., down syndrome), or a pre-morbid psychiatric condition resulting in admission to a hospital. Participants with claustrophobia or with metallic implants making it unsafe to enter the MRI environment were also excluded. Study protocols were approved by the institutional review board of Casa Colina Hospital and Centers for Healthcare. Participants provided written informed consent.
Data acquisition and analyses
In this prospective longitudinal study, the following data were collected, at monthly intervals (every 4 to 6 week) across 6 time-points (study duration: 6.43 ± 0.84 months):
Data collection. Single whole blood samples were collected from participants. After collection of the whole blood, the blood was allowed to clot by leaving it undisturbed at room temperature for 15 minutes. The clot was removed by centrifuging at 1,000–2,000 x g for 10 minutes. The resulting supernatant was designated serum. Following centrifugation, the serum was aliquoted directly to microcentifuge tubes (500ul). The samples were maintained at 2–8°C while handling and stored at –20°C until being transferred on dry ice to a -80°C freezer for long term storage and processing at the Biology Department of the University of La Verne.
Data analyses. A quantitative sandwich enzyme-linked immunosorbent assay (ELISA) technique was used to analyze serum levels of glial fibrillary acidic protein (GFAP) (Millipore, Burlington, MA), neuron-specific enolase (NSE) (R&D Systems, Minneapolis, MN), human alpha II-spectrin breakdown product (SBDP145) (MyBioSource, San Diego, CA), S100 calcium-binding protein A12 (S100A12) (Millipore, Burlington, MA), Ubiquitin C-terminal hydrolase-L1 (UCH-L1) (MyBioSource, San Diego, CA) as well as T-tau (Abcam, Cambridge, UK), P-tau (MyBioSource, San Diego, CA), and P-tau/T-tau ratio. Samples for each participant were run in duplicate on the same plate. Protein levels were determined based on the standard curve using the average of the duplicate values.
The behavioral assessment was performed (in English or Spanish) by trained clinicians and included cognitive measures (i.e., the Montreal Cognitive Assessment and subtests of the Repeatable Battery for the Assessment of Neuropsychological Status), a mood measure (i.e., the Beck Depression Inventory) and a measure of functional outcome (i.e., the Disability Rating Sale).
Montreal Cognitive Assessment (MOCA). The MOCA is a widely used screening assessment for detecting gross cognitive impairment in a variety of diseases including traumatic brain injuries. It assesses short-term and long-term memory, attention/working memory, executive functions, language and orientation to time and space. Its total score ranges from 0 to 30. It exists in around 50 languages and it has 3 alternative versions. In this study, either the English or the Spanish version was used.[10-11] The alternative versions were used twice across assessments (in the following order: 1, 2, 3, 1, 2, 3).
Repeatable Battery for the Assessment of Neuropsychological Status (RBANS). The RBANS is a neuropsychological assessment which tests five cognitive domains (short-term and long-term memory, visuospatial/constructional abilities, language, attention/working memory). It was originally introduced in the screening for dementia, but has also found application to traumatic brain injury.[12] It has four alternative versions but, since these alternative versions do not exist in Spanish, we decided to administer the following non-verbal subtests to both English and Spanish speakers: digit span (DS) (score range: 0-16), figure copy (FC) (score range: 0-20), and coding (C) (score range: 0-89). The alternative versions were administered in the following order across assessments: a, b, c, d, a, b.
Beck Depression Inventory (BDI-II). The BDI-II is a 21-question multiple-choice self-report inventory, and is one of the most widely used psychometric tests for measuring depression. It includes items related to hopelessness, irritability, guilt, as well as fatigue, weight loss, and lack of sex drive. Its total score ranges from 0 to 63. Both its English and Spanish version were used.[13-14]
Disability Rating Scale (DRS). The DRS was developed to assess functional recovery of patients after a traumatic brain injury.[15] It assesses eye opening, communication ability, feeding, toileting, grooming, level of functioning and employability. It does not request an interview with the patient but it relies on the rater’s behavioral observation and knowledge of the patient. Its total score ranges from 29 (worst outcome) to 0 (best outcome).
- Magnetic resonance imaging (MRI)
Data collection. MRI data were collected on a Siemens Magnetom Verio 3T system at the Casa Colina Diagnostic Imaging Center. Each participant underwent a conventional structural MRI T1-weighted 3-dimensional magnetic-preparation rapid gradient echo scan (MPRAGE) with the following parameters: TR = 2300 ms, TE = 2 ms, flip angle = 9°, FOV = 23 cm, slice thickness = 1 mm with no gap, number of slices = 160, matrix size = 224 x 224.
Subcortical volumes were calculated using FSL's FIRST.[19] Subcortical morphometry (shape statistics) measurements were calculated using the Metric Optimization for Computational Anatomy toolbox.[20] We modified FSL’s vertex analysis to allow for longitudinal analysis within subject across multiple timepoints and regressed against the independent variables (i.e., the blood-based biomarkers).[19] First, 15 subcortical structures were segmented for each participant for each time point. This novel workflow creates a cohort average subcortex over all time points and participants, as well as an image per subcortical region of interest that represents the localized longitudinal shape change within each subject compared with the cohort. Group-level significance was assessed using a general linear model with a non-parametric permutation test at a level of p < 0.05 corrected for multiple comparisons using family-wise cluster correction and threshold free cluster enhancement (TFCE) as implemented in FSL randomize.[21]
- Resting state electroencephalogram (EEG)
Data collection. EEG-data were acquired at the Casa Colina research institute in a faraday room. A 24 electrode skull cap was connected to a wireless portable digital EEG amplifier (B-Alert wireless EEG system; advancedbrainmonitoring.com). Fp1 and Fp2 were used to monitor eye artifacts. A mastoid reference was applied. EEG-data were recorded on a laptop computer. The impedances were kept below 40 kW. A sampling rate of 250 Hz was used. EEG recordings were performed while the participants were in a wakeful state with eyes open, sitting in a setting with minimal ambient noise. Patients were instructed to relax and look at a fixed object in the room. Duration of recording was 5 minutes.
Data analyses. EEG data were analyzed using the GUI (graphic user interface) of EEGLAB (version 14.0.0b). Preprocessing included channel location, band-pass filtered between 1 and 40 Hz (FIR filter), Independent Component Analyses (ICA; with jader decomposition algorithm]) to identify and exclude ocular and motor artifacts related components (maximum of 3), and a final visual inspection of the EEG data in order to manually remove any residual artifacts. For the spectral analysis, the continuous signal of each recording was segmented into 1s epochs. Each epoch was fast Fourier transformed (FFT) with a Hanning-tapered window. Subsequently, all epochs were averaged, and the mean of the power spectral density in different frequency bands were exported for statistical analysis. Frequency bands were chosen at target electrode sites based on previous literature pinpointing sites of maximal amplitude for each band: alpha (8-13 Hz) at the occipital electrodes (O1, Oz, O2), theta (4-7 Hz) at the frontal electrodes (F7, F3, Fz, F4, F8), and delta (1-3 Hz) at the central electrodes (C3, Cz, C4). The averaged power spectral density extracted for each site (F, C, O) was expressed in absolute power (μV2).[22]
Statistical analyses
For blood, behavioral and EEG data, LME analyses were performed on our 6 timepoints (repeated effect) per patient and 48 observations (repeated covariance type: diagonal, severity of lesions as a fixed factor).[23] First, time since injury was covaried (as a fixed factor) to each blood-based biomarker (i.e., GFAP, NSE, SDBP145, UCH-L1, T-tau, P-tau, and P-tau/T-tau ratio). Time since injury was adjusted across patients by subtracting the mean of the first assessment from each observation for this variable. Then, the blood-based biomarkers showing significant changes were covaried with 1) the outcome measures (as fixed factors: MOCA, RBANS DS, RBANS FC, RBANS C, BDI, DRS) (Bonferroni correction applied: p≤.008), and 2) the resting state EEG (as fixed factors: alpha, delta, theta) (Bonferroni correction applied: p≤.02). Missing data were discarded in our analyses. Values for NSE and GFAP were not obtained in 2 out of 8 patients (due to deficient ELISA kits) and analyses were performed on 6 patients. Statistical analyses performed for the MRI data have been described above [see “Magnetic resonance imaging (MRI)” section].We only report results that remain significant after corrections for multiple comparisons.