Changes in synaptic density in the subacute phase after ischemic stroke: A 11 C-UCB-J PET/MR study

Abstract

Functional alterations after ischemic stroke have been described with Magnetic Resonance Imaging (MRI) and perfusion Positron Emission Tomography (PET), but no data on in vivo synaptic changes exist. Recently, imaging of synaptic density became available by targeting synaptic vesicle protein 2 A, a protein ubiquitously expressed in all presynaptic nerve terminals. We hypothesized that in subacute ischemic stroke loss of synaptic density can be evaluated with ¹¹C-UCB-J PET in the ischemic tissue and that alterations in synaptic density can be present in brain regions beyond the ischemic core. We recruited ischemic stroke patients to undergo ¹¹C-UCB-J PET/MR imaging 21 ± 8 days after stroke onset to investigate regional ¹¹C-UCB-J SUVR (standardized uptake value ratio). There was a decrease (but residual signal) of ¹¹C-UCB-J SUVR within the lesion of 16 stroke patients compared to 40 healthy controls (ratio_{lesion/controls} = 0.67 ± 0.28, p = 0.00023). Moreover, ¹¹C-UCB-J SUVR was lower in the non-lesioned tissue of the affected hemisphere compared to the unaffected hemisphere (ΔSUVR = −0.17, p = 0.0035). The contralesional cerebellar hemisphere showed a lower ¹¹C-UCB-J SUVR compared to the ipsilesional cerebellar hemisphere (ΔSUVR = −0.14, p = 0.0048). In 8 out of 16 patients, the asymmetry index suggested crossed cerebellar diaschisis. Future research is required to longitudinally study these changes in synaptic density and their association with outcome.

Introduction

Predicting outcome in patients with subacute ischemic stroke is challenging. Most of the developed prediction tools focus on motor recovery and use baseline characteristics obtained in the acute phase that are associated with outcome. Age and stroke severity, measured by a variety of clinical scales, are the most frequently used clinical characteristics and a number of Magnetic Resonance Imaging (MRI) parameters and results of Transcranial Magnetic Stimulation (TMS) have been investigated. MRI and TMS are mostly used to assess the integrity of the corticospinal tract and if the corticospinal tract is too severely damaged, the proportional recovery rule (which states that 3 months after stroke, patients will have regained about 70% of their maximum recovery potential) is not valid anymore. ¹ For most of these prediction tools, however, external validation and implementation in clinical practice is lacking.^2–4 Therefore, there is a need to identify subclinical, preferentially modifiable biomarkers of motor recovery during the acute and subacute phase to improve prognostication and potentially evaluate the response to treatment aimed at improving recovery. ⁵ However, as long as the processes underlying recovery remain obscure, identification of such a biomarker is hard. Preclinical studies investigating the mechanisms underlying stroke recovery identified changes of synaptic density in the peri-ischemic area and remote regions, suggesting plasticity of the brain.^6–8 Functional MRI and perfusion Positron Emission Tomography (PET) in humans during the recovery period revealed alterations in activity with evidence of both ipsilesional as well as contralesional and bilateral cerebral (over)activation. ⁹ Functional changes also occur in the cerebellum after supratentorial stroke, known as crossed cerebellar diaschisis. Crossed cerebellar diaschisis is characterized by a loss of function of the contralesional cerebellum due to disruption of cerebro-cerebellar pathways.^10–15 Recently, imaging of synaptic density has become available by targeting synaptic vesicle protein 2 A (SV2A), a protein ubiquitously expressed in presynaptic nerve terminals.^16–19 This novel PET tracer offers a unique opportunity to gain insight in the spatial and temporal dynamics of synaptic alterations after ischemic stroke in humans, which may be a key factor in unraveling the processes underlying recovery.

In this study, we investigated the feasibility of SV2A PET imaging in patients with subacute ischemic stroke. We hypothesized that ischemic brain tissue is characterized by loss of synaptic density.^20,21 Furthermore, we investigated if SV2A imaging could also identify alterations in synaptic density beyond the ischemic core possibly in the context of plasticity or injury beyond the ischemic lesion.^6–9 In addition, we explored if crossed cerebellar diaschisis could be detected by SV2A PET as disruption of cerebro-cerebellar connections could translate into decrease of presynaptic vesicles.

Methods

This study was approved by the local University Hospital Ethics committee (UZ Leuven/KU Leuven) and was conducted in accordance with the latest version of the Declaration of Helsinki. All participants (or relatives in case the patient was unable to give written informed consent) signed written informed consent prior to participation. This study is part of a larger, multimodal, multitracer, longitudinal imaging study with different patient populations (ClinicalTrials.gov: NCT03514524).

Participants

We recruited patients between 18–85 years with a first ischemic stroke. Inclusion criteria were: first ischemic stroke within 2 weeks prior to inclusion; significant motor deficit of the upper extremity, defined as MRC (Medical Research Council) ²² grade 3 or less for finger extension and/or elbow flexion and/or shoulder abduction and sufficient cooperation/comprehension to cooperate with clinical evaluation. Major exclusion criteria were: pre-stroke Barthel Index <19 ²³ ; history of traumatic brain injury, Alzheimer disease or other important neurological disorders; contra-indications for MRI; use of levetiracetam/brivaracetam and/or previous participation in other research studies involving ionizing radiation with >1 mSv over the past 12 months. Patients underwent an extensive test battery at inclusion, including NIHSS (National Institutes of Health Stroke Scale),^24,25 Fugl-Meyer upper extremity (FMUE) and lower extremity (FMLE) ²⁶ and Barthel Index (BI). ²³ Patients underwent ¹¹C-UCB-J PET/MR imaging within the first month after the event.

We also included healthy controls between 18–85 years, with main exclusion criteria: history of important neurological, psychiatric or internal medical disorder; chronic use of medication with effects on the central nervous system; pregnancy; history or active use of alcohol or other drug abuse; contra-indications for MR; brain white matter disease Fazekas 3 (for subjects ≥60 years) or other relevant MR abnormalities and previous participation in other research studies involving ionizing radiation with >1 mSv over the past 12 months.

Image acquisition and reconstruction

All participants underwent simultaneous ¹¹C-UCB-J PET/MR imaging 60–90 minutes post tracer injection on a General Electric (GE, Milwaukee, WI, USA) Signa PET/MR (time-of-flight (TOF) PET, 3-Tesla MR) scanner. Median injected ¹¹C-UCB-J activity was 214 MBq (IQR: 204-231 MBq) with a mean specific activity of 244 ± 147 GBq/µmol for patients. Median injected ¹¹C-UCB-J activity for healthy controls was 246 MBq (IQR: 216-294 MBq) with a median specific activity of 222 GBq/µmol (IQR: 159–318 GBq/µmol). PET acquisitions were performed in list mode, rebinned into 6 frames of 5 minutes and reconstructed using vendor software. As reconstruction algorithm, an ordered subset expectation maximization (28 subsets, 4 iterations) was used including TOF information and decay, scatter, attenuation, deadtime and random correction. Isotropic Gaussian post smoothing was performed to reduce the noise (full width at half maximum (FWHM) of 4 mm; final image resolution 5 mm). Individual attenuation correction was performed using a validated zero echo time (ZTE) approach. ²⁷

We acquired 3 D T1-weighted MR images with an 8-channel head coil (plane: sagittal; TE: 3.2 ms; TR: 8.5 ms; TI: 450 ms; flip angle: 12; receiver bandwidth: 31.25 kHz; NEX: 1; voxelsize: 1 × 1 × 1 mm) as well as 3 D T2-weighted FLAIR images (plane: sagittal; TE: 136/137 ms; TR: 8500 ms; TI: 2298 ms; receiver bandwidth: 31.25 kHz; NEX: 1; voxelsize: 0.7 × 1 × 1 mm).

Image processing

We preprocessed the ¹¹C-UCB-J PET images using PMODv3.9 (PMOD technologies, Zurich, Switzerland): first, motion correction was applied by rigidly realigning each frame to the first frame, then averaging all frames. Second, we co-registered the PET images to the corresponding T1-weighted MRI data and calculated standardized uptake value ratios (SUVR) using the centrum semiovale (CS) as reference region (see Supplementary Methods for detailed description). Representative images of 3 patients and 1 healthy control are shown in Figure 1.

Figure 1.

Representative images of different stroke patients and one healthy control, shown in radiological convention. Top row: ¹¹C-UCB-J SUVR images overlaid on T1-weighted MR images; bottom row: T2 FLAIR images showing the lesions of the patients. Patient 1 was a 82-year old female patient with a large ischemic stroke in the right middle cerebral artery territory. Patient 2 was a 78-year old female patient with a subcortical stroke in the left corona radiata, internal capsule and lentiform nucleus. Patient 3, a 51-year old male, suffered an ischemic stroke in the watershed area of the left hemisphere. The last column shows a ¹¹C-UCB-J SUVR image and T2 FLAIR image of a 58-year-old, female healthy control. SUVR = standardized uptake value ratio.

For volume-of-interest (VOI)-based analyses, we delineated VOI’s with FreeSurfer (Laboratory for Computational Neuroimaging v6.0, Boston, USA; cortical parcellation based on the Desikan-Killiany atlas,^28,29 subcortical structures as described in Fischl et al. ³⁰ ). As the large cortical stroke lesions could interfere with proper VOI delineation, we used virtual brain grafting (version 4) for all patients with cortical involvement of the stroke. Virtual brain grafting is a processing pipeline that replaces a delineated lesion by healthy tissue, based on the contralesional hemisphere. This healthy T1-weighted image is then used as input for VOI delineation in FreeSurfer. ³¹

To correct for partial volume effects (PVE), we applied a region-based voxel-wise (RBV) partial volume correction (point spread function 5 mm).^32,33 For optimization of this partial volume correction, the lesion itself was not considered as a homogeneous entity, but subdivision in grey and white matter regions was preserved as SV2A uptake is different in grey and white matter. The peri-ischemic area was not handled as separate region for RBV as this area is not homogeneous either. However, in contrast to the lesion, subdividing the peri-ischemic area in grey and white matter regions would have created too many too small regions resulting in less robust RBV estimation.

We defined the stroke lesions semi-automatically using the 3 D T2-weighted FLAIR data and created a peri-ischemic area by applying a 4 mm simple binary dilation of the lesion mask. We adapted the peri-ischemic area manually to exclude voxels not consisting of brain parenchyma. Representative images of these VOI’s are shown in Figure 2. To be able to compare values of the stroke lesions and peri-ischemic areas with healthy controls, we transformed all VOI’s of the lesions and peri-ischemic areas to MNI space (Montreal Neurological Institute; see Supplementary Methods for more information about the non-linear deformation fields). We created an average ¹¹C-UCB-J SUVR image of the healthy controls in MNI space (=AV HC) to extract an average normal value for all lesions and peri-ischemic areas. Moreover, to compare the values of the lesions and peri-ischemic areas to the homologous areas in the unaffected hemisphere (=mirror regions), we flipped the ¹¹C-UCB-J SUVR images of all patients in MNI space over the sagittal axis and we co-registered them to the original MNI space. We extracted the average values of the lesion/peri-ischemic mirror areas.

Figure 2.

Representative images of four patients with stroke, shown in radiological convention. Top row: Volumes-of-interest, overlaid on T1-weighted MR images. The lesion is visualized in red, the peri-ischemic area in yellow, the non-lesioned tissue of the affected hemisphere in blue and the non-lesioned tissue of the unaffected hemisphere in green. Bottom row: corresponding T2 FLAIR images. VOI = volume-of-interest.

We defined the non-lesioned tissue as all grey matter brain tissue of the affected and unaffected hemisphere that was not part of the ischemic lesion or peri-ischemic area. For each hemisphere, we extracted an average ¹¹C-UCB-J SUVR value of the non-lesioned tissue (see Figure 2 for visualization of the non-lesioned tissue).

Simulation study

To assess the impact of changes in blood flow on SUVR quantification, we performed a simulation study using the ¹¹C-UCB-J input function of one healthy volunteer (historical data), by applying a one-tissue compartment model. ³⁴ We simulated region-specific changes in tracer delivery by altering K₁ and adapting k₂ accordingly to preserve the K₁/k₂ ratio (=V_T (volume of distribution)) constant (PMODv3.9).^35,36 We applied following baseline values: K₁ =0.45 ml.cm⁻³.min⁻¹ and k₂ = 0.015 min⁻¹ (V_T = 30.0 ml.cm⁻³) for the cerebral cortex, K₁ = 0.45 ml.cm⁻³.min⁻¹ and k₂ = 0.020 min⁻¹ (V_T = 22.5 ml.cm⁻³) for the cerebellum and K₁ = 0.123 ml.cm⁻³.min⁻¹ and k₂ = 0.023 min⁻¹ (V_T = 5.3 ml.cm⁻³) for the CS. We used the simulation of the cerebral cortex relative to the CS to assess the impact of perfusion asymmetries in non-lesioned tissue and the simulation of the cerebellum relative to the CS to assess crossed cerebellar diaschisis.

Based on available literature data, we assumed regional perfusion asymmetry in the non-lesioned tissue ranging from −15% to +15%^37–39 and for crossed cerebellar diaschisis, we assumed asymmetry in perfusion ranging from −35% to +15%.^11,40

We assessed the possible role of perfusion changes on ΔSUVR by reporting the results of the simulation and the measured values of the non-lesioned tissue (measured ΔSUVR_non-lesioned) and cerebellar hemispheres (measured ΔSUVR_cerebellum).

$$\begin{array}{l}\mathit{measured}\mathit{\hspace{0.17em}}\Delta {SU\mathit{V}R}_{\mathit{non}-\mathit{lesioned}}\\ =\frac{{\mathit{SUVR}}_{\mathit{ipsilesional}\mathit{\hspace{0.17em}}\mathit{non}-\mathit{lesioned}}-S{\mathit{UVR}}_{\mathit{contralesional}\mathit{\hspace{0.17em}}\mathit{non}-\mathit{lesioned}}}{{\mathit{SUVR}}_{\mathit{contralesional}\mathit{\hspace{0.17em}}\mathit{non}-\mathit{lesioned}}}\end{array}$$

$$\begin{array}{l}\mathit{measured}\mathit{ }\Delta {\mathit{SUVR}}_{\mathit{cerebellum}}\\ =\frac{{\mathit{SUVR}}_{\mathit{contralesional}\mathit{ }\mathit{cerebellum}}-S{\mathit{UVR}}_{\mathit{ipsilesional}\mathit{ }\mathit{cerebellum}}}{{\mathit{SUVR}}_{\mathit{ipsilesional}\mathit{ }\mathit{cerebellum}}}\end{array}$$

Statistics

We performed general statistics in RStudio (v1.1.463. RStudio, Inc., Boston, MA) and we presented the data as mean ± standard deviation if normally distributed and as median (interquartile range (IQR)) if not normally distributed. We verified normality of the distributions with Shapiro-Wilk test (α = 0.05).

We tested differences in baseline characteristics with tests for proportions or Wilcoxon rank sum tests as appropriate.

To assess the changes in synaptic density within the lesion and the peri-ischemic area, we calculated the ratios for the lesion and peri-ischemic area for each patient compared to the set of healthy controls:

$${\mathit{ratio}}_{\mathit{lesion}/\mathit{controls}}=\frac{{\mathit{SUVR}}_{\mathit{lesion},\mathit{\hspace{0.17em}}\mathit{patient}}}{{\mathit{SUVR}}_{\mathit{lesion},\mathit{\hspace{0.17em}}\mathit{A}\mathit{V}\mathit{\hspace{0.17em}}\mathit{H}\mathit{C}}}$$

$${\mathit{ratio}}_{\mathit{peri}/\mathit{controls}}=\frac{{\mathit{SUVR}}_{\mathit{peri},\mathit{\hspace{0.17em}}\mathit{patient}}}{{\mathit{SUVR}}_{\mathit{peri},\mathit{\hspace{0.17em}}\mathit{A}\mathit{V}\mathit{\hspace{0.17em}}\mathit{H}\mathit{C}}}$$

We conducted one sample t-tests to assess whether the ratios were different from 1.

As the normal inter-subject variation of ¹¹C-UCB-J SUVR appeared rather high, we calculated the ratios for the lesions and peri-ischemic areas compared to the lesion/peri-ischemic mirror as well, as this may be more sensitive to detect a decrease:

$${\mathit{ratio}}_{\mathit{lesion}/\mathit{mirror}}=\frac{{\mathit{SUVR}}_{\mathit{lesion},\mathit{\hspace{0.17em}}\mathit{patient}}}{{\mathit{SUVR}}_{\mathit{lesion}\mathit{\hspace{0.17em}}\mathit{mirror},\mathit{patient}}}$$

$${\mathit{ratio}}_{\mathit{peri}/\mathit{mirror}}=\frac{{\mathit{SUVR}}_{\mathit{peri},\mathit{\hspace{0.17em}}\mathit{pat}\mathit{i}\mathit{ent}}}{{\mathit{SUVR}}_{\mathit{peri}\mathit{\hspace{0.17em}}\mathit{mirror},\mathit{\hspace{0.17em}patient}}}$$

We conducted one sample t-tests to assess whether the ratios were different from 1. We also conducted a one sample t-tests to assess whether the mean SUVR in the lesion area was different from 1.

We investigated the relationship between ratio_{lesion/controls} and ratio_{lesion/mirror} with the timing of scan (= days after stroke) with a Pearson correlation.

To detect differences in the non-lesioned tissue, we used Welch’s t-tests for group comparisons of VOI’s between patients and healthy controls. We conducted a paired t-test to detect differences between both hemispheres within the patients with ischemic stroke.

To investigate crossed cerebellar diaschisis on a group level, we used Welch’s t-test. We performed a paired t-test between ipsilesional and contralesional cerebellar hemispheres for the patients. We calculated individual asymmetry indices (AI’s) for all subjects to define crossed cerebellar diaschisis, as described in literature:^11–13

$$\mathit{A}\mathit{I}=\frac{{\mathit{SUVR}}_{\mathit{ipsilesional}\mathit{\hspace{0.17em}}\mathit{cerebellum}}-{\mathit{SUVR}}_{\mathit{contralesional}\mathit{\hspace{0.17em}}\mathit{cerebellum}}}{{\mathit{SUVR}}_{\mathit{ipsilesional}\mathit{\hspace{0.17em}}\mathit{cerebellum}}}\mathrm{*}100$$

AI’s of patients were considered significantly asymmetric if they exceeded the mean ± 2 standard deviation boundaries (= z-score >2 of <-2) of the AI’s of the controls. A significant AI in disfavor of the contralesional cerebellar hemisphere (= positive AI value) is the hallmark of crossed cerebellar diaschisis.

We also explored the effect of lesion size on ¹¹C-UCB-J SUVR by correlating lesion size with SUVR, with lesion ratios and peri-ischemic ratios and with cerebellar AI (Spearman correlation).

For all group level analyses with healthy controls, the right hemisphere was assigned as the equivalent of the affected hemisphere and the left hemisphere the equivalent of the unaffected hemisphere as the majority of the patients had a stroke lesion in the right hemisphere.

As this study was exploratory, significance level for all statistical analyses was set to α = 0.05 without correction for multiple comparisons. 95% confidence intervals (95% CI) are reported as well.

We performed all analyses without and with correction for PVE. All reported results are without correction for PVE unless otherwise specified. All PVE-corrected results are shown in the Supplementary Results (including Supplementary Table 1). In addition, we decided to analyze the data in the subgroup of patients without vascular lesions in the reference region (CS, n = 9) considering the possibility that variability in nonspecific uptake could increase as a result of white matter pathology. We reported these results only if they differed in reaching statistical significance from the main analysis. The full results of the subanalysis are provided in Supplementary Table 2.

Results

Demographics

Baseline characteristics are shown in Table 1 (and Supplementary Table 3). We included 20 patients with a first ischemic stroke of whom 19 underwent ¹¹C-UCB-J PET/MR imaging. We excluded 3 patients due to missing or insufficient quality of the MR data resulting in a dataset of 16 patients (5 male, 11 female). Median age was 74 years (IQR: 66–82 years) and patients were scanned at a mean interval of 21 ± 8 days after the event. Mean NIHSS at inclusion was 12 ± 5, median FMUE was 7.5 (IQR: 4–13), mean FMLE was 12 ± 8 and mean Barthel Index was 4 ± 3. Median stroke volume was 26 ml (IQR: 7.1–105 ml) which was present in the right hemisphere in 11 out of 16 patients.

Table 1.

Demographics and clinical characteristics.

	Patients with stroke	Healthy controls	p-value
Number of subjects (n)	16	40
Male / Female (n)	5/11	19/21	0.25
Age (years)	74 (66–82)	59.5 (43–72)	0.0085*
Injected 11C-UCB-J activity (MBq; median (IQR))	214 (204–231)	246 (216–294)	0.016*
Interval between event and scan (days; mean ± SD)	21 ± 8	NA
NIHSS (mean ± SD)	12 ± 5	NA
Fugl-Meyer – upper extremity (median (IQR))	7.5 (4–13)	NA
Fugl-Meyer – lower extremity (mean ± SD)	12 ± 8	NA
Barthel index (mean ± SD)	4 ± 3	NA
Lesioned hemisphere (n)	11 right/5 left	NA
Lesion volume (ml; median (IQR))	26 (7.1–105)	NA

Significant p-values are marked with *. For NIHSS, higher values represent more severe stroke (maximum score 42); for Barthel index, higher values represent more independence in activities of daily living (maximum score 20); for Fugl-Meyer, higher values represent better performance (maximum score upper extremity 66, maximum score lower extremity 34). IQR = interquartile range; NA = not applicable; NIHSS = National Institutes of Health Stroke Scale; SD = standard deviation.

The set of healthy controls consisted of 40 subjects (19 male, 21 female). These controls were younger compared to the patients (median age 59.5 years (IQR 43–72 years); 95% CI = 3 – 23, p = 0.0085).

¹¹C-UCB-J SUVR within the lesion and correlation with time of imaging

The mean ¹¹C-UCB-J SUVR within the lesion was 1.94 ± 0.71, which is significantly higher than 1 (95% CI = 1.6 – 2.3, p = 9.3*10⁻⁵). The ratios_{lesion/controls} of ¹¹C-UCB-J SUVR in the lesion were lower than 1 (mean = 0.67 ± 0.28, 95% CI = 0.52 – 0.82, p = 0.00023). The ratios_{lesion/mirror} also differed from 1 (mean = 0.66 ± 0.24, 95% CI = 0.53 – 0.78, p = 3.3*10⁻⁵; Table 2). There was no correlation between the decrease of SUVR (expressed as ratios_{lesion/controls} or ratios_{lesion/mirror}) and time between the index event and imaging (ratios_{lesion/controls}: r_P = −0.30, 95% CI = −0.69 – 0.23, p = 0.26; ratios_{lesion/mirror}: r_P = −0.35, 95% CI = −0.72 − 0.18, p = 0.19; Supplementary Figure 1).

Table 2.

Results of ¹¹C-UCB-J SUVR in the different regions for patients and healthy controls.

	Patients with stroke	Healthy controls	One sample t-test (µ=1)
Lesion	1.94 ± 0.71	NA	p = 9.3*10-5 *
ratiolesion/controls	0.67 ± 0.28		p = 0.00023 *
ratiolesion/mirror	0.66 ± 0.24		p = 3.3*10-5 *
Peri-ischemic area	2.66 ± 0.58	NA
ratioperi/controls	0.86 ± 0.18		p = 0.0076 *
ratioperi/mirror	0.85 ± 0.13		p = 0.00035 *
			Two sample t-test
Non-lesioned tissue
affected hemisphere	4.44 ± 0.81	4.61 ± 0.41	p = 0.43
unaffected hemisphere	4.61 ± 0.85	4.59 ± 0.41	p = 0.95
paired t-test	p = 0.0035 *
Cerebellar hemisphere
ipsilesional	3.98 ± 0.77	3.83 ± 0.37	p = 0.44
contralesional	3.84 ± 0.77	3.82 ± 0.38	p = 0.91
paired t-test	p = 0.0048 *

Data are presented as mean ± standard deviation. Significant p-values are marked with *. NA = not applicable; SUVR = Standardized Uptake Value Ratio.

¹¹C-UCB-J SUVR within the peri-ischemic area

The mean ¹¹C-UCB-J SUVR within the peri-ischemic area was 2.66 ± 0.58. The ratios_{peri/controls} of ¹¹C-UCB-J SUVR in the peri-ischemic area were lower than 1 (mean = 0.86 ± 0.18, 95% CI = 0.76 – 0.96, p = 0.0076) as were the ratios_peri/mirror (mean = 0.85 ± 0.13, 95% CI = 0.79 – 0.92, p = 0.00035; Table 2).

We performed partial volume correction (PVC) as this region is located at the border between the uninjured brain parenchyma and the ischemic lesion with decreased uptake. After PVC the ratios_{peri/controls} were no longer different from 1 (mean = 0.93 ± 0.19; 95% CI = 0.83 – 1.0; p = 0.16), but the ratios_peri/mirror of ¹¹C-UCB-J SUVR still differed from 1 (mean = 0.92 ± 0.12, 95% CI = 0.86 – 0.99; p = 0.026).

¹¹C-UCB-J SUVR in the non-lesioned tissue

¹¹C-UCB-J SUVR in the non-lesioned grey matter part of the affected hemisphere (mean = 4.44 ± 0.81) did not differ from ¹¹C-UCB-J SUVR of healthy controls in the right hemisphere (mean = 4.61 ± 0.41, 95% CI = −0.61 – 0.27, p = 0.43). ¹¹C-UCB-J SUVR of the grey matter of the unaffected hemisphere (mean = 4.61 ± 0.85) did not differ from ¹¹C-UCB-J SUVR of healthy controls in the left hemisphere (mean = 4.59 ± 0.41, 95% CI = −0.45 – 0.48, p = 0.95). In an intra-individual analysis, we did identify decreased ¹¹C-UCB-J SUVR when comparing the non-lesioned tissue grey matter of the affected hemisphere with the unaffected hemisphere (mean difference = −0.17, 95% CI = −0.27 – −0.064, p = 0.0035; Table 2).

In the subgroup analysis of patients without vascular lesions in the CS (n = 9), we also detected a decrease of ¹¹C-UCB-J SUVR in the non-lesioned grey matter of the affected hemisphere compared to the hemispheric grey matter of healthy controls (mean patients = 4.04 ± 0.63, 95% CI = −1.1 – −0.073, p = 0.029).

¹¹C-UCB-J SUVR in the cerebellar hemispheres

On a group level, ¹¹C-UCB-J SUVR did not differ in the ipsilesional (mean = 3.98 ± 0.77), nor contralesional cerebellar hemisphere (mean = 3.84 ± 0.77) of the patients with stroke compared to healthy controls (right: mean = 3.83 ± 0.37; left: mean = 3.82 ± 0.38; ipsilesional vs right: 95% CI = −0.27 – 0.58, p = 0.44; contralesional vs left: 95% CI = −0.40 – 0.45, p = 0.91). In stroke patients we assessed the presence of an asymmetry by comparing the ipsilesional cerebellar hemisphere with the contralesional cerebellar hemisphere. The contralesional cerebellar hemisphere showed a decrease of ¹¹C-UCB-J SUVR compared to the ipsilesional cerebellar hemisphere (mean difference = −0.14, 95% CI = −0.24 – −0.051, p = 0.0048; Table 2).

The mean AI in the cerebellar hemispheres of healthy controls was 0.18 ± 1.6. The mean AI of patients was 3.6 ± 4.9 and in 9 patients (56%) the AI exceeded the 2 standard deviation boundaries of the AI in healthy controls, indicative of cerebellar diaschisis. In most of those patients (8; 89%) the asymmetry was in disfavor of the contralesional cerebellum.

Effect of lesion size on ¹¹C-UCB-J SUVR

There was no correlation between lesion size and SUVR in any of the investigated regions (lesion, peri-ischemic area, non-lesioned tissue and cerebellar hemispheres). There was no correlation between lesion size and lesion ratios or peri-ischemic ratios, nor between lesion size and cerebellar AI.

Simulation of perfusion effects

The results of the simulation study are shown in Supplementary Figure 2 and Supplementary Table 4. The measured ΔSUVR_non-lesioned ranged from -13% to +5.5% (mean = −3.6% ± 4.1%), whereas the simulated ΔSUVR was −4.3% for ΔK₁ −15% and +2.5% for ΔK₁ +15% for cerebral cortical regions. The measured ΔSUVR_cerebellum ranged from −16% to +6.6% (mean = −4.0% ± 5.4%), whereas the simulated ΔSUVR was −7.2% for ΔK₁ −35% and −0.27% for ΔK₁ +15% for the cerebellum.

Discussion

In this study, we identified a decrease, but no total loss, of synaptic density with ¹¹C-UCB-J PET within the ischemic stroke lesion several weeks after the index event. Evaluating synaptic density in the non-ischemic brain parenchyma revealed a reduction in the ipsilesional hemisphere compared to the contralesional hemisphere and we detected crossed cerebellar diaschisis.

The strong decrease of synaptic density signal within the FLAIR-based ischemic lesion is a largely expected finding as ischemia results in cell death and tissue loss. However, we did not observe a total loss of synaptic density in the first weeks after the ischemic stroke, as reflected by values above 1 in most lesions. Potentially, synaptic loss could progress over weeks which we studied by analyzing the correlation between time after onset and synaptic density. Likely due to the low sample size, rather narrow range of time after stroke onset and imaging between patients and heterogeneous population we could not detect this correlation.

In the peri-ischemic area we could not demonstrate any increases in synaptic density in the tissue surrounding the lesion. This is in contrast to preclinical research that described increases in synaptic density and metabolism in the peri-ischemic tissue.^6,41 Heterogeneity within the peri-ischemic area, the early time window between the index event and imaging as well as insufficient resolution of PET to detect subtle changes, may account for the lack of identifying this in our cohort. However, dissimilarities between different species or non-representative preclinical stroke models may also contribute to the contrasting results. ⁴² Our results are in line with a ¹⁸ F-flumazenil study, showing a decreased signal in the peri-infarct cortex compared to contralesional cortex. ^{43 18}F-flumazenil binds to postsynaptic GABA_A (gamma-aminobutyric acid) receptors and is therefore considered as a proxy for neuronal density.

In the non-lesioned tissue of the affected hemisphere, we found a lower synaptic density compared to the unaffected hemisphere. Stroke lesions and underlying white matter disease were quite heterogeneous as reflected in the wide range in SUVR, which may have hampered detection of changes in the total cohort of patients compared to controls. In the subgroup of patients without white matter pathology in the CS, we did observe a reduction in the ipsilesional hemisphere compared to controls. This may suggest reduced synaptic density to exist in the ipsilesional hemisphere outside the ischemic lesion. This reduction may be explained by Wallerian degeneration, a phenomenon in which distal axons of damaged neurons degenerate and therefore cause a decrease of presynaptic vesicles. Although most studies investigating Wallerian degeneration after stroke focused on the corticospinal tract,^44,45 reduced intrahemispheric connections are also described.^46,47 Likewise, transcallosal Wallerian degeneration leading to reduced interhemispheric connections has been reported as well,^48,49 but this was not (yet) seen in our population. Follow-up on outcome of this cohort and further research in larger samples will have to explore if this reduced synaptic density within and/or outside the ischemic lesion is linked to outcome and may have a predictive value.

Neural plasticity, including reorganization and strengthening of existing connections, could translate into functional recovery. Underlying processes such as axonal sprouting and synaptic alterations may be at play. Functional compensation after ischemic stroke has been explored in patients by measuring relative cerebral blood flow using ¹⁵O tracers and resulted in a variety of patterns of ipsilesional, contralesional or bilateral (over)activation.^50–56 The populations studied previously were always relatively small and heterogeneous, both in terms of included patients as well as timing of scans, possibly responsible for rather ambiguous results. We explored changes in synaptic density by measuring SV2A, which is a marker of structural presynaptic vesicles, as a potential biomarker of neuronal loss or neural plasticity. So far, ¹¹C-UCB-J PET has only been used to detect losses of synaptic density^57–64 and we have no data on the relationship between changes in SV2A density and synaptic function, making it uncertain if upregulation of SV2A is related to plasticity or recovery. It is unclear which time window is most optimal to identify loss of synapses, axonal reorganization or strengthening of existing connections. Most studies exploring functional recovery after stroke so far, focused on later time windows. Considering the short time interval between the ischemic insult and PET imaging in the presented data, we assumed to be able to capture synaptic loss but possibly not yet increases in the SUVR, as potential indicator of synaptic plasticity. In the first month after ischemic stroke, we did not detect any increases in synaptic density. It is possible that the neural plasticity has only been initiated and will further develop over weeks and months which may become more apparent during the follow-up of our cohort with imaging at 6 months and 2 years after the index event.

Crossed cerebellar diaschisis is a known phenomenon in the first weeks after ischemic stroke and reflects a decrease in function, metabolism or perfusion of the contralesional cerebellum caused by a supratentorial, remotely connected lesion. Perfusion PET/SPECT (single photon emission computed tomography),^{10–13,40,65} FDG (fluorodeoxyglucose) PET¹⁴ and ASL (arterial spin labeled) MRI ¹⁵ have revealed these changes. We demonstrated that crossed cerebellar diaschisis can be assessed with ¹¹C-UCB-J PET and the rate as measured with ¹¹C-UCB-J in our population (50%) was in line with previously reported percentages as measured with perfusion SPECT/PET in patients with supratentorial lesions (46-62.5%).^11,40,65 One patient had a significant cerebellar asymmetry reversed to the expected crossed cerebellar diaschisis, which is previously described as well. ¹¹ Thus far, the exact interpretation of the crossed cerebellar diaschisis that we found is unclear. Although it can reflect a loss of cortical input connections, it is not excluded that the results are (partially) driven by changes in perfusion as quantification with SUVR may be influenced by perfusion effects. However, the results of the simulation study suggests that perfusion changes are likely not sufficient to fully explain our findings as some individual differences in cerebellar asymmetry were more pronounced than could be explained by perfusion changes alone. A ¹²³I-iomazenil SPECT study could not identify crossed cerebellar diaschisis in two patients scanned at an earlier time point. ⁶⁶ This might suggest that the presynaptic SV2A tracer and postsynaptic GABA_A receptor tracer behave differently, but differences in time window after the event and the fact that not all stroke patients manifest crossed cerebellar diaschisis, may be other explanations.

Our study has some limitations. The first relates to the quantification method which made use of a reference tissue normalization with 60-90 minutes SUVR data normalized on the CS as proxy to the distribution volume. Although this reference region is validated by different research groups^19,34,67,68 and used in several patient populations,^{57–59,61–64} validation in a cohort of ischemic stroke patients is lacking. The main concern is the potential influence of white matter disease in the CS as it is unsure whether the nonspecific binding in the CS is influenced by white matter pathology. There was no difference in SUV of the CS between patients and healthy controls in our study (See Supplementary Material), but since the role of white matter changes on the SUV in the CS is unknown, we conducted a subgroup analysis in patients without white matter pathology in the CS (n = 9). This approach reduced the variation in SUVR compared to the total cohort which is suggestive of a potential influence of white matter disease in the CS. This needs further follow-up and study in a larger cohort of patients with white matter pathology to ensure the validity of the CS as a reference region or evaluate a different reference region. However, the reproducibility of the results in the within subjects analyses, which are not influenced by the reference region, suggests that our findings are valid, even despite the potential influence of white matter pathology. Moreover, the results of a post hoc sensitivity analysis with pons as normalization region were similar to the results of the subanalysis with CS as reference region (data not shown). The pons is no validated reference region for ¹¹C-UCB-J PET quantification, but the similarities between all different methods strengthen the validity of our results. Another limitation of the quantification method is the possible bias in SUVR that is introduced by perfusion changes and/or changes in blood brain barrier permeability due to blood brain barrier breakdown in the subacute phase after ischemic stroke. Based on our simulation study we cannot exclude that perfusion changes after stroke may have influenced the SUVR. This could have been further explored by dynamic scan acquisition with full kinetic modelling which however was considered unfeasible in this patient population. Another limitation is the difference in age between the two groups as the patients were older compared to the controls. We have previously demonstrated that there is no major influence of age on synaptic density with SUVR as measure and therefore we included the full set of healthy controls instead of limiting the control cohort to age-matched subjects as this would decrease the power of the study. ⁶⁹ To confirm this hypothesis we conducted a post hoc analysis with an age-matched control cohort (data not shown) which revealed comparable results. Another difference we observed between patients and controls was the lower injected ¹¹C-UCB-J activity which was inevitable as patients will receive follow-up scans, necessitating a lower injected dose to not exceed the legally allowed yearly radiation dose (as this study is part of a multimodal, multitracer imaging study). Last, ¹¹C-UCB-J PET is a marker of presynaptic SV2A vesicles and the approximation of SV2A PET as proxy for synaptic density is under the assumption that there is a correlation between density of SV2A vesicles and synaptic density and that this correlation is not influenced by the pathology under investigation. Moreover, although there is no evidence of off-target binding to inflammatory cells so far, this possibility should be kept in mind as subacute stroke is accompanied by an inflammatory process. Future in vitro assays should assess the possibility of off-target binding to inflammatory or other cells.

A strength of our study is the inclusion of healthy controls. Neuroimaging studies of patients with stroke frequently define the contralateral hemisphere as control region instead of hemispheres of healthy controls. This hampers the detection of potential changes occurring in the unaffected hemisphere. This may become of even greater relevance in the longitudinal follow-up of patients which will enable making comparisons with the contralateral hemisphere within each subject, but also with healthy controls. With this follow-up we intend to provide more insights in changes in synaptic density during the recovery phase after an ischemic insult.

In conclusion, ¹¹C-UCB-J PET was able to demonstrate a decrease, but not total loss, of synaptic density in ischemic stroke lesions in the first month after the index event. When comparing non-lesioned tissue in the two cerebral hemispheres, an asymmetry in synaptic density was present likely caused by a decrease in the ipsilesional hemisphere. ¹¹C-UCB-J PET was able to detect crossed cerebellar diaschisis in a proportion of patients with ischemic stroke. Further follow-up of patients with ¹¹C-UCB-J PET is needed to gain more insight in the temporal and spatial changes of synaptic density after ischemic stroke and their relation to outcome parameters.

Supplemental Material

sj-pdf-1-jcb-10.1177_0271678X211047759 - Supplemental material for Changes in synaptic density in the subacute phase after ischemic stroke: A ¹¹C-UCB-J PET/MR study

Supplemental material, sj-pdf-1-jcb-10.1177_0271678X211047759 for Changes in synaptic density in the subacute phase after ischemic stroke: A ¹¹C-UCB-J PET/MR study by Laura Michiels, Nathalie Mertens, Liselot Thijs, Ahmed Radwan, Stefan Sunaert, Mathieu Vandenbulcke, Geert Verheyden, Michel Koole, Koen Van Laere and Robin Lemmens in Journal of Cerebral Blood Flow & Metabolism