Faster Monte Carlo simulation of radiotherapy geometries

Using Monte Carlo simulation for radiotherapy dose calculation can provide more accurate results when compared to the analytical methods usually found in modern treatment planning systems, especially in regions with a high degree of inhomogeneity. These more accurate results acquired using Monte Carlo simulation however, often require orders of magnitude more calculation time so as to attain high precision, thereby reducing its utility within the clinical environment. This work aims to improve the utility of Monte Carlo simulation within the clinical environment by developing techniques which enable faster Monte Carlo simulation of radiotherapy geometries. This is achieved principally through the use new high performance computing environments and simpler alternative, yet equivalent representations of complex geometries. Firstly the use of cloud computing technology and it application to radiotherapy dose calculation is demonstrated. As with other super-computer like environments, the time to complete a simulation decreases as 1=n with increasing n cloud based computers performing the calculation in parallel. Unlike traditional super computer infrastructure however, there is no initial outlay of cost, only modest ongoing usage fees; the simulations described in the following are performed using this cloud computing technology. The definition of geometry within the chosen Monte Carlo simulation environment - Geometry & Tracking 4 (GEANT4) in this case - is also addressed in this work. At the simulation implementation level, a new computer aided design interface is presented for use with GEANT4 enabling direct coupling between manufactured parts and their equivalent in the simulation environment, which is of particular importance when defining linear accelerator treatment head geometry. Further, a new technique for navigating tessellated or meshed geometries is described, allowing for up to 3 orders of magnitude performance improvement with the use of tetrahedral meshes in place of complex triangular surface meshes. The technique has application in the definition of both mechanical parts in a geometry as well as patient geometry. Static patient CT datasets like those found in typical radiotherapy treatment plans are often very large and present a significant performance penalty on a Monte Carlo simulation. By extracting the regions of interest in a radiotherapy treatment plan, and representing them in a mesh based form similar to those used in computer aided design, the above mentioned optimisation techniques can be used so as to reduce the time required to navigation the patient geometry in the simulation environment. Results presented in this work show that these equivalent yet much simplified patient geometry representations enable significant performance improvements over simulations that consider raw CT datasets alone. Furthermore, this mesh based representation allows for direct manipulation of the geometry enabling motion augmentation for time dependant dose calculation for example. Finally, an experimental dosimetry technique is described which allows the validation of time dependant Monte Carlo simulation, like the ones made possible by the afore mentioned patient geometry definition. A bespoke organic plastic scintillator dose rate meter is embedded in a gel dosimeter thereby enabling simultaneous 3D dose distribution and dose rate measurement. This work demonstrates the effectiveness of applying alternative and equivalent geometry definitions to complex geometries for the purposes of Monte Carlo simulation performance improvement. Additionally, these alternative geometry definitions allow for manipulations to be performed on otherwise static and rigid geometry.

Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.

Segmental torso masses and coronal plane joint torques in the adolescent scoliotic spine

Introduction. Calculating segmental (vertebral level-by-level) torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be determined. This study used CT scans of AIS patients to measure segmental torso masses and explores how joint moments in the coronal plane are affected by changes in the position of the intervertebral joint’s axis of rotation; particularly at the apex of a scoliotic major curve. Methods. Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint torques occurring in the spine for a group of 20 female AIS patients (mean age 15.0 ± 2.7 years, mean Cobb angle 53 ± 7.1°). Image processing software, ImageJ (v1.45 NIH USA) was used to threshold the T1 to L5 CT images and calculate the segmental torso volume and mass corresponding to each vertebral level. Body segment masses for the head, neck and arms were taken from published anthropometric data. Intervertebral (IV) joint torques at each vertebral level were found using principles of static equilibrium together with the segmental body mass data. Summing the torque contributions for each level above the required joint, allowed the cumulative joint torque at a particular level to be found. Since there is some uncertainty in the position of the coronal plane Instantaneous Axis of Rotation (IAR) for scoliosis patients, it was assumed the IAR was located in the centre of the IV disc. A sensitivity analysis was performed to see what effect the IAR had on the joint torques by moving it laterally 10mm in both directions. Results. The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint torques during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm being found in patient 13. Shifting the assumed IAR by 10mm towards the convexity of the spine, increased the joint torque at that level by a mean 9.0%, showing that calculated joint torques were moderately sensitive to the assumed IAR location. When the IAR midline position was moved 10mm away from the convexity of the spine, the joint torque reduced by a mean 8.9%. Conclusion. Coronal plane joint torques as high as 7Nm can occur during relaxed standing in scoliosis patients, which may help to explain the mechanics of AIS progression. This study provides new anthropometric reference data on vertebral level-by-level torso mass in AIS patients which will be useful for biomechanical models of scoliosis progression and treatment. However, the CT scans were performed in supine (no gravitational load on spine) and curve magnitudes are known to be smaller than those measured in standing.

Measurement of a PAGAT gel dosimeter by ultrasound computed tomography

In this work we used a 3D quantitative CT ultrasound imaging system to characterise polymer gel dosimeters. The system comprised of two identical 5 MHz 128 element phased-array ultrasound transducers co-axially aligned and submerged in water as a coupling agent. Rotational and translational movement of the gel dosimeter sample between the transducers were performed using a robotic arm. Ultrasound signals were generated and received using an Olympus Omniscan unit. Dose sensitivity of attenuation and time of flight ultrasonic parameters were assessed using this system.

The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

The development of an effective vaccine against Chlamydia : utilisation of a non-toxic mucosal adjuvant to generate a protective mucosal response

Chlamydia is responsible for a wide range of diseases with enormous global economic and health burden. As the majority of chlamydial infections are asymptomatic, a vaccine has greatest potential to reduce infection and disease prevalence. Protective immunity against Chlamydia requires the induction of a mucosal immune response, ideally, at the multiple sites in the body where an infection can be established. Mucosal immunity is most effectively stimulated by targeting vaccination to the epithelium, which is best accomplished by direct vaccine application to mucosal surfaces rather than by injection. The efficacy of needle-free vaccines however is reliant on a powerful adjuvant to overcome mucosal tolerance. As very few adjuvants have proven able to elicit mucosal immunity without harmful side effects, there is a need to develop non-toxic adjuvants or safer ways to administered pre-existing toxic adjuvants. In the present study we investigated the novel non-toxic mucosal adjuvant CTA1-DD. The immunogenicity of CTA1-DD was compared to our "gold-standard" mucosal adjuvant combination of cholera toxin (CT) and cytosine-phosphate-guanosine oligodeoxynucleotide (CpG-ODN). We also utilised different needle-free immunisation routes, intranasal (IN), sublingual (SL) and transcutaneous (TC), to stimulate the induction of immunity at multiple mucosal surfaces in the body where Chlamydia are known to infect. Moreover, administering each adjuvant by different routes may also limit the toxicity of the CT/CpG adjuvant, currently restricted from use in humans. Mice were immunised with either adjuvant together with the chlamydial major outer membrane protein (MOMP) to evaluate vaccine safety and quantify the induction of antigen-specific mucosal immune responses. The level of protection against infection and disease was also assessed in vaccinated animals following a live genital or respiratory tract infectious challenge. The non-toxic CTA1-DD was found to be safe and immunogenic when delivered via the IN route in mice, inducing a comparable mucosal response and level of protective immunity against chlamydial challenge to its toxic CT/CpG counterpart administered by the same route. The utilisation of different routes of immunisation strongly influenced the distribution of antigen-specific responses to distant mucosal surfaces and also abrogated the toxicity of CT/CpG. The CT/CpG-adjuvanted vaccine was safe when administered by the SL and TC routes and conferred partial immunity against infection and pathology in both challenge models. This protection was attributed to the induction of antigen-specific pro-inflammatory cellular responses in the lymph nodes regional to the site of infection and rather than in the spleen. Development of non-toxic adjuvants and effective ways to reduce the side effects of toxic adjuvants has profound implications for vaccine development, particularly against mucosal pathogens like Chlamydia. Interestingly, we also identified two contrasting vaccines in both infection models capable of preventing infection or pathology exclusively. This indicated that the development of pathology following an infection of vaccinated animals was independent of bacterial load and was instead the result of immunopathology, potentially driven by the adaptive immune response generated following immunisation. While both vaccines expressed high levels of interleukin (IL)-17 cytokines, the pathology protected group displayed significantly reduced expression of corresponding IL-17 receptors and hence an inhibition of signalling. This indicated that the balance of IL-17-mediated responses defines the degree of protection against infection and tissue damage generated following vaccination. This study has enabled us to better understand the immune basis of pathology and protection, necessary to design more effective vaccines.

Segmental torso masses and gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Introduction: Calculating segmental (vertebral level-by-level) torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be estimated. This study used supine CT scans of AIS patients to measure segmental torso masses and explored the joint moments in the coronal plane, particularly at the apex of a scoliotic major curve. Methods: Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint moments occurring in the spine for a group of 20 female AIS patients with right sided thoracic curves. The mean age was 15.0 ± 2.7 years and all curves were classified Lenke Type 1 with a mean Cobb angle 52 ± 5.9°. Image processing software, ImageJ (v1.45 NIH USA) was used to create reformatted coronal plane images, reconstruct vertebral level-by-level torso segments and subsequently measure the torso volume corresponding to each vertebral level. Segment mass was then determined by assuming a tissue density of 1.04x103 kg/m3. Body segment masses for the head, neck and arms were taken from published anthropometric data (Winter 2009). Intervertebral joint moments in the coronal plane at each vertebral level were found from the position of the centroid of the segment masses relative to the joint centres with the segmental body mass data. Results and Discussion: The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint moments during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm. The CT scans were performed in the supine position and curve magnitudes are known to be 7-10° smaller than those measured in standing, due to the absence of gravity acting on the spine. Hence, it can be expected that the moments produced by gravity in the standing individual will be greater than those calculated here.

Analysis of the MTHFR C677T variant with migraine phenotypes

Background The methylenetetrahydrofolate reductase (MTHFR) gene variant C677T has been implicated as a genetic risk factor in migraine susceptibility, particularly in Migraine with Aura. Migraine, with and without aura (MA and MO) have many diagnostic characteristics in common. It is postulated that migraine symptomatic characteristics might themselves be influenced by MTHFR. Here we analysed the clinical profile, migraine symptoms, triggers and treatments of 267 migraineurs previously genotyped for the MTHFR C677T variant. The chi-square test was used to analyse all potential relationships between genotype and migraine clinical variables. Regression analyses were performed to assess the association of C677T with all migraine clinical variables after adjusting for gender. Findings The homozygous TT genotype was significantly associated with MA (P < 0.0001) and unilateral head pain (P = 0.002). While the CT genotype was significantly associated with physical activity discomfort (P < 0.001) and stress as a migraine trigger (P = 0.002). Females with the TT genotype were significantly associated with unilateral head pain (P < 0.001) and females with the CT genotype were significantly associated with nausea (P < 0.001), osmophobia (P = 0.002), and the use of natural remedy for migraine treatment (P = 0.003). Conversely, male migraineurs with the TT genotype experienced higher incidences of bilateral head pain (63% vs 34%) and were less likely to use a natural remedy as a migraine treatment compared to female migraineurs (5% vs 20%). Conclusions MTHFR genotype is associated with specific clinical variables of migraine including unilateral head pain, physical activity discomfort and stress.

Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated in-bred human populations. An understanding of the population history of domestic dog breeds suggests that LD mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LD mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LD was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LD in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.

Identify disorders in health records using Conditional Random Fields and Metamap: AEHRC at ShARe/CLEF 2013 eHealth Evaluation Lab Task 1

The Australian e-Health Research Centre (AEHRC) recently participated in the ShARe/CLEF eHealth Evaluation Lab Task 1. The goal of this task is to individuate mentions of disorders in free-text electronic health records and map disorders to SNOMED CT concepts in the UMLS metathesaurus. This paper details our participation to this ShARe/CLEF task. Our approaches are based on using the clinical natural language processing tool Metamap and Conditional Random Fields (CRF) to individuate mentions of disorders and then to map those to SNOMED CT concepts. Empirical results obtained on the 2013 ShARe/CLEF task highlight that our instance of Metamap (after ltering irrelevant semantic types), although achieving a high level of precision, is only able to identify a small amount of disorders (about 21% to 28%) from free-text health records. On the other hand, the addition of the CRF models allows for a much higher recall (57% to 79%) of disorders from free-text, without sensible detriment in precision. When evaluating the accuracy of the mapping of disorders to SNOMED CT concepts in the UMLS, we observe that the mapping obtained by our ltered instance of Metamap delivers state-of-the-art e ectiveness if only spans individuated by our system are considered (`relaxed' accuracy).

G-protein β3 subunit gene (GNB3) variant in causation of essential hypertension

Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein β3 subunit gene (GNB3) induces formation of a splice variant (Gβ3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (χ2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105±7, 109±16, and 128±28 mm Hg (mean±SD) for CC, CT, and 7T, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein β3 subunit gene in the causation of essential hypertension.

Fabrication and in vitro characterisation of bioactive glass composite scaffolds for bone regeneration

Here we fabricate and characterise bioactive composite scaffolds for bone tissue engineering applications. 45S5 Bioglass® (45S5) or strontium-substituted bioactive glass (SrBG) were incorporated into polycaprolactone (PCL) and fabricated into 3D bioactive composite scaffolds utilising additive manufacturing technology. We show that composite scaffolds (PCL/45S5 and PCL/SrBG) can be reproducibly manufactured with a scaffold morphology highly resembling that of PCL scaffolds. Additionally, micro-CT analysis reveals BG particles were homogeneously distributed throughout the scaffolds. Mechanical data suggested that PCL/45S5 and PCL/SrBG composite scaffolds have higher compressive Young’s modulus compared to PCL scaffolds at similar porosity (~75%). After 1 day in accelerated degradation conditions using 5M NaOH, PCL/SrBG, PCL/45S5 and PCL lost 48.6 ±3.8%, 12.1 ±1% and 1.6 ±1% of its original mass, respectively. In vitro studies were conducted using MC3T3 cells under normal and osteogenic conditions. All scaffolds were shown to be non-cytotoxic, and supported cell attachment and proliferation. Our results also indicate that the inclusion of bioactive glass (BG) promotes precipitation of calcium phosphate on the scaffold surfaces which leads to earlier cell differentiation and matrix mineralisation when compared to PCL scaffolds. However, as indicated by ALP activity, no significant difference in osteoblast differentiation was found between PCL/45S5 and PCL/SrBG scaffolds. These results suggest that PCL/45S5 and PCL/SrBG composite scaffold shows potential as a next generation bone scaffold.

Association of HincII RFLP of low density lipoprotein receptor gene with obesity in essential hypertensives

Obese (BMI ≥ 26 kg/m 2; n = 51) and lean (BMI <26 kg/m 2; n = 61) Caucasian patients with severe, familial essential hypertension, were compared with respect to genotype and allele frequencies of a HincII RFLP of the low density lipoprotein receptor gene (LDLR). A similar analysis was performed in obese (n = 28) and lean (n = 68) normotensives. A significant association of the C allele of the T→C variant responsible for this RFLP was seen with obesity (χ 2 = 4.6, P = 0.029) in the hypertensive, but not in the normotensive, group (odds ratio = 3.0 for the CC genotype and 2.7 for CT). Furthermore, BMI tracked with genotypes of this allele in the hypertensives (P = 0.046). No significant genotypic relationship was apparent for plasma lipids. Significant linkage disequilibrium was, moreover, noted between the HincII RFLP and an ApaLI RFLP (χ 2 = 33, P<0.0005) that has previously shown even stronger association with obesity (odds ratio 19.6 for cases homozygous for the susceptibility allele and 15.2 for het-erozygotes). The present study therefore adds to our previous evidence implicating LDLR as a locus for obesity in patients with essential hypertension.

Effects of hyperbaric oxygen and inducible nitric oxide synthase inhibitor treatment on femoral head osteonecrosis in a rat model

Introduction: Apoptosis is the final destiny of many cells in the body, though this process has been observed in some pathological processes. One of these pathological processes is femoral head non-traumatic osteonecrosis. Among many pro/anti-apoptotic factors, nitric oxide has recently been an area of further interest. Osteocyte apoptosis and its relation to pro-apoptotic action invite further research, and the inducible form of nitric oxide synthase (iNOS)—which produces a high concentration of nitric oxide—has been flagged. The aim of this study was to investigate the effect of hyperbaric oxygen (HBO) and inducible NOS suppressor (Aminoguanidine) in prevention of femoral head osteonecrosis in an experimental model of osteonecrosis in spontaneous hypertensive rats (SHRs). Methods: After animal ethic approval 34 SHR rats were divided into four groups. Ten rats were allocated to the control group without any treatment, and eight rats were allocated to three treatment groups namely: HBO, Aminoguanidine (AMG), and the combination of HBO and AMG treatments (HBO+AMG). The HBO group received 250 kPa of oxygen via hyperbaric chamber for 30 days started at their 5th week of life; the AMG group received 1mg/ml of AMG in drinking water from the fifth week till the 17th week of life; and the last group received a combination of these treatments. Rats were sacrificed at the end of the 17th week of life and both femurs were analysed for evidence of osteonecrosis using Micro CT scan and H&E staining. Also, osteocyte apoptosis and the presence of two different forms of NOS (inducible (iNOS) and endothelial (eNOS)) were analysed by immunostaining and apoptosis staining (Hoechst and TUNEL). Results: Bone morphology of metaphyseal and epiphyseal area of all rats were investigated and analysed. Micro CT findings revealed significantly higher mean fractional trabecular bone volume (FBV) of metaphyseal area in untreated SHRs compared with all other treatments (HBO, P<0.05, HBO+AMG, P<0.005, and AMG P<0.001). Bone surface to volume ratio also significantly increased with HBO+AMG and AMG treatments when compared with the control group (18.7 Vs 20.8, P<0.05, and 18.7 Vs 21.1, P<0.05). Epiphyseal mean FBV did not change significantly among groups. In the metaphyseal area, trabecular thickness and numbers significantly decreased with AMG treatment, while trabecular separation significantly increased with both AMG and HBO+AMG treatment. Histological ratio of no ossification and osteonecrosis was 37.5%, 43.7%, 18.7% and 6.2% of control, HBO, HBO+AMG and AMG groups respectively with only significant difference observed between HBO and AMG treatment (P<0.01). High concentration of iNOS was observed in the region of osteonecrosis while there was no evidence of eNOS activity around that region. In comparison with the control group, the ratio of osteocyte apoptosis significantly reduced in AMG treatment (P<0.005). We also observed significantly fewer apoptotic osteocytes in AMG group comparing with HBO treatment (P<0.05). Conclusion: None of our treatments prevents osteonecrosis at the histological or micro CT scan level. High concentration of iNOS in the region of osteonecrosis and significant reduction of osteocyte apoptosis with AMG treatment were supportive of iNOS modulating osteocyte apoptosis in SHRs.

Can vasculature changes be characterised morphologically after closed soft tissue trauma?

INTRODUCTION Closed soft tissue trauma (CSTT) can be the result of a blunt impact, or a prolonged crush injury and involves damage to the skin, muscles and the neurovascular system. It causes a variety of symptoms such as haematoma and in severe cases may result in hypoxia and necrosis. There is evidence that early vasculature changes following the injury delays the tissue healing [1]. However, a precise qualitative and quantitative morphological assessment of vasculature changes after trauma and the effect of this on CSTT healing is currently missing. Research aims: Developing an experimental rat model to characterise the structural changes to the vasculature after trauma qualitatively and quantitatively using micro CT. MATERIAL AND METHODS An impact device was developed to apply a controlled reproducible CSTT to the left thigh (Biceps Femoris) of anaesthetised rats [3]. After euthanizing the animals at 6 hours after trauma, CSTT was qualitatively evaluated by macroscopic observations of the skin and muscles. For vasculature visualisation, the blood vessels of sacrificed rats were flushed with heparinised saline and then perfused with a radio-opaque contrast agent (Microfil) using an infusion pump (Figure 4). The overall changes to the vasculature as a result of impact trauma were characterised qualitatively based on the 3D reconstructed images of the vasculature (Figure 5). For a smaller region of interest, the morphological parameters such as vessel thickness (diameter), spacing, and average number per volume were quantified using the scanner’s software. RESULTS AND DISCUSSION Visual observation of CSTT has revealed a haematoma in some animals (Figure 3). Micro CT images indicate good perfusion of the vasculature with contrast agent, allowing the major vessels to be identified (Figure 5). Qualitatively and quantitatively, no differences between injured and non-injured legs were observed at 6 h after trauma. Further time points of 12h, 24h, 3 days and 14 days after trauma will be characterised for identifying temporal changes of the vasculature during healing. Histomorphometical studies are required for validation of the results derived from the micro CT imaging. CONCLUSION AND FUTURE DIRECTION Findings of this research may contribute towards the establishment of a fundamental basis for the quantitative assessment and monitoring of CSTT based on microvasculature changes after trauma, which will ultimately allow for optimising the clinical treatment and improve patient outcomes.