Research Article Volume 10 Issue 6
Hackensack Meridian JFK University Medical Center, USA
Correspondence: Dr. Farah Fourcand, Stroke and Neurovascular Center, 65 James Street, Edison, New Jersey
Received: October 28, 2020 | Published: November 24, 2020
Citation: Fourcand F, Kulhair A, Mehta S, et al. Virtual reality adds value to analysis of aneurysms post-flow diversion, coiling, and clipping. J Neurol Stroke. 2020;10(6):195-201. DOI: 10.15406/jnsk.2020.10.00438
Purpose: The purpose of this study was to measure sensitivity of virtual reality (VR) in detecting biomarkers of neurovascular remodeling suboptimally evaluated in digital subtraction angiography (DSA) of treated unruptured, intracranial aneurysms.
Methods: The sensitivity of virtual reality and digital subtraction angiography in detection of neurovascular biomarkers in aneurysms treated with flow diversion, coiling, and clipping were evaluated. Validated grading scales were integrated into a standardized rating platform. The respective novel and conceptual measures of minimal imaging important difference (MIID) and number needed to image (NNI) were calculated for each biomarker.
Results: In flow diversion, coiling, and clipping, minimal imaging important difference and number needed to image were associated with virtual reality in detection of abnormal biomarkers, with the exception of stasis phase associated with digital subtraction angiography. Number needed to image was associated with flow diversion stent stenosis (RR: 7.00, 95% CI 0.37 to 131.97; OR: 7.46, 95% CI 0.38 to 148.49). Minimal imaging important difference was greatest in residual aneurysm filling (25%±66, 95% CI) in flow diversion and Meyer score in coiling (42.5%±17.69, 95% CI) and clipping (22.2%±13.58, 95% CI). Regression models demonstrated minimal imaging important difference and number needed to image shared a significant correlation (R20.99, 95% CI, p<0.001).
Conclusion: Virtual reality adds value to digital subtraction angiography in evaluation of aneurysms treated with flow diversion, coiling, and clipping. Larger, prospective studies are warranted to increase statistical power and validate clinical significance.
Keywords: sphenoid bone, craniofacial abnormalities, malocclusion, megalencephaly, facial asymmetry, exophthalmos
Unruptured, intracranial aneurysms, have a prevalence of 3% worldwide, global incidence of 700,000 person-years, and rupture rate of 0.3% per year.1 As neuroimaging advances, more asymptomatic aneurysms are detected incidentally, which does not correlate with increased risk of rupture.1 Why then, is management of aneurysms a highly involved process warranting a multidisciplinary approach?
The reason is that while aneurysmal SAH has a mortality of 35%, treatment carries risk of poor neurological outcome in <1% to 10% of cases depending on treatment type and device generation, which may be higher than risk of rupture.2,3 In addition to patient-specific modifiable and non-modifiable risk factors, the International Study of Unruptured Intracranial Aneurysms,4 American Heart Association/American Stroke Association guidelines,5 validated grading scales.5,6 clinical trials,7 and meta-analyses6 have also focused on aneurysm-specific neuroimaging biometrics, such as location, size, and aspects ratios (dome: neck, height:neck). These biometrics are macroscopic results of microscopic neurovascular changes, or biomarkers, mainly studied in basic science research, but important in clinical decision-making.8–11 And so, we looked at what translational tools exist to fill this unmet need in the clinical research setting using virtual reality (VR).12 Currently, some neurovascular biomarkers in flow diverted, coiled, and clipped aneurysms are suboptimally evaluated with standard imaging and diagnostic tools making their validated clinical significance unknown.13–19 The use of VR in treated aneurysms has not been explored.
DSA, the technology used to extract images from diagnostic cerebral angiogram, is the gold standard for surveillance of treated aneurysms. The purpose is to evaluate evolving characteristics of aneurysms at prospective intervals over time.8,9 The clinician is then more adept at predicting clinical outcomes that guide management and understanding an aneurysm’s natural history, the cumulative aim of which is to develop novel diagnostic tools and biomarkers of prognostication.10,11,20–26
Diagnostic limitations
In flow diversion, limitations of DSA are suboptimal visualization of stent structure and alignment, free-hanging edges, and stent stenosis.27 In coiled aneurysms, complications may arise from aberrancies in coil-thrombus complex.25,28 Predictive biomarkers include thrombus organization29 coil packing30 and coil orientation.31 In clipping of extradural, petrous or cavernous segment aneurysms, DSA is unable to assess bony erosion or peripheral nerve pathology.32,33 In CTA, which is part of the surveillance algorithm, clip and bone artifact are limiting factors.27,33
All data and materials have been made publicly available in Open Science Framework and can be accessed at https://osf.io/6kwyc/. VR reconstructions of unruptured, intracranial aneurysms treated with flow diversion,34 coiling, and clipping from March 2010 to September 2018 at our university-affiliated, community-based comprehensive stroke center were retrospectively built using Surgical Theater35 (Precision VR® by Surgical Theater Inc., Cleveland, OH). Four institutional endovascular surgical neuroradiologists independently evaluated cases. Variables were inputted into a standardized rating platform to facilitate comparison of VR and DSA sensitivity in detecting biomarkers of neurovascular remodeling. Validated grading scales were integrated to increase intra- and interrater agreement and generalizability.
Validated grading of neurovascular biomarkers
Flow diversion cases were evaluated through adaptation of variables from the following validated grading scales: O’Kelly-Marotta,20 SMART,36 Kamran-Byrne,37 and flow-diverting stent scale.38 In our study, parameters of residual aneurysm filling were: none 0%, entry remnant <5%, subtotal 5-95%, and total >95%. Parameters of stasis phase were: arterial clearance before capillary phase, clearance before venous phase, or stagnation in venous phase. Parameters of branch filling defined by TICI39 were: complete (TICI 3), partial I (TICI 2b >50%), partial II (TICI 2a <50%), or insufficient. Endoleak was identified as: absent or present. Free-hanging edge was identified as: absent, present without thrombus, or present with thrombus. Parameters of stent stenosis adapted from NASCET criteria40 were: minimal 0-29%, mild 30-49%, moderate 50-69%, severe 70-94%, critical 95-99%, or occluded 100%. Coiling cases were evaluated using variables adapted from modified Raymond Roy classification41 (class I: complete obliteration, class II: residual neck, class IIIa: residual aneurysm with contrast opacification within coil, class IIIb: contrast opacification outside coil) and Meyer scale42 (grade 1 > or = 90%, grade 2 70–89%, grade 3 50–69%, grade 4 26–49%, or grade 5 < or = 25% occlusion; modifier I: coil opacification; modifier G: independent aneurysm growth). Clipping cases were evaluated using variables adapted from Raymond Roy occlusion classification41 (class I: complete obliteration, class II: residual neck, class III: residual aneurysm) and Meyer scale. Qualitative information was included, if applicable.
360° VR model rendering and manipulation
The VR model can be manipulated in any axis or magnitude to examine global neuroanatomy and neuropathology. VR software controls scan intensities to visualize selected systems and facilitate treatment planning for neurovascular conditions.12 This technology is also able to highlight specific components in flow diverters, coils, and clips as well as clip-in through any aspect of aneurysm or device.
Each aneurysm was reconstructed as a 360°VR model using patient-specific CT volumetric scans (3DCT) derived from DSA. Pre-contrast images, or mask,43,44 and contrast-enhanced images that subtract the mask, or sub, were the two sequences used. These sequences are established protocols in our institution’s neuroimaging department. There are specific resolution parameters to optimize scans for VR reconstruction. The dimensions are 512, 512, 502 on an x, y, z plane with 0.25mm spacing per slice collecting 498 to 502 files per sequence. Sub and mask 3DCTs were loaded into VR software to render the 360°VR model, which was later edited to highlight critical components. For each case, sub and mask 3DCTs from the most recent DSA were merged.
The sub sequence only showed vessels from DSA. Specific digital presets were created to show parent vessels and distal branches around aneurysms. The mask sequence showed bone and other radiopaque elements, such as metal in flow diverters (cobalt-chromium and platinum-tungsten wires), coils (platinum), and clips (titanium or Phynox® alloy45). A mesh generation tool produced a 3D object of each device using high-intensity variations per voxel from mask sequence. This extraction allowed raters to interact with devices dependently and independently of sub sequence.
Data analysis
Relative risk (RR)46,47 (95% CI) and OR46, 47 (95% CI) for each neurovascular biomarker was determined. Haldane-Anscombe correction was used for zero values in contingency tables,48 if applicable. Pearson’s coefficient was used to measure association of RR and OR size with CI width.49
The novel concept of minimal imaging important difference was used to define the smallest amount a neuroimaging finding changed to be meaningful. Sensitivity of VR in detecting a neuroimaging abnormality to a greater degree than DSA was defined as a meaningful difference. The diagnostic relevance of minimal imaging important difference was based on principles of clinical relevance in minimal clinically important difference.50,51 The conceptual measure of number needed to image52 (NNI) was defined as number of aneurysms with VR needed to detect an abnormality to a greater degree than DSA. The diagnostic relevance of NNI was based on principles of clinical relevance in number needed to treat.52,53 Minimal imaging important difference was expressed as a percentage with associated margin of error (95% CI).
In NNI, to account for potentially insignificant CI given dependence on significant absolute risk reduction, CI were expressed as: NNI for an additional harmful outcome to ∞ to NNI for an additional beneficial outcome. In this method, CI ranges from NNI (harm)=-1 to NNI (benefit)=1 via infinity, whereby these are theoretically unattainable values (absolute risk reduction of -100% and 100%, respectively) and NNI=∞ represents midpoint (absolute risk reduction=0%).54
All neurovascular biomarkers were analyzed weighted and unweighted. For example, in flow diversion cases, normal designation (i.e. TICI 3) was assigned weight of 1, TICI 2b weight of 2, TICI 2a weight of 3, and insufficient weight of 4. Cases were measured unweighted in a dichotomous manner (TICI 3 versus non-TICI 3) and weighted (TICI 3 versus summation of non-TICI 3) for minimal imaging important difference and NNI. Two-tailed t-test was performed to evaluate for significant differences between unweighted and weighted samples.46,47,49
Correlation scatter models were made using ordinary least-squares regression55 (quadratic regression equation: y = a + bx +cx2, 95% confidence and probability) where R, or coefficient of determination, determines strength of association between x (minimal imaging important difference) and y (NNI) and variability in y (NNI) explained by model in weighted and unweighted analysis. Outlier detection was performed prior to regression analysis to account for heterogeneity and limited sample size. Statistical power was calculated using dichotomous endpoint, two-independent sample study model49. Statistical software used was MedCalc Software Ltd.56
From March 2010 to September 2018, 86 flow diversion, 62 coiling, and 96 clipping cases for treatment of unruptured, intracranial aneurysms were performed. Of these cases, 48 flow diversion, 15 coiling, and 18 clipping cases met technical parameters required to build a VR model.
Measure of association between VR and neuroimaging abnormality, or OR (Table 1), and probability of VR detecting neuroimaging abnormality, or RR (Table 2), were greater than 1 in all neurovascular biomarkers (n=10) and greatest in stent stenosis (RR: 7.00, 95% CI 0.40 to 132.00; OR: 7.50, 95% CI 0.40 to 149.00). Wide CI shared a positive correlation with high RR (Pearson correlation coefficient 0.98, p<0.0001, 95% CI 0.90 to 0.99) and high OR (Pearson correlation coefficient 0.97, p<0.0001, 95% CI 0.89 to 0.99).
Procedure |
Biomarkers |
OR 95%CI |
Minimal imaging important difference1 |
|
Unweighted 95%CI |
Weighted 95%CI |
|||
Flow divert |
Residual |
1.82, 0.80 to 4.13 |
14.6±7.06 |
25.0±8.66 |
|
Stasis Phase2 |
0.04, 0.002 to 0.67 |
-20.8±8.12 |
-52.1±9.99 |
|
Endoleak |
2.04, 0.18 to 23.32 |
2.1±2.87 |
4.1±3.97 |
|
Hanging Edge |
2.04, 0.18 to 23.32 |
2.1±2.87 |
4.1±3.97 |
|
Stent Stenosis |
7.46, 0.38 to 148.5 |
6.25±4.84 |
12.5±6.62 |
|
Branch Filling |
1.72, 0.52 to 5.70 |
6.3±4.86 |
17.1±7.53 |
Coil |
Raymond Roy3 |
1.71, 0.40 to 7.29 |
13.3±12.15 |
33.3±16.87 |
|
Meyer Scale |
1.75, 0.39 to 7.86 |
13.8±12.34 |
42.5±17.69 |
Clip |
Raymond Roy |
1.35, 0.30 to 6.13 |
5.6±7.51 |
11.2±10.30 |
|
Meyer Scale |
1.35, 0.30 to 6.13 |
5.6±7.51 |
22.2±13.58 |
Table 1 Association of virtual reality and neurovascular biomarkers
Footrnotes:
1Expressed as percentage with margin of error
2OR <1 and negative minimal imaging important difference represent greater association with DSA
3Modified Raymond Roy classification
Procedure |
Biomarkers |
RR 95%CI |
Number needed to image1 |
|
Unweighted 95%CI |
Weighted 95%CI |
|||
Flow Divert |
Residual |
1.41, 0.88 to 2.27 |
7, 2.93 ∞ 20.04 |
10, 4.08 ∞ 26.24 |
|
Stasis Phase2 |
0.05, 0.003 to 0.79 |
-5, 3.10 to 11.66 |
-3, 1.90 to 4.02 |
|
Endoleak |
2.00, 0.19 to 21.33 |
48, 11.07 ∞ 20.55 |
25, 7.53 ∞ 18.35 |
|
Hanging Edge |
2.00, 0.19 to 21.33 |
48, 11.07 ∞ 20.55 |
25, 7.53 ∞ 18.35 |
|
Stent Stenosis |
7.00, 0.37 to 131.97 |
16, 7.21 ∞ 61.80 |
9, 4.73 to 54.80 |
|
Branch Filling |
1.60, 0.56 to 4.54 |
16, 5.03 ∞ 13.55 |
11, 3.87 ∞ 14.28 |
Coil |
Raymond Roy3 |
1.33, 0.61 to 2.91 |
8, 2.05 ∞ 4.56 |
9, 2.71 ∞ 7.28 |
|
Meyer Scale |
1.36, 0.61 to 2.99 |
7, 1.98 ∞ 4.38 |
9, 2.96 ∞ 8.12 |
Clip |
Raymond Roy |
1.25, 0.40 to 3.91 |
18, 2.96 ∞ 4.41 |
14, 2.80 ∞ 4.66 |
|
Meyer Scale |
1.25, 0.40 to 3.91 |
18, 2.96 ∞ 4.41 |
10, 2.67 ∞ 5.84 |
Table 2 Probability of virtual reality detecting neurovascular biomarkers
Footnotes:
1Rounded to nearest whole number
2RR <1 and negative NNI represent higher probability in DSA
3Modified Raymond Roy classification
Figures 1 and 2 provided visual representation of aberrant neurovascular biomarkers in VR, which were quantified by minimal imaging important difference (Table 1) and NNI (Table 2). In NNI, CI included infinity with the exception of weighted stent stenosis in VR that had significant absolute risk reduction (11.76%, 95% CI 2.35-21.18). There were significant differences between unweighted and weighted minimal imaging important difference (test statistic t=3.70, P=0.006) and NNI (test statistic t=-2.29, P=0.05).
Figure 1 Biomarker detection in flow diversion. VR model of treated paraophthalmic aneurysm with endoleak (A) is not visualized on DSA (B). VR shows ophthalmic artery branch filling (C) not visualized on DSA (D).
Figure 2 Device abnormalities detected in VR. Irregularly-shaped coil of treated posterior communicating artery aneurysm (A) protrudes into parent vessel increasing coil herniation risk. Stent-stenosis of flow diverted paraophthalmic aneurysm (B) is shown. Residual neck of clipped PICA aneurysm in VR (C) is difficult to visualize on DSA (D).
Scatterplot regression model (Figure 3) demonstrated NNI (n=9) and minimal imaging important difference (n=9) had significant correlation (unweighted: R2 0.99, 95% CI, p<0.001) and variability in y (NNI) was explained by model. In weighted analysis (Figure 4), minimal imaging important difference was greater and NNI lower in the same biomarkers (weighted: R2 0.92, 95% CI, p<0.001). Outlier detection did not show significant outliers (NNI: coefficient of skewness 1.27, P=0.08; coefficient of Kurtosis 0.20, P=0.74) (minimal imaging important difference: coefficient of skewness 0.46, P=0.51; coefficient of Kurtosis -1.50, P=0.21). Subgroup analysis of flow diversion, coiling, and clipping was not performed to prevent bias given small sample size and study design limitations. Of note, stasis phase was the only biomarker more strongly associated with DSA (OR 0.04, 95% CI 0.002 to 0.70; RR 0.05, 95% CI 0.002 to 0.80). Negative minimal imaging important difference and NNI for an additional harmful outcome were reported. Values were included in all analyses except regression model.
To our knowledge, this is the first study comparing VR to DSA in evaluation of unruptured, intracranial aneurysms following flow diversion, coiling, and clipping. Results suggested a conservative degree of internal validity in VR conferring an advantage in detection of neuroimaging abnormalities representative of aberrant neurovascular remodeling. These limitations are the result of small sample size with suboptimal precision and generalizability. A modest estimation of minimal sample size (alpha type I error: 0.05, beta type II error: power 80%) assuming small effect size (5% increased incidence) of VR in detecting a neuroimaging abnormality was 304. Prospective, multicenter studies and concurrent advances in existing technologies, standardization protocols, and resource allocation to reliably meet parameters to build a VR model from existing cases are needed.
Clinical relevance of statistical significance
Our study proposed a novel concept: minimal imaging important difference based on minimal clinically important difference, which prioritizes objective measures of beneficence, rather than statistical distribution-based approaches57 The conceptual measure of NNI52 was evaluated based on number needed to treat, a type of effect size58 vital to interpretation of results associated with, yet independent of, statistical significance.
There were challenges in determining clinical significance in neuroimaging-dependent measures of minimal imaging important difference and NNI as compared to minimal clinically important difference and number needed to treat. As such, we attempted to test if a statistically significant association existed between minimal imaging important difference and NNI that would strengthen their independent reliability and synergistically validate their integration into diagnostic and predictive tools with external validity. Correlation and regression models suggested a significant association existed between minimal imaging important difference and NNI (unweighted: R2 0.99, weighted: R2 0.92), whereby values nearing R=1 or 100% of variance theoretically presume fitted values always equal observed values and all data points fall on a fitted regression line. However, residual plots of regression models showed marginal bias in unweighted analysis and marked bias in weighted analysis. We posit that further validation and rigorous critical appraisal of minimal imaging important difference and NNI as interdependent variables is warranted to potentially reveal a symbiotic relationship that achieves statistical significance. The redistribution of weight assignments in weighted analysis is also needed.
It follows that although neuroimaging biomarkers, both aberrant and adaptive, reflect a discrete stage of the neurovascular remodeling process, their dynamic interplay is unknown. Validated grading scales in this study were adapted as static measures of clinical significance for abnormalities that reflect these dynamic processes. In comparison to coiling and clipping, flow diversion is a more recent technology undergoing rapid evolution of device design, indications for use, and, for our purposes, grading scale development.20,36–38 As such, grading scales that integrate principles of hemodynamics and endothelialization in aneurysmal thrombosis into standard biometrics will continue to be redefined to accurately and reliably reflect clinical significance.59
Future directions of virtual reality
Although biomarkers suboptimally evaluated on DSA were appreciable in VR, the technology was unable to identify other fundamental biomarkers like organization of thrombus-coil complex. Presently, neither neuroscience-based VR nor standard tools are designed to evaluate histological and molecular components that are, in effect, the primordial soup of the remodeling process.60–62 Pathology-based augmented-reality models63 in the pipeline facilitate high-resolution imaging, microscopic evaluation, and real-time pathology-radiology correlations. A merger of these models may be the interface between mechanistic, proof-of-concept data and technologies oriented towards imaging, diagnostic, and clinical domains. Although promising, clinical equipoise is necessary as these technologies advance in breadth and scope. Moreover, from a practical standpoint, it is important to be mindful of logistics and cost-effectiveness of VR.52,64 The further development of validated measures of clinical significance, such as those presented, are of value in justifying widespread implementation of VR.
Virtual reality serves an adjunctive role in enhancing the sensitivity of DSA in detecting aberrant biomarkers of neurovascular remodeling in treated unruptured, intracranial aneurysms. The perspective from which our study results are examined should be prefaced by a closing remark. In domains that are clinically unexplored such as the role of VR in the present study, compelling but principled arguments that aim to explain rather than quantify observations are part of the natural history of clinically relevant, statistically significant findings.65–69 By building upon recent preliminary findings,66–69 a limited technology investigating evidence-based, yet clinically non-validated, neurovascular biomarkers with a proposed standardized methodology has been able to provide an exploratory answer to an interesting clinical question. Inherent statistical limitations in magnitude of effect, articulation of precision, and generalizability of results were partially mitigated by applied importance and provisional credibility based on theoretical coherence, methodological soundness, and results of limited statistical and yet to be determined clinical significance.70–73
Conflict of interest: The authors declare that they have no conflict of interest.
Informed consent: Informed consent was obtained from all individual participants included in the study.
Ethical approval: All procedures performed in the studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.
Acknowledgments: None
Funding: No funding was received for this study.
©2020 Fourcand, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.