How to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index.

Cortical folding (gyrification) is determined during the first months of life, so that adverse events occurring during this period leave traces that will be identifiable at any age. As recently reviewed by Mangin and colleagues(2), several methods exist to quantify different characteristics of gyrification. For instance, sulcal morphometry can be used to measure shape descriptors such as the depth, length or indices of inter-hemispheric asymmetry(3). These geometrical properties have the advantage of being easy to interpret. However, sulcal morphometry tightly relies on the accurate identification of a given set of sulci and hence provides a fragmented description of gyrification. A more fine-grained quantification of gyrification can be achieved with curvature-based measurements, where smoothed absolute mean curvature is typically computed at thousands of points over the cortical surface(4). The curvature is however not straightforward to comprehend, as it remains unclear if there is any direct relationship between the curvedness and a biologically meaningful correlate such as cortical volume or surface. To address the diverse issues raised by the measurement of cortical folding, we previously developed an algorithm to quantify local gyrification with an exquisite spatial resolution and of simple interpretation. Our method is inspired of the Gyrification Index(5), a method originally used in comparative neuroanatomy to evaluate the cortical folding differences across species. In our implementation, which we name local Gyrification Index (lGI(1)), we measure the amount of cortex buried within the sulcal folds as compared with the amount of visible cortex in circular regions of interest. Given that the cortex grows primarily through radial expansion(6), our method was specifically designed to identify early defects of cortical development. In this article, we detail the computation of local Gyrification Index, which is now freely distributed as a part of the FreeSurfer Software (http://surfer.nmr.mgh.harvard.edu/, Martinos Center for Biomedical Imaging, Massachusetts General Hospital). FreeSurfer provides a set of automated reconstruction tools of the brain's cortical surface from structural MRI data. The cortical surface extracted in the native space of the images with sub-millimeter accuracy is then further used for the creation of an outer surface, which will serve as a basis for the lGI calculation. A circular region of interest is then delineated on the outer surface, and its corresponding region of interest on the cortical surface is identified using a matching algorithm as described in our validation study(1). This process is repeatedly iterated with largely overlapping regions of interest, resulting in cortical maps of gyrification for subsequent statistical comparisons (Fig. 1). Of note, another measurement of local gyrification with a similar inspiration was proposed by Toro and colleagues(7), where the folding index at each point is computed as the ratio of the cortical area contained in a sphere divided by the area of a disc with the same radius. The two implementations differ in that the one by Toro et al. is based on Euclidian distances and thus considers discontinuous patches of cortical area, whereas ours uses a strict geodesic algorithm and include only the continuous patch of cortical area opening at the brain surface in a circular region of interest.

Design, data management, and population baseline characteristics of the PERFORM magnetic resonance imaging project.

Quantitative information from magnetic resonance imaging (MRI) may substantiate clinical findings and provide additional insight into the mechanism of clinical interventions in therapeutic stroke trials. The PERFORM study is exploring the efficacy of terutroban versus aspirin for secondary prevention in patients with a history of ischemic stroke. We report on the design of an exploratory longitudinal MRI follow-up study that was performed in a subgroup of the PERFORM trial. An international multi-centre longitudinal follow-up MRI study was designed for different MR systems employing safety and efficacy readouts: new T2 lesions, new DWI lesions, whole brain volume change, hippocampal volume change, changes in tissue microstructure as depicted by mean diffusivity and fractional anisotropy, vessel patency on MR angiography, and the presence of and development of new microbleeds. A total of 1,056 patients (men and women ≥ 55 years) were included. The data analysis included 3D reformation, image registration of different contrasts, tissue segmentation, and automated lesion detection. This large international multi-centre study demonstrates how new MRI readouts can be used to provide key information on the evolution of cerebral tissue lesions and within the macrovasculature after atherothrombotic stroke in a large sample of patients.

Identification of different heart tissues from MRI C-SENC images using an unsupervised multi-stage fuzzy clustering technique.

PURPOSE: To objectively characterize different heart tissues from functional and viability images provided by composite-strain-encoding (C-SENC) MRI. MATERIALS AND METHODS: C-SENC is a new MRI technique for simultaneously acquiring cardiac functional and viability images. In this work, an unsupervised multi-stage fuzzy clustering method is proposed to identify different heart tissues in the C-SENC images. The method is based on sequential application of the fuzzy c-means (FCM) and iterative self-organizing data (ISODATA) clustering algorithms. The proposed method is tested on simulated heart images and on images from nine patients with and without myocardial infarction (MI). The resulting clustered images are compared with MRI delayed-enhancement (DE) viability images for determining MI. Also, Bland-Altman analysis is conducted between the two methods. RESULTS: Normal myocardium, infarcted myocardium, and blood are correctly identified using the proposed method. The clustered images correctly identified 90 +/- 4% of the pixels defined as infarct in the DE images. In addition, 89 +/- 5% of the pixels defined as infarct in the clustered images were also defined as infarct in DE images. The Bland-Altman results show no bias between the two methods in identifying MI. CONCLUSION: The proposed technique allows for objectively identifying divergent heart tissues, which would be potentially important for clinical decision-making in patients with MI.

ZHARP: three-dimensional motion tracking from a single image plane.

Three-dimensional imaging and quantification of myocardial function are essential steps in the evaluation of cardiac disease. We propose a tagged magnetic resonance imaging methodology called zHARP that encodes and automatically tracks myocardial displacement in three dimensions. Unlike other motion encoding techniques, zHARP encodes both in-plane and through-plane motion in a single image plane without affecting the acquisition speed. Postprocessing unravels this encoding in order to directly track the 3-D displacement of every point within the image plane throughout an entire image sequence. Experimental results include a phantom validation experiment, which compares zHARP to phase contrast imaging, and an in vivo study of a normal human volunteer. Results demonstrate that the simultaneous extraction of in-plane and through-plane displacements from tagged images is feasible.

Ligamentous and tendinous anatomy of the intermetacarpal and common carpometacarpal joints: evaluation with MR imaging and MR arthrography.

PURPOSE: The purpose of this work was to demonstrate the normal ligamentous and tendinous anatomy of the intermetacarpal (IMC) and common carpometacarpal (CCMC) joints with MRI and MR arthrography. METHOD: MR images of 22 wrists derived from fresh human cadavers were obtained before and after arthrography. The MR imaging features of the ligaments and tendons about the CCMC and IMC joints and the joints themselves were analyzed in a randomized fashion and correlated with those seen on anatomic sections. RESULTS: Six CCMC ligaments were visualized. The dorsal and palmar CCMC ligaments and the pisometacarpal ligament were best visualized in the sagittal plane. The radial and ulnar CCMC collateral ligaments and the capito-third metacarpal ligament were best visualized in the coronal plane. Three main IMC ligaments were observed: a dorsal and a palmar ligament and an interosseous ligament complex. All three ligaments were best visualized in the axial plane. Four tendinous insertions to the metacarpal bases were evident. CONCLUSION: The anatomy of the ligaments and tendinous insertions about the second to fifth IMC and the CCMC joints is well demonstrated by MR imaging and MR arthrography. MR arthrography does not significantly improve the visualization of these complex structures.

SA2RAGE: a new sequence for fast B1+ -mapping.

At high magnetic field strengths (≥ 3T), the radiofrequency wavelength used in MRI is of the same order of magnitude of (or smaller than) the typical sample size, making transmit magnetic field (B1+) inhomogeneities more prominent. Methods such as radiofrequency-shimming and transmit SENSE have been proposed to mitigate these undesirable effects. A prerequisite for such approaches is an accurate and rapid characterization of the B1+ field in the organ of interest. In this work, a new phase-sensitive three-dimensional B1+-mapping technique is introduced that allows the acquisition of a 64 × 64 × 8 B1+-map in ≈ 20 s, yielding an accurate mapping of the relative B1+ with a 10-fold dynamic range (0.2-2 times the nominal B1+). Moreover, the predominant use of low flip angle excitations in the presented sequence minimizes specific absorption rate, which is an important asset for in vivo B1+-shimming procedures at high magnetic fields. The proposed methodology was validated in phantom experiments and demonstrated good results in phantom and human B1+-shimming using an 8-channel transmit-receive array.

Is there histomorphological evidence of plantar metatarsal fat pad atrophy in patients with diabetes?

The etiology of diabetic foot ulceration remains incompletely understood. Among other factors such as foot deformity in the presence of neuropathy, plantar fat pad atrophy has been identified as a contributory factor in diabetic foot ulceration. An association between fat pad atrophy and diabetic foot ulceration has been documented by imaging and histomorphological analysis of the calcaneal fat pad. However, histomorphological analysis of the metatarsal fat pad has not been performed to date. The present study entailed 14 patients with diabetes and 14 nondiabetic controls and was aimed at documenting histomorphological evidence for presumed plantar metatarsal fat pad atrophy in patients with diabetes. Histological stains and computer-assisted planimetry were performed on samples of metatarsal fat obtained during forefoot surgery. The histomorphological and planimetric analyses of adipocyte cross-sectional area and nuclear density demonstrated no differences between patients with diabetes and control patients. Our findings demonstrate that systemic atrophy of the metatarsal fat pad is not present in the diabetic foot and may not explain the structural changes previously proposed by noninvasive imaging. Level of Clinical Evidence: 3.

Official Positions for FRAX® clinical regarding smoking from Joint Official Positions Development Conference of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX®.

The worldwide prevalence of smoking has been estimated at about 50% in men, and 10% in women, with larger variations among different populations studied. Smoking has been shown to affect many organ systems resulting in severe morbidity and increased mortality. In addition, smoking has been identified as a predictor of ten-year fracture risk in men and women, largely independent of an individual's bone mineral density. This finding has eventually lead to incorporation of this risk factor into FRAX®, an algorithm that has been developed to calculate an individual's ten-year fracture risk. However, only little, or conflicting data is available on a possible association between smoking dose, duration, length of time after cessation, type of tobacco and fracture risk, limiting this risk factor's applicability in the context of FRAX®.

A comparison of in vivo 13C MR brain glycogen quantification at 9.4 and 14.1 T.

The high molecular weight and low concentration of brain glycogen render its noninvasive quantification challenging. Therefore, the precision increase of the quantification by localized (13) C MR at 9.4 to 14.1 T was investigated. Signal-to-noise ratio increased by 66%, slightly offset by a T(1) increase of 332 ± 15 to 521 ± 34 ms. Isotopic enrichment after long-term (13) C administration was comparable (≈ 40%) as was the nominal linewidth of glycogen C1 (≈ 50 Hz). Among the factors that contributed to the 66% observed increase in signal-to-noise ratio, the T(1) relaxation time impacted the effective signal-to-noise ratio by only 10% at a repetition time = 1 s. The signal-to-noise ratio increase together with the larger spectral dispersion at 14.1 T resulted in a better defined baseline, which allowed for more accurate fitting. Quantified glycogen concentrations were 5.8 ± 0.9 mM at 9.4 T and 6.0 ± 0.4 mM at 14.1 T; the decreased standard deviation demonstrates the compounded effect of increased magnetization and improved baseline on the precision of glycogen quantification.

Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy.

We have developed a digital holographic microscope (DHM), in a transmission mode, especially dedicated to the quantitative visualization of phase objects such as living cells. The method is based on an original numerical algorithm presented in detail elsewhere [Cuche et al., Appl. Opt. 38, 6994 (1999)]. DHM images of living cells in culture are shown for what is to our knowledge the first time. They represent the distribution of the optical path length over the cell, which has been measured with subwavelength accuracy. These DHM images are compared with those obtained by use of the widely used phase contrast and Nomarski differential interference contrast techniques.

Preliminary report on in vivo coronary MRA at 3 Tesla in humans.

Current limitations of coronary magnetic resonance angiography (MRA) include a suboptimal signal-to-noise ratio (SNR), which limits spatial resolution and the ability to visualize distal and branch vessel coronary segments. Improved SNR is expected at higher field strengths, which may provide improved spatial resolution. However, a number of potential adverse effects on image quality have been reported at higher field strengths. The limited availability of high-field systems equipped with cardiac-specific hardware and software has previously precluded successful in vivo human high-field coronary MRA data acquisition. In the present study we investigated the feasibility of human coronary MRA at 3.0 T in vivo. The first results obtained in nine healthy adult subjects are presented.

Joint Official Positions of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX(®). Executive Summary of the 2010 Position Development Conference on Interpretation and use of FRAX® in clinical practice.

The International Society for Clinical Densitometry (ISCD) and the International Osteoporosis Foundation (IOF) convened the FRAX(®) Position Development Conference (PDC) in Bucharest, Romania, on November 14, 2010, following a two-day joint meeting of the ISCD and IOF on the "Interpretation and Use of FRAX(®) in Clinical Practice." These three days of critical discussion and debate, led by a panel of international experts from the ISCD, IOF and dedicated task forces, have clarified a number of important issues pertaining to the interpretation and implementation of FRAX(®) in clinical practice. The Official Positions resulting from the PDC are intended to enhance the quality and clinical utility of fracture risk assessment worldwide. Since the field of skeletal assessment is still evolving rapidly, some clinically important issues addressed at the PDCs are not associated with robust medical evidence. Accordingly, some Official Positions are based largely on expert opinion. Despite limitations inherent in such a process, the ISCD and IOF believe it is important to provide clinicians and technologists with the best distillation of current knowledge in the discipline of bone densitometry and provide an important focus for the scientific community to consider. This report describes the methodology and results of the ISCD-IOF PDC dedicated to FRAX(®).

Relaxivity of Gd-based contrast agents on X nuclei with long intrinsic relaxation times in aqueous solutions.

The relaxivity of commercially available gadolinium (Gd)-based contrast agents was studied for X-nuclei resonances with long intrinsic relaxation times ranging from 6 s to several hundred seconds. Omniscan in pure 13C formic acid had a relaxivity of 2.9 mM(-1) s(-1), whereas its relaxivity on glutamate C1 and C5 in aqueous solution was approximately 0.5 mM(-1) s(-1). Both relaxivities allow the preparation of solutions with a predetermined short T1 and suggest that in vitro substantial sensitivity gains in their measurement can be achieved. 6Li has a long intrinsic relaxation time, on the order of several minutes, which was strongly affected by the contrast agents. Relaxivity ranged from approximately 0.1 mM(-1) s(-1) for Omniscan to 0.3 for Magnevist, whereas the relaxivity of Gd-DOTP was at 11 mM(-1) s(-1), which is two orders of magnitude higher. Overall, these experiments suggest that the presence of 0.1- to 10-microM contrast agents should be detectable, provided sufficient sensitivity is available, such as that afforded by hyperpolarization, recently introduced to in vivo imaging.

Prothèses totales de genoux. Analyse clinique et numérique des micromouvements de l'implant tibial [Total knee prosthesis. Clinical and numerical study of micromovements of the tibial implant].

INTRODUCTION: The importance of the micromovements in the mechanism of aseptic loosening is clinically difficult to evaluate. To complete the analysis of a series of total knee arthroplasties (TKA), we used a tridimensional numerical model to study the micromovements of the tibial implant. MATERIAL AND METHODS: Fifty one patients (with 57 cemented Porous Coated Anatomic TKAs) were reviewed (mean follow-up 4.5 year). Radiolucency at the tibial bone-cement interface was sought on the AP radiographs and divided in 7 areas. The distribution of the radiolucency was then correlated with the axis of the lower limb as measured on the orthoradiograms. The tridimensional numerical model is based on the finite element method. It allowed the measurement of the cemented prosthetic tibial implant's displacements and the micromovements generated at bone-ciment interface. A total load (2000 Newton) was applied at first vertically and asymetrically on the tibial plateau, thereby simulating an axial deviation of the lower limbs. The vector's posterior inclination then permitted the addition of a tangential component to the axial load. This type of effort is generated by complex biomechanical phenomena such as knee flexion. RESULTS: 81 per cent of the 57 knees had a radiolucent line of at least 1 mm, at one or more of the tibial cement-epiphysis jonctional areas. The distribution of these lucent lines showed that they came out more frequently at the periphery of the implant. The lucent lines appeared most often under the unloaded margin of the tibial plateau, when axial deviation of lower limbs was present. Numerical simulations showed that asymetrical loading on the tibial plateau induced a subsidence of the loaded margin (0-100 microns) and lifting off at the opposite border (0-70 microns). The postero-anterior tangential component induced an anterior displacement of the tibial implant (160-220 microns), and horizontal micromovements with non homogenous distribution at the bone-ciment interface (28-54 microns). DISCUSSION: Comparison of clinical and numerical results showed a relation between the development of radiolucent lines and the unloading of the tibial implant's margin. The deleterious effect of lower limbs' axial deviation is thereby proven. The irregular distribution of lucent lines under the tibial plateau was similar of the micromovements' repartition at the bone-cement interface when tangential forces were present. A causative relation between the two phenomenaes could not however be established. Numerical simulation is a truly useful method of study; it permits to calculate micromovements which are relative, non homogenous and of very low amplitude. However, comparative clinical studies remain as essential to ensure the credibility of results.

Evidence for segregated and integrative connectivity patterns in the human Basal Ganglia.

Detailed knowledge of the anatomy and connectivity pattern of cortico-basal ganglia circuits is essential to an understanding of abnormal cortical function and pathophysiology associated with a wide range of neurological and neuropsychiatric diseases. We aim to study the spatial extent and topography of human basal ganglia connectivity in vivo. Additionally, we explore at an anatomical level the hypothesis of coexistent segregated and integrative cortico-basal ganglia loops. We use probabilistic tractography on magnetic resonance diffusion weighted imaging data to segment basal ganglia and thalamus in 30 healthy subjects based on their cortical and subcortical projections. We introduce a novel method to define voxel-based connectivity profiles that allow representation of projections from a source to more than one target region. Using this method, we localize specific relay nuclei within predefined functional circuits. We find strong correlation between tractography-based basal ganglia parcellation and anatomical data from previously reported invasive tracing studies in nonhuman primates. Additionally, we show in vivo the anatomical basis of segregated loops and the extent of their overlap in prefrontal, premotor, and motor networks. Our findings in healthy humans support the notion that probabilistic diffusion tractography can be used to parcellate subcortical gray matter structures on the basis of their connectivity patterns. The coexistence of clearly segregated and also overlapping connections from cortical sites to basal ganglia subregions is a neuroanatomical correlate of both parallel and integrative networks within them. We believe that this method can be used to examine pathophysiological concepts in a number of basal ganglia-related disorders.