Dietary polyunsaturated fatty acids modulate adipose secretome and is associated with changes in mammary epithelial stem cell self-renewal

Chronic low-grade adipose inflammation, characterized by aberrant adipokine production and pro-inflammatory macrophage activation/polarization is associated with increased risk of breast cancer. Adipocyte fatty acid composition is influenced by dietary availability and may regulate adipokine secretion and adipose inflammation. After feeding F344 rats for 20 weeks with a Western diet or a fish oil-supplemented diet, we cultured primary rat adipose tissue in a three-dimensional explant culture and collected the conditioned medium. The rat adipose tissue secretome was assayed using the Proteome Profiler Cytokine XL Array, and adipose tissue macrophage polarization (M1/M2 ratio) was assessed using the iNOS/ARG1 ratio. We then assessed the adipokine's effects upon stem cell self-renewal using primary human mammospheres from normal breast mammoplasty tissue. Adipose from rats fed the fish oil diet had an ω-3:ω-6 fatty acid ratio of 0.28 compared to 0.04 in Western diet rats. The adipokine profile from the fish oil-fed rats was shifted toward adipokines associated with reduced inflammation compared to the rats fed the Western diet. The M1/M2 macrophage ratio decreased by 50% in adipose of fish oil-fed rats compared to that from rats fed the Western diet. Conditioned media from rats fed the high ω-6 Western diet increased stem cell self-renewal by 62%±9% (X¯%±SD) above baseline compared to only an 11%±11% increase with the fish oil rat adipose. Modulating the adipokine secretome with dietary interventions therefore may alter stromal-epithelial signaling that plays a role in controlling mammary stem cell self-renewal.

Bone marrow adipocytes resist lipolysis and remodeling in response to β-adrenergic stimulation

Bone marrow adipose tissue (BMAT) is preserved or increased in states of caloric restriction. Similarly, we found that BMAT in the tail vertebrae, but not the red marrow in the tibia, resists loss of neutral lipid with acute, 48-hour fasting in rats. The mechanisms underlying this phenomenon and its seemingly distinct regulation from peripheral white adipose tissue (WAT) remain unknown. To test the role of β-adrenergic stimulation, a major regulator of adipose tissue lipolysis, we examined the responses of BMAT to β-adrenergic agonists. Relative to inguinal WAT, BMAT had reduced phosphorylation of hormone sensitive lipase (HSL) after treatment with pan-β-adrenergic agonist isoproterenol. Phosphorylation of HSL in response to β3-adrenergic agonist CL316,243 was decreased by an additional ~90% (distal tibia BMAT) or could not be detected (tail vertebrae). Ex vivo, adrenergic stimulation of lipolysis in purified BMAT adipocytes was also substantially less than iWAT adipocytes and had site-specific properties. Specifically, regulated bone marrow adipocytes (rBMAs) from proximal tibia and femur underwent lipolysis in response to both CL316,243 and forskolin, while constitutive BMAs from the tail responded only to forskolin. This occurred independently of changes in gene expression of β-adrenergic receptors, which were similar between adipocytes from iWAT and BMAT, and could not be explained by defective coupling of β-adrenergic receptors to lipolytic machinery through caveolin 1. Specifically, we found that whereas caveolin 1 was necessary to mediate maximal stimulation of lipolysis in iWAT, overexpression of caveolin 1 was insufficient to rescue impaired BMAT signaling. Lastly, we tested the ability of BMAT to respond to 72-hour treatment with CL316,243 in vivo. This was sufficient to cause beiging of iWAT adipocytes and a decrease in iWAT adipocyte cell size. By contrast, adipocyte size in the tail BMAT and distal tibia remained unchanged. However, within the distal femur, we identified a subpopulation of BMAT adipocytes that underwent lipid droplet remodeling. This response was more pronounced in females than in males and resembled lipolysis-induced lipid partitioning rather than traditional beiging. In summary, BMAT has the capacity to respond to β-adrenergic stimuli, however, its responses are muted and BMAT generally resists lipid hydrolysis and remodeling relative to iWAT. This resistance is more pronounced in distal regions of the skeleton where the BMAT adipocytes are larger with little intervening hematopoiesis, suggesting that there may be a role for both cell-autonomous and microenvironmental determinants. Resistance to β-adrenergic stimuli further separates BMAT from known regulators of energy partitioning and contributes to our understanding of why BMAT is preserved in states of fasting and caloric restriction.

Abstract 5256: Intracellular lipid droplet quantity as a biomarker for obesity associated resistance to omega-3 fatty acid colon cancer prevention

Rationale: Obesity induces resistance to Omega-3 (ω-3) FA anticarcinogenesis effect. Obesity associated excessive fatty acid intracellular concentration results in excess triglycerides stored as intracellular lipid droplets, saturated acyl-co-A transport of FA to triglycerides resulting in excess substrate for cyclooxygenase metabolism with increased formation of PGE2 and other inflammatory eicosanoids.

Methods: Five weeks old male F344 rats were randomized equally into two diet groups. One cohort fed a western diet of ω-6 fatty acids (EPA:ω6 ratio of 0) and other fed an ω-3 fatty acid supplemented diet (EPA:ω6 ratio of 0.4). One half from each diet group were further treated with azoxymethane/dextran sulfate sodium (AOM/DSS) and then euthanized at 21 weeks after the carcinogen induction. Snap frozen colon tissue were assayed for Prostaglandins and fatty acids. OCT fixed colonic tissue stained with oil red O (ORO) to detect lipid droplets which were quantified using Aperio image analysis. From the second set, rat colons were opened longitudinally and fixed in 10% Formalin. Tumors were identified and enumerated using stereomicroscope to get tumor count, incidence and volume. These tissue were subjected to hematoxylin & eosin staining for histological evaluation of colonic tumors.

Results: The ω-3 FA diet significantly increased the percent of eicosapentaenoic acid (EPA) in total fatty acids in both untreated and AOM/DSS-treated rats. EPA in the colonic tissue was increased with ω-3 diet as compared to ω-6 from 1.4±0.8 to 2.6±1.0* in untreated rats and from 1.3±1.3 to 3.3±1.7* in AOM/DSS-treated rats (*p<0.005). PGE2 increased with AOM/DSS treatment. The ω-3 diet decreased PGE2 in colonic mucosa of untreated rats from 132±64 in ω-6 to 78±17 in ω-3 and 150±74 in ω-6 to 92±43 in ω-3 with AOM/DSS treatment. Most tumors were found in the distal colon (DC) with ω-3 diet significantly reducing the number of tumors (p value=0.0011). The ω-3 diet significantly reduced colonic adenoma volume and multiplicity. When tumors were stratified by size and location, there was a significant effect of diet in the DC, with ω-3 diet reducing the number of tumors in the 1-3mm range (p value=0.0018), 3-5mm (p value=0.0041) and >5mm (p value=0.0055). ω-3 also reduced the number of 3-5mm size tumors in the cecum (p value=0.0377). On preliminary quantification of lipid droplets, there was a reduction in the positive pixels divided by the total pixels of ORO stained colonic tissue of ω-3 rats (0.004 ±0.01 as compared to 0.013±0.07 in ω-6 rats; p value=0.0121)

Conclusions: Dietary ω-3 supplementation that increases the EPA:AA ratio in colonic tissue phospholipids reduces excess intracellular FA concentrations reflected by a reduction in size and number of lipid droplets. Lipid droplet number and size may be useful as a biomarker for effective reduction of obesity associated inflammation and carcinogenesis risk.

Citation Format: Muhammad N. Aslam, Nithya Sridhar, Maheen Nadeem, Becky R. Simon, Jianwei Ren, Ananda Sen, James Varani, William L. Smith, Zora Djuric, Dean E. Brenner. Intracellular lipid droplet quantity as a biomarker for obesity associated resistance to omega-3 fatty acid colon cancer prevention [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5256. doi:10.1158/1538-7445.AM2017-5256

Fatty acid and lipidomic data in normal and tumor colon tissues of rats fed diets with and without fish oil

Data is provided to show the detailed fatty acid and lipidomic composition of normal and tumor rat colon tissues. Rats were fed either a Western fat diet or a fish oil diet, and half the rats from each diet group were treated with chemical carcinogens that induce colon cancer (azoxymethane and dextran sodium sulfate). The data show total fatty acid profiles of sera and of all the colon tissues, namely normal tissue from control rats and both normal and tumor tissues from carcinogen-treated rats, as obtained by gas chromatography with mass spectral detection. Data from lipidomic analyses of a representative subset of the colon tissue samples is also shown in heat maps generated from hierarchical cluster analysis. These data display the utility lipidomic analyses to enhance the interpretation of dietary feeding studies aimed at cancer prevention and support the findings published in the companion paper (Effects of fish oil supplementation on prostaglandins in normal and tumor colon tissue: modulation by the lipogenic phenotype of colon tumors, Djuric et al., 2017 [1]).

Effects of fish oil supplementation on prostaglandins in normal and tumor colon tissue: modulation by the lipogenic phenotype of colon tumors

Dietary fish oils have potential for prevention of colon cancer, and yet the mechanisms of action in normal and tumor colon tissues are not well defined. Here we evaluated the impact of the colonic fatty acid milieu on the formation of prostaglandins and other eicosanoids. Distal tumors in rats were chemically induced to model inflammatory colonic carcinogenesis. After 21 weeks of feeding with either a fish oil diet containing an eicosapentaenoic acid/ω-6 fatty acid ratio of 0.4 or a Western fat diet, the relationships between colon fatty acids and prostaglandin E₂ (PGE₂) concentrations were evaluated. PGE₂ is a key proinflammatory mediator in the colon tightly linked with the initiation and progression of colon cancer. The fish oil vs. the Western fat diet resulted in reduced total fatty acid concentrations in serum but not in colon. In the colon, the effects of the fish oil on fatty acids differed in normal and tumor tissue. There were distinct lipodomic patterns consistent with a lipogenic phenotype in tumors. In tumor tissue, the eicosapentaenoic acid/arachidonic acid ratio, cyclooxygenase-2 expression and the mole percent of saturated fatty acids were significant predictors of inter-animal variability in colon PGE₂ after accounting for diet. In normal tissues from either control rats or carcinogen-treated rats, only diet was a significant predictor of colon PGE₂. These results show that the fatty acid milieu can modulate the efficacy of dietary fish oils for colon cancer prevention, and this could extend to other preventive agents that function by reducing inflammatory stress.

Expansion of Bone Marrow Adipose Tissue During Caloric Restriction Is Associated With Increased Circulating Glucocorticoids and Not With Hypoleptinemia

Bone marrow adipose tissue (MAT) accounts for up to 70% of bone marrow volume in healthy adults and increases further in clinical conditions of altered skeletal or metabolic function. Perhaps most strikingly, and in stark contrast to white adipose tissue, MAT has been found to increase during caloric restriction (CR) in humans and many other species. Hypoleptinemia may drive MAT expansion during CR but this has not been demonstrated conclusively. Indeed, MAT formation and function are poorly understood; hence, the physiological and pathological roles of MAT remain elusive. We recently revealed that MAT contributes to hyperadiponectinemia and systemic adaptations to CR. To further these observations, we have now performed CR studies in rabbits to determine whether CR affects adiponectin production by MAT. Moderate or extensive CR decreased bone mass, white adipose tissue mass, and circulating leptin but, surprisingly, did not cause hyperadiponectinemia or MAT expansion. Although this unexpected finding limited our subsequent MAT characterization, it demonstrates that during CR, bone loss can occur independently of MAT expansion; increased MAT may be required for hyperadiponectinemia; and hypoleptinemia is not sufficient for MAT expansion. We further investigated this relationship in mice. In females, CR increased MAT without decreasing circulating leptin, suggesting that hypoleptinemia is also not necessary for MAT expansion. Finally, circulating glucocorticoids increased during CR in mice but not rabbits, suggesting that glucocorticoids might drive MAT expansion during CR. These observations provide insights into the causes and consequences of CR-associated MAT expansion, knowledge with potential relevance to health and disease.

Increased Circulating Adiponectin in Response to Thiazolidinediones: Investigating the Role of Bone Marrow Adipose Tissue

BACKGROUND

Bone marrow adipose tissue (MAT) contributes to increased circulating adiponectin, an insulin-sensitizing hormone, during caloric restriction (CR), but whether this occurs in other contexts remains unknown. The antidiabetic thiazolidinediones (TZDs) also promote MAT expansion and hyperadiponectinemia, even without increasing adiponectin expression in white adipose tissue (WAT).

OBJECTIVES

To test the hypothesis that MAT expansion contributes to TZD-associated hyperadiponectinemia, we investigated the effects of rosiglitazone, a prototypical TZD, in wild-type (WT) or Ocn-Wnt10b mice. The latter resist MAT expansion during CR, leading us to postulate that they would also resist this effect of rosiglitazone.

DESIGN

Male and female WT or Ocn-Wnt10b mice (C57BL/6J) were treated with or without rosiglitazone for 2, 4, or 8 weeks, up to 30 weeks of age. MAT content was assessed by osmium tetroxide staining and adipocyte marker expression. Circulating adiponectin was determined by ELISA.

RESULTS

In WT mice, rosiglitazone caused hyperadiponectinemia and MAT expansion. Compared to WT mice, Ocn-Wnt10b mice had significantly less MAT in distal tibiae and sometimes in proximal tibiae; however, interpretation was complicated by the leakage of osmium tetroxide from ruptures in some tibiae, highlighting an important technical consideration for osmium-based MAT analysis. Despite decreased MAT in Ocn-Wnt10b mice, circulating adiponectin was generally similar between WT and Ocn-Wnt10b mice; however, in females receiving rosiglitazone for 4 weeks, hyperadiponectinemia was significantly blunted in Ocn-Wnt10b compared to WT mice. Notably, this was also the only group in which tibial adiponectin expression was lower than in WT mice, suggesting a close association between MAT adiponectin production and circulating adiponectin. However, rosiglitazone significantly increased adiponectin protein expression in WAT, suggesting that WAT contributes to hyperadiponectinemia in this context. Finally, rosiglitazone upregulated uncoupling protein 1 in brown adipose tissue (BAT), but this protein was undetectable in tibiae, suggesting that MAT is unlikely to share thermogenic properties of BAT.

CONCLUSION

TZD-induced hyperadiponectinemia is closely associated with increased adiponectin production in MAT but is not prevented by the partial loss of MAT that occurs in Ocn-Wnt10b mice. Thus, more robust loss-of-MAT models are required for future studies to better establish MAT's elusive functions, both on an endocrine level and beyond.

Bone marrow adipose tissue is an endocrine organ that contributes to increased circulating adiponectin during caloric restriction

The adipocyte-derived hormone adiponectin promotes metabolic and cardiovascular health. Circulating adiponectin increases in lean states such as caloric restriction (CR), but the reasons for this paradox remain unclear. Unlike white adipose tissue (WAT), bone marrow adipose tissue (MAT) increases during CR, and both MAT and serum adiponectin increase in many other clinical conditions. Thus, we investigated whether MAT contributes to circulating adiponectin. We find that adiponectin secretion is greater from MAT than WAT. Notably, specific inhibition of MAT formation in mice results in decreased circulating adiponectin during CR despite unaltered adiponectin expression in WAT. Inhibiting MAT formation also alters skeletal muscle adaptation to CR, suggesting that MAT exerts systemic effects. Finally, we reveal that both MAT and serum adiponectin increase during cancer therapy in humans. These observations identify MAT as an endocrine organ that contributes significantly to increased serum adiponectin during CR and perhaps in other adverse states.

Administration of saccharin to neonatal mice influences body composition of adult males and reduces body weight of females

Nutritional or pharmacological perturbations during perinatal growth can cause persistent effects on the function of white adipose tissue, altering susceptibility to obesity later in life. Previous studies have established that saccharin, a nonnutritive sweetener, inhibits lipolysis in mature adipocytes and stimulates adipogenesis. Thus, the current study tested whether neonatal exposure to saccharin via maternal lactation increased susceptibility of mice to diet-induced obesity. Saccharin decreased body weight of female mice beginning postnatal week 3. Decreased liver weights on week 14 corroborated this diminished body weight. Initially, saccharin also reduced male mouse body weight. By week 5, weights transiently rebounded above controls, and by week 14, male body weights did not differ. Body composition analysis revealed that saccharin increased lean and decreased fat mass of male mice, the latter due to decreased adipocyte size and epididymal, perirenal, and sc adipose weights. A mild improvement in glucose tolerance without a change in insulin sensitivity or secretion aligned with this leaner phenotype. Interestingly, microcomputed tomography analysis indicated that saccharin also increased cortical and trabecular bone mass of male mice and modified cortical bone alone in female mice. A modest increase in circulating testosterone may contribute to the leaner phenotype in male mice. Accordingly, the current study established a developmental period in which saccharin at high concentrations reduces adiposity and increases lean and bone mass in male mice while decreasing generalized growth in female mice.

Sweet taste receptor deficient mice have decreased adiposity and increased bone mass

Functional expression of sweet taste receptors (T1R2 and T1R3) has been reported in numerous metabolic tissues, including the gut, pancreas, and, more recently, in adipose tissue. It has been suggested that sweet taste receptors in these non-gustatory tissues may play a role in systemic energy balance and metabolism. Smaller adipose depots have been reported in T1R3 knockout mice on a high carbohydrate diet, and sweet taste receptors have been reported to regulate adipogenesis in vitro. To assess the potential contribution of sweet taste receptors to adipose tissue biology, we investigated the adipose tissue phenotypes of T1R2 and T1R3 knockout mice. Here we provide data to demonstrate that when fed an obesogenic diet, both T1R2 and T1R3 knockout mice have reduced adiposity and smaller adipocytes. Although a mild glucose intolerance was observed with T1R3 deficiency, other metabolic variables analyzed were similar between genotypes. In addition, food intake, respiratory quotient, oxygen consumption, and physical activity were unchanged in T1R2 knockout mice. Although T1R2 deficiency did not affect adipocyte number in peripheral adipose depots, the number of bone marrow adipocytes is significantly reduced in these knockout animals. Finally, we present data demonstrating that T1R2 and T1R3 knockout mice have increased cortical bone mass and trabecular remodeling. This report identifies novel functions for sweet taste receptors in the regulation of adipose and bone biology, and suggests that in these contexts, T1R2 and T1R3 are either dependent on each other for activity or have common independent effects in vivo.

Artificial sweeteners stimulate adipogenesis and suppress lipolysis independently of sweet taste receptors

G protein-coupled receptors mediate responses to a myriad of ligands, some of which regulate adipocyte differentiation and metabolism. The sweet taste receptors T1R2 and T1R3 are G protein-coupled receptors that function as carbohydrate sensors in taste buds, gut, and pancreas. Here we report that sweet taste receptors T1R2 and T1R3 are expressed throughout adipogenesis and in adipose tissues. Treatment of mouse and human precursor cells with artificial sweeteners, saccharin and acesulfame potassium, enhanced adipogenesis. Saccharin treatment of 3T3-L1 cells and primary mesenchymal stem cells rapidly stimulated phosphorylation of Akt and downstream targets with functions in adipogenesis such as cAMP-response element-binding protein and FOXO1; however, increased expression of peroxisome proliferator-activated receptor γ and CCAAT/enhancer-binding protein α was not observed until relatively late in differentiation. Saccharin-stimulated Akt phosphorylation at Thr-308 occurred within 5 min, was phosphatidylinositol 3-kinase-dependent, and occurred in the presence of high concentrations of insulin and dexamethasone; phosphorylation of Ser-473 occurred more gradually. Surprisingly, neither saccharin-stimulated adipogenesis nor Thr-308 phosphorylation was dependent on expression of T1R2 and/or T1R3, although Ser-473 phosphorylation was impaired in T1R2/T1R3 double knock-out precursors. In mature adipocytes, artificial sweetener treatment suppressed lipolysis even in the presence of forskolin, and lipolytic responses were correlated with phosphorylation of hormone-sensitive lipase. Suppression of lipolysis by saccharin in adipocytes was also independent of T1R2 and T1R3. These results suggest that some artificial sweeteners have previously uncharacterized metabolic effects on adipocyte differentiation and metabolism and that effects of artificial sweeteners on adipose tissue biology may be largely independent of the classical sweet taste receptors, T1R2 and T1R3.

A one-dimensional mixed porohyperelastic transport swelling finite element model with growth

A one-dimensional, large-strain, mixed porohyperelastic transport and swelling (MPHETS) finite element model was developed in MATLAB and incorporated with a well-known growth model for soft tissues to allow the model to grow (increase in length) or shrink (decrease in length) at constant material density. By using the finite element model to determine the deformation and stress state, it is possible to implement different growth laws in the program in the future to simulate how soft tissues grow and behave when exposed to various stimuli (e.g. mechanical, chemical, or electrical). The essential assumptions needed to use the MPHETS model with growth are clearly identified and explained in this paper. The primary assumption in this work, however, is that the stress upon which growth acts is the stress in the solid skeleton, i.e. the effective stress, S(eff). It is shown that significantly different amounts of growth are experienced for the same loading conditions when using a porohyperelastic model as compared to a purely solid model. In one particular example, approximately 51% less total growth occurred in the MPHETS model than in the solid model even though both problems were subjected to the same external loading. This work represents a first step in developing more sophisticated models capable of capturing the complex mechanical and biochemical environment in growing and remodeling tissues.

Coupled porohyperelastic mass transport (PHEXPT) finite element models for soft tissues using ABAQUS

Finite element models (FEMs) including characteristic large deformations in highly nonlinear materials (hyperelasticity and coupled diffusive/convective transport of neutral mobile species) will allow quantitative study of in vivo tissues. Such FEMs will provide basic understanding of normal and pathological tissue responses and lead to optimization of local drug delivery strategies. We present a coupled porohyperelastic mass transport (PHEXPT) finite element approach developed using a commercially available ABAQUS finite element software. The PHEXPT transient simulations are based on sequential solution of the porohyperelastic (PHE) and mass transport (XPT) problems where an Eulerian PHE FEM is coupled to a Lagrangian XPT FEM using a custom-written FORTRAN program. The PHEXPT theoretical background is derived in the context of porous media transport theory and extended to ABAQUS finite element formulations. The essential assumptions needed in order to use ABAQUS are clearly identified in the derivation. Representative benchmark finite element simulations are provided along with analytical solutions (when appropriate). These simulations demonstrate the differences in transient and steady state responses including finite deformations, total stress, fluid pressure, relative fluid, and mobile species flux. A detailed description of important model considerations (e.g., material property functions and jump discontinuities at material interfaces) is also presented in the context of finite deformations. The ABAQUS-based PHEXPT approach enables the use of the available ABAQUS capabilities (interactive FEM mesh generation, finite element libraries, nonlinear material laws, pre- and postprocessing, etc.). PHEXPT FEMs can be used to simulate the transport of a relatively large neutral species (negligible osmotic fluid flux) in highly deformable hydrated soft tissues and tissue-engineered materials.

Identification and determination of material properties for porohyperelastic analysis of large arteries

A "porohyperelastic" (PHE) material model is described and the theoretical framework presented that allows identification of the necessary material properties functions for soft arterial tissues. A generalized Fung form is proposed for the PHE constitutive law in which the two fundamental Lagrangian material properties are the effective strain energy density function, W(e), and the hydraulic permeability, kij. The PHE model is based on isotropic forms using W(e) = Ue (phi) = 1/2C0(e phi - 1) and the radial component of permeability, kRR = kRR(phi), with phi = C1'(I1 - 3) + C2'(I2 - 3) + K'(J - 1)2. The methods for determination of these material properties are illustrated using experimental data from in situ rabbit aortas. Three experiments are described to determine parameters in Ue and kRR for the intima and media of the aortas, i.e., (1) undrained tests to determine C0, C1', and C2'; (2) drained tests to determine K'; and (3) steady-state pressurization tests of intact and de-endothelialized vessels to determine intimal and medial permeability (adventitia removed in these models). Data-reduction procedures are presented that allow determination of kRR for the intima and media and Ue for the media using experimental data. The effectiveness and accuracy of these procedures are studied using input "data" from finite element models generated with the ABAQUS program. The isotropic theory and data-reduction methods give good approximations for the PHE properties of in situ aortas. These methods can be extended to include arterial tissue remodeling and anisotropic behavior when appropriate experimental data are available.

A poroelastic finite element formulation including transport and swelling in soft tissue structures

A field theory is presented for the study of swelling in soft tissue structures that are modeled as poroelastic materials. As a first approximation, soft tissues are assumed to be linear isotropic materials undergoing infinitesimal strains. Material properties are identified that are necessary for the solution of initial boundary value problems where swelling and convection are significant. A finite element model is developed that includes the solid displacements, the relative fluid displacements, and a representative concentration as the primary unknowns. A numerical example is presented based on a triphasic model. The finite model simulates a typical experimental protocol for soft tissue testing and demonstrates the interaction and coupling associated with relative fluid motion and swelling in a deforming poroelastic materiaL The theory and finite element model provide a starting point for nonlinear porohyperelastic transport-swelling analyses of soft tissue structures that include finite strains in anisotropic materials.

A dynamic material parameter estimation procedure for soft tissue using a poroelastic finite element model

A three-dimensional finite element model for a poroelastic medium has been coupled with a least squares parameter estimation method for the purpose of assessing material properties based on intradiscal displacement and reactive forces. Parameter optimization may be based on either load or displacement control experiments. In this paper we present the basis of the finite element model and the parameter estimation process. The method is then applied to a test problem and the computational behavior is discussed. Sequential optimization on different parameter groups was found to have superior convergence properties. Some guidelines for choosing the starting parameter values for optimization were deduced by considering the form of the objective function. For load control experiments, in which displacement data is used for the optimization, the starting values for the elastic modulus should be lower in magnitude than an "anticipated" modulus. The permeability starting values should be higher than an anticipated permeability. For displacement control experiments, the reverse is true. The optimization scheme was also tested on data with random variations.

Finite element models for arterial wall mechanics

Arterial wall mechanics has been studied for nearly 200 years. This subject is of importance if we are to gain a fundamental understanding of this complex biological structure, as well as information needed to design prosthetics. Biomechanical arterial models continue to play an important role in the study of atherosclerosis, a disease of the arterial wall that is the chief cause of mortality and morbidity in the United States and the Western World. Over the past 20 years, the finite element model (FEM) has been used in a variety of ways to simulate the structural response of large arteries. Our purpose is to summarize the uses of FEMs in arterial mechanics. We will also indicate directions for future research in this area. A specialized FEM was described in the literature for the study of transport in the arterial wall, however the convection was not directly linked to arterial wall mechanics. In this paper special attention will be given to the development of FEMs based on the poroelastic view of arterial tissues which couple wall deformation, free tissue fluid motion, and associated transport phenomena in the arterial wall. In the future such models should provide fundamental quantitative information relating arterial wall mechanics and transport which may lead to a better understanding of both normal arterial physiology and atherogenesis.

A poroelastic-swelling finite element model with application to the intervertebral disc

The swelling process that occurs in soft tissue is incorporated into a poroelastic finite element model. The model is applied to a spinal segment consisting of two vertebrae and a single intervertebral disc. The theory is an extension of the poroelastic theory developed by Biot and the model is an adaptation of an axisymmetric poroelastic finite element model of the intervertebral disc by Simon. The model is completely three-dimensional although the results presented here assume symmetry about the sagittal plane. The theory is presented in two stages. First the development of the poroelastic theory. Following this, the effects of swelling caused by osmotic pressure are developed and expressed as a modification of the constitutive law and initial stresses. In the case of the disc, this pressure is produced mainly by the fixed negative charges on the proteoglycans within the disc. In this development we assume that the number of fixed charges remains constant over time and that the distribution of mobile ions has reached equilibrium. The variations over time in osmotic pressure, and thus in swelling effects are therefore only dependent on the initial state and the change in water content. Variations of the swelling effects caused by changes in mobile ion concentrations will be the subject of a future paper. The results reported in this article illustrate the dramatic effect of swelling on the load carrying mechanisms in the disc. The authors believe it is likely that this will have important useful implications for our understanding not only of normal disc function, but also of abnormal function, such as disc degeneration, herniation, and others.

Arterial mechanics in spontaneously hypertensive rats. Mechanical properties, hydraulic conductivity, and two-phase (solid/fluid) finite element models

To characterize the interaction between mechanical and fluid transport properties in hypertension, we measured in vivo elastic material constants and hydraulic conductivity in intact segments of carotid arteries in normal and spontaneously hypertensive rats (SHR). With the use of a finite element model, the arterial wall was modeled as a large-deformation, two-phase (solid/fluid) medium, which accounts for the existence and motion of the tissue fluid. Measurements of internal diameter and transmural pressures were obtained during continuous increases in pressure from 0 to 200 mm Hg. Strain and stress components were calculated based on a pseudostrain exponential energy density function. To measure the hydraulic conductivity, segments of the carotid artery were isolated, filled with a 4% oxygenated albumin-Tyrode's solution, and connected to a capillary tube. The movement of the meniscus of the capillary tube represented the fluid filtration across the artery. To study the influence of transmural pressure on hydraulic conductivity, measurement of fluid filtration across the arterial wall was obtained at transmural pressures of 50 and 100 mm Hg. The material constants in the SHR (n = 9) were higher (p less than 0.05 for all variables) than in normal rats (n = 10): c = 1,343 +/- 96 versus 1,158 +/- 65 mm Hg, b1 = 1.84 +/- 0.24 versus 1.22 +/- 0.22, b2 = 0.769 +/- 0.114 versus 0.616 +/- 0.11, b3 = 0.017 +/- 0.005 versus 0.0065 +/- 0.002, b4 = 0.206 +/- 0.04 versus 0.083 +/- 0.03, b5 = 0.0594 +/- 0.007 versus 0.0217 +/- 0.006, and b6 = 0.22 +/- 0.09 versus 0.123 +/- 0.02, respectively. The hydraulic conductivity of the total wall, calculated from the filtration data, was lower (p less than 0.05) at both 50 and 100 mm Hg in the SHR (n = 6) compared with normal rats (n = 7): 1.12 +/- 0.31 x 10(-8) and 0.72 +/- 0.23 x 10(-8) versus 1.95 +/- 0.53 x 10(-8) and 1.35 +/- 0.47 x 10(-8) cm/(sec.mm Hg), respectively. The intergroup comparisons between 50 and 100 mm Hg in both SHR and normal rats were also different (p less than 0.05). The finite element model was used to predict tissue fluid pressure distribution, tissue fluid velocity distribution, and total Cauchy stress gradients developed in the arterial wall during fluid pressurization in both species.(ABSTRACT TRUNCATED AT 400 WORDS)

Effect of pressure on aortic hydraulic conductance

This study was performed to determine whether the transmural hydraulic conductance (Lp) of the rabbit aortic wall depends on its distension. In 19 rabbits, the aorta was cannulated in situ and perfused at a given pressure with a physiologically buffered solution containing 4% bovine serum albumin. The output cannula was then occluded to limit fluid flow to that traversing the artery wall. External diameter and transmural fluid flow were measured at three pressures (eight rabbits, group 1) or at four pressures (12 rabbits, group 2) in each vessel. Transmural fluid flow was determined by monitoring the velocity of an air bubble within a buffer-filled tube leading to the input cannula. From group 1 measurements, Lp values (mean +/- SD) at 50, 100, and 150 mm Hg were calculated to be 3.8 +/- 2.8, 3.5 +/- 1.3, and 4.1 +/- 1.2 x 10(-8) cm/sec/mm Hg, respectively. Group 2 measurements gave values of 4.2 +/- 1.6, 3.8 +/- 1.1, 3.8 +/- 1.1, and 4.2 +/- 1.1 x 10(-8) cm/sec/mm Hg at 75, 100, 125, and 150 mm Hg, respectively. Paired Student's t tests indicated no significant change in Lp with pressure. However, linear regression analysis demonstrated a weak correlation between Lp values obtained at 50 and 100 mm Hg (r2 = 0.30) and at 75 and 100 mm Hg (r2 = 0.36). Values of Lp at 100 and 150 mm Hg and at 125 and 150 mm Hg were closely correlated in each case. These results suggest that between 50 and 100 mm Hg the structural properties of the aortic wall change so as to alter Lp but not in the same way in each vessel. Lp may increase or decrease depending on which structural change predominates in a particular vessel.

Structural models for human spinal motion segments based on a poroelastic view of the intervertebral disk

Analytical and finite element models (FEMs) were used to quantify poroelastic material properties for a human intervertebral disk. An axisymmetric FEM based on a poroelastic view of disk constituents was developed for a representative human spinal motion segment (SMS). Creep and steady-state response predicted by FEMs agreed with experimental observations, i.e., long-time creep occurs with flow in the SMS, whereas for rapid steady-state loading an "undrained," nearly incompressible response is evident. A relatively low value was determined for discal permeability. Transient and long-term creep FE analyses included the study of deformation, pore fluid flow, stress, and pore fluid pressure. Relative fluid motion associated with transient creep is related to nuclear nutrition and the overall mechanical response in the normal disk. Degeneration of the disk may be associated with an increase in permeability.

Poroelastic dynamic structural models of rhesus spinal motion segments

Finite element models (FEMs) and analytical and experimental models based on poroelastic constitutive laws were developed for rhesus spinal motion segments (SMSs). Long-time creep, transient creep, and impact were studied for SMSs with normal and simulated degenerated discs. The results suggested that long-time creep observed in excised SMSs may be reduced in the in vivo SMS. The fluid phase included in these FEMs was shown to play a significant role in the mechanical response of SMSs. Relative fluid motion fields predicted in the SMS could be related to nutritional paths to the avascular interior of the disc and were found to be very sensitive to changes in discal stiffness. Reduced disc height, increased discal bulge, altered fluid motion, and stresses were quantified and may be related to mechanical failure, disc degeneration, and low-back pain.

Relaxation and creep quasilinear viscoelastic models for normal articular cartilage

A quasilinear viscoelastic model was used to develop relaxation and creep forms for a constitutive law for soft tissues. Combined relaxation and cyclic test data as well as preconditioned and nonpreconditioned creep data were used to demonstrate the approach for normal bovine articular cartilage. Values for mechanical parameters in the analytical models were determined using a generalized least squares method.

Less rigid internal fixation plates: historical perspectives and new concepts

In the application of "rigid" plates for diaphyseal fractures, lack of callus healing and overprotection of the underlying bone are viewed by many investigators as undesirable consequences. Potential solutions offered to overcome these deficiencies include modification of the timing of plate removal, use of biologically degradable materials for plates so that stress-shielding can be minimized, and use of less rigid plate fixation systems. This study emphasizes the selection of appropriate design criteria for less rigid plate-fixation systems. To accomplish this goal, the axial, bending, and torsional stiffness parameters are considered in place of the oversimplified terms such as "flexible plate" or "elastic fixation." With the aid of finite element modeling and simplified bench experiments, we performed parametric studies and singled out the plate axial stiffness as the dominant factor in altering the bone stresses. As a result, we designed two experimental plates, one a thin Ti-6Al-4V (titanium with 6% aluminum and 4% vanadium) alloy plate with low stiffness in axial and bending directions, and the other a tubular stainless steel plate with low stiffness in the axial direction but moderate stiffness in bending and torsional directions. The low-stiffness Ti-6Al-4V alloy plate was first tested in a demanding bilateral canine midshaft osteotomy, and proved to be inadequate. Both experimental plates were successful in the unilateral osteotomies, with the tubular plate yielding the best results. After 6 months of plating, the bones beneath the tubular plate had superior mechanical and structural properties as compared to those of the control "rigid" stainless steel and the Ti-6Al-4V alloy plates. Application of this plate prolonged for 9 months did not cause reduction in bone properties and strength. The success of the tubular plate is due to its moderate bending and torsional stiffnesses, which provide adequate fixation to achieve callus union, while its low axial stiffness permits the underlying bone to share the physiological stresses needed for bone remodeling. These drastic changes in mechanical demands on the internal fixation plate during the early healing phase and the postunion remodeling phase are discussed.

A new approach to the design of internal fixation plates

Mathematical analysis by finite element modeling methods was used in conjunction with laboratory bench experiments in selecting appropriate stiffness parameters for internal fixation plate designs. Bench experiments performed included tests of an idealized plated magnesium tube and of a plated canine femur. Finite element modeling of the plate-tube and plate-bone structures was also performed, and the computed results were compared with those obtained experimentally. Three types of internal fixation plates--one "rigid" control plate and two less rigid experimental plates--were examined using both the finite element and experimental models. These plates were further examined by finite element modeling of a plated human femur during single stance loading. In all cases, the "rigid" control plate was found to shield the underlying bone from stress, while both experimental plates were found to have significantly less bone stress shielding. Our findings suggest that an improved plate design should have a low axial stiffness but moderate bending and torsional stiffnesses to facilitate fracture healing and bone remodeling without causing osteopenia.

Quasi-linear viscoelastic properties of normal articular cartilage

A combined experimental and analytical approach was used to determine the history-dependent viscoelastic properties of normal articular cartilage in tension. Specimens along the surface split line direction, taken from the middle zone of articular cartilage were subjected to relaxation and cyclic tests. A quasi-linear viscoelastic theory proposed by Fung was used in combination with the experimental results to determine the nonlinear viscoelastic properties and the elastic stress-strain relationship of normal articular cartilage.

Reevaluation of arterial constitutive relations. A finite-deformation approach

The purpose of this investigation was to use the finite-deformation theory of elasticity to interpret pressure-diameter data for in situ canine aortas and other arterial response data reported in the literature. A meaningful mechanical property for arterial tissue was identified as ∂ W 1 /∂ I , the partial derivative of the strain-energy density function with respect to the first strain invariant. An exponential function was found to characterize the mechanical property ∂W 1 /∂I for all arteries considered. Thin-walled tube stress approximations were found to result in inaccurate values for arterial stresses and incremental elastic mechanical properties. Wave speeds calculated using ∂ W 1 /∂ I for these arterial tissues agreed well with experimental measurements of wave speeds reported in the literature. Elevated values for strain-energy density were found in the inner arterial tissue layers. These high values for strain energy may contribute to atherogenesis in relatively straight arteries (e.g., the abdominal aorta) subjected to hypertension.