Best individual template selection from deformation tensor minimization

We study the influence of the choice of template in tensor-based morphometry. Using 3D brain MR images from 10 monozygotic twin pairs, we defined a tensor-based distance in the log-Euclidean framework [1] between each image pair in the study. Relative to this metric, twin pairs were found to be closer to each other on average than random pairings, consistent with evidence that brain structure is under strong genetic control. We also computed the intraclass correlation and associated permutation p-value at each voxel for the determinant of the Jacobian matrix of the transformation. The cumulative distribution function (cdf) of the p-values was found at each voxel for each of the templates and compared to the null distribution. Surprisingly, there was very little difference between CDFs of statistics computed from analyses using different templates. As the brain with least log-Euclidean deformation cost, the mean template defined here avoids the blurring caused by creating a synthetic image from a population, and when selected from a large population, avoids bias by being geometrically centered, in a metric that is sensitive enough to anatomical similarity that it can even detect genetic affinity among anatomies.

Corrigendum to "Mapping hippocampal and ventricular change in Alzheimer disease" [NeuroImage 22 (2004) 1754-1766] (DOI:10.1016/j.neuroimage.2004.03.040)

We recently noticed an error in the demographic data in this article. The validity of the findings and the conclusions of the paper is not affected. However, there is an error in the reported sample size and in the means and standard deviations of the subjects’ ages and MMSE scores. We would like to correct this error, which came to light when we were re-analyzing the data for a meta-analysis. The error occurred because an older version of a spreadsheet was incorrectly used when reporting the sample composition. Instead of examining 12 Alzheimer's disease patients and 14 healthy elderly controls, we in fact examined 17 Alzheimer’s disease patients and 14 healthy elderly controls. All maps and morphometric data reported in the paper are correct, except that the sample size was in fact slightly higher than that originally reported, and the maps computed in the paper were based on the larger sample (which included five more subjects in the Alzheimer’s disease group). All of the maps and figures in the paper are correct, and the conclusions of the paper are unchanged. We apologize for this error, which falls under the sole responsibility of the first author. The corrected demographic information appears below.

Brain network efficiency and topology depend on the fiber tracking method: 11 tractography algorithms compared in 536 subjects

As connectivity analyses become more popular, claims are often made about how the brain's anatomical networks depend on age, sex, or disease. It is unclear how results depend on tractography methods used to compute fiber networks. We applied 11 tractography methods to high angular resolution diffusion images of the brain (4-Tesla 105-gradient HARDI) from 536 healthy young adults. We parcellated 70 cortical regions, yielding 70×70 connectivity matrices, encoding fiber density. We computed popular graph theory metrics, including network efficiency, and characteristic path lengths. Both metrics were robust to the number of spherical harmonics used to model diffusion (4th-8th order). Age effects were detected only for networks computed with the probabilistic Hough transform method, which excludes smaller fibers. Sex and total brain volume affected networks measured with deterministic, tensor-based fiber tracking but not with the Hough method. Each tractography method includes different fibers, which affects inferences made about the reconstructed networks.

Analyzing multi-fiber reconstruction in high angular resolution diffusion imaging using the tensor distribution function

High-angular resolution diffusion imaging (HARDI) can reconstruct fiber pathways in the brain with extraordinary detail, identifying anatomical features and connections not seen with conventional MRI. HARDI overcomes several limitations of standard diffusion tensor imaging, which fails to model diffusion correctly in regions where fibers cross or mix. As HARDI can accurately resolve sharp signal peaks in angular space where fibers cross, we studied how many gradients are required in practice to compute accurate orientation density functions, to better understand the tradeoff between longer scanning times and more angular precision. We computed orientation density functions analytically from tensor distribution functions (TDFs) which model the HARDI signal at each point as a unit-mass probability density on the 6D manifold of symmetric positive definite tensors. In simulated two-fiber systems with varying Rician noise, we assessed how many diffusionsensitized gradients were sufficient to (1) accurately resolve the diffusion profile, and (2) measure the exponential isotropy (EI), a TDF-derived measure of fiber integrity that exploits the full multidirectional HARDI signal. At lower SNR, the reconstruction accuracy, measured using the Kullback-Leibler divergence, rapidly increased with additional gradients, and EI estimation accuracy plateaued at around 70 gradients.

The tensor distribution function

Diffusion weighted magnetic resonance imaging is a powerful tool that can be employed to study white matter microstructure by examining the 3D displacement profile of water molecules in brain tissue. By applying diffusion-sensitized gradients along a minimum of six directions, second-order tensors (represented by three-by-three positive definite matrices) can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve more complicated white matter configurations, e.g., crossing fiber tracts. Recently, a number of high-angular resolution schemes with more than six gradient directions have been employed to address this issue. In this article, we introduce the tensor distribution function (TDF), a probability function defined on the space of symmetric positive definite matrices. Using the calculus of variations, we solve the TDF that optimally describes the observed data. Here, fiber crossing is modeled as an ensemble of Gaussian diffusion processes with weights specified by the TDF. Once this optimal TDF is determined, the orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function. Moreover, a tensor orientation distribution function (TOD) may also be derived from the TDF, allowing for the estimation of principal fiber directions and their corresponding eigenvalues.

Now you're talking

Undergraduate medical imaging students at an urban Australian university attend their first clinical placement towards the end of Semester two in their first year. While they often have a great deal of contextual knowledge and targeted theoretical learning prior to their first clinical placement, many of them are unfamiliar with the dynamics of a hospital environment. With the increase in online communication platforms, there was a concern that face-to-face clinical communication skills were becoming a source of anxiety for students, especially prior to their first clinical placements. Given the impact of anxiety on communication, reducing communication-related apprehension was considered a programme priority. A pilot series of scenario-based tutorials was implemented in the general radiography unit for the first year students.

Comparison of silicate-substituted calcium phosphate (Actifuse®) with rhBMP-2 (Infuse®) in posterolateral instrumented lumbar fusion

Study Design This was a randomised controlled trial in patients with degenerative disc disease (DDD) who underwent instrumented posterolateral lumbar fusion (PLF) surgery. Objective The aim of this study was to assess the efficacy of the bone grafting substitute, silicate-substituted calcium phosphate (SiCaP) compared with bone morphogenetic protein (rhBMP-2) and to evaluate clinical outcomes over a period of two years. Methods Patients undergoing PLF surgery for DDD at a single centre were recruited and randomised to one of two groups; SiCaP (n=9) or rhBMP-2 (n=10). One patient withdrew prior to randomisation and another from the rhBMP-2 group after randomisation. The radiological and clinical outcomes were examined and compared. Fusion was assessed at 12 months with computed tomography (CT) and plain radiographs. Clinical outcomes were evaluated by recording measures of pain, quality of life, disability and neurological status from six weeks to two years postoperatively. Results In the SiCaP and rhBMP-2 groups, fusion was observed in 9/9 and 8/9 patients respectively. Pain and disability scores were reduced and quality of life increased in both groups. Leg pain, disability and satisfaction scores were similar between the groups at each postoperative time point, however, back pain was less at six weeks and quality of life was higher at six months in the SiCaP group than the rhBMP-2 group. Conclusions SiCaP and rhBMP-2 were comparable in terms of achieving successful bone growth and fusion. Both groups similarly alleviated pain and improved quality of life, neurological, satisfaction and return to work outcomes following PLF surgery.

A semiautomatic method to identify vertebral end plate lesions (Schmorl's nodes)

Background Context There are differences in definitions of end plate lesions (EPLs), often referred to as Schmorl’s nodes, that may, to some extent, account for the large range of reported prevalence (3.8 to 76%). Purpose To develop a technique to measure the size, prevalence and location of EPLs in a consistent manner. Study Design/Setting This study proposed a method using a detection algorithm which was applied to five adolescent females (average age 15.1 years, range 13.0 to 19.2 years) with idiopathic scoliosis (average major Cobb angle 60°, range 55 to 67°). Methods Existing low-dose, computed tomography scans were segmented semi-automatically to extract 3D morphology of each vertebral endplate. Any remaining attachments to the posterior elements of adjacent vertebrae or endplates were then manually sectioned. An automatic algorithm was used to determine the presence and position of EPLs. Results EPLs were identified in 15 of the 170 (8.8%) endplates analysed with an average depth of 3.1mm. 11/15 of the EPLs were seen in the lumbar spine. The algorithm was found to be most sensitive to changes in the minimum EPL gradient at the edges of the EPL. Conclusions This study describes an imaging analysis technique for consistent measurement of the prevalence, location and size of EPLs. The technique can be used to analyse large populations without observer errors in EPL definitions.

A virtual radiation therapy workflow training simulation

Aim Simulation forms an increasingly vital component of clinical skills development in a wide range of professional disciplines. Simulation of clinical techniques and equipment is designed to better prepare students for placement by providing an opportunity to learn technical skills in a “safe” academic environment. In radiotherapy training over the last decade or so this has predominantly comprised treatment planning software and small ancillary equipment such as mould room apparatus. Recent virtual reality developments have dramatically changed this approach. Innovative new simulation applications and file processing and interrogation software have helped to fill in the gaps to provide a streamlined virtual workflow solution. This paper outlines the innovations that have enabled this, along with an evaluation of the impact on students and educators. Method Virtual reality software and workflow applications have been developed to enable the following steps of radiation therapy to be simulated in an academic environment: CT scanning using a 3D virtual CT scanner simulation; batch CT duplication; treatment planning; 3D plan evaluation using a virtual linear accelerator; quantitative plan assessment, patient setup with lasers; and image guided radiotherapy software. Results Evaluation of the impact of the virtual reality workflow system highlighted substantial time saving for academic staff as well as positive feedback from students relating to preparation for clinical placements. Students valued practice in the “safe” environment and the opportunity to understand the clinical workflow ahead of clinical department experience. Conclusion Simulation of most of the radiation therapy workflow and tasks is feasible using a raft of virtual reality simulation applications and supporting software. Benefits of this approach include time-saving, embedding of a case-study based approach, increased student confidence, and optimal use of the clinical environment. Ongoing work seeks to determine the impact of simulation on clinical skills.

Practice patterns of radiation therapy technology in Australia: Results of a national audit

This article presents the results of a single-day census of radiation therapy (RT) treatment and technology use in Australia. The primary aim of the study was to ascertain patterns of RT practice and technology in use across Australia. These data were primarily collated to inform curriculum development of academic programs, thereby ensuring that training is matched to workforce patterns of practice. Methods: The study design was a census method with all 59 RT centres in Australia being invited to provide quantitative summary data relating to patient case mix and technology use on a randomly selected but common date. Anonymous and demographic-free data were analysed using descriptive statistics. Results: Overall data were provided across all six Australian States by 29 centres of a possible 59, yielding a response rate of 49% and representing a total of 2743 patients. Findings from this study indicate the increasing use of emerging intensity-modulated radiotherapy (IMRT), image fusion and image-guided radiation therapy (IGRT) technology in Australian RT planning and delivery phases. IMRT in particular was used for 37% of patients, indicating a high uptake of the technology in Australia when compared to other published data. The results also highlight the resource-intensive nature of benign tumour radiotherapy. Conclusions: In the absence of routine national data collection, the single-day census method offers a relatively convenient means of measuring and tracking RT resource utilisation. Wider use of this tool has the potential to not only track trends in technology implementation but also inform evidence-based guidelines for referral and resource planning.

Identifying vertebral end plate lesions (schmorl's nodes) using a semi-automatic method

INTRODUCTION There is a large range in the reported prevalence of end plate lesions (EPLs), sometimes referred to as Schmorl's nodes in the general population (3.8-76%). One possible reason for this large range is the differences in definitions used by authors. Previous research has suggested that EPLs may potentially be a primary disturbance of growth plates that leads to the onset of scoliosis. The aim of this study was to develop a technique to measure the size, prevalence and location of EPLs on Computed Tomography (CT) images of scoliosis patients in a consistent manner. METHODS A detection algorithm was developed and applied to measure EPLs for five adolescent females with idiopathic scoliosis (average age 15.1 years, average major Cobb 60°). In this algorithm, the EPL definition was based on the lesion depth, the distance from the edge of the vertebral body and the gradient of the lesion edge. Existing low-dose, CT scans of the patients' spines were segmented semi-automatically to extract 3D vertebral endplate morphology. Manual sectioning of any attachments between posterior elements of adjacent vertebrae and, if necessary, endplates was carried out before the automatic algorithm was used to determine the presence and position of EPLs. RESULTS EPLs were identified in 15 of the 170 (8.8%) endplates analysed with an average depth of 3.1mm. 73% of the EPLs were seen in the lumbar spines (11/15). A sensitivity study demonstrated that the algorithm was most sensitive to changes in the minimum gradient required at the lesion edge. CONCLUSION An imaging analysis technique for consistent measurement of the prevalence, location and size of EPLs on CT images has been developed. Although the technique was tested on scoliosis patients, it can be used to analyse other populations without observer errors in EPL definitions.

Factors influencing radiation therapy student clinical placement satisfaction

Introduction Radiation therapy students at Queensland University of Technology (QUT) attend clinical placements at five different clinical departments with varying resources and support strategies. This study aimed to determine the relative availability and perceived importance of different factors affecting student support while on clinical placement. The purpose of the research was to inform development of future support mechanisms to enhance radiation therapy students’ experience on clinical placement. Methods This study used anonymous Likert-style surveys to gather data from years 1 and 2 radiation therapy students from QUT and clinical educators from Queensland relating to availability and importance of support mechanisms during clinical placements in a semester. Results The study findings demonstrated student satisfaction with clinical support and suggested that level of support on placement influenced student employment choices. Staff support was perceived as more important than physical resources; particularly access to a named mentor, a clinical educator and weekly formative feedback. Both students and educators highlighted the impact of time pressures. Conclusions The support offered to radiation therapy students by clinical staff is more highly valued than physical resources or models of placement support. Protected time and acknowledgement of the importance of clinical education roles are both invaluable. Joint investment in mentor support by both universities and clinical departments is crucial for facilitation of effective clinical learning.

Morphometric analysis of the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis using low dose computerized tomography

Introduction. The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). Methods. In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. Results. The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). Conclusion. Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

Using low dose computerized tomography to morphometrically analyse the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis

INTRODUCTION The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). METHODS In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. RESULTS The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). CONCLUSION Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/ H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p-0.037), but the number of heterozygotes was not significantly increased. At a cost per test of £1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only £64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocatcinosis is recommended.