Artificial intelligence of imaging and clinical neurological data for predictive, preventive and personalized (P3) medicine for Parkinson Disease: The NeuroArtP3 protocol for a multi-center research study.

BACKGROUND: The burden of Parkinson Disease (PD) represents a key public health issue and it is essential to develop innovative and cost-effective approaches to promote sustainable diagnostic and therapeutic interventions. In this perspective the adoption of a P3 (predictive, preventive and personalized) medicine approach seems to be pivotal. The NeuroArtP3 (NET-2018-12366666) is a four-year multi-site project co-funded by the Italian Ministry of Health, bringing together clinical and computational centers operating in the field of neurology, including PD. OBJECTIVE: The core objectives of the project are: i) to harmonize the collection of data across the participating centers, ii) to structure standardized disease-specific datasets and iii) to advance knowledge on disease's trajectories through machine learning analysis. METHODS: The 4-years study combines two consecutive research components: i) a multi-center retrospective observational phase; ii) a multi-center prospective observational phase. The retrospective phase aims at collecting data of the patients admitted at the participating clinical centers. Whereas the prospective phase aims at collecting the same variables of the retrospective study in newly diagnosed patients who will be enrolled at the same centers. RESULTS: The participating clinical centers are the Provincial Health Services (APSS) of Trento (Italy) as the center responsible for the PD study and the IRCCS San Martino Hospital of Genoa (Italy) as the promoter center of the NeuroartP3 project. The computational centers responsible for data analysis are the Bruno Kessler Foundation of Trento (Italy) with TrentinoSalute4.0 -Competence Center for Digital Health of the Province of Trento (Italy) and the LISCOMPlab University of Genoa (Italy). CONCLUSIONS: The work behind this observational study protocol shows how it is possible and viable to systematize data collection procedures in order to feed research and to advance the implementation of a P3 approach into the clinical practice through the use of AI models.

Data-driven modelling of neurodegenerative disease progression: thinking outside the black box.

Data-driven disease progression models are an emerging set of computational tools that reconstruct disease timelines for long-term chronic diseases, providing unique insights into disease processes and their underlying mechanisms. Such methods combine a priori human knowledge and assumptions with large-scale data processing and parameter estimation to infer long-term disease trajectories from short-term data. In contrast to 'black box' machine learning tools, data-driven disease progression models typically require fewer data and are inherently interpretable, thereby aiding disease understanding in addition to enabling classification, prediction and stratification. In this Review, we place the current landscape of data-driven disease progression models in a general framework and discuss their enhanced utility for constructing a disease timeline compared with wider machine learning tools that construct static disease profiles. We review the insights they have enabled across multiple neurodegenerative diseases, notably Alzheimer disease, for applications such as determining temporal trajectories of disease biomarkers, testing hypotheses about disease mechanisms and uncovering disease subtypes. We outline key areas for technological development and translation to a broader range of neuroscience and non-neuroscience applications. Finally, we discuss potential pathways and barriers to integrating disease progression models into clinical practice and trial settings.

Longitudinal Cluster Analysis of Hemodialysis Patients with COVID-19 in the Pre-Vaccination Era.

Coronavirus disease 2019 (COVID-19) in hemodialysis patients (HD) is characterized by heterogeneity of clinical presentation and outcomes. To stratify patients, we collected clinical and laboratory data in two cohorts of HD patients at COVID-19 diagnosis and during the following 4 weeks. Baseline and longitudinal values were used to build a linear mixed effect model (LME) and define different clusters. The development of the LME model in the derivation cohort of 17 HD patients (66.7 ± 12.3 years, eight males) allowed the characterization of two clusters (cl1 and cl2). Patients in cl1 presented a prevalence of females, higher lymphocyte count, and lower levels of lactate dehydrogenase, C-reactive protein, and CD8 + T memory stem cells as a possible result of a milder inflammation. Then, this model was tested in an independent validation cohort of 30 HD patients (73.3 ± 16.3 years, 16 males) assigned to cl1 or cl2 (16 and 14 patients, respectively). The cluster comparison confirmed that cl1 presented a milder form of COVID-19 associated with reduced disease activity, hospitalization, mortality rate, and oxygen requirement. Clustering analysis on longitudinal data allowed patient stratification and identification of the patients at high risk of complications. This strategy could be suitable in different clinical settings.

First Hard X-Ray Imaging Results by Solar Orbiter STIX.

The Spectrometer/Telescope for Imaging X-rays (STIX) is one of six remote sensing instruments on-board Solar Orbiter. The telescope applies an indirect imaging technique that uses the measurement of 30 visibilities, i.e., angular Fourier components of the solar flare X-ray source. Hence, the imaging problem for STIX consists of the Fourier inversion of the data measured by the instrument. In this work, we show that the visibility amplitude and phase calibration of 24 out of 30 STIX sub-collimators has reached a satisfactory level for scientific data exploitation and that a set of imaging methods is able to provide the first hard X-ray images of solar flares from Solar Orbiter. Four visibility-based image reconstruction methods and one count-based are applied to calibrated STIX observations of six events with GOES class between C4 and M4 that occurred in May 2021. The resulting reconstructions are compared to those provided by an optimization algorithm used for fitting the amplitudes of STIX visibilities. We show that the five imaging methods produce results morphologically consistent with the ones provided by the Atmospheric Imaging Assembly on board the Solar Dynamic Observatory (SDO/AIA) in UV wavelengths. The χ 2 values and the parameters of the reconstructed sources are comparable between methods, thus confirming their robustness.

Neurological long-COVID in the outpatient clinic: Two subtypes, two courses.

INTRODUCTION: Symptoms referable to central and peripheral nervous system involvement are often evident both during the acute phase of COVID-19 infection and during long-COVID. In this study, we evaluated a population of patients with prior COVID-19 infection who showed signs and symptoms consistent with neurological long-COVID. METHODS: We prospectively collected demographic and acute phase course data from patients with prior COVID-19 infection who showed symptoms related to neurological involvement in the long-COVID phase. Firstly, we performed a multivariate logistic linear regression analysis to investigate the impact of demographic and clinical data, the severity of the acute COVID-19 infection and hospitalization course, on the post-COVID neurological symptoms at three months follow-up. Secondly, we performed an unsupervised clustering analysis to investigate whether there was evidence of different subtypes of neurological long COVID-19. RESULTS: One hundred and nine patients referred to the neurological post-COVID outpatient clinic. Clustering analysis on the most common neurological symptoms returned two well-separated and well-balanced clusters: long-COVID type 1 contains the subjects with memory disturbances, psychological impairment, headache, anosmia and ageusia, while long-COVID type 2 contains all the subjects with reported symptoms related to PNS involvement. The analysis of potential risk-factors among the demographic, clinical presentation, COVID 19 severity and hospitalization course variables showed that the number of comorbidities at onset, the BMI, the number of COVID-19 symptoms, the number of non-neurological complications and a more severe course of the acute infection were all, on average, higher for the cluster of subjects with reported symptoms related to PNS involvement. CONCLUSION: We analyzed the characteristics of neurological long-COVID and presented a method to identify well-defined patient groups with distinct symptoms and risk factors. The proposed method could potentially enable treatment deployment by identifying the optimal interventions and services for well-defined patient groups, so alleviating long-COVID and easing recovery.

Revealing the Timeline of Structural MRI Changes in Premanifest to Manifest Huntington Disease.

BACKGROUND AND OBJECTIVES: Longitudinal measurements of brain atrophy using structural MRI (sMRI) can provide powerful markers for tracking disease progression in neurodegenerative diseases. In this study, we use a disease progression model to learn individual-level disease times and hence reveal a new timeline of sMRI changes in Huntington disease (HD). METHODS: We use data from the 2 largest cohort imaging studies in HD-284 participants from TRACK-HD (100 control, 104 premanifest, and 80 manifest) and 159 participants from PREDICT-HD (36 control and 128 premanifest)-to train and test the model. We longitudinally register T1-weighted sMRI scans from 3 consecutive time points to reduce intraindividual variability and calculate regional brain volumes using an automated segmentation tool with rigorous manual quality control. RESULTS: Our model reveals, for the first time, the relative magnitude and timescale of subcortical and cortical atrophy changes in HD. We find that the largest (∼20% average change in magnitude) and earliest (∼2 years before average abnormality) changes occur in the subcortex (pallidum, putamen, and caudate), followed by a cascade of changes across other subcortical and cortical regions over a period of ∼11 years. We also show that sMRI, when combined with our disease progression model, provides improved prediction of onset over the current best method (root mean square error = 4.5 years and maximum error = 7.9 years vs root mean square error = 6.6 years and maximum error = 18.2 years). DISCUSSION: Our findings support the use of disease progression modeling to reveal new information from sMRI, which can potentially inform imaging marker selection for clinical trials.

Positron emission tomography and single photon emission computed tomography imaging of tertiary lymphoid structures during the development of lupus nephritis.

Lymphoid neogenesis occurs in tissues targeted by chronic inflammatory processes, such as infection and autoimmunity. In systemic lupus erythematosus (SLE), such structures develop within the kidneys of lupus-prone mice ((NZBXNZW)F1) and are observed in kidney biopsies taken from SLE patients with lupus nephritis (LN). The purpose of this prospective longitudinal animal study was to detect early kidney changes and tertiary lymphoid structures (TLS) using in vivo imaging. Positron emission tomography (PET) by tail vein injection of 18-F-fluoro-2-deoxy-D-glucose (18F-FDG)(PET/FDG) combined with computed tomography (CT) for anatomical localization and single photon emission computed tomography (SPECT) by intraperitoneal injection of 99mTC labeled Albumin Nanocoll (99mTC-Nanocoll) were performed on different disease stages of NZB/W mice (n = 40) and on aged matched control mice (BALB/c) (n = 20). By using one-way ANOVA analyses, we compared two different compartmental models for the quantitative measure of 18F-FDG uptake within the kidneys. Using a new five-compartment model, we observed that glomerular filtration of 18FFDG in lupus-prone mice decreased significantly by disease progression measured by anti-dsDNA Ab production and before onset of proteinuria. We could not visualize TLS within the kidneys, but we were able to visualize pancreatic TLS using 99mTC Nanocoll SPECT. Based on our findings, we conclude that the five-compartment model can be used to measure changes of FDG uptake within the kidney. However, new optimal PET/SPECT tracer administration sites together with more specific tracers in combination with magnetic resonance imaging (MRI) may make it possible to detect formation of TLS and LN before clinical manifestations.

Investigating hypotheses of neurodegeneration by learning dynamical systems of protein propagation in the brain.

We introduce a theoretical framework for estimating, comparing and interpreting mechanistic hypotheses on long term protein propagation across brain networks in neurodegenerative disorders (ND). The model is expressed within a Bayesian non-parametric regression setting, where mechanisms of protein dynamics are inferred by means of gradient matching on dynamical systems (DS). The Bayesian formalism, combined with stochastic variational inference, naturally allows for model comparison via assessment of model evidence, while providing uncertainty quantification of causal relationship underlying protein progressions. When applied to in-vivo AV45-PET brain imaging data measuring topographic amyloid deposition in Alzheimer's disease (AD), our model identified the mechanisms of accumulation, clearance and propagation as the best suited DS for bio-mechanical description of amyloid dynamics in AD, enabling realistic and accurate personalized simulation of amyloidosis.

Prion propagation estimated from brain diffusion MRI is subtype dependent in sporadic Creutzfeldt-Jakob disease.

Sporadic Creutzfeldt-Jakob disease (sCJD) is a transmissible brain proteinopathy. Five main clinicopathological subtypes (sCJD-MM(V)1, -MM(V)2C, -MV2K, -VV1, and -VV2) are currently distinguished. Histopathological evidence suggests that the localisation of prion aggregates and spongiform lesions varies among subtypes. Establishing whether there is an initial site with detectable imaging abnormalities (epicentre) and an order of lesion propagation would be informative for disease early diagnosis, patient staging, management and recruitment in clinical trials. Diffusion magnetic resonance imaging (MRI) is the most-used and most-sensitive test to detect spongiform degeneration. This study was designed to identify, in vivo and for the first time, subtype-dependent epicentre and lesion propagation in the brain using diffusion-weighted images (DWI), in the largest known cross-sectional dataset of autopsy-proven subjects with sCJD. We estimate lesion propagation by cross-sectional DWI using event-based modelling, a well-established data-driven technique. DWI abnormalities of 594 autopsy-diagnosed subjects (448 patients with sCJD) were scored in 12 brain regions by 1 neuroradiologist blind to the diagnosis. We used the event-based model to reconstruct sequential orderings of lesion propagation in each of five pure subtypes. Follow-up data from 151 patients validated the estimated sequences. Results showed that epicentre and ordering of lesion propagation are subtype specific. The two most common subtypes (-MM1 and -VV2) showed opposite ordering of DWI abnormality appearance: from the neocortex to subcortical regions, and vice versa, respectively. The precuneus was the most likely epicentre also in -MM2 and -VV1 although at variance with -MM1, abnormal signal was also detected early in cingulate and insular cortices. The caudal-rostral sequence of lesion propagation that characterises -VV2 was replicated in -MV2K. Combined, these data-driven models provide unprecedented dynamic insights into subtype-specific epicentre at onset and propagation of the pathologic process, which may also enhance early diagnosis and enable disease staging in sCJD.

Differences in topological progression profile among neurodegenerative diseases from imaging data.

The spatial distribution of atrophy in neurodegenerative diseases suggests that brain connectivity mediates disease propagation. Different descriptors of the connectivity graph potentially relate to different underlying mechanisms of propagation. Previous approaches for evaluating the influence of connectivity on neurodegeneration consider each descriptor in isolation and match predictions against late-stage atrophy patterns. We introduce the notion of a topological profile - a characteristic combination of topological descriptors that best describes the propagation of pathology in a particular disease. By drawing on recent advances in disease progression modeling, we estimate topological profiles from the full course of pathology accumulation, at both cohort and individual levels. Experimental results comparing topological profiles for Alzheimer's disease, multiple sclerosis and normal ageing show that topological profiles explain the observed data better than single descriptors. Within each condition, most individual profiles cluster around the cohort-level profile, and individuals whose profiles align more closely with other cohort-level profiles show features of that cohort. The cohort-level profiles suggest new insights into the biological mechanisms underlying pathology propagation in each disease.

DIVE: A spatiotemporal progression model of brain pathology in neurodegenerative disorders.

Current models of progression in neurodegenerative diseases use neuroimaging measures that are averaged across pre-defined regions of interest (ROIs). Such models are unable to recover fine details of atrophy patterns; they tend to impose an assumption of strong spatial correlation within each ROI and no correlation among ROIs. Such assumptions may be violated by the influence of underlying brain network connectivity on pathology propagation - a strong hypothesis e.g. in Alzheimer's Disease. Here we present DIVE: Data-driven Inference of Vertexwise Evolution. DIVE is an image-based disease progression model with single-vertex resolution, designed to reconstruct long-term patterns of brain pathology from short-term longitudinal data sets. DIVE clusters vertex-wise (i.e. point-wise) biomarker measurements on the cortical surface that have similar temporal dynamics across a patient population, and concurrently estimates an average trajectory of vertex measurements in each cluster. DIVE uniquely outputs a parcellation of the cortex into areas with common progression patterns, leading to a new signature for individual diseases. DIVE further estimates the disease stage and progression speed for every visit of every subject, potentially enhancing stratification for clinical trials or management. On simulated data, DIVE can recover ground truth clusters and their underlying trajectory, provided the average trajectories are sufficiently different between clusters. We demonstrate DIVE on data from two cohorts: the Alzheimer's Disease Neuroimaging Initiative (ADNI) and the Dementia Research Centre (DRC), UK. The DRC cohort contains patients with Posterior Cortical Atrophy (PCA) as well as typical Alzheimer's disease (tAD). DIVE finds similar spatial patterns of atrophy for tAD subjects in the two independent datasets (ADNI and DRC), and further reveals distinct patterns of pathology in different diseases (tAD vs PCA) and for distinct types of biomarker data - cortical thickness from Magnetic Resonance Imaging (MRI) vs amyloid load from Positron Emission Tomography (PET). We demonstrate that DIVE stages have potential clinical relevance, despite being based only on imaging data, by showing that the stages correlate with cognitive test scores. Finally, DIVE can be used to estimate a fine-grained spatial distribution of pathology in the brain using any kind of voxelwise or vertexwise measures including Jacobian compression maps, fractional anisotropy (FA) maps from diffusion tensor imaging (DTI) or other PET measures.

Data-Driven Sequence of Changes to Anatomical Brain Connectivity in Sporadic Alzheimer's Disease.

Model-based investigations of transneuronal spreading mechanisms in neurodegenerative diseases relate the pattern of pathology severity to the brain's connectivity matrix, which reveals information about how pathology propagates through the connectivity network. Such network models typically use networks based on functional or structural connectivity in young and healthy individuals, and only end-stage patterns of pathology, thereby ignoring/excluding the effects of normal aging and disease progression. Here, we examine the sequence of changes in the elderly brain's anatomical connectivity over the course of a neurodegenerative disease. We do this in a data-driven manner that is not dependent upon clinical disease stage, by using event-based disease progression modeling. Using data from the Alzheimer's Disease Neuroimaging Initiative dataset, we sequence the progressive decline of anatomical connectivity, as quantified by graph-theory metrics, in the Alzheimer's disease brain. Ours is the first single model to contribute to understanding all three of the nature, the location, and the sequence of changes to anatomical connectivity in the human brain due to Alzheimer's disease. Our experimental results reveal new insights into Alzheimer's disease: that degeneration of anatomical connectivity in the brain may be a viable, even early, biomarker and should be considered when studying such neurodegenerative diseases.

Expectation maximization and the retrieval of the atmospheric extinction coefficients by inversion of Raman lidar data.

We consider the problem of retrieving the aerosol extinction coefficient from Raman lidar measurements. This is an ill-posed inverse problem that needs regularization, and we propose to use the Expectation-Maximization (EM) algorithm to provide stable solutions. Indeed, EM is an iterative algorithm that imposes a positivity constraint on the solution, and provides regularization if iterations are stopped early enough. We describe the algorithm and propose a stopping criterion inspired by a statistical principle. We then discuss its properties concerning the spatial resolution. Finally, we validate the proposed approach by using both synthetic data and experimental measurements; we compare the reconstructions obtained by EM with those obtained by the Tikhonov method, by the Levenberg-Marquardt method, as well as those obtained by combining data smoothing and numerical derivation.

A new compartmental method for the analysis of liver FDG kinetics in small animal models.

BACKGROUND: Compartmental analysis is a standard method to quantify metabolic processes using fluorodeoxyglucose-positron emission tomography (FDG-PET). For liver studies, this analysis is complex due to the hepatocyte capability to dephosphorylate and release glucose and FDG into the blood. Moreover, a tracer is supplied to the liver by both the hepatic artery and the portal vein, which is not visible in PET images. This study developed an innovative computational approach accounting for the reversible nature of FDG in the liver and directly computing the portal vein tracer concentration by means of gut radioactivity measurements. METHODS: Twenty-one mice were subdivided into three groups: the control group 'CTR' (n = 7) received no treatment, the short-term starvation group 'STS' (n = 7) was submitted to food deprivation with free access to water within 48 h before imaging, and the metformin group 'MTF' (n = 7) was treated with metformin (750 mg/Kg per day) for 1 month. All mice underwent a dynamic micro-PET study for 50 min after an (18)F-FDG injection. The compartmental analysis considered two FDG pools (phosphorylated and free) in both the gut and liver. A tracer was carried into the liver by the hepatic artery and the portal vein, and tracer delivery from the gut was considered as the sole input for portal vein tracer concentration. Accordingly, both the liver and gut were characterized by two compartments and two exchange coefficients. Each one of the two two-compartment models was mathematically described by a system of differential equations, and data optimization was performed by applying a Newton algorithm to the inverse problems associated to these differential systems. RESULTS: All rate constants were stable in each group. The tracer coefficient from the free to the metabolized compartment in the liver was increased by STS, while it was unaltered by MTF. By contrast, the tracer coefficient from the metabolized to the free compartment was reduced by MTF and increased by STS. CONCLUSIONS: Data demonstrated that our method was able to analyze FDG kinetics under pharmacological or pathophysiological stimulation, quantifying the fraction of the tracer trapped in the liver or dephosphorylated and released into the bloodstream.

Estimate of FDG excretion by means of compartmental analysis and ant colony optimization of nuclear medicine data.

[(18)F]fluoro-2-deoxy-D-glucose (FDG) is one of the most utilized tracers for positron emission tomography (PET) applications in oncology. FDG-PET relies on higher glycolytic activity in tumors compared to normal structures as the basis of image contrast. As a glucose analog, FDG is transported into malignant cells which typically exhibit an increased radioactivity. However, different from glucose, FDG is not reabsorbed by the renal system and is excreted to the bladder. The present paper describes a novel computational method for the quantitative assessment of this excretion process. The method is based on a compartmental analysis of FDG-PET data in which the excretion process is explicitly accounted for by the bladder compartment and on the application of an ant colony optimization (ACO) algorithm for the determination of the tracer coefficients describing the FDG transport effectiveness. The validation of this approach is performed by means of both synthetic data and real measurements acquired by a PET device for small animals (micro-PET). Possible oncological applications of the results are discussed in the final section.