Current accepted cerebrovascular reactivity indices suffer from the need of high frequency data capture and export for post-acquisition processing. The role for minute-by-minute data in cerebrovascular reactivity monitoring remains uncertain. The goal was to explore the statistical time-series relationships between intra-cranial pressure (ICP), mean arterial pressure (MAP) and pressure reactivity index (PRx) using both 10-s and minute data update frequency in TBI. Prospective data from 31 patients from 3 centers with moderate/severe TBI and high-frequency archived physiology were reviewed. Both 10-s by 10-s and minute-by-minute mean values were derived for ICP and MAP for each patient. Similarly, PRx was derived using 30 consecutive 10-s data points, updated every minute. While long-PRx (L-PRx) was derived via similar methodology using minute-by-minute data, with L-PRx derived using various window lengths (5, 10, 20, 30, 40, and 60 min; denoted L-PRx_5, etc.). Time-series autoregressive integrative moving average (ARIMA) and vector autoregressive integrative moving average (VARIMA) models were created to analyze the relationship of these parameters over time. ARIMA modelling, Granger causality testing and VARIMA impulse response function (IRF) plotting demonstrated that similar information is carried in minute mean ICP and MAP data, compared to 10-s mean slow-wave ICP and MAP data. Shorter window L-PRx variants, such as L-PRx_5, appear to have a similar ARIMA structure, have a linear association with PRx and display moderate-to-strong correlations (r ~ 0.700, p < 0.0001 for each patient). Thus, these particular L-PRx variants appear closest in nature to standard PRx. ICP and MAP derived via 10-s or minute based averaging display similar statistical time-series structure and co-variance patterns. PRx and L-PRx based on shorter windows also behave similarly over time. These results imply certain L-PRx variants may carry similar information to PRx in TBI.

Background Chordomas are rare malignant bone cancers of the skull-base and spine. Patient survival is variable and not reliably predicted using clinical factors or molecular features. This study identifies prognostic epigenetic chordoma subtypes that are detected non-invasively using plasma methylomes. Methods Methylation profiles of 68 chordoma surgical samples were obtained between 1996-2018 across three international centres along with matched plasma methylomes where available. Results Consensus clustering identified two stable tissue clusters with a disease-specific survival difference that was independent of clinical factors in a multivariate Cox analysis (HR=14.2, 95%CI: 2.1–94.8, p=0.0063). Immune-related pathways with genes hypomethylated at promoters and increased immune cell abundance were observed in the poor-performing “Immune-infiltrated” subtype. Cell-to-cell interaction plus extracellular matrix pathway hypomethylation and higher tumor purity was observed in the better-performing “Cellular” subtype. The findings were validated in additional DNA methylation and RNA sequencing datasets as well as with immunohistochemical staining. Plasma methylomes distinguished chordomas from other clinical differential diagnoses by applying fifty chordoma-versus-other binomial generalized linear models in random 20% testing sets (mean AUROC=0.84, 95%CI: 0.52-1.00). Tissue-based and plasma-based methylation signals were highly correlated in both prognostic clusters. Additionally, leave-one-out models accurately classified all tumors into their correct cluster based on plasma methylation data. Conclusions Here, we show the first identification of prognostic epigenetic chordoma subtypes and first use of plasma methylome-based biomarkers to non-invasively diagnose and subtype chordomas. These results may transform patient management by allowing treatment aggressiveness to be balanced with patient risk according to prognosis.