A comparison of calculated and experimental HRTEM images of twinned boron carbide

High resolution TEM images of boron carbide (B13C2) have been recorded and compared with images calculated using the multislice method as implemented by M. A. O'Keefe in the SHRLI programs. Images calculated for the [010] zone, using machine parameters for the JEOL 2000FX AEM operating at 200 keV, indicate that for the structure model of Will et al., the optimum defocus image can be interpreted such that white spots correspond to B12 icosahedra for thin specimens and to low density channels through the structure adjacent to the direct inter-icosahedral bonds for specimens of intermediate thickness (-40 > t > -100 nm). With this information, and from the symmetry observed in the TEM images, it is likely that the (101) twin plane passes through the center of icosahedron located at the origin. This model was tested using the method of periodic continuation. Resulting images compare favorably with experimental images, thus supporting the structural model. The introduction of a (101) twin plane through the origin creates distortions to the icosahedral linkages as well as to the intra-icosahedral bonding. This increases the inequivalence of adjacent icosahedral sites along the twin plane, and thereby increases the likelihood of bipolaron hopping.

The development of a Monte Carlo system to verify radiotherapy treatment dose calculations

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 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. In this study, software has been developed that enables the transfer of treatment plan information from the treatment planning system to a Monte Carlo dose calculation engine. A database of commissioned linear accelerator models (Elekta Precise and Varian 2100CD at various energies) has been developed using the EGSnrc/BEAMnrc Monte Carlo suite. Planned beam descriptions and CT images can be exported from the treatment planning system using the DICOM framework. The information in these files is combined with an appropriate linear accelerator model to allow the accurate calculation of the radiation field incident on a modelled patient geometry. The Monte Carlo dose calculation results are combined according to the monitor units specified in the exported plan. The result is a 3D dose distribution that could be used to verify treatment planning system calculations. The software, MCDTK (Monte Carlo Dicom ToolKit), has been developed in the Java programming language and produces BEAMnrc and DOSXYZnrc input files, ready for submission on a high-performance computing cluster. The code has been tested with the Eclipse (Varian Medical Systems), Oncentra MasterPlan (Nucletron B.V.) and Pinnacle3 (Philips Medical Systems) planning systems. In this study the software was validated against measurements in homogenous and heterogeneous phantoms. Monte Carlo models are commissioned through comparison with quality assurance measurements made using a large square field incident on a homogenous volume of water. This study aims to provide a valuable confirmation that Monte Carlo calculations match experimental measurements for complex fields and heterogeneous media.

Anatomical MR imaging of long bones : comparative performance of MRI at 1.5 T and 3 T

The current gold standard for the design of orthopaedic implants is 3D models of long bones obtained using computed tomography (CT). However, high-resolution CT imaging involves high radiation exposure, which limits its use in healthy human volunteers. Magnetic resonance imaging (MRI) is an attractive alternative for the scanning of healthy human volunteers for research purposes. Current limitations of MRI include difficulties of tissue segmentation within joints and long scanning times. In this work, we explore the possibility of overcoming these limitations through the use of MRI scanners operating at a higher field strength. We quantitatively compare the quality of anatomical MR images of long bones obtained at 1.5 T and 3 T and optimise the scanning protocol of 3 T MRI. FLASH images of the right leg of five human volunteers acquired at 1.5 T and 3 T were compared in terms of signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). The comparison showed a relatively high CNR and SNR at 3 T for most regions of the femur and tibia, with the exception of the distal diaphyseal region of the femur and the mid diaphyseal region of the tibia. This was accompanied by an ~65% increase in the longitudinal spin relaxation time (T1) of the muscle at 3 T compared to 1.5 T. The results suggest that MRI at 3 T may be able to enhance the segmentability and potentially improve the accuracy of 3D anatomical models of long bones, compared to 1.5 T. We discuss how the total imaging times at 3 T can be kept short while maximising the CNR and SNR of the images obtained.

Contested notions of civics and citizenship education as national education in the Australian curriculum

Civics and Citizenship (CC) education is a contested concept and a learning area that creates curriculum and implementation challenges for schools in many nations. The current development of the first national curriculum to be implemented in Australia, the Australian Curriculum, provides a national opportunity for educators to rethink curriculum priorities and to decide on new emphases for learning in schools, in response to policy that emphasizes the importance of CC for all young Australians. In this paper, we discuss the contested notion of citizenship education as ‘national education’ in Australia, the development of this learning area and some challenges schools will encounter implementing CC in the Australian Curriculum.

Desperately seeking Asia through China: Reading ‘China’ in the Australian Curriculum: History through children’s literature

This paper considers how Asia can be meaningfully studied and understood in the first national history curriculum to be implemented in Australia. Its focus is on how empathy might be conceptualised as part of the process of becoming ‘Asia literate’ and the ways in which an empathetic understanding can be developed in the Australian Curriculum: History by engaging students with children’s literature. We argue that stories about Chinese experiences in Australia from particular episodes in the nation’s past can be utilised for their potential to prompt historical inquiry and empathetic engagement in the classroom.

An FE investigation simulating intra-operative corrective forces applied to correct scoliosis deformity

Background: Adolescent idiopathic scoliosis (AIS) is a deformity of the spine, which may 34 require surgical correction by attaching a rod to the patient’s spine using screws 35 implanted in the vertebral bodies. Surgeons achieve an intra-operative reduction in the 36 deformity by applying compressive forces across the intervertebral disc spaces while 37 they secure the rod to the vertebra. We were interested to understand how the 38 deformity correction is influenced by increasing magnitudes of surgical corrective forces 39 and what tissue level stresses are predicted at the vertebral endplates due to the 40 surgical correction. 41 Methods: Patient-specific finite element models of the osseoligamentous spine and 42 ribcage of eight AIS patients who underwent single rod anterior scoliosis surgery were 43 created using pre-operative computed tomography (CT) scans. The surgically altered 44 spine, including titanium rod and vertebral screws, was simulated. The models were 45 analysed using data for intra-operatively measured compressive forces – three load 46 profiles representing the mean and upper and lower standard deviation of this data 47 were analysed. Data for the clinically observed deformity correction (Cobb angle) were 48 compared with the model-predicted correction and the model results investigated to 49 better understand the influence of increased compressive forces on the biomechanics of 50 the instrumented joints. 51 Results: The predicted corrected Cobb angle for seven of the eight FE models were 52 within the 5° clinical Cobb measurement variability for at least one of the force profiles. 53 The largest portion of overall correction was predicted at or near the apical 54 intervertebral disc for all load profiles. Model predictions for four of the eight patients 55 showed endplate-to-endplate contact was occurring on adjacent endplates of one or 56 more intervertebral disc spaces in the instrumented curve following the surgical loading 57 steps. 58 Conclusion: This study demonstrated there is a direct relationship between intra-59 operative joint compressive forces and the degree of deformity correction achieved. The 60 majority of the deformity correction will occur at or in adjacent spinal levels to the apex 61 of the deformity. This study highlighted the importance of the intervertebral disc space 62 anatomy in governing the coronal plane deformity correction and the limit of this 63 correction will be when bone-to-bone contact of the opposing vertebral endplates 64 occurs.

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.