Expression of Fas, FasL, caspase-8 and other factors of the extrinsic apoptotic pathway during the onset of interdigital tissue elimination.

Elimination of the interdigital web is considered to be the classical model for assessing apoptosis. So far, most of the molecules described in the process have been connected to the intrinsic (mitochondrial) pathway. The extrinsic (receptor mediated) apoptotic pathway has been rather neglected, although it is important in development, immunomodulation and cancer therapy. This work aimed to investigate factors of the extrinsic apoptotic machinery during interdigital regression with a focus on three crucial initiators: Fas, Fas ligand and caspase-8. Immunofluorescent analysis of mouse forelimb histological sections revealed abundant expression of these molecules prior to digit separation. Subsequent PCR Array analyses indicated the expression of several markers engaged in the extrinsic pathway. Between embryonic days 11 and 13, statistically significant increases in the expression of Fas and caspase-8 were observed, along with other molecules involved in the extrinsic apoptotic pathway such as Dapk1, Traf3, Tnsf12, Tnfrsf1A and Ripk1. These results demonstrate for the first time the presence of extrinsic apoptotic components in mouse limb development and indicate novel candidates in the molecular network accompanying the regression of interdigital tissue during digitalisation.

Cellular and molecular drivers of differential organ growth: insights from the limbs of Monodelphis domestica.

A fundamental question in biology is "how is growth differentially regulated during development to produce organs of particular sizes?" We used a new model system for the study of differential organ growth, the limbs of the opossum (Monodelphis domestica), to investigate the cellular and molecular basis of differential organ growth in mammals. Opossum forelimbs grow much faster than hindlimbs, making opossum limbs an exceptional system with which to study differential growth. We first used the great differences in opossum forelimb and hindlimb growth to identify cellular processes and molecular signals that underlie differential limb growth. We then used organ culture and pharmacological addition of FGF ligands and inhibitors to test the role of the Fgf/Mitogen-activated protein kinases (MAPK) signaling pathway in driving these cellular processes. We found that molecular signals from within the limb drive differences in cell proliferation that contribute to the differential growth of the forelimb and hindlimbs of opossums. We also found that alterations in the Fgf/MAPK pathway can generate differences in cell proliferation that mirror those observed between wild-type forelimb and hindlimbs of opossums and that manipulation of Fgf/MAPK signaling affects downstream focal adhesion-extracellular matrix (FA-ECM) and Wnt signaling in opossum limbs. Taken together, these findings suggest that evolutionary changes in the Fgf/MAPK pathway could help drive the observed differences in cell behaviors and growth in opossum forelimb and hindlimbs.

A Hoxa13:Cre mouse strain for conditional gene manipulation in developing limb, hindgut, and urogenital system.

The developing limb is a useful model for studying organogenesis and developmental processes. Although Cre alleles exist for conditional loss- or gain-of-function in limbs, Cre alleles targeting specific limb subdomains are desirable. Here we report on the generation of the Hoxa13:Cre line, in which the Cre gene is inserted in the endogenous Hoxa13 gene. We provide evidence that the Cre is active in embryonic tissues/regions where the endogenous Hoxa13 gene is expressed. Our results show that cells expressing Hoxa13 in developing limb buds contribute to the entire autopod (hand/feet) skeleton and validate Hoxa13 as a distal limb marker as far as the skeleton is concerned. In contrast, in the limb musculature, Cre-based fate mapping shows that almost all muscle masses of the zeugopod (forearm) and part of the triceps contain Hoxa13-expressing cells and/or their descendants. Besides the limb, the activity of the Cre is detectable in the urogenital system and the hindgut, primarily in the epithelium and smooth muscles. Together our data show that the Hoxa13:Cre allele is a useful tool for conditional gene manipulation in the urogenital system, posterior digestive tract, autopod and part of the limb musculature.

Limb bud colonization by somite-derived angioblasts is a crucial step for myoblast emigration.

We have combined the use of mouse genetic strains and the mouse-into-chicken chimera system to determine precisely the sequence of forelimb colonization by presomitic mesoderm (PSM)-derived myoblasts and angioblasts, and the possible role of this latter cell type in myoblast guidance. By creating a new Flk1/Pax3 double reporter mouse line, we have established the precise timetable for angioblast and myoblast delamination/migration from the somite to the limb bud. This timetable was conserved when mouse PSM was grafted into a chicken host, which further validates the experimental model. The use of Pax3(GFP/GFP) knockout mice showed that establishment of vascular endothelial and smooth muscle cells (SMCs) is not compromised by the absence of Pax3. Of note, Pax3(GFP/GFP) knockout mouse PSM-derived cells can contribute to aortic, but not to limb, SMCs that are derived from the somatopleure. Finally, using the Flk1(lacZ)(/)(lacZ) knockout mouse, we show that, in the absence of angioblast and vascular network formation, myoblasts are prevented from migrating into the limb. Taken together, our study establishes for the first time the time schedule for endothelial and skeletal muscle cell colonization in the mouse limb bud and establishes the absolute requirement of endothelial cells for myoblast delamination and migration to the limb. It also reveals that cells delaminating from the somites display marked differentiation traits, suggesting that if a common progenitor exists, its lifespan is extremely short and restricted to the somite.

Microstructural analysis of collagen and elastin fibres in the kangaroo articular cartilage reveals a structural divergence depending on its local mechanical environment.

OBJECTIVE: To assess the microstructure of the collagen and elastin fibres in articular cartilage under different natural mechanical loading conditions and determine the relationship between the microstructure of collagen and its mechanical environment. METHOD: Articular cartilage specimens were collected from the load bearing regions of the medial femoral condyle and the medial distal humerus of adult kangaroos. The microstructure of collagen and elastin fibres of these specimens was studied using laser scanning confocal microscopy (LSCM) and the orientation and texture features of the collagen were analysed using ImageJ. RESULTS: A zonal arrangement of collagen was found in kangaroo articular cartilage: the collagen fibres aligned parallel to the surface in the superficial zone and ran perpendicular in the deep zone. Compared with the distal humerus, the collagen in the femoral condyle was less isotropic and more clearly oriented, especially in the superficial and deep zones. The collagen in the femoral condyle was highly heterogeneous, less linear and more complex. Elastin fibres were found mainly in the superficial zone of the articular cartilage of both femoral condyle and distal humerus. CONCLUSIONS: The present study demonstrates that the collagen structure and texture of kangaroo articular cartilage is joint-dependent. This finding emphasizes the effects of loading on collagen development and suggests that articular cartilage with high biochemical and biomechanical qualities could be achieved by optimizing joint loading, which may benefit cartilage tissue engineering and prevention of joint injury. The existence of elastin fibres in articular cartilage could have important functional implications.

Cellular and molecular investigations into the development of the pectoral girdle.

The forelimbs of higher vertebrates are composed of two portions: the appendicular region (stylopod, zeugopod and autopod) and the less prominent proximal girdle elements (scapula and clavicle) that brace the limb to the main trunk axis. We show that the formation of the muscles of the proximal limb occurs through two distinct mechanisms. The more superficial girdle muscles (pectoral and latissimus dorsi) develop by the "In-Out" mechanism whereby migration of myogenic cells from the somites into the limb bud is followed by their extension from the proximal limb bud out onto the thorax. In contrast, the deeper girdle muscles (e.g. rhomboideus profundus and serratus anterior) are induced by the forelimb field which promotes myotomal extension directly from the somites. Tbx5 inactivation demonstrated its requirement for the development of all forelimb elements which include the skeletal elements, proximal and distal muscles as well as the sternum in mammals and the cleithrum of fish. Intriguingly, the formation of the diaphragm musculature is also dependent on the Tbx5 programme. These observations challenge our classical views of the boundary between limb and trunk tissues. We suggest that significant structures located in the body should be considered as components of the forelimb.

Sox9 Is Required for Nail-Bed Differentiation and Digit-Tip Regeneration.

The nail organ is a specialized appendage in which several ectodermal tissues coordinately function to sustain nail growth, a process that is coupled to digit regeneration. In this study, we show that the transcription factor Sox9 is expressed in several cell populations in the mouse digit tip. We found a SOX9+ cell population in the nail bed, and genetic lineage tracing showed that this is a transient cell population differentiated from matrix nail stem cells. In the absence of Sox9, nail matrix stem cells fail to differentiate into epithelial nail-bed cells and proliferate, thus expanding distally and following the corneocyte fate, which results in outlandishly large fingernails. In addition, the tip of the underlying terminal phalanx undergoes bone regression. Sox9-lineage tracing also revealed the existence of a continuous cell supply from a Sox9-expressing population residing in the basal layers to the entire hyponychium epidermis. Furthermore, digit-tip regeneration is compromised in Sox9-knockout mice, revealing an essential role for the gene during this process. These results will contribute to understand the cellular and molecular basis of mammalian nail organ homeostasis and disease and digit-tip regeneration and will help to design new treatment strategies for patients with nail diseases or amputation.

Immunoreactivities of AR, ERα, ERß and aromatase in the nuptial pad of Chinese brown frog (Rana dybowskii) during pre-hibernation and the breeding period.

There is a prominent local raised pad called nuptial pad on the forelimb of Chinese brown frog (Rana dybowskii), which is hypothetically concluded as an enhancement of the grip and a spreader of pheromone during the amplexus. In this study, we investigated the immunolocalization and protein expression levels of AR, ERα, ERß and aromatase in the nuptial pad of R. dybowskii during pre-hibernation and the breeding period. Histologically, the annual development of the nuptial pad in R. dybowskii is manifested as the larger area of specialized mucous gland and the longer length of papillary epidermal projection during the breeding period. AR, ERα, ERß and aromatase are present in the stratum granulosum, stratum spinosum, stratum basale and the secretory portion of specialized mucous glands during both periods. Western blotting results confirmed that AR, ERα and ERß protein levels are higher during pre-hibernation than those during the breeding season. These results suggest that nuptial pad is the direct target organ of androgen and estrogen. Androgen may participate in the regulation of annual development and glandular function of nuptial pad, and estrogen may play an endocrine, autocrine or paracrine role during pre-hibernation and the breeding period.

Conditional bursting enhances resonant firing in neocortical layer 2-3 pyramidal neurons.

The frequency response properties of neurons are critical for signal transmission and control of network oscillations. At subthreshold membrane potential, some neurons show resonance caused by voltage-gated channels. During action potential firing, resonance of the spike output may arise from subthreshold mechanisms and/or spike-dependent currents that cause afterhyperpolarizations (AHPs) and afterdepolarizations (ADPs). Layer 2-3 pyramidal neurons (L2-3 PNs) have a fast ADP that can trigger bursts. The present study investigated what stimuli elicit bursting in these cells and whether bursts transmit specific frequency components of the synaptic input, leading to resonance at particular frequencies. We found that two-spike bursts are triggered by step onsets, sine waves in two frequency bands, and noise. Using noise adjusted to elicit firing at approximately 10 Hz, we measured the gain for modulation of the time-varying firing rate as a function of stimulus frequency, finding a primary peak (7-16 Hz) and a high-frequency resonance (250-450 Hz). Gain was also measured separately for single and burst spikes. For a given spike rate, bursts provided higher gain at the primary peak and lower gain at intermediate frequencies, sharpening the high-frequency resonance. Suppression of bursting using automated current feedback weakened the primary and high-frequency resonances. The primary resonance was also influenced by the SK channel-mediated medium AHP (mAHP), because the SK blocker apamin reduced the sharpness of the primary peak. Our results suggest that resonance in L2-3 PNs depends on burst firing and the mAHP. Bursting enhances resonance in two distinct frequency bands.

The Role of Retinoic Acid in Establishing the Early Limb Bud.

Retinoic acid (RA) was one of the first molecules in the modern era of experimental embryology to be shown capable of generating profound effects on limb development. In this review, we focus on the earliest events of limb development and specifically on the role of RA in establishing the domain of cells that will go on to form the limb itself. Although there is some consensus on the role of RA during the earliest stages of limb formation, some controversy remains on the mechanism of RA action and the requirement for RA signaling in forming the hindlimb buds.

PI3K/Akt signaling as a key regulatory pathway for chondrocyte terminal differentiation.

Chondrogenesis is a well-coordinated multi-step differentiation process in which resting chondrocytes produce terminally differentiated hypertrophic chondrocytes through a proliferative stage. Here we show that phosphoinositide-3 kinase (PI3K) and its major downstream molecule, Akt, a serine-threonine kinase, play pivotal roles in this process. Akt signaling was activated in resting and proliferative chondrocytes but was reduced during terminal differentiation. We adopted two chondrocyte differentiation systems to investigate the roles of PI3K/Akt signaling in chondrogenesis. First, we employed an embryonic forelimb organ culture of transgenic mice expressing an Akt-Mer (a ligand-binding domain of a mutated estrogen receptor) fusion protein whose kinase activity was conditionally activated by treatment with 4-hydroxytamoxifen (4OHT). Activation of Akt signaling in embryonic chondrogenesis enhanced chondrocyte proliferation and inhibited hypertrophic differentiation, presumably due to the suppressed expression of Runx2, a transcription factor critical for chondrocyte terminal differentiation. Conversely, inhibition of PI3K by its inhibitor accelerated terminal hypertrophic differentiation, resulting in a shorter bone. Essentially the same results were obtained in a second line of experiments using human synovial stromal cells (hSSCs), which are mesenchymal progenitor cells isolated from adult joints. These findings demonstrate that PI3K/Akt signaling is a key regulator in terminal chondrocyte differentiation in both embryonic and adult chondrogenesis.

Do morphogenetic tissue rearrangements require active cell movements? The reversible inhibition of cell sorting and tissue spreading by cytochalasin B.

Previous studies have indicated that cell sorting and tissue spreading are caused by cell combination-specific differences in intercellular adhesive energies, acting in a system of motile cells. We wished to determine whether these adhesive energies could drive cell rearrangements as well as guide them. Accordingly, aggregates of intermixed embryonic cells were cultured in solutions of the drug cytochalasin B (CCB) at a concentration shown to inhibit the locomotion of cells on a solid surface. In addition, spherical aggregates of several kinds were cultured in mutual contact under similar conditions. Both cell sorting and tissue spreading were found to be inhibited. The prompt release of this inhibition upon removal of the CCB showed that the inhibited cells were not merely injured. Moreover, aggregation experiments showed that CCB did not prevent cells of several kinds from initiating mutual adhesions. In fact, heart cell aggregation was enhanced by CCB. We conclude that interfacial forces, originating outside the cell, act together with forces originating inside it in bringing about the morphogenetic movements of cell sorting and tissue spreading. We propose the term "cooperative cell locomotion" to describe translational movements of cells arising from such a combination of intrinsic and extrinsic forces.

Pectoral fins of Micropogonias furnieri: a histochemical and ultrastructural study.

The myotomal fibres of the pectoral fins of white croaker (Micropogonias furnieri) have been studied using histochemical techniques and transmission electron microscopy (TEM). Succinic dehydrogenase (SDH) for mitochondria, periodic acid Schiff (PAS) for glycogen, Sudan Black for lipids and myosin-adenosintriphosphatase (mATPase) pre-incubated at alkaline and acid pHs were used to visualize the contraction velocity. Three zones were determined: superficial (SZ), medium (MZ) and deep (DZ). Staining for SDH, PAS and Sudan Black was positive only in the SZ. The level of alkaline mATPase was the highest in fibres from the DZ, intermediate in the MZ and low in the SZ; at an acid pH, the reverse was obtained. Fibres from the SZ were small with large quantities of subsarcolemmal mitochondria, scarce intermyofibrilar mitochondria and a well-developed sarcoplasmic reticulum; the myofibrils displayed a polygonal distribution along the entire length of the fibre. Fibres in the MZ were larger than those in the SZ, the myofibrils were densely packed, mitochondria prevailed under the sarcolemma and the sarcoplasmic reticulum was not abundant. Fibres from the DZ were the largest, with ribbon-shaped myofibrils and scarce mitochondria. The intercellular space was abundant and nervous endings were frequently observed.

Hyaluronic acid, CD44 and RHAMM regulate myoblast behavior during embryogenesis.

Hyaluronic acid (HA) is an extracellular matrix (ECM) component that has been shown to play a significant role in regulating muscle cell behavior during repair and regeneration. For instance, ECM remodeling after muscle injury involves an upregulation in HA expression that is coupled with skeletal muscle precursor cell recruitment. However, little is known about the role of HA during skeletal muscle development. To gain insight into the way in which HA mediates embryonic myogenesis, we first determined the spatial distribution and gene expression of CD44, RHAMM and other HA related proteins in embryonic day (E)10.5 to E12.5 murine forelimbs. While HA and CD44 expression remained high, RHAMM decreased at both the protein (via immunohistochemistry) and RNA (via qPCR) levels. Next, we determined that 4-methylumbelliferone-mediated knockdown of HA synthesis inhibited the migration and proliferation of E11.5/E12.5 forelimb-derived cells. Then, the influence of CD44 and RHAMM on myoblast and connective tissue cell behavior was investigated using antibodies against these receptors. Anti-RHAMM, but not anti-CD44, significantly decreased the total distance myogenic progenitors migrated over 24 h, whereas both inhibited connective tissue cell migration. In contrast, anti-CD44 inhibited the proliferation of connective tissue cells and muscle progenitors, but anti-RHAMM had no effect. However, when myoblasts and connective tissue cells were depleted of CD44 and RHAMM by shRNA, motility and proliferation were significantly inhibited in both cells indicating that blocking cell surface-localized CD44 and RHAMM does not have as pronounced effect as global shRNA-mediated depletion of these receptors. These results show, for the first time, the distribution and activity of RHAMM in the context of skeletal muscle. Furthermore, our data indicate that HA, through interactions with CD44 and RHAMM, promotes myogenic progenitor migration and proliferation. Confirmation of the role of HA and its receptors in directing myogenesis will be useful for the design of regenerative therapies that aim to promote the restoration of damaged or diseased muscle.

Developmental 'awakening' of primary motor cortex to the sensory consequences of movement.

Before primary motor cortex (M1) develops its motor functions, it functions like a somatosensory area. Here, by recording from neurons in the forelimb representation of M1 in postnatal day (P) 8-12 rats, we demonstrate a rapid shift in its sensory responses. At P8-10, M1 neurons respond overwhelmingly to feedback from sleep-related twitches of the forelimb, but the same neurons do not respond to wake-related movements. By P12, M1 neurons suddenly respond to wake movements, a transition that results from opening the sensory gate in the external cuneate nucleus. Also at P12, fewer M1 neurons respond to individual twitches, but the full complement of twitch-related feedback observed at P8 is unmasked through local disinhibition. Finally, through P12, M1 sensory responses originate in the deep thalamorecipient layers, not primary somatosensory cortex. These findings demonstrate that M1 initially establishes a sensory framework upon which its later-emerging role in motor control is built.

Visualizing muscle cell migration in situ.

BACKGROUND: Cell migration has been studied extensively by manipulating and observing cells bathed in putative chemotactic or chemokinetic agents on planar substrates. This environment differs from that in vivo and, consequently, the cells can behave abnormally. Embryo slices provide an optically accessible system for studying cellular navigation pathways during development. We extended this system to observe the migration of muscle precursors from the somite into the forelimb, their cellular morphology, and the localization of green fluorescent protein (GFP)-tagged adhesion-related molecules under normal and perturbed conditions. RESULTS: Muscle precursors initiated migration synchronously and migrated in broad, rather than highly defined, regions. Bursts of directed migration were followed by periods of meandering or extension and retraction of cell protrusions. Although paxillin did not localize to discernible intracellular structures, we found that alpha-actinin localized to linear, punctate structures, and the alpha5 integrin to some focal complexes and/or vesicle-like concentrations. Alterations in the expression of adhesion molecules inhibited migration. The muscle precursors migrating in situ formed unusually large, long-lived protrusions that were polarized in the direction of migration. Unlike wild-type Rac, a constitutively active Rac localized continuously around the cell surface and promoted random protrusive activity and migration. CONCLUSIONS: The observation of cellular migration and the dynamics of molecular organization at high temporal and spatial resolution in situ is feasible. Migration from the somite to the wing bud is discontinuous and not highly stereotyped. In situ, local activation of Rac appears to produce large protrusions, which in turn, leads to directed migration. Adhesion can also regulate migration.

Pitx1 determines characteristic hindlimb morphologies in cartilage micromass culture.

The shapes of homologous skeletal elements in the vertebrate forelimb and hindlimb are distinct, with each element exquisitely adapted to their divergent functions. Many of the signals and signalling pathways responsible for patterning the developing limb bud are common to both forelimb and hindlimb. How disparate morphologies are generated from common signalling inputs during limb development remains poorly understood. We show that, similar to what has been shown in the chick, characteristic differences in mouse forelimb and hindlimb cartilage morphology are maintained when chondrogenesis proceeds in vitro away from the endogenous limb bud environment. Chondrogenic nodules that form in high-density micromass cultures derived from forelimb and hindlimb buds are consistently different in size and shape. We described analytical tools we have developed to quantify these differences in nodule morphology and demonstrate that characteristic hindlimb nodule morphology is lost in the absence of the hindlimb-restricted limb modifier gene Pitx1. Furthermore, we show that ectopic expression of Pitx1 in the forelimb is sufficient to generate nodule patterns characteristic of the hindlimb. We also demonstrate that hindlimb cells are less adhesive to the tissue culture substrate and, within the limb environment, to the extracellular matrix and to each other. These results reveal autonomously programmed differences in forelimb and hindlimb cartilage precursors of the limb skeleton are controlled, at least in part, by Pitx1 and suggest this has an important role in generating distinct limb-type morphologies. Our results demonstrate that the micromass culture system is ideally suited to study cues governing morphogenesis of limb skeletal elements in a simple and experimentally tractable in vitro system that reflects in vivo potential.

Survival of developing motor neurons mediated by Rho GTPase signaling pathway through Rho-kinase.

A variety of neurons generated during embryonic development survive or undergo programmed cell death (PCD) at later developmental stages. Survival or death of developing neurons is generally considered to depend on trophic support from various target tissues. The small GTPase Rho regulates diverse cellular processes such as cell morphology, cell adhesion, cell motility, and apoptosis. Rho-dependent serine-threonine protein kinase (Rho-kinase-ROK-ROCK), one of the effector proteins, transmits signals for some Rho-mediated processes. Here, we report the in vivo role of the Rho signaling pathway through Rho-kinase during development of motor neurons (MNs) in the spinal cord. We performed conditional expression of a dominant-negative form for RhoA (RhoA DN) or for Rho-kinase (Rho-K DN) in transgenic mice by using the Cre-loxP system to suppress the activity of these signaling molecules in developing MNs. Expression of RhoA DN reduced the number of MNs in the spinal cord because of increased apoptosis while preserving the gross patterning of motor axons. Expression of Rho-K DN produced developmental defects similar to those observed in RhoA DN expression. In addition, analysis of transgenic mice expressing Rho-K DN showed that the increased apoptosis of MNs was induced at the early embryonic stages before the initiation of PCD, and that MN death at the late embryonic stages corresponding to the period of PCD was moderately enhanced in the transgenic mice. These findings indicate that the Rho signaling pathway, primarily through Rho-kinase, plays a crucial role in survival of spinal MNs during embryogenesis, particularly at the early developmental stages.

Bone quality in swine composite tissue allografts: effects of combination immunotherapy.

BACKGROUND: Tacrolimus (FK506)/mycophenolate mofetil (MMF)/prednisone combination immunosuppression therapy has been found to effectively prevent composite tissue allograft (CTA) rejection with minimal toxicity in a preclinical porcine model. These findings have been reproduced in 24 human hands transplanted in 18 patients. In CTAs containing bone, adequate bone quality and healing are essential for long-term functional success. The purpose of this study was to determine the effect FK506/MMF/prednisone immunotherapy has on bone quality and healing. METHODS: Forelimb CTA-flaps were transplanted in nine pigs. Recipient animals received FK506/MMF/prednisone therapy for 3 months. Bone quality was studied pre- and posttransplant by measuring acoustic velocity and density and by calculating elastic coefficients. Additional bone quality analyses were performed on unoperated limbs, and in bone grafts from two pigs that had autograft procedures performed. Bone healing was assessed using radiographic analysis. RESULTS: Three animals were lost to immunosuppression-related complications before the endpoint of the study. The bone component of all six CTA-flaps showed normal healing. Although results of the bone density measurements were not significantly different when comparing pre- to posttransplant values, acoustic velocity and elastic coefficient measurements showed a significant decrease posttransplant indicating a decrease in bone quality. CONCLUSIONS: FK506/MMF/prednisone combination therapy prevented rejection, did not adversely affect bone quality, and showed normal bone healing. The transplant procedure itself decreased bone quality more than the immunosuppression regimen did over the observation period in this study. Based on these findings, we conclude to prevent CTA failure it is important to monitor bone quality posttransplant.

A framework for three-dimensional simulation of morphogenesis.

We present COMPUCELL3D, a software framework for three-dimensional simulation of morphogenesis in different organisms. COMPUCELL3D employs biologically relevant models for cell clustering, growth, and interaction with chemical fields. COMPUCELL3D uses design patterns for speed, efficient memory management, extensibility, and flexibility to allow an almost unlimited variety of simulations. We have verified COMPUCELL3D by building a model of growth and skeletal pattern formation in the avian (chicken) limb bud. Binaries and source code are available, along with documentation and input files for sample simulations, at http:// compucell.sourceforge.net.