Monte Carlo modelling of X-ray absorptiometry and its application to body composition studies

This thesis applies Monte Carlo techniques to the study of X-ray absorptiometric methods of bone mineral measurement. These studies seek to obtain information that can be used in efforts to improve the accuracy of the bone mineral measurements. A Monte Carlo computer code for X-ray photon transport at diagnostic energies has been developed from first principles. This development was undertaken as there was no readily available code which included electron binding energy corrections for incoherent scattering and one of the objectives of the project was to study the effects of inclusion of these corrections in Monte Carlo models. The code includes the main Monte Carlo program plus utilities for dealing with input data. A number of geometrical subroutines which can be used to construct complex geometries have also been written. The accuracy of the Monte Carlo code has been evaluated against the predictions of theory and the results of experiments. The results show a high correlation with theoretical predictions. In comparisons of model results with those of direct experimental measurements, agreement to within the model and experimental variances is obtained. The code is an accurate and valid modelling tool. A study of the significance of inclusion of electron binding energy corrections for incoherent scatter in the Monte Carlo code has been made. The results show this significance to be very dependent upon the type of application. The most significant effect is a reduction of low angle scatter flux for high atomic number scatterers. To effectively apply the Monte Carlo code to the study of bone mineral density measurement by photon absorptiometry the results must be considered in the context of a theoretical framework for the extraction of energy dependent information from planar X-ray beams. Such a theoretical framework is developed and the two-dimensional nature of tissue decomposition based on attenuation measurements alone is explained. This theoretical framework forms the basis for analytical models of bone mineral measurement by dual energy X-ray photon absorptiometry techniques. Monte Carlo models of dual energy X-ray absorptiometry (DEXA) have been established. These models have been used to study the contribution of scattered radiation to the measurements. It has been demonstrated that the measurement geometry has a significant effect upon the scatter contribution to the detected signal. For the geometry of the models studied in this work the scatter has no significant effect upon the results of the measurements. The model has also been used to study a proposed technique which involves dual energy X-ray transmission measurements plus a linear measurement of the distance along the ray path. This is designated as the DPA( +) technique. The addition of the linear measurement enables the tissue decomposition to be extended to three components. Bone mineral, fat and lean soft tissue are the components considered here. The results of the model demonstrate that the measurement of bone mineral using this technique is stable over a wide range of soft tissue compositions and hence would indicate the potential to overcome a major problem of the two component DEXA technique. However, the results also show that the accuracy of the DPA( +) technique is highly dependent upon the composition of the non-mineral components of bone and has poorer precision (approximately twice the coefficient of variation) than the standard DEXA measurements. These factors may limit the usefulness of the technique. These studies illustrate the value of Monte Carlo computer modelling of quantitative X-ray measurement techniques. The Monte Carlo models of bone densitometry measurement have:- 1. demonstrated the significant effects of the measurement geometry upon the contribution of scattered radiation to the measurements, 2. demonstrated that the statistical precision of the proposed DPA( +) three tissue component technique is poorer than that of the standard DEXA two tissue component technique, 3. demonstrated that the proposed DPA(+) technique has difficulty providing accurate simultaneous measurement of body composition in terms of a three component model of fat, lean soft tissue and bone mineral,4. and provided a knowledge base for input to decisions about development (or otherwise) of a physical prototype DPA( +) imaging system. The Monte Carlo computer code, data, utilities and associated models represent a set of significant, accurate and valid modelling tools for quantitative studies of physical problems in the fields of diagnostic radiology and radiography.

The differential impact of holocaust trauma across three generations

In the current thesis, the reasons for the differential impact of Holocaust trauma on Holocaust survivors, and the differential intergenerational transmission of this trauma to survivors’ children and grandchildren were explored. A model specifically related to Holocaust trauma and its transmission was developed based on trauma, family systems and attachment theories as well as theoretical and anecdotal conjecture in the Holocaust literature. The Model of the Differential Impact of Holocaust Trauma across Three Generations was tested firstly by extensive meta-analyses of the literature pertaining to the psychological health of Holocaust survivors and their descendants and secondly via analysis of empirical study data. The meta-analyses reported in this thesis represent the first conducted with research pertaining to Holocaust survivors and grandchildren of Holocaust survivors. The meta-analysis of research conducted with children of survivors is the first to include both published and unpublished research. Meta-analytic techniques such as meta-regression and sub-set meta-analyses provided new information regarding the influence of a number of unmeasured demographic variables on the psychological health of Holocaust survivors and descendants. Based on the results of the meta-analyses it was concluded that Holocaust survivors and their children and grandchildren suffer from a statistically significantly higher level or greater severity of psychological symptoms than the general population. However it was also concluded that there is statistically significant variation in psychological health within the Holocaust survivor and descendant populations. Demographic variables which may explain a substantial amount of this variation have been largely under-assessed in the literature and so an empirical study was needed to clarify the role of demographics in determining survivor and descendant mental health. A total of 124 participants took part in the empirical study conducted for this thesis with 27 Holocaust survivors, 69 children of survivors and 28 grandchildren of survivors. A worldwide recruitment process was used to obtain these participants. Among the demographic variables assessed in the empirical study, aspects of the survivors’ Holocaust trauma (namely the exact nature of their Holocaust experiences, the extent of family bereavement and their country of origin) were found to be particularly potent predictors of not only their own psychological health but continue to be strongly influential in determining the psychological health of their descendants. Further highlighting the continuing influence of the Holocaust was the finding that number of Holocaust affected ancestors was the strongest demographic predictor of grandchild of survivor psychological health. Apart from demographic variables, the current thesis considered family environment dimensions which have been hypothesised to play a role in the transmission of the traumatic impact of the Holocaust from survivors to their descendants. Within the empirical study, parent-child attachment was found to be a key determinant in the transmission of Holocaust trauma from survivors to their children and insecure parent-child attachment continues to reverberate through the generations. In addition, survivors’ communication about the Holocaust and their Holocaust experiences to their children was found to be more influential than general communication within the family. Ten case studies (derived from the empirical study data set) are also provided; five Holocaust survivors, three children of survivors and two grandchildren of survivors. These cases add further to the picture of heterogeneity of the survivor and descendant populations in both experiences and adaptations. It is concluded that the legacy of the Holocaust continues to leave its mark on both its direct survivors and their descendants. Even two generations removed, the direct and indirect effects of the Holocaust have yet to be completely nullified. Research with Holocaust survivor families serves to highlight the differential impacts of state-based trauma and the ways in which its effects continue to be felt for generations. The revised and empirically tested Model of the Differential Impact of Holocaust Trauma across Three Generations presented at the conclusion of this thesis represents a further clarification of existing trauma theories as well as the first attempt at determining the relative importance of both cognitive, interpersonal/interfamilial interaction processes and demographic variables in post-trauma psychological health and transmission of traumatic impact.

Repair of calvarial defects with customized tissue-engineered bone grafts - I. Evaluation of osteogenesis in a three-dimensional culture system

Bone generation by autogenous cell transplantation in combination with a biodegradable scaffold is one of the most promising techniques being developed in craniofacial surgery. The objective of this combined in vitro and in vivo study was to evaluate the morphology and osteogenic differentiation of bone marrow derived mesenchymal progenitor cells and calvarial osteoblasts in a two-dimensional (2-D) and three-dimensional (3-D) culture environment (Part I of this study) and their potential in combination with a biodegradable scaffold to reconstruct critical-size calvarial defects in an autologous animal model [Part II of this study; see Schantz, J.T., et al. Tissue Eng. 2003;9(Suppl. 1):S-127-S-139; this issue]. New Zealand White rabbits were used to isolate osteoblasts from calvarial bone chips and bone marrow stromal cells from iliac crest bone marrow aspirates. Multilineage differentiation potential was evaluated in a 2-D culture setting. After amplification, the cells were seeded within a fibrin matrix into a 3-D polycaprolactone (PCL) scaffold system. The constructs were cultured for up to 3 weeks in vitro and assayed for cell attachment and proliferation using phase-contrast light, confocal laser, and scanning electron microscopy and the MTS cell metabolic assay. Osteogenic differentiation was analyzed by determining the expression of alkaline phosphatase (ALP) and osteocalcin. The bone marrow-derived progenitor cells demonstrated the potential to be induced to the osteogenic, adipogenic, and chondrogenic pathways. In a 3-D environment, cell-seeded PCL scaffolds evaluated by confocal laser microscopy revealed continuous cell proliferation and homogeneous cell distribution within the PCL scaffolds. On osteogenic induction mesenchymal progenitor cells (12 U/L) produce significantly higher (p < 0.05) ALP activity than do osteoblasts (2 U/L); however, no significant differences were found in osteocalcin expression. In conclusion, this study showed that the combination of a mechanically stable synthetic framework (PCL scaffolds) and a biomimetic hydrogel (fibrin glue) provides a potential matrix for bone tissue-engineering applications. Comparison of osteogenic differentiation between the two mesenchymal cell sources revealed a similar pattern.

Scaffold development using 3D printing with a starch-based polymer

Rapid prototyping (RP) techniques have been utilised by tissue engineers to produce three-dimensional (3D) porous scaffolds. RP technologies allow the design and fabrication of complex scaffold geometries with a fully interconnected pore network. Three-dimensional printing (3DP) technique was used to fabricate scaffolds with a novel micro- and macro-architecture. In this study, a unique blend of starch-based polymer powders (cornstarch, dextran and gelatin) was developed for the 3DP process. Cylindrical scaffolds of five different designs were fabricated and post-processed to enhance the mechanical and chemical properties. The scaffold properties were characterised by scanning electron microscopy (SEM), differential scanning calorimetry (DSC), porosity analysis and compression tests

A comparative study of the effect of calibration conditions on the water equivalence of a range of gel dosimeters

Gel dosimeters are of increasing interest in the field of radiation oncology as the only truly three-dimensional integrating radiation dosimeter. There are a range of ferrous-sulphate and polymer gel dosimeters. To be of use, they must be water-equivalent. On their own, this relates to their radiological properties as determined by their composition. In the context of calibration of gel dosimeters, there is the added complexity of the calibration geometry; the presence of containment vessels may influence the dose absorbed. Five such methods of calibration are modelled here using the Monte Carlo method. It is found that the Fricke gel best matches water for most of the calibration methods, and that the best calibration method involves the use of a large tub into which multiple fields of different dose are directed. The least accurate calibration method involves the use of a long test tube along which a depth dose curve yields multiple calibration points.

A comparative study of the effect of calibration conditions on the water equivalence of a range of gel dosimeters

Gel dosimeters are of increasing interest in the field of radiation oncology as the only truly three-dimensional integrating radiation dosimeter. There are a range of ferrous-sulphate and polymer gel dosimeters. To be of use, they must be water-equivalent. On their own, this relates to their radiological properties as determined by their composition. In the context of calibration of gel dosimeters, there is the added complexity of the calibration geometry; the presence of containment vessels may influence the dose absorbed. Five such methods of calibration are modelled here using the Monte Carlo method. It is found that the Fricke gel best matches water for most of the calibration methods, and that the best calibration method involves the use of a large tub into which multiple fields of different dose are directed. The least accurate calibration method involves the use of a long test tube along which a depth dose curve yields multiple calibration points.

The feasibility of vibration analysis as a technique to detect osseointegration of transfemoral implants

One of the main causes of above knee or transfemoral amputation (TFA) in the developed world is trauma to the limb. The number of people undergoing TFA due to limb trauma, particularly due to war injuries, has been increasing. Typically the trauma amputee population, including war-related amputees, are otherwise healthy, active and desire to return to employment and their usual lifestyle. Consequently there is a growing need to restore long-term mobility and limb function to this population. Traditionally transfemoral amputees are provided with an artificial or prosthetic leg that consists of a fabricated socket, knee joint mechanism and a prosthetic foot. Amputees have reported several problems related to the socket of their prosthetic limb. These include pain in the residual limb, poor socket fit, discomfort and poor mobility. Removing the socket from the prosthetic limb could eliminate or reduce these problems. A solution to this is the direct attachment of the prosthesis to the residual bone (femur) inside the residual limb. This technique has been used on a small population of transfemoral amputees since 1990. A threaded titanium implant is screwed in to the shaft of the femur and a second component connects between the implant and the prosthesis. A period of time is required to allow the implant to become fully attached to the bone, called osseointegration (OI), and be able to withstand applied load; then the prosthesis can be attached. The advantages of transfemoral osseointegration (TFOI) over conventional prosthetic sockets include better hip mobility, sitting comfort and prosthetic retention and fewer skin problems on the residual limb. However, due to the length of time required for OI to progress and to complete the rehabilitation exercises, it can take up to twelve months after implant insertion for an amputee to be able to load bear and to walk unaided. The long rehabilitation time is a significant disadvantage of TFOI and may be impeding the wider adoption of the technique. There is a need for a non-invasive method of assessing the degree of osseointegration between the bone and the implant. If such a method was capable of determining the progression of TFOI and assessing when the implant was able to withstand physiological load it could reduce the overall rehabilitation time. Vibration analysis has been suggested as a potential technique: it is a non destructive method of assessing the dynamic properties of a structure. Changes in the physical properties of a structure can be identified from changes in its dynamic properties. Consequently vibration analysis, both experimental and computational, has been used to assess bone fracture healing, prosthetic hip loosening and dental implant OI with varying degrees of success. More recently experimental vibration analysis has been used in TFOI. However further work is needed to assess the potential of the technique and fully characterise the femur-implant system. The overall aim of this study was to develop physical and computational models of the TFOI femur-implant system and use these models to investigate the feasibility of vibration analysis to detect the process of OI. Femur-implant physical models were developed and manufactured using synthetic materials to represent four key stages of OI development (identified from a physiological model), simulated using different interface conditions between the implant and femur. Experimental vibration analysis (modal analysis) was then conducted using the physical models. The femur-implant models, representing stage one to stage four of OI development, were excited and the modal parameters obtained over the range 0-5kHz. The results indicated the technique had limited capability in distinguishing between different interface conditions. The fundamental bending mode did not alter with interfacial changes. However higher modes were able to track chronological changes in interface condition by the change in natural frequency, although no one modal parameter could uniquely distinguish between each interface condition. The importance of the model boundary condition (how the model is constrained) was the key finding; variations in the boundary condition altered the modal parameters obtained. Therefore the boundary conditions need to be held constant between tests in order for the detected modal parameter changes to be attributed to interface condition changes. A three dimensional Finite Element (FE) model of the femur-implant model was then developed and used to explore the sensitivity of the modal parameters to more subtle interfacial and boundary condition changes. The FE model was created using the synthetic femur geometry and an approximation of the implant geometry. The natural frequencies of the FE model were found to match the experimental frequencies within 20% and the FE and experimental mode shapes were similar. Therefore the FE model was shown to successfully capture the dynamic response of the physical system. As was found with the experimental modal analysis, the fundamental bending mode of the FE model did not alter due to changes in interface elastic modulus. Axial and torsional modes were identified by the FE model that were not detected experimentally; the torsional mode exhibited the largest frequency change due to interfacial changes (103% between the lower and upper limits of the interface modulus range). Therefore the FE model provided additional information on the dynamic response of the system and was complementary to the experimental model. The small changes in natural frequency over a large range of interface region elastic moduli indicated the method may only be able to distinguish between early and late OI progression. The boundary conditions applied to the FE model influenced the modal parameters to a far greater extent than the interface condition variations. Therefore the FE model, as well as the experimental modal analysis, indicated that the boundary conditions need to be held constant between tests in order for the detected changes in modal parameters to be attributed to interface condition changes alone. The results of this study suggest that in a clinical setting it is unlikely that the in vivo boundary conditions of the amputated femur could be adequately controlled or replicated over time and consequently it is unlikely that any longitudinal change in frequency detected by the modal analysis technique could be attributed exclusively to changes at the femur-implant interface. Therefore further development of the modal analysis technique would require significant consideration of the clinical boundary conditions and investigation of modes other than the bending modes.

Three-dimensional hydrogen-bonded structures in the guanidinium salts of the monocyclic dicarboxylic acids rac-trans-cyclohexane-1,2-dicarboxylic acid (2:1, anhydrous), isophthalic acid (1:1, monohydrate) and terephthalic acid (2:1, trihydrate).

The structures of bis(guanidinium)rac-trans-cyclohexane-1,2-dicarboxylate, 2(CH6N3+) C8H10O4- (I), guanidinium 3-carboxybenzoate monohydrate CH6N3+ C8H5O4- . H2O (II) and bis(guanidinium) benzene-1,4-dicarboxylate trihydrate, 2(CH6N3+) C8H4O4^2- . 3H2O (III) have been determined and the hydrogen bonding in each examined. All three compounds form three-dimensional hydrogen-bonded framework structures. In anhydrous (I), both guanidinium cations give classic cyclic R2/2(8) N--H...O,O'(carboxyl) and asymmetric cyclic R1/2(6) hydrogen-bonding interactions while one cation gives an unusual enlarged cyclic interaction with O acceptors of separate ortho-related carboxyl groups [graph set R2/2(11)]. Cations and anions also associate across inversion centres giving cyclic R2/4(8) motifs. In the 1:1 guanidinium salt (II), the cation gives two separate cyclic R1/2(6) interactions, one with a carboxyl O-acceptor, the other with the water molecule of solvation. The structure is unusual in that both carboxyl groups give short inter-anion O...H...O contacts, one across a crystallographic inversion centre [2.483(2)\%A], the other about a two-fold axis of rotation [2.462(2)\%A] with a half-occupancy hydrogen delocalized on the symmetry element in each. The water molecule links the cation--anion ribbon structures into a three-dimensional framework. In (III), the repeating molecular unit comprises a benzene-1,4-dicarboxylate dianion which lies across a crystallographic inversion centre, two guanidinium cations and two water molecules of solvation (each set related by two-fold rotational symmetry), and a single water molecule which lies on a two-fold axis. Each guanidinium cation gives three types of cyclic interactions with the dianions: one R^1^~2~(6), the others R2/3(8) and R3/3(10) (both of these involving the water molecules), giving a three-dimensional structure through bridges down the b cell direction. The water molecule at the general site also forms an unusual cyclic R2/2(4) homodimeric association across an inversion centre [O--H...O, 2.875(2)\%A]. The work described here provides further examples of the common cyclic guanidinium cation...carboxylate anion hydrogen-bonding associations as well as featuring other less common cyclic motifs.

Differential roles of the pRb and Arf/p53 pathways in murine naevus and melanoma genesis

We report on a systematic analysis of genotype-specific melanocyte (MC) UVR responses in transgenic mouse melanoma models along with tumour penetrance and comparative histopathology. pRb or p53 pathway mutations cooperated with NrasQ61K to transform MCs. We previously reported that MCs migrate from the follicular outer root sheath into the epidermis after neonatal UVR. Here, we found that Arf or p53 loss markedly diminished this response. Despite this, mice carrying these mutations developed melanoma with very early age of onset after neonatal UVR. Cdk4R24C did not affect the MC migration. Instead, independent of UVR exposure, interfollicular dermal MCs were more prevalent in Cdk4R24C mice. Subsequently, in adulthood, these mutants developed dermal MC proliferations reminiscent of superficial congenital naevi. Two types of melanoma were observed in this model. The location and growth pattern of the first was consistent with derivation from the naevi, while the second appeared to be of deep dermal origin. In animals carrying the Arf or p53 defects, no naevi were detected, with all tumours ostensibly skipping the benign precursor stage in progression.