Bayesian computational methods for spatial analysis of images

This thesis introduces a new way of using prior information in a spatial model and develops scalable algorithms for fitting this model to large imaging datasets. These methods are employed for image-guided radiation therapy and satellite based classification of land use and water quality. This study has utilized a pre-computation step to achieve a hundredfold improvement in the elapsed runtime for model fitting. This makes it much more feasible to apply these models to real-world problems, and enables full Bayesian inference for images with a million or more pixels.

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.

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.

Ectopic calcification among families in the Azores: Clinical and radiologic manifestations in families with diffuse idiopathic skeletal hyperostosis and chondrocalcinosis

Objective. Twelve families that were multiply affected with diffuse idiopathic skeletal hyperostosis (DISH) and/or chondrocalcinosis, were identified on the island of Terceira, The Azores, potentially supporting the hypothesis that the 2 disorders share common etiopathogenic factors. The present study was undertaken to investigate this hypothesis. Methods. One hundred three individuals from 12 unrelated families were assessed. Probands were identified from patients attending the Rheumatic Diseases Clinic, Hospital de Santo Espirito, in The Azores. Family members were assessed by rheumatologists and radiologists. Radiographs of all family members were obtained, including radiographs of the dorsolumbar spine, pelvis, knees, elbows, and wrists, and all cases were screened for known features of chondrocalcinosis. Results. Ectopic calcifications were identified in 70 patients. The most frequent symptoms or findings were as follows: axial pain, elbow, knee and metacarpophalangeal (MCP) joint pain, swelling, and/or deformity, and radiographic enthesopathic changes. Elbow and MCP joint periarticular calcifications were observed in 35 and 5 patients, respectively, and chondrocalcinosis was identified in 12 patients. Fifteen patients had sacroiliac disease (ankylosis or sclerosis) on computed tomography scans. Fifty-two patients could be classified as having definite (17%), probable (26%), or possible (31%) DISH. Concomitant DISH and chondrocalcinosis was diagnosed in 12 patients. Pyrophosphate crystals were identified from knee effusions in 13 patients. The pattern of disease transmission was compatible with an autosomal-dominant monogenic disease. The mean age at which symptoms developed was 38 years. Conclusion. These families may represent a familial type of pyrophosphate arthropathy with a phenotype that includes peripheral and axial enthesopathic calcifications. The concurrence of DISH and chondrocalcinosis suggests a shared pathogenic mechanism in the 2 conditions.

Biofabrication of customized bone grafts by combination of additive manufacturing and bioreactor knowhow

This study reports on an original concept of additive manufacturing for the fabrication of tissue engineered constructs (TEC), offering the possibility of concomitantly manufacturing a customized scaffold and a bioreactor chamber to any size and shape. As a proof of concept towards the development of anatomically relevant TECs, this concept was utilized for the design and fabrication of a highly porous sheep tibia scaffold around which a bioreactor chamber of similar shape was simultaneously built. The morphology of the bioreactor/scaffold device was investigated by micro-computed tomography and scanning electron microscopy confirming the porous architecture of the sheep tibiae as opposed to the non-porous nature of the bioreactor chamber. Additionally, this study demonstrates that both the shape, as well as the inner architecture of the device can significantly impact the perfusion of fluid within the scaffold architecture. Indeed, fluid flow modelling revealed that this was of significant importance for controlling the nutrition flow pattern within the scaffold and the bioreactor chamber, avoiding the formation of stagnant flow regions detrimental for in vitro tissue development. The bioreactor/scaffold device was dynamically seeded with human primary osteoblasts and cultured under bi-directional perfusion for two and six weeks. Primary human osteoblasts were observed homogenously distributed throughout the scaffold, and were viable for the six week culture period. This work demonstrates a novel application for additive manufacturing in the development of scaffolds and bioreactors. Given the intrinsic flexibility of the additive manufacturing technology platform developed, more complex culture systems can be fabricated which would contribute to the advances in customized and patient-specific tissue engineering strategies for a wide range of applications.

A phase 1b study of placenta-derived mesenchymal stromal cells in patients with idiopathic pulmonary fibrosis

BACKGROUND AND OBJECTIVE: Idiopathic pulmonary fibrosis (IPF) is a degenerative disease characterized by fibrosis following failed epithelial repair. Mesenchymal stromal cells (MSC), a key component of the stem cell niche in bone marrow and possibly other organs including lung, have been shown to enhance epithelial repair and are effective in preclinical models of inflammation-induced pulmonary fibrosis, but may be profibrotic in some circumstances. METHODS: In this single centre, non-randomized, dose escalation phase 1b trial, patients with moderately severe IPF (diffusing capacity for carbon monoxide (DLCO ) ≥ 25% and forced vital capacity (FVC) ≥ 50%) received either 1 × 10(6) (n = 4) or 2 × 10(6) (n = 4) unrelated-donor, placenta-derived MSC/kg via a peripheral vein and were followed for 6 months with lung function (FVC and DLCO ), 6-min walk distance (6MWD) and computed tomography (CT) chest. RESULTS: Eight patients (4 female, aged 63.5 (57-75) years) with median (interquartile range) FVC 60 (52.5-74.5)% and DLCO 34.5 (29.5-40)% predicted were treated. Both dose schedules were well tolerated with only minor and transient acute adverse effects. MSC infusion was associated with a transient (1% (0-2%)) fall in SaO2 after 15 min, but no changes in haemodynamics. At 6 months FVC, DLCO , 6MWD and CT fibrosis score were unchanged compared with baseline. There was no evidence of worsening fibrosis. CONCLUSIONS: Intravenous MSC administration is feasible and has a good short-term safety profile in patients with moderately severe IPF.

Synthesis of a multimodal molecular imaging probe based on a hyperbranched polymer architecture

Molecular imaging is utilised in modern medicine to aid in the diagnosis and treatment of disease by allowing its spatiotemporal state to be examined in vivo. This study focuses on the development of novel multimodal molecular imaging agents based on hyperbranched polymers that combine the complementary capabilities of optical fluorescence imaging and positron emission tomography-computed tomography (PET/CT) into one construct. RAFT-mediated polymerisation was used to prepare two hydrophilic hyperbranched polymers that were differentiated by their size and level of branching. The multiple functional end-groups facilitated covalent attachment of both near infrared fluorescent dyes for optical imaging, as well as a copper chelator allowing binding of 64Cu as a PET radio nuclei. In vivo multimodal imaging of mice using PET/CT and planar optical imaging was first used to assess the biodistribution of the polymeric materials and it was shown that the larger and more branched polymer had a significantly longer circulation time. The larger constructs were also shown to exhibit enhanced accumulation in solid tumours in a murine B16 melanoma model. Importantly, it was demonstrated that the PET modality gave rise to high sensitivity immediately after injection of the agent, while the optical modality facilitated extended longitudinal studies, thus highlighting how the complementary capabilities of the molecular imaging agents can be useful for studying various diseases, including cancer.

A mouse model for spondyloepiphyseal dysplasia congenita with secondary osteoarthritis due to a Col2a1 mutation

Progeny of mice treated with the mutagen N-ethyl-N-nitrosourea (ENU) revealed a mouse, designated Longpockets (Lpk), with short humeri, abnormal vertebrae, and disorganized growth plates, features consistent with spondyloepiphyseal dysplasia congenita (SEDC). The Lpk phenotype was inherited as an autosomal dominant trait. Lpk/+ mice were viable and fertile and Lpk/Lpk mice died perinatally. Lpk was mapped to chromosome 15 and mutational analysis of likely candidates from the interval revealed a Col2a1 missense Ser1386Pro mutation. Transient transfection of wild-type and Ser1386Pro mutant Col2a1 c-Myc constructs in COS-7 cells and CH8 chondrocytes demonstrated abnormal processing and endoplasmic reticulum retention of the mutant protein. Histology revealed growth plate disorganization in 14-day-old Lpk/+ mice and embryonic cartilage from Lpk/+ and Lpk/Lpk mice had reduced safranin-O and type-II collagen staining in the extracellular matrix. The wild-type and Lpk/+ embryos had vertical columns of proliferating chondrocytes, whereas those in Lpk/Lpk mice were perpendicular to the direction of bone growth. Electron microscopy of cartilage from 18.5 dpc wild-type, Lpk/+, and Lpk/Lpk embryos revealed fewer and less elaborate collagen fibrils in the mutants, with enlarged vacuoles in the endoplasmic reticulum that contained amorphous inclusions. Micro-computed tomography (CT) scans of 12-week-old Lpk/+ mice revealed them to have decreased bone mineral density, and total bone volume, with erosions and osteophytes at the joints. Thus, an ENU mouse model with a Ser1386Pro mutation of the Col2a1 C-propeptide domain that results in abnormal collagen processing and phenotypic features consistent with SEDC and secondary osteoarthritis has been established.

Association between aneurysm shoulder stress and abdominal aortic aneurysm expansion: A longitudinal follow-up study

Background Aneurysm expansion rate is an important indicator of the potential risk of abdominal aortic aneurysm (AAA) rupture. Stress within the AAA wall is also thought to be a trigger for its rupture. However, the association between aneurysm wall stresses and expansion of AAA is unclear. Methods and Results Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computed tomography scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up computed tomography images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A nonlinear large-strain finite element method was used to compute the wall-stress distribution. The relationship between wall stresses and expansion rate was investigated. Slowly and rapidly expanding aneurysms had comparable baseline maximum diameters (median, 4.35 cm [interquartile range, 4.12 to 5.0 cm] versus 4.6 cm [interquartile range, 4.2 to 5.0 cm]; P=0.32). Rapidly expanding AAAs had significantly higher shoulder stresses than slowly expanding AAAs (median, 300 kPa [interquartile range, 280 to 320 kPa] versus 225 kPa [interquartile range, 211 to 249 kPa]; P=0.0001). A good correlation between shoulder stress at baseline and expansion rate was found (r=0.71; P=0.0001). Conclusion A higher shoulder stress was found to have an association with a rapidly expanding AAA. Therefore, it may be useful for estimating the expansion of AAAs and improve risk stratification of patients with AAAs.

Invasive and non-invasive modalities of imaging carotid stenosis

Despite recent therapeutic advances, acute ischemic complications of atherosclerosis remain the primary cause of morbidity and mortality in Western countries, with carotid atherosclerotic disease one of the major preventable causes of stroke. As the impact of this disease challenges our healthcare systems, we are becoming aware that factors influencing this disease are more complex than previously realized. In current clinical practice, risk stratification relies primarily on evaluation of the degree of luminal stenosis and patient symptomatology. Adequate investigation and optimal imaging are important factors that affect the quality of a carotid endarterectomy (CEA) service and are fundamental to patient selection. Digital subtraction angiography is still perceived as the most accurate imaging modality for carotid stenosis and historically has been the cornerstone of most of the major CEA trials but concerns regarding potential neurological complications have generated substantial interest in non-invasive modalities, such as contrast-enhanced magnetic resonance angiography. The purpose of this review is to give an overview to the vascular specialist of the current imaging modalities in clinical practice to identify patients with carotid stenosis. Advantages and disadvantages of each technique are outlined. Finally, limitations of assessing luminal stenosis in general are discussed. This article will not cover imaging of carotid atheroma morphology, function and other emerging imaging modalities of assessing plaque risk, which look beyond simple luminal measurements.

Comparison of active-set method deconvolution and matched-filtering for derivation of an ultrasound transit time spectrum

The quality of ultrasound computed tomography imaging is primarily determined by the accuracy of ultrasound transit time measurement. A major problem in analysis is the overlap of signals making it difficult to detect the correct transit time. The current standard is to apply a matched-filtering approach to the input and output signals. This study compares the matched-filtering technique with active set deconvolution to derive a transit time spectrum from a coded excitation chirp signal and the measured output signal. The ultrasound wave travels in a direct and a reflected path to the receiver, resulting in an overlap in the recorded output signal. The matched-filtering and deconvolution techniques were applied to determine the transit times associated with the two signal paths. Both techniques were able to detect the two different transit times; while matched-filtering has a better accuracy (0.13 μs vs. 0.18 μs standard deviation), deconvolution has a 3.5 times improved side-lobe to main-lobe ratio. A higher side-lobe suppression is important to further improve image fidelity. These results suggest that a future combination of both techniques would provide improved signal detection and hence improved image fidelity.

CBCT image guidance in head and neck irradiation: The impact of daily and weekly imaging protocols

Purpose This study evaluated the impact of a daily and weekly image-guided radiotherapy protocols in reducing setup errors and setting of appropriate margins in head and neck cancer patients. Materials and methods Interfraction and systematic shifts for the hypothetical day 1–3 plus weekly imaging were extrapolated from daily imaging data from 31 patients (964 cone beam computed tomography (CBCT) scans). In addition, residual setup errors were calculated by taking the average shifts in each direction for each patient based on the first three shifts and were presumed to represent systematic setup error. The clinical target volume (CTV) to planning target volume (PTV) margins were calculated using van Herk formula and analysed for each protocol. Results The mean interfraction shifts for daily imaging were 0·8, 0·3 and 0·5 mm in the S-I (superior-inferior), L-R (left-right) and A-P (anterior-posterior) direction, respectively. On the other hand the mean shifts for day 1–3 plus weekly imaging were 0·9, 1·8 and 0·5 mm in the S-I, L-R and A-P direction, respectively. The mean day 1–3 residual shifts were 1·5, 2·1 and 0·7 mm in the S-I, L-R and A-P direction, respectively. No significant difference was found in the mean setup error for the daily and hypothetical day 1–3 plus weekly protocol. However, the calculated CTV to PTV margins for the daily interfraction imaging data were 1·6, 3·8 and 1·4 mm in the S-I, L-R and A-P directions, respectively. Hypothetical day 1–3 plus weekly resulted in CTV–PTV margins of 5, 4·2 and 5 mm in the S-I, L-R and A-P direction. Conclusions The results of this study show that a daily CBCT protocol reduces setup errors and allows setup margin reduction in head and neck radiotherapy compared to a weekly imaging protocol.

Micro-CT based characterisation of changes to the vascular network following closed soft tissue trauma and cryotherapy

This project has investigated how the architecture of the blood vessels supplying nutrients to skeletal muscles is affected by muscle contusion injuries, and how it changes during healing with or without initial treatment of the injury by icing. In order to do this, we used contrast agents to visualise blood vessels in 3D with micro-computed tomography imaging. This research significantly contributes to the fields of orthopaedics, traumatology and sports medicine, as it improves our understanding of muscle contusion injuries. Furthermore, the methods developed in this thesis may help to improve the diagnosis and monitoring of these injuries.

Ultrasound of paediatric appendicitis and its secondary sonographic signs: Providing a more meaningful finding

Sonography is an important clinical tool in diagnosing appendicitis in children as it can obviate both exposure to potentially harmful ionising radiation from computed tomography scans and the need for unnecessary appendicectomies. This review examines the diagnostic accuracy of ultrasound in the identification of acute appendicitis, with a particular focus on the the utility of secondary sonographic signs as an adjunct or corollary to traditionally examined criteria. These secondary signs can be important in cases where the appendix cannot be identified with ultrasound and a more meaningful finding may be made by incorporating the presence or absence of secondary sonographic signs. There is evidence that integrating these secondary signs into the final ultrasound diagnosis can improve the utility of ultrasound in cases where appendicitis is expected, though there remains some conjecture about whether they play a more important role in negative or positive prediction in the absence of an identifiable appendix.