Neuroconductor: an R platform for medical imaging analysis

Neuroconductor (https://neuroconductor.org) is an open-source platform for rapid testing and dissemination of reproducible computational imaging software. The goals of the project are to: (i) provide a centralized repository of R software dedicated to image analysis, (ii) disseminate software updates quickly, (iii) train a large, diverse community of scientists using detailed tutorials and short courses, (iv) increase software quality via automatic and manual quality controls, and (v) promote reproducibility of image data analysis. Based on the programming language R (https://www.r-project.org/), Neuroconductor starts with 51 inter-operable packages that cover multiple areas of imaging including visualization, data processing and storage, and statistical inference. Neuroconductor accepts new R package submissions, which are subject to a formal review and continuous automated testing. We provide a description of the purpose of Neuroconductor and the user and developer experience

Sub‐phenotyping Metabolic Disorders Using Body Composition: An Individualized, Nonparametric Approach Utilizing Large Data Sets

Objective This study performed individual-centric, data-driven calculations of propensity for coronary heart disease (CHD) and type 2 diabetes (T2D), utilizing magnetic resonance imaging-acquired body composition measurements, for sub-phenotyping of obesity and nonalcoholic fatty liver disease (NAFLD). Methods A total of 10,019 participants from the UK Biobank imaging substudy were included and analyzed for visceral and abdominal subcutaneous adipose tissue, muscle fat infiltration, and liver fat. An adaption of the k-nearest neighbors algorithm was applied to the imaging variable space to calculate individualized CHD and T2D propensity and explore metabolic sub-phenotyping within obesity and NAFLD. Results The ranges of CHD and T2D propensity for the whole cohort were 1.3% to 58.0% and 0.6% to 42.0%, respectively. The diagnostic performance, area under the receiver operating characteristic curve (95% CI), using disease propensities for CHD and T2D detection was 0.75 (0.73-0.77) and 0.79 (0.77-0.81). Exploring individualized disease propensity, CHD phenotypes, T2D phenotypes, comorbid phenotypes, and metabolically healthy phenotypes were found within obesity and NAFLD. Conclusions The adaptive k-nearest neighbors algorithm allowed an individual-centric assessment of each individual’s metabolic phenotype moving beyond discrete categorizations of body composition. Within obesity and NAFLD, this may help in identifying which comorbidities a patient may develop and consequently enable optimization of treatment

A Multiresolution Census Algorithm for Calculating Vortex Statistics in Turbulent Flows

The fundamental equations that model turbulent flow do not provide much insight into the size and shape of observed turbulent structures. We investigate the efficient and accurate representation of structures in two-dimensional turbulence by applying statistical models directly to the simulated vorticity field. Rather than extract the coherent portion of the image from the background variation, as in the classical signal-plus-noise model, we present a model for individual vortices using the non-decimated discrete wavelet transform. A template image, supplied by the user, provides the features to be extracted from the vorticity field. By transforming the vortex template into the wavelet domain, specific characteristics present in the template, such as size and symmetry, are broken down into components associated with spatial frequencies. Multivariate multiple linear regression is used to fit the vortex template to the vorticity field in the wavelet domain. Since all levels of the template decomposition may be used to model each level in the field decomposition, the resulting model need not be identical to the template. Application to a vortex census algorithm that records quantities of interest (such as size, peak amplitude, circulation, etc.) as the vorticity field evolves is given. The multiresolution census algorithm extracts coherent structures of all shapes and sizes in simulated vorticity fields and is able to reproduce known physical scaling laws when processing a set of voriticity fields that evolve over time

Large‑scale analysis of iliopsoas muscle volumes in the UK Biobank

Psoas muscle measurements are frequently used as markers of sarcopenia and predictors of health. Manually measured cross-sectional areas are most commonly used, but there is a lack of consistency regarding the position of the measurement and manual annotations are not practical for large population studies. We have developed a fully automated method to measure iliopsoas muscle volume (comprised of the psoas and iliacus muscles) using a convolutional neural network. Magnetic resonance images were obtained from the UK Biobank for 5000 participants, balanced for age, gender and BMI. Ninety manual annotations were available for model training and validation. The model showed excellent performance against out-of-sample data (average dice score coefficient of 0.9046 ± 0.0058 for six-fold cross-validation). Iliopsoas muscle volumes were successfully measured in all 5000 participants. Iliopsoas volume was greater in male compared with female subjects. There was a small but significant asymmetry between left and right iliopsoas muscle volumes. We also found that iliopsoas volume was significantly related to height, BMI and age, and that there was an acceleration in muscle volume decrease in men with age. Our method provides a robust technique for measuring iliopsoas muscle volume that can be applied to large cohorts

Traditional eye medicine use in microbial keratitis in Uganda : a mixed methods study [version 2; peer review: 2 approved]

Background: Traditional eye medicine (TEM) is frequently used to treat microbial keratitis (MK) in many parts of Africa. Few reports have suggested that this is associated with a worse outcome. We undertook this large prospective study to determine how TEM use impacts presentation and outcome of MK and to explore reasons why people use TEM for treatment in Uganda. Methods: In a mixed method prospective cohort study, we enrolled patients presenting with MK at the two main eye units in Southern Uganda between December 2016 and March 2018 and collected information on history, TEM use, microbiology and 3-month outcomes. We conducted qualitative interviews with patients, carers traditional healers on reasons why people use TEM. Outcome measures included presenting vision and at 3-months, comparing TEM Users versus Non-Users. A thematic coding framework was deployed to explore reasons for use of TEM. Results: Out of 313 participants enrolled, 188 reported TEM use. TEM Users had a delayed presentation; median presenting time 18 days versus 14 days, p= 0.005; had larger ulcers 5.6 mm versus 4.3 mm p=0.0005; a worse presenting visual acuity median logarithm of the minimum angle of resolution (Log MAR) 1.5 versus 0.6, p=0.005; and, a worse visual acuity at 3 months median Log MAR 0.6 versus 0.2, p=0.010. In a multivariable logistic regression model, distance from the eye hospital and delayed presentation were associated with TEM use. Reasons for TEM use included lack of confidence in conventional medicine, health system breakdown, poverty, fear of the eye hospital, cultural belief in TEM, influence from traditional healers, personal circumstances and ignorance. Conclusion: TEM users had poorer clinical presentation and outcomes. Capacity building of the primary health centres to improve access to eye care and community behavioural change initiatives against TEM use should be encouraged

Normalized indices derived from visceral adipose mass assessed by MRI and their correlation with markers for insulin resistance and prediabetes

Visceral adipose tissue (VAT) plays an important role in the pathogenesis of insulin resistance (IR), prediabetes and type 2 diabetes. However, VAT volume alone might not be the best marker for insulin resistance and prediabetes or diabetes, as a given VAT volume may differently impact on these metabolic traits based on body height, gender, age and ethnicity. In a cohort of 1295 subjects from the Tübingen Diabetes Family Study (TDFS) and in 9978 subjects from the UK Biobank (UKBB), undergoing magnetic resonance imaging for quantification of VAT volume, total adipose tissue (TAT, in the TDFS), total abdominal adipose tissue (TAAT) in the UKBB, and total lean tissue (TLT), VAT volume and several VAT-indices were investigated for their relationships with insulin resistance and glycemic traits. VAT-related indices were calculated by correcting for body height (VAT/m: VAT/body height; VAT/m²: VAT/(body height)², and VAT/m³: VAT/(body height)³), TAT (%VAT), TLT (VAT/TLT) and weight (VAT/WEI), with closest equivalents used within the UKBB dataset. Prognostic values of VAT and VAT-related indices for insulin sensitivity, HbA1c levels and prediabetes/diabetes were analyzed for males and females. Males had higher VAT volume and VAT-related indices than females in both cohorts (p < 0.0001) and VAT volume has shown to be a stronger determinant for insulin sensitivity than anthropometric variables. Among the parameters uncorrected VAT and derived indices, VAT/m³ most strongly correlated negatively with insulin sensitivity and positively with HbA1c levels and prediabetes/diabetes in the TDFS (R² = 0.375/0.305 for females/males for insulin sensitivity, 0.178/0.148 for HbA1c levels vs. – e.g. – 0.355/0.293 and 0.144/0.133 for VAT, respectively) and positively with HbA1c (R² = 0.046/0.042) in the UKBB for females and males. Furthermore, VAT/m³ was found to be a significantly better determinant of insulin resistance or prediabetes than uncorrected VAT volume (p < 0.001/0.019 for females/males regarding insulin sensitivity, p < 0.001/< 0.001 for females/males regarding HbA1c). Evaluation of several indices derived from VAT volume identified VAT/m³ to most strongly correlate with insulin sensitivity and glucose metabolism. Thus, VAT/m³ appears to provide better indications of metabolic characteristics (insulin sensitivity and pre-diabetes/diabetes) than VAT volume alone

Verbal working memory and functional large-scale networks in schizophrenia

The aim of this study was to test whether bilinear and nonlinear effective connectivity (EC) measures of working memory fMRI data can differentiate between patients with schizophrenia (SZ) and healthy controls (HC). We applied bilinear and nonlinear Dynamic Causal Modeling (DCM) for the analysis of verbal working memory in 16 SZ and 21 HC. The connection strengths with nonlinear modulation between the dorsolateral prefrontal cortex (DLPFC) and the ventral tegmental area/substantia nigra (VTA/SN) were evaluated. We used Bayesian Model Selection at the group and family levels to compare the optimal bilinear and nonlinear models. Bayesian Model Averaging was used to assess the connection strengths with nonlinear modulation. The DCM analyses revealed that SZ and HC used different bilinear networks despite comparable behavioral performance. In addition, the connection strengths with nonlinear modulation between the DLPFC and the VTA/SN area showed differences between SZ and HC. The adoption of different functional networks in SZ and HC indicated neurobiological alterations underlying working memory performance, including different connection strengths with nonlinear modulation between the DLPFC and the VTA/SN area. These novel findings may increase our understanding of connectivity in working memory in schizophrenia

Neuroconductor: an R platform for medical imaging analysis

Neuroconductor (https://neuroconductor.org) is an open-source platform for rapid testing and dissemination of reproducible computational imaging software. The goals of the project are to: (i) provide a centralized repository of R software dedicated to image analysis, (ii) disseminate software updates quickly, (iii) train a large, diverse community of scientists using detailed tutorials and short courses, (iv) increase software quality via automatic and manual quality controls, and (v) promote reproducibility of image data analysis. Based on the programming language R (https://www.r-project.org/), Neuroconductor starts with 51 inter-operable packages that cover multiple areas of imaging including visualization, data processing and storage, and statistical inference. Neuroconductor accepts new R package submissions, which are subject to a formal review and continuous automated testing. We provide a description of the purpose of Neuroconductor and the user and developer experience

Genetic evidence for different adiposity phenotypes and their opposing influence on ectopic fat and risk of cardiometabolic disease

To understand the causal role of adiposity and ectopic fat in type 2 diabetes and cardiometabolic diseases, we aimed to identify two clusters of adiposity genetic variants, one with ‘adverse’ metabolic effects (UFA) and the other with, paradoxically, ‘favourable’ metabolic effects (FA). We performed a multivariate genome-wide association study using body fat percentage and metabolic biomarkers from UK Biobank and identified 38 UFA and 36 FA variants. Adiposity-increasing alleles were associated with an adverse metabolic profile, higher risk of disease, higher CRP, higher fat in subcutaneous and visceral adipose tissue, liver and pancreas for UFA; and a favourable metabolic profile, lower risk of disease, higher CRP, higher subcutaneous adipose tissue but lower liver fat for FA. We detected no sexual dimorphism. The Mendelian randomization studies provided evidence for risk-increasing effect of UFA and protective effect of FA on type 2 diabetes, heart disease, hypertension, stroke, non-alcoholic fatty liver disease and polycystic ovary syndrome. FA is distinct from UFA by its association with lower liver fat, and protection from cardiometabolic diseases; it was not associated with visceral or pancreatic fat. Understanding the difference in FA and UFA may lead to new insights in preventing, predicting and treating of cardiometabolic diseases

Effects of a balanced translocation between chromosomes 1 and 11 disrupting the DISC1 locus on white matter integrity

Objective Individuals carrying rare, but biologically informative genetic variants provide a unique opportunity to model major mental illness and inform understanding of disease mechanisms. The rarity of such variations means that their study involves small group numbers, however they are amongst the strongest known genetic risk factors for major mental illness and are likely to have large neural effects. DISC1 (Disrupted in Schizophrenia 1) is a gene containing one such risk variant, identified in a single Scottish family through its disruption by a balanced translocation of chromosomes 1 and 11; t(1;11) (q42.1;q14.3). Method Within the original pedigree, we examined the effects of the t(1;11) translocation on white matter integrity, measured by fractional anisotropy (FA). This included family members with (n = 7) and without (n = 13) the translocation, along with a clinical control sample of patients with psychosis (n = 34), and a group of healthy controls (n = 33). Results We report decreased white matter integrity in five clusters in the genu of the corpus callosum, the right inferior fronto-occipital fasciculus, acoustic radiation and fornix. Analysis of the mixed psychosis group also demonstrated decreased white matter integrity in the above regions. FA values within the corpus callosum correlated significantly with positive psychotic symptom severity. Conclusions We demonstrate that the t(1;11) translocation is associated with reduced white matter integrity in frontal commissural and association fibre tracts. These findings overlap with those shown in affected patients with psychosis and in DISC1 animal models and highlight the value of rare but biologically informative mutations in modeling psychosis