Lipid distributions in the Global Diagnostics Network across five continents.

AIMS: Lipids are central in the development of cardiovascular disease, and the present study aimed to characterize variation in lipid profiles across different countries to improve understanding of cardiovascular risk and opportunities for risk-reducing interventions. METHODS AND RESULTS: This first collaborative report of the Global Diagnostics Network (GDN) evaluated lipid distributions from nine laboratory organizations providing clinical laboratory testing in 17 countries on five continents. This cross-sectional study assessed aggregated lipid results from patients aged 20-89 years, tested at GDN laboratories, from 2018 through 2020. In addition to mean levels, the World Health Organization total cholesterol risk target (<5.00 mmol/L, <193 mg/dL) and proportions in guideline-based low-density lipoprotein cholesterol (LDL-C) categories were assessed. This study of 461 888 753 lipid results found wide variation by country/region, sex, and age. In most countries, total cholesterol and LDL-C peaked at 50-59 years in females and 40-49 years in males. Sex- and age-group adjusted mean total cholesterol levels ranged from 4.58 mmol/L (177.1 mg/dL) in the Republic of Korea to 5.40 mmol/L (208.8 mg/dL) in Austria. Mean total cholesterol levels exceeded the World Health Organization target in Japan, Australia, North Macedonia, Switzerland, Germany, Slovakia, and Austria. Considering LDL-C categories, North Macedonia had the highest proportions of LDL-C results >4.91 mmol/L (>190 mg/dL) for both females (9.9%) and males (8.7%). LDL-C levels <1.55 mmol/L (<60 mg/dL) were most common among females in Canada (10.7%) and males in the UK (17.3%). CONCLUSION: With nearly a half billion lipid results, this study sheds light on the worldwide variability in lipid levels, which may reflect inter-country differences in genetics, lipid testing, lifestyle habits, and pharmacologic treatment. Despite variability, elevated atherogenic lipid levels are a common global problem, and these results can help inform national policies and health system approaches to mitigate lipid-mediated risk of cardiovascular disease.

Translating AI to Clinical Practice: Overcoming Data Shift with Explainability.

To translate artificial intelligence (AI) algorithms into clinical practice requires generalizability of models to real-world data. One of the main obstacles to generalizability is data shift, a data distribution mismatch between model training and real environments. Explainable AI techniques offer tools to detect and mitigate the data shift problem and develop reliable AI for clinical practice. Most medical AI is trained with datasets gathered from limited environments, such as restricted disease populations and center-dependent acquisition conditions. The data shift that commonly exists in the limited training set often causes a significant performance decrease in the deployment environment. To develop a medical application, it is important to detect potential data shift and its impact on clinical translation. During AI training stages, from premodel analysis to in-model and post hoc explanations, explainability can play a key role in detecting model susceptibility to data shift, which is otherwise hidden because the test data have the same biased distribution as the training data. Performance-based model assessments cannot effectively distinguish the model overfitting to training data bias without enriched test sets from external environments. In the absence of such external data, explainability techniques can aid in translating AI to clinical practice as a tool to detect and mitigate potential failures due to data shift. ©RSNA, 2023 Quiz questions for this article are available in the supplemental material.

Do Microsurgical Outcomes Differ Based on Which Specialty Does the Operation? A NSQIP Analysis.

BACKGROUND: Because plastic surgeons do not "own" a specific anatomic region, other surgical specialties have increasingly assumed procedures historically performed by plastic surgery. Decreased case volume is postulated to be associated with higher complication rates. Herein, we investigate whether volume and surgical specialty have an impact on microsurgical complications, specifically surgical site infection (SSI) and reoperation rates. METHODS: The 2005-2015 National Surgical Quality Improvement Program participant use file was queried by Current Procedural Terminology code for breast and head/neck microsurgeries. Multivariate logistic regression was performed to compare the outcomes between surgical specialties. A cumulative frequency variable was introduced to investigate the effect of case volume on complication rates. RESULTS: We captured 6,617 microsurgical cases. Multivariate logistic regression revealed that although the rate of SSI was lower in plastic surgery compared with otolaryngology for head and neck reconstructions (13.3% versus 10.5%) and compared with general surgery for breast reconstructions (5.4% versus 4.7%), there was no significant difference between specialties (P = 0.13; P = 0.96). Increased case volume is negatively correlated with complications. CONCLUSIONS: Plastic surgery is at risk given case cannibalization by other specialties. We conclude that surgical specialty does not affect the rates of SSI and reoperation. We demonstrate a correlation between lower volumes and increased complications, implying that, once a specialty has amassed critical case experience, complication rates may decrease, and outcomes can be equivalent or superior. Case breadth and volumes should be maintained to preserve skills, optimize outcomes, and maintain the specialty as it currently exists.

High-Risk Plastic Surgery: An Analysis of 108,303 Cases From the American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP).

BACKGROUND: There is a lack of large-scale data that examine complications in plastic surgery. A description of baseline rates and patient outcomes allows better understanding of ways to improve patient care and cost-savings for health systems. Herein, we determine the most frequent complications in plastic surgery, identify procedures with high complication rates, and examine predictive risk factors. METHODS: A retrospective analysis of the 2012 to 2016 American College of Surgeons National Surgical Quality Improvement Program plastic surgery data set was conducted. Complication rates were calculated for the entire cohort and each procedure therein. Microsurgical procedures were analyzed as a subgroup, where multivariate logistic regression models determined the risk factors for surgical site infection (SSI) and related reoperation. RESULTS: We identified 108 303 patients undergoing a plastic surgery procedure of which 6 264 (5.78%) experienced ≥1 complication. The outcome with the highest incidence was related reoperation (3.31%), followed by SSI (3.11%). Microsurgical cases comprised 6 148 (5.68%) of all cases, and 1211 (19.33%) experienced ≥1 complication. Similar to the entire cohort, the related reoperation (12.83%) and SSI (5.66%) were common complications. Increased operative time was a common independent risk factor predictive of a related reoperation or development of an SSI (P < 001). Of all microsurgeries, 23.3% had an operative time larger than 10 hours which lead to faster increase in reoperation likelihood. CONCLUSIONS: The complication rate in plastic surgery remains relatively low but is significantly increased for microsurgery. Increased operative time is a common risk factor. Two-team approaches and staged operations could be explored, as a large portion of microsurgeries are vulnerable to increased complications.

HISTORIQUE: Les données à grande échelle sur les complications de la chirurgie plastique font défaut. Une description des taux de référence et des résultats cliniques des patients permettrait de mieux déterminer comment améliorer les soins aux patients et réaliser des économies dans les systèmes de santé. Dans le présent article, les chercheurs recensent les complications les plus fréquentes en chirurgie plastique, dégagent les interventions aux taux de complication élevés et examinent les facteurs de risque prédictifs. MÉTHODOLOGIE: Les chercheurs ont réalisé une analyse rétrospective des données de chirurgie plastique tirées du programme national d'amélioration de la qualité chirurgicale de l'American College of Surgeons entre 2012 et 2016. Ils ont calculé les taux de complications de toute la cohorte et de chaque intervention recensée. Ils ont analysé les interventions microchirurgicales en sous-groupe, où ils ont utilisé des modèles de régression logistique multivariée pour déterminer les facteurs de risque d'infection des plaies opératoires (IPO) et de réopérations s'y rapportant. RÉSULTATS: Les chercheurs ont dénombré 108 303 patients qui avaient subi une intervention en chirurgie plastique, dont 6 264 (5,78 %) avaient souffert d'au moins une complication. Les réopérations (3,31 %), suivies des IPO (3,11 %) étaient les résultats à la plus forte incidence. Les cas de microchirurgie représentaient 6 148 (5,68 %) de toutes les occurrences, et 1211 (19,33 %) ont souffert d'au moins une complication. Tout comme dans l'ensemble de la cohorte, les réopérations (12,83 %) et les IPO (5,66 %) étaient des complications courantes. La plus longue durée de l'opération était un facteur de risque indépendant fréquent, prédicteur d'une réopération ou d'une IPO (p<0,001). Ainsi, 23,3 % des microchirurgies duraient plus de dix heures, ce qui s'associait à une plus forte augmentation du risque de réopération. CONCLUSIONS: Le taux de complications demeure relativement faible en chirurgie plastique, mais est significativement plus élevé en microchirurgie. La longue durée des opérations représente un facteur de risque courant. On pourrait explorer les approches à deux équipes et les opérations échelonnées, car une forte proportion des microchirurgies sont vulnérables à un accroissement des complications.

Association of intensive care unit occupancy during admission and inpatient mortality: a retrospective cohort study.

PURPOSE: There is conflicting evidence regarding the influence of intensive care unit (ICU) occupancy at the time of admission on important patient outcomes such as mortality. The objective of this analysis was to characterize the association between ICU occupancy at the time of ICU admission and subsequent mortality. METHODS: This single-centre, retrospective cohort study included all patients admitted to the ICU at the Vancouver General Hospital between 4 January 2010 and 8 October 2017. Intensive care unit occupancy was defined as the number of ICU bed hours utilized in a day divided by the total amount of ICU bed hours available for that day. We constructed mixed-effects logistic regression models controlling for relevant covariates to assess the impact of admission occupancy quintiles on total inpatient (ICU and ward) and early (72-hr) ICU mortality. RESULTS: This analysis included 10,365 ICU admissions by 8,562 unique patients. Compared with ICU admissions in the median occupancy quintile, admissions in the highest and second highest occupancy quintile were associated with a significant increase in the odds of inpatient mortality (highest: odds ratio [OR], 1.33; 95% confidence interval [CI], 1.12 to 1.59; P value < 0.001; second highest: OR, 1.21; 95% CI, 1.02 to 1.44; P value < 0.03). No association between admission occupancy and 72-hr ICU mortality was observed. CONCLUSIONS: Admission to the ICU on days of high occupancy was associated with increased inpatient mortality, but not with increased 72-hr ICU mortality. Capacity strain on the ICU may result in significant negative consequences for patients, but further research is needed to fully characterize the complex effects of capacity strain.

Deep learning and radiomics: the utility of Google TensorFlow™ Inception in classifying clear cell renal cell carcinoma and oncocytoma on multiphasic CT.

PURPOSE: Currently, all solid enhancing renal masses without microscopic fat are considered malignant until proven otherwise and there is substantial overlap in the imaging findings of benign and malignant renal masses, particularly between clear cell RCC (ccRCC) and benign oncocytoma (ONC). Radiomics has attracted increased attention for its utility in pre-operative work-up on routine clinical images. Radiomics based approaches have converted medical images into mineable data and identified prognostic imaging signatures that machine learning algorithms can use to construct predictive models by learning the decision boundaries of the underlying data distribution. The TensorFlow™ framework from Google is a state-of-the-art open-source software library that can be used for training deep learning neural networks for performing machine learning tasks. The purpose of this study was to investigate the diagnostic value and feasibility of a deep learning-based renal lesion classifier using open-source Google TensorFlow™ Inception in differentiating ccRCC from ONC on routine four-phase MDCT in patients with pathologically confirmed renal masses. METHODS: With institutional review board approval for this 1996 Health Insurance Portability and Accountability Act compliant retrospective study and a waiver of informed consent, we queried our institution's pathology, clinical, and radiology databases for histologically proven cases of ccRCC and ONC obtained between January 2000 and January 2016 scanned with a an intravenous contrast-enhanced four-phase renal mass protocol (unenhanced (UN), corticomedullary (CM), nephrographic (NP), and excretory (EX) phases). To extract features to be used for the machine learning model, the entire renal mass was contoured in the axial plane in each of the four phases, resulting in a 3D volume of interest (VOI) representative of the entire renal mass. We investigated thirteen different approaches to convert the acquired VOI data into a set of images that adequately represented each tumor which was used to train the final layer of the neural network model. Training was performed over 4000 iterations. In each iteration, 90% of the data were designated as training data and the remaining 10% served as validation data and a leave-one-out cross-validation scheme was implemented. Accuracy, sensitivity, specificity, positive (PPV) and negative predictive (NPV) values, and CIs were calculated for the classification of the thirteen processing modes. RESULTS: We analyzed 179 consecutive patients with 179 lesions (128 ccRCC and 51 ONC). The ccRCC cohort had a mean size of 3.8 cm (range 0.8-14.6 cm) and the ONC cohort had a mean lesion size of 3.9 cm (range 1.0-13.1 cm). The highest specificity and PPV (52.9% and 80.3%, respectively) were achieved in the EX phase when we analyzed the single mid-slice of the tumor in the axial, coronal and sagittal plane, and when we increased the number of mid-slices of the tumor to three, with an accuracy of 75.4%, which also increased the sensitivity to 88.3% and the PPV to 79.6%. Using the entire tumor volume also showed that classification performance was best in the EX phase with an accuracy of 74.4%, a sensitivity of 85.8% and a PPV of 80.1%. When the entire tumor volume, plus mid-slices from all phases and all planes presented as tiled images, were submitted to the final layer of the neural network we achieved a PPV of 82.5%. CONCLUSIONS: The best classification result was obtained in the EX phase among the thirteen classification methods tested. Our proof of concept study is the first step towards understanding the utility of machine learning in the differentiation of ccRCC from ONC on routine CT images. We hope this could lead to future investigation into the development of a multivariate machine learning model which may augment our ability to accurately predict renal lesion histology on imaging.

Alteration of brain network topology in HIV-associated neurocognitive disorder: A novel functional connectivity perspective.

HIV is capable of invading the brain soon after seroconversion. This ultimately can lead to deficits in multiple cognitive domains commonly referred to as HIV-associated neurocognitive disorders (HAND). Clinical diagnosis of such deficits requires detailed neuropsychological assessment but clinical signs may be difficult to detect during asymptomatic injury of the central nervous system (CNS). Therefore neuroimaging biomarkers are of particular interest in HAND. In this study, we constructed brain connectivity profiles of 40 subjects (20 HIV positive subjects and 20 age-matched seronegative controls) using two different methods: a non-linear mutual connectivity analysis approach and a conventional method based on Pearson's correlation. These profiles were then summarized using graph-theoretic methods characterizing their topological network properties. Standard clinical and laboratory assessments were performed and a battery of neuropsychological (NP) tests was administered for all participating subjects. Based on NP testing, 14 of the seropositive subjects exhibited mild neurologic impairment. Subsequently, we analyzed associations between the network derived measures and neuropsychological assessment scores as well as common clinical laboratory plasma markers (CD4 cell count, HIV RNA) after adjusting for age and gender. Mutual connectivity analysis derived graph-theoretic measures, Modularity and Small Worldness, were significantly (p < 0.05, FDR adjusted) associated with the Executive as well as Overall z-score of NP performance. In contrast, network measures derived from conventional correlation-based connectivity did not yield any significant results. Thus, changes in connectivity can be captured using advanced time-series analysis techniques. The demonstrated associations between imaging-derived graph-theoretic properties of brain networks with neuropsychological performance, provides opportunities to further investigate the evolution of HAND in larger, longitudinal studies. Our analysis approach, involving non-linear time-series analysis in conjunction with graph theory, is promising and it may prove to be useful not only in HAND but also in other neurodegenerative disorders.

Building a high-resolution T2-weighted MR-based probabilistic model of tumor occurrence in the prostate.

PURPOSE: We present a method for generating a T2 MR-based probabilistic model of tumor occurrence in the prostate to guide the selection of anatomical sites for targeted biopsies and serve as a diagnostic tool to aid radiological evaluation of prostate cancer. MATERIALS AND METHODS: In our study, the prostate and any radiological findings within were segmented retrospectively on 3D T2-weighted MR images of 266 subjects who underwent radical prostatectomy. Subsequent histopathological analysis determined both the ground truth and the Gleason grade of the tumors. A randomly chosen subset of 19 subjects was used to generate a multi-subject-derived prostate template. Subsequently, a cascading registration algorithm involving both affine and non-rigid B-spline transforms was used to register the prostate of every subject to the template. Corresponding transformation of radiological findings yielded a population-based probabilistic model of tumor occurrence. The quality of our probabilistic model building approach was statistically evaluated by measuring the proportion of correct placements of tumors in the prostate template, i.e., the number of tumors that maintained their anatomical location within the prostate after their transformation into the prostate template space. RESULTS: Probabilistic model built with tumors deemed clinically significant demonstrated a heterogeneous distribution of tumors, with higher likelihood of tumor occurrence at the mid-gland anterior transition zone and the base-to-mid-gland posterior peripheral zones. Of 250 MR lesions analyzed, 248 maintained their original anatomical location with respect to the prostate zones after transformation to the prostate. CONCLUSION: We present a robust method for generating a probabilistic model of tumor occurrence in the prostate that could aid clinical decision making, such as selection of anatomical sites for MR-guided prostate biopsies.

Deep transfer learning for characterizing chondrocyte patterns in phase contrast X-Ray computed tomography images of the human patellar cartilage.

Phase contrast X-ray computed tomography (PCI-CT) has been demonstrated to be effective for visualization of the human cartilage matrix at micrometer resolution, thereby capturing osteoarthritis induced changes to chondrocyte organization. This study aims to systematically assess the efficacy of deep transfer learning methods for classifying between healthy and diseased tissue patterns. We extracted features from two different convolutional neural network architectures, CaffeNet and Inception-v3 for characterizing such patterns. These features were quantitatively evaluated in a classification task measured by the area (AUC) under the Receiver Operating Characteristic (ROC) curve as well as qualitative visualization through a dimension reduction approach t-Distributed Stochastic Neighbor Embedding (t-SNE). The best classification performance, for CaffeNet, was observed when using features from the last convolutional layer and the last fully connected layer (AUCs >0.91). Meanwhile, off-the-shelf features from Inception-v3 produced similar classification performance (AUC >0.95). Visualization of features from these layers further confirmed adequate characterization of chondrocyte patterns for reliably distinguishing between healthy and osteoarthritic tissue classes. Such techniques, can be potentially used for detecting the presence of osteoarthritis related changes in the human patellar cartilage.

Detecting Altered connectivity patterns in HIV associated neurocognitive impairment using Mutual Connectivity Analysis.

The use of functional Magnetic Resonance Imaging (fMRI) has provided interesting insights into our understanding of the brain. In clinical setups these scans have been used to detect and study changes in the brain network properties in various neurological disorders. A large percentage of subjects infected with HIV present cognitive deficits, which are known as HIV associated neurocognitive disorder (HAND). In this study we propose to use our novel technique named Mutual Connectivity Analysis (MCA) to detect differences in brain networks in subjects with and without HIV infection. Resting state functional MRI scans acquired from 10 subjects (5 HIV+ and 5 HIV-) were subject to standard pre-processing routines. Subsequently, the average time-series for each brain region of the Automated Anatomic Labeling (AAL) atlas are extracted and used with the MCA framework to obtain a graph characterizing the interactions between them. The network graphs obtained for different subjects are then compared using Network-Based Statistics (NBS), which is an approach to detect differences between graphs edges while controlling for the family-wise error rate when mass univariate testing is performed. Applying this approach on the graphs obtained yields a single network encompassing 42 nodes and 65 edges, which is significantly different between the two subject groups. Specifically connections to the regions in and around the basal ganglia are significantly decreased. Also some nodes corresponding to the posterior cingulate cortex are affected. These results are inline with our current understanding of pathophysiological mechanisms of HIV associated neurocognitive disease (HAND) and other HIV based fMRI connectivity studies. Hence, we illustrate the applicability of our novel approach with network-based statistics in a clinical case-control study to detect differences connectivity patterns.

Assessing vertebral fracture risk on volumetric quantitative computed tomography by geometric characterization of trabecular bone structure.

The current clinical standard for measuring Bone Mineral Density (BMD) is dual X-ray absorptiometry, however more recently BMD derived from volumetric quantitative computed tomography has been shown to demonstrate a high association with spinal fracture susceptibility. In this study, we propose a method of fracture risk assessment using structural properties of trabecular bone in spinal vertebrae. Experimental data was acquired via axial multi-detector CT (MDCT) from 12 spinal vertebrae specimens using a whole-body 256-row CT scanner with a dedicated calibration phantom. Common image processing methods were used to annotate the trabecular compartment in the vertebral slices creating a circular region of interest (ROI) that excluded cortical bone for each slice. The pixels inside the ROI were converted to values indicative of BMD. High dimensional geometrical features were derived using the scaling index method (SIM) at different radii and scaling factors (SF). The mean BMD values within the ROI were then extracted and used in conjunction with a support vector machine to predict the failure load of the specimens. Prediction performance was measured using the root-mean-square error (RMSE) metric and determined that SIM combined with mean BMD features (RMSE = 0.82 ± 0.37) outperformed MDCT-measured mean BMD (RMSE = 1.11 ± 0.33) (p < 10-4). These results demonstrate that biomechanical strength prediction in vertebrae can be significantly improved through the use of SIM-derived texture features from trabecular bone.

Investigating Changes in Brain Network Properties in HIV-Associated Neurocognitive Disease (HAND) using Mutual Connectivity Analysis (MCA).

About 50% of subjects infected with HIV present deficits in cognitive domains, which are known collectively as HIV associated neurocognitive disorder (HAND). The underlying synaptodendritic damage can be captured using resting state functional MRI, as has been demonstrated by a few earlier studies. Such damage may induce topological changes of brain connectivity networks. We test this hypothesis by capturing the functional interdependence of 90 brain network nodes using a Mutual Connectivity Analysis (MCA) framework with non-linear time series modeling based on Generalized Radial Basis function (GRBF) neural networks. The network nodes are selected based on the regions defined in the Automated Anatomic Labeling (AAL) atlas. Each node is represented by the average time series of the voxels of that region. The resulting networks are then characterized using graph-theoretic measures that quantify various network topology properties at a global as well as at a local level. We tested for differences in these properties in network graphs obtained for 10 subjects (6 male and 4 female, 5 HIV+ and 5 HIV-). Global network properties captured some differences between these subject cohorts, though significant differences were seen only with the clustering coefficient measure. Local network properties, such as local efficiency and the degree of connections, captured significant differences in regions of the frontal lobe, precentral and cingulate cortex amongst a few others. These results suggest that our method can be used to effectively capture differences occurring in brain network connectivity properties revealed by resting-state functional MRI in neurological disease states, such as HAND.

Mutual Connectivity Analysis (MCA) Using Generalized Radial Basis Function Neural Networks for Nonlinear Functional Connectivity Network Recovery in Resting-State Functional MRI.

We investigate the applicability of a computational framework, called mutual connectivity analysis (MCA), for directed functional connectivity analysis in both synthetic and resting-state functional MRI data. This framework comprises of first evaluating non-linear cross-predictability between every pair of time series prior to recovering the underlying network structure using community detection algorithms. We obtain the non-linear cross-prediction score between time series using Generalized Radial Basis Functions (GRBF) neural networks. These cross-prediction scores characterize the underlying functionally connected networks within the resting brain, which can be extracted using non-metric clustering approaches, such as the Louvain method. We first test our approach on synthetic models with known directional influence and network structure. Our method is able to capture the directional relationships between time series (with an area under the ROC curve = 0.92 ± 0.037) as well as the underlying network structure (Rand index = 0.87 ± 0.063) with high accuracy. Furthermore, we test this method for network recovery on resting-state fMRI data, where results are compared to the motor cortex network recovered from a motor stimulation sequence, resulting in a strong agreement between the two (Dice coefficient = 0.45). We conclude that our MCA approach is effective in analyzing non-linear directed functional connectivity and in revealing underlying functional network structure in complex systems.

Volumetric quantitative characterization of human patellar cartilage with topological and geometrical features on phase-contrast X-ray computed tomography.

Phase-contrast X-ray computed tomography (PCI-CT) has attracted significant interest in recent years for its ability to provide significantly improved image contrast in low absorbing materials such as soft biological tissue. In the research context of cartilage imaging, previous studies have demonstrated the ability of PCI-CT to visualize structural details of human patellar cartilage matrix and capture changes to chondrocyte organization induced by osteoarthritis. This study evaluates the use of geometrical and topological features for volumetric characterization of such chondrocyte patterns in the presence (or absence) of osteoarthritic damage. Geometrical features derived from the scaling index method (SIM) and topological features derived from Minkowski Functionals were extracted from 1392 volumes of interest (VOI) annotated on PCI-CT images of ex vivo human patellar cartilage specimens. These features were subsequently used in a machine learning task with support vector regression to classify VOIs as healthy or osteoarthritic; classification performance was evaluated using the area under the receiver operating characteristic curve (AUC). Our results show that the classification performance of SIM-derived geometrical features (AUC: 0.90 ± 0.09) is significantly better than Minkowski Functionals volume (AUC: 0.54 ± 0.02), surface (AUC: 0.72 ± 0.06), mean breadth (AUC: 0.74 ± 0.06) and Euler characteristic (AUC: 0.78 ± 0.04) (p < 10(-4)). These results suggest that such geometrical features can provide a detailed characterization of the chondrocyte organization in the cartilage matrix in an automated manner, while also enabling classification of cartilage as healthy or osteoarthritic with high accuracy. Such features could potentially serve as diagnostic imaging markers for evaluating osteoarthritis progression and its response to different therapeutic intervention strategies.

Integrating dimension reduction and out-of-sample extension in automated classification of ex vivo human patellar cartilage on phase contrast X-ray computed tomography.

Phase contrast X-ray computed tomography (PCI-CT) has been demonstrated as a novel imaging technique that can visualize human cartilage with high spatial resolution and soft tissue contrast. Different textural approaches have been previously investigated for characterizing chondrocyte organization on PCI-CT to enable classification of healthy and osteoarthritic cartilage. However, the large size of feature sets extracted in such studies motivates an investigation into algorithmic feature reduction for computing efficient feature representations without compromising their discriminatory power. For this purpose, geometrical feature sets derived from the scaling index method (SIM) were extracted from 1392 volumes of interest (VOI) annotated on PCI-CT images of ex vivo human patellar cartilage specimens. The extracted feature sets were subject to linear and non-linear dimension reduction techniques as well as feature selection based on evaluation of mutual information criteria. The reduced feature set was subsequently used in a machine learning task with support vector regression to classify VOIs as healthy or osteoarthritic; classification performance was evaluated using the area under the receiver-operating characteristic (ROC) curve (AUC). Our results show that the classification performance achieved by 9-D SIM-derived geometric feature sets (AUC: 0.96 ± 0.02) can be maintained with 2-D representations computed from both dimension reduction and feature selection (AUC values as high as 0.97 ± 0.02). Thus, such feature reduction techniques can offer a high degree of compaction to large feature sets extracted from PCI-CT images while maintaining their ability to characterize the underlying chondrocyte patterns.

Volumetric Characterization of Human Patellar Cartilage Matrix on Phase Contrast X-Ray Computed Tomography.

Phase contrast X-ray computed tomography (PCI-CT) has recently emerged as a novel imaging technique that allows visualization of cartilage soft tissue, subsequent examination of chondrocyte patterns, and their correlation to osteoarthritis. Previous studies have shown that 2D texture features are effective at distinguishing between healthy and osteoarthritic regions of interest annotated in the radial zone of cartilage matrix on PCI-CT images. In this study, we further extend the texture analysis to 3D and investigate the ability of volumetric texture features at characterizing chondrocyte patterns in the cartilage matrix for purposes of classification. Here, we extracted volumetric texture features derived from Minkowski Functionals and gray-level co-occurrence matrices (GLCM) from 496 volumes of interest (VOI) annotated on PCI-CT images of human patellar cartilage specimens. The extracted features were then used in a machine-learning task involving support vector regression to classify ROIs as healthy or osteoarthritic. Classification performance was evaluated using the area under the receiver operating characteristic (ROC) curve (AUC). The best classification performance was observed with GLCM features correlation (AUC = 0.83 ± 0.06) and homogeneity (AUC = 0.82 ± 0.07), which significantly outperformed all Minkowski Functionals (p < 0.05). These results suggest that such quantitative analysis of chondrocyte patterns in human patellar cartilage matrix involving GLCM-derived statistical features can distinguish between healthy and osteoarthritic tissue with high accuracy.

Functional Connectivity Analysis in Resting State fMRI with Echo-State Networks and Non-Metric Clustering for Network Structure Recovery.

Echo state networks (ESN) are recurrent neural networks where the hidden layer is replaced with a fixed reservoir of neurons. Unlike feed-forward networks, neuron training in ESN is restricted to the output neurons alone thereby providing a computational advantage. We demonstrate the use of such ESNs in our mutual connectivity analysis (MCA) framework for recovering the primary motor cortex network associated with hand movement from resting state functional MRI (fMRI) data. Such a framework consists of two steps - (1) defining a pair-wise affinity matrix between different pixel time series within the brain to characterize network activity and (2) recovering network components from the affinity matrix with non-metric clustering. Here, ESNs are used to evaluate pair-wise cross-estimation performance between pixel time series to create the affinity matrix, which is subsequently subject to non-metric clustering with the Louvain method. For comparison, the ground truth of the motor cortex network structure is established with a task-based fMRI sequence. Overlap between the primary motor cortex network recovered with our model free MCA approach and the ground truth was measured with the Dice coefficient. Our results show that network recovery with our proposed MCA approach is in close agreement with the ground truth. Such network recovery is achieved without requiring low-pass filtering of the time series ensembles prior to analysis, an fMRI preprocessing step that has courted controversy in recent years. Thus, we conclude our MCA framework can allow recovery and visualization of the underlying functionally connected networks in the brain on resting state fMRI.

Characterizing Trabecular Bone structure for Assessing Vertebral Fracture Risk on Volumetric Quantitative Computed Tomography.

While the proximal femur is preferred for measuring bone mineral density (BMD) in fracture risk estimation, the introduction of volumetric quantitative computed tomography has revealed stronger associations between BMD and spinal fracture status. In this study, we propose to capture properties of trabecular bone structure in spinal vertebrae with advanced second-order statistical features for purposes of fracture risk assessment. For this purpose, axial multi-detector CT (MDCT) images were acquired from 28 spinal vertebrae specimens using a whole-body 256-row CT scanner with a dedicated calibration phantom. A semi-automated method was used to annotate the trabecular compartment in the central vertebral slice with a circular region of interest (ROI) to exclude cortical bone; pixels within were converted to values indicative of BMD. Six second-order statistical features derived from gray-level co-occurrence matrices (GLCM) and the mean BMD within the ROI were then extracted and used in conjunction with a generalized radial basis functions (GRBF) neural network to predict the failure load of the specimens; true failure load was measured through biomechanical testing. Prediction performance was evaluated with a root-mean-square error (RMSE) metric. The best prediction performance was observed with GLCM feature 'correlation' (RMSE = 1.02 ± 0.18), which significantly outperformed all other GLCM features (p < 0.01). GLCM feature correlation also significantly outperformed MDCT-measured mean BMD (RMSE = 1.11 ± 0.17) (p < 10-4). These results suggest that biomechanical strength prediction in spinal vertebrae can be significantly improved through characterization of trabecular bone structure with GLCM-derived texture features.

Investigating the use of mutual information and non-metric clustering for functional connectivity analysis on resting-state functional MRI.

Functional MRI (fMRI) is currently used to investigate structural and functional connectivity in human brain networks. To this end, previous studies have proposed computational methods that involve assumptions that can induce information loss, such as assumed linear coupling of the fMRI signals or requiring dimension reduction. This study presents a new computational framework for investigating the functional connectivity in the brain and recovering network structure while reducing the information loss inherent in previous methods. For this purpose, pair-wise mutual information (MI) was extracted from all pixel time series within the brain on resting-state fMRI data. Non-metric topographic mapping of proximity (TMP) data was subsequently applied to recover network structure from the pair-wise MI analysis. Our computational framework is demonstrated in the task of identifying regions of the primary motor cortex network on resting state fMRI data. For ground truth comparison, we also localized regions of the primary motor cortex associated with hand movement in a task-based fMRI sequence with a finger-tapping stimulus function. The similarity between our pair-wise MI clustering results and the ground truth is evaluated using the dice coefficient. Our results show that non-metric clustering with the TMP algorithm, as performed on pair-wise MI analysis, was able to detect the primary motor cortex network and achieved a dice coefficient of 0.53 in terms of overlap with the ground truth. Thus, we conclude that our computational framework can extract and visualize valuable information concerning the underlying network structure between different regions of the brain in resting state fMRI.

Nonlinear Functional Connectivity Network Recovery in the Human Brain with Mutual Connectivity Analysis (MCA): Convergent Cross-Mapping and Non-Metric Clustering.

We explore a computational framework for functional connectivity analysis in resting-state functional MRI (fMRI) data acquired from the human brain for recovering the underlying network structure and understanding causality between network components. Termed mutual connectivity analysis (MCA), this framework involves two steps, the first of which is to evaluate the pair-wise cross-prediction performance between fMRI pixel time series within the brain. In a second step, the underlying network structure is subsequently recovered from the affinity matrix using non-metric network clustering approaches, such as the so-called Louvain method. Finally, we use convergent cross-mapping (CCM) to study causality between different network components. We demonstrate our MCA framework in the problem of recovering the motor cortex network associated with hand movement from resting state fMRI data. Results are compared with a ground truth of active motor cortex regions as identified by a task-based fMRI sequence involving a finger-tapping stimulation experiment. Our results regarding causation between regions of the motor cortex revealed a significant directional variability and were not readily interpretable in a consistent manner across subjects. However, our results on whole-slice fMRI analysis demonstrate that MCA-based model-free recovery of regions associated with the primary motor cortex and supplementary motor area are in close agreement with localization of similar regions achieved with a task-based fMRI acquisition. Thus, we conclude that our MCA methodology can extract and visualize valuable information concerning the underlying network structure between different regions of the brain in resting state fMRI.