Additional related information may be found at: |
Neuropsychopharmacology: The Fifth Generation of Progress |
In Vivo Structural Brain Assessment
Kelvin O. Lim, Margaret Rosenbloom, and Adolf Pfefferbaum
Over the past few years, technical improvements in the acquisition and processing of structural neuroimaging data, particularly with magnetic resonance imaging (MRI), have provided vivid visual representations of the human brain, both as two-dimensional (2D) slices through the brain, and as three-dimensional (3D) views of the external surface and internal structures of the brain. Computerized tomography (CT) was the first imaging modality to provide in vivo evidence of gross brain morphological abnormalities in schizophrenia (39), with many CT reports of increase in cerebrospinal fluid (CSF)-filled spaces, both centrally (ventricles), and peripherally (sulci) in a variety of psychiatric patients (80). MRI, with its superior tissue contrast for differentiation of white matter from gray matter, and greater flexibility in acquiring, or re-slicing images in orientations best suited for viewing and measuring specific structures, has provided additional capabilities for uncovering both generalized and regionally specific morphometric brain abnormalities associated with schizophrenia, and other mental disorders (46, 88). Patients with psychiatric disorders generally manifest brain abnormalities of a subtle and/or diffuse nature, best characterized by systematic and reliable quantification techniques rather than qualitative assessments.
In addition to characterizing ever more specific brain dysmorphology (the size, tissue composition, and shape of different cortical or subcortical neuroanatomic structures), new advances have also made it possible to use MRI to visualize and quantify the directional coherence of white matter fibers (diffusion tensor imaging, DTI) (9), an application of potential relevance for investigating levels of connectivity and dysconnectivity between different brain regions.
All measures derived from neuroimaging technologies need to be carefully interpreted in the context of other brain and biological features, as well as the clinical and demographic characteristics of the patient, before they can provide insight into the possible etiology, pathophysiology and course of different mental disorders. The goal of this updated chapter is to briefly recapitulate prior information, describe some new advances in neuroimaging and outline their capabilities and limitations and reiterate some of the important technical and methodological issues connected with acquiring, measuring and interpreting these types of data.
ACQUISITION OF MAGNETIC RESONANCE IMAGES
Mechanisms Of Tissue Contrast: Pulse Sequences
One of the strengths of MRI is that investigators can manipulate the amount of contrast between different biological tissues by varying elements of the image acquisition sequence. A highly simplified outline of the principle of MRI is provided here to orient the reader to how these acquisition elements can be varied. More detailed descriptions can be found elsewhere (68). Atomic nuclei with an odd number of nucleons (protons plus neutrons) behave as spinning tops with an electric charge, and thus they produce a magnetic field, oriented along the axis of rotation of the nucleus. Under normal circumstances, the spin axes of these nuclei are randomly oriented so that the sum of all magnetic fields is zero. However, when these nuclei are exposed to a strong external magnetic field, they align either parallel or anti-parallel with the field in an equilibrium state, that is either in a high or low energy state, with a slight preponderance of the high energy state. These aligned atomic nuclei rotate around their axes at a characteristic frequency, also dependent on the strength of the magnetic field. Electromagnetic radiofrequency (RF) pulses, applied at the appropriate frequency perpendicular to the main magnetic field, excite the spinning nuclei, knocking them out of equilibrium. When the radio waves are turned off, the nuclei return to equilibrium. MRI employs appropriately tuned radiofrequency coils to detect the radio signals given off as the nuclei return to equilibrium, and localize them using magnetic gradients in the three orthogonal axes.
Hydrogen atoms (protons) are the most abundant element in biological tissue, primarily in the form of water. Proton MRI thus records signals arising predominantly from free water (i.e., water not biochemically bound in complex molecules), but there is also a measurable but much lesser contribution from protons in fat. Proton density (PD), based on the number of nuclei stimulated, and the relaxation times, T1 and T2, which reflect the chemical environment of the tissue stimulated, each determine in different ways the intensity at which different tissues will appear in an acquired image. T1 is an exponential time constant which represents the time taken for excited nuclei to return to equilibrium after the RF pulse has been turned off. T2 is an exponential time constant describing the time it takes for the excited nuclei to lose signal, mainly due to dephasing in the transverse plane. The time between radiofrequency pulses (TR), and the amount of time after the pulse at which the signal is acquired (echo time, or TE), determine the influence each of these parameters (PD, T1 and T2) have on the acquired image. For example a long TR (2-3 sec), spin echo sequence with a long TE (80 msec) provides a heavily "T2-weighted" image in which CSF remains bright, because CSF has a long T2, while signal from brain parenchyma,which has a shorter T2, has declined, rendering this tisue darker and homogeneous. In contrast, a relatively short TR (200-400 msec) spin echo or inversion recovery sequence with a short TE (20 msec) provides a "T1-weighted" image which optimizes contrast between gray matter and white matter. In T1-weighted images, tissue with shorter T1s are brighter, whereas in T2 weighted images, tissues with longer T2s are brighter. For brain imaging, the choice of T1 vs. T2 weighted images presents a tradeoff between obtaining high contrast between gray matter and white matter or obtaining high contrast between CSF and brain parenchyma. While T1-weighted spin-echo and inversion recovery sequences provide excellent white-gray contrast, and have been used extensively for morphometric studies, they are limited by their poor definition of CSF/skull margins for reliably measuring intracranial volume. Thus, there is no single correct MR brain sequence, but rather a range of options, with the optimal selection determined by the questions being asked.
A major limitation in scanning patients has been the length of time needed to obtain images, resulting in image degradation due to subject motion. A major thrust for development of clinical imaging techniques has been to reduce imaging time, either by reducing TR and TE while incorporating additional pulses to provide contrast, or by acquiring more data after each excitation pulse. Techniques that acquire the same amount of data in less time produce image contrasts similar, but not identical, to those of traditional spin-echo sequences. Sequences that acquire all of the data needed for reconstructing an image after a single excitation pulse (single shot techniques) include Echo Planar Imaging (EPI) (53) and spiral imaging (61) techniques. While fast and convenient, these images often suffer from lower signal-to-noise ratio (SNR) and susceptibility induced distortions. Work is constantly under way by the manufacturers and research groups to produce further enhancements and overcome limitations. Approaches that minimize the influence of susceptibility artifacts of EPI based systems by acquiring only part of the data after each excitation include Rapid Acquisition with Relaxation Enhancement (RARE) and Fast Spin Echo (FSE) sequences.
Sources of contrast other than that based on manipulation of T1, T2, and proton density have been exploited to extend the information available through MRI. One of these is based on magnetization transfer, which exploits the difference between free and restricted (typically in myelin) hydrogen molecules. Off-resonance RF pulses are used to selectively saturate the restricted pool, inducing an exchange of magnetization between free and restricted hydrogen protons and promoting T1 relaxation. The resulting contrast is believed to reflect the structural integrity of the tissue being imaged (27, 58) particularly white matter.
Another source of contrast is based on the self-diffusion or the random Brownian movement of water molecules. When diffusion is unrestricted, e.g. as in a glass of water, the movement of water molecules has no directional preference and is termed isotropic. However, when the motion is restricted, e.g., by fibers such as constitute the white matter tracts, diffusion is greater parallel to the fiber than perpendicular to it, and is termed anisotropic. Anisotropic diffusion can be imaged, and by sequentially applying magnetic gradients in six non-collinear directions, its directionality and magnitude can be quantified. This form of imaging has been called diffusion tensor imaging (DTI) (9). Subsequent analysis applies a mathematical construct, a tensor, to describe the multidimensional vector systems and quantify the diffusion. Vector measures have been used to map white matter fiber tracts in the brain. Scalar measures of the diffusion, such as fractional anisotropy and relative anisotropy, provide an overall index of tissue anisotropy (52). This approach has recently been used to reveal possibly compromised white matter tracts in patients with schizophrenia (13, 50).
Two-dimensional Multi-slice and Three-dimensional Imaging
Two-dimensional (2D)
images are acquired in a defined orthogonal plane, typically in axial, sagittal
and coronal orientations (Figure
1). Image orientation is defined
by use of selective RF pulse excitations and the appropriate magnetic gradients
in the three orthogonal axes. Each orientation
provides a different view of the brain with optimal visualization of different
structures. For example, mid-sagittal
images provide clear delineation of the prefrontal cortex and the corpus callosum,
coronal images provide views of the hippocampus and limbic structures over several
sections or slices and axial images provide a good view of basal ganglia structures
such as the putamen, globus pallidus, caudate and substantia nigra, as well
as the lateral ventricular system. One limitation of 2D image acquisition is that
only selected slices, of specified thickness, are excited and imaged and, unless
these slices have been collected using standardized internal neuroanatomic landmarks,
they may not be directly comparable across subjects. Comparability is enhanced by orienting each
slice acquisition relative to a specific anatomic plane such as that passing
through the anterior and posterior commissures, perpendicular to the sagittal
plane (95).
Three-dimensional (3D) volume acquisition protocols encompass the entire brain and can be reformatted, post acquisition, into any plane subject to the limitations imposed by the resolution of the original data. This post-acquisition flexibility in defining 3D planes for analysis can reduce measurement errors related to head misalignment between subjects, or within a subject in a longitudinal study. 3D images are acquired by exciting a broad volume of tissue with a non-slice selective RF pulse (5). Additional magnetic gradients are then used to spatially encode within the excited region. After the raw image data have been collected, a Fourier transform decodes slice positions from the entire data set. These slice data are further reconstructed with a 2D Fourier transform into the final image. Limitations in post-scan analysis of 3D data sets are discussed in more detail in the section on image analysis below. Among the 3D acquisition protocols which are widely used in psychiatric neuroimaging are spoiled gradient recall acquisition (SPGR) and magnetization prepared rapid gradient echo (MPRAGE) (63) each of which are T1-weighted, provide reasonable gray/white differentiation, and can sample the entire brain with 1.5mm or thinner slices in 10 minutes or less.
Image Resolution, Signal-To-Noise-Ratio, and Slice Thickness
An image is composed of a grid (usually 256 x 256) of 2D picture elements (pixels). The resolution of the image is affected both by in-plane resolution and by slice thickness. An in-plane resolution of 1 mm means that each pixel in the image matrix represents 1 mm2. However if the slice is as thick as 5 mm, once standard for 2D MRI protocols, each image pixel actually represent a voxel volume of 5 mm3. The larger the voxel, the greater the likelihood it will be "partially volumed", i.e., contain some combination of CSF, white matter or gray matter rather than being homogenous for brain tissue type. To reduce partial voluming, high in-plane resolution and thin slices are desired. However, technical factors such as gradient strength and the amount of RF energy used to excite the object, limit the minimum thickness effectively available to approximately 1-2 mm. In addition the SNR per voxel will be lower because there is less material producing the signal in a thin slice. Thus, while the resolution may be "finer" there will be more noise and less signal per pixel. Thinner slices may also require an increased number of acquisitions and thus longer scanning time to encompass an equivalent volume of the brain, with the attendant disadvantages of requiring subjects to remain still in the scanner.
Artifact
MR images are susceptible to various types of artifact, some of which can be minimized during image acquisition and others which can only be remedied during data analysis.
Hardware Induced Artifact
The RF coil's sensitivity is not homogeneous across all dimensions of the imaging plane which often introduces a low frequency gradient of signal intensity across a given image. While the human visual system can maintain contrast detection in the face of this artifact, most automated thresholding techniques assume a constant baseline level and require the RF inhomogeneity to be corrected before the image can be processed. A variety of solutions have been proposed to this problem, including filtering (51) bi-feature thresholding (44), and line by line inversion of the slope of the pixel values within each slice (17) and adaptive filtering (107).
Effects Of CSF And Blood Pulsation, And Subject Movement
Physiological sources of artifact include pulsation of blood and CSF through regions of the brain being imaged. The limbic system is particularly vulnerable to this artifact because of its vascularization pattern and relationship to the ventricular system. Cardiac-gating and flow compensated pulse sequences help reduce these sources of movement artifact, although not without cost. Because more gradients are used with flow compensation, the minimum echo time is lengthened and the maximum number of slices that can be acquired is reduced. Voluntary or involuntary gross head movements are another source of artifact to which 3D, functional MRI and DTI are particularly susceptible. Head restraints or bitebars (41, 60) provide one solution A complementary approach, adopted by some research groups for structural imaging, involves administering a sedative hypnotic, adequate to calm the subject or induce sleep.
Need For Phantoms
The accuracy and reproducibility of image data is directly related to the stability of the imaging hardware, especially the magnetic field gradient systems, the main magnetic field, and the RF pulse system. Routine imaging of a quality assurance phantom permits monitoring of the scanner's stability over time, stability between hardware upgrades, and establishment of standardized techniques for multi-center studies (101).
Tradeoffs
A recurring theme in any discussion of MRI data acquisition for clinical
studies is the constant tradeoff an investigator must make between resolution
of the image, its SNR, and time of the subject in the scanner.
Even with cooperative and relaxed subjects, involuntary head movements
during a lengthy scan session can be a problem. Optimal image acquisition parameters for specific studies need to
be determined empirically and will vary depending on the structure being imaged,
scanner hardware, and the clinical characteristics of the subjects being studied.
The high cost of MRI scanning time, along with limitations in patient
tolerance, combine to favor short acquisition sessions.
However, reductions in scanning time may be at the cost of poorer image
resolution and contrast, or sampling less of the brain.
A more detailed technical treatment of pulse design issues and their
relative strengths and weaknesses is available elsewhere (e.g., 85).
ANALYSIS OF MAGNETIC RESONANCE BRAIN IMAGES
Automated computerized analysis of many elements or characteristics of an image is becoming increasingly feasible and accessible to investigators. These approaches offer considerable advantages of flexibility and efficiency, though even automated analysis approaches may not overcome all limitations of image resolution, and distortions due to artifact in data collection. Most analysis approaches combine both automated and interactive elements (e.g., 24). Clinical reading or qualitative scoring of MRI brain data, still widely used by neuroradiologists, serves as a useful complement to quantification, particularly for assessment of image features that until recently have defied reliable quantification such as white matter lucencies or gross "abnormalities" in brain morphology. However recent advances in characterizing tissue quality (e.g., 96) and shape (e.g., 10, 19) of defined structures are making even these abnormalities amenable to quantification.
Resolution Issues
Resolution is limited by parameters of the data collection protocol, including
in-plane resolution and slice thickness. Generally
2D data are analyzed as acquired, with a known in-plane resolution and slice
thickness. In contrast, 3D data sets
offer the opportunity for reformatting images in user-defined planes in any
orientation and slice thickness. This provides greater flexibility in optimizing
the view of specific structures. However,
re-slicing is bound, not only by the hardware and software constraints of the
acquisition parameters, but also by constraints of the object space ( i.e.,
the defined plane and thickness) of the original images.
The effective thickness of a user-defined slice will depend on the orientation
of the slice in the 3D volume. For example,
with an isotropic data set, (i.e., each voxel is a true cube with side length
= L) where resolution in the x, y, and z axes is equal, the maximum distance
between opposite corners of the cube will be sqrt(3)*L. Hence a slice reconstructed normal to a line connecting these two
corners would have an effective slice sampling of sqrt(3)*L. This orientation-dependent effective slice-thickness
needs to be kept in mind when analyzing such data. This problem is compounded further when the
data set is anisotropic, as are most.
Measurement Of
Gross Morphology/Issues Of Reliability And Validity
MRI provides great flexibility in acquiring a number of different visual
representations of the human brain, each designed to enhance contrast between
different tissue types, highlight abnormalities within a tissue type, or provide
an optimal view of a specific lobar region or subcortical neuroanatomic structure. Techniques for morphometric analysis rely
to varying degrees on user interaction with a visual representation of the brain
image, displayed on film or a computer console, and automated algorithms applied
to the matrix of numerical values representing each pixel in the image.
Validity must be established for both interactive and automated approaches,
though in the absence of satisfactory in vivo validation models, automated approaches
are frequently validated against the "gold standard" of a human expert,
by establishing consistency with known effects such as changes with age, by
comparisons with pathological samples, or by use of phantoms simulating different
tissue types (1). Given the complexity
and individual variability of each brain, issues of reliability acquire a particular
salience when human input is involved, but should also be documented for automated
elements of any analysis procedure (e.g., 8, 105)
Volumetric Vs.
Area Measures
The size of structures of the brain can be estimated using linear, area,
or volumetric measurements. In the early
CT literature, planimetric (area) measures were most commonly used and applied
to one or two "best view" slices.
This approach was accessible to investigators without computerized image
processing facilities, but was vulnerable to sampling error, and sacrificed
volumetric information. With the advent
of MR, and increasing availability of computerized analysis, volumetric analysis
has become the norm. A volumetric approach,
ideally encompassing the entire structure in each hemisphere, is particularly
important for studying left-right asymmetry of structures that are not evenly
aligned relative to the plane of view on a single slice of a given image (112).
This approach can compensate for any bias in selection of imaging plane, or
apparent asymmetries due to head placement.
Segmentation of
CSF, Gray Matter and White Matter
Before brain tissue can be segmented, its extent must be defined. In vivo estimates of brain size are usually based on estimates of
the total volume of brain tissue and CSF, including both ventricular and subarachnoid
fluid. Non-brain material must be identified or removed from the image before
the brain as a whole can be measured and its external features identified.
Correct identification of the outer brain/dural boundary involves certain
assumptions and estimations and is complicated by the fact that bone can easily
be confused with CSF and tissue in many MRI acquisition sequences, especially
T1 sequences designed to enhance contrast between white matter and gray matter.
Such ambiguities may underlie the recently observed variation in estimated
brain size between a number of MRI studies (71). One solution to the problem is to use late
echo images in which CSF, with its long T2, can be sharply differentiated from
bone, for identifying subarachnoid fluid margins and then transfer information
about these boundary locations to images obtained with sequence parameters on
which CSF/tissue and white/gray matter contrasts are optimized (51). Although brain tissue volume is the principal
object of study in the effort to identify neuronal substrates of dysfunction,
the size of CSF-filled spaces relative to the size of the intracranial vault
can also be instructive for interpretation regarding degeneration of tissue
and progression of disease.
Tissue segmentation routines to differentiate CSF from tissue, and subdivide tissue into gray matter and white matter all have limitations because there is no discrete image intensity signature, even across echo times, for each tissue type. Segmentation is also complicated both by the presence of areas of white matter hyperintensity (WMHI), and by the fact that individual voxels may be partially volumed and thus contain mixtures of CSF and gray matter or gray matter and white matter. Approaches to this analysis include edge detection, application of various forms of filters, and fuzzy clustering, each of which can be implemented with varying degrees of operator interaction. Edge detection, for example, of the outer margin of the brain, the ventricles, white matter tracts, or specific subcortical structures such as the hippocampus or striatal structures (89) can be done by a human tracing outlines on a computer display or by using algorithms which identify pixels whose value differs sufficiently from neighboring pixels to qualify them as a boundary pixels (7). Such automated edge detection works well when contrast is high and structures to be outlined are of a relatively simple shape. Operator input frequently supplements automated outlining.
Thresholding involves classifying pixels according to image intensity as representative of CSF or different tissue types. Some approaches take advantage of the different contrasts obtained by different pulse sequences or echo times, either combining data from T1, T2 and proton density weighted images (e.g., 1) or subtracting one image from another to enhance CSF-tissue contrast (32, 51). More recently, "fuzzy" cluster or classifier approaches, which adopt a probabilistic model towards each voxel and thus are better able to account for partial voluming effects have been found to have high reliability, accuracy and validity (33, 82). Automated assessments may employ discriminant functions based on ranges of pixel values representing known tissue and fluid samples (e.g., 14, 17).
One interactive approach requires an operator to adjust a 2-bit image until a pixel value is identified which correctly differentiates CSF from tissue (51). Another takes advantage of the fact that while losses in signal differentiation, for whatever cause, make automated segmentation more difficult, the human visual system, which is quite robust to local intensity changes, can often compensate for them—people can see the gray-white tissue differences even when the segmentation routines fail. In this situation, stereological point counting methods (11, 28, 42) might be superior. This approach involves the random placement of a grid with sufficient resolution (in two or three dimensions) over the structure of interest and counting the points overlying the region or structure of interest. The only requirements are that the grid completely encompass the region or structure, that it be placed on the structure randomly, and that an adequate number of points (150-200) be counted on an adequate number of slices (~10). The approach can be very efficient and is statistically sound. The stereological point counting approach has the additional advantage of providing a coefficient of error of the measurement of the volume of the structure of interest which estimates the magnitude of difference that can be determined between subject groups.
Models for tissue segmentation which work well on healthy young brains may fail when applied to aging subjects or populations with greater tissue pathologies, particularly of the white matter (e.g., 73). Furthermore, assessments of subcortical gray matter volume using thresholding algorithms can be affected by variations in tissue iron content which shortens T2, resulting in signal loss. Iron content-related changes in T2 occur with normal aging and in basal ganglia disorders, particularly in the substantia nigra and globus pallidus (78). The misclassification of WMHI as gray matter can be handled by visual checking and manual correction (e.g. 75) in samples where this occurs infrequently. Clinical populations where focal abnormalities such as gliosis, edema, demyelination, and ischemia, processes that alter the properties of tissue water and produce changes in white and gray matter signal intensity, are common require new approaches to tissue segmentation. Magnetization transfer ratio is one such approach, specifically designed to characterize white matter pathology (96, 103).
Defining Anatomical
Structures And Regions Of Interest For Measurement
As in image acquisition, where tradeoffs need to be made between time in scanner
and optimizing image resolution, tissue contrast, and range of structures visualized,
so in image analysis there are tradeoffs to be made between labor intensive,
interactive delineation of boundaries, and automated approaches, whether using
geometrically based boundary drawing algorithms (e.g., 111), parcellation against
a standardized set of coordinates (e.g., 4), or template matching (98).
As research laboratories collect ever increasing numbers of brains for
analysis, automated approaches which can be implemented quickly and with less
ongoing neuroanatomic expertise become very attractive.
The jury is still out on whether the advantages of speed and reliability
of such systems offset inaccuracies introduced by the deformations involved
in a "one size fits all" measurement approach.
Approaches combining automation with expert supervision attempt to get
the best of both worlds.
One approach to automated measurement of major brain regions uses an MRI-derived anatomic atlas, based on a stereotactic coordinate system (95). Individual 3D brains are linearly transformed to fit neuroanatomic anchor points from the atlas and volumes of ROIs are computed by summing brain voxels coinciding with the atlas boxes assigned to each ROI (e.g., 4). This approach has been extended by thresholding brains into CSF, white, and gray matter to delineate subcortical gray matter structures as well as overall cerebrum, cerebellum and ventricular volumes (40) and validated against interactively defined gold standards and by demonstrating predicted gender, age and group effects. For a more comprehensive treatment the variety of brain atlases available, and different approaches to deformations the reader is referred to a recent specialized publication (102).
While much can be inferred from global measures of overall brain volume and volume of different tissue compartments, specific regions particularly the frontal and temporal lobes as well as limbic and striatal structures are believed to be of particular interest in characterizing the pathophysiology of psychiatric disorders (12). New capabilities for 3D display of brain images have greatly expanded options for delineating brain regions according to sulcal and gyral landmarks, but also can confound the investigator with the complexity and variability of the human brain. For example, while the central sulcus provides the boundary separating frontal lobes from parietal lobes, its configuration varies between individuals and between each cerebral hemisphere. One of the earliest MRI estimates of frontal lobe volume was based on manual tracing of a single midsaggital slice (2), yielding a single measure representing a fraction of the entire structure. Now investigators are combining surface sulcal landmarks (79) with multiple, orthogonal views of brain tissue to identify and trace boundaries of specific cortical gyri (11, 108), although these are not the same as functionally homogenous cytoarchitectonic brain regions. Among the subregions defined and segmented into white and gray matter are superior, medial, inferior, orbital, insular and cingulate regions. This painstaking approach requires interactive neuroanatomic judgments, considerably speeded by visual enhancement of sulcal boundaries, and simultaneous display of internal images in different orthogonal planes. A recent application of this technique to the prefrontal cortex in schizophrenia found a gray matter deficit only in the inferior region (11).
Rules for defining
boundaries for temporal lobe as a whole, and its subsections such as the superior
temporal gyrus (STG) and planum temporale have been constrained in many cases
by the resolution and plane of view available to the investigator, and have
not been comparable across laboratories. Increasing availability of high resolution 3D acquisition protocols,
and interactive 3D image processing systems will enable closer approximation
to established anatomical criteria, and hopefully the adoption of standardized
criteria. Absence of standard definitions
for mesial temporal lobe structures, particularly hippocampus and amygdala,
has been particularly evident in MRI studies of these regions in patients with
schizophrenia (for review, see 35). Meta analysis, which provides greater power
to detect effects submerged in varying sources of error in individual studies,
indicates that mesial temporal structures are reduced in volume in schizophrenia,
whether or not the amygdala is included in the measure (66). The intensely demanding requirements for interactive
delineation of hippocampal structures, as well as the error inherent in subjective
assessment of structure boundaries motivated development of an automated approach,
that measures not only volume (31) but also characterizes shape (30).
Application of multi-dimensional transformations based on pattern matching
has identified specific shape deformations of the hippocampus in patients with
schizophrenia localized to regions of this structure that send projections to
the frontal cortex (19). Such analysis
was not only more sensitive in differentiating patients from controls than simple
volumetric measurements, but also identified an abnormality consistent with
the hypothesized disconnection of hippocampal-prefrontal circuitry.
The development and application of laboratory-specific boundaries, criteria, and procedures for defining neuroanatomic regions will lead to inconsistencies rather than convergence in the literature.
Assessment of
Shape and Complex Morphological Features.
New analyses incorporating advances in physics and mathematics are now enabling
examination of both focal and global characteristics of a structure that are
independent of variations in size, and can be performed with varying degrees
of automaticity. Until recently, a 3D representation of brain
structures could only be visualized in the investigator's imagination. Now 3D displays of the external surfaces, with
ability to slice through selected planes are routinely available in most brain
imaging laboratories. Creation of 3D
images employs a range of rapidly evolving computer graphics software and hardware
methods (23), surface and volume rendering
techniques. While visually compelling,
rendered images may contain distortions introduced by the viewing model and
image preprocessing that must be taken into account before quantitative analysis
is made of the rendered image.
Neurosurgeons and
neurologists have lead the way in applying these techniques for localizing lesions
and vasculature before neurosurgery (e.g., 20), developing neuroanatomic teaching
tools enabling students to explore the structure and function of the brain in
virtual reality (e.g., 99), characterizing the normal surface anatomy of the
brain (104), and identifying specific neuroanatomic regions based on specific
sulcal and gyral landmarks (79). Fluid deformations of one exterior brain surface
against another (template) followed by encoding and mapping of probability of
displacement are being used to assess individual and group differences in the
immensely complex and variable gyral and sulcal topography (98).
These probability statistics expand on earlier concepts of anatomic atlases
by simultaneously characterizing both locational and geometric variation and
open up a new dimension in characterizing and quantifying regional abnormality
relative to age and gender norms. Other structures for which shape analyses have
uncovered specific dysmorphologies in schizophrenia include the hippocampus
(19) and the corpus callosum (16, 97).
Co-Registration
With Other In Vivo Neuroimages
The relatively poor in-plane resolution and wide slices of images from functional
imaging techniques such as positron emission tomography (PET) and proton or
phosphorus magnetic spectroscopic imaging (MRSI) limits anatomic localization. Structural information from MR images is now
routinely used for localizing and displaying high intensity areas of activation
derived from functional imaging modalities including PET (67, 91) and functional
MRI (fMRI) (15). Indeed, this localization
of functional data is important not only for more precisely characterizing the
structures involved, but also to take regional or group differences in CSF,
gray and white matter composition of the partially volumed functional voxels
into account in interpreting results. Since voxels obtained in MRS or PET tend to be larger than those
now available with MRI, co-registration of MR and PET or MRS images allows voxel
composition information from structural images to be applied to functional images.
This approach has been used for PET (59, 64), and recently also for MRSI
(34, 49, 75). Such corrections are computationally challenging,
but eventually will allow much greater precision in assessing functional activity.
The simplest case of co-registration is when both data sets have been collected on the same imaging hardware, without moving the subject, as is typically the case with fMRI and MRSI scanning (Figure 3). For PET scanning, differences in hardware and time of scan are unavoidable. Several different approaches to image acquisition to facilitate subsequent matching of acquired PET and MR images have been developed. These include use of external landmarks, or fiducials, consistently placed and of a size and property to be adequately imaged by both modalities (22), and head holders to stabilize and standardize head position between scanners (41). Varying levels of sophistication have been applied to compute a transformation matrix between two sets of images, including iterative matching of brain surfaces (48, 70) projections (93) or overall brain tissue (25, 110). The automated image registration (AIR) (110) and statistical parametric mapping (SPM) (25) systems are now the most commonly used. They have recently been cross-validated and their reliability established (43)
A particular application for merging of data sets within an imaging modality is the generation of averaged brains (3, 109). This approach involves aligning and matching successive images relative to stereotactic coordinates (21) and provides a useful visual summary of group differences. Regions of statistically significant differences can also be identified, but should be verified by a separate case by case quantification of the structure or region in question (72). Image averaging works better for assessing subcortical and midline structures, less affected by deformations necessary to standardize brains, than cortical structures.
go to Figure 4
INTERPRETATION OF MAGNETIC RESONANCE BRAIN IMAGES
Neuroimaging findings in some mental disorders can be elusive. The true effect is small, and brain morphology is strongly affected by several variables which are relatively independent of disease, such as age, sex and somatic size, as well as those which might be associated with mental disorder, e.g., alcohol consumption, socio-economic status and handedness. These variables must be taken into consideration and appropriate controls implemented before inferences about the pathophysiology of various mental disorders can be properly drawn.
Selection of Appropriate
Controls
When CT was the prevailing technology, medical patients with "normal"
scans were frequently used as controls in imaging studies of psychiatric patients
- a practice motivated, to some extent, by reluctance to expose healthy subjects
to unnecessary radiation. Whether or
not this practice contributed to observed group differences was never fully
resolved (81, 92). With MRI scanning,
this radiation concern is not relevant. However,
the prospective recruitment of control samples for any neuroimaging studies
of psychiatric disorder requires attention not only to excluding cases with
psychiatric disorders (90), but also to including subjects who match the parental
socioeconomic status, educational and ethnic status of the patients (84).
In addition to comparisons with normative data, comparisons between pathologic
groups are also valuable, especially between different diagnostic groups with
documented abnormalities in common anatomic region(s).
For example, comparisons between patients with temporal lobe epilepsy
and schizophrenia are appropriate for studies investigating the role of limbic
structures in these disorders (55).
Accounting for
Normal Developmental and Aging Effects
Schizophrenia is typically a disease with onset in late adolescence or early
adulthood, a period during which cortical gray matter undergoes significant
reduction, but cortical white continues to increase in volume (74).
Its course coincides with a period of stability of white matter volume,
but gradual reduction in gray matter and increase in CSF (37, 65, 74, 83) (Figure
5). These normal developmental and
aging effects need to be taken into account when interpreting observed brain
dysmorphology in schizophrenia. Studies
of Alzheimer's disease (AD), Parkinson's disease (PD) and other diseases which
typically make their appearance after middle age, likewise need to take into
account not only normal age-related changes occurring at this point in the life
span, but also the greater variability in brain tissue volumes found among the
healthy elderly than the healthy young (37, 38).
Accounting For
Normal Variations In Somatic Size
Given the existence of a wide range of body and head sizes among the population,
it is important to demonstrate that any size difference of a specific brain
structure is independent of such non-disease-relevant differences. In general, large people will have larger heads
and larger brains than small people, although such relationships have been easier
to demonstrate across gender and between racial groups than within them (71).
In most research contexts, this variability constitutes a form of noise
which investigators need to minimize in order to compare the volume of a particular
brain region or structure between individuals and groups.
Thus head size correction seems particularly appropriate when gender
differences in brain morphology are under investigation (18, 45).
One approach for head-size correction has been use of a ratio, such as the ventricle brain ratio (VBR), which was developed to provide an index of ventricular enlargement relative to the size of the total brain. Other correction approaches have included proportions, where the size of a structure, or the amount of tissue or CSF in a ROI, is divided by an estimate of the intracranial volume (e.g., 37); or regression analysis which partials out the components of quantitative brain measures attributable to headsize variation (e.g., 74).
Concern that the ratio or proportion approaches to head size correction could introduce error or obscure the complexity of underlying relationships has lead some investigators to warn against their use, and advocate regression or multifactorial approaches, as necessary, to parse out contribution of non-specific variables (6). Others (57) have demonstrated that a regression approach to head-size correction removes irrelevant true-score variance which reduces reliability and improves the correlation with validity criteria such as age and diagnostic status. Linear regression of specific ROIs against estimated head-size in a sample of normal control subjects highlighted the fact that different parts of the brain bear different relationships to overall brain size and suggests that ROI-specific correction factors should be applied (57).
While head-size correction is appropriate in some contexts, in others, such as developmental studies (74), it may not be. For such studies, absolute volumes of particular brain regions rather than head-size corrected values should be used. Because investigators have described reduced intracranial volume in patients with schizophrenia (69), the issue has become more complicated. If schizophrenia affects the brain before it has completed its growth, leading to diminished brain size and as a result diminished skull size, the correction of all structural brain measures for intracranial volume could obscure pathologically relevant information. Regardless of possible disease-related reduction in head size, however, an investigator asking whether particular brain structures are disproportionately reduced, must examine and take into account differences in intracranial volume.
Accounting For
Widespread Tissue Deficits
Generalized morphological differences in the brains of psychiatric patients
compared to controls should be taken into account before assessing smaller structures. For instance, schizophrenia appears to be characterized
by generalized ventricular and cortical sulcal enlargement as well as a widespread
loss of cortical gray matter (88). Therefore, it is necessary to consider whether
loss of tissue from the target region exceeds that found from adjacent or control
regions before concluding that a specific brain structure or region is uniquely
affected by the disease (e.g., 112). Using such an analysis, recent studies have indeed identified certain
regions of frontal and parieto-temporal cortices (86, 94) as particularly affected
in schizophrenia.
As neuroimaging matures and becomes a standard element of psychiatric clinical research, longitudinal studies of change over time will increasingly replace the inferences made from cross-sectional studies. Direct longitudinal comparisons over time are hampered by major advances in technology (e.g. replacement of CT with MRI) as well as system upgrades within a technology. However, the advantages of direct longitudinal evaluation of the progression of a disease over inference from cross sectional studies motivates the development and application of procedures both to standardize collection parameters over time and to estimate and correct for the error introduced by uncontrolled factors.
Within psychiatry, longitudinal neuroimaging studies are beginning to show promise for studying the progression of the disease (29, 56, 76), ideally with a matched control group also followed for comparable intervals. Neuroimaging studies of multiple sclerosis, where the waxing and waning of white matter lesions are regularly monitored provide some important pointers for longitudinal studies of psychiatric patients (62). Studies of multiple sclerosis were among the first to alert investigators to the influence of measurement error due to misalignment between scans over time on measures of real change (26). Other subtle sources of error between tests identified by MS researchers include variation in scanner sensitivity (106) and voxel dimensions (47). The use of phantom standards in all scans greatly facilitates application of post hoc adjustments (101).
Standardizing
Measures In Face Of Changing Technology
As neuroimaging methods advance and new imaging techniques are used, it is necessary
to understand how morphometric methods are affected by differences in imaging
methods. A major transition in the
past 15 years has been the replacement of CT scanning with MRI scanning. While MRI adds much information, particularly
from its enhanced soft-tissue resolution, to that provided by CT, it also provides
similar information on CSF-filled spaces, particularly ventricular enlargement,
which has been a hallmark of CT studies. Thus MRI values can potentially offer continuity
for longitudinal studies in which a transition was made from CT to MRI.
A comparison of CT and MRI-derived estimates of ventricular and sulcal
CSF volume obtained within two weeks of each other in patients meeting RDC criteria
for alcohol dependence (77) found that while CT and MRI estimates of absolute
ventricular volume were equivalent, CT considerably underestimated cortical
sulcal volume relative to MRI. This
underestimation could be due to poorer resolution and spectral shift artifact
on CT. Thus longitudinal studies of subjects initially scanned with CT and subsequently
with MRI may, with suitable caution, assess changes in ventricular volume, but
should not compare sulcal volume.
The increasing replacement of conventional 2D image acquisition sequences with fast 2D and 3D acquisition sequences makes issues of comparability across acquisition sequence within MRI extremely relevant. While comparisons between MRI sequences for sensitivity to detecting various types of lesions are common (e.g., 100) few comparisons between regional volumes of CSF, white and gray matter obtained by different sequences have been published. One pilot study, tested the comparability of CSF, GM, and WM measures made from 3D SPGR images to those acquired using a T2-weighted axial spin-echo protocol (54) and found that although there was a high degree of correlation between the CSF, white and gray matter values from both sets of images, there was relatively greater white matter and less CSF and gray matter in SPGR compared to spin echo images, suggesting a systematic difference between the methods. Thus longitudinal comparisons of data collected using different acquisition methods should be made with caution since the differences in tissue contrast, signal-to-noise ratios, and slice profiles will affect the results of morphometric measures.
Longitudinal Measurement
Error
Even when software protocols are identical, differences in head positioning
between scans can contribute measurement error in studies involving repeat scans
(26), particularly when 2D, rather than 3D images are being compared. A longitudinal CT study on patients with AD
identified apparent intra-individual change in total cranial volume between
first and second scans ranging from -60.5 to 89.5 cc, which correlated significantly,
but differentially with several of the ROI change scores. Since total cranial
volume (CSF plus tissue) should not change between scans in adults, this apparent
change was used to derive a specific correction factor for each ROI change (87).
This procedure has since been applied to longitudinal studies of patients
with alcoholism (76) and schizophrenia (56).
Other approaches to correcting for differences in head positioning between
scans include realignment of scans in three axes and reslicing along the AC-PC
line (29). These efforts to correct
measurement error due to differences in head placement have only recently been
applied in longitudinal studies, and the lack of such corrections in earlier
longitudinal studies of brain morphological change in patients with schizophrenia
may account for inconsistent results.
Other sources of error include unreliability in measurement techniques (36) as well as changing approaches to measurement. Problems limiting the value of earlier neuroimaging follow-up studies can now be substantially reduced by applying the same measurement approach to baseline as well as follow-up data sets, assessing inter and intra-rater reliabilities, and retesting control groups along with patient groups.
The superior resolution and flexibility of MR imaging has enabled the examination of more detailed and specific anatomical questions in psychiatric disorders. New developments in acquisition and image processing software continue to expand the possibilities available to investigators, not only for visual presentation, but also for quantitative analyses. While there continue to be important methodological and technical issues to attend to, investigators are becoming increasingly sophisticated in using appropriate techniques, and adopting strategies to ensure the reliability and validity of the data being produced. Standardization of techniques between laboratories is becoming increasingly viable though the dissemination of image analysis software and posting of brain atlases in web sites. It continues to be important that data are collected and analyzed in such a way that specific hypotheses about particular regions of the brain can be tested within the context of the brain as a whole.
published 2000