Sequential MRI reveals individual level deformities in the growing scoliotic spine that are clinically masked by the Cobb angle

Introduction Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. Methods A 0.49mm isotropic 3D MRI technique was developed to investigate the level-by-level changes that occur in the growing spine of a group of Adolescent Idiopathic Scoliosis (AIS) patients, who received two to four sequential scans (spaced 3-12 months apart). The coronal plane wedge angles of each vertebra and disc in the major curve were measured to capture any changes that occurred during their adolescent growth phase. Results Seventeen patients had at least two scans. Mean patient age was 12.9 years (SD 1.5 years). Sixteen were classified as right-sided major thoracic Lenke Type 1 (one left sided). Mean standing Cobb angle at initial presentation was 31° (SD 12°). Six received two scans, nine three scans and two four scans, with 65% showing a Cobb angle progression of 5° or more between scans. Overall, there was no clear pattern of deformity progression of individual vertebrae and discs, nor between patients who progressed and those who didn’t. There were measurable changes in the wedging of the vertebrae and discs in all patients. In sequential scans, change in direction of wedging was also seen. In several patients there was reverse wedging in the discs that counteracted increased wedging of the vertebrae such that no change in overall Cobb angle was seen. Conclusion Sequential MRI data showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

The vertebral venous system in healthy and scoliotic adolescent spines - a 3D MRI investigation

Introduction. The venous drainage system within vertebral bodies (VBs) has been well documented previously in cadaveric specimens. Advances in 3D imaging and image processing now allow for in vivo quantification of larger venous vessels, such as the basivertebral vein. Differences between healthy and scoliotic VB veins can therefore be investigated. Methods. 20 healthy adolescent controls and 21 AIS patients were recruited (with ethics approval) to undergo 3D MRI, using a 3 Tesla, T1-weighted 3D gradient echo sequence, resulting in 512 slices across the thoraco-lumbar spine, with a voxel size of 0.5x0.5x0.5mm. Using Amira Filament Editor, five transverse slices through the VB were examined simultaneously and the resulting observable vascular network traced. Each VB was assessed, and a vascular network recorded when observable. A local coordinate system was created in the centre of each VB and the vascular networks aligned to this. The length of the vascular network on the left and right sides (with a small central region) of the VB was calculated, and the spatial patterning of the networks assessed level-by-level within each subject. Results. An average of 6 (range 4-10) vascular networks, consistent with descriptions of the basivertebral vein, were identifiable within each subject, most commonly between T10-L1. Differences were seen in the left/right distribution of vessels in the control and AIS subjects. Healthy controls saw a percentage distribution of 29:18:53 across the left:centre:right regions respectively, whereas the AIS subjects had a slightly shifted distribution of 33:25:42. The control group showed consistent spatial patterning of the vascular networks across most levels, but this was not seen in the AIS group. Conclusion. Observation and quantification of the basivertebral vein in vivo is possible using 3D MRI. The AIS group lacked the spatial pattern repetition seen in the control group and minor differences were seen in the left/right distribution of vessels.

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

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

Sequential MRI of growing scoliotic spines reveals individual level deformities that are clinically masked by the cobb angle

INTRODUCTION. Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. METHODS. A 0.49mm isotropic 3D MRI technique was developed to investigate the level-by-level changes that occur in the growing spine of a group of Adolescent Idiopathic Scoliosis (AIS) patients, who received two to four sequential scans (spaced 3-12 months apart). The coronal plane wedge angles of each vertebra and disc in the major curve were measured to capture any changes that occurred during their adolescent growth phase. RESULTS. Seventeen patients had at least two scans. Mean patient age was 12.9 years (SD 1.5 years). Sixteen were classified as right-sided major thoracic Lenke Type 1 (one left sided). Mean standing Cobb angle at initial presentation was 31° (SD 12°). Six received two scans, nine three scans and two four scans, with 65% showing a Cobb angle progression of 5° or more between scans. Overall, there was no clear pattern of deformity progression of individual vertebrae and discs, nor between patients who progressed and those who didn’t. There were measurable changes in the wedging of the vertebrae and discs in all patients. In sequential scans, change in direction of wedging was also seen. In several patients there was reverse wedging in the discs that counteracted increased wedging of the vertebrae such that no change in overall Cobb angle was seen. CONCLUSION. Sequential MRI data showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

On the biomechanical significance of inter-lamellar interfaces in the intervertebral disc

The intervertebral disc (IVD) is a unique soft tissue structure which provides structural support and flexibility in the axial skeleton of vertebrates. From a structural perspective, the disc behaves somewhat like a thick walled pressure vessel, where the walls are comprised of a series of composite annular rings (lamellae). However, a prior study (Marchand and Ahmed, 1990) found a high proportion of circumferentially discontinuous lamellae in human lumbar IVDs. The presence of these discontinuities raises important structural questions, because discontinuous lamellae cannot withstand high nucleus pressures via the generation of circumferential (hoop) stress. A possible alternative mechanism may be that inter-lamellar cohesion allows shear stress transfer between adjacent annular layers. The aim of the present study was therefore to investigate the importance of inter-lamellar shear resistance in the intervertebral disc. This work found that inter-lamellar shear resistance has a strong influence on the compressive stiffness of the intervertebral disc, with a change in interface condition from tied (no slip) to frictionless (no shear resistance) reducing disc compressive stiffness by 40%. However, it appears that substantial inter-lamellar shear resistance is present in the bovine tail disc. Decreases in inter-lamellar shear resistance due to degradation of bridging collagenous or elastic fibre structures could therefore be an important part of the process of disc degeneration.

Predicting surgical outcomes for paediatric spinal deformity patients using subject specific finite element models

INTRODUCTION Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. Our novel patient specific modelling software creates individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We have recently applied the model to investigate the influence of increasing magnitudes of surgically applied corrective force on predicted deformity correction...

Sequential MRI reveals individual level deformities in the growing scoliotic spine that are not seen clinically by the cobb angle

Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. Sequential MRI data in this project showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

Association study of genes related to bone formation and resorption and the extent of radiographic change in ankylosing spondylitis

Objective: To identify genetic associations with severity of radiographic damage in ankylosing spondylitis (AS). Method: We studied 1537 AS cases of European descent; all fulfilled the modified New York Criteria. Radiographic severity was assessed from digitised lateral radiographs of the cervical and lumbar spine using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). A two-phase genotyping design was used. In phase 1, 498 single nucleotide polymorphisms (SNPs) were genotyped in 688 cases; these were selected to capture >90% of the common haplotypic variation in the exons, exon-intron boundaries, and 5 kb flanking DNA in the 5' and 3' UTR of 74 genes involved in anabolic or catabolic bone pathways. In phase 2, 15 SNPs exhibiting p<0.05 were genotyped in a further cohort of 830 AS cases; results were analysed both separately and in combination with the discovery phase data. Association was tested by contingency tables after separating the samples into 'mild' and 'severe' groups, defined as the bottom and top 40% by mSASSS, adjusted for gender and disease duration. Results: Experiment-wise association was observed with the SNP rs8092336 (combined OR 0.32, p=1.2×10-5), which lies within RANK (receptor activator of NF?B), a gene involved in osteoclastogenesis, and in the interaction between T cells and dendritic cells. Association was also found with the SNP rs1236913 in PTGS1 (prostaglandin-endoperoxide synthase 1, cyclooxygenase 1), giving an OR of 0.53 (p=2.6×10-3). There was no observed association between radiographic severity and HLA-B*27. Conclusions: These findings support roles for bone resorption and prostaglandins pathways in the osteoproliferative changes in AS.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Life-course determinants of bone mass in young adults from a transitional rural community in India: The Andhra Pradesh Children and Parents Study (APCAPS)

Background: Undernutrition and physical inactivity are both associated with lower bone mass. Objective: This study aimed to investigate the combined effects of early-life undernutrition and urbanized lifestyles in later life on bone mass accrual in young adults from a rural community in India that is undergoing rapid socioeconomic development. Design: This was a prospective cohort study of participants of the Hyderabad Nutrition Trial (1987–1990), which offered balanced protein-calorie supplementation to pregnant women and preschool children younger than 6 y in the intervention villages. The 2009–2010 follow-up study collected data on current anthropometric measures, bone mineral density (BMD) measured by dual-energy X-ray absorptiometry, blood samples, diet, physical activity, and living standards of the trial participants (n = 1446, aged 18–23 y). Results: Participants were generally lean and had low BMD [mean hip BMD: 0.83 (women), 0.95 (men) g/cm2; lumbar spine: 0.86 (women), 0.93 (men) g/cm2]. In models adjusted for current risk factors, no strong evidence of a positive association was found between BMD and early-life supplementation. On the other hand, current lean mass and weight-bearing physical activity were positively associated with BMD. No strong evidence of an association was found between BMD and current serum 25-hydroxyvitamin D or dietary intake of calcium, protein, or calories. Conclusions: Current lean mass and weight-bearing physical activity were more important determinants of bone mass than was early-life undernutrition in this population. In transitional rural communities from low-income countries, promotion of physical activity may help to mitigate any potential adverse effects of early nutritional disadvantage.

Histopathological features of bone regeneration in a canine segmental ulnar defect model

Background Today, finding an ideal biomaterial to treat the large bone defects, delayed unions and non-unions remains a challenge for orthopaedic surgeions and researchers. Several studies have been carried out on the subject of bone regeneration, each having its own advantages. The present study has been designed in vivo to evaluate the effects of cellular auto-transplantation of tail vertebrae on healing of experimental critical bone defect in a dog model. Methods Six indigenous breeds of dog with 32 ± 3.6 kg average weight from both sexes (5 males and 1 female) received bilateral critical-sized ulnar segmental defects. After determining the health condition, divided to 2 groups: The Group I were kept as control I (n = 1) while in Group II (experimental group; n = 5) bioactive bone implants were inserted. The defects were implanted with either autogeneic coccygeal bone grafts in dogs with 3-4 cm diaphyseal defects in the ulna. Defects were stabilized with internal plate fixation, and the control defects were not stabilized. Animals were euthanized at 16 weeks and analyzed by histopathology. Results Histological evaluation of this new bone at sixteen weeks postoperatively revealed primarily lamellar bone, with the formation of new cortices and normal-appearing marrow elements. And also reformation cortical compartment and reconstitution of marrow space were observed at the graft-host interface together with graft resorption and necrosis responses. Finally, our data were consistent with the osteoconducting function of the tail autograft. Conclusions Our results suggested that the tail vertebrae autograft seemed to be a new source of autogenous cortical bone in order to supporting segmental long bone defects in dogs. Furthermore, cellular autotransplantation was found to be a successful replacement for the tail vertebrae allograft bone at 3-4 cm segmental defects in the canine mid- ulna. Clinical application using graft expanders or bone autotransplantation should be used carefully and requires further investigation.

Altered nociceptive, endocrine, and dorsal horn neuron responses in rats following a neonatal immune challenge

Summary The neonatal period is characterized by significant plasticity where the immune, endocrine, and nociceptive systems undergo fine-tuning and maturation. Painful experiences during this period can result in long-term alterations in the neurocircuitry underlying nociception, including increased sensitivity to mechanical or thermal stimuli. Less is known about the impact of neonatal exposure to mild inflammatory stimuli, such as lipopolysaccharide (LPS), on subsequent inflammatory pain responses. Here we examine the impact of neonatal LPS exposure on inflammatory pain sensitivity and HPA axis activity during the first three postnatal weeks. Wistar rats were injected with LPS (0.05 mg/kg IP, Salmonella enteritidis) or saline on postnatal days (PNDs) 3 and 5 and later subjected to the formalin test at PNDs 7, 13, and 22. One hour after formalin injection, blood was collected to assess corticosterone responses. Transverse spinal cord slices were also prepared for whole-cell patch clamp recording from lumbar superficial dorsal horn neurons (SDH). Brains were obtained at PND 22 and the hypothalamus was isolated to measure glucocorticoid (GR) and mineralocorticoid receptor (MR) transcript expression using qRT-PCR. Behavioural analyses indicate that at PND 7, no significant differences were observed between saline- or LPS-challenged rats. At PND 13, LPS-challenged rats exhibited enhanced licking (p < .01), and at PND 22, increased flinching in response to formalin injection (p < .05). LPS-challenged rats also displayed increased plasma corticosterone at PND 7 and PND 22 (p < .001) but not at PND 13 following formalin administration. Furthermore, at PND 22 neonatal LPS exposure induced decreased levels of GR mRNA and increased levels of MR mRNA in the hypothalamus. The intrinsic properties of SDH neurons were similar at PND 7 and PND 13. However, at PND 22, ipsilateral SDH neurons in LPS-challenged rats had a lower input resistance compared to their saline-challenged counterparts (p < .05). These data suggest neonatal LPS exposure produces developmentally regulated changes in formalin-induced behavioural responses, corticosterone levels, and dorsal horn neuron properties following noxious stimulation later in life. These findings highlight the importance of immune activation during the neonatal period in shaping pain sensitivity later in life. This programming involves both spinal cord neurons and the HPA axis.

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

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

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.