Estimation and testing of Wilcoxon–Mann–Whitney effects in factorial clustered data designs

Clustered data arise frequently in many practical applications whenever units are repeatedly observed under a certain condition. One typical example for clustered data are animal experiments, where several animals share the same cage and should not be assumed to be completely independent. Standard methods for the analysis of such data are Linear Mixed Models and Generalized Estimating Equations—however, checking their assumptions is not easy, especially in scenarios with small sample sizes, highly skewed, count, and ordinal or binary data. In such situations, Wilcoxon–Mann–Whitney type effects are suitable alternatives to mean-based or other distributional approaches. Hence, no specific data distribution, symmetric or asymmetric, is required. Within this work, we will present different estimation techniques of such effects in clustered factorial designs and discuss quadratic- and multiple contrast type-testing procedures for hypotheses formulated in terms of Wilcoxon–Mann–Whitney effects. Additionally, the framework allows for the occurrence of missing data: estimation and testing hypotheses are based on all-available data instead of complete-cases. An extensive simulation study investigates the precision of the estimators and the behavior of the test procedures in terms of their type-I error control. One real world dataset exemplifies the applicability of the newly proposed procedures

Validation of computer-adaptive contrast sensitivity as a tool to assess visual impairment in multiple sclerosis patients

BACKGROUND: Impairment of visual function is one of the major symptoms of people with multiple sclerosis (pwMS). A multitude of disease effects including inflammation and neurodegeneration lead to structural impairment in the visual system. However, the gold standard of disability quantification, the expanded disability status scale (EDSS), relies on visual assessment charts. A more comprehensive assessment of visual function is the full contrast sensitivity function (CSF), but most tools are time consuming and not feasible in clinical routine. The quantitative CSF (qCSF) test is a computerized test to assess the full CSF. We have already shown a better correlation with visual quality of life (QoL) than for classical high and low contrast charts in multiple sclerosis (MS). OBJECTIVE: To study the precision, test duration, and repeatability of the qCSF in pwMS. In order to evaluate the discrimination ability, we compared the data of pwMS to healthy controls. METHODS: We recruited two independent cohorts of MS patients. Within the precision cohort (n = 54), we analyzed the benefit of running 50 instead of 25 qCSF trials. The repeatability cohort (n = 44) was assessed by high contrast vision charts and qCSF assessments twice and we computed repeatability metrics. For the discrimination ability we used the data from all pwMS without any previous optic neuritis and compared the area under the log CSF (AULCSF) to an age-matched healthy control data set. RESULTS: We identified 25 trials of the qCSF algorithm as a sufficient amount for a precise estimate of the CSF. The median test duration for one eye was 185 s (range 129-373 s). The AULCSF had better test-retest repeatability (Mean Average Precision, MAP) than visual acuity measured by standard high contrast visual acuity charts or CSF acuity measured with the qCSF (0.18 vs. 0.11 and 0.17, respectively). Even better repeatability (MAP = 0.19) was demonstrated by a CSF-derived feature that was inspired by low-contrast acuity charts, i.e., the highest spatial frequency at 25% contrast. When compared to healthy controls, the MS patients showed reduced CSF (average AULCSF 1.21 vs. 1.42, p < 0.01). CONCLUSION: High precision, usability, repeatability, and discrimination support the qCSF as a tool to assess contrast vision in pwMS

Fingolimod after a first unilateral episode of acute optic neuritis (MOVING) - preliminary results from a randomized, rater-blind, active-controlled, phase 2 trial

BACKGROUND: Neuroprotection and promotion of remyelination represent important therapeutic gaps in multiple sclerosis (MS). Acute optic neuritis (ON) is a frequent MS manifestation. Based on the presence and properties of sphingosine-1-phosphate receptors (S1PR) on astrocytes and oligodendrocytes, we hypothesized that remyelination can be enhanced by treatment with fingolimod, a S1PR modulator currently licensed for relapsing-remitting MS. METHODS: MOVING was an investigator-driven, rater-blind, randomized clinical trial. Patients with acute unilateral ON, occurring as a clinically isolated syndrome or MS relapse, were randomized to 6 months of treatment with 0.5 mg oral fingolimod or subcutaneous IFN-β 1b 250 μg every other day. The change in multifocal visual evoked potential (mfVEP) latency of the qualifying eye was examined as the primary (month 6 vs. baseline) and secondary (months 3, 6 and 12 vs. baseline) outcome. In addition, full field visual evoked potentials, visual acuity, optical coherence tomography as well as clinical relapses and measures of disability, cerebral MRI, and self-reported visual quality of life were obtained for follow-up. The study was halted due to insufficient recruitment (n = 15), and available results are reported. RESULTS: Per protocol analysis of the primary endpoint revealed a significantly larger reduction of mfVEP latency at 6 months compared to baseline with fingolimod treatment (n = 5; median decrease, 15.7 ms) than with IFN-β 1b treatment (n = 4; median increase, 8.15 ms) (p <  0.001 for interaction). Statistical significance was maintained in the secondary endpoint analysis. Descriptive results are reported for other endpoints. CONCLUSION: Preliminary results of the MOVING trial argue in support of a beneficial effect of fingolimod on optic nerve remyelination when compared to IFN-β treatment. Interpretation is limited by the small number of complete observations, an unexpected deterioration of the control group and a difference in baseline mfVEP latencies. The findings need to be confirmed in larger studies. TRIAL REGISTRATION: The trial was registered as EUDRA-CT 2011-004787-30 on October 26, 2012 and as NCT01647880 on July 24, 2012

Altered fovea in AQP4-IgG-seropositive neuromyelitis optica spectrum disorders

OBJECTIVE: To investigate disease-specific foveal shape changes in patients with neuromyelitis optica spectrum disorders (NMOSDs) using foveal morphometry. METHODS: This cross-sectional study included macular spectral domain optical coherence tomography scans of 52 eyes from 28 patients with aquaporin-4 immunoglobulin G (AQP4-IgG)-seropositive NMOSD, 116 eyes from 60 patients with MS, and 123 eyes from 62 healthy controls (HCs), retrospectively, and an independent confirmatory cohort comprised 33/33 patients with NMOSD/MS. The fovea was characterized using 3D foveal morphometry. We included peripapillary retinal nerve fiber layer (pRNFL) thickness and combined macular ganglion cell and inner plexiform layer (GCIPL) volume to account for optic neuritis (ON)-related neuroaxonal damage. RESULTS: Group comparison showed significant differences compared with HC in the majority of foveal shape parameters in NMOSD, but not MS. Pit flat disk area, average pit flat disk diameter, inner rim volume, and major slope disk length, as selected parameters, showed differences between NMOSD and MS (p value = 0.017, 0.002, 0.005, and 0.033, respectively). This effect was independent of ON. Area under the curve was between 0.7 and 0.8 (receiver operating characteristic curve) for discriminating between NMOSD and MS. Pit flat disk area and average pit flat disk diameter changes independent of ON were confirmed in an independent cohort. CONCLUSIONS: Foveal morphometry reveals a wider and flatter fovea in NMOSD in comparison to MS and HC. Comparison to MS and accounting for ON suggest this effect to be at least in part independent of ON. This supports a primary retinopathy in AQP4-IgG–seropositive NMOSD

Prior optic neuritis detection on peripapillary ring scans using deep learning

BACKGROUND: The diagnosis of multiple sclerosis (MS) requires demyelinating events that are disseminated in time and space. Peripapillary retinal nerve fiber layer (pRNFL) thickness as measured by optical coherence tomography (OCT) distinguishes eyes with a prior history of acute optic neuritis (ON) and may provide evidence to support a demyelinating attack. OBJECTIVE: To investigate whether a deep learning (DL)-based network can distinguish between eyes with prior ON and healthy control (HC) eyes using peripapillary ring scans. METHODS: We included 1033 OCT scans from 415 healthy eyes (213 HC subjects) and 510 peripapillary ring scans from 164 eyes with prior acute ON (140 patients with MS). Data were split into 70% training, 15% validation, and 15% test data. We included 102 OCT scans from 80 healthy eyes (40 HC) and 61 scans from 40 ON eyes (31 MS patients) from an independent second center. Receiver operating characteristic curve analyses with area under the curve (AUC) were used to investigate performance. RESULTS: We used a dilated residual convolutional neural network for the classification. The final network had an accuracy of 0.85 and an AUC of 0.86, whereas pRNFL only had an AUC of 0.77 in recognizing ON eyes. Using data from a second center, the network achieved an accuracy of 0.77 and an AUC of 0.90 compared to pRNFL, which had an AUC of 0.84. INTERPRETATION: DL-based disease classification of prior ON is feasible and has the potential to outperform thickness-based classification of eyes with and without history of prior ON

CD4+ T cell mitochondrial genotype in Multiple Sclerosis: a cross-sectional and longitudinal analysis

Multiple Sclerosis (MS) is a chronic autoimmune demyelinating disease of the central nervous system (CNS), with a largely unknown etiology, where mitochondrial dysfunction likely contributes to neuroaxonal loss and brain atrophy. Mirroring the CNS, peripheral immune cells from patients with MS, particularly CD4+ T cells, show inappropriate mitochondrial phenotypes and/or oxidative phosphorylation (OxPhos) insufficiency, with a still unknown contribution of mitochondrial DNA (mtDNA). We hypothesized that mitochondrial genotype in CD4+ T cells might influence MS disease activity and progression. Thus, we performed a retrospective cross-sectional and longitudinal study on patients with a recent diagnosis of either Clinically Isolated Syndrome (CIS) or Relapsing–Remitting MS (RRMS) at two timepoints: 6 months (VIS1) and 36 months (VIS2) after disease onset. Our primary outcomes were the differences in mtDNA extracted from CD4+ T cells between: (I) patients with CIS/RRMS (PwMS) at VIS1 and age- and sex-matched healthy controls (HC), in the cross-sectional analysis, and (II) different diagnostic evolutions in PwMS from VIS1 to VIS2, in the longitudinal analysis. We successfully performed mtDNA whole genome sequencing (mean coverage: 2055.77 reads/base pair) in 183 samples (61 triplets). Nonetheless, mitochondrial genotype was not associated with a diagnosis of CIS/RRMS, nor with longitudinal diagnostic evolution

Interactions of optic radiation lesions with retinal and brain atrophy in early multiple sclerosis

OBJECTIVE: Retrograde trans-synaptic neuroaxonal degeneration is considered a key pathological factor of subclinical retinal neuroaxonal damage in multiple sclerosis (MS). We aim to evaluate the longitudinal association of optic radiation (OR) lesion activity with retinal neuroaxonal damage and its role in correlations between retinal and brain atrophy in people with clinically isolated syndrome and early MS (pweMS). METHODS: Eighty-five pweMS were retrospectively screened from a prospective cohort (Berlin CIS cohort). Participants underwent 3T magnetic resonance imaging (MRI) for OR lesion volume and brain atrophy measurements and optical coherence tomography (OCT) for retinal layer thickness measurements. All pweMS were followed with serial OCT and MRI over a median follow-up of 2.9 (interquartile range: 2.6-3.4) years. Eyes with a history of optic neuritis prior to study enrollment were excluded. Linear mixed models were used to analyze the association of retinal layer thinning with changes in OR lesion volume and brain atrophy. RESULTS: Macular ganglion cell-inner plexiform layer (GCIPL) thinning was more pronounced in pweMS with OR lesion volume increase during follow-up compared to those without (Difference: -0.82 μm [95% CI:-1.49 to -0.15], p = 0.018). Furthermore, GCIPL thinning correlated with both OR lesion volume increase (β [95% CI] = -0.27 [-0.50 to -0.03], p = 0.028) and brain atrophy (β [95% CI] = 0.47 [0.25 to 0.70], p < 0.001). Correlations of GCIPL changes with brain atrophy did not differ between pweMS with or without OR lesion increase (η(2)(p) = 5.92e(-7), p = 0.762). INTERPRETATION: Faster GCIPL thinning rate is associated with increased OR lesion load. Our results support the value of GCIPL as a sensitive biomarker reflecting both posterior visual pathway pathology and global brain neurodegeneration

Impaired motion perception is associated with functional and structural visual pathway damage in multiple sclerosis and neuromyelitis optica spectrum disorders

BACKGROUND: Decreased motion perception has been suggested as a marker for visual pathway demyelination in optic neuritis (ON) and/or multiple sclerosis (MS). OBJECTIVES: To examine the influence of neuro-axonal damage on motion perception in MS and neuromyelitis optica spectrum disorders (NMOSD). METHODS: We analysed motion perception with numbers-from-motion (NFM), visual acuity, (multifocal (mf)) VEP, optical coherence tomography in patients with MS (n = 38, confirmatory cohort n = 43), NMOSD (n = 13) and healthy controls (n = 33). RESULTS: NFM was lower compared with controls in MS (B = −12.37, p < 0.001) and NMOSD (B = −34.5, p < 0.001). NFM was lower in ON than in non-ON eyes (B = −30.95, p = 0.041) in NMOSD, but not MS. In MS and NMOSD, lower NFM was associated with worse visual acuity (B = −139.4, p < 0.001/B = −77.2, p < 0.001) and low contrast letter acuity (B = 0.99, p = 0.002/B = 1.6, p < 0.001), thinner peripapillary retinal nerve fibre layer (B = 1.0, p < 0.001/ B = 0.92, p = 0.016) and ganglion cell/inner plexiform layer (B = 64.8, p < 0.001/B = 79.5, p = 0.006), but not with VEP P100 latencies. In the confirmatory MS cohort, lower NFM was associated with thinner retinal nerve fibre layer (B = 1.351, p < 0.001) and increased mfVEP P100 latencies (B = −1.159, p < 0.001). CONCLUSIONS: Structural neuro-axonal visual pathway damage is an important driver of motion perception impairment in MS and NMOSD

Structure-function correlates of vision loss in neuromyelitis optica spectrum disorders

Optic neuritis (ON) in neuromyelitis optica spectrum disorders (NMOSD) regularly leads to more profound vision loss compared to multiple sclerosis (MS) and myelin-oligodendrocyte-glycoprotein-antibody associated disease (MOGAD). Here we investigate ON-related vision loss in NMOSD compared to MS and MOGAD in order to identify neuroaxonal and retinal contributors to visual dysfunction. In this retrospective study we included patients with aquaporin-4-antibody seropositive NMOSD (n = 28), MOGAD (n = 14), MS (n = 29) and controls (n = 14). We assessed optic nerve damage and fovea morphometry by optical coherence tomography. Visual function was assessed as high (HCVA) and low contrast visual acuity (LCVA), and visual fields' mean deviation (MD). In all diseases, lower visual function was associated with peripapillary retinal nerve fiber layer (pRNFL) and ganglion cell and inner plexiform layer (GCIP) thinning following a broken stick model, with pRNFL and GCIP cutoff point at ca. 60 µm. HCVA loss per µm pRNFL and GCIP thinning was stronger in NMOSD compared with MOGAD. Foveal inner rim volume contributed to MD and LCVA in NMOSD eyes, only. Together these data supports that visual dysfunction in NMOSD is associated with neuroaxonal damage beyond the effect seen in MS and MOGAD. A primary retinopathy, respectively Müller cell pathology, may contribute to this effect