Independent and combined effects of obesity and traumatic joint injury to the structure and composition of rat knee cartilage.

Osteoarthritis (OA) is a multifactorial joint disease characterized by articular cartilage degradation. Risk factors for OA include joint trauma, obesity, and inflammation, each of which can affect joint health independently, but their interaction and the associated consequences of such interaction were largely unexplored. Here, we studied compositional and structural alterations in knee joint cartilages of Sprague-Dawley rats exposed to two OA risk factors: joint injury and diet-induced obesity. Joint injury was imposed by surgical transection of anterior cruciate ligaments (ACLx), and obesity was induced by a high fat/high sucrose diet. Depth-dependent proteoglycan (PG) content and collagen structural network of cartilage were measured from histological sections collected previously in Collins et al.. (2015). We found that ACLx primarily affected the superficial cartilages. Compositionally, ACLx led to reduced PG content in lean animals, but increased PG content in obese rats. Structurally, ACLx caused disorganization of collagenous network in both lean and obese animals through increased collagen orientation in the superficial tissues and a change in the degree of fibrous alignment. However, the cartilage degradation attributed to joint injury and obesity was not necessarily additive when the two risk factors were present simultaneously, particularly for PG content and collagen orientation in the superficial tissues. Interestingly, sham surgeries caused a through-thickness disorganization of collagen network in lean and obese animals. We conclude that the interactions of multiple OA risk factors are complex and their combined effects cannot be understood by superposition principle. Further research is required to elucidate the interactive mechanism between OA subtypes.

Correction: A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

[This corrects the article DOI: 10.1371/journal.pcbi.1009398.].

Dependence of crack shape in loaded articular cartilage on the collagenous structure.

Cartilage cracks disrupt tissue mechanics, alter cell mechanobiology, and often trigger tissue degeneration. Yet, some tissue cracks heal spontaneously. A primary factor determining the fate of tissue cracks is the compression-induced mechanics, specifically whether a crack opens or closes when loaded. Crack deformation is thought to be affected by tissue structure, which can be probed by quantitative polarized light microscopy (PLM). It is unclear how the PLM measures are related to deformed crack morphology. Here, we investigated the relationship between PLM-derived cartilage structure and mechanical behavior of tissue cracks by testing if PLM-derived structural measures correlated with crack morphology in mechanically indented cartilages. METHODS: Knee joint cartilages harvested from mature and immature animals were used for their distinct collagenous fibrous structure and composition. The cartilages were cut through thickness, indented over the cracked region, and processed histologically. Sample-specific birefringence was quantified as two-dimensional (2D) maps of azimuth and retardance, two measures related to local orientation and degree of alignment of the collagen fibers, respectively. The shape of mechanically indented tissue cracks, measured as depth-dependent crack opening, were compared with azimuth, retardance, or "PLM index," a new parameter derived by combining azimuth and retardance. RESULTS: Of the three parameters, only the PLM index consistently correlated with the crack shape in immature and mature tissues. CONCLUSION: In conclusion, we identified the relative roles of azimuth and retardance on the deformation of tissue cracks, with azimuth playing the dominant role. The applicability of the PLM index should be tested in future studies using naturally-occurring tissue cracks.

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.

Dynamic Deformation Calculation of Articular Cartilage and Cells Using Resonance-Driven Laser Scanning Microscopy.

The deformation of articular cartilage and its cells at the micro-scale during dynamic activities such as gait has high mechanoregulatory importance. Measuring the cellular geometries during such dynamics has been limited by the rate of microscopic image acquisition. The introduction of resonating mirrors for image rasterization (resonant scanning), rather than the conventional servo control (galvano scanning), has significantly improved the scanning rate by more than 100×. However, the high scanning rate comes at the cost of image quality, thereby posing challenges in image processing. Here, resonance-driven 3-D laser microscopy is used to observe the transient, micro-scale deformation of articular cartilage and its cells under osmotic challenge conditions. Custom image segmentation and deformable registration software were implemented for analysis of the resonance-scanned microscopy data. The software exhibited robust and accurate performance on the osmotic swelling measurements, as well as quantitative validation testing. The resonance-scanning protocol and developed analysis software allow for simultaneous strain calculation of both the local tissue and cells, and are thus a valuable tool for real-time probing of the cell-matrix interactions that are highly relevant in the fields of orthopedic biomechanics, cell mechanobiology, and functional tissue engineering.

Chondrocyte morphology as an indicator of collagen network integrity.

Osteochondral allograft (OCA) transplantation offers an attractive treatment option as it can be used to repair large cartilage defects that otherwise would not heal. The currently accepted criterion for OCA selection for joint reconstruction is the percentage of viable chondrocytes, but this criterion alone may not be sufficient to ensure structural integrity and functional performance of allografts following transplantation. We sought to determine an additional parameter that indicates matrix integrity. We used multi-photon microscopy to quantitatively assess chondrocyte viability, chondrocyte shape, and collagen structure of articular cartilage of OCAs. Chondrocyte shape varied considerably in otherwise macroscopically healthy-looking OCAs with good (>90%) cell viability. Shape varied from the expected ellipsoidal form found in healthy cartilage, to excessively elongated and flattened cells that often contained multiple cytoplasmic processes reminiscent of those observed in fibroblasts. Chondrocytes with abnormal morphology were associated with degradation of their pericellular matrix and disruption of the collagen fiber orientation, reflected by an increase in heterogeneity of second harmonic signal intensity. Cell shape may be an important marker for collagen network integrity in articular cartilage in general and OCAs specifically. We propose that, aside from cell viability, cell shape may be used as an additional criterion measure for the selection of OCAs. OCAs selected for transplantation based on these criteria showed good graft-host integration post-operation. In view of the rapid and nondestructive nature of the current approach, it may be suitable for clinical application in the future.

Effect of cells on spatial quantification of proteoglycans in articular cartilage of small animals.

OBJECTIVE: Histochemical characterization of proteoglycan content in articular cartilage is important for the understanding of osteoarthritis pathogenesis. However, cartilage cells may interfere with the measurement of matrix proteoglycan content in small animal models (e.g. mice and rats) due to the high cell volume fraction (38%) in mice compared to human tissue (~1%). We investigated whether excluding the cells from image analysis affects the histochemically measured proteoglycan content of rat knee joint cartilage and assessed the effectiveness of a deep learning algorithm-based tool named U-Net in cell segmentation. DESIGN: Histological sections were stained with Safranin-O, after which optical densities were measured using digital densitometry to estimate proteoglycan content. U-Net was trained with 600 annotated Safranin-O cartilage images for exclusion of cells from the cartilage extracellular matrix. Optical densities of the ECM with and without cells were compared as a function of normalized tissue depth. RESULTS: U-Net cell segmentation was accurate, with the measured cell area fraction following largely that of ground-truth images (average difference: 4.3%). Cell area fraction varied as a function of tissue depth and took up 8-21% of the tissue area. The exclusion of cells from the analysis led to an increase in the analyzed depth-dependent optical density of cartilage by approximately 0.6-1.8% (p < 0.01). CONCLUSIONS: Although the effect of cells on the analyzed proteoglycan content is small, it should be considered for improved sensitivity, especially at the onset of the disease during which cells may proliferate in small animals.

The sarcomere force-length relationship in an intact muscle-tendon unit.

The periodic striation pattern in skeletal muscle reflects the length of the basic contractile unit: the sarcomere. More than half a century ago, Gordon, Huxley and Julian provided strong support for the 'sliding filament' theory through experiments with single muscle fibres. The sarcomere force-length (FL) relationship has since been extrapolated to whole muscles in an attempt to unravel in vivo muscle function. However, these extrapolations were frequently associated with non-trivial assumptions, such as muscle length changes corresponding linearly to SL changes. Here, we determined the in situ sarcomere FL relationship in a whole muscle preparation by simultaneously measuring muscle force and individual SLs in an intact muscle-tendon unit (MTU) using state-of-the-art multi-photon excitation microscopy. We found that despite great SL non-uniformity, the mean value of SLs measured from a minute volume of the mid-belly, equivalent to about 5×10-6% of the total muscle volume, agrees well with the theoretically predicted FL relationship, but only if the precise contractile filament lengths are known, and if passive forces from parallel elastic components and activation-associated sarcomere shortening are considered properly. As SLs are not uniformly distributed across the whole muscle and changes in SL with muscle length are location dependent, our results may not be valid for the proximal or distal parts of the muscle. The approach described here, and our findings, may encourage future studies to determine the role of SL non-uniformity in influencing sarcomere FL properties in different muscles and for different locations within single muscles.

On sarcomere length stability during isometric contractions before and after active stretching.

Sarcomere length (SL) instability and SL non-uniformity have been used to explain fundamental properties of skeletal muscles, such as creep, force depression following active muscle shortening and residual force enhancement following active stretching of muscles. Regarding residual force enhancement, it has been argued that active muscle stretching causes SL instability, thereby increasing SL non-uniformity. However, we recently showed that SL non-uniformity is not increased by active muscle stretching, but it remains unclear if SL stability is affected by active stretching. Here, we used single myofibrils of rabbit psoas muscle and measured SL non-uniformity and SL instability during isometric contractions and for isometric contractions following active stretching at average SLs corresponding to the descending limb of the force-length relationship. We defined isometric contractions as contractions during which mean SL remained constant. SL instability was quantified by the rate of change of individual SLs over the course of steady-state isometric force and SL non-uniformity was defined as deviations of SLs from the mean SL at an instant of time. We found that whereas the mean SL remained constant during isometric contraction, by definition, individual SLs did not. SLs were more stable in the force-enhanced, isometric state following active stretching compared with the isometric reference state. We also found that SL instability was not correlated with the rate of change of SL non-uniformity. Also, SL non-uniformity was not different in the isometric and the post-stretch isometric contractions. We conclude that since SL is more stable but similarly non-uniform in the force-enhanced compared with the corresponding isometric reference contraction, it appears unlikely that either SL instability or SL non-uniformity contribute to the residual force enhancement property of skeletal muscle.

Titin force enhancement following active stretch of skinned skeletal muscle fibres.

In actively stretched skeletal muscle sarcomeres, titin-based force is enhanced, increasing the stiffness of active sarcomeres. Titin force enhancement in sarcomeres is vastly reduced in mdm, a genetic mutation with a deletion in titin. Whether loss of titin force enhancement is associated with compensatory mechanisms at higher structural levels of organization, such as single fibres or entire muscles, is unclear. The aim of this study was to determine whether mechanical deficiencies in titin force enhancement are also observed at the fibre level, and whether mechanisms compensate for the loss of titin force enhancement. Single skinned fibres from control and mutant mice were stretched actively and passively beyond filament overlap to observe titin-based force. Mutant fibres generated lower contractile stress (force divided by cross-sectional area) than control fibres. Titin force enhancement was observed in control fibres stretched beyond filament overlap, but was overshadowed in mutant fibres by an abundance of collagen and high variability in mechanics. However, titin force enhancement could be measured in all control fibres and most mutant fibres following short stretches, accounting for ∼25% of the total stress following active stretch. Our results show that the partial loss of titin force enhancement in myofibrils is not preserved in all mutant fibres and this mutation likely affects fibres differentially within a muscle. An increase in collagen helps to reestablish total force at long sarcomere lengths with the loss in titin force enhancement in some mutant fibres, increasing the overall strength of mutant fibres.

Microscale investigation of the anisotropic swelling of cartilage tissue and cells in response to hypo-osmotic challenges.

Tissue swelling represents an early sign of osteoarthritis, reflecting osmolarity changes from iso- to hypo-osmotic in the diseased joints. Increased tissue hydration may drive cell swelling. The opposing cartilages in a joint may swell differently, thereby predisposing the more swollen cartilage and cells to mechanical injuries. However, our understanding of the tissue-cell interdependence in osmotically loaded joints is limited as tissue and cell swellings have been studied separately. Here, we measured tissue and cell responses of opposing patellar (PAT) and femoral groove (FG) cartilages in lapine knees exposed to an extreme hypo-osmotic challenge. We found that the tissue matrix and most cells swelled during the hypo-osmotic challenge, but to a different extent (tissue: <3%, cells: 11%-15%). Swelling-induced tissue strains were anisotropic, showing 2%-4% stretch and 1%-2% compression along the first and third principal directions, respectively. These strains were amplified by 5-8 times in the cells. Interestingly, the first principal strains of tissue and cells occurred in different directions (60-61° for tissue vs. 8-13° for cells), suggesting different mechanisms causing volume expansion in the tissue and the cells. Instead of the continuous swelling observed in the tissue matrix, >88% of cells underwent regulatory volume decrease to return to their pre-osmotic challenge volumes. Cell shapes changed in the early phase of swelling but stayed constant thereafter. Kinematic changes to tissue and cells were larger for PAT cartilage than for FG cartilage. We conclude that the swelling-induced deformation of tissue and cells is anisotropic. Cells actively restored volume independent of the surrounding tissues and seemed to prioritize volume restoration over shape restoration. Our findings shed light on tissue-cell interdependence in changing osmotic environments that is crucial for cell mechano-transduction in swollen/diseased tissues.

Functional Assessment of Human Articular Cartilage Using Second Harmonic Generation (SHG) Imaging: A Feasibility Study.

Many arthroscopic tools developed for knee joint assessment are contact-based, which is challenging for in vivo application in narrow joint spaces. Second harmonic generation (SHG) laser imaging is a non-invasive and non-contact method, thus presenting an attractive alternative. However, the association between SHG-based measures and cartilage quality has not been established systematically. Here, we investigated the feasibility of using image-based measures derived from SHG microscopy for objective evaluation of cartilage quality as assessed by mechanical testing. Human tibial plateaus harvested from nine patients were used. Cartilage mechanical properties were determined using indentation stiffness (Einst) and streaming potential-based quantitative parameters (QP). The correspondence of the cartilage electromechanical properties (Einst and QP) and the image-based measures derived from SHG imaging, tissue thickness and cell viability were evaluated using correlation and logistic regression analyses. The SHG-related parameters included the newly developed volumetric fraction of organised collagenous network (Φcol) and the coefficient of variation of the SHG intensity (CVSHG). We found that Φcol correlated strongly with Einst and QP (ρ = 0.97 and - 0.89, respectively). CVSHG also correlated, albeit weakly, with QP and Einst, (|ρ| = 0.52-0.58). Einst and Φcol were the most sensitive predictors of cartilage quality whereas CVSHG only showed moderate sensitivity. Cell viability and tissue thickness, often used as measures of cartilage health, predicted the cartilage quality poorly. We present a simple, objective, yet effective image-based approach for assessment of cartilage quality. Φcol correlated strongly with electromechanical properties of cartilage and could fuel the continuous development of SHG-based arthroscopy.

The properties of chondrocyte membrane reservoirs and their role in impact-induced cell death.

Impact loading of articular cartilage causes extensive chondrocyte death. Cell membranes have a limited elastic range of 3-4% strain but are protected from direct stretch during physiological loading by their membrane reservoir, an intricate pattern of membrane folds. Using a finite-element model, we suggested previously that access to the membrane reservoir is strain-rate-dependent and that during impact loading, the accessible membrane reservoir is drastically decreased, so that strains applied to chondrocytes are directly transferred to cell membranes, which fail when strains exceed 3-4%. However, experimental support for this proposal is lacking. The purpose of this study was to measure the accessible membrane reservoir size for different membrane strain rates using membrane tethering techniques with atomic force microscopy. We conducted atomic force spectroscopy on isolated chondrocytes (n = 87). A micron-sized cantilever was used to extract membrane tethers from cell surfaces at constant pulling rates. Membrane tethers could be identified as force plateaus in the resulting force-displacement curves. Six pulling rates were tested (1, 5, 10, 20, 40, and 80 µm/s). The size of the membrane reservoir, represented by the membrane tether surface areas, decreased exponentially with increasing pulling rates. The current results support our theoretical findings that chondrocytes exposed to impact loading die because of membrane ruptures caused by high tensile membrane strain rates.

Load-induced fluid pressurisation in hydrogel systems before and after reinforcement by melt-electrowritten fibrous meshes.

Fluid pressure develops transiently within mechanically-loaded, cell-embedding hydrogels, but its magnitude depends on the intrinsic material properties of the hydrogel and cannot be easily altered. The recently developed melt-electrowriting (MEW) technique enables three-dimensional printing of structured fibrous mesh with small fibre diameter (20 µm). The MEW mesh with 20 µm fibre diameter can synergistically increase the instantaneous mechanical stiffness of soft hydrogels. However, the reinforcing mechanism of the MEW meshes is not well understood, and may involve load-induced fluid pressurisation. Here, we examined the reinforcing effect of MEW meshes in three hydrogels: gelatin methacryloyl (GelMA), agarose and alginate, and the role of load-induced fluid pressurisation in the MEW reinforcement. We tested the hydrogels with and without MEW mesh (i.e., hydrogel alone, and MEW-hydrogel composite) using micro-indentation and unconfined compression, and analysed the mechanical data using biphasic Hertz and mixture models. We found that the MEW mesh altered the tension-to-compression modulus ratio differently for hydrogels that are cross-linked differently, which led to a variable change to their load-induced fluid pressurisation. MEW meshes only enhanced the fluid pressurisation for GelMA, but not for agarose or alginate. We speculate that only covalently cross-linked hydrogels (GelMA) can effectively tense the MEW meshes, thereby enhancing the fluid pressure developed during compressive loading. In conclusion, load-induced fluid pressurisation in selected hydrogels was enhanced by MEW fibrous mesh, and may be controlled by MEW mesh of different designs in the future, thereby making fluid pressure a tunable cell growth stimulus for tissue engineering involving mechanical stimulation.

The Protective Function of Directed Asymmetry in the Pericellular Matrix Enveloping Chondrocytes.

The specialized pericellular matrix (PCM) surrounding chondrocytes within articular cartilage is critical to the tissue's health and longevity. Growing evidence suggests that PCM alterations are ubiquitous across all trajectories of osteoarthritis, a crippling and prevalent joint disease. The PCM geometry is of particular interest as it influences the cellular mechanical environment. Observations of asymmetrical PCM thickness have been reported, but a quantified characterization is lacking. To this end, a novel microscopy protocol was developed and applied to acquire images of the PCM surrounding live cells. Morphometric analysis indicated a statistical bias towards thicker PCM on the inferior cellular surface. The mechanical effects of this bias were investigated with multiscale modelling, which revealed potentially damaging, high tensile strains in the direction perpendicular to the membrane and localized on the inferior surface. These strains varied substantially between PCM asymmetry cases. Simulations with a thicker inferior PCM, representative of the observed geometry, resulted in strain magnitudes approximately half of those calculated for a symmetric geometry, and a third of those with a thin inferior PCM. This strain attenuation suggests that synthesis of a thicker inferior PCM may be a protective adaptation. PCM asymmetry may thus be important in cartilage development, pathology, and engineering.

The intrinsic quality of proteoglycans, but not collagen fibres, degrades in osteoarthritic cartilage.

The function of articular cartilage as a load-bearing connective tissue is derived primarily from a balanced interaction between the swelling proteoglycan (PG) matrix and tension-resistant collagen fibrous network. Such balance is compromised during joint disease such as osteoarthritis (OA) due to degradation to PGs and/or collagens. While the PG degradation is generally thought to be related to a loss of protein abundance, the collagenous degradation is more complex as it can be caused independently by a decrease of collagen content, disorganisation of fibrous structure and softening of individual collagen fibrils. A comprehensive understanding of the initial trajectories of degradation of PGs and collagen network can improve our chance of finding potential therapeutic solutions for OA. Here, we developed geometrically, structurally, and compositionally realistic and sample-specific Finite Element (FE) models under the framework of multiphasic mixture theory, from which the elastic moduli of collagen fibres and the PG load-bearing quality in healthy and diseased cartilages were estimated by numerical optimisation of the multi-step indentation stress relaxation force-time curves. We found the intrinsic quality of collagen fibres, measured by their elastic moduli, to stay constant for healthy and diseased cartilages. Combining with previous findings which show unaltered collagen content during early stages of OA, our results suggest the disorganisation of collagen fibrous network as the first form of collagenous degradation in osteoarthritic cartilage. We also found that PG degradation involves not only a loss of protein abundance, but also the quality of the remaining PGs in generating sufficient osmotic pressure for load bearing. This study sheds light on the mechanism of OA pathogenesis and highlights the restoration of collageneous organisation in cartilage as key medical intervention for OA. STATEMENT OF SIGNIFICANCE: Collagen network in articular cartilage consists of individual fibres that are organised into depth-dependent structure specialised for joint load-bearing and lubrication. During osteoarthritis, the collagen network undergoes mechanical degradation, but it is unclear if a loss of content, disorganisation of fibrous structure, or softening of individual fibres causes this degeneration. Using mechanical indentation, Finite Element modelling, and numerical optimisation methods, we determined that individual fibres did not soften in early disease stage. Together with previous findings showing unaltered collagen content, our results pinpoint the disorganisation of collagen structure as the main culprit for early collagenous degradation in osteoarthritic cartilage. Thus, early restoration in cartilage of collagen organisation, instead of individual fibre quality, may be key to slow osteoarthritis development.

Deformation behaviors and mechanical impairments of tissue cracks in immature and mature cartilages.

Degeneration of articular cartilage is often triggered by a small tissue crack. As cartilage structure and composition change with age, the mechanics of cracked cartilage may depend on the tissue age, but this relationship is poorly understood. Here, we investigated cartilage mechanics and crack deformation in immature and mature cartilage exposed to a full-thickness tissue crack using indentation testing and histology, respectively. When a cut was introduced, tissue cracks opened wider in the mature cartilage compared to the immature cartilage. However, the opposite occurred upon mechanical indentation over the cracked region. Functionally, the immature-cracked cartilages stress-relaxed faster, experienced increased tissue strain, and had reduced instantaneous stiffness, compared to the mature-cracked cartilages. Taken together, mature cartilage appears to withstand surface cracks and maintains its mechanical properties better than immature cartilage and these superior properties can be explained by the structure of their collagen fibrous network.

Collagen fibres determine the crack morphology in articular cartilage.

Cracks in articular cartilage compromise tissue integrity and mechanical properties and lead to chondral lesions if untreated. An understanding of the mechanics of cracked cartilage may help in the prevention of cartilage deterioration and the development of tissue-engineered substitutes. The degeneration of cartilage in the presence of cracks may depend on the ultrastructure and composition of the tissue, which changes with aging, disease and habitual loading. It is unknown if the structural and compositional differences between immature and mature cartilage affect the mechanics of cartilage cracks, possibly predisposing one to a greater risk of degeneration than the other. We used a fibre-reinforced poro-viscoelastic swelling material model that accounts for large deformations and tension-compression non-linearity, and the finite element method to investigate the role of cartilage structure and composition on crack morphology and tissue mechanics. We demonstrate that the crack morphology predicted by our theoretical model agrees well with the histo-morphometric images of young and mature cracked cartilages under indentation loading. We also determined that the crack morphology was primarily dependent on collagen fibre orientation which differs as a function of cartilage depth and tissue maturity. The arcade-like collagen fibre orientation, first discussed by Benninghoff in his classical 1925 paper, appears to be beneficial for slowing the progression of tissue cracks by 'sealing' the crack and partially preserving fluid pressure during loading. Preservation of the natural load distribution between solid and fluid constituents of cartilage may be a key factor in slowing or preventing the propagation of tissue cracks and associated tissue matrix damage. STATEMENT OF SIGNIFICANCE: Cracks in articular cartilage can be detrimental to joint health if not treated, but it is not clear how they propagate and lead to tissue degradation. We used an advanced numerical model to determine the role of cartilage structure and composition on crack morphology under loading. Based on the structure and composition found in immature and mature cartilages, our model successfully predicts the crack morphology in these cartilages and determines that collagen fibre as the major determinant of crack morphology. The arcade-like Benninghoff collagen fibre orientation appears to be crucial in 'sealing' the tissue crack and preserves normal fluid-solid load distribution in cartilage. Inclusion of the arcade-like fibre orientation in tissue-engineered construct may help improve its integration within the host tissue.

Sarcomere Lengths Become More Uniform Over Time in Intact Muscle-Tendon Unit During Isometric Contractions.

The seemingly uniform striation pattern of skeletal muscles, quantified in terms of sarcomere lengths (SLs), is inherently non-uniform across all hierarchical levels. The SL non-uniformity theory has been used to explain the force creep in isometric contractions, force depression following shortening of activated muscle, and residual force enhancement following lengthening of activated muscle. Our understanding of sarcomere contraction dynamics has been derived primarily from in vitro experiments using regular bright-field light microscopy or laser diffraction techniques to measure striation/diffraction patterns in isolated muscle fibers or myofibrils. However, the collagenous extracellular matrices present around the muscle fibers, as well as the complex architecture in the whole muscles may lead to different contraction dynamics of sarcomeres than seen in the in vitro studies. Here, we used multi-photon excitation microscopy to visualize in situ individual sarcomeres in intact muscle tendon units (MTUs) of mouse tibialis anterior (TA), and quantified the temporal changes of SL distribution as a function of SLs in relaxed and maximally activated muscles for quasi-steady state, fixed-end isometric conditions. The corresponding muscle forces were simultaneously measured using a force transducer. We found that SL non-uniformity, quantified by the coefficient of variation (CV) of SLs, decreased at a rate of 1.9-3.1%/s in the activated muscles, but remained constant in the relaxed muscles. The force loss during the quasi-steady state likely did not play a role in the decrease of SL non-uniformity, as similar force losses were found in the activated and relaxed muscles, but the CV of SLs in the relaxed muscles underwent negligible change over time. We conclude that sarcomeres in the mid-belly of maximally contracting whole muscles constantly re-organize their lengths into a more uniform pattern over time. The molecular mechanisms accounting for SL non-uniformity appear to differ in active and passive muscles, and need further elucidation, as do the functional implications of the SL non-uniformity.

Single sarcomere contraction dynamics in a whole muscle.

The instantaneous sarcomere length (SL) is regarded as an important indicator of the functional properties of striated muscle. Previously, we found greater sarcomere elongations at the distal end compared to the mid-portion in the mouse tibialis anterior (TA) when the muscle was stretched passively. Here, we wanted to see if SL dispersions increase with activation, as has been observed in single myofibrils, and if SL dispersions differ for different locations in a muscle. Sarcomere lengths were measured at a mid- and a distal location of the TA in live mice using second harmonic generation imaging. Muscle force was measured using a tendon force transducer. We found that SL dispersions increased substantially from the passive to the active state, and were the same for the mid- and distal portions of TA. Sarcomere length non-uniformities within a segment of ~30 serial sarcomeres were up to 1.0 µm. We conclude from these findings that passive, mean SLs obtained from a single location are not necessarily representative of the distribution of SL in active muscle, and thus may be misinterpreted when deriving muscle mechanical properties, such as the force-length relationship. In view of these findings, it seems crucial to determine how SL distributions within a muscle relate to the most fundamental properties of muscle, such as the maximal isometric force.