Alexandre S. Stephens

Temporal pattern of gene expression and histology of stress fracture healing

Loading of the rat ulna is an ideal model to examine stress fracture healing. The aim of this study was to undertake a detailed examination of the histology, histomorphometry and gene expression of the healing and remodelling process initiated by fatigue loading of the rat ulna. Ulnae were harvested 1, 2, 4, 6, 8, and 10 weeks following creation of a stress fracture. Stress fracture healing involved direct remodelling that progressed along the fracture line as well as woven bone proliferation at the site of the fracture. Histomorphometry demonstrated rapid progression of basic multicellular units from 1 to 4 weeks with significant slowing down of healing by 10 weeks after loading. Quantitative PCR was performed at 4 hours, 24 hours, 4 days, 7 days, and 14 days after loading. Gene expression was compared to an unloaded control group. At 4 hours after fracture, there was a marked 220-fold increase (P<0.0001) in expression of IL-6. There were also prominent peak increases in mRNA expression for OPG, COX-2, and VEGF (all P<0.0001). At 24 hours, there was a peak increase in mRNA expression for IL-11 (73-fold increase, P<0.0001). At 4 days, there was a significant increase in mRNA expression for Bcl-2, COX-1, IGF-1, OPN, and SDF-1. At 7 days, there was significantly increased mRNA expression of RANKL and OPN. Prominent, upregulation of COX-2, VEGF, OPG, SDF-1, BMP-2, and SOST prior to peak expression of RANKL indicates the importance of these factors in mediating directed remodelling of the fracture line. Dramatic, early upregulation of IL-6 and IL-11 demonstrate their central role in initiating signalling events for remodelling and stress fracture healing.

Internal control genes for quantitative RT-PCR expression analysis in mouse osteoblasts, osteoclasts and macrophages

BackgroundReal-time quantitative RT-PCR (qPCR) is a powerful technique capable of accurately quantitating mRNA expression levels over a large dynamic range. This makes qPCR the most widely used method for studying quantitative gene expression. An important aspect of qPCR is selecting appropriate controls or normalization factors to account for any differences in starting cDNA quantities between samples during expression studies. Here, we report on the selection of a concise set of housekeeper genes for the accurate normalization of quantitative gene expression data in differentiating osteoblasts, osteoclasts and macrophages. We implemented the use of geNorm, an algorithm that determines the suitability of genes to function as housekeepers by assessing expression stabilities. We evaluated the expression stabilities of 18S, ACTB, B2M, GAPDH, HMBS and HPRT1 genes.FindingsOur analyses revealed that 18S and GAPDH were regulated during osteoblast differentiation and are not suitable for use as reference genes. The most stably expressed genes in osteoblasts were ACTB, HMBS and HPRT1 and their geometric average constitutes a suitable normalization factor upon which gene expression data can be normalized. In macrophages, 18S and GAPDH were the most variable genes while HMBS and B2M were the most stably expressed genes. The geometric average of HMBS and B2M expression levels forms a suitable normalization factor to account for potential differences in starting cDNA quantities during gene expression analysis in macrophages. The expression stabilities of the six candidate reference genes in osteoclasts were, on average, more variable than that observed in macrophages but slightly less variable than those seen in osteoblasts. The two most stably expressed genes in osteoclasts were HMBS and B2M and the genes displaying the greatest levels of variability were 18S and GAPDH. Notably, 18S and GAPDH were the two most variably expressed control genes in all three cell types. The geometric average of HMBS, B2M and ACTB creates an appropriate normalization factor for gene expression studies in osteoclasts.ConclusionWe have identified concise sets of genes suitable to use as normalization factors for quantitative real-time RT-PCR gene expression studies in osteoblasts, osteoclasts and macrophages.

Association of Functionally Different RUNX2 P2 Promoter Alleles with BMD

RUNX2 gene SNPs were genotyped in subjects from the upper and lower deciles of age‐ and weight‐adjusted femoral neck BMD. Of 16 SNPs in RUNX2 and its two promoters (P1 and P2), only SNPs in the P2 promoter were significantly associated with BMD. These P2 promoter SNPs were functionally different in gel‐shift and promoter activity assays.

Glutamine repeat variants in human RUNX2 associated with decreased femoral neck BMD, broadband ultrasound attenuation and target gene transactivation.

RUNX2 is an essential transcription factor required for skeletal development and cartilage formation. Haploinsufficiency of RUNX2 leads to cleidocranial displaysia (CCD) a skeletal disorder characterised by gross dysgenesis of bones particularly those derived from intramembranous bone formation. A notable feature of the RUNX2 protein is the polyglutamine and polyalanine (23Q/17A) domain coded by a repeat sequence. Since none of the known mutations causing CCD characterised to date map in the glutamine repeat region, we hypothesised that Q-repeat mutations may be related to a more subtle bone phenotype. We screened subjects derived from four normal populations for Q-repeat variants. A total of 22 subjects were identified who were heterozygous for a wild type allele and a Q-repeat variant allele: (15Q, 16Q, 18Q and 30Q). Although not every subject had data for all measures, Q-repeat variants had a significant deficit in BMD with an average decrease of 0.7SD measured over 12 BMD-related parameters (p = 0.005). Femoral neck BMD was measured in all subjects (−0.6SD, p = 0.0007). The transactivation function of RUNX2 was determined for 16Q and 30Q alleles using a reporter gene assay. 16Q and 30Q alleles displayed significantly lower transactivation function compared to wild type (23Q). Our analysis has identified novel Q-repeat mutations that occur at a collective frequency of about 0.4%. These mutations significantly alter BMD and display impaired transactivation function, introducing a new class of functionally relevant RUNX2 mutants.

Polyalanine Repeat Polymorphism in RUNX2 Is Associated with Site-Specific Fracture in Post-Menopausal Females

Runt related transcription factor 2 (RUNX2) is a key regulator of osteoblast differentiation. Several variations within the RUNX2 gene have been found to be associated with significant changes in BMD, which is a major risk factor for fracture. In this study we report that an 18 bp deletion within the polyalanine tract (17A>11A) of RUNX2 is significantly associated with fracture. Carriers of the 11A allele were found to be nearly twice as likely to have sustained fracture. Within the fracture category, there was a significant tendency of 11A carriers to present with fractures of distal radius and bones of intramembranous origin compared to bones of endochondral origin (p = 0.0001). In a population of random subjects, the 11A allele was associated with decreased levels of serum collagen cross links (CTx, p = 0.01), suggesting decreased bone turnover. The transactivation function of the 11A allele showed a minor quantitative decrease. Interestingly, we found no effect of the 11A allele on BMD at multiple skeletal sites. These findings suggest that the 11A allele is a biologically relevant polymorphism that influences serum CTx and confers enhanced fracture risk in a site-selective manner related to intramembranous bone ossification.

Myocyte Enhancer Factor 2C, an Osteoblast Transcription Factor Identified by Dimethyl Sulfoxide (DMSO)-enhanced Mineralization

Rapid mineralization of cultured osteoblasts could be a useful characteristic in stem cell-mediated therapies for fracture and other orthopedic problems. Dimethyl sulfoxide (DMSO) is a small amphipathic solvent molecule capable of stimulating cell differentiation. We report that, in primary human osteoblasts, DMSO dose-dependently enhanced the expression of osteoblast differentiation markers alkaline phosphatase activity and extracellular matrix mineralization. Furthermore, similar DMSO-mediated mineralization enhancement was observed in primary osteoblast-like cells differentiated from mouse mesenchymal cells derived from fat, a promising source of starter cells for cell-based therapy. Using a convenient mouse pre-osteoblast model cell line MC3T3-E1, we further investigated this phenomenon showing that numerous osteoblast-expressed genes were elevated in response to DMSO treatment and correlated with enhanced mineralization. Myocyte enhancer factor 2c (Mef2c) was identified as the transcription factor most induced by DMSO, among the numerous DMSO-induced genes, suggesting a role for Mef2c in osteoblast gene regulation. Immunohistochemistry confirmed expression of Mef2c in osteoblast-like cells in mouse mandible, cortical, and trabecular bone. shRNAi-mediated Mef2c gene silencing resulted in defective osteoblast differentiation, decreased alkaline phosphatase activity, and matrix mineralization and knockdown of osteoblast specific gene expression, including osteocalcin and bone sialoprotein. A flow on knockdown of bone-specific transcription factors, Runx2 and osterix by shRNAi knockdown of Mef2c, suggests that Mef2c lies upstream of these two important factors in the cascade of gene expression in osteoblasts.

Novel Target Genes of RUNX2 Transcription Factor and 1,25‐Dihydroxyvitamin D3

The RUNX2 transcription factor is indispensable for skeletal development and controls bone formation by acting as a signaling hub and transcriptional regulator to coordinate target gene expression. A signaling partner of RUNX2 is the nuclear vitamin D receptor (VDR) that becomes active when bound by its ligand 1,25‐dihydroxyvitamin D3 (VD3). RUNX2 and VDR unite to cooperatively regulate the expression of numerous genes. In this study, we overexpressed RUNX2 in NIH3T3 fibroblasts concomitantly treated with VD3 and show that RUNX2 alone, or in combination with VD3, failed to promote an osteoblastic phenotype in NIH3T3 cells. However, the expression of numerous osteoblast‐related genes was up‐regulated by RUNX2 and large‐scale gene expression profiling using microarrays identified over 800 transcripts that displayed a twofold of greater change in expression in response to RUNX2 overexpression or VD3 treatment. Functional analysis using gene ontology (GO) revealed GO terms for ossification, cellular motility, biological adhesion, and chromosome organization were enriched in the pool of genes regulated by RUNX2. For the set of genes whose expression was modulated by VD3, the GO terms response to hormone stimulus, chemotaxis, and metalloendopeptidase activity where overrepresented. Our study provides a functional insight into the consequences of RUNX2 overexpression and VD3 treatment in NIH3T3 cells in addition to identifying candidate genes whose expression is controlled by either factor individually or through their functional cooperation. J. Cell. Biochem. 115: 1594–1608, 2014.

Characterisation of CMP-Sialic Acid Transporter Substrate Recognition

CMP-sialic acid transporter: We report an in-depth, multidisciplinary, structural study that has identified the amino acid residues intimately involved in CMP-sialic acid transporter (CST) substrate specificity. Our data provide a significant contribution towards a better understanding the structure-function relationship of this important family of transporters and the rational design of CST inhibitors.

Life expectancy estimation in small administrative areas with non-uniform population sizes: application to Australian New South Wales local government areas

Objective To determine a practical approach for deriving life expectancy estimates in Australian New South Wales local government areas which display a large diversity in population sizes. Design Population-based study utilising mortality and estimated residential population data. Setting 153 local government areas in New South Wales, Australia. Outcome measures Key performance measures of Chiang II, Silcocks, adjusted Chiang II and Bayesian random effects model methodologies of life expectancy estimation including agreement analysis of life expectancy estimates and comparison of estimate SEs. Results Chiang II and Silcocks methods produced almost identical life expectancy estimates across a large range of population sizes but calculation failures and excessively large SEs limited their use in small populations. A population of 25 000 or greater was required to estimate life expectancy with SE of 1 year or less using adjusted Chiang II (a composite of Chiang II and Silcocks methods). Data aggregation offered some remedy for extending the use of adjusted Chiang II in small populations but reduced estimate currency. A recently developed Bayesian random effects model utilising the correlation in mortality rates between genders, age groups and geographical areas markedly improved the precision of life expectancy estimates in small populations. Conclusions We propose a hybrid approach for the calculation of life expectancy using the Bayesian random effects model in populations of 25 000 or lower permitting the precise derivation of life expectancy in small populations. In populations above 25 000, we propose the use of adjusted Chiang II to guard against violations of spatial correlation, to benefit from a widely accepted method that is simpler to communicate to local health authorities and where its slight inferior performance compared with the Bayesian approach is of minor practical significance.

Diagnosis of fetal growth restriction in perinatal deaths using brain to liver weight ratios.

Summary We determined brain to liver weight ratio (BLWR) thresholds for fetal growth restriction (FGR) using autopsy information on 395 perinatal deaths comprising stillborn babies who died during labour and neonatal deaths. FGR was defined using two methods: (1) birth weight for gestational age (WGA) less than the 10th percentile; and (2) WGA less than the 10th percentile or discordant birth weight/length. The association between BLWR and FGR was investigated using odds ratios, and classification statistics were calculated for a range of BLWR thresholds. Using WGA, 84 cases (21.3%) were FGR and a further 15 cases (n = 99, 25%) had discordant birth weight/length. The BLWR ranged from 1.02 to 7.30 and was positively associated with FGR. BLWR was not associated with FGR for babies with congenital central nervous system or chromosomal abnormalities. Excluding these, for FGR defined using WGA and discordant birth weight/length, a BLWR threshold of 5.0 was 100% predictive of FGR. A BLWR threshold of 3.0 for babies over 28 weeks gestation and 3.7 for more preterm babies optimised case detection while minimising missed and false positive cases. Additional evidence of FGR should be sought for babies with a BLWR of less than 5.0 to confirm FGR.