Using multidimensional topological data analysis to identify traits of hip osteoarthritis

Development of a Systems Medicine Approach to Spinal Cord Injury

Traumatic spinal cord injury (SCI) causes a sudden onset multi-system disease, permanently altering homeostasis with multiple complications. Consequences include aberrant neuronal circuits, multiple organ system dysfunctions, and chronic phenotypes such as neuropathic pain and metabolic syndrome. Reductionist approaches are used to classify SCI patients based on residual neurological function. Still, recovery varies due to interacting variables, including individual biology, comorbidities, complications, therapeutic side effects, and socioeconomic influences for which data integration methods are lacking. Infections, pressure sores, and heterotopic ossification are known recovery modifiers. However, the molecular pathobiology of the disease-modifying factors altering the neurological recovery-chronic syndrome trajectory is mainly unknown, with significant data gaps between intensive early treatment and chronic phases. Changes in organ function such as gut dysbiosis, adrenal dysregulation, fatty liver, muscle loss, and autonomic dysregulation disrupt homeostasis, generating progression-driving allostatic load. Interactions between interdependent systems produce emergent effects, such as resilience, that preclude single mechanism interpretations. Due to many interacting variables in individuals, substantiating the effects of treatments to improve neurological outcomes is difficult. Acute injury outcome predictors, including blood and cerebrospinal fluid biomarkers, neuroimaging signal changes, and autonomic system abnormalities, often do not predict chronic SCI syndrome phenotypes. In systems medicine, network analysis of bioinformatics data is used to derive molecular control modules. To better understand the evolution from acute SCI to chronic SCI multi-system states, we propose a topological phenotype framework integrating bioinformatics, physiological data, and allostatic load tested against accepted established recovery metrics. This form of correlational phenotyping may reveal critical nodal points for intervention to improve recovery trajectories. This study examines the limitations of current classifications of SCI and how these can evolve through systems medicine.

The evolution of Big Data in neuroscience and neurology

Neurological diseases are on the rise worldwide, leading to increased healthcare costs and diminished quality of life in patients. In recent years, Big Data has started to transform the fields of Neuroscience and Neurology. Scientists and clinicians are collaborating in global alliances, combining diverse datasets on a massive scale, and solving complex computational problems that demand the utilization of increasingly powerful computational resources. This Big Data revolution is opening new avenues for developing innovative treatments for neurological diseases. Our paper surveys Big Data's impact on neurological patient care, as exemplified through work done in a comprehensive selection of areas, including Connectomics, Alzheimer's Disease, Stroke, Depression, Parkinson's Disease, Pain, and Addiction (e.g., Opioid Use Disorder). We present an overview of research and the methodologies utilizing Big Data in each area, as well as their current limitations and technical challenges. Despite the potential benefits, the full potential of Big Data in these fields currently remains unrealized. We close with recommendations for future research aimed at optimizing the use of Big Data in Neuroscience and Neurology for improved patient outcomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40537-023-00751-2.

A progression analysis of motor features in Parkinson's disease based on the mapper algorithm

Background

Parkinson's disease (PD) is a neurodegenerative disease with a broad spectrum of motor and non-motor symptoms. The great heterogeneity of clinical symptoms, biomarkers, and neuroimaging and lack of reliable progression markers present a significant challenge in predicting disease progression and prognoses.

Methods

We propose a new approach to disease progression analysis based on the mapper algorithm, a tool from topological data analysis. In this paper, we apply this method to the data from the Parkinson's Progression Markers Initiative (PPMI). We then construct a Markov chain on the mapper output graphs.

Results

The resulting progression model yields a quantitative comparison of patients' disease progression under different usage of medications. We also obtain an algorithm to predict patients' UPDRS III scores.

Conclusions

By using mapper algorithm and routinely gathered clinical assessments, we developed a new dynamic models to predict the following year's motor progression in the early stage of PD. The use of this model can predict motor evaluations at the individual level, assisting clinicians to adjust intervention strategy for each patient and identifying at-risk patients for future disease-modifying therapy clinical trials.

Topological data analysis in biomedicine: A review

Significant technological advances made in recent years have shepherded a dramatic increase in utilization of digital technologies for biomedicine- everything from the widespread use of electronic health records to improved medical imaging capabilities and the rising ubiquity of genomic sequencing contribute to a "digitization" of biomedical research and clinical care. With this shift toward computerized tools comes a dramatic increase in the amount of available data, and current tools for data analysis capable of extracting meaningful knowledge from this wealth of information have yet to catch up. This article seeks to provide an overview of emerging mathematical methods with the potential to improve the abilities of clinicians and researchers to analyze biomedical data, but may be hindered from doing so by a lack of conceptual accessibility and awareness in the life sciences research community. In particular, we focus on topological data analysis (TDA), a set of methods grounded in the mathematical field of algebraic topology that seeks to describe and harness features related to the "shape" of data. We aim to make such techniques more approachable to non-mathematicians by providing a conceptual discussion of their theoretical foundations followed by a survey of their published applications to scientific research. Finally, we discuss the limitations of these methods and suggest potential avenues for future work integrating mathematical tools into clinical care and biomedical informatics.

Use of machine learning in osteoarthritis research: a systematic literature review

Objective

The aim of this systematic literature review was to provide a comprehensive and exhaustive overview of the use of machine learning (ML) in the clinical care of osteoarthritis (OA).

Methods

A systematic literature review was performed in July 2021 using MEDLINE PubMed with key words and MeSH terms. For each selected article, the number of patients, ML algorithms used, type of data analysed, validation methods and data availability were collected.

Results

From 1148 screened articles, 46 were selected and analysed; most were published after 2017. Twelve articles were related to diagnosis, 7 to prediction, 4 to phenotyping, 12 to severity and 11 to progression. The number of patients included ranged from 18 to 5749. Overall, 35% of the articles described the use of deep learning And 74% imaging analyses. A total of 85% of the articles involved knee OA and 15% hip OA. No study investigated hand OA. Most of the studies involved the same cohort, with data from the OA initiative described in 46% of the articles and the MOST and Cohort Hip and Cohort Knee cohorts in 11% and 7%. Data and source codes were described as publicly available respectively in 54% and 22% of the articles. External validation was provided in only 7% of the articles.

Conclusion

This review proposes an up-to-date overview of ML approaches used in clinical OA research and will help to enhance its application in this field.

[18 F]-Sodium Fluoride PET/MR Imaging for Bone-Cartilage Interactions in Hip Osteoarthritis: A Feasibility Study

This study characterized the distribution of [¹⁸ F]-sodium fluoride (NaF) uptake and blood flow in the femur and acetabulum in hip osteoarthritis (OA) patients to find associations between bone remodeling and cartilage composition in the presence of morphological abnormalities using simultaneous positron emission tomography and magnetic resonance imaging (PET/MR), quantitative magnetic resonance imaging (MRI) and femur shape modeling. Ten patients underwent a [¹⁸ F]-NaF PET/MR dynamic scan of the hip simultaneously with: (i) fast spin-echo CUBE for morphology grading and (ii) T_1ρ /T₂ magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots for cartilage, bone segmentation, bone shape modeling, and T_1ρ /T₂ quantification. The standardized uptake values (SUVs) and Patlak kinetic parameter (K_pat ) were calculated for each patient as PET outcomes, using an automated post-processing pipeline. Shape modeling was performed to extract the variations in bone shapes in the patients. Pearson's correlation coefficients were used to study the associations between bone shapes, PET outcomes, and patient reported pain. Direct associations between quantitative MR and PET evidence of bone remodeling were established in the acetabulum and femur. Associations of shaft thickness with SUV in the femur (p = 0.07) and K_pat in the acetabulum (p = 0.02), cam deformity with acetabular score (p = 0.09), osteophytic growth on the femur head with K_pat (p = 0.01) were observed. Pain had increased correlations with SUV in the acetabulum (p = 0.14) and femur (p = 0.09) when shaft thickness was accounted for. This study demonstrated the ability of [¹⁸ F]-NaF PET-MRI, 3D shape modeling, and quantitative MRI to investigate cartilage-bone interactions and bone shape features in hip OA, providing potential investigative tools to diagnose OA. © 2019 The Authors. Journal of Orthopaedic Research^® published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society J Orthop Res 37:2671-2680, 2019.

Generating dynamical neuroimaging spatiotemporal representations (DyNeuSR) using topological data analysis

In this article, we present an open source neuroinformatics platform for exploring, analyzing, and validating distilled graphical representations of high-dimensional neuroimaging data extracted using topological data analysis (TDA). TDA techniques like Mapper have been recently applied to examine the brain's dynamical organization during ongoing cognition without averaging data in space, in time, or across participants at the outset. Such TDA-based approaches mark an important deviation from standard neuroimaging analyses by distilling complex high-dimensional neuroimaging data into simple-yet neurophysiologically valid and behaviorally relevant-representations that can be interactively explored at the single-participant level. To facilitate wider use of such techniques within neuroimaging and general neuroscience communities, our work provides several tools for visualizing, interacting with, and grounding TDA-generated graphical representations in neurophysiology. Through Python-based Jupyter notebooks and open datasets, we provide a platform to assess and visualize different intermittent stages of Mapper and examine the influence of Mapper parameters on the generated representations. We hope this platform could enable researchers and clinicians alike to explore topological representations of neuroimaging data and generate biological insights underlying complex mental disorders.

Diagnosing osteoarthritis from T2 maps using deep learning: an analysis of the entire Osteoarthritis Initiative baseline cohort

OBJECTIVE

We aim to study to what extent conventional and deep-learning-based T2 relaxometry patterns are able to distinguish between knees with and without radiographic osteoarthritis (OA).

METHODS

T2 relaxation time maps were analyzed for 4,384 subjects from the baseline Osteoarthritis Initiative (OAI) Dataset. Voxel Based Relaxometry (VBR) was used for automatic quantification and voxel-based analysis of the differences in T2 between subjects with and without radiographic OA. A Densely Connected Convolutional Neural Network (DenseNet) was trained to diagnose OA from T2 data. For comparison, more classical feature extraction techniques and shallow classifiers were used to benchmark the performance of our algorithm's results. Deep and shallow models were evaluated with and without the inclusion of risk factors. Sensitivity and Specificity values and McNemar test were used to compare the performance of the different classifiers.

RESULTS

The best shallow model was obtained when the first ten Principal Components, demographics and pain score were included as features (AUC = 77.77%, Sensitivity = 67.01%, Specificity = 71.79%). In comparison, DenseNet trained on raw T2 data obtained AUC = 83.44%, Sensitivity = 76.99%, Specificity = 77.94%. McNemar test on two misclassified proportions form the shallow and deep model showed that the boost in performance was statistically significant (McNemar's chi-squared = 10.33, degree of freedom (DF) = 1, P-value = 0.0013).

CONCLUSION

In this study, we presented a Magnetic Resonance Imaging (MRI)-based data-driven platform using T2 measurements to characterize radiographic OA. Our results showed that feature learning from T2 maps has potential in uncovering information that can potentially better diagnose OA than simple averages or linear patterns decomposition.

Translation of morphological and functional musculoskeletal imaging

In an effort to develop quantitative biomarkers for degenerative joint disease and fill the void that exists for diagnosing, monitoring, and assessing the extent of whole joint degeneration, the past decade has been marked by a greatly increased role of noninvasive imaging. This coupled with recent advances in image processing and deep learning opens new possibilities for promising quantitative techniques. The clinical translation of quantitative imaging was previously hampered by tedious non-scalable and subjective image analysis. Osteoarthritis (OA) diagnosis using X-rays can be automated by the use of deep learning models and pilot studies showed feasibility of using similar techniques to reliably segment multiple musculoskeletal tissues and detect and stage the severity of morphological abnormalities in magnetic resonance imaging (MRI). Automation and more advanced feature extraction techniques have applications on larger more heterogeneous samples. Analyses based on voxel based relaxometry have shown local patterns in relaxation time elevations and local correlations with outcome variables. Bone cartilage interactions are also enhanced by the analysis of three-dimensional bone morphology and the potential for the assessment of metabolic activity with simultaneous Positron Emission Tomography (PET)/MR systems. Novel techniques in image processing and deep learning are augmenting imaging to be a source of quantitative and reliable data and new multidimensional analytics allow us to exploit the interactions of data from various sources. In this review, we aim to summarize recent advances in quantitative imaging, the application of image processing and deep learning techniques to study knee and hip OA. ©2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res XX:XX-XX, 2018.