- 1Department of Neurology, University Hospital, Ludwig-Maximilians-Universität, Munich, Germany
- 2German Center for Vertigo and Balance Disorders-IFB (Integriertes Forschungs- und Behandlungszentrum), University Hospital, Ludwig-Maximilians-Universität, Munich, Germany
- 3Graduate School of Systemic Neuroscience (GSN), Ludwig-Maximilians-Universität, Munich, Germany
- 4Department of Radiology, The Hospital for Sick Children, University of Toronto, Toronto, ON, Canada
- 5Department of Radiology, University Hospital, Ludwig-Maximilians-Universität, Munich, Germany
- 6Munich Cluster for Systems Neurology (SyNergy), Munich, Germany
In-vivo non-invasive verification of endolymphatic hydrops (ELH) by means of intravenous delayed gadolinium (Gd) enhanced magnetic resonance imaging of the inner ear (iMRI) is rapidly developing into a standard clinical tool to investigate peripheral vestibulo-cochlear syndromes. In this context, methodological comparative studies providing standardization and comparability between labs seem even more important, but so far very few are available. One hundred eight participants [75 patients with Meniere's disease (MD; 55.2 ± 14.9 years) and 33 vestibular healthy controls (HC; 46.4 ± 15.6 years)] were examined. The aim was to understand (i) how variations in acquisition protocols influence endolymphatic space (ELS) MR-signals; (ii) how ELS quantification methods correlate to each other or clinical data; and finally, (iii) how ELS extent influences MR-signals. Diagnostics included neuro-otological assessment, video-oculography during caloric stimulation, head-impulse test, audiometry, and iMRI. Data analysis provided semi-quantitative (SQ) visual grading and automatic algorithmic quantitative segmentation of ELS area [2D, mm2] and volume [3D, mm3] using deep learning-based segmentation and volumetric local thresholding. Within the range of 0.1–0.2 mmol/kg Gd dosage and a 4 h ± 30 min time delay, SQ grading and 2D- or 3D-quantifications were independent of signal intensity (SI) and signal-to-noise ratio (SNR; FWE corrected, p < 0.05). The ELS quantification methods used were highly reproducible across raters or thresholds and correlated strongly (0.3–0.8). However, 3D-quantifications showed the least variability. Asymmetry indices and normalized ELH proved the most useful for predicting quantitative clinical data. ELH size influenced SI (cochlear basal turn p < 0.001), but not SNR. SI could not predict the presence of ELH. In conclusion, (1) Gd dosage of 0.1–0.2 mmol/kg after 4 h ± 30 min time delay suffices for ELS quantification. (2) A consensus is needed on a clinical SQ grading classification including a standardized level of evaluation reconstructed to anatomical fixpoints. (3) 3D-quantification methods of the ELS are best suited for correlations with clinical variables and should include both ears and ELS values reported relative or normalized to size. (4) The presence of ELH increases signal intensity in the basal cochlear turn weakly, but cannot predict the presence of ELH.
Introduction
In-vivo non-invasive verification of endolymphatic hydrops (ELH) by means of delayed gadolinium (Gd) enhanced magnetic resonance imaging of the inner ear (iMRI) is rapidly developing into a standard clinical tool to investigate episodic vertigo (1–3). This is due to iMRI allowing pre-mortem detection of ELH for the first time (4, 5), demonstrating that ELH is not pathognomonic to Menière's disease (MD) (6–8), but rather a concomitant that can be found in various etiologies of episodic vertigo (9–13). Consequently, the clinical prevalence and pathophysiological significance of ELH has yet to be conclusively clarified. Understanding the underpinnings of the ELH syndrome requires a systematic investigation of pathologies involving endolymphatic space (ELS) changes as well as its base physiological condition.
Data acquisition protocols have undergone a continuous optimization of MR sequences (14, 15), as well as a steady minimization of procedural invasiveness (via a shift from intratympanic to intravenous application), duration and Gd dosage (16–18). A variety of cochlear and vestibular ELH quantification conventions have been suggested, including ELS semi-quantitative visual grading (19–25), manual measurement (26–28), semi-automatic (29, 30), and automatic algorithmic area ratio (AR), and volumetric segmentation (31, 32).
Given the plurality of approaches (some manual, some algorithmic), not all published results are inherently comparable. ELH features may vary greatly depending on ELS classification (for an overview, see Table 1) and data analysis choices. In this context, methodological comparative studies providing normalization and standard values between the methods and classifications used seem all the more important but still remain rare. On this note, this study aims to investigate the following questions:
(i) How variations in data acquisition protocols, such as Gd dosage or time delay, influence signal intensity (SI) and signal-to-noise ratio (SNR) within the ELS.
(ii) How ELH measures correlate with each other, as well as with clinical symptoms or neurophysiological testing.
(iii) How ELH influences SNR and SI within the ELS.
Materials and Methods
Setting and Institutional Review Board Approval
All data was acquired at the Interdisciplinary German Center for Vertigo and Balance Disorders (DSGZ) and the Department of Neurology of Munich University Hospital (LMU) between 2016 and 2019. Institutional Review Board approval was obtained before the initiation of the study (no. 641-15). All participants provided informed oral and written consent in accordance with the Declaration of Helsinki before inclusion in the study.
Study Population
One hundred eight consecutive participants [75 patients with Meniere's disease (MD) and 33 vestibular healthy controls (HC)] underwent delayed intravenous gadolinium-enhanced magnetic resonance imaging (iMRI) for exclusion or verification of ELH. The diagnosis of Meniere's disease (MD) was based on the Classification Committee of the Bárány Society 2015 (33). HC were inpatients of the Department of Neurology without symptoms or underlying pathologies of the peripheral and central vestibular and auditory system that underwent MRI with a contrast agent as part of their diagnostic workup and agreed to undergo iMRI sequences after 4 h. HC underwent audio-vestibular testing to confirm the soundness of their peripheral end organs. The reasons for their admission to the clinic included movement disorders (n = 6), epilepsy (n = 5), optic neuritis (n = 4), trigeminal neuralgia (n = 4), headache (n = 4), idiopathic facial nerve palsy (n = 3), viral meningitis (n = 3), subdural hematoma (n = 2), spinal inflammatory lesion (n = 1), and decompensated esophoria (n = 1). The laterality quotient for right-handedness was assessed with the 10-item inventory of the Edinburgh test (34, 35). The inclusion criterion was age between 18 and 85 years. The exclusion criteria were other neurological or psychiatric disorders, as well as any MR-related contraindications (36), poor image quality, or missing MR sequences.
Nomenclature
In the following, “ipsilateral” refers to the clinically leading side (or affected side) and “contralateral” to the opposite side (or non-affected side). In the case of patients presenting without a leading clinical side, a pseudorandom number generator [“Mersenne Twister” algorithm (37), uniform distribution] was used to generate a random number between 1 (=minimum value) and 9 (=maximum value). Even numbers meant “left side = ipsilateral side” and uneven numbers indicated “right = ipsilateral side.” “Vegetative symptoms” refers to nausea and/or vomiting due to the episodic vertigo attack. “Ear symptoms” includes attack-associated tinnitus, hearing loss, ear pressure, and/or ear pain both uni- and bilaterally that fit the criteria for MD. “Other ear symptoms” refers to non-MD ear symptoms.
Measurement of the Auditory, Semicircular Canal, and Otolith Functions
Diagnostic workup included a thorough neurological workup (e.g., history-taking, clinical examination), neuro-orthoptic assessment [e.g., Frenzel glasses, fundus photography, and adjustments of the subjective visual vertical (SVV)], video-oculography (VOG) during caloric stimulation and head impulse test (HIT), as well as ocular (o) and cervical (c) vestibular evoked myogenic potentials (VEMPs) and pure tone audiometry (PTA).
A tilt of the SVV is a sensitive sign of a graviceptive vestibular tone imbalance. SVV was assessed with the subject sitting in an upright position in front of a half-spherical dome with the head fixed on a chin rest (38). A mean deviation of >2.5° from the true vertical was considered a pathological tilt of SVV.
The impairment of the vestibulo-ocular reflex (VOR) in higher frequencies was measured by HIT (39) using high-frame-rate VOG with EyeSeeCam [(40), EyeSeeTech, Munich, Germany]. A median gain during head impulses <0.6 (eye velocity in °/s divided by head velocity in °/s) was considered a pathological VOR (41). Furthermore, canal responsiveness in lower frequencies was assessed by caloric testing with VOG, which was performed for both ears with 30°C cold and 44°C warm water. Vestibular paresis was defined as >25% asymmetry between the right- and left-sided responses (42). The caloric asymmetry index (AIC) was calculated based on the slow-phase velocity of the caloric nystagmus: .
Vestibular evoked myogenic potentials (VEMPs) are short-latency, mainly otolith-driven vestibular reflexes elicited by air-conducted sound (ACS), or bone-conducted vibration (BCV) and recorded from the inferior oblique eye muscle (ocular or oVEMPs) or the sternocleidomastoid muscle (cervical or cVEMPs). VEMPs were recorded with the Eclipse platform (Interacoustics, Middelfart, Denmark), as described previously (43, 44). Only those VEMP responses that were clearly discernible from background noise were included in the analysis. To avoid bias due to examiners, only the asymmetry index (AIo/cV) of VEMP amplitudes and latencies was analyzed in detail (45).
Delayed Intravenous Gadolinium-Enhanced MRI of the Inner Ear
Data Acquisition
Four hours after intravenous injection of a standard dose (0.1–0.2 mmol/kg body weight, i.e., 0.1 −0.1 mmol/kg body weight) of Gadobutrol (Gadovist®, Bayer, Leverkusen, Germany), MR imaging (MRI) data were acquired in a whole-body 3 Tesla MRI scanner (Magnetom Skyra, Siemens Healthcare, Erlangen, Germany) with a 20-channel head coil. We used a 3D-FLAIR sequence to differentiate endolymph from perilymph and bone, and a CISS sequence to delineate the total inner ear fluid space from the surrounding bone. The T2-weighted, three-dimensional, fluid-attenuated inversion recovery sequence (3D-FLAIR) had the following parameters: TR 6,000 ms, TE 134 ms, TI 2,240 ms, FA 180°, FOV 160 × 160 mm2, 36 slices, base resolution 320, averages 1, acceleration factor of 2 using a parallel imaging technique with a generalized auto-calibrating partially parallel acquisition (GRAPPA) algorithm, slice thickness 0.5 mm, acquisition time 15:08 min. The high-resolution, strongly T2-weighted, 3D constructive interference steady state (CISS) sequence of the temporal bones was performed to evaluate the anatomy of the whole-fluid-filled labyrinthine spaces and had the following parameters: TR 1,000 ms, TE 133 ms, FA 100°, FOV 192 × 192 mm2, 56 slices, base resolution 384, averages 4, acceleration factor of 2 using GRAPPA algorithm, slice thickness of 0.5 mm and acquisition time 8:36 min. The presence of ELH was observed on the 3D-FLAIR images as enlarged negative-signal spaces inside the labyrinth, according to a previously reported method (18, 46).
Signal Quality Assessment
Signal quality was validated using signal-to-noise ratio (SNR) and signal-homogeneity (SH) in different regions of interest (ROIs). ROIs were labeled in the left and right inner ear within the “endolymph” and “perilymph” fluid, “cochlear basal turn,” as well as in the surrounding tissue or subject matter, such as the “petrous bone,” “cerebellum,” “medulla,” and “air.”
In detail, the endolymph ROI consisted of 0.6 mm2 circular 2D-selections of the left/right utricle. The perilymph ROIs consisted of multiple 0.6 mm2 circular 2-D selections in the perilymphatic space (PLS) on both sides and were spread within the inner ear to obtain a signal intensity map. Said selections were placed in the vestibulum, twice inside the basal cochlea turn, the apex cochleae, the horizontal semicircular canal (hSCC) as well as the posterior SCC (pSCC). ROIs in the surrounding tissue or subject matter (“petrous bone,” “cerebellum,” “medulla oblongata,” and “air”) consisted of 60.8 mm2 circular selections. Signal intensity extraction (mean, minimum, and maximum) was performed on axial slices of the FLAIR raw images via the “Analyze Regions” plugin of the “MorphoLibJ toolbox” (47) within ImageJ (48).
SNR was calculated in each ROI as , i.e., the fraction of mean signal intensity in an ROI S(ROI), and the standard deviation (STD) of the region labeled “air,” std(air). The label “air” was defined as “MRI signal measure of background variations in the signal devoid of fluid.” In other words, a region's SNR was calculated as a mean signal relative to the extent of the background variation.
The signal's statistical homogeneity was examined between ROIs for each group, and between groups for each ROI. SH was defined as the identical distribution of two samples except for shifts and scaling of the overall distribution. The median of each sample was removed and the interquartile range was scaled to the value of one. The two samples were then compared using the minimum statistical energy [minEn; (49)] and the maximum mean discrepancy [MMD; (50)], whilst adding 10,000 permutations with a threshold of maximally one failed test to reach statistical significance. Consequently, two samples were deemed to have different distributions if they diverged in shape, either due to kurtosis, skewness or the extent, and number of outliers. Note that no correction for multiple testing was applied in these tests in order to be more sensitive toward violations of SH, i.e., significant differences.
Semi-quantitative Grading of the Endolymphatic Space
Semi-quantitative (SQ) grading of the endolymphatic space (ELS) was performed independently by three experienced head and neck radiologists or neurologists (BE-W, VK, and JG) who were blinded to the clinical patient data. Rater statistical homogeneity was calculated just as the signal's statistical homogeneity. The ELS's characterization in the vestibulum and cochlea was based on criteria previously described (22) and can be viewed in Table 1 and is described in further detail in Figures 1A,B, grade 0–3.
Figure 1. Semi-quantitative (SQ) grading used for the endolymphatic hydrops. The vertical columns show the different semi-quantitative (SQ) grades from 0 to 3 used [cf. Table 1; according to a classification first described in Kirsch et al. (9)]. The horizontal rows each give an overview of how each grade looks (A), in FLAIR raw data (B), after VOLT processing (C), as used for 2D quantification (D), and as used for 3D quantification (E). A detailed description of each grade is given in the second paragraph of subsection ‘Semi-quantitative Grading of the Endolymphatic Space'.
The characterization describes a 4-point grading for the cochlear and vestibular ELH. The cochlear grading is done on the midmodiolar level (19) and the vestibular grading on the inferior part of the vestibulum, where the left semicircular canal (L-SCC) is still visible (19). The cochlear grading can be thought of as a fusion of previously described grading suggestions (21–24). Grade 0 (no vestibular ELH) can be reduced to “X-mas tree built from circles that are divided by “very thin, clear, hypointense lines” [cp. also Figure 1 of (23), Figure 1A in (21)] that represent the non-enhanced ELS (scala media) between the enhanced PLS (scala vestibuli and tympani). Grade 1 (mild cochlear ELH) can be reduced to “X-mas tree with lights,” where the ELS is slightly enlarged and indirectly visible as a nodular black cut out of the scala vestibuli [cp. further Figure 2 in (24), Figure 1B in (21)]. Grade 2 (marked cochlear ELH) can be reduced to “X-mas tree with X-mas balls,” where the ELS is bulging into the scala tympani whilst giving the PLS a semicircular appearance [cp. Figure 3 in (24)]. Grade 3 (severe cochlear ELH) can be reduced to “X-mas tree with garlands,” where the severely distended ELS is causes a flattened appearance of the PLS [cp. also Figure 4 in (24), Figure 5A in (21), Figure 1C in (20)]. The vestibular grading is a fusion of previously described grading suggestions (21–23). Grade 0 (no vestibular ELH) can be reduced to “sacculus < utriculus,” where the otolith organs are distinguishable and the sacculus is smaller than the utriculus [cp. also Figure 2A in (21), Figure 6 in (23)]. Grade 1 (mild vestibular ELH) can be reduced to “sacculus≥utriculus,” where the sacculus is as large or larger than the utriculus [cp. also Figure 2B in (21), Figure 9 in (23)]. Grade 2 (marked vestibular ELH) can be reduced to “sacculus and utriculus are confluent,” where the otoliths organs are no longer distinguishable with a surrounding PLS rim [cp. also Figure 2C in (21), Figure 7 in (23)]. Grade 3 (severe vestibular ELH) can be reduced to “otolith organs not distinguishable” with no PLS visible [cp. also Figure 2D in (21), Figure 8 in (23)].
2D- and 3D-Quantification of the Endolymphatic Space
Segmentation of the total fluid space (TFS) was based on a recently proposed (Ahmadi et al., under review) and pre-trained volumetric deep convolutional neural network (CNN) with V-net architecture (51) that was deployed via the TOMAAT module (52) in 3D–Slicer toolbox [version 4.11 (53)]. ELS and PLS were differentiated within the TFS using Volumetric Local Thresholding [VOLT; (31)] using ImageJ Fiji (48) with the “Fuzzy and artificial neural networks image processing toolbox” (54) and the “MorphoLibJ Toolbox” (47).
The resulting 3D volume can be regarded as a probabilistic map of the inner ear, which includes the classification into its two different compartments (ELS and PLS). The final classification strongly depends on the chosen cutoff. Based on empirical observations (31), 2D- and 3D-quantifications were examined at three cutoff variations (c6, c8, and c10). Each cutoff matches a percentage of positive classifications. For example, cutoff 6 (c6) corresponds to 79.2%, cutoff 8 (c8) to 70.8%, and cutoff 10 (c10) to 62.5% classifications into endolymphatic space. Examples of the pipeline outputs can be viewed in Figure 1C.
2D-quantification was done on axial slices of the VOLT volume. The mid-modiolar level was chosen for the cochlea and the inferior part of the vestibulum where the lateral semicircular canal (L-SCC) is still visible was selected for the vestibulum. However, the majority of volumes allowed for both a cochlear and vestibular measurement on the same slice. Easier visual selection was enabled by a look-up-table (LUT, “phase”) included in ImageJ that was applied to the VOLT volumes. An example can be seen in Figure 1D. Areas were then measured using the “Analyze Regions” plugin which is part of the “MorpholibJ Toolbox” (47).
3D-quantification was done on the VOLT volume that included the entire inner ear. The cochlear volume was cropped using a cylindrical volumetric selection and applied to the VOLT volumes. The volume of the vestibulum including otolith organs and semicircular canals arose from subtracting the cochlear volume from the inner ear VOLT volume. Measurements were performed using the “Analyze Regions (3D)” plugin of the “MorpholibJ Toolbox” (47). A visualization can be seen in Figure 1E.
Parameters Derived From Endolymphatic Space Measures
The ELS ratio, , was calculated for 2D- and 3D-quantification of the ELS analogous to the area ratio (AR) in previous classification conventions (19–21). ER indicates the relative size of the ELH to the TFS and as such is independent of the absolute size which might differ between subjects (for example due to body size).
ELS symmetry between both inner ears was assessed via the ratio of ELS side differences Diff − ER [%] = ERi − ERc, where ERi and ERcare the respective ipsilateral and contralateral ELS ratios in percent relative to the TFS. Another parameter was the asymmetry index, , where ELSi is the semi-, 2D- or 3D-quantification of the ipsilateral ELS and ELSc of the contralateral ELS. The asymmetry index can be interpreted as a normalized difference and as such is also independent of the individual TFS.
Areas and volumes were normalized according to their TFS,
if c/v/a > c/v/a + 2.5 × std (TFS), where “e” is the individual value and “mean” is the mean of the respective group (HC or MD). For an overview of TFS, see Figure 6.
Statistics and Validation Parameters
All statistics were implemented with self-written scripts in MATLAB version 7.19.0 (R2019b) using the “Statistics and Machine Learning” toolbox provided with MATLAB (Natick, Massachusetts: the MathWorks Inc.). ELS quantification measures were validated and compared using parameters describing different characteristics on different levels (i.e., between groups, ELS analysis methods, and diagnostic methods) and between different entities (i.e., inter-rater, inter-threshold, and inter-ROI). Parameters considered the ordering of subjects between samples (concordance), Spearman correlations between samples (rank-correlation), the form of the distribution of samples via “minimum statistical energy” and “maximum mean discrepancy” (statistical homogeneity), and covariance between samples via ANCOVA (analysis of covariance). All statistical tests used multiple comparison correction, if multiple tests (e.g., more than two regions or two thresholds) were compared independently with each other. The FWE level was set at p = 0.05/N with N being the number of tests (e.g., regions, thresholds), i.e., Bonferroni correction.
Influence of Gd Dosage, Gd Time Delay on SNR and via SNR on SQ Grading, 2D- or 3D-Quantitatification
The influence of Gd dosage and time delay (from Gd injection to MR measurement) on the SNR and signal intensity (SI), as well as SNR, Gd dosage and time delay on SQ grading, 2D- or 3D- quantification measures was evaluated using ANCOVA modeling. The model included interaction of the group with each individual variable as well as the interaction of group, dosage and time delay variables. Additionally, covariates of no interest, such as age and BMI, were included. In other words, we checked whether SNR, Gd dosage, and Gd time delay each had an influence on the ELH measure in question, as well as the interaction of Gd dosage and Gd time delay, allowing for the possibility that the relationship might be different for each group.
Interrelations Between SQ Gradings and 2D- or 3D Quantification
Statistical Homogeneity
Statistical homogeneity between SQ grading, and 2D- or 3D- quantification methods between groups was, in principle, calculated in the same way as the signal statistical homogeneity (cf. signal quality validation). First, the median of each group was removed and the interquartile range was scaled to the value of one. The two groups were then compared using minEn and MMD test statistics whilst using 10,000 permutations between groups. Any instance of a random permutation with a higher test-statistic than the unpermuted groups was considered a failure. The groups were deemed statistically homogeneous if at most one test failed, otherwise the groups were deemed inhomogeneous, as they could be distinguished based solely on their distribution shape (kurtosis and skewness or the extent and number of outliers). Note that no correction for multiple comparisons was performed here in order to be more sensitive to violations of SH, i.e., significant differences.
Rater Repeatability and Reliability
Repeatability and reliability of the three different raters for SQ grading, as well as of the three different thresholds (c6, c8, and c10) for 2D- or 3D-quantification were measured using rank-based correlations and Kendall's W measure for concordance (55). This assessment shows whether the ordering of subjects between raters is similar and therefore can be assumed to be repeatable over the raters. Furthermore, we compared ratings by subtracting the SQ grading scores between raters to see if the extent of differences in rating values differed. Correction for multiple comparisons, i.e., multiple tests was done over data types (SQ, 2D and 3D), therefore p(FWE) = 0.05 was set to p = 0.05/N with N = 3 for the three data types.
Interrelations Between SQ Grading, 2D- and 3D-Quantification
Interrelations between SQ grading, 2D- and 3D-quantification were examined via Spearman, i.e., rank-based correlations. Significant rank correlations indicated that the ordering of subjects was very similar or concordant across these measures. Rank-correlation was used so that linear as well as non-linear relationships could be examined and the gradings (ordinal measures) could be related to the quantitative measures. Correction for multiple comparisons p(FWE) = 0.05/N, i.e., Bonferroni correction, was done over all pairs of correlations in each correlation matrix, i.e., for SQ- × -2D quantification and SQ- × -3D quantification N = 12- × 12 = 144 and for the correlation of asymmetry indices N = 6- × 6 = 36.
Influence of Thresholds on Quantitative Measures
The influence of VOLT thresholds (c6, c8, and c10) on group differences was assessed using general linear model (GLM) based two-sample t-tests (including age as a covariate of no interest). The resulting slopes for the effects of thresholds and their standard errors were used to calculate t-statistic values for each group comparison at each threshold. Furthermore, a slope difference test (56) was used to check if group differences depended on the VOLT thresholds. A slope difference test compares differences in slopes with standard errors for the group differences across thresholds to determine if group differences depended on the cutoff-threshold. Correction for multiple comparisons p(FWE) = 0.05/N, i.e., Bonferroni correction, was done for three tests of between threshold comparisons resulting from three thresholds (c6,c8,c10), i.e., c6-vs-c8, c6-vs-c10, and c8-vs-c10, and therefore N = 3.
Covariance of Clinical Measures and iMRI
Clinical (e.g., disease duration, number of attacks) and diagnostic measures (e.g., HIT, calorics, and VEMPs), as well as parameters derived from ELS measures in SQ gradings and 2D- and 3-D quantifications (ER, Diff-ER, and AI) were included in an analysis of covariance (ANCOVA). An overview of clinical symptoms and diagnostic measures can be viewed in Table 2. Furthermore, the analysis accounted for categorical variables, such as symptoms like headache, and continuous covariates, such as body mass index (BMI) and the age of the patients. For detection of diverging trends between MD or HC, parameters derived from ELS measures were allowed interactions with the group. That means each group was allowed to have a different trend in the model. We used Bonferroni-correction for the post-hoc assessment of the individual factors in the ANCOVA.
Influence of ELH Presence on SNR and SI
The influence of the presence and extent of ELH on SNR and SI was examined with two approaches. First, SNR and SI data were investigated using classifications derived from SQ grades and 3D-quantification measures. The SQ grades were used to distinguish between “no ELH” and “ELH present,” while the 3D-quantification was used to distinguish between “low/small ELH” and “high/large ELH.” For the classification using SQ grades, all grades equal to zero (SQ grade = = 0) were allocated to “no ELH” and the rest to “definite ELH present.” For the classification using 3D-quantification, data values below the median were in the “low ELH” class and data values above the median in the “high ELH” class. The SNR and SI data were analyzed using two-sample t-tests and Wilcoxon rank-sum tests for differences from these classifications. The two tests (i.e., a parametric and non-parametric test), were used to ensure that any of the significant differences found were not purely dependent on the assumed distribution. Correction for multiple comparisons p(FWE) = 0.05/N, i.e., Bonferroni correction, was done for five tests between regions (split by ELH) comparisons, i.e., N = 5 (see Figures 5A,B).
Then, the inverse question was asked. This time SQ and 3D-quantification values were compared following SI or SNR value classification and then analyzed accordingly for differences with two-sample t-tests and Wilcoxon rank-sum tests. For both SI and SNR classification, “low SI or SNR class” was defined by their values below the respective median, and “high SI or SNR class” by their values above the respective median. Correction for multiple comparisons p(FWE) = 0.05/N, i.e., Bonferroni correction, was done separately for the test between SNR (split by ELH 3D-quantification) and 3D-quantification (split by SI). The number of tests for the SNR comparison was N = 2, and the number of tests for ELH 3D-quantification was N = 4 (see Figures 5C,D).
Results
Descriptive Statistics
Seventy-five MD patients (35 females; aged 22–81 years, mean age 56.6 ± 14.9 years; 97% RH) and 33 HC participants (20 females; aged 20–84 years, mean age 42.1 ± 18.9 years; 94% RH) were included in the study. An overview of the most important clinical features in MD compared to HC can be seen in Table 2. An overview of the ELS grading for HC and MD can be viewed in Table 3.
Influences of Signal Quality on ELS Quantification Methods (i)
• The signal intensity (SI) of each region of interest significantly (FWE-corrected, p < 0.05) depended on Gd dosage (range: 0.08–0.28 ml/kg; mean ± std: 0.17 ± 0.05 ml/kg; 48% of subjects got 1 dose, 12% 1.5 doses; 40% got 2 doses of 0.1 ml/kg) and Gd time delay (range: 2 h and 51 min to 5 h and 20 min; mean ± std: 4 h and 24 ± 25 min; 25% are between 3 and 4 h, 50% are between 4 and 4 h and 31 min and the remaining 25% were longer). However, the effect sizes (eta-squared) were small (5–12%).
• The mean SNR (range: 24.8–130.49; mean ± std: 64.82 ± 20.64) was significantly related to Gd dosage and time delay (FWE-corrected, p = 0.006), but only 6.8% of the total variance (r-squared) could be explained. If looked at separately, Gd dosage (4.2%, FWE-corrected, p = 0.03) and Gd time delay (5.3%, FWE-corrected, p = 0.02) explained even less of the variance. See Figure 2 for an overview of the minor influence of the iMRI acquisition parameters on SNR.
• The mean SNR was significantly different between the MD and HC group (p < 0.05). SNR asymmetry between left and right ear was not significantly related to Gd dosage, Gd time delay, or Gd dosage × Gd time delay.
• SQ gradings and 2D- or 3D- quantifications were not significantly related to Gd dosage, Gd time delay, Gd dosage × Gd time delay interaction, SI or SNR (FWE corrected, p ≤ 0.05). There were some simple significant relationships (p < 0.05 uncorrected) for the iMRI variables with Gd dosage and SNR, but all these relationships were small in effect size (around 0.5–5% omega squared).
Figure 2. Influence of gadolinium (Gd) dosage and Gd time delay on the signal-to-noise ratio (SNR) (belonging to question i). The scatter plot of the signal-to-noise ratio (SNR) over the gadolinium (Gd) dosage (in ml, plot in the top left), of the SNR over the time delay (in minutes, plot in the top right) and of the SNR over the interaction of Gd dosage × time delay (ml × minutes, plot in the bottom). SNR data points are plotted as black crosses and the trend lines of the fitted model are plotted as red lines. The SNR depended significantly on dosage, time delay as well as the interaction of dosage × time delay. However, the effect size was very small (4.2, 5.3, and 6.8% of explained variance, respectively). The SNR displayed here “SNR mean” is the mean of the SNRs calculated for each region of interest, i.e., the two inside the basal cochlear turn (CBT), the apex cochleae (AC), the horizontal semicircular canal (hSCC) as well as the posterior SCC (pSCC). For each region, the SNR is the mean signal divided by the standard deviation of the region labeled “air.” The left and the right ear are averaged for each region, before regions are averaged to form “SNR mean”.
Interrelations Between ELS Quantification Methods (ii)
• Inter-rater SQ gradings (R1-3) were statistically homogeneous, as were 2D- and 3D- quantification values including ipsi- and contralateral or cochlea and vestibulum.
• Inter-rater SQ gradings (R1-3) were highly matched in the vestibular (v) and cochlear (c) part of the inner ear for all subjects (HC, MD) and slightly less for MD only. The results can be viewed in Table 4 (column SQ) and suggest a high but imperfect reproducibility due to remaining variability.
• Inter-threshold (c6, c8, and c10) 2D- and 3D-quantification was highly concordant. These results can be viewed in Table 4 (column 2D and 3D) and indicate an almost perfect agreement over VOLT thresholds with a basically perfect reproducibility.
• SQ grades correlated strongly with 2D-quantification values (range of correlation from 0.3 to 0.7) and 3D-quantification values (range of correlation from 0.3 to 0.7). The correlations of 2D- and 3D-quantification values with SQ grades was mainly driven by the MD group, due to the higher variability within the group, compared to HC group which did not vary much in grades or 2D- and 3D-quantification values (cp. Figure 3, plots on the left and in the middle).
• 2D- and 3D-quantification correlated substantially (range of correlation from 0.3 to 0.8) for the total inner ear, cochlea, and vestibulum on both the ipsilateral and contralateral side. However, there were no significant correlations of the ipsilateral with the contralateral sides (cp. Figure 3, plots on the right). AISQ (asymmetry-index of SQ quantification) correlated significantly (range of correlation from 0.3 to 0.7) with AI2D and AI3D (asymmetry-indices of 2D- and 3D-quantification) except for the cochlear AI in the 2D- and 3D-quantifications in the c6-cutoff ( and , cp. Figure 3).
• Inter-rater SQ grading differences did not differ strongly between R1-3. Figure 4 shows the results in more detail. For the vestibular part, the percentage of ratings that agreed, i.e., showed zero differences, was 54.6% (R2-R1), 50% (R3-R1), and 67.6% (R3-R2), while the percentage of differences of maximally one grade apart was 85.2% (R2-R1), 90.7% (R3-R1), and 98.2% (R3-R2). For the cochlear part, the percentage of ratings that agreed was 53.7% (R2-R1), 53.7% (R3-R1), and 71.3% (R3-R2), while the percentage of differences of maximally one grade apart was 88.9% (R2-R1), 90.7% (R3-R1), and 97.2% (R3-R2).
• Inter-threshold 2D- and 3D-quantification measures were statistically homogenous and showed group differences for each threshold.
• Clinical variables correlated with symmetry parameters derived from SQ grading and 2D- or 3D-quantification values such as the asymmetry index (AI) or the plain ELH difference between ipsilateral and contralateral side for the inner ear, vestibulum, and cochlea. The AI worked for un-normalized data and its results were comparable to the normalized data, while the pure differences between ipsilateral and contralateral sides were only useful when the data was first normalized to the fraction ER [%] of the total fluid space (TFS, cf. legend on standard values). This indicates that relative proportions of both ears, and the relative size of ELH are most useful for predicting quantitative clinical data from iMRI measures. Fittingly, vestibular AI for the 3D-quantification data explained 35% of the variance of the number of attacks in the 3 months prior to the examination and another 16% of variance could be explained by the AI for the 3D-quantification data of the whole of the inner ear (vestibular and cochlear parts combined). A more detailed clinical study and discussion can be found in another work (57).
Figure 3. Significant correlations between endolymphatic hydrops (ELH) quantification methods (belonging to question ii). The top row shows the correlations between semi-quantitative (SQ) gradings (x-axis) and 2D- or 3D- quantification values (y-axis) for vestibulum (v) or cochlea (c), and for the ipsilateral (ipsi) or contralateral side (contra). In addition, SQ gradings are rater-specific (R1-R3), and 2D- or 3D- quantification values are cutoff-specific (c6, c8, and c10). The higher the significant results (p < 0.05 FWE-corrected) correlated, the more they are colored in yellow (thresholded to 0–0.8). Overall, SQ gradings and 2D- or 3D-quantification values correlated with the respective other method, although to a higher extent on the ipsilateral side (on the contralateral side in c10) and vestibular part of the inner ear. The bottom row shows the corresponding correlations between methods (SQ gradings and 2D- or 3D- quantification values) for the respective asymmetry-index (AI).
Figure 4. Inter-rater and -threshold between ELH quantification methods (belonging to question ii). Shown are differences between the three raters. The differences between raters are shown as percentages of the total number of subjects rated. Most grades between raters agree (no difference; in blue), and the next largest difference was by 1 grade (in green), the remaining differences between raters were mostly 2 grades (in yellow) and rarely 3 grades (in red) apart.
Influences on Signal Quality (iii)
• There were significant differences in SI due to the presence of ELH in the following ROIs: cochlear basal turn [p(t−test) = 0.0009 and p(rank−sum) = 0.0003], apex cochleae [p(t−test) = 0.002 and p(rank−sum) = 0.001], hSCC [p(t−test) = 0.038 and p(rank−sum) = 0.022], and pSCC [p(t−test) = 0.046 and p(rank−sum) = 0.018]. Generally, higher ELH 3D-quantification values had higher SI values. However, due to a significant spread, SI could not distinguish the presence of ELH from the absence of ELH (tested by means of the split of SI values based on defining “absence of ELH” as an SQ grading equal to zero and “presence of ELH” as all grades higher than zero). Fittingly, the opposite approach (splitting ELH values by SI brightness) did not show significant differences. For an overview, please see Figure 5. The signal intensity SI was significantly different between the MD and HC group for both ROIs in the cochlear basal turn, but not for the cochlear apex, hSCC, or pSCC (p < 0.05 FWE). The group differences in iMRI variables between the MD and HC groups persisted after removing effects of Gd dosage, time delay, and SNR, indicating that iMRI assessment was not significantly affected by the differences in Gd dosage, time delay, and SNR in the present dataset.
• SNR was not influenced significantly by the presence or absence of an ELH. Selecting SNR values for all SQ grades = 0 (“absence of ELH”) and comparing them with the remaining SNR values (where SQ grades >0, “presence of ELH”) led to two-sample t-test p = 0.99 and two-sample rank-sum test p = 0.94. Furthermore, comparing the SNR values for low ELH values (3D-quantification values below the median) with SNR values for high ELH values (3D-quantification values above the median) did not show any significant differences in SNR (two-sample t-test p = 0.31 and two-sample rank-sum test p = 0.45). Analog to this, splitting ELH values due to low SNR values vs. high SNR values did not result in significant differences [p(t−test) = 0.66 and p(rank−sum) = 0.47 on the ipsilateral side and p(t−testtest) = 0.2 and p(rank−sum) = 0.16 for the contralateral side]. For an overview, see Figure 5.
Figure 5. Influence of endolymphatic hydrops (ELH) on signal intensity (SI) or signal-to-noise ratio (SNR) and vice versa (belonging to question iii). The top two rows (A,B) show the signal intensity (SI) in various areas of the inner ear. The SI values are split by either (A) 3D-quantification values (below or above median, indicated by a downward-arrow or an upward-arrow) or (B) endolymphatic hydrops (ELH) being absent [“–”] or present [“+”]. Here, “ELH absence” is defined as the average semi-quantitative (SQ) rating being zero and “ELH presence” is defined as the average rating being non-zero. Significant differences are indicated by black lines and p-values for the group differences evaluated with two-sample t-test or rank-sum test (**p < 0.001 and *p < 0.05). The third row (C) shows the mean SNR split by 3D-quantification values (below or above median; indicated by a downward-arrow or an upward-arrow), as well as the mean SNR split by ELH being absent [“–”] or present [“+”], where presence is defined as ELH being non-zero and absence as ELH rating being zero. The group differences were evaluated with the two-sample t-test and the rank-sum test; neither of these was significant. The fourth row (D) shows the 3D-quantification values split by mean SI values (below or above median; indicated by a downward-arrow or an upward-arrow), as well as the 3D-quantification values split by SNR values (below or above median; indicated by a downward-arrow or an upward-arrow). The group differences were evaluated with the two-sample t-test and the rank-sum test; neither of these was significant.
Standard Values
• Areas and volumes were normalized according to their TFS (total fluid space/surface) and can be viewed in Figure 6.
• Our calculations showed that the chosen threshold did not change the group differences between MD and HC. The grading-specific 2D- and 3D-quantification values, the TFS values and resulting ratios can be seen in Figures 6, 7. Furthermore, we show the relationship of 2D- and 3D-quantification for the vestibular and cochlear part broken down by SQ grades in Figure 8 and Table 5. While grades increase, one can observe that 2D- as well as 3D-quantification increased.
• ELH 3D-quantification values (see also Figure 6): The medians (ipsilateral, contralateral) of the vestibular data were (15 mm3, 11 mm3) for the MD group and (12 mm3, 12 mm3) for the vestibular healthy control (HC) group. The medians of the cochlear data were (5.4 mm3, 4.6 mm3) for the MD group and (4.6 mm3, 4.7 mm3) for the HC group.
• Total fluid space (TFS) 3D-quantification values of the vestibular and cochlear part of the inner ear [in mm3] that were used for normalizing the data of each individual to generate ELS ratio (ER) [%], the percentage of the TFS occupied by the ELH. The medians (ipsilateral, contralateral) of the vestibular TFS data were (182 mm3, 183 mm3) for the MD group and (179 mm3, 181 mm3) for the HC group. The medians of the cochlear TFS data were (93 mm3, 93 mm3) for the MD group and (88 mm3, 90 mm3) for the HC group.
• Therefore, the medians (ipsilateral, contralateral) of the vestibular ER [%] data were (8.8%, 6.3%) for the MD group and (6.3%, 6.5%) for the HC group. The medians of the cochlear ER [%] data were (6.0%, 5.1%) for the MD group and (5.1%, 5.0%) for the HC group.
• ELH 2D-quantification values (see also Figure 7): The medians (ipsilateral, contralateral) of the vestibular data were (3.6 mm2, 2.2 mm2) for the MD group and (2.3 mm2, 2.4 mm2) for the vestibular healthy control (HC) group. The medians of the cochlear data were (1.2 mm2, 0.8 mm2) for the MD group and (1.1 mm2, 1.0 mm2) for the HC group.
• Total fluid surface (TFS) 2D-quantification values of the vestibular and cochlear part of the inner ear, in [mm2] that were used for normalizing the data of each individual to generate ELS ratio (ER) [%], the percentage of the TFS occupied by the ELH. The medians (ipsilateral, contralateral) of the vestibular TFS data were (20.1 mm2, 20.8 mm2) for the MD group and (20.9 mm2, 20.4 mm2) for the HC group. The medians of the cochlear TFS data were (16.4 mm2, 16.1 mm2) for the MD group and (15.9 mm2, 15.8 mm2) for the HC group.
• Therefore, the medians (ipsilateral, contralateral) of the vestibular ER [%] data were (19.2%, 10.5%) for the MD group and (10.7%, 11.7%) for the HC group. The medians of the cochlear ER [%] data were (7.9%, 5.3%) for the MD group and (6.3%, 6.3%) for the HC group.
Figure 6. Normalization and group differences for 3D-quantification of endolymphatic hydrops (ELH). 3D-quantification values of endolymphatic hydrops (ELH) in the vestibular and cochlear part of the inner ear [in mm3] (top). 3D-quantification values of the total fluid space (TFS) of the vestibular and cochlear part of the inner ear [in mm3] (middle row) that were used for normalizing the data of each individual to ELS ratio (ER) [%], the percentage of the TFS occupied by the ELH (bottom row).
Figure 7. Normalization and group differences for 2D-quantification of endolymphatic hydrops (ELH). 2D-quantification values of endolymphatic hydrops (ELH) in the vestibular and cochlear part of the inner ear, in [mm2] (top row). 2D-quantification values of the total fluid surface (TFS) of the vestibular and cochlear part of the inner ear, in [mm2] (middle row) that were used for normalizing the data of each individual to ELS ratio (ER) [%], the percentage of the TFS occupied by the ELH (bottom row).
Figure 8. Semi-quantitative (SQ) grade-specific 2D- and 3D-quantification values. These four plots show the quantitative endolymphatic hydrops (ELH) measures (3D-quantification in the top row and 2D-quantification in the bottom row) split up along the semi-quantitative visual grading steps (from grade 0 to grade 3 in steps of 1). This means each boxplot with data points overlaid shows the distribution of quantitative measures that correspond to the respective grading. The grading here is of the average over the three raters, rounded to the nearest integer to preserve gradings going from 0 to 3 in steps of 1. The top row shows the distributions for the 3D-quantification of ELH for vestibular and cochlear parts of the inner ear, left and right, respectively, while the bottom row shows the distributions for the 2D-quantification of ELH for vestibular and cochlear parts of the inner ear.
Discussion
This methodological study with 108 participants (75 MD, 33 HC) focused on comparability and parametrization of different ELS quantification methods (SQ grading of three raters, 2D- or 3D-quantification of three cutoffs) used in iMRI and their (i) interrelations with subtle variations in data acquisition protocols (that influence SNR or SI); (ii) correlations to each other, clinical symptoms, or neurophysiological testing; and (iii) the influence of ELH on signal quality. The results were as follows: (i) Within the range of 0.1–0.2 mmol/kg (mean ± std: 0.16 ± 0.05 mmol/kg) Gd dosage and a 3 h 41 min to 5 h 19 min (mean ± std: 4 h 39 min ± 25 min) time delay, SQ gradings, and 2D- or 3D-quantifications were independent of signal intensity (SI) and signal-to-noise ratio (SNR), but they were found to be significantly related to Gd dosage and time delay themselves. (ii) The ELS quantification methods used were highly reproducible across raters (SQ gradings) or thresholds (2D- and 3D-quantification), although 3D-quantifications showed least variability in comparison to 2D-quantifications and SQ gradings. The relative proportions of both ears, and the relative size of ELH proved to be most useful for predicting quantitative clinical data from iMRI measures. (iii) ELH size significantly influenced SI but not SNR. In contrast, SI could not predict ELH size. In the following, results (i–iii) will be discussed.
Within a Specific Dosage and Time Delay Range ELS Quantification Methods Remain Independent of Signal Intensity (i)
The 3D fluid-attenuated inversion recovery (3D- FLAIR) imaging used has high sensitivity to low concentrations of Gd-based contrast agents (GBCA) in fluid compared with conventional T1-weighted imaging (58). In particular, the heavily T2-weighted 3D-FLAIR imaging with a long effective echo time is very sensitive to subtle T1 shortening and can detect low concentrations of GBCAs in the perilymphatic space after intravenous administration of a single dose of GBCA (18, 59, 60). In the tested Gd dosage and Gd time delay range (see above), at most weak influences on SNR and no influence on ELS quantification methods were found. It can therefore be assumed that, although Gd dosage and Gd time delay certainly have an influence on iMRI quality parameters, the sweet spot for ELH quantification by iMRI is within the range of the tested parameters. These results tie in well with earlier studies that showed strongest enhancement in 3D FLAIR sequences between 3 and 6 h (61), and optimally 4 h (62) after intravenous administration of 0.1 mmol/kg gadolinium diethylenetriaminepentaacetic acid (standard dose Gd-DTPA, Magnevist®), or 4.5 h (63) after intravenous administration of a standard dose gadoteridol (0.1 mmol/kg Gd-HP-DO3A, ProHance®), double dose (0.2 mmol/kg) Gd-HP-DO3A (64) or triple dose (0.3 mmol/kg) of Gd-DTPA-BMA, Omniscan® (65).
Another feature of the good performance within the chosen ranges may be the homogeneous distribution of the contrast agent in the entire volume of the inner ear (66, 67).
Further improvement of SNR and visualization in terms of rapid, morphological enhancement for analysis of the temporal and spatial distribution in the PLS of the inner ear can be achieved through careful selection of MR sequences (59, 68), combination (69, 70), and post-processing (14) of MR sequences, MR Gd complex (71), MR coil, and MR field strength (72).
Is There a Hierarchy Within ELS Quantification Methods? (ii)
In line with the only comparative methodological study of ELS quantification methods published to date in 11 participants (9 patients and 2 healthy controls) (26), SQ grading and 2D- or 3D- quantification methods were found to be reliable and useful for the diagnosis of endolymphatic hydrops. However, the degree of reliability based on comparisons between raters or thresholds increased from SQ grading to 2D- and again to 3D-quantification methods. The increase in repeatability corresponds to the decrease in dependency of human decision (visually > specific slice in 2D > whole volume in 3D) and increase of automatization and data points (semi-quantitative < area < volume).
Another aspect that makes relying solely on SQ grading tricky is the comparability of methods between different research groups, besides inter-rater disparities. SQ grading conventions (cf. Table 1) vary in grading resolution from three [in cochlea (19–21, 23) and vestibulum (19, 24, 25)] to four steps [in cochlea (22, 24, 25) and vestibulum (21–23)]. Accordingly, not all ELH grade results in cochlea or vestibulum correspond to each other due to the usage of different conventions [as an example grade 1 in (19, 20, 24)], or not at all [as an example (73)]. Based on either manually drawn (28, 74) regions of interest (ROIs) or a convolutional neural network (CNN) segmentation (32), 2D quantification methods already offer an increased comparability and variability of information. However, the comparability of the results remains limited by the slice selection for the calculation of the ratio and the differing slices emerging from slice planning or MRI setup (sequence type, slice thickness, slice resolution). Concerning these issues, 3D-quantification can be a solution (no slice selection, independent of slice planning) or at least an improvement (sequence type, slice thickness, slice resolution). In addition, more information (data points) enables better fitting of diagnostic and clinical parameters (75). Yet here, too, methodological variations affect reproducibility and availability of results. The critical points are the segmentation of the inner ear from the background [manually (29), via atlas (76, 77), or CNN (31); (Ahmadi et al., under review)] and the ELS and PLS from the TFS [manually (26), semi-automatic (29), automatic (31)], as well as the availability of the software solutions [commercial (26, 28, 29, 78) vs. open source (31)]. The less human-dependent and the more automated, the more reproducible the method in most cases. Therefore, the usefulness of the available quantification methods depends on its intended application. While visual SQ grading is highly useful in a clinical setting, automated 3D-quantification seems most suitable for research.
ELS Patterns in MD, HC (ii)
Significant differences between groups could especially be found for the ipsilateral (or affected) side of the MD group vs. HC group, as was already shown for 4-point [cochlea: (24, 25)] and 3-point [cochlea: (79); vestibulum: (80); sacculum (81)] ELS SQ grading, 2D-quantification [cochlea: (82); vestibulum: (82)] and 3D-quantification [cochlea: (30, 57); vestibulum: (30, 57)] results.
Clinical variables correlated highest with symmetry parameters derived from SQ grading and 2D- or 3D-quantification values such as the asymmetry index (AI) or the plain ELH difference between ipsilateral and contralateral side for the inner ear, vestibulum, and cochlea. Recent studies using ELS asymmetry indexes confirm this inclination (57). A more detailed clinical study and discussion can be found in another work (76). To date, correlations were found for SQ grades 3-point [electrocochleography (EcochG) (83, 84)] and 4-point ordinal cochlear scale [PTA (24–26, 79); auditory symptoms (85, 86); disease duration (24, 79); but not for the glycerol test (25)]. Furthermore, correlations were found for 3-point ordinal vestibular scale [cVEMP-side difference (SD) (24); PTA (11, 25, 87); oVEMP-amplitudes (88), but not with VOG during caloric stimulation SD (24, 89) or the glycerol test (25)]. SQ correlations coincided with 2D-quantification [cochlea: PTA (82); vestibulum: SP/AP ratio of ECoG (82)] and 3D-quantification [cochlea: PTA (26); vestibulum: duration of illness >30 months (26), side difference in response to caloric irrigation (57)] correlation results.
ELH Extent Influences Signal Intensity in the Basal Cochlear Turn (iii)
Zhang et al. (90) investigated 19 MD patients following double-dose iMRI and found that the signal intensity ratio of the cochlear basal turns in the affected ear was significantly higher than in the unaffected ear and that there was a positive correlation between the signal intensity ratio of the cochlear basal turn and the grades of cochlear and vestibular hydrops in the affected ear. The SNR was assessed and calculated manually according to (91) using the signal in perilymph of both cochlear basal turns and noise in coplanar circular 50 mm2 ROIs in the cerebellum. The interpretation of these findings was that increased permeability of the blood-labyrinth barrier (higher SNR) may play a role in the process of endolymphatic hydrops in MD.
The results of the current study suggest, however, a general pathophysiological effect tied to the extent of the ELH and not MD as a pathology, since higher ELH 3D-quantification values had higher signal intensity (SI) values in the cochlear basal turn, apex cochlea, and hSCC ROIs. Within MD, SI only in the cochlear basal turn was significantly higher on the ipsilateral side when compared to the contralateral side. The SI was generally different between the MD and HC groups, indicating an effect of ELH also on signal presentation. SNR differed between the MD and HC groups; however, the effect was small and the group differences in ELH were not significantly affected by SNR, indicating that the group differences are a persistent effect of the underlying condition and not related to the imaging settings that were used in the current study.
Normalization and Standardization of ELS Values
Clinical variables correlated better and more correctly with relative (AI) or normalized values [to the fraction ER [%] of the total fluid space]. This indicates that relative proportions of both ears, and the relative size of ELH are most useful for predicting quantitative clinical data from iMRI measures (57).
However, to date not many iMRI ELS values have been published in absolute (26, 30, 92) and relative sizes (26, 28, 85, 93–95); those that have been published were mostly group-specific but not grade-specific, and one grade-specific but relative (26). In Table 5, 2D- and 3D-quantifications, relative, and normalized to TFS are presented.
Recommendations for Future iMRI Studies
The following methodological recommendations for future studies can be derived from the present work and the current available literature:
• MR setup: Improved hybrid of reversed image of positive endolymph signal and native image of positive perilymph signal (iHYDROPS-Mi2) (15) or 3D-real inversion recovery (3D-real IR) (28, 96), highest possible MRI field strength (72), smallest possible isotropic voxel size, deep learning reconstruction denoising (14) if applicable.
• MR measurement: 4 h ±30 min after single-dose (0.1 mmol/kg) intravenous application (64) of Dotarem (Gd-DOTA, 100% morphological enhancement) or Gadovist (Gd-Do3A, 88% morphological enhancement) (71).
• SQ grading: 4-point ordinal scale (0: no hydrops, 1: mild hydrops, 2: marked hydrops, 3: severe hydrops) for cochlear and vestibular SQ grading (cp. Table 1), preferably with a level of evaluation reconstructed to distinctive anatomical fixpoints.
• Scalar ELS values: 3D-quantification, optimally using algorithm-based segmentation of both the TFS (Seyed-Ahmad et al., under review); (76, 77) and ELS (31, 32) should be included. 2D-quantification if 3D-quantification is not available. Reported values should be normalized for TFS size.
• Correlations with clinical variables should include both ears and are most promising in symmetry parameters, such as asymmetry-indices for un-normalized data and relative size ELS for normalized data.
Methodological Limitations and Outlook
There are methodological limitations in the current study that need to be considered in the interpretation of the data. First, despite the comparatively wide range of contrast agent dosage and delay time within this study, the results should be (to some degree) considered specific to the study's MR settings (MR sequence, MR contrast agent, intravenous application). Second, despite the extensive analyses within this study, it was not possible to try all, but only representations of the methods used in this study [SQ following Figure 1 and (22), 2D- and 3D-quantification using VOLT (31)]. Third, the study lacks histological confirmation of endolymphatic hydrops. However, the in-vivo acquisition of histological specimens in Menière's disease is currently not possible. Fourth, the size of the control group (n = 33) was small in comparison to the MD group (n = 105). However, due to findings of signal intensity in the dendate nucleus and globus pallidus on unenhanced T1-weighted MR images (97–99) that are still under investigation, measurements were restricted to inpatients of the Department of Neurology that underwent MRI with a contrast agent as part of their diagnostic workup and agreed to undergo iMRI sequences after 4 h.
Conclusion
The current comparative methodological study has shown that: (1) A Gd dosage of 0.1–0.2 mmol/kg after 4 h ± 30 min Gd time delay will provide sufficient SNR when using recommended MR sequences and contrast agents. (2) An agreed upon clinical SQ grading classification including a standardized level of evaluation reconstructed to anatomical fixpoints is needed to provide unambiguous comparability between labs. (3) 3D-quantification methods of the ELS using algorithm-based segmentation of the TFS and ELS seem to be best suited for research purposes. Correlations with clinical variables should include both ears and ELS values reported relative or normalized to size. (4) The presence of ELH increases signal intensity in the basal cochlear turn weakly, but cannot predict the presence of ELH.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Medical Faculty of the Ludwig-Maximilians-Universität, Munich, Germany. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
RB: conception and design of the study, analysis of the data, and drafting the manuscript. JG, EK, SB-B, and BE-W: acquisition, analysis of the data, and drafting of the manuscript. MD: conception and design of the study, drafting the manuscript, and providing funding. VK: conception and design of the study, acquisition, analysis of the data, drafting the manuscript, and providing funding. All authors contributed to the article and approved the submitted version.
Funding
This work was partially funded by the Friedrich-Baur-Stiftung (FBS), the German Federal Ministry of Education and Research (BMBF) in connection with the foundation of the German Center for Vertigo and Balance Disorders (DSGZ; grant number 01 EO 0901) and the German Research Foundation (DFG) in connection with the Munich Cluster for Systems Neurology (SyNergy; grant number: EXC2145, project ID: 390857198), and the German Foundation for Neurology (DSN).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This was part of the dissertation of JG. We thank K. Göttlinger for copyediting the manuscript.
Abbreviations
±, standard deviation; 2D, two-dimensional; 3D, three-dimensional; AR, area ratio; AI, asymmetry index; B, bilateral; c, cochlea; contra, contralateral; CEMD, central eye movement disorder; CISS, constructive interference in steady-state; d, definite; Diff-ER, ELS ratio of the side difference; DL, deep learning; ELH, endolymphatic hydrops; ELS, endolymphatic space; ER, ELS ratio; FLAIR, fluid-attenuated inversion recovery; Gd, gadolinium; GBCA, Gd-based contrast agents; Gd-DOTA, gadoteric acid (trade name: Dotarem®); Gd-Do3A, gadobutrol (trade name: Gadovist®); Gd-DTPA-BM, gadodiamide (trade name: Omniscan®); Gd-DTPA, gadopentetic acid (trade name: Magnevist®); Gd-HP-DO3A, gadoteridol (trade name: ProHance®); GLM, general linear model; GRAPPA, generalized auto-calibrating partially parallel acquisition; HC, healthy controls; HIT, head-impulse test; iMRI, delayed intravenous gadolinium-enhanced MRI of the inner ear; ipsi, ipsilateral; iv, intravenous; L, left; R, right; MD, Meniere's disease; minEn, minimum energy statistic; MMD, maximum mean discrepancy; MRI, magnetic resonance imaging; n, number; p, possible; ROI, region-of-interest; std, standard deviation; SNR, signal-to-noise ratio; SVV, subjective visual vertical; TFS, total fluid space; v, vestibulum; vHIT, video-oculography during head-impulse test; VM, vestibular migraine; VOG, video-oculography.
References
1. Strupp M, Brandt T, Dieterich M. Vertigo - Leitsymptom Schwindel. 3rd ed. Berlin Heidelberg: Springer-Verlag (2021). doi: 10.1007/978-3-662-61397-9
2. Pyykkö I, Zou J, Gürkov R, Naganawa S, Nakashima T. Imaging of temporal bone. Adv Otorhinolaryngol. (2019) 82:12–31. doi: 10.1159/000490268
3. van der Lubbe MFJA, Vaidyanathan A, Van Rompaey V, Postma AA, Bruintjes TD, Kimenai DM, et al. The “hype” of hydrops in classifying vestibular disorders: a narrative review. J Neurol. (2020) 267:197–211. doi: 10.1007/s00415-020-10278-8
4. Rauch SD, Merchant SN, Thedinger BA. Meniere's syndrome and endolymphatic hydrops. Double-blind temporal bone study. Ann Otol Rhinol Laryngol. (1989) 98:873–83. doi: 10.1177/000348948909801108
5. Nakashima T, Pyykkö I, Arroll MA, Casselbrant ML, Foster CA, Manzoor NF, et al. Meniere's disease. Nat Rev Dis Prim. (2016) 2:16028. doi: 10.1038/nrdp.2016.28
6. Nakada T, Yoshida T, Suga K, Kato M, Otake H, Kato K, et al. Endolymphatic space size in patients with vestibular migraine and Ménière's disease. J Neurol. (2014) 261:2079–84. doi: 10.1007/s00415-014-7458-9
7. Suga K, Kato M, Yoshida T, Nishio N, Nakada T, Sugiura S, et al. Changes in endolymphatic hydrops in patients with Ménière's disease treated conservatively for more than 1 year. Acta Oto-Laryngologica. (2015) 135:866–70. doi: 10.3109/00016489.2015.1015607
8. Attyé A, Eliezer M, Galloux A, Pietras J, Tropres I, Schmerber S, et al. Endolymphatic hydrops imaging: differential diagnosis in patients with Meniere disease symptoms. Diagn Interv Imaging. (2017) 98:699–706. doi: 10.1016/j.diii.2017.06.002
9. Kirsch V, Ertl-Wagner B, Berman A, Gerb J, Dieterich M, Becker-Bense S. High-resolution MRI of the inner ear enables syndrome differentiation and specific treatment of cerebellar downbeat nystagmus and secondary endolymphatic hydrops in a postoperative ELST patient. J Neurol. (2018) 265:48–50. doi: 10.1007/s00415-018-8858-z
10. Zwergal A, Kirsch V, Gerb J, Dlugaiczyk J, Becker-Bense S, Dieterich M. Neurootologie: Grenzfälle zwischen Ohr und Gehirn. Nervenarzt. (2018) 89:1106–14. doi: 10.1007/s00115-018-0598-x
11. Attyé A, Eliezer M, Medici M, Tropres I, Dumas G, Krainik A, et al. In vivo imaging of saccular hydrops in humans reflects sensorineural hearing loss rather than Meniere's disease symptoms. Eur Radiol. (2018) 28:2916–22. doi: 10.1007/s00330-017-5260-7
12. Eliezer M, Hautefort C, Van Nechel C, Duquesne U, Guichard J-P, Herman P, et al. Electrophysiological and inner ear MRI findings in patients with bilateral vestibulopathy. Eur Arch Otorhinolaryngol. (2020) 277:1305–14. doi: 10.1007/s00405-020-05829-8
13. Jerin C, Krause E, Ertl-Wagner B, Gürkov R. Clinical features of delayed endolymphatic hydrops and intralabyrinthine schwannoma : an imaging-confirmed comparative case series. English version. HNO. (2017) 65:41–5. doi: 10.1007/s00106-016-0199-6
14. Naganawa S, Nakamichi R, Ichikawa K, Kawamura M, Kawai H, Yoshida T, et al. MR imaging of endolymphatic hydrops: utility of iHYDROPS-Mi2 combined with deep learning reconstruction denoising. Magn Reson Med Sci. (2020) doi: 10.2463/mrms.mp.2020-0082. [Epub ahead of print].
15. Naganawa S, Kawai H, Taoka T, Sone M. Improved HYDROPS: Imaging of endolymphatic hydrops after intravenous administration of gadolinium. Magn Reson Med Sci. (2017) 16:357–61. doi: 10.2463/mrms.tn.2016-0126
16. Yamazaki M, Naganawa S, Tagaya M, Kawai H, Ikeda M, Sone M, et al. Comparison of contrast effect on the cochlear perilymph after intratympanic and intravenous gadolinium injection. Am J Neuroradiol. (2012) 33:773–8. doi: 10.3174/ajnr.A2821
17. Naganawa S, Satake H, Iwano S, Fukatsu H, Sone M, Nakashima T. Imaging endolymphatic hydrops at 3 tesla using 3D-FLAIR with intratympanic Gd-DTPA administration. Magn Reson Med Sci. (2008) 7:85–91. doi: 10.2463/mrms.7.85
18. Naganawa S, Yamazaki M, Kawai H, Bokura K, Sone M, Nakashima T. Visualization of endolymphatic hydrops in Ménière's disease with single-dose intravenous gadolinium-based contrast media using heavily T(2)-weighted 3D-FLAIR. Magn Reson Med Sci. (2010) 9:237–42. doi: 10.2463/mrms.9.237
19. Nakashima T, Naganawa S, Pyykko I, Gibson WPR, Sone M, Nakata S, et al. Grading of endolymphatic hydrops using magnetic resonance imaging. Acta Otolaryngol Suppl. (2009) 560:5–8. doi: 10.1080/00016480902729827
20. Baráth K, Schuknecht B, Naldi AM, Schrepfer T, Bockisch CJ, Hegemann SCA. Detection and grading of endolymphatic hydrops in Menière disease using MR imaging. Am J Neuroradiol. (2014) 35:1387–92. doi: 10.3174/ajnr.A3856
21. Bernaerts A, Vanspauwen R, Blaivie C, van Dinther J, Zarowski A, Wuyts FL, et al. The value of four stage vestibular hydrops grading and asymmetric perilymphatic enhancement in the diagnosis of Menière's disease on MRI. Neuroradiology. (2019) 61:421–9. doi: 10.1007/s00234-019-02155-7
22. Kirsch V, Becker-Bense S, Berman A, Kierig E, Ertl-Wagner B, Dieterich M. Transient endolymphatic hydrops after an attack of vestibular migraine: a longitudinal single case study. J Neurol. (2018) 265(Suppl. 1):51–3. doi: 10.1007/s00415-018-8870-3
23. Bernaerts A, De Foer B. Imaging of Ménière disease. Neuroimaging Clin North Am. (2019) 29:19–28. doi: 10.1016/j.nic.2018.09.002
24. Gürkov R, Flatz W, Louza J, Strupp M, Ertl-Wagner B, Krause E. In vivo visualized endolymphatic hydrops and inner ear functions in patients with electrocochleographically confirmed Ménière's disease. Otol Neurotol. (2012) 33:1040–5. doi: 10.1097/MAO.0b013e31825d9a95
25. Yang S, Zhu H, Zhu B, Wang H, Chen Z, Wu Y, et al. Correlations between the degree of endolymphatic hydrops and symptoms and audiological test results in patients with Menière's disease: a reevaluation. Otol Neurotol. (2018) 39:351–6. doi: 10.1097/MAO.0000000000001675
26. Homann G, Vieth V, Weiss D, Nikolaou K, Heindel W, Notohamiprodjo M, et al. Semi-quantitative vs. volumetric determination of endolymphatic space in Menière's disease using endolymphatic hydrops 3T-HR-MRI after intravenous gadolinium injection. PLoS ONE. (2015) 10:e0120357. doi: 10.1371/journal.pone.0120357
27. Naganawa S, Suzuki K, Nakamichi R, Bokura K, Yoshida T, Sone M, et al. Semi-quantification of endolymphatic size on MR imaging after intravenous injection of single-dose gadodiamide: comparison between two types of processing strategies. Magn Reson Med Sci. (2013) 12:261–9. doi: 10.2463/mrms.2013-0019
28. Ohashi T, Naganawa S, Takeuchi A, Katagiri T, Kuno K. Quantification of endolymphatic space volume after intravenous administration of a single dose of gadolinium-based contrast agent: 3D-real inversion recovery versus HYDROPS-Mi2. Magn Reson Med Sci. (2019) 19:119–24. doi: 10.2463/mrms.mp.2019-0013
29. Gürkov R, Berman A, Dietrich O, Flatz W, Jerin C, Krause E, et al. MR volumetric assessment of endolymphatic hydrops. Eur Radiol. (2015) 25:585–95. doi: 10.1007/s00330-014-3414-4
30. Inui H, Sakamoto T, Ito T, Kitahara T. Magnetic resonance-based volumetric measurement of the endolymphatic space in patients with Meniere's disease and other endolymphatic hydrops-related diseases. Auris Nasus Larynx. (2019) 46:493–7. doi: 10.1016/j.anl.2018.11.008
31. Gerb J, Ahmadi SA, Kierig E, Ertl-Wagner B, Dieterich M, Kirsch V. VOLT: a novel open-source pipeline for automatic segmentation of endolymphatic space in inner ear MRI. J Neurol. (2020) 267:185–96. doi: 10.1007/s00415-020-10062-8
32. Cho YS, Cho K, Park CJ, Chung MJ, Kim JH, Kim K, et al. Automated measurement of hydrops ratio from MRI in patients with Ménière's disease using CNN-based segmentation. Sci Rep. (2020) 10:7003. doi: 10.1038/s41598-020-63887-8
33. Lopez-Escamez JA, Carey J, Chung W-H, Goebel JA, Magnusson M, Mandalà M, et al. Diagnostic criteria for Menière's disease. J Vestib Res. (2015) 25:1–7. doi: 10.3233/VES-150549
34. Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. (1971) 9:97–113. doi: 10.1016/0028-3932(71)90067-4
35. Salmaso D, Longoni AM. Problems in the assessment of hand preference. Cortex. (1985) 21:533–49. doi: 10.1016/s0010-9452(58)80003-9
36. Dill T. Contraindications to magnetic resonance imaging. Heart. (2008) 94:943–8. doi: 10.1136/hrt.2007.125039
37. Matsumoto M, Nishimura T. Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator. ACM Trans Model Comput Simul. (1998) 8:3–30. doi: 10.1145/272991.272995
38. Dieterich M, Brandt T. Ocular torsion and tilt of subjective visual vertical are sensitive brainstem signs. Ann Neurol. (1993) 33:292–9. doi: 10.1002/ana.410330311
39. Halmagyi GM Curthoys IS. A clinical sign of canal paresis. Arch Neurol. (1988) 45:737–9. doi: 10.1001/archneur.1988.00520310043015
40. Schneider E, Villgrattner T, Vockeroth J, Bartl K, Kohlbecher S, Bardins S, et al. EyeSeeCam: an eye movement-driven head camera for the examination of natural visual exploration. Ann N Y Acad Sci. (2009) 1164:461–7. doi: 10.1111/j.1749-6632.2009.03858.x
41. Strupp M, Kim J-S, Murofushi T, Straumann D, Jen JC, Rosengren SM, et al. Bilateral vestibulopathy: diagnostic criteria consensus document of the classification committee of the Bárány society. J Vestib Res. (2017) 27:177–89. doi: 10.3233/VES-170619
42. Jongkees LB, Maas JP, Philipszoon AJ. Clinical nystagmography. A detailed study of electro-nystagmography in 341 patients with vertigo. Pract Otorhinolaryngol. (1962) 24:65–93.
43. Dlugaiczyk J, Habs M, Dieterich M. Vestibular evoked myogenic potentials in vestibular migraine and Menière's disease: cVEMPs make the difference. J Neurol. (2020) 267(Suppl. 1):169–80. doi: 10.1007/s00415-020-09902-4
44. Curthoys IS Dlugaiczyk J. Physiology, clinical evidence and diagnostic relevance of sound-induced and vibration-induced vestibular stimulation. Curr Opin Neurol. (2020) 33:126–35. doi: 10.1097/WCO.0000000000000770
45. Ertl M, Boegle R, Kirsch V, Dieterich M. On the impact of examiners on latencies and amplitudes in cervical and ocular vestibular-evoked myogenic potentials evaluated over a large sample (N = 1,038). Eur Arch Otorhinolaryngol. (2016) 273:317–23. doi: 10.1007/s00405-015-3510-3
46. Nahmani S, Vaussy A, Hautefort C, Guichard J-P, Guillonet A, Houdart E, et al. Comparison of enhancement of the vestibular perilymph between variable and constant flip angle-delayed 3D-FLAIR sequences in Menière disease. Am J Neuroradiol. (2020) doi: 10.3174/ajnr.A6483
47. Legland D, Arganda-Carreras I, Andrey P. MorphoLibJ: integrated library and plugins for mathematical morphology with ImageJ. Bioinformatics. (2016) 32:3532–4. doi: 10.1093/bioinformatics/btw413
48. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. (2012) 9:676–82. doi: 10.1038/nmeth.2019
49. Székely GJ, Rizzo ML. Energy statistics: a class of statistics based on distances. J Statist Plan Inference. (2013) 143:1249–72. doi: 10.1016/j.jspi.2013.03.018
50. Gretton A, Borgwardt K, Rasch M, Schölkopf B, Smola A. A kernel two-sample test. J Mach Learn Res. (2012) 13:723–73. doi: 10.5555/2503308.2188410
51. Milletari F, Navab N, Ahmadi S-A. V-Net: fully convolutional neural networks for volumetric medical image segmentation. arXiv:160604797. (2016). Available online at: http://arxiv.org/abs/1606.04797
52. Milletari F, Frei J, Aboulatta M, Vivar G, Ahmadi S-A. Cloud deployment of high-resolution medical image analysis with TOMAAT. IEEE J Biomed Health Inform. (2019) 23:969–77. doi: 10.1109/JBHI.2018.2885214
53. Fedorov A, Beichel R, Kalpathy-Cramer J, Finet J, Fillion-Robin J-C, Pujol S, et al. 3D slicer as an image computing platform for the quantitative imaging network. Magn Reson Imaging. (2012) 30:1323–41. doi: 10.1016/j.mri.2012.05.001
54. Nizamutdinov V. Neural Network and Fuzzy Logic Based Plugins for ImageJ. (2015). Available online at: https://github.com/astartes91/imagej-neural-fuzzy-plugins
55. Kendall MG, Smith BB. The problem of m rankings. Ann Math Statist. (1939) 10:275–87. doi: 10.1214/aoms/1177732186
56. Armitage P, Berry G, Matthews JNS. Modelling continuous data. In: Statistical Methods in Medical Research. 4th edn. Oxford: John Wiley & Sons, Ltd. (2002). p. 312–77. doi: 10.1002/9780470773666.ch11
57. Oh S-Y, Dieterich M, Lee BN, Boegle R, Kang J-J, Lee N-R, et al. Endolymphatic hydrops in patients with vestibular migraine and concurrent Meniere's disease. Front. Neurol. (2021) 12:594481. doi: 10.3389/fneur.2021.594481
58. Fukuoka H, Takumi Y, Tsukada K, Miyagawa M, Oguchi T, Ueda H, et al. Comparison of the diagnostic value of 3 T MRI after intratympanic injection of GBCA, electrocochleography, and the glycerol test in patients with Meniere's disease. Acta Oto-Laryngol. (2012) 132:141–5. doi: 10.3109/00016489.2011.635383
59. Naganawa S, Kawai H, Sone M, Nakashima T. Increased sensitivity to low concentration gadolinium contrast by optimized heavily T2-weighted 3D-FLAIR to visualize endolymphatic space. Magn Reson Med Sci. (2010) 9:73–80. doi: 10.2463/mrms.9.73
60. Naganawa S. The technical and clinical features of 3D-FLAIR in neuroimaging. Magn Reson Med Sci. (2015) 14:93–106. doi: 10.2463/mrms.2014-0132
61. Hirai T, Ando Y, Yamura M, Kitajima M, Hayashida Y, Korogi Y, et al. Transthyretin-related familial amyloid polyneuropathy: evaluation of CSF enhancement on serial T1-weighted and fluid- attenuated inversion recovery images following intravenous contrast administration. Am J Neuroradiol. (2005) 26:2043–8.
62. Naganawa S, Komada T, Fukatsu H, Ishigaki T, Takizawa O. Observation of contrast enhancement in the cochlear fluid space of healthy subjects using a 3D-FLAIR sequence at 3 Tesla. Eur Radiol. (2006) 16:733–7. doi: 10.1007/s00330-005-0046-8
63. Naganawa S, Suzuki K, Yamazaki M, Sakurai Y. Serial scans in healthy volunteers following intravenous administration of gadoteridol: time course of contrast enhancement in various cranial fluid spaces. Magn Reson Med Sci. (2014) 13:7–13. doi: 10.2463/mrms.2013-0056
64. Suzuki H, Teranishi M, Sone M, Yamazaki M, Naganawa S, Nakashima T. Contrast enhancement of the inner ear after intravenous administration of a standard or double dose of gadolinium contrast agents. Acta Otolaryngol. (2011) 131:1025–31. doi: 10.3109/00016489.2011.598552
65. Naganawa S, Koshikawa T, Nakamura T, Fukatsu H, Ishigaki T, Aoki I. High-resolution T1-weighted 3D real IR imaging of the temporal bone using triple-dose contrast material. Eur Radiol. (2003) 13:2650–8. doi: 10.1007/s00330-003-1922-8
66. Yamazaki M, Naganawa S, Kawai H, Sone M, Nakashima T. Gadolinium distribution in cochlear perilymph: differences between intratympanic and intravenous gadolinium injection. Neuroradiology. (2012) 54:1161–9. doi: 10.1007/s00234-012-1078-9
67. Iida T, Teranishi M, Yoshida T, Otake H, Sone M, Kato M, et al. Magnetic resonance imaging of the inner ear after both intratympanic and intravenous gadolinium injections. Acta Oto-Laryngol. (2013) 133:434–8. doi: 10.3109/00016489.2012.753640
68. Mamourian AC, Hoopes PJ, Lewis LD. Visualization of intravenously administered contrast material in the CSF on fluid-attenuated inversion-recovery MR images: an in vitro and animal-model investigation. Am J Neuroradiol. (2000) 21:105–111.
69. Naganawa S, Yamazaki M, Kawai H, Bokura K, Sone M, Nakashima T. Imaging of Ménière's disease by subtraction of MR cisternography from positive perilymph image. Magn Reson Med Sci. (2012) 11:303–9. doi: 10.2463/mrms.11.303
70. Naganawa S, Yamazaki M, Kawai H, Bokura K, Sone M, Nakashima T. Imaging of Ménière's disease after intravenous administration of single-dose gadodiamide: utility of multiplication of MR cisternography and HYDROPS image. Magn Reson Med Sci. (2013) 12:63–8. doi: 10.2463/mrms.2012-0027
71. Counter SA, Nikkhou-Aski S, Damberg P, Berglin CE, Laurell G. Ultra-high-field (9.4 T) MRI analysis of contrast agent transport across the blood-perilymph barrier and intrastrial fluid-blood barrier in the mouse inner ear. Otol Neurotol. (2017) 38:1052–9. doi: 10.1097/MAO.0000000000001458
72. Counter SA. MRI evidence of endolymphatic impermeability to the gadolinium molecule in the in vivo mouse inner ear at 9.4 Tesla. Open Neuroimag J. (2013) 7:27–31. doi: 10.2174/1874440001307010027
73. Naganawa S, Yamazaki M, Kawai H, Bokura K, Sone M, Nakashima T. Imaging of Ménière's disease after intravenous administration of single-dose gadodiamide: utility of subtraction images with different inversion time. Magn Reson Med Sci. (2012) 11:213–9. doi: 10.2463/mrms.11.213
74. Naganawa S, Kanou M, Ohashi T, Kuno K, Sone M. Simple estimation of the endolymphatic volume ratio after intravenous administration of a single-dose of gadolinium contrast. Magn Reson Med Sci. (2016) 15:379–85. doi: 10.2463/mrms.mp.2015-0175
75. Sepahdari AR, Ishiyama G, Vorasubin N, Peng KA, Linetsky M, Ishiyama A. Delayed intravenous contrast-enhanced 3D FLAIR MRI in Meniere's disease: correlation of quantitative measures of endolymphatic hydrops with hearing. Clinical Imaging. (2015) 39:26–31. doi: 10.1016/j.clinimag.2014.09.014
76. Kirsch V, Nejatbakhshesfahani F, Ahmadi S-A, Dieterich M, Ertl-Wagner B. A probabilistic atlas of the human inner ear's bony labyrinth enables reliable atlas-based segmentation of the total fluid space. J Neurol. (2019) 266(Suppl. 1):52–61. doi: 10.1007/s00415-019-09488-6
77. Ahmadi S-A, Raiser T, Ruehl RM, Flanagin VL, zu Eulenburg P. IE-map: a novel in-vivo atlas template of the human inner ear. Sci Rep. (2021) 11:3293. doi: 10.1038/s41598-021-82716-0
78. Ito T, Inui H, Miyasaka T, Shiozaki T, Matsuyama S, Yamanaka T, et al. Three-Dimensional magnetic resonance imaging reveals the relationship between the control of vertigo and decreases in endolymphatic hydrops after endolymphatic sac drainage with steroids for Meniere's disease. Front Neurol. (2019) 10:46. doi: 10.3389/fneur.2019.00046
79. Shi S, Guo P, Li W, Wang W. Clinical features and endolymphatic hydrops in patients with MRI evidence of hydrops. Ann Otol Rhinol Laryngol. (2019) 128:286–92. doi: 10.1177/0003489418819551
80. Shimono M, Teranishi M, Yoshida T, Kato M, Sano R, Otake H, et al. Endolymphatic hydrops revealed by magnetic resonance imaging in patients with acute low-tone sensorineural hearing loss. Otol Neurotol. (2013) 34:1241–6. doi: 10.1097/MAO.0b013e3182990e81
81. Attyé A, Eliezer M, Boudiaf N, Tropres I, Chechin D, Schmerber S, et al. MRI of endolymphatic hydrops in patients with Meniere's disease: a case-controlled study with a simplified classification based on saccular morphology. Eur Radiol. (2017) 27:3138–46. doi: 10.1007/s00330-016-4701-z
82. Cho YS, Ahn JM, Choi JE, Park HW, Kim Y-K, Kim H-J, et al. Usefulness of intravenous gadolinium inner ear MR imaging in diagnosis of Ménière's disease. Sci Rep. (2018) 8:17562. doi: 10.1038/s41598-018-35709-5
83. Hornibrook J, Flook E, Greig S, Babbage M, Goh T, Coates M, et al. MRI inner ear imaging and tone burst electrocochleography in the diagnosis of Ménière's disease. Otol Neurotol. (2015) 36:1109–14. doi: 10.1097/MAO.0000000000000782
84. He B, Zhang F, Zheng H, Sun X, Chen J, Chen J, et al. The correlation of a 2D volume-referencing endolymphatic-hydrops grading system with extra-tympanic electrocochleography in patients with definite Ménière's disease. Front Neurol. (2021) 11:595038. doi: 10.3389/fneur.2020.595038
85. Gürkov R, Kantner C, Strupp M, Flatz W, Krause E, Ertl-Wagner B. Endolymphatic hydrops in patients with vestibular migraine and auditory symptoms. Eur Arch Otorhinolaryngol. (2014) 271:2661–7. doi: 10.1007/s00405-013-2751-2
86. Jerin C, Floerke S, Maxwell R, Gürkov R. Relationship between the extent of endolymphatic hydrops and the severity and fluctuation of audiovestibular symptoms in patients with Menière's disease and MRI evidence of hydrops. Otol Neurotol. (2018) 39:e123–30. doi: 10.1097/MAO.0000000000001681
87. Eliezer M, Poillon G, Maquet C, Gillibert A, Horion J, Marie J-P, et al. Sensorineural hearing loss in patients with vestibular schwannoma correlates with the presence of utricular hydrops as diagnosed on heavily T2-weighted MRI. Diagn Interv Imaging. (2019) 100:259–68. doi: 10.1016/j.diii.2019.01.006
88. Katayama N, Yamamoto M, Teranishi M, Naganawa S, Nakata S, Sone M, et al. Relationship between endolymphatic hydrops and vestibular-evoked myogenic potential. Acta Oto-Laryngol. (2010) 130:917–23. doi: 10.3109/00016480903573187
89. Kato M, Teranishi M, Katayama N, Sone M, Naganawa S, Nakashima T. Association between endolymphatic hydrops as revealed by magnetic resonance imaging and caloric response. Otol Neurotol. (2011) 32:1480–5. doi: 10.1097/MAO.0b013e318235568d
90. Zhang W, Hui L, Zhang B, Ren L, Zhu J, Wang F, et al. The correlation between endolymphatic hydrops and clinical features of Meniere disease. Laryngoscope. (2020) 131:E144–50. doi: 10.1002/lary.28576
91. Shi S, Guo P, Wang W. Magnetic resonance imaging of Ménière's disease after intravenous administration of gadolinium. Ann Otol Rhinol Laryngol. (2018) 127:777–82. doi: 10.1177/0003489418794699
92. Inui H, Sakamoto T, Ito T, Kitahara T. Volumetric measurements of the inner ear in patients with Meniere's disease using three-dimensional magnetic resonance imaging. Acta Oto-Laryngol. (2016) 136:888–93. doi: 10.3109/00016489.2016.1168940
93. Naganawa S, Ohashi T, Kanou M, Kuno K, Sone M, Ikeda M. Volume quantification of endolymph after intravenous administration of a single dose of gadolinium contrast agent: comparison of 18- versus 8-minute imaging protocols. Magn Reson Med Sci. (2015) 14:257–62. doi: 10.2463/mrms.2014-0118
94. Ito T, Inui H, Miyasaka T, Shiozaki T, Hasukawa A, Yamanaka T, et al. Endolymphatic volume in patients with Meniere's disease and healthy controls: three-dimensional analysis with magnetic resonance imaging. Laryngoscope Investig Otolaryngol. (2019) 4:653–8. doi: 10.1002/lio2.313
95. Ito T, Kitahara T, Inui H, Miyasaka T, Kichikawa K, Ota I, et al. Endolymphatic space size in patients with Meniere's disease and healthy controls. Acta Oto-Laryngol. (2016) 136:879–82. doi: 10.3109/00016489.2016.1169556
96. Naganawa S, Satake H, Kawamura M, Fukatsu H, Sone M, Nakashima T. Separate visualization of endolymphatic space, perilymphatic space and bone by a single pulse sequence; 3D-inversion recovery imaging utilizing real reconstruction after intratympanic Gd-DTPA administration at 3 Tesla. Eur Radiol. (2008) 18:920–4. doi: 10.1007/s00330-008-0854-8
97. Adin ME, Kleinberg L, Vaidya D, Zan E, Mirbagheri S, Yousem DM. Hyperintense dentate nuclei on T1-weighted MRI: relation to repeat gadolinium administration. Am J Neuroradiol. (2015) 36:1859–65. doi: 10.3174/ajnr.A4378
98. Barbieri S, Schroeder C, Froehlich JM, Pasch A, Thoeny HC. High signal intensity in dentate nucleus and globus pallidus on unenhanced T1-weighted MR images in three patients with impaired renal function and vascular calcification. Contrast Media Mol Imaging. (2016) 11:245–50. doi: 10.1002/cmmi.1683
99. Errante Y, Cirimele V, Mallio CA, Di Lazzaro V, Zobel BB, Quattrocchi CC. Progressive increase of T1 signal intensity of the dentate nucleus on unenhanced magnetic resonance images is associated with cumulative doses of intravenously administered gadodiamide in patients with normal renal function, suggesting dechelation. Investigat Radiol. (2014) 49:685–90. doi: 10.1097/RLI.0000000000000072
Keywords: endolymphatic hydrops, endolymphatic space, inner ear imaging, gadolinium based contrast agent, intravenous, convolutional neural network, deep learning, volumetric local thresholding
Citation: Boegle R, Gerb J, Kierig E, Becker-Bense S, Ertl-Wagner B, Dieterich M and Kirsch V (2021) Intravenous Delayed Gadolinium-Enhanced MR Imaging of the Endolymphatic Space: A Methodological Comparative Study. Front. Neurol. 12:647296. doi: 10.3389/fneur.2021.647296
Received: 29 December 2020; Accepted: 24 February 2021;
Published: 22 April 2021.
Edited by:
Stefan K. Plontke, Martin Luther University of Halle-Wittenberg, GermanyReviewed by:
Alexander Gussew, Universitätsklinikum Halle, GermanyArne Ernst, Unfallkrankenhaus Berlin, Germany
Copyright © 2021 Boegle, Gerb, Kierig, Becker-Bense, Ertl-Wagner, Dieterich and Kirsch. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Valerie Kirsch, dmFsZXJpZS5raXJzY2gmI3gwMDA0MDttZWQubG11LmRl