Correlating flow induced shear stress and chondrocytes activity in micro-porous scaffold using computational fluid dynamics and rapid prototyping

Flow induced shear stress plays an important role in regulating cell growth and distribution in scaffolds. This study sought to correlate wall shear stress and chondrocytes activity for engineering design of micro-porous osteochondral grafts based on the hypothesis that it is possible to capture and discriminate between the transmitted force and cell response at the inner irregularities. Unlike common tissue engineering therapies with perfusion bioreactors in which flow-mediated stress is the controlling parameter, this work assigned the associated stress as a function of porosity to influence in vitro proliferation of chondrocytes. D-optimality criterion was used to accommodate three pore characteristics for appraisal in a mixed level fractional design of experiment (DOE); namely, pore size (4 levels), distribution pattern (2 levels) and density (3 levels). Micro-porous scaffolds (n=12) were fabricated according to the DOE using rapid prototyping of an acrylic-based bio-photopolymer. Computational fluid dynamics (CFD) models were created correspondingly and used on an idealized boundary condition with a Newtonian fluid domain to simulate the dynamic microenvironment inside the pores. In vitro condition was reproduced for the 3D printed constructs seeded by high pellet densities of human chondrocytes and cultured for 72 hours. The results showed that cell proliferation was significantly different in the constructs (p<0.05). Inlet fluid velocity of 3×10-2mms-1 and average shear stress of 5.65×10-2 Pa corresponded with increased cell proliferation for scaffolds with smaller pores in hexagonal pattern and lower densities. Although the analytical solution of a Poiseuille flow inside the pores was found insufficient for the description of the flow profile probably due to the outside flow induced turbulence, it showed that the shear stress would increase with cell growth and decrease with pore size. This correlation demonstrated the basis for determining the relation between the induced stress and chondrocyte activity to optimize microfabrication of engineered cartilaginous constructs.

A novel method using intranasal delivery of EdU demonstrates that accessory olfactory ensheathing cells respond to injury by proliferation

Olfactory ensheathing cells (OECs) play an important role in the continuous regeneration of the primary olfactory nervous system throughout life and for regeneration of olfactory neurons after injury. While it is known that several individual OEC subpopulations with distinct properties exist in different anatomical locations, it remains unclear how these different subpopulations respond to a major injury. We have examined the proliferation of OECs from one distinct location, the peripheral accessory olfactory nervous system, following large-scale injury (bulbectomy) in mice. We used crosses of two transgenic reporter mouse lines, S100ß-DsRed and OMP-ZsGreen, to visualise OECs, and main/accessory olfactory neurons, respectively. We surgically removed one olfactory bulb including the accessory olfactory bulb to induce degeneration, and found that accessory OECs in the nerve bundles that terminate in the accessory olfactory bulb responded by increased proliferation with a peak occurring 2 days after the injury. To label proliferating cells we used the thymidine analogue ethynyl deoxyuridine (EdU) using intranasal delivery instead of intraperitoneal injection. We compared and quantified the number of proliferating cells at different regions at one and four days after EdU labelling by the two different methods and found that intranasal delivery method was as effective as intrapeitoneal injection. We demonstrated that accessory OECs actively respond to widespread degeneration of accessory olfactory axons by proliferating. These results have important implications for selecting the source of OECs for neural regeneration therapies and show that intranasal delivery of EdU is an efficient and reliable method for assessing proliferation of olfactory glia.

Crucial roles for olfactory ensheathing cells and olfactory mucosal cells in the repair of damaged neural tracts

Olfactory ensheathing cells, the glial cells of the olfactory nervous system, exhibit unique growth-promoting and migratory properties that make them interesting candidates for cell therapies targeting neuronal injuries such as spinal cord injury. Transplantation of olfactory cells is feasible and safe in humans; however, functional outcomes are highly variable with some studies showing dramatic improvements and some no improvements at all. We propose that the reason for this is that the identity and purity of the cells is different in each individual study. We have shown that olfactory ensheathing cells are not a uniform cell population and that individual subpopulations of OECs are present in different regions of the olfactory nervous system, with strikingly different behaviors. Furthermore, the presence of fibroblasts and other cell types in the transplant can dramatically alter the behavior of the transplanted glial cells. Thus, a thorough characterization of the differences between olfactory ensheathing cell subpopulations and how the behavior of these cells is affected by the presence of other cell types is highly warranted.

Spatial and developmental heterogeneity of calcium signaling in olfactory ensheathing cells

Olfactory ensheathing cells (OECs) are specialized glial cells in the mammalian olfactory system supporting growth of axons from the olfactory epithelium into the olfactory bulb. OECs in the olfactory bulb can be subdivided into OECs of the outer nerve layer and the inner nerve layer according to the expression of marker proteins and their location in the nerve layer. In the present study, we have used confocal calcium imaging of OECs in acute mouse brain slices and olfactory bulbs in toto to investigate physiological differences between OEC subpopulations. OECs in the outer nerve layer, but not the inner nerve layer, responded to glutamate, ATP, serotonin, dopamine, carbachol, and phenylephrine with increases in the cytosolic calcium concentration. The calcium responses consisted of a transient and a tonic component, the latter being mediated by store-operated calcium entry. Calcium measurements in OECs during the first three postnatal weeks revealed a downregulation of mGluR(1) and P2Y(1) receptor-mediated calcium signaling within the first 2 weeks, suggesting that the expression of these receptors is developmentally controlled. In addition, electrical stimulation of sensory axons evoked calcium signaling via mGluR(1) and P2Y(1) only in outer nerve layer OECs. Downregulation of the receptor-mediated calcium responses in postnatal animals is reflected by a decrease in amplitude of stimulation-evoked calcium transients in OECs from postnatal days 3 to 21. In summary, the results presented reveal striking differences in receptor responses during development and in axon-OEC communication between the two subpopulations of OECs in the olfactory bulb.

Mesenchymal stromal cells and the repair of cartilage tissue

Articular cartilage has a limited intrinsic repair capacity, and thus defects are more likely to further degrade rather than undergo spontaneous self-repair. Whilst a number of surgical techniques have been developed to repair cartilage defects, their efficacy is generally poor and total joint replacement remains the gold standard, albeit last resort, treatment option. Cell-based therapies hold the greatest promise, as they appear uniquely capable of generating de novo cartilage tissue. Two approved therapies (ACI and MACI) are based on the premise that the transplantation of ex vivo expanded autologous chondrocyte populations, harvested from a non-load bearing region of the same joint, could be utilized to effectively regenerate cartilage tissue in the primary defect site. These therapeutic strategies are partially limited by our inability to harvest and expand adequate numbers of autologous chondrocytes that retain the appropriate phenotype. By contrast, the harvest and expansion of large numbers of mesenchymal stem/stromal cells (MSC) derived from tissues such as bone marrow and adipose is comparatively straightforward and has become routine in laboratories worldwide. Additionally, our understanding of the biochemical and biophysical signals required to drive the chondrogenic differentiation of MSC is rapidly increasing. It is conceivable that in the near future MSC expansion and differentiation technologies will offer a means to generate sufficient cell numbers, of an appropriate phenotype, for use in cartilage defect repair. In this chapter we review the relative potential of MSC and their likely contribution to cartilage regeneration.

Multiphasic construct studied in an ectopic osteochondral defect model

In vivo osteochondral defect models predominantly consist of small animals, such as rabbits. Although they have an advantage of low cost and manageability, their joints are smaller and more easily healed compared with larger animals or humans. We hypothesized that osteochondral cores from large animals can be implanted subcutaneously in rats to create an ectopic osteochondral defect model for routine and high-throughput screening of multiphasic scaffold designs and/or tissue-engineered constructs (TECs). Bovine osteochondral plugs with 4 mm diameter osteochondral defect were fitted with novel multiphasic osteochondral grafts composed of chondrocyte-seeded alginate gels and osteoblast-seeded polycaprolactone scaffolds, prior to being implanted in rats subcutaneously with bone morphogenic protein-7. After 12 weeks of in vivo implantation, histological and micro-computed tomography analyses demonstrated that TECs are susceptible to mineralization. Additionally, there was limited bone formation in the scaffold. These results suggest that the current model requires optimization to facilitate robust bone regeneration and vascular infiltration into the defect site. Taken together, this study provides a proof-of-concept for a high-throughput osteochondral defect model. With further optimization, the presented hybrid in vivo model may address the growing need for a cost-effective way to screen osteochondral repair strategies before moving to large animal preclinical trials.

Bone sialoprotein supports breast cancer cell adhesion proliferation and migration through differential usage of the αvβ3 and αvβ5 integrins

Bone sialoprotein (BSP), a secreted glycoprotein found in bone matrix, has been implicated in the formation of mammary microcalcifications and osteotropic metastasis of human breast cancer (HBC). BSP possesses an integrin-binding RGD (Arg-Gly-Asp) domain, which may promote interactions between HBC cells and bone extracellular matrix. Purified BSP, recombinant human BSP fragments and BSP-derived RGD peptides are shown to elicit migratory, adhesive, and proliferative responses in the MDA-MB-231 HBC cell line. Recombinant BSP fragment analysis localized a significant component of these activities to the RGD domain of the protein, and synthetic RGD peptides with BSP flanking sequences (BSPRGD) also conferred these responses. The fibronectin-derived RGD counterpart, GRGDSP (Gly-Arg-Gly-Asp-Ser-Pro), could not support these cellular responses, emphasizing specificity of the BSP configuration. Although most of the proliferative and adhesive responses could be attributed to RGD interactions, these interactions were only partly responsible for the migrational responses. Experiments with integrin-blocking antibodies demonstrated that BSP-RGD-induced migration utilizes the αvβ3 vitronectin receptor, whereas adhesion and proliferation responses were αvβ5-mediated. Using fluorescence activated cell sorting, we selected two separate subpopulations of MDA-MB-231 cells enriched for αvβ3 or αvβ5 respectively. Although some expression of the alternate αv integrin was still retained, the αvβ5-enriched MDA-MB-231 cells showed enhanced proliferative and adhesive responses, whereas the αvβ3-enriched subpopulation was suppressed for proliferation and adhesion, but showed enhanced migratory responses to BSP-RGD. In addition, similar analysis of two other HBC cell lines showed less marked, but similar RGD-dependent trends in adhesion and proliferation to the BSP fragments. Collectively, these data demonstrate BSP effects on proliferative, migratory, and adhesive functions in HBC cells and that the RGD-mediated component differentially employs αvβ3 and αvβ5 integrin receptors.

Laminin adhesion-selected primary human colon-cancer cells are more tumorigenic than the parenteral and nonadherent cells

Laminin has been shown to promote the malignant phenotype and the level of the 32/67 Kd laminin receptor has been found to correlate with Dukes' staging of colon cancer. A biopsy of a Dukes' stage B2 human colon carcinoma formed a tumor in a nude mouse after coinjection with Matrigel. The parental tumor and the murine tumor appeared identical at the histological level. A cell line LCC-C1 was established from the murine tumor. The cell line appeared moderately differentiated although it did not produce mucin in vitro; however, the xenograft in vivo did produce low levels of mucin. Laminin adherent and non-adherent cell lines were selected. The parental and the laminin-selected cell subclones adhered equally well to plastic and to fibronectin and showed similar growth rates on plastic. When injected subcutaneously into nude mice, the laminin-adherent cells formed relatively undifferentiated tumors that were twice as large as the parental cell tumors whereas the laminin non adherent cells formed very small, but highly differentiated tumors. These data demonstrate that subpopulations of tumor cells which differ in their tumorigenic properties can be selected based on their adhesion to laminin and thus provide models for studying the mechanisms of tumor growth.

Role of dentin matrix protein 1 in cartilage redifferentiation and osteoarthritis

Objective The aim of this study was to test the possible involvement, relevance and significance of dentin matrix protein 1 (DMP1) in chondrocyte redifferentiation and OA. Methods To examine the function of DMP1 in vitro, bone marrow stromal cells (BMSCs) and articular chondrocytes (ACs) were isolated and differentiated in micromasses in the presence or absence of DMP1 small interfering RNA and analysed for chondrogenic phenotype. The association of DMP1 expression with OA progression was analysed time dependently in the OA menisectomy rat model and in grade-specific OA human samples. Results It was found that DMP1 was strongly related to chondrogenesis, which was evidenced by the strong expression of DMP1 in the 14.5-day mouse embryonic cartilage development stage and in femoral heads of post-natal days 0 and 4. In vitro chondrogenesis in BMSCs and ACs was accompanied by a gradual increase in DMP1 expression at both the gene and protein levels. In addition, knockdown of DMP1 expression led to decreased chondrocyte marker genes, such as COL2A1, ACAN and SOX9, and an increase in the expression of COL10A and MMP13 in ACs. Moreover, treatment with IL-1β, a well-known catabolic culprit of proteoglycan matrix loss, significantly reduced the expression of DMP1. Furthermore, we also observed the suppression of DMP1 protein in a grade-specific manner in knee joint samples from patients with OA. In the menisectomy-induced OA model, an increase in the Mankin score was accompanied by the gradual loss of DMP1 expression. Conclusion Observations from this study suggest that DMP1 may play an important role in maintaining the chondrogenic phenotype and its possible involvement in altered cartilage matrix remodelling and degradation in disease conditions like OA.

Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.

Collective motion of dimers

We consider a discrete agent-based model on a one-dimensional lattice and a two-dimensional square lattice, where each agent is a dimer occupying two sites. Agents move by vacating one occupied site in favor of a nearest-neighbor site and obey either a strict simple exclusion rule or a weaker constraint that permits partial overlaps between dimers. Using indicator variables and careful probability arguments, a discrete-time master equation for these processes is derived systematically within a mean-field approximation. In the continuum limit, nonlinear diffusion equations that describe the average agent occupancy of the dimer population are obtained. In addition, we show that multiple species of interacting subpopulations give rise to advection-diffusion equations. Averaged discrete simulation data compares very well with the solution to the continuum partial differential equation models. Since many cell types are elongated rather than circular, this work offers insight into population-level behavior of collective cellular motion.

Building macroscale models from microscale probabilistic models : a general probabilistic approach for nonlinear diffusion and multispecies phenomena

A discrete agent-based model on a periodic lattice of arbitrary dimension is considered. Agents move to nearest-neighbor sites by a motility mechanism accounting for general interactions, which may include volume exclusion. The partial differential equation describing the average occupancy of the agent population is derived systematically. A diffusion equation arises for all types of interactions and is nonlinear except for the simplest interactions. In addition, multiple species of interacting subpopulations give rise to an advection-diffusion equation for each subpopulation. This work extends and generalizes previous specific results, providing a construction method for determining the transport coefficients in terms of a single conditional transition probability, which depends on the occupancy of sites in an influence region. These coefficients characterize the diffusion of agents in a crowded environment in biological and physical processes.

The glycosphingolipid P1 is an ovarian cancer-associated carbohydrate antigen involved in migration

Background The level of plasma-derived naturally circulating anti-glycan antibodies (AGA) to P1 trisaccharide has previously been shown to significantly discriminate between ovarian cancer patients and healthy women. Here we aim to identify the Ig class that causes this discrimination, to identify on cancer cells the corresponding P1 antigen recognised by circulating anti-P1 antibodies and to shed light into the possible function of this glycosphingolipid. Method An independent Australian cohort was assessed for the presence of anti-P1 IgG and IgM class antibodies using suspension array. Monoclonal and human derived anti-glycan antibodies were verified using three independent glycan-based immunoassays and flow cytometry-based inhibition assay. The P1 antigen was detected by LC-MS/MS and flow cytometry. FACS-sorted cell lines were studied on the cellular migration by colorimetric assay and real-time measurement using xCELLigence system. Results Here we show in a second independent cohort (n=155) that the discrimination of cancer patients is mediated by the IgM class of anti-P1 antibodies (P=0.0002). The presence of corresponding antigen P1 and structurally related epitopes in fresh tissue specimens and cultured cancer cells is demonstrated. We further link the antibody and antigen (P1) by showing that human naturally circulating and affinity-purified anti-P1 IgM isolated from patients ascites can bind to naturally expressed P1 on the cell surface of ovarian cancer cells. Cell-sorted IGROV1 was used to obtain two study subpopulations (P1-high, 66.1%; and P1-low, 33.3%) and observed that cells expressing high P1-levels migrate significantly faster than those with low P1-levels. Conclusions This is the first report showing that P1 antigen, known to be expressed on erythrocytes only, is also present on ovarian cancer cells. This suggests that P1 is a novel tumour-associated carbohydrate antigen recognised by the immune system in patients and may have a role in cell migration. The clinical value of our data may be both diagnostic and prognostic; patients with low anti-P1 IgM antibodies present with a more aggressive phenotype and earlier relapse.

Choosing an appropriate modelling framework for analysing multispecies co-culture cell biology experiments

In vitro cell biology assays play a crucial role in informing our understanding of the migratory, proliferative and invasive properties of many cell types in different biological contexts. While mono-culture assays involve the study of a population of cells composed of a single cell type, co-culture assays study a population of cells composed of multiple cell types (or subpopulations of cells). Such co-culture assays can provide more realistic insights into many biological processes including tissue repair, tissue regeneration and malignant spreading. Typically, system parameters, such as motility and proliferation rates, are estimated by calibrating a mathematical or computational model to the observed experimental data. However, parameter estimates can be highly sensitive to the choice of model and modelling framework. This observation motivates us to consider the fundamental question of how we can best choose a model to facilitate accurate parameter estimation for a particular assay. In this work we describe three mathematical models of mono-culture and co-culture assays that include different levels of spatial detail. We study various spatial summary statistics to explore if they can be used to distinguish between the suitability of each model over a range of parameter space. Our results for mono-culture experiments are promising, in that we suggest two spatial statistics that can be used to direct model choice. However, co-culture experiments are far more challenging: we show that these same spatial statistics which provide useful insight into mono-culture systems are insuffcient for co-culture systems. Therefore, we conclude that great care ought to be exercised when estimating the parameters of co-culture assays.

Making robust decisions for conservation with restricted money and knowledge

1. In conservation decision-making, we operate within the confines of limited funding. Furthermore, we often assume particular relationships between management impact and our investment in management. The structure of these relationships, however, is rarely known with certainty - there is model uncertainty. We investigate how these two fundamentally limiting factors in conservation management, money and knowledge, impact optimal decision-making. 2. We use information-gap decision theory to find strategies for maximizing the number of extant subpopulations of a threatened species that are most immune to failure due to model uncertainty. We thus find a robust framework for exploring optimal decision-making. 3. The performance of every strategy decreases as model uncertainty increases. 4. The strategy most robust to model uncertainty depends not only on what performance is perceived to be acceptable but also on available funding and the time horizon over which extinction is considered. 5. Synthesis and applications. We investigate the impact of model uncertainty on robust decision-making in conservation and how this is affected by available conservation funding. We show that subpopulation triage can be a natural consequence of robust decision-making. We highlight the need for managers to consider triage not as merely giving up, but as a tool for ensuring species persistence in light of the urgency of most conservation requirements, uncertainty and the poor state of conservation funding. We illustrate this theory by a specific application to allocation of funding to reduce poaching impact on the Sumatran tiger Panthera tigris sumatrae in Kerinci Seblat National Park. © 2008 The Authors.