Paul M. Matthews

Advances in functional and structural MR image analysis and implementation as FSL.

The techniques available for the interrogation and analysis of neuroimaging data have a large influence in determining the flexibility, sensitivity, and scope of neuroimaging experiments. The development of such methodologies has allowed investigators to address scientific questions that could not previously be answered and, as such, has become an important research area in its own right. In this paper, we present a review of the research carried out by the Analysis Group at the Oxford Centre for Functional MRI of the Brain (FMRIB). This research has focussed on the development of new methodologies for the analysis of both structural and functional magnetic resonance imaging data. The majority of the research laid out in this paper has been implemented as freely available software tools within FMRIBs Software Library (FSL).

Characterization and propagation of uncertainty in diffusion-weighted MR imaging.

A fully probabilistic framework is presented for estimating local probability density functions on parameters of interest in a model of diffusion. This technique is applied to the estimation of parameters in the diffusion tensor model, and also to a simple partial volume model of diffusion. In both cases the parameters of interest include parameters defining local fiber direction. A technique is then presented for using these density functions to estimate global connectivity (i.e., the probability of the existence of a connection through the data field, between any two distant points), allowing for the quantification of belief in tractography results. This technique is then applied to the estimation of the cortical connectivity of the human thalamus. The resulting connectivity distributions correspond well with predictions from invasive tracer methods in nonhuman primate. Magn Reson Med 50:1077–1088, 2003.

Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging

Evidence concerning anatomical connectivities in the human brain is sparse and based largely on limited post-mortem observations. Diffusion tensor imaging has previously been used to define large white-matter tracts in the living human brain, but this technique has had limited success in tracing pathways into gray matter. Here we identified specific connections between human thalamus and cortex using a novel probabilistic tractography algorithm with diffusion imaging data. Classification of thalamic gray matter based on cortical connectivity patterns revealed distinct subregions whose locations correspond to nuclei described previously in histological studies. The connections that we found between thalamus and cortex were similar to those reported for non-human primates and were reproducible between individuals. Our results provide the first quantitative demonstration of reliable inference of anatomical connectivity between human gray matter structures using diffusion data and the first connectivity-based segmentation of gray matter.

Large recurrent microdeletions associated with schizophrenia.

Reduced fecundity, associated with severe mental disorders, places negative selection pressure on risk alleles and may explain, in part, why common variants have not been found that confer risk of disorders such as autism, schizophrenia and mental retardation. Thus, rare variants may account for a larger fraction of the overall genetic risk than previously assumed. In contrast to rare single nucleotide mutations, rare copy number variations (CNVs) can be detected using genome-wide single nucleotide polymorphism arrays. This has led to the identification of CNVs associated with mental retardation and autism. In a genome-wide search for CNVs associating with schizophrenia, we used a population-based sample to identify de novo CNVs by analysing 9,878 transmissions from parents to offspring. The 66 de novo CNVs identified were tested for association in a sample of 1,433 schizophrenia cases and 33,250 controls. Three deletions at 1q21.1, 15q11.2 and 15q13.3 showing nominal association with schizophrenia in the first sample (phase I) were followed up in a second sample of 3,285 cases and 7,951 controls (phase II). All three deletions significantly associate with schizophrenia and related psychoses in the combined sample. The identification of these rare, recurrent risk variants, having occurred independently in multiple founders and being subject to negative selection, is important in itself. CNV analysis may also point the way to the identification of additional and more prevalent risk variants in genes and pathways involved in schizophrenia.

Accurate, robust, and automated longitudinal and cross-sectional brain change analysis.

Quantitative measurement of brain size, shape, and temporal change (for example, in order to estimate atrophy) is increasingly important in biomedical image analysis applications. New methods of structural analysis attempt to improve robustness, accuracy, and extent of automation. A fully automated method of longitudinal (temporal change) analysis, SIENA, was presented previously. In this paper, improvements to this method are described, and also an extension of SIENA to a new method for cross-sectional (single time point) analysis. The methods are fully automated, robust, and accurate: 0.15% brain volume change error (longitudinal): 0.5-1% brain volume accuracy for single-time point (cross-sectional). A particular advantage is the relative insensitivity to differences in scanning parameters. The methods provide easy manual review of their output by the automatic production of summary images which show the results of the brain extraction, registration, tissue segmentation, and final atrophy estimation.

Common variants conferring risk of schizophrenia

Schizophrenia is a complex disorder, caused by both genetic and environmental factors and their interactions. Research on pathogenesis has traditionally focused on neurotransmitter systems in the brain, particularly those involving dopamine. Schizophrenia has been considered a separate disease for over a century, but in the absence of clear biological markers, diagnosis has historically been based on signs and symptoms. A fundamental message emerging from genome-wide association studies of copy number variations (CNVs) associated with the disease is that its genetic basis does not necessarily conform to classical nosological disease boundaries. Certain CNVs confer not only high relative risk of schizophrenia but also of other psychiatric disorders. The structural variations associated with schizophrenia can involve several genes and the phenotypic syndromes, or the ‘genomic disorders’, have not yet been characterized. Single nucleotide polymorphism (SNP)-based genome-wide association studies with the potential to implicate individual genes in complex diseases may reveal underlying biological pathways. Here we combined SNP data from several large genome-wide scans and followed up the most significant association signals. We found significant association with several markers spanning the major histocompatibility complex (MHC) region on chromosome 6p21.3-22.1, a marker located upstream of the neurogranin gene (NRGN) on 11q24.2 and a marker in intron four of transcription factor 4 (TCF4) on 18q21.2. Our findings implicating the MHC region are consistent with an immune component to schizophrenia risk, whereas the association with NRGN and TCF4 points to perturbation of pathways involved in brain development, memory and cognition.

Distinct patterns of brain activity in young carriers of the APOE-epsilon4 allele.

The APOE ε4 allele is a risk factor for late-life pathological changes that is also associated with anatomical and functional brain changes in middle-aged and elderly healthy subjects. We investigated structural and functional effects of the APOE polymorphism in 18 young healthy APOE ε4-carriers and 18 matched noncarriers (age range: 20–35 years). Brain activity was studied both at rest and during an encoding memory paradigm using blood oxygen level-dependent fMRI. Resting fMRI revealed increased “default mode network” (involving retrosplenial, medial temporal, and medial-prefrontal cortical areas) coactivation in ε4-carriers relative to noncarriers. The encoding task produced greater hippocampal activation in ε4-carriers relative to noncarriers. Neither result could be explained by differences in memory performance, brain morphology, or resting cerebral blood flow. The APOE ε4 allele modulates brain function decades before any clinical or neurophysiological expression of neurodegenerative processes.

Oxygenation dependence of the transverse relaxation time of water protons in whole blood at high field.

At high and medium magnetic field, the transverse NMR relaxation rate (T-1(2)) of water proteins in blood is determined predominantly by the oxygenation state of haemoglobin. T-1(2) depends quadratically on the field strength and on the proportion of haemoglobin that is deoxygenated. Deoxygenation increases the volume magnetic susceptibility within the erythrocytes and thus creates local field gradients around these cells. From volume susceptibility measurements and the dependence of T-1(2) on the pulse rate in the Carr-Purcell-Meiboom-Gill experiment, we show that the increase in T-1(2) with increasing blood deoxygenation arises from diffusion of water through these field gradients.

Meta-analysis of genome scans and replication identify CD6, IRF8 and TNFRSF1A as new multiple sclerosis susceptibility loci

We report the results of a meta-analysis of genome-wide association scans for multiple sclerosis (MS) susceptibility that includes 2,624 subjects with MS and 7,220 control subjects. Replication in an independent set of 2,215 subjects with MS and 2,116 control subjects validates new MS susceptibility loci at TNFRSF1A (combined P = 1.59 × 10−11), IRF8 (P = 3.73 × 10−9) and CD6 (P = 3.79 × 10−9). TNFRSF1A harbors two independent susceptibility alleles: rs1800693 is a common variant with modest effect (odds ratio = 1.2), whereas rs4149584 is a nonsynonymous coding polymorphism of low frequency but with stronger effect (allele frequency = 0.02; odds ratio = 1.6). We also report that the susceptibility allele near IRF8, which encodes a transcription factor known to function in type I interferon signaling, is associated with higher mRNA expression of interferon-response pathway genes in subjects with MS.

The role of ipsilateral premotor cortex in hand movement after stroke

Movement of an affected hand after stroke is associated with increased activation of ipsilateral motor cortical areas, suggesting that these motor areas in the undamaged hemisphere may adaptively compensate for damaged or disconnected regions. However, this adaptive compensation has not yet been demonstrated directly. Here we used transcranial magnetic stimulation (TMS) to interfere transiently with processing in the ipsilateral primary motor or dorsal premotor cortex (PMd) during finger movements. TMS had a greater effect on patients than controls in a manner that depended on the site, hemisphere, and time of stimulation. In patients with right hemiparesis (but not in healthy controls), TMS applied to PMd early (100 ms) after the cue to move slowed simple reaction-time finger movements by 12% compared with controls. The relative slowing of movements with ipsilateral PMd stimulation in patients correlated with the degree of motor impairment, suggesting that functional recruitment of ipsilateral motor areas was greatest in the more impaired patients. We also used functional magnetic resonance imaging to monitor brain activity in these subjects as they performed the same movements. Slowing of reaction time after premotor cortex TMS in the patients correlated inversely with the relative hemispheric lateralization of functional magnetic resonance imaging activation in PMd. This inverse correlation suggests that the increased activation in ipsilateral cortical motor areas during movements of a paretic hand, shown in this and previous functional imaging studies, represents a functionally relevant, adaptive response to the associated brain injury.