Mechano-regulation of mesenchymal stem cell differentiation and collagen organisation during skeletal tissue repair

A number of mechano-regulation theories have been proposed that relate the differentiation pathway of mesenchymal stem cells (MSCs) to their local biomechanical environment. During spontaneous repair processes in skeletal tissues, the organisation of the extracellular matrix is a key determinant of its mechanical fitness. In this paper, we extend the mechano-regulation theory proposed by Prendergast et al. (J Biomech 30(6):539–548, 1997) to include the role of the mechanical environment on the collagen architecture in regenerating soft tissues. A large strain anisotropic poroelastic material model is used in a simulation of tissue differentiation in a fracture subject to cyclic bending (Cullinane et al. in J Orthop Res 20(3):579–586, 2002). The model predicts non-union with cartilage and fibrous tissue formation in the defect. Predicted collagen fibre angles, as determined by the principal decomposition of strain- and stress-type tensors, are similar to the architecture seen in native articular cartilage and neoarthroses induced by bending of mid-femoral defects in rats. Both stress and strain-based remodelling stimuli successfully predicted the general patterns of collagen fibre organisation observed in vivo. This provides further evidence that collagen organisation during tissue differentiation is determined by the mechanical environment. It is envisioned that such predictive models can play a key role in optimising MSC-based skeletal repair therapies where recapitulation of the normal tissue architecture is critical to successful repair.     

Bone tissue,lyophilized and stored at room temperature for 15 days or more,is not capable of transmitting HIV,HCV or HBV

Over the past 57 years, 17 recipients of frozen bone have been infected with: HIV (Centers for Disease Control and Prevention in Morb Mortal Wkly Rep MMWR 37(39):597–599, 1988; Li et al. in J Formos Med Assoc 100(5):350–351, 2001; Simonds et al. in NEJM 326(11):726–732, 1992; Schratt et al. in Unfallchirurg 99(9):679–684, 1996); HCV (Eggen and Nordbo in NEJM 326(6):411, 1992; Conrad et al. in J Bone Joint Surg Am 77:214–224, 1995; Trotter in J Bone Joint Surg Am 851(11):2215–2217, 2003; Tugwell et al. in Ann of Internal Med 143(9):648–654, 2005); or HBV (Shutkin in J Bone Joint Surg Am 36:160–162, 1954). However, bone, lyophilized and stored at room temperature, has never transmitted these viral diseases. A literature review was undertaken to determine whether there is any evidence that lyophilized bone is capable of transmitting HIV, HCV and HBV.     

Evolving mechanical properties of a model of abdominal aortic aneurysm

The novel three-dimensional (3D) mathematical model for the development of abdominal aortic aneurysm (AAA) of Watton et al. Biomech Model Mechanobiol 3(2): 98–113, (2004) describes how changes in the micro-structure of the arterial wall lead to the development of AAA, during which collagen remodels to compensate for loss of elastin. In this paper, we examine the influence of several of the model’s material and remodelling parameters on growth rates of the AAA and compare with clinical data. Furthermore, we calculate the dynamic properties of the AAA at different stages in its development and examine the evolution of clinically measurable mechanical properties. The model predicts that the maximum diameter of the aneurysm increases exponentially and that the ratio of systolic to diastolic diameter decreases from 1.13 to 1.02 as the aneurysm develops; these predictions are consistent with physiological observations of Vardulaki et al. Br J Surg 85:1674–1680 (1998) and Lanne et al. Eur J Vasc Surg 6:178–184 (1992), respectively. We conclude that mathematical models of aneurysm growth have the potential to be useful, noninvasive diagnostic tools and thus merit further development.     

Growth factor regulation of proliferation in primary cultures of small intestinal epithelium

Summary Although the intestinal epithelium is one of the most rapidly renewing tissues, little is known about the major growth factors that control the rate of cell replacement and migration. Recently, a primary culture model has been described for the developing rat small intestinal epithelium, which permits epithelial growth while maintaining interactions with associated stromal cells, thereby possessing several contextual advantages over established cell lines (Evans et al., 1992). We have used this model to begin to determine the factors that may be involved in controlling intestinal epithelial cell proliferation. Under the conditions examined, no single growth factor promoted exclusive proliferation of epithelial cells; stromal cell proliferation was also apparent. The most potent stimulators of epithelial proliferation were insulin and insulin-like growth factor 1 (IGF-1). These factors also appeared to inhibit migration of the epithelial cells. 5–10 ng/ml EGF, 5–20 ng/ml TGFα, and 10–20 ng/ml PDGF also slightly increased epithelial cell numbers. Cell proliferation was inhibited by 0.1 ng/ml TGFβ-1. In Dulbecco’s modified Eagle’s medium (DMEM) containing 0.25 IU/ml insulin, glucose levels of 2–3 g/liter permitted epithelial growth with limited expansion of the stromal cell population. Higher levels of glucose further stimulated the nonepithelial cell types. Transferrin was also a potent stimulator of both cell types.     

CCN1 contributes to skin connective tissue aging by inducing age-associated secretory phenotype in human skin dermal fibroblasts

Dermal connective tissue collagen is the major structural protein in skin. Fibroblasts within the dermis are largely responsible for collagen production and turnover. We have previously reported that dermal fibroblasts, in aged human skin in vivo, express elevated levels of CCN1, and that CCN1 negatively regulates collagen homeostasis by suppressing collagen synthesis and increasing collagen degradation (Quan et al. Am J Pathol 169:482–90, 2006, J Invest Dermatol 130:1697–706, 2010). In further investigations of CCN1 actions, we find that CCN1 alters collagen homeostasis by promoting expression of specific secreted proteins, which include matrix metalloproteinases and proinflammatory cytokines. We also find that CCN1-induced secretory proteins are elevated in aged human skin in vivo. We propose that CCN1 induces an “Age-Associated Secretory Phenotype”, in dermal fibroblasts, which mediates collagen reduction and fragmentation in aged human skin.     

ROS-Mediated ABA Signaling

In plants, reactive oxygen species (ROS) are short-lived molecules produced through various cellular mechanisms in response to biotic and abiotic stimuli. ROS function as second messengers for hormone signaling, development, oxygen deprivation, programmed cell death, and plant–pathogen interactions. Recent research on ROS-mediated responses has produced stimulating findings such as the specific sources of ROS production, molecular elements that work in ROS-mediated signaling and homeostasis, and a ROS-regulated gene network (Neill et al., Curr Opin Plant Biol 5:388–395, 2002a; Apel and Hirt, Annu Rev Plant Biol 55:373–399, 2004; Mittler et al., Trends Plant Sci 9:490–498, 2004; Mori and Schroeder, Plant Physiol 135:702–708, 2004; Kwak et al., Plant Physiol 141:323–329, 2006; Torres et al., Plant Physiol 141:373–378, 2006; Miller et al., Physiol Plant 133:481–489, 2008). In this review, we highlight new discoveries in ROS-mediated abscisic acid (ABA) signaling. Drs. Daeshik Cho and June M. Kwak are the corresponding authors for this paper.     

Development of a hemicornea from human primary cell cultures for pharmacotoxicology testing

We report the reconstruction and characterization of a hemicornea (epithelialized stroma), using primary human cells, for use in research and as an alternative to the use of animals in pharmacotoxicology testing. To create a stromal equivalent, keratocytes from human corneas were cultured in collagen–glycosaminoglycan–chitosan foams. Limbal stem cell-derived epithelial cells were seeded on top of these, giving rise to hemi-corneas. The epithelium appeared morphologically similar to its physiological counterpart, as shown by the basal cell expression of p63 isoforms including, in some cases, the stem cell marker p63ΔNα, and the expression of keratin 3 and 14-3-3σ in the upper cell layers. In addition, the cuboidal basal epithelial cells were anchored to a basement membrane containing collagen IV, laminin 5, and hemidesmosomes. In the stromal part, the keratocytes colonized the porous scaffold, formed a network of interconnecting cells, and synthesized an ultrastructurally organized extracellular matrix (ECM) containing collagen types I, V, and VI. Electron microscopy showed the newly synthesized collagen fibrils to have characteristic periodic striations, with diameters and interfibril spacings similar to those found in natural corneas. Compared to existing models for corneal pharmacotoxicology testing, this new model more closely approaches physiological conditions by including the inducing effects of mesenchyme and cell–matrix interactions on epithelial cell morphogenesis.     

Temperature controls a latitudinal gradient in the proportion of watershed nitrogen exported to coastal ecosystems

Increased export of biologically available nitrogen (N) to the coastal zone is strongly linked to eutrophication, which is a major problem in coastal marine ecosystems (NRC (2000) Clean Coastal Waters: Understanding and Reducing the Effects of Nutrient Pollution. National Academy Press, Washington, DC; Bricker et al. (1999) National Estuarine Eutrophication Assessment. Effects of nutrient enrichment in the nation’s estuaries. NOAA-NOS Special Projects Office, Silver Spring, MD). However, not all of the nitrogen input to a watershed is exported to the coast (Howarth et al. (1996) Biogeochemistry 35:75–139; Jordan and Weller (1996) Bioscience 46:655–664). Global estimates of nitrogen export to coasts have been taken to be 25% of watershed input, based largely on northeastern U.S. observations (Galloway et al. (2004) Biogeochemistry 70:153–226; Boyer et al. (2006) Global Biogeochem Cycle 20:Art. No. GB1S91). We applied the N budgeting methodology developed for the International SCOPE Nitrogen project (Howarth et al. (1996) Biogeochemistry 35:75–139; Boyer et al. (2002) Biogeochemistry 57:137–169) to 12 watersheds in the southeastern U.S., and compared them with estimates of N export for 16 watersheds in the northeastern U.S. (Boyer et al. (2002) Biogeochemistry 57:137–169). In southeastern watersheds, average N export was only 9% of input, suggesting the need for downward revision of global estimates. The difference between northern and southern watersheds is not a function of the absolute value of N inputs, which spanned a comparable range and were positively related to export in both cases. Rather, the proportion of N exported was significantly related to average watershed temperature (% N export = 58.41 e^{−0.11 * temperature}; R ² = 0.76), with lower proportionate nitrogen export in warmer watersheds. In addition, we identified a threshold in proportionate N export at 38°N latitude that corresponds to a reported breakpoint in the rate of denitrification at 10–12°C. We hypothesize that temperature, by regulating denitrification, results in increased proportionate N export at higher latitudes. Regardless of the mechanism, these observations suggest that temperature increases associated with future climate change may well reduce the amount of nitrogen that reaches estuaries, which will have implications for coastal eutrophication.     

Estimating material parameters of a structurally based constitutive relation for skin mechanics

This paper presents a structurally based modeling framework to characterize the structure–function relation in skin tissues, based upon biaxial tensile experiments performed in vitro on porcine skin. Equi-axial deformations were imposed by stretching circular skin specimens uniformly along twelve directions, and the resultant loads at the membrane attachment points were measured. Displacement fields at each deformation step were tracked using an image 2D cross-correlation technique. A modeling framework was developed to simulate the experiments, whereby measured forces were applied to finite element models that were created to represent the geometry and structure of the tissue samples. Parameters of a structurally based constitutive relation were then identified using nonlinear optimization. Results showed that the ground matrix stiffness ranged from 5 to 32 kPa, fiber orientation mean from 2 to 13^◦ from the torso midline, fiber undulation mean from 1.04 to 1.34 and collagen fiber stiffness from 48 to 366 MPa. It was concluded that the objective function was highly sensitive to the mean orientation and that a priori information about fiber orientation mean was important for the reliable identification of constitutive parameters.     

Keratan sulfate expression during avian craniofacial morphogenesis

Summary Using the monoclonal antibody MZ15 in immunocytochemical and ultrastructural studies we have been able to determine the spatiotemporal pattern of keratan sulfate (KS) distribution during quail craniofacial morphogenesis. KS-containing proteoglycans are found associated with invaginating placodes (olfactory, lens and otic), in developing pronephric tubules, notochord, pharynx and endocardium, and display developmental regulation. The appearance of such proteoglycans (PGs) during placode morphogenesis is particularly striking and we suggest that they may be an important component of the extracellular matrix which has been previously implicated in mediating the morphogenetic interactions and cell movements occurring at these sites. The otic vesicle during stage 18–22 displays a notable asymmetric distribution of KS-containing PGs. The role that these molecules may play and the reasons for this regionalization are, as yet, unclear but it is conceivable that the distribution of proteoglycans at this stage reflects subsequent differentiative events during otocyst development. Furthermore, our ultrastructural observations indicate that over the developmental period studied (H & H stages 8–22) keratan sulfate exists in at least two proteoglycan forms. Some spatiotemporal correlation has been found to exist between the distributions of KS-containing PGs and type II collagen as previously reported by Thorogood et al. (1986). We suggest that the proteoglycan detected at such sites is cartilage-specific proteoglycan and that it plays an important role, together with type II collagen, in the “signalling” mechanism which specifies the subsequent pattern of the chondrocranium. It is proposed that this interaction at epithelio-mesenchymal interfaces in the developing head parallels the matrix-mediated tissue interaction between notochord and somites which results in the formation of the cartilaginous primordia of the vertebrae from the sclerotomes as reported by Lash and Vasan (1978).     

Modeling interlamellar interactions in angle-ply biologic laminates for annulus fibrosus tissue engineering

Mechanical function of the annulus fibrosus of the intervertebral disc is dictated by the composition and microstructure of its highly ordered extracellular matrix. Recent work on engineered angle-ply laminates formed from mesenchymal stem cell (MSC)-seeded nanofibrous scaffolds indicates that the organization of collagen fibers into planes of alternating alignment may play an important role in annulus fibrosus tissue function. Specifically, these engineered tissues can resist tensile deformation through shearing of the interlamellar matrix as layers of collagen differentially reorient under load. In the present work, a hyperelastic constitutive model was developed to describe the role of interlamellar shearing in reinforcing the tensile response of biologic laminates, and was applied to experimental results from engineered annulus constructs formed from MSC-seeded nanofibrous scaffolds. By applying the constitutive model to uniaxial tensile stress–strain data for bilayers with three different fiber orientations, material parameters were generated that characterize the contributions of extrafibrillar matrix, fibers, and interlamellar shearing interactions. By 10 weeks of in vitro culture, interlamellar shearing accounted for nearly 50% of the total stress associated with uniaxial extension in the anatomic range of ply angle. The model successfully captured changes in function with extracellular matrix deposition through variations in the magnitude of model parameters with culture duration. This work illustrates the value of engineered tissues as tools to further our understanding of structure–function relations in native tissues and as a test-bed for the development of constitutive models to describe them.     

Surface-bound selectin–ligand binding is regulated by carrier diffusion

Two-dimensional (2D) kinetics of receptor–ligand interactions governs cell adhesion in many biological processes. While the dissociation kinetics of receptor–ligand bond is extensively investigated, the association kinetics has much less been quantified. Recently receptor–ligand interactions between two surfaces were investigated using a thermal fluctuation assay upon biomembrane force probe technique (Chen et al. in Biophys J 94:694–701, 2008). The regulating factors on association kinetics, however, are not well characterized. Here we developed an alternative thermal fluctuation assay using optical trap technique, which enables to visualize consecutive binding–unbinding transition and to quantify the impact of microbead diffusion on receptor–ligand binding. Three selectin constructs (sLs, sPs, and PLE) and their ligand P-selectin glycoprotein ligand 1 were used to conduct the measurements. It was indicated that bond formation was reduced by enhancing the diffusivity of selectin-coupled carrier, suggesting that carrier diffusion is crucial to determine receptor–ligand binding. It was also found that 2D forward rate predicted upon first-order kinetics was in the order of sPs > sLs > PLE and bond formation was history-dependent. These results further the understandings in regulating association kinetics of surface-bound receptor–ligand interactions.     

Simulation of planar soft tissues using a structural constitutive model: Finite element implementation and validation

Computational implementation of physical and physiologically realistic constitutive models is critical for numerical simulation of soft biological tissues in a variety of biomedical applications. It is well established that the highly nonlinear and anisotropic mechanical behaviors of soft tissues are an emergent behavior of the underlying tissue microstructure. In the present study, we have implemented a structural constitutive model into a finite element framework specialized for membrane tissues. We noted that starting with a single element subjected to uniaxial tension, the non-fibrous tissue matrix must be present to prevent unrealistic tissue deformations. Flexural simulations were used to set the non-fibrous matrix modulus because fibers have little effects on tissue deformation under three-point bending. Multiple deformation modes were simulated, including strip biaxial, planar biaxial with two attachment methods, and membrane inflation. Detailed comparisons with experimental data were undertaken to insure faithful simulations of both the macro-level stress–strain insights into adaptations of the fiber architecture under stress, such as fiber reorientation and fiber recruitment. Results indicated a high degree of fidelity and demonstrated interesting microstructural adaptions to stress and the important role of the underlying tissue matrix. Moreover, we apparently resolve a discrepancy in our 1997 study (Billiar and Sacks, 1997. J. Biomech. 30 (7), 753–756) where we observed that under strip biaxial stretch the simulated fiber splay responses were not in good agreement with the experimental results, suggesting non-affine deformations may have occurred. However, by correctly accounting for the isotropic phase of the measured fiber splay, good agreement was obtained. While not the final word, these simulations suggest that affine fiber kinematics for planar collagenous tissues is a reasonable assumption at the macro level. Simulation tools such as these are imperative in the design and simulation of native and engineered tissues.     

Use of weighted reference panels based on empirical estimates of ancestry for capturing untyped variation

Many association methods use a subset of genotyped single nucleotide polymorphisms (SNPs) to capture or infer genotypes at other untyped SNPs. We and others previously showed that tag SNPs selected to capture common variation using data from The International HapMap Consortium (Nature 437:1299–1320, 2005), The International HapMap Consortium (Nature 449:851–861, 2007) could also capture variation in populations of similar ancestry to HapMap reference populations (de Bakker et al. in Nat Genet 38:1298–1303, 2006; González-Neira et al. in Genome Res 16:323–330, 2006; Montpetit et al. in PLoS Genet 2:282–290, 2006; Mueller et al. in Am J Hum Genet 76:387–398, 2005). To capture variation in admixed populations or populations less similar to HapMap panels, a “cosmopolitan approach,” in which all samples from HapMap are used as a single reference panel, was proposed. Here we refine this suggestion and show that use of a “weighted reference panel,” constructed based on empirical estimates of ancestry in the target population (relative to available reference panels), is more efficient than the cosmopolitan approach. Weighted reference panels capture, on average, only slightly fewer common variants (minor allele frequency > 5%) than the cosmopolitan approach (mean r ² = 0.977 vs. 0.989, 94.5% variation captured vs. 96.8% at r ² > 0.8), across the five populations of the Multiethnic Cohort, but entail approximately 25% fewer tag SNPs per panel (average 538 vs. 718). These results extend a recent study in two Indian populations (Pemberton et al. in Ann Hum Genet 72:535–546, 2008). Weighted reference panels are potentially useful for both the selection of tag SNPs in diverse populations and perhaps in the design of reference panels for imputation of untyped genotypes in genome-wide association studies in admixed populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome-wide meta-analysis for rheumatoid arthritis

Meta-analysis is being increasingly used as a tool for integrating data from different studies of complex phenotypes, because the power of any one study to identify causal loci is limited. We applied a novel meta-analytical approach (Loesgen et al. in Genet Epidemiol 21(Suppl 1):S142–S147, 2001) in compiling results from four studies of rheumatoid arthritis in Caucasians including two studies from NARAC (Jawaheer et al. in Am J Hum Genet 68:927–936, 2001; Jawaheer et al. in Arthritis Rheum 48:906–916, 2003), one study from the UK (MacKay et al. in Arthritis Rheum 46:632–639, 2001) and one from France (Cornelis et al. in Proc Natl Acad Sci USA 95:10746–10750, 1998). For each study, we obtained NPL scores by performing interval mapping (2 cM intervals) using GeneHunter2 (Kruglyak et al. in Am J Hum Genet 58:1347–1363, 1996; Markianos et al. in Am J Hum Genet 68:963–977, 2001). The marker maps differed among the three consortium groups, therefore, the marker maps were aligned after the interval mapping was completed and the NPL scores that were within 1 cM of each other were combined using the method of Loesgen et al. (Genet Epidemiol 21(Suppl 1):S142–S147, 2001) by calculating the weighted average of the NPL score. This approach avoids some problems in analysis encountered by using GeneHunter2 when some markers in the sample are not genotyped. This procedure provided marginal evidence (P<0.05) of linkage on chromosome 1, 2, 5 and 18, strong evidence (P<0.01) on chromosomes 8 and 16, and overwhelming evidence in the HLA region of chromosome 6.     

The Morphostatic Limit for a Model of Skeletal Pattern Formation in the Vertebrate Limb

A recently proposed mathematical model of a “core” set of cellular and molecular interactions present in the developing vertebrate limb was shown to exhibit pattern-forming instabilities and limb skeleton-like patterns under certain restrictive conditions, suggesting that it may authentically represent the underlying embryonic process (Hentschel et al., Proc. R. Soc. B 271, 1713–1722, 2004). The model, an eight-equation system of partial differential equations, incorporates the behavior of mesenchymal cells as “reactors,” both participating in the generation of morphogen patterns and changing their state and position in response to them. The full system, which has smooth solutions that exist globally in time, is nonetheless highly complex and difficult to handle analytically or numerically. According to a recent classification of developmental mechanisms (Salazar-Ciudad et al., Development 130, 2027–2037, 2003), the limb model of Hentschel et al. is “morphodynamic,” since differentiation of new cell types occurs simultaneously with cell rearrangement. This contrasts with “morphostatic” mechanisms, in which cell identity becomes established independently of cell rearrangement. Under the hypothesis that development of some vertebrate limbs employs the core mechanism in a morphostatic fashion, we derive in an analytically rigorous fashion a pair of equations representing the spatiotemporal evolution of the morphogen fields under the assumption that cell differentiation relaxes faster than the evolution of the overall cell density (i.e., the morphostatic limit of the full system). This simple reaction–diffusion system is unique in having been derived analytically from a substantially more complex system involving multiple morphogens, extracellular matrix deposition, haptotaxis, and cell translocation. We identify regions in the parameter space of the reduced system where Turing-type pattern formation is possible, which we refer to as its “Turing space.” Obtained values of the parameters are used in numerical simulations of the reduced system, using a new Galerkin finite element method, in tissue domains with nonstandard geometry. The reduced system exhibits patterns of spots and stripes like those seen in developing limbs, indicating its potential utility in hybrid continuum-discrete stochastic modeling of limb development. Lastly, we discuss the possible role in limb evolution of selection for increasingly morphostatic developmental mechanisms.     

The Repair Response to Osteochondral Implant Types in a Rabbit Model

Current treatments for damaged articular cartilage (i.e., shaving the articular surface, perforation or abrasion of the subchondral bone, and resurfacing with periosteal and perichondrial resurfacing) often produce fibrocartilage, or hyaline-appearing repair that is not sustained over time (Henche 1967, Ligament and Articular Cartilage Injuries. Springer-Verlag, New York, NY, pp. 157–164; Insall 1974, Clin. Orthop. 101: 61–67; Mitchell and Shepard 1976, J. Bone Joint Surg. [Am.] 58: 230–233; O’Driscoll et al. 1986, J. Bone Joint Surg. [Am.] 68: 1017–1035; 1989, Trans. Orthop. Res. Soc. 14: 145; Kim et al. 1991, J. Bone Joint Surg. [Am.] 73: 1301–1315). Autologous chondrocyte transplantation, although promising, requires two surgeries, has site-dependent and patient age limitations, and has unknown long-term donor site morbidity (Brittberg et al. 1994, N Engl. J. Med. 331: 889–895; Minas 2003, Orthopedics 26: 945–947; Peterson et al. 2003, J. Bone Joint Surg. Am. 85-A(Suppl. 2): S17–S24). Osteochondral allografts remain a widely used method of articular resurfacing to delay arthritic progression. The present study compared the histological response to four types of osteochondral implants in a rabbit model: autograft, frozen, freeze-dried, and fresh implants. Specimens implanted in the femoral groove were harvested at 6 and 12 weeks. Results showed similar restoration of the joint surface regardless of implant type, with a trend toward better repair at the later timepoint. As has been observed in other studies (Frenkel et al. 1997, J. Bone Joint Surg. 79B: 281–286; Toolan et al. 1998, J. Biomed. Mater. Res. 41: 244–250), each group in this study had at least one specimen in which a healthy-appearing surface on the implant was not well-integrated with host tissues. Although the differences were not statistically significant, freeze-dried implants at both timepoints had the best histological scores. The osteochondral grafts tested successfully restored the gross joint surface and congruity. At 12 weeks, no significant differences were observed between the various allografts and autologous osteochondral grafts.     

Development of a biologically-based controlled growth and differentiation model for developmental toxicology

 A mathematical model is developed with a highly controlled birth and death process for precursor cells. This model is both biologically- and statistically-based. The controlled growth and differentiation (CGD) model limits the number of replications allowed in the development of a tissue or organ and thus, more closely reflects the presence of a true stem cell population. Leroux et al. (1996) presented a biologically-based dose-response model for developmental toxicology that was derived from a partial differential equation for the generating function. This formulation limits further expansion into more realistic models of mammalian development. The same formulae for the probability of a defect (a system of ordinary differential equations) can be derived through the Kolmogorov forward equations due to the nature of this Markov process. This modified approach is easily amenable to the expansion of more complicated models of the developmental process such as the one presented here. Comparisons between the Leroux et al. (1996) model and the controlled growth and differentiation (CGD) model as developed in this paper are also discussed. Received: 8 June 2001 / Revised version: 15 June 2002 / Published online: 26 September 2002 Keywords or phrases: Teratology – Multistate process – Cellular kinetics – Numerical simulation