Multiorgan MRI findings after hospitalisation with COVID-19 in the UK (C-MORE): a prospective, multicentre, observational cohort study

INTRODUCTION

The multiorgan impact of moderate to severe coronavirus infections in the post-acute phase is still poorly understood. We aimed to evaluate the excess burden of multiorgan abnormalities after hospitalisation with COVID-19, evaluate their determinants, and explore associations with patient-related outcome measures.

METHODS

In a prospective, UK-wide, multicentre MRI follow-up study (C-MORE), adults (aged ≥18 years) discharged from hospital following COVID-19 who were included in Tier 2 of the Post-hospitalisation COVID-19 study (PHOSP-COVID) and contemporary controls with no evidence of previous COVID-19 (SARS-CoV-2 nucleocapsid antibody negative) underwent multiorgan MRI (lungs, heart, brain, liver, and kidneys) with quantitative and qualitative assessment of images and clinical adjudication when relevant. Individuals with end-stage renal failure or contraindications to MRI were excluded. Participants also underwent detailed recording of symptoms, and physiological and biochemical tests. The primary outcome was the excess burden of multiorgan abnormalities (two or more organs) relative to controls, with further adjustments for potential confounders. The C-MORE study is ongoing and is registered with ClinicalTrials.gov, NCT04510025.

FINDINGS

Of 2710 participants in Tier 2 of PHOSP-COVID, 531 were recruited across 13 UK-wide C-MORE sites. After exclusions, 259 C-MORE patients (mean age 57 years [SD 12]; 158 [61%] male and 101 [39%] female) who were discharged from hospital with PCR-confirmed or clinically diagnosed COVID-19 between March 1, 2020, and Nov 1, 2021, and 52 non-COVID-19 controls from the community (mean age 49 years [SD 14]; 30 [58%] male and 22 [42%] female) were included in the analysis. Patients were assessed at a median of 5·0 months (IQR 4·2-6·3) after hospital discharge. Compared with non-COVID-19 controls, patients were older, living with more obesity, and had more comorbidities. Multiorgan abnormalities on MRI were more frequent in patients than in controls (157 [61%] of 259 vs 14 [27%] of 52; p<0·0001) and independently associated with COVID-19 status (odds ratio [OR] 2·9 [95% CI 1·5-5·8]; padjusted=0·0023) after adjusting for relevant confounders. Compared with controls, patients were more likely to have MRI evidence of lung abnormalities (p=0·0001; parenchymal abnormalities), brain abnormalities (p<0·0001; more white matter hyperintensities and regional brain volume reduction), and kidney abnormalities (p=0·014; lower medullary T1 and loss of corticomedullary differentiation), whereas cardiac and liver MRI abnormalities were similar between patients and controls. Patients with multiorgan abnormalities were older (difference in mean age 7 years [95% CI 4-10]; mean age of 59·8 years [SD 11·7] with multiorgan abnormalities vs mean age of 52·8 years [11·9] without multiorgan abnormalities; p<0·0001), more likely to have three or more comorbidities (OR 2·47 [1·32-4·82]; padjusted=0·0059), and more likely to have a more severe acute infection (acute CRP >5mg/L, OR 3·55 [1·23-11·88]; padjusted=0·025) than those without multiorgan abnormalities. Presence of lung MRI abnormalities was associated with a two-fold higher risk of chest tightness, and multiorgan MRI abnormalities were associated with severe and very severe persistent physical and mental health impairment (PHOSP-COVID symptom clusters) after hospitalisation.

INTERPRETATION

After hospitalisation for COVID-19, people are at risk of multiorgan abnormalities in the medium term. Our findings emphasise the need for proactive multidisciplinary care pathways, with the potential for imaging to guide surveillance frequency and therapeutic stratification.

FUNDING

UK Research and Innovation and National Institute for Health Research.

Automated detection of cerebral microbleeds on MR images using knowledge distillation framework

Introduction

Cerebral microbleeds (CMBs) are associated with white matter damage, and various neurodegenerative and cerebrovascular diseases. CMBs occur as small, circular hypointense lesions on T2*-weighted gradient recalled echo (GRE) and susceptibility-weighted imaging (SWI) images, and hyperintense on quantitative susceptibility mapping (QSM) images due to their paramagnetic nature. Accurate automated detection of CMBs would help to determine quantitative imaging biomarkers (e.g., CMB count) on large datasets. In this work, we propose a fully automated, deep learning-based, 3-step algorithm, using structural and anatomical properties of CMBs from any single input image modality (e.g., GRE/SWI/QSM) for their accurate detections.

Methods

In our method, the first step consists of an initial candidate detection step that detects CMBs with high sensitivity. In the second step, candidate discrimination step is performed using a knowledge distillation framework, with a multi-tasking teacher network that guides the student network to classify CMB and non-CMB instances in an offline manner. Finally, a morphological clean-up step further reduces false positives using anatomical constraints. We used four datasets consisting of different modalities specified above, acquired using various protocols and with a variety of pathological and demographic characteristics.

Results

On cross-validation within datasets, our method achieved a cluster-wise true positive rate (TPR) of over 90% with an average of <2 false positives per subject. The knowledge distillation framework improves the cluster-wise TPR of the student model by 15%. Our method is flexible in terms of the input modality and provides comparable cluster-wise TPR and better cluster-wise precision compared to existing state-of-the-art methods. When evaluating across different datasets, our method showed good generalizability with a cluster-wise TPR >80 % with different modalities. The python implementation of the proposed method is openly available.

Association of gout with brain reserve and vulnerability to neurodegenerative disease

Studies of neurodegenerative disease risk in gout are contradictory. Relationships with neuroimaging markers of brain structure, which may offer insights, are uncertain. Here we investigated associations between gout, brain structure, and neurodegenerative disease incidence. Gout patients had smaller global and regional brain volumes and markers of higher brain iron, using both observational and genetic approaches. Participants with gout also had higher incidence of all-cause dementia, Parkinson's disease, and probable essential tremor. Risks were strongly time dependent, whereby associations with incident dementia were highest in the first 3 years after gout diagnosis. These findings suggest gout is causally related to several measures of brain structure. Lower brain reserve amongst gout patients may explain their higher vulnerability to multiple neurodegenerative diseases. Motor and cognitive impairments may affect gout patients, particularly in early years after diagnosis.

ICAM-Reg: Interpretable Classification and Regression With Feature Attribution for Mapping Neurological Phenotypes in Individual Scans

An important goal of medical imaging is to be able to precisely detect patterns of disease specific to individual scans; however, this is challenged in brain imaging by the degree of heterogeneity of shape and appearance. Traditional methods, based on image registration, historically fail to detect variable features of disease, as they utilise population-based analyses, suited primarily to studying group-average effects. In this paper we therefore take advantage of recent developments in generative deep learning to develop a method for simultaneous classification, or regression, and feature attribution (FA). Specifically, we explore the use of a VAE-GAN (variational autoencoder - general adversarial network) for translation called ICAM, to explicitly disentangle class relevant features, from background confounds, for improved interpretability and regression of neurological phenotypes. We validate our method on the tasks of Mini-Mental State Examination (MMSE) cognitive test score prediction for the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort, as well as brain age prediction, for both neurodevelopment and neurodegeneration, using the developing Human Connectome Project (dHCP) and UK Biobank datasets. We show that the generated FA maps can be used to explain outlier predictions and demonstrate that the inclusion of a regression module improves the disentanglement of the latent space. Our code is freely available on GitHub https://github.com/CherBass/ICAM.

Telomere length and brain imaging phenotypes in UK Biobank

Telomeres form protective caps at the ends of chromosomes, and their attrition is a marker of biological aging. Short telomeres are associated with an increased risk of neurological and psychiatric disorders including dementia. The mechanism underlying this risk is unclear, and may involve brain structure and function. However, the relationship between telomere length and neuroimaging markers is poorly characterized. Here we show that leucocyte telomere length (LTL) is associated with multi-modal MRI phenotypes in 31,661 UK Biobank participants. Longer LTL is associated with: i) larger global and subcortical grey matter volumes including the hippocampus, ii) lower T1-weighted grey-white tissue contrast in sensory cortices, iii) white-matter microstructure measures in corpus callosum and association fibres, iv) lower volume of white matter hyperintensities, and v) lower basal ganglia iron. Longer LTL was protective against certain related clinical manifestations, namely all-cause dementia (HR 0.93, 95% CI: 0.91-0.96), but not stroke or Parkinson's disease. LTL is associated with multiple MRI endophenotypes of neurodegenerative disease, suggesting a pathway by which longer LTL may confer protective against dementia.

Phenotypic and genetic associations of quantitative magnetic susceptibility in UK Biobank brain imaging

A key aim in epidemiological neuroscience is identification of markers to assess brain health and monitor therapeutic interventions. Quantitative susceptibility mapping (QSM) is an emerging magnetic resonance imaging technique that measures tissue magnetic susceptibility and has been shown to detect pathological changes in tissue iron, myelin and calcification. We present an open resource of QSM-based imaging measures of multiple brain structures in 35,273 individuals from the UK Biobank prospective epidemiological study. We identify statistically significant associations of 251 phenotypes with magnetic susceptibility that include body iron, disease, diet and alcohol consumption. Genome-wide associations relate magnetic susceptibility to 76 replicating clusters of genetic variants with biological functions involving iron, calcium, myelin and extracellular matrix. These patterns of associations include relationships that are unique to QSM, in particular being complementary to T2* signal decay time measures. These new imaging phenotypes are being integrated into the core UK Biobank measures provided to researchers worldwide, creating the potential to discover new, non-invasive markers of brain health.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence of brain-related abnormalities in COVID-191–13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51–81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans—with 141 days on average separating their diagnosis and the second scan—as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence of brain-related abnormalities in COVID-191-13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51-81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans-with 141 days on average separating their diagnosis and the second scan-as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence for brain-related abnormalities in COVID-19 1-13 . It remains unknown however whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here, we investigated brain changes in 785 UK Biobank participants (aged 51-81) imaged twice, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans, with 141 days on average separating their diagnosis and second scan, and 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including: (i) greater reduction in grey matter thickness and tissue-contrast in the orbitofrontal cortex and parahippocampal gyrus, (ii) greater changes in markers of tissue damage in regions functionally-connected to the primary olfactory cortex, and (iii) greater reduction in global brain size. The infected participants also showed on average larger cognitive decline between the two timepoints. Importantly, these imaging and cognitive longitudinal effects were still seen after excluding the 15 cases who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease via olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious impact can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow up.

Automated Detection of Candidate Subjects With Cerebral Microbleeds Using Machine Learning

Cerebral microbleeds (CMBs) appear as small, circular, well defined hypointense lesions of a few mm in size on T2^*-weighted gradient recalled echo (T2^*-GRE) images and appear enhanced on susceptibility weighted images (SWI). Due to their small size, contrast variations and other mimics (e.g., blood vessels), CMBs are highly challenging to detect automatically. In large datasets (e.g., the UK Biobank dataset), exhaustively labelling CMBs manually is difficult and time consuming. Hence it would be useful to preselect candidate CMB subjects in order to focus on those for manual labelling, which is essential for training and testing automated CMB detection tools on these datasets. In this work, we aim to detect CMB candidate subjects from a larger dataset, UK Biobank, using a machine learning-based, computationally light pipeline. For our evaluation, we used 3 different datasets, with different intensity characteristics, acquired with different scanners. They include the UK Biobank dataset and two clinical datasets with different pathological conditions. We developed and evaluated our pipelines on different types of images, consisting of SWI or GRE images. We also used the UK Biobank dataset to compare our approach with alternative CMB preselection methods using non-imaging factors and/or imaging data. Finally, we evaluated the pipeline's generalisability across datasets. Our method provided subject-level detection accuracy > 80% on all the datasets (within-dataset results), and showed good generalisability across datasets, providing a consistent accuracy of over 80%, even when evaluated across different modalities.

Adapting UK Biobank imaging for use in a routine memory clinic setting: The Oxford Brain Health Clinic

The Oxford Brain Health Clinic (BHC) is a joint clinical-research service that provides memory clinic patients and clinicians access to high-quality assessments not routinely available, including brain MRI aligned with the UK Biobank imaging study (UKB). In this work we present how we 1) adapted the UKB MRI acquisition protocol to be suitable for memory clinic patients, 2) modified the imaging analysis pipeline to extract measures that are in line with radiology reports and 3) explored the alignment of measures from BHC patients to the largest brain MRI study in the world (ultimately 100,000 participants). Adaptations of the UKB acquisition protocol for BHC patients include dividing the scan into core and optional sequences (i.e., additional imaging modalities) to improve patients' tolerance for the MRI assessment. We adapted the UKB structural MRI analysis pipeline to take into account the characteristics of a memory clinic population (e.g., high amount of white matter hyperintensities and hippocampal atrophy). We then compared the imaging derived phenotypes (IDPs) extracted from the structural scans to visual ratings from radiology reports, non-imaging factors (age, cognition) and to reference distributions derived from UKB data. Of the first 108 BHC attendees (August 2020-November 2021), 92.5 % completed the clinical scans, 88.0 % consented to use of data for research, and 43.5 % completed the additional research sequences, demonstrating that the protocol is well tolerated. The high rates of consent to research makes this a valuable real-world quality research dataset routinely captured in a clinical service. Modified tissue-type segmentation with lesion masking greatly improved grey matter volume estimation. CSF-masking marginally improved hippocampal segmentation. The IDPs were in line with radiology reports and showed significant associations with age and cognitive performance, in line with the literature. Due to the age difference between memory clinic patients of the BHC (age range 65-101 years, average 78.3 years) and UKB participants (44-82 years, average 64 years), additional scans on elderly healthy controls are needed to improve reference distributions. Current and future work aims to integrate automated quantitative measures in the radiology reports and evaluate their clinical utility.

Adapting the UK Biobank Brain Imaging Protocol and Analysis Pipeline for the C-MORE Multi-Organ Study of COVID-19 Survivors

SARS-CoV-2 infection has been shown to damage multiple organs, including the brain. Multiorgan MRI can provide further insight on the repercussions of COVID-19 on organ health but requires a balance between richness and quality of data acquisition and total scan duration. We adapted the UK Biobank brain MRI protocol to produce high-quality images while being suitable as part of a post-COVID-19 multiorgan MRI exam. The analysis pipeline, also adapted from UK Biobank, includes new imaging-derived phenotypes (IDPs) designed to assess the possible effects of COVID-19. A first application of the protocol and pipeline was performed in 51 COVID-19 patients post-hospital discharge and 25 controls participating in the Oxford C-MORE study. The protocol acquires high resolution T1, T2-FLAIR, diffusion weighted images, susceptibility weighted images, and arterial spin labelling data in 17 min. The automated imaging pipeline derives 1,575 IDPs, assessing brain anatomy (including olfactory bulb volume and intensity) and tissue perfusion, hyperintensities, diffusivity, and susceptibility. In the C-MORE data, IDPs related to atrophy, small vessel disease and olfactory bulbs were consistent with clinical radiology reports. Our exploratory analysis tentatively revealed some group differences between recovered COVID-19 patients and controls, across severity groups, but not across anosmia groups. Follow-up imaging in the C-MORE study is currently ongoing, and this protocol is now being used in other large-scale studies. The protocol, pipeline code and data are openly available and will further contribute to the understanding of the medium to long-term effects of COVID-19.

An expanded set of genome-wide association studies of brain imaging phenotypes in UK Biobank

UK Biobank is a major prospective epidemiological study, including multimodal brain imaging, genetics and ongoing health outcomes. Previously, we published genome-wide associations of 3,144 brain imaging-derived phenotypes, with a discovery sample of 8,428 individuals. Here we present a new open resource of genome-wide association study summary statistics, using the 2020 data release, almost tripling the discovery sample size. We now include the X chromosome and new classes of imaging-derived phenotypes (subcortical volumes and tissue contrast). Previously, we found 148 replicated clusters of associations between genetic variants and imaging phenotypes; in this study, we found 692, including 12 on the X chromosome. We describe some of the newly found associations, focusing on the X chromosome and autosomal associations involving the new classes of imaging-derived phenotypes. Our novel associations implicate, for example, pathways involved in the rare X-linked STAR (syndactyly, telecanthus and anogenital and renal malformations) syndrome, Alzheimer's disease and mitochondrial disorders.

Medium-term effects of SARS-CoV-2 infection on multiple vital organs, exercise capacity, cognition, quality of life and mental health, post-hospital discharge

BACKGROUND

The medium-term effects of Coronavirus disease (COVID-19) on organ health, exercise capacity, cognition, quality of life and mental health are poorly understood.

METHODS

Fifty-eight COVID-19 patients post-hospital discharge and 30 age, sex, body mass index comorbidity-matched controls were enrolled for multiorgan (brain, lungs, heart, liver and kidneys) magnetic resonance imaging (MRI), spirometry, six-minute walk test, cardiopulmonary exercise test (CPET), quality of life, cognitive and mental health assessments.

FINDINGS

At 2-3 months from disease-onset, 64% of patients experienced breathlessness and 55% reported fatigue. On MRI, abnormalities were seen in lungs (60%), heart (26%), liver (10%) and kidneys (29%). Patients exhibited changes in the thalamus, posterior thalamic radiations and sagittal stratum on brain MRI and demonstrated impaired cognitive performance, specifically in the executive and visuospatial domains. Exercise tolerance (maximal oxygen consumption and ventilatory efficiency on CPET) and six-minute walk distance were significantly reduced. The extent of extra-pulmonary MRI abnormalities and exercise intolerance correlated with serum markers of inflammation and acute illness severity. Patients had a higher burden of self-reported symptoms of depression and experienced significant impairment in all domains of quality of life compared to controls (p<0.0001 to 0.044).

INTERPRETATION

A significant proportion of patients discharged from hospital reported symptoms of breathlessness, fatigue, depression and had limited exercise capacity. Persistent lung and extra-pulmonary organ MRI findings are common in patients and linked to inflammation and severity of acute illness.

FUNDING

NIHR Oxford and Oxford Health Biomedical Research Centres, British Heart Foundation Centre for Research Excellence, UKRI, Wellcome Trust, British Heart Foundation.

The human hippocampus and its subfield volumes across age, sex and APOE e4 status

Female sex, age and carriage of the apolipoprotein E e4 allele are the greatest risk factors for sporadic Alzheimer's disease. The hippocampus has a selective vulnerability to atrophy in ageing that may be accelerated in Alzheimer's disease, including in those with increased genetic risk of the disease, years before onset. Within the hippocampal complex, subfields represent cytoarchitectonic and connectivity based divisions. Variation in global hippocampal and subfield volume associated with sex, age and apolipoprotein E e4 status has the potential to provide a sensitive biomarker of future vulnerability to Alzheimer's disease. Here, we examined non-linear age, sex and apolipoprotein E effects, and their interactions, on hippocampal and subfield volumes across several decades spanning mid-life to old age in 36 653 healthy ageing individuals. FMRIB Software Library derived estimates of total hippocampal volume and Freesurfer derived estimates hippocampal subfield volume were estimated. A model-free, sliding-window approach was implemented that does not assume a linear relationship between age and subfield volume. The annualized percentage of subfield volume change was calculated to investigate associations with age, sex and apolipoprotein E e4 homozygosity. Hippocampal volume showed a marked reduction in apolipoprotein E e4/e4 female carriers after age 65. Volume was lower in homozygous e4 individuals in specific subfields including the presubiculum, subiculum head, cornu ammonis 1 body, cornu ammonis 3 head and cornu ammonis 4. Nearby brain structures in medial temporal and subcortical regions did not show the same age, sex and apolipoprotein E interactions, suggesting selective vulnerability of the hippocampus and its subfields. The findings demonstrate that in healthy ageing, two factors-female sex and apolipoprotein E e4 status-confer selective vulnerability of specific hippocampal subfields to volume loss.

Discovering correlates of age-related decline in a healthy late-midlife male birth cohort

Studies exploring age-related brain and cognitive change have identified substantial heterogeneity among individuals, but the underlying reasons for the differential trajectories remain largely unknown. We investigated cross-sectional and longitudinal associations between brain-imaging phenotypes (IDPs) and cognitive ability, and how these relations may be modified by common risk and protective factors. Participants were recruited from the 1953 Danish Male Birth Cohort (N=123), a longitudinal study of cognitive and brain ageing. Childhood IQ and socio-demographic factors are available for these participants who have been assessed regularly on multiple IDPs and behavioural factors in midlife. Using Pearson correlations and canonical correlation analysis (CCA), we explored the relation between 454 IDPs and 114 behavioural variables. CCA identified a single mode of population covariation coupling cross-subject longitudinal changes in brain structure to changes in cognitive performance and to a range of age-related covariates (r=0.92, Pcorrected < 0.001). Specifically, this CCA-mode indicated that; decreases in IQ and speed assessed tasks, higher rates of familial myocardial infarct, less physical activity, and poorer mental health are associated with larger decreases in whole brain grey matter and white matter. We found no evidence supporting the role of baseline scores as predictors of impending brain and behavioural change in late-midlife.

Common Genetic Variation Indicates Separate Causes for Periventricular and Deep White Matter Hyperintensities

BACKGROUND AND PURPOSE

Periventricular white matter hyperintensities (WMH; PVWMH) and deep WMH (DWMH) are regional classifications of WMH and reflect proposed differences in cause. In the first study, to date, we undertook genome-wide association analyses of DWMH and PVWMH to show that these phenotypes have different genetic underpinnings.

METHODS

Participants were aged 45 years and older, free of stroke and dementia. We conducted genome-wide association analyses of PVWMH and DWMH in 26,654 participants from CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology), ENIGMA (Enhancing Neuro-Imaging Genetics Through Meta-Analysis), and the UKB (UK Biobank). Regional correlations were investigated using the genome-wide association analyses -pairwise method. Cross-trait genetic correlations between PVWMH, DWMH, stroke, and dementia were estimated using LDSC.

RESULTS

In the discovery and replication analysis, for PVWMH only, we found associations on chromosomes 2 (NBEAL), 10q23.1 (TSPAN14/FAM231A), and 10q24.33 (SH3PXD2A). In the much larger combined meta-analysis of all cohorts, we identified ten significant regions for PVWMH: chromosomes 2 (3 regions), 6, 7, 10 (2 regions), 13, 16, and 17q23.1. New loci of interest include 7q36.1 (NOS3) and 16q24.2. In both the discovery/replication and combined analysis, we found genome-wide significant associations for the 17q25.1 locus for both DWMH and PVWMH. Using gene-based association analysis, 19 genes across all regions were identified for PVWMH only, including the new genes: CALCRL (2q32.1), KLHL24 (3q27.1), VCAN (5q27.1), and POLR2F (22q13.1). Thirteen genes in the 17q25.1 locus were significant for both phenotypes. More extensive genetic correlations were observed for PVWMH with small vessel ischemic stroke. There were no associations with dementia for either phenotype.

CONCLUSIONS

Our study confirms these phenotypes have distinct and also shared genetic architectures. Genetic analyses indicated PVWMH was more associated with ischemic stroke whilst DWMH loci were implicated in vascular, astrocyte, and neuronal function. Our study confirms these phenotypes are distinct neuroimaging classifications and identifies new candidate genes associated with PVWMH only.

Confound modelling in UK Biobank brain imaging

Dealing with confounds is an essential step in large cohort studies to address problems such as unexplained variance and spurious correlations. UK Biobank is a powerful resource for studying associations between imaging and non-imaging measures such as lifestyle factors and health outcomes, in part because of the large subject numbers. However, the resulting high statistical power also raises the sensitivity to confound effects, which therefore have to be carefully considered. In this work we describe a set of possible confounds (including nonlinear effects and interactions that researchers may wish to consider for their studies using such data). We include descriptions of how we can estimate the confounds, and study the extent to which each of these confounds affects the data, and the spurious correlations that may arise if they are not controlled. Finally, we discuss several issues that future studies should consider when dealing with confounds.

The UK Biobank imaging enhancement of 100,000 participants: rationale, data collection, management and future directions

UK Biobank is a population-based cohort of half a million participants aged 40-69 years recruited between 2006 and 2010. In 2014, UK Biobank started the world's largest multi-modal imaging study, with the aim of re-inviting 100,000 participants to undergo brain, cardiac and abdominal magnetic resonance imaging, dual-energy X-ray absorptiometry and carotid ultrasound. The combination of large-scale multi-modal imaging with extensive phenotypic and genetic data offers an unprecedented resource for scientists to conduct health-related research. This article provides an in-depth overview of the imaging enhancement, including the data collected, how it is managed and processed, and future directions.

Brain aging comprises many modes of structural and functional change with distinct genetic and biophysical associations

Brain imaging can be used to study how individuals' brains are aging, compared against population norms. This can inform on aspects of brain health; for example, smoking and blood pressure can be seen to accelerate brain aging. Typically, a single 'brain age' is estimated per subject, whereas here we identified 62 modes of subject variability, from 21,407 subjects' multimodal brain imaging data in UK Biobank. The modes represent different aspects of brain aging, showing distinct patterns of functional and structural brain change, and distinct patterns of association with genetics, lifestyle, cognition, physical measures and disease. While conventional brain-age modelling found no genetic associations, 34 modes had genetic associations. We suggest that it is important not to treat brain aging as a single homogeneous process, and that modelling of distinct patterns of structural and functional change will reveal more biologically meaningful markers of brain aging in health and disease.

Handedness, language areas and neuropsychiatric diseases: insights from brain imaging and genetics

Ninety per cent of the human population has been right-handed since the Paleolithic, yet the brain signature and genetic basis of handedness remain poorly characterized. Here, we correlated brain imaging phenotypes from ∼9000 UK Biobank participants with handedness, and with loci found significantly associated with handedness after we performed genome-wide association studies (GWAS) in ∼400 000 of these participants. Our imaging-handedness analysis revealed an increase in functional connectivity between left and right language networks in left-handers. GWAS of handedness uncovered four significant loci (rs199512, rs45608532, rs13017199, and rs3094128), three of which are in-or expression quantitative trait loci of-genes encoding proteins involved in brain development and patterning. These included microtubule-related MAP2 and MAPT, as well as WNT3 and MICB, all implicated in the pathogenesis of diseases such as Parkinson's, Alzheimer's and schizophrenia. In particular, with rs199512, we identified a common genetic influence on handedness, psychiatric phenotypes, Parkinson's disease, and the integrity of white matter tracts connecting the same language-related regions identified in the handedness-imaging analysis. This study has identified in the general population genome-wide significant loci for human handedness in, and expression quantitative trait loci of, genes associated with brain development, microtubules and patterning. We suggest that these genetic variants contribute to neurodevelopmental lateralization of brain organization, which in turn influences both the handedness phenotype and the predisposition to develop certain neurological and psychiatric diseases.

Discovering markers of healthy aging: a prospective study in a Danish male birth cohort

There is a pressing need to identify markers of cognitive and neural decline in healthy late-midlife participants. We explored the relationship between cross-sectional structural brain-imaging derived phenotypes (IDPs) and cognitive ability, demographic, health and lifestyle factors (non-IDPs). Participants were recruited from the 1953 Danish Male Birth Cohort (N=193). Applying an extreme group design, members were selected in 2 groups based on cognitive change between IQ at age ~20y (IQ-20) and age ~57y (IQ-57). Subjects showing the highest (n=95) and lowest (n=98) change were selected (at age ~57) for assessments on multiple IDPs and non-IDPs. We investigated the relationship between 453 IDPs and 70 non-IDPs through pairwise correlation and multivariate canonical correlation analysis (CCA) models. Significant pairwise associations included positive associations between IQ-20 and gray-matter volume of the temporal pole. CCA identified a richer pattern - a single "positive-negative" mode of population co-variation coupling individual cross-subject variations in IDPs to an extensive range of non-IDP measures (r = 0.75, P_corrected < 0.01). Specifically, this mode linked higher cognitive performance, positive early-life social factors, and mental health to a larger brain volume of several brain structures, overall volume, and microstructural properties of some white matter tracts. Interestingly, both statistical models identified IQ-20 and gray-matter volume of the temporal pole as important contributors to the inter-individual variation observed. The converging patterns provide novel insight into the importance of early adulthood intelligence as a significant marker of late-midlife neural decline and motivates additional study.

Hippocampal volume across age: Nomograms derived from over 19,700 people in UK Biobank

Measurement of hippocampal volume has proven useful to diagnose and track progression in several brain disorders, most notably in Alzheimer's disease (AD). For example, an objective evaluation of a patient's hippocampal volume status may provide important information that can assist diagnosis or risk stratification of AD. However, clinicians and researchers require access to age-related normative percentiles to reliably categorise a patient's hippocampal volume as being pathologically small. Here we analysed effects of age, sex, and hemisphere on the hippocampus and neighbouring temporal lobe volumes, in 19,793 generally healthy participants in the UK Biobank. A key finding of the current study is a significant acceleration in the rate of hippocampal volume loss in middle age, more pronounced in females than in males. In this report, we provide normative values for hippocampal and total grey matter volume as a function of age for reference in clinical and research settings. These normative values may be used in combination with our online, automated percentile estimation tool to provide a rapid, objective evaluation of an individual's hippocampal volume status. The data provide a large-scale normative database to facilitate easy age-adjusted determination of where an individual hippocampal and temporal lobe volume lies within the normal distribution.

•

Largest normative database for hippocampal volume across age to date.

•

Hippocampal volume loss accelerates in middle age.

•

Acceleration of hippocampal volume loss is more pronounced in women than in men.

•

With the online tool the database provides a useful resource for research and clinical studies.

Estimation of brain age delta from brain imaging

It is of increasing interest to study “brain age” - the apparent age of a subject, as inferred from brain imaging data. The difference between brain age and actual age (the “delta”) is typically computed, reflecting deviation from the population norm. This therefore may reflect accelerated aging (positive delta) or resilience (negative delta) and has been found to be a useful correlate with factors such as disease and cognitive decline. However, although there has been a range of methods proposed for estimating brain age, there has been little study of the optimal ways of computing the delta. In this technical note we describe problems with the most common current approach, and present potential improvements. We evaluate different estimation methods on simulated and real data. We also find the strongest correlations of corrected brain age delta with 5,792 non-imaging variables (non-brain physical measures, life-factor measures, cognitive test scores, etc.), and also with 2,641 multimodal brain imaging-derived phenotypes, with data from 19,000 participants in UK Biobank.

•

It is of interest to study "brain age'', as inferred from brain imaging data.

•

The delta between brain age and actual age is typically computed.

•

We describe problems with the most common current approach.

•

We present potential improvements.

•

We evaluate methods on simulated data and data from UK Biobank.

Modelling the distribution of white matter hyperintensities due to ageing on MRI images using Bayesian inference

White matter hyperintensities (WMH), also known as white matter lesions, are localised white matter areas that appear hyperintense on MRI scans. WMH commonly occur in the ageing population, and are often associated with several factors such as cognitive disorders, cardiovascular risk factors, cerebrovascular and neurodegenerative diseases. Despite the fact that some links between lesion location and parametric factors such as age have already been established, the relationship between voxel-wise spatial distribution of lesions and these factors is not yet well understood. Hence, it would be of clinical importance to model the distribution of lesions at the population-level and quantitatively analyse the effect of various factors on the lesion distribution model. In this work we compare various methods, including our proposed method, to generate voxel-wise distributions of WMH within a population with respect to various factors. Our proposed Bayesian spline method models the spatio-temporal distribution of WMH with respect to a parametric factor of interest, in this case age, within a population. Our probabilistic model takes as input the lesion segmentation binary maps of subjects belonging to various age groups and provides a population-level parametric lesion probability map as output. We used a spline representation to ensure a degree of smoothness in space and the dimension associated with the parameter, and formulated our model using a Bayesian framework. We tested our algorithm output on simulated data and compared our results with those obtained using various existing methods with different levels of algorithmic and computational complexity. We then compared the better performing methods on a real dataset, consisting of 1000 subjects of the UK Biobank, divided in two groups based on hypertension diagnosis. Finally, we applied our method on a clinical dataset of patients with vascular disease. On simulated dataset, the results from our algorithm showed a mean square error (MSE) value of 7.27×10-5, which was lower than the MSE value reported in the literature, with the advantage of being robust and computationally efficient. In the UK Biobank data, we found that the lesion probabilities are higher for the hypertension group compared to the non-hypertension group and further verified this finding using a statistical t-test. Finally, when applying our method on patients with vascular disease, we observed that the overall probability of lesions is significantly higher in later age groups, which is in line with the current literature.

Discovering dynamic brain networks from big data in rest and task

Brain activity is a dynamic combination of the responses to sensory inputs and its own spontaneous processing. Consequently, such brain activity is continuously changing whether or not one is focusing on an externally imposed task. Previously, we have introduced an analysis method that allows us, using Hidden Markov Models (HMM), to model task or rest brain activity as a dynamic sequence of distinct brain networks, overcoming many of the limitations posed by sliding window approaches. Here, we present an advance that enables the HMM to handle very large amounts of data, making possible the inference of very reproducible and interpretable dynamic brain networks in a range of different datasets, including task, rest, MEG and fMRI, with potentially thousands of subjects. We anticipate that the generation of large and publicly available datasets from initiatives such as the Human Connectome Project and UK Biobank, in combination with computational methods that can work at this scale, will bring a breakthrough in our understanding of brain function in both health and disease.

Genome-wide association studies of brain imaging phenotypes in UK Biobank

The genetic architecture of brain structure and function is largely unknown. To investigate this, we carried out genome-wide association studies of 3,144 functional and structural brain imaging phenotypes from UK Biobank (discovery dataset 8,428 subjects). Here we show that many of these phenotypes are heritable. We identify 148 clusters of associations between single nucleotide polymorphisms and imaging phenotypes that replicate at P < 0.05, when we would expect 21 to replicate by chance. Notable significant, interpretable associations include: iron transport and storage genes, related to magnetic susceptibility of subcortical brain tissue; extracellular matrix and epidermal growth factor genes, associated with white matter micro-structure and lesions; genes that regulate mid-line axon development, associated with organization of the pontine crossing tract; and overall 17 genes involved in development, pathway signalling and plasticity. Our results provide insights into the genetic architecture of the brain that are relevant to neurological and psychiatric disorders, brain development and ageing.

Automated quality control for within and between studies diffusion MRI data using a non-parametric framework for movement and distortion correction

Diffusion MRI data can be affected by hardware and subject-related artefacts that can adversely affect downstream analyses. Therefore, automated quality control (QC) is of great importance, especially in large population studies where visual QC is not practical. In this work, we introduce an automated diffusion MRI QC framework for single subject and group studies. The QC is based on a comprehensive, non-parametric approach for movement and distortion correction: FSL EDDY, which allows us to extract a rich set of QC metrics that are both sensitive and specific to different types of artefacts. Two different tools are presented: QUAD (QUality Assessment for DMRI), for single subject QC and SQUAD (Study-wise QUality Assessment for DMRI), which is designed to enable group QC and facilitate cross-studies harmonisation efforts.

•

Two tools to automatically perform QC of diffusion MRI data.

•

Automated generation of single subject reports for visual inspection and database.

•

Group databases and reports allow to compare subjects within and between studies.

•

Categorical and continuous variables can be used to update the reports.

An empirical, 21st century evaluation of phrenology

Phrenology was a nineteenth century endeavour to link personality traits with scalp morphology, which has been both influential and fiercely criticised, not least because of the assumption that scalp morphology can be informative of underlying brain function. Here we test the idea empirically rather than dismissing it out of hand. Whereas nineteenth century phrenologists had access to coarse measurement tools (digital technology referring then to fingers), we were able to re-examine phrenology using 21st century methods and thousands of subjects drawn from the largest neuroimaging study to date. High-quality structural MRI was used to quantify local scalp curvature. The resulting curvature statistics were compared against lifestyle measures acquired from the same cohort of subjects, being careful to match a subset of lifestyle measures to phrenological ideas of brain organisation, in an effort to evoke the character of Victorian times. The results represent the most rigorous evaluation of phrenological claims to date.

Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank

UK Biobank is a large-scale prospective epidemiological study with all data accessible to researchers worldwide. It is currently in the process of bringing back 100,000 of the original participants for brain, heart and body MRI, carotid ultrasound and low-dose bone/fat x-ray. The brain imaging component covers 6 modalities (T1, T2 FLAIR, susceptibility weighted MRI, Resting fMRI, Task fMRI and Diffusion MRI). Raw and processed data from the first 10,000 imaged subjects has recently been released for general research access. To help convert this data into useful summary information we have developed an automated processing and QC (Quality Control) pipeline that is available for use by other researchers. In this paper we describe the pipeline in detail, following a brief overview of UK Biobank brain imaging and the acquisition protocol. We also describe several quantitative investigations carried out as part of the development of both the imaging protocol and the processing pipeline.

Investigations into within- and between-subject resting-state amplitude variations

The amplitudes of spontaneous fluctuations in brain activity may be a significant source of within-subject and between-subject variability, and this variability is likely to be carried through into functional connectivity (FC) estimates (whether directly or indirectly). Therefore, improving our understanding of amplitude fluctuations over the course of a resting state scan and variation in amplitude across individuals is of great relevance to the interpretation of FC findings. We investigate resting state amplitudes in two large-scale studies (HCP and UK Biobank), with the aim of determining between-subject and within-subject variability. Between-subject clustering distinguished between two groups of brain networks whose amplitude variation across subjects were highly correlated with each other, revealing a clear distinction between primary sensory and motor regions ('primary sensory/motor cluster') and cognitive networks. Within subjects, all networks in the primary sensory/motor cluster showed a consistent increase in amplitudes from the start to the end of the scan. In addition to the strong increases in primary sensory/motor amplitude, a large number of changes in FC were found when comparing the two scans acquired on the same day (HCP data). Additive signal change analysis confirmed that all of the observed FC changes could be fully explained by changes in amplitude. Between-subject correlations in UK Biobank data showed a negative correlation between primary sensory/motor amplitude and average sleep duration, suggesting a role of arousal. Our findings additionally reveal complex relationships between amplitude and head motion. These results suggest that network amplitude is a source of significant variability both across subjects, and within subjects on a within-session timescale. Future rfMRI studies may benefit from obtaining arousal-related (self report) measures, and may wish to consider the influence of amplitude changes on measures of (dynamic) functional connectivity.

BIDS apps: Improving ease of use, accessibility, and reproducibility of neuroimaging data analysis methods

The rate of progress in human neurosciences is limited by the inability to easily apply a wide range of analysis methods to the plethora of different datasets acquired in labs around the world. In this work, we introduce a framework for creating, testing, versioning and archiving portable applications for analyzing neuroimaging data organized and described in compliance with the Brain Imaging Data Structure (BIDS). The portability of these applications (BIDS Apps) is achieved by using container technologies that encapsulate all binary and other dependencies in one convenient package. BIDS Apps run on all three major operating systems with no need for complex setup and configuration and thanks to the comprehensiveness of the BIDS standard they require little manual user input. Previous containerized data processing solutions were limited to single user environments and not compatible with most multi-tenant High Performance Computing systems. BIDS Apps overcome this limitation by taking advantage of the Singularity container technology. As a proof of concept, this work is accompanied by 22 ready to use BIDS Apps, packaging a diverse set of commonly used neuroimaging algorithms.

Multimodal population brain imaging in the UK Biobank prospective epidemiological study

Medical imaging has enormous potential for early disease prediction, but is impeded by the difficulty and expense of acquiring data sets before symptom onset. UK Biobank aims to address this problem directly by acquiring high-quality, consistently acquired imaging data from 100,000 predominantly healthy participants, with health outcomes being tracked over the coming decades. The brain imaging includes structural, diffusion and functional modalities. Along with body and cardiac imaging, genetics, lifestyle measures, biological phenotyping and health records, this imaging is expected to enable discovery of imaging markers of a broad range of diseases at their earliest stages, as well as provide unique insight into disease mechanisms. We describe UK Biobank brain imaging and present results derived from the first 5,000 participants' data release. Although this covers just 5% of the ultimate cohort, it has already yielded a rich range of associations between brain imaging and other measures collected by UK Biobank.

Cortical morphology of adolescents with bipolar disorder and with schizophrenia

INTRODUCTION

Recent evidence points to overlapping decreases in cortical thickness and gyrification in the frontal lobe of patients with adult-onset schizophrenia and bipolar disorder with psychotic symptoms, but it is not clear if these findings generalize to patients with a disease onset during adolescence and what may be the mechanisms underlying a decrease in gyrification.

METHOD

This study analyzed cortical morphology using surface-based morphometry in 92 subjects (age range 11-18 years, 52 healthy controls and 40 adolescents with early-onset first-episode psychosis diagnosed with schizophrenia (n=20) or bipolar disorder with psychotic symptoms (n=20) based on a two year clinical follow up). Average lobar cortical thickness, surface area, gyrification index (GI) and sulcal width were compared between groups, and the relationship between the GI and sulcal width was assessed in the patient group.

RESULTS

Both patients groups showed decreased cortical thickness and increased sulcal width in the frontal cortex when compared to healthy controls. The schizophrenia subgroup also had increased sulcal width in all other lobes. In the frontal cortex of the combined patient group sulcal width was negatively correlated (r=-0.58, p<0.001) with the GI.

CONCLUSIONS

In adolescents with schizophrenia and bipolar disorder with psychotic symptoms there is cortical thinning, decreased GI and increased sulcal width of the frontal cortex present at the time of the first psychotic episode. Decreased frontal GI is associated with the widening of the frontal sulci which may reduce sulcal surface area. These results suggest that abnormal growth (or more pronounced shrinkage during adolescence) of the frontal cortex represents a shared endophenotype for psychosis.