Trabecular bone microcrack accumulation in patients treated with bisphosphonates for durations up to 16 years.

This study quantified the length, number, and density of microcracks in bone from patients treated with bisphosphonates as a function of duration. Anterior iliac crest bone samples from 51 osteoporotic Caucasian females continuously treated with oral bisphosphonates for 1-16 years were obtained by bone biopsy. Samples were histologically processed and analyzed for bone area, microcrack number, and microcrack length. The analyses used statistical modeling and considered patient age, bone mineral density, bone volume/total volume, trabecular thickness, and bone turnover as potential covariates. Microcrack density (number of microcracks/total examined bone area) was linearly related (p = 0.018) to bisphosphonate treatment duration. None of the analyzed covariates contributed significantly to the observed relationship between microcrack density and bisphosphonate treatment duration. Observed increases in microcrack density with increasing bisphosphonate treatment duration is important because increasing levels of microcracks may not only affect bone remodeling but also reduce elastic modulus and are suspected to adversely affect other mechanical properties that may influence fracture risk. The present findings add to our prior results showing changes in bone material properties and modulus with bisphosphonate treatment duration and thereby provide a more comprehensive assessment of the relationship between bisphosphonate treatment duration and bone quality. Statement of Clinical Significance: The present findings provide information guiding clinical use of oral bisphosphonates for post-menopausal osteoporosis therapy.

Effect of Skeletal Paracrine Signals on the Proliferation of Interzone Cells.

OBJECTIVE: Articular cartilage in mammals has limited intrinsic capacity to repair structural defects, a fact that contributes to the chronic and progressive nature of osteoarthritis. In contrast, Mexican axolotl salamanders have demonstrated the remarkable ability to spontaneously and completely repair large joint cartilage lesions, a healing process that involves interzone cells in the intraarticular space. Furthermore, interzone tissue transplanted into skeletal defects in the axolotl salamander demonstrates a multi-differentiation potential. Cellular and molecular mechanisms of this repair process remain unclear. The objective of this study was to examine whether paracrine mitogenic signals are an important variable in the interaction between interzone cells and the skeletal microenvironment. DESIGN: The paracrine regulation of the proliferation of equine interzone cells was evaluated in an in vitro co-culture system. Cell viability and proliferation were measured in equine fetal interzone cells after exposure to conditioned medium from skeletal and nonskeletal primary cell lines. Steady-state expression was determined for genes encoding 37 putative mitogens secreted by cells that generated the conditioned medium. RESULTS: All experimental groups of conditioned media elicited a mitogenic response in interzone cells. Fetal anlage chondrocytes (P < 0.0001) and dermal fibroblasts (P < 0.0001) conditioned medium showed a significantly higher mitogenic potential compared with interzone cells. Conditioned medium from bone marrow-derived cells elicited a significantly higher proliferative response relative to that from young adult articular chondrocytes (P < 0.0001) or dermal fibroblasts (P < 0.0001). Sixteen genes had expression patterns consistent with the functional proliferation assays. CONCLUSIONS: The results indicate a mitogenic effect of skeletal paracrine signals on interzone cells.

Tall fescue endophyte effects on tolerance to water-deficit stress.

BACKGROUND: The endophytic fungus, Neotyphodium coenophialum, can enhance drought tolerance of its host grass, tall fescue. To investigate endophyte effects on plant responses to acute water deficit stress, we did comprehensive profiling of plant metabolite levels in both shoot and root tissues of genetically identical clone pairs of tall fescue with endophyte (E+) and without endophyte (E-) in response to direct water deficit stress. The E- clones were generated by treating E+ plants with fungicide and selectively propagating single tillers. In time course studies on the E+ and E- clones, water was withheld from 0 to 5 days, during which levels of free sugars, sugar alcohols, and amino acids were determined, as were levels of some major fungal metabolites. RESULTS: After 2-3 days of withholding water, survival and tillering of re-watered plants was significantly greater for E+ than E- clones. Within two to three days of withholding water, significant endophyte effects on metabolites manifested as higher levels of free glucose, fructose, trehalose, sugar alcohols, proline and glutamic acid in shoots and roots. The fungal metabolites, mannitol and loline alkaloids, also significantly increased with water deficit. CONCLUSIONS: Our results suggest that symbiotic N. coenophialum aids in survival and recovery of tall fescue plants from water deficit, and acts in part by inducing rapid accumulation of these compatible solutes soon after imposition of stress.

Identification of gene expression patterns using planned linear contrasts.

BACKGROUND: In gene networks, the timing of significant changes in the expression level of each gene may be the most critical information in time course expression profiles. With the same timing of the initial change, genes which share similar patterns of expression for any number of sampling intervals from the beginning should be considered co-expressed at certain level(s) in the gene networks. In addition, multiple testing problems are complicated in experiments with multi-level treatments when thousands of genes are involved. RESULTS: To address these issues, we first performed an ANOVA F test to identify significantly regulated genes. The Benjamini and Hochberg (BH) procedure of controlling false discovery rate (FDR) at 5% was applied to the P values of the F test. We then categorized the genes with a significant F test into 4 classes based on the timing of their initial responses by sequentially testing a complete set of orthogonal contrasts, the reverse Helmert series. For genes within each class, specific sequences of contrasts were performed to characterize their general 'fluctuation' shapes of expression along the subsequent sampling time points. To be consistent with the BH procedure, each contrast was examined using a stepwise Studentized Maximum Modulus test to control the gene based maximum family-wise error rate (MFWER) at the level of alphanew determined by the BH procedure. We demonstrated our method on the analysis of microarray data from murine olfactory sensory epithelia at five different time points after target ablation. CONCLUSION: In this manuscript, we used planned linear contrasts to analyze time-course microarray experiments. This analysis allowed us to characterize gene expression patterns based on the temporal order in the data, the timing of a gene's initial response, and the general shapes of gene expression patterns along the subsequent sampling time points. Our method is particularly suitable for analysis of microarray experiments in which it is often difficult to take sufficiently frequent measurements and/or the sampling intervals are non-uniform.

Analysis of oligonucleotide array experiments with repeated measures using mixed models.

BACKGROUND: Two or more factor mixed factorial experiments are becoming increasingly common in microarray data analysis. In this case study, the two factors are presence (Patients with Alzheimer's disease) or absence (Control) of the disease, and brain regions including olfactory bulb (OB) or cerebellum (CER). In the design considered in this manuscript, OB and CER are repeated measurements from the same subject and, hence, are correlated. It is critical to identify sources of variability in the analysis of oligonucleotide array experiments with repeated measures and correlations among data points have to be considered. In addition, multiple testing problems are more complicated in experiments with multi-level treatments or treatment combinations. RESULTS: In this study we adopted a linear mixed model to analyze oligonucleotide array experiments with repeated measures. We first construct a generalized F test to select differentially expressed genes. The Benjamini and Hochberg (BH) procedure of controlling false discovery rate (FDR) at 5% was applied to the P values of the generalized F test. For those genes with significant generalized F test, we then categorize them based on whether the interaction terms were significant or not at the alpha-level (alphanew = 0.0033) determined by the FDR procedure. Since simple effects may be examined for the genes with significant interaction effect, we adopt the protected Fisher's least significant difference test (LSD) procedure at the level of alphanew to control the family-wise error rate (FWER) for each gene examined. CONCLUSIONS: A linear mixed model is appropriate for analysis of oligonucleotide array experiments with repeated measures. We constructed a generalized F test to select differentially expressed genes, and then applied a specific sequence of tests to identify factorial effects. This sequence of tests applied was designed to control for gene based FWER.

Statistical implications of pooling RNA samples for microarray experiments.

BACKGROUND: Microarray technology has become a very important tool for studying gene expression profiles under various conditions. Biologists often pool RNA samples extracted from different subjects onto a single microarray chip to help defray the cost of microarray experiments as well as to correct for the technical difficulty in getting sufficient RNA from a single subject. However, the statistical, technical and financial implications of pooling have not been explicitly investigated. RESULTS: Modeling the resulting gene expression from sample pooling as a mixture of individual responses, we derived expressions for the experimental error and provided both upper and lower bounds for its value in terms of the variability among individuals and the number of RNA samples pooled. Using "virtual" pooling of data from real experiments and computer simulations, we investigated the statistical properties of RNA sample pooling. Our study reveals that pooling biological samples appropriately is statistically valid and efficient for microarray experiments. Furthermore, optimal pooling design(s) can be found to meet statistical requirements while minimizing total cost. CONCLUSIONS: Appropriate RNA pooling can provide equivalent power and improve efficiency and cost-effectiveness for microarray experiments with a modest increase in total number of subjects. Pooling schemes in terms of replicates of subjects and arrays can be compared before experiments are conducted.