A Functional-Based Distribution Diagnostic for a Linear Model with Correlated Outcomes: Technical Report

Despite the widespread popularity of linear models for correlated outcomes (e.g. linear mixed modesl and time series models), distribution diagnostic methodology remains relatively underdeveloped in this context. In this paper we present an easy-to-implement approach that lends itself to graphical displays of model fit. Our approach involves multiplying the estimated marginal residual vector by the Cholesky decomposition of the inverse of the estimated marginal variance matrix. Linear functions or the resulting "rotated" residuals are used to construct an empirical cumulative distribution function (ECDF), whose stochastic limit is characterized. We describe a resampling technique that serves as a computationally efficient parametric bootstrap for generating representatives of the stochastic limit of the ECDF. Through functionals, such representatives are used to construct global tests for the hypothesis of normal margional errors. In addition, we demonstrate that the ECDF of the predicted random effects, as described by Lange and Ryan (1989), can be formulated as a special case of our approach. Thus, our method supports both omnibus and directed tests. Our method works well in a variety of circumstances, including models having independent units of sampling (clustered data) and models for which all observations are correlated (e.g., a single time series).

Establishment of a novel intervertebral disc/endplate culture model: analysis of an ex vivo in vitro whole-organ rabbit culture system

STUDY DESIGN: Ex vivo in vitro study evaluating a novel intervertebral disc/endplate culture system. OBJECTIVES: To establish a whole-organ intervertebral disc culture model for the study of disc degeneration in vitro, including the characterization of basic cell and organ function. SUMMARY OF BACKGROUND DATA: With current in vivo models for the study of disc and endplate degeneration, it remains difficult to investigate the complex disc metabolism and signaling cascades. In contrast, more controlled but simplified in vitro systems using isolated cells or disc fragments are difficult to culture due to the unconstrained conditions, with often-observed cell death or cell dedifferentiation. Therefore, there is a demand for a controlled culture model with preserved cell function that offers the possibility to investigate disc and endplate pathologies in a structurally intact organ. METHODS: Naturally constrained intervertebral disc/endplate units from rabbits were cultured in multi-well plates. Cell viability, metabolic activity, matrix composition, and matrix gene expression profile were monitored using the Live/Dead cell viability test (Invitrogen, Basel, Switzerland), tetrazolium salt reduction (WST-8), proteoglycan and deoxyribonucleic acid quantification assays, and quantitative polymerase chain reaction. RESULTS: Viability and organ integrity were preserved for at least 4 weeks, while proteoglycan and deoxyribonucleic acid content decreased slightly, and matrix genes exhibited a degenerative profile with up-regulation of type I collagen and suppression of collagen type II and aggrecan genes. Additionally, cell metabolic activity was reduced to one third of the initial value. CONCLUSIONS: Naturally constrained intervertebral rabbit discs could be cultured for several weeks without losing cell viability. Structural integrity and matrix composition were retained. However, the organ responded to the artificial environment with a degenerative gene expression pattern and decreased metabolic rate. Therefore, the described system serves as a promising in vitro model to study disc degeneration in a whole organ.

Cyclooxygenase-2 inhibition selectively attenuates bone morphogenetic protein-6 synthesis and bone formation during guided tissue regeneration in a rat model

OBJECTIVES: Bone formation during guided tissue regeneration is a tightly regulated process involving cells, extracellular matrix and growth factors. The aims of this study were (i) to examine the expression of cyclooxygenase-2 (COX-2) during bone regeneration and (ii) the effects of selective COX-2 inhibition on osseous regeneration and growth factor expression in the rodent femur model. MATERIAL AND METHODS: A standardized transcortical defect of 5 x 1.5 mm was prepared in the femur of 12 male rats and a closed half-cylindrical titanium chamber was placed over the defect. The expression of COX-2 and of platelet-derived growth factor-B (PDGF-B), bone morphogenetic protein-6 (BMP-6) and insulin-like growth factor-I/II (IGF-I/II) was analyzed at Days 3, 7, 21 and 28 semiquantitatively by reverse transcriptase-polymerase chain reaction and immunohistochemistry. The effects of COX-2 inhibition by intraperitoneal injection of NS-398 (3 mg/kg/day) were analyzed in five additional animals sacrificed at Day 14. RESULTS: Histomorphometry revealed that new bone formation occurred in the cortical defect area as well as in the supracortical region, i.e. region within the chamber by Day 7 and increased through Day 28. Immunohistochemical evidence of COX-2 and PDGF-B levels were observed early (i.e. Day 3) and decreased rapidly by Day 7. BMP-6 expression was maximal at Day 3 and slowly declined by Day 28. In contrast, IGF-I/II expression gradually increased during the 28-day period. Systemic administration NS-398 caused a statistically significant reduction (P<0.05) in new bone formation (25-30%) and was associated with a statistically significant reduction in BMP-6 protein and mRNA expression (50% and 65% at P<0.05 and P<0.01, respectively). PDGF-B mRNA or protein expression was not affected by NS-398 treatment. CONCLUSION: COX-2 inhibition resulted in reduced BMP-6 expression and impaired osseous regeneration suggesting an important role for COX-2-induced signaling in BMP synthesis and new bone formation.

Cholesky Residuals for Assessing Normal Errors in a Linear Model with Correlated Outcomes: Technical Report

Despite the widespread popularity of linear models for correlated outcomes (e.g. linear mixed models and time series models), distribution diagnostic methodology remains relatively underdeveloped in this context. In this paper we present an easy-to-implement approach that lends itself to graphical displays of model fit. Our approach involves multiplying the estimated margional residual vector by the Cholesky decomposition of the inverse of the estimated margional variance matrix. The resulting "rotated" residuals are used to construct an empirical cumulative distribution function and pointwise standard errors. The theoretical framework, including conditions and asymptotic properties, involves technical details that are motivated by Lange and Ryan (1989), Pierce (1982), and Randles (1982). Our method appears to work well in a variety of circumstances, including models having independent units of sampling (clustered data) and models for which all observations are correlated (e.g., a single time series). Our methods can produce satisfactory results even for models that do not satisfy all of the technical conditions stated in our theory.

On the Accelerated Failure Time Model for Current Status and Interval Censored Data

This paper introduces a novel approach to making inference about the regression parameters in the accelerated failure time (AFT) model for current status and interval censored data. The estimator is constructed by inverting a Wald type test for testing a null proportional hazards model. A numerically efficient Markov chain Monte Carlo (MCMC) based resampling method is proposed to simultaneously obtain the point estimator and a consistent estimator of its variance-covariance matrix. We illustrate our approach with interval censored data sets from two clinical studies. Extensive numerical studies are conducted to evaluate the finite sample performance of the new estimators.

A BAYESIAN HIERARCHICAL MODEL FOR CONSTRAINED DISTRIBUTED LAG FUNCTIONS: ESTIMATING THE TIME COURSE OF HOSPITALIZATION ASSOCIATED WITH AIR POLLUTION EXPOSURE

Numerous time series studies have provided strong evidence of an association between increased levels of ambient air pollution and increased levels of hospital admissions, typically at 0, 1, or 2 days after an air pollution episode. An important research aim is to extend existing statistical models so that a more detailed understanding of the time course of hospitalization after exposure to air pollution can be obtained. Information about this time course, combined with prior knowledge about biological mechanisms, could provide the basis for hypotheses concerning the mechanism by which air pollution causes disease. Previous studies have identified two important methodological questions: (1) How can we estimate the shape of the distributed lag between increased air pollution exposure and increased mortality or morbidity? and (2) How should we estimate the cumulative population health risk from short-term exposure to air pollution? Distributed lag models are appropriate tools for estimating air pollution health effects that may be spread over several days. However, estimation for distributed lag models in air pollution and health applications is hampered by the substantial noise in the data and the inherently weak signal that is the target of investigation. We introduce an hierarchical Bayesian distributed lag model that incorporates prior information about the time course of pollution effects and combines information across multiple locations. The model has a connection to penalized spline smoothing using a special type of penalty matrix. We apply the model to estimating the distributed lag between exposure to particulate matter air pollution and hospitalization for cardiovascular and respiratory disease using data from a large United States air pollution and hospitalization database of Medicare enrollees in 94 counties covering the years 1999-2002.

Enhanced early chondrogenesis in articular defects following arthroscopic mesenchymal stem cell implantation in an equine model

Mesenchymal stem cells (MSCs) provide an important source of pluripotent cells for musculoskeletal tissue repair. This study examined the impact of MSC implantation on cartilage healing characteristics in a large animal model. Twelve full-thickness 15-mm cartilage lesions in the femoropatellar articulations of six young mature horses were repaired by injection of a self-polymerizing autogenous fibrin vehicle containing mesenchymal stem cells, or autogenous fibrin alone in control joints. Arthroscopic second look and defect biopsy was obtained at 30 days, and all animals were euthanized 8 months after repair. Cartilage repair tissue and surrounding cartilage were assessed by histology, histochemistry, collagen type I and type II immunohistochemistry, collagen type II in situ hybridization, and matrix biochemical assays. Arthroscopic scores for MSC-implanted defects were significantly improved at the 30-day arthroscopic assessment. Biopsy showed MSC-implanted defects contained increased fibrous tissue with several defects containing predominantly type II collagen. Long-term assessment revealed repair tissue filled grafted and control lesions at 8 months, with no significant difference between stem cell-treated and control defects. Collagen type II and proteoglycan content in MSC-implanted and control defects were similar. Mesenchymal stem cell grafts improved the early healing response, but did not significantly enhance the long-term histologic appearance or biochemical composition of full-thickness cartilage lesions.

Essential role of choline for pneumococcal virulence in an experimental model of meningitis

Objectives: The goal of the present study was to elucidate the contribution of the newly recognized virulence factor choline to the pathogenesis of Streptococcus pneumoniae in an animal model of meningitis. Results: The choline containing strain D39Cho(-) and its isogenic choline-free derivative D39Cho(-)licA64 -each expressing the capsule polysaccharide 2 - were introduced intracisternally at an inoculum size of 10(3) CFU into 11 days old Wistar rats. During the first 8 h post infection both strains multiplied and stimulated a similar immune response that involved expression of high levels of proinflammatory cytokines, the matrix metalloproteinase 9 (MMP-9), IL-10, and the influx of white blood cells into the CSF. Virtually identical immune response was also elicited by intracisternal inoculation of 10(7) CFU equivalents of either choline-containing or choline-free cell walls. At sampling times past 8 h strain D39Cho(-) continued to replicate accompanied by an intense inflammatory response and strong granulocytic pleiocytosis. Animals infected with D39Cho(-) died within 20 h and histopathology revealed brain damage in the cerebral cortex and hippocampus. In contrast, the initial immune response generated by the choline-free strain D39Cho(-)licA64 began to decline after the first 8 h accompanied by elimination of the bacteria from the CSF in parallel with a strong WBC response peaking at 8 h after infection. All animals survived and there was no evidence for brain damage. Conclusion: Choline in the cell wall is essential for pneumococci to remain highly virulent and survive within the host and establish pneumococcal meningitis.

Matrix metalloproteinases contribute to brain damage in experimental pneumococcal meningitis

The present study was performed to evaluate the role of matrix metalloproteinases (MMP) in the pathogenesis of the inflammatory reaction and the development of neuronal injury in a rat model of bacterial meningitis. mRNA encoding specific MMPs (MMP-3, MMP-7, MMP-8, and MMP-9) and the inflammatory cytokine tumor necrosis factor alpha (TNF-alpha) were significantly (P < 0.04) upregulated, compared to the beta-actin housekeeping gene, in cortical homogenates at 20 h after infection. In parallel, concentrations of MMP-9 and TNF-alpha in cerebrospinal fluid (CSF) were significantly increased in rats with bacterial meningitis compared to uninfected animals (P = 0.002) and showed a close correlation (r = 0.76; P < 0. 001). Treatment with a hydroxamic acid-type MMP inhibitor (GM6001; 65 mg/kg intraperitoneally every 12 h) beginning at the time of infection significantly lowered the MMP-9 (P < 0.02) and TNF-alpha (P < 0.02) levels in CSF. Histopathology at 25.5 +/- 5.7 h after infection showed neuronal injury (median [range], 3.5% [0 to 17.5%] of the cortex), which was significantly (P < 0.01) reduced to 0% (0 to 10.8%) by GM6001. This is the first report to demonstrate that MMPs contribute to the development of neuronal injury in bacterial meningitis and that inhibition of MMPs may be an effective approach to prevent brain damage as a consequence of the disease.

Experimental model to evaluate in vivo and in vitro cartilage MR imaging by means of histological analyses

OBJECTIVES: Implementation of an experimental model to compare cartilage MR imaging by means of histological analyses. MATERIAL AND METHODS: MRI was obtained from 4 patients expecting total knee replacement at 1.5 and/or 3T prior surgery. The timeframe between pre-op MRI and knee replacement was within two days. Resected cartilage-bone samples were tagged with Ethi((R))-pins to reproduce the histological cutting course. Pre-operative scanning at 1.5T included following parameters for fast low angle shot (FLASH: TR/TE/FA=33ms/6ms/30 degrees , BW=110kHz, 120mmx120mm FOV, 256x256 matrix, 0.65mm slice-thickness) and double echo steady state (DESS: TR/TE/FA=23.7ms/6.9ms/40 degrees , BW=130kHz, 120x120mm FOV, 256x256 matrix, 0.65mm slice-thickness). At 3T, scan parameters were: FLASH (TR/TE/FA=12.2ms/5.1ms/10 degrees , BW=130kHz, 170x170mm FOV, 320x320, 0.5mm slice-thickness) and DESS (TR/TE/FA=15.6ms/4.5ms/25 degrees , BW=200kHz, 135mmx150mm FOV, 288x320matrix, 0.5mm slice-thickness). Imaging of the specimens was done the same day at 1.5T. MRI (Noyes) and histological (Mankin) score scales were correlated using the paired t-test. Sensitivity and specificity for the detection of different grades of cartilage degeneration were assessed. Inter-reader and intra-reader reliability was determined using Kappa analysis. RESULTS: Low correlation (sensitivity, specificity) was found for both sequences in normal to mild Mankin grades. Only moderate to severe changes were diagnosed with higher significance and specificity. The use of higher field-strengths was advantageous for both protocols with sensitivity values ranging from 13.6% to 93.3% (FLASH) and 20.5% to 96.2% (DESS). Kappa values ranged from 0.488 to 0.944. CONCLUSIONS: Correlating MR images with continuous histological slices was feasible by using three-dimensional imaging, multi-planar-reformat and marker pins. The capability of diagnosing early cartilage changes with high accuracy could not be proven for both FLASH and DESS.

Origin of renal myofibroblasts in the model of unilateral ureter obstruction in the rat

Tubulo-interstitial fibrosis is a constant feature of chronic renal failure and it is suspected to contribute importantly to the deterioration of renal function. In the fibrotic kidney there exists, besides normal fibroblasts, a large population of myofibroblasts, which are supposedly responsible for the increased production of intercellular matrix. It has been proposed that myofibroblasts in chronic renal failure originate from the transformation of tubular cells via epithelial-mesenchymal transition (EMT) or from infiltration by bone marrow-derived precursors. Little attention has been paid to the possibility of a transformation of resident fibroblasts into myofibroblasts in renal fibrosis. Therefore we examined the fate of resident fibroblasts in the initial phase of renal fibrosis in the classical model of unilateral ureter obstruction (UUO) in the rat. Rats were perfusion-fixed on days 1, 2, 3 and 4 after ligature of the right ureter. Starting from 1 day of UUO an increasing expression of alpha-smooth muscle actin (alphaSMA) in resident fibroblasts was revealed by immunofluorescence and confirmed by the observation of bundles of microfilaments and webs of intermediate filaments in the electron microscope. Inversely, there was a decreased expression of 5'-nucleotidase (5'NT), a marker of renal cortical fibroblasts. The RER became more voluminous, suggesting an increased synthesis of matrix. Intercellular junctions, a characteristic feature of myofibroblasts, became more frequent. The mitotic activity in fibroblasts was strongly increased. Renal tubules underwent severe regressive changes but the cells retained their epithelial characteristics and there was no sign of EMT. In conclusion, after ureter ligature, resident peritubular fibroblasts proliferated and they showed progressive alterations, suggesting a transformation in myofibroblasts. Thus the resident fibroblasts likely play a central role in fibrosis in that model.

High temperature, low cycle fatigue of a hybrid particulate/fiber aluminum metal-matrix composite

The effect of shot particles on the high temperature, low cycle fatigue of a hybrid fiber/particulate metal-matrix composite (MMC) was studied. Two hybrid composites with the general composition A356/35%SiC particle/5%Fiber (one without shot) were tested. It was found that shot particles acting as stress concentrators had little effect on the fatigue performance. It appears that fibers with a high silica content were more likely to debond from the matrix. Final failure of the composite was found to occur preferentially in the matrix. SiC particles fracture progressively during fatigue testing, leading to higher stress in the matrix, and final failure by matrix overload. A continuum mechanics based model was developed to predict failure in fatigue based on the tensile properties of the matrix and particles. By accounting for matrix yielding and recovery, composite creep and particle strength distribution, failure of the composite was predicted.

Model to describe the mode I fracture of steel fiber reinforced ultra-high performance concrete

Ultra-high performance fiber reinforced concrete (UHPFRC) has arisen from the implementation of a variety of concrete engineering and materials science concepts developed over the last century. This material offers superior strength, serviceability, and durability over its conventional counterparts. One of the most important differences for UHPFRC over other concrete materials is its ability to resist fracture through the use of randomly dispersed discontinuous fibers and improvements to the fiber-matrix bond. Of particular interest is the materials ability to achieve higher loads after first crack, as well as its high fracture toughness. In this research, a study of the fracture behavior of UHPFRC with steel fibers was conducted to look at the effect of several parameters related to the fracture behavior and to develop a fracture model based on a non-linear curve fit of the data. To determine this, a series of three-point bending tests were performed on various single edge notched prisms (SENPs). Compression tests were also performed for quality assurance. Testing was conducted on specimens of different cross-sections, span/depth (S/D) ratios, curing regimes, ages, and fiber contents. By comparing the results from prisms of different sizes this study examines the weakening mechanism due to the size effect. Furthermore, by employing the concept of fracture energy it was possible to obtain a comparison of the fracture toughness and ductility. The model was determined based on a fit to P-w fracture curves, which was cross referenced for comparability to the results. Once obtained the model was then compared to the models proposed by the AFGC in the 2003 and to the ACI 544 model for conventional fiber reinforced concretes.

Selective matrix remodeling to increase interstitial flow in lymphedema

Lymphedema is a disease characterized by swelling resulting from the accumulation of fluid in the extracellular matrix (ECM) of the skin. In order to alleviate this swelling, the authors sought to selectively degrade certain hydrophilic molecules in the ECM called glycosaminoglycans (GAGs). GAGs are long unbranched sugar molecules present in the ECM that attract water to their numerous negative charges. The authors hypothesized that the density of GAGs would increase in lymphedema and inhibit fluid from leaving the tissue. An existing mouse tail model of experimental lymphedema that reproduced important features of the human condition was used to evaluate GAG content in swollen tissue. In this model, a surgical excision of tissue was made circumferentially around the tail that caused swelling distal to the wound site. Tissue distal to the wound site was analyzed via two assays; one that measured hyaluronan (an unsulfated GAG) and another that measured sulfated GAGs (including Dermatan Sulfate and Chondroitin Sulfate), at various timepoints post surgical intervention. Hyaluronan (HA) levels were significantly higher than control (tissues with no surgical intervention) by day 5 (p

Sirolimus ameliorates the enhanced expression of metalloproteinases in a rat model of autosomal dominant polycystic kidney disease

BACKGROUND: Remodelling of matrix and tubular basement membranes (TBM) is a characteristic of polycystic kidney disease. We hypothesized that matrix and TBM degradation by metalloproteinases (MMPs) could promote cyst formation. We therefore investigated the renal expression of MMPs in the Han:SPRD rat model of autosomal dominant polycystic kidney disease (ADPKD) and examined the effect of sirolimus treatment on MMPs. METHODS: 5-week-old male heterozygous (Cy/+) and wild-type normal (+/+) rats were treated with sirolimus (2 mg/kg/day) through drinking water for 3 months. RESULTS: The mRNA and protein levels of MMP-2 and MMP-14 were markedly increased in the kidneys of heterozygous Cy/+ animals compared to wild-type +/+ as shown by RT-PCR and Western blot analyses for MMP-2 and MMP-14, and by zymography for MMP-2. Strong MMP-2 expression was detected by immunoperoxidase staining in cystic epithelial cells that also displayed an altered, thickened TBM. Tissue inhibitor of metalloproteinases-2 (TIMP-2) expression was not changed in Cy/+ kidneys. Sirolimus treatment leads to decreased protein expression of MMP-2 and MMP-14 in Cy/+, whereas MMP-2 and MMP-14 mRNA levels and TIMP-2 protein levels were not affected by sirolimus. CONCLUSION: In summary, in kidneys of the Han:SPRD rat model of ADPKD, there is a marked upregulation of MMP-2 and MMP-14. Sirolimus treatment was associated with a marked improvement of MMP-2 and MMP-14 overexpression, and this correlated also with less matrix and TBM alterations and milder cystic disease.