- Research article
- Open Access
- Open Peer Review
Articular cartilage gene expression patterns in the tissue surrounding the impact site following applications of shear and axial loads
- R. S. McCulloch1, 2Email authorView ORCID ID profile,
- P. L. Mente^1,
- A. T. O’Nan3 and
- M. S. Ashwell3
© The Author(s). 2018
- Received: 24 April 2018
- Accepted: 6 December 2018
- Published: 22 December 2018
Osteoarthritis is a degradative joint disease found in humans and commercial swine which can develop from a number of factors, including prior joint trauma. An impact injury model was developed to deliver in vitro loads to disease-free porcine patellae in a model of OA.
Axial impactions (2000 N normal) and shear impactions (500 N normal with induced shear forces) were delivered to 48 randomly assigned patellae. The patellae were then cultured for 0, 3, 7, or 14 days following the impact. Specimens in the tissue surrounding the loading site were harvested and expression of 18 OA related genes was studied via quantitative PCR. The selected genes were previously identified from published work and fell into four categories: cartilage matrix, degradative enzymes, inflammatory response, and apoptosis.
Type II collagen (Col2a1) showed significantly lower expression in shear vs. axial adjacent tissue at day 0 and 7 (fold changes of 0.40 & 0.19, respectively). In addition, higher expression of degradative enzymes and Fas, an apoptosis gene, was observed in the shear specimens.
The results suggest that a more physiologically valid shear load may induce more damage to surrounding articular cartilage than a normal load alone.
- Articular injury
Although osteoarthritis (OA) is the most common joint disorder in humans , the multi-factorial pathogenesis of OA is still not completely understood [2, 3]. Joint trauma, however, is a known causative factor in the development of OA [2, 4]. Controlled experimental injury models in an in vitro setting [5–11] provide the ability to precisely control loading to study the early stages of articular cartilage degradation in OA. Most impact studies use loading normal to the cartilage surface, however a realistic physiological trauma will most likely generate large shear forces. An axial load will cause a compression of the cartilage surface, and subsequent movement of fluid through the matrix following compression of the tissue. However, our work with three dimension load cells indicates that the while some shear loading is present at the surface, these forces are not large. A realistic trauma that may lead to OA, such as a fall, sports injury, car crash, would likely involve multi-axial loading with much higher shear forces. Thus, it was our aim to intentionally deliver higher shear loads to better model what may happen in a physiological trauma.
ACL transection models are used to study in vivo progression of OA , however these models induce general instability and make it hard to precisely control the degree of load changes. Other models include disuse models of OA , however these likely result in different pathways than that induced by a traumatic injury. Many models use cartilage explants that are removed from the surrounding tissue and then subjected to loading . However, these models make it hard to differentiate whether the changes were due to the harvesting of the explant or the intentionally induced trauma.
Degenerative joint disease is not limited to only humans. Lameness is one of the primary reasons for culling in commercial swine, accounting for 22.5% of the culled sows in Southern China . Structural lameness and leg weakness have been issues in commercial swine for decades and continues to contribute to sow longevity issues . Kirk and coworkers found many of the locomotive problems found in sows were due to arthritis . In our previous work (unpublished), we found that nearly 75% of the femoral heads from culled sows examined in a slaughterhouse have naturally occurring osteoarthritis at different stages of degeneration. Even in growing pigs, degenerative joint disease is a major cause of lameness , where trauma is a common cause of that lameness . Because of the occurrence and similarities of OA in pigs, porcine models of OA have been used in both in vitro and in in vivo models to study OA and its treatments. [20–23]
In our previous work, we developed a shear injury model of OA and evaluated gene expression changes in the articular cartilage directly below the impact injury site . In the present work, expression changes in the tissue adjacent to the loading site were evaluated. It was hypothesized that shear impacts would generate more degradative changes to adjacent tissue than normal impacts.
Gene primer information
Sequence (5′ - > 3′)
Collagen, Type I, Alpha 1 (Col1a1)
Collagen, Type II, Alpha 1 (Col2a1)
SRY (sex determining gene region Y) box-9 (Sox9)
Cartilage oligometric matrix protein (Comp)
Degradative Enzymes & Inhibitors
Matrix metalloprotease-1 (Mmp1)
Matrix metalloprotease-3 (Mmp3)
Matrix metalloprotease-13 (Mmp13)
TIMP Metallopeptidase Inhibitor-1 (Timp1)
TIMP Metallopeptidase Inhibitor-2 (Timp2)
ADAM Metallopeptidase with Thrombospondin Type 1 Motif 5 (Adamts5)
Indian Hedgehog (Ihh)
Transforming growth factor β (Tgfb)
nitric oxide synthase 2, inducible (Inos)
Chitinase-3-like protein 1 (Chi3l1)
Cell Proliferation and Apoptosis
Fas (TNF receptor superfamily, member 6) (Fas)
The relative gene expression levels were compared between the ADJ specimens for axial and shear impacts, and between ADJ and their associated AOI specimens. Differences in AOI specimens were reported previously . A linear mixed model (SAS, SAS Institute, Cary, NC) was used to evaluate differences in fold changes, following the methods of Steibel et al. . The family-wise error rate was controlled using the false discovery rate method (FDR)  to generate a q-value. Because each time point/impact combination had a small number of samples, the q-value significance threshold was set at q < 0.2. This allowed detection of changes without being overly restrictive, and due to how FDR controls for error in results deemed significant, a threshold of up to 0.5 may be acceptable .
Differential gene expression between the tissue adjacent to the impact site (ADJ) and the area of impact (AOI) at each time point for each type of impact (axial and shear)
Axial ADJ vs. AOI
Shear ADJ vs. AOI
Degradative Enzymes & Inhibitors
Inflammatory Response & Signaling
Cell Proliferation & Apoptosis
The goal of this study was to compare two in vitro loading scenarios, one with primarily normal loading, and another with elevated shear loading. The elevated shear loading is likely more representative of true physiological loading. It was believed that shear loading would induce more deleterious effects in the adjacent non-impacted tissue. Our results suggest that this was the case, with Col2a1 downregulated in shear ADJ tissue, Mmps showing higher expression, and elevated levels of Fas expression.
Though Col1a1 levels were elevated early in shear vs axial, the difference was not significant, and the expression did not remain elevated. However, Col2a1 remained lower at all time points in shear ADJ tissue, and tended to be lower than shear AOI. Similarly, Col1a1 was more highly expressed in shear specimens vs. control, and Col2a1 showed lower expression in the shear specimens vs. control. This finding is notable, as Col2a1 is the primary collagen found in articular cartilage, and its depression may indicate that the chondrocytes are not effectively repairing the damage, especially when taken in light of the more highly expressed levels of Col1a1. This could suggest that the chondrocytes are reverting to a dedifferentiated state following the loading [8, 28, 29]. The elevation of Sox9 and Acan indicates a repair effort is being mounted . Taken together with the collagen expression differences, the cells may be attempting recovery from the trauma, but with responses that are ineffective for proper tissue repair.
The chondrocytes show elevated levels of degradative enzymes, and at later time points, elevated levels of an associated inhibitor, consistent with the belief that a repair attempt is being mounted. Early elevation of degradative enzymes (Mmps) in shear vs axial ADJ and in shear ADJ vs AOI, and in shear ADJ vs. control indicates a chondrocyte response to the trauma and attempt to effect repairs, and suggests that shear loads are more damaging, likely due to additional forceful fluid movement through the matrix. This is consistent with other work that has found elevated levels of Mmps in early OA [10, 31]. Also, Adamts5, an aggrecanase, was elevated more in the shear ADJ vs AOI tissue at the last time point, which suggests more matrix breakdown in the shear adjacent tissue. Following the early elevated degradative enzyme expression, the increase of timp1 (an inhibitor of Mmps) at the last time point may suggest the chondrocytes are attempting to limit matrix breakdown, consistent with the idea that the cells are attempting repairs via initial breakdown that is stemmed at later points following the injury.
Ihh has been associated with chondrocyte hypertrophy and has been found to be elevated in OA , and its lower expression in our shear ADJ vs AOI tissue may be the result of the aforementioned dedifferentiation of the chondrocytes, therefore generating inconsistent and ineffective repair efforts. Fas is associated with apoptosis, and its elevation in shear vs axial ADJ tissue may indicate more cell death, in-line with work showing its increase at OA lesions .
If we can gain a better understanding of the very early degenerative changes that precede full blown OA we will be better able to identify potential targets for therapeutic intervention, treat the disease, and prevent the debilitating changes that occur further down the road, having dual benefit in human medicine as well as sustainable pork production. The results presented here show the effects on surrounding tissue when a more realistic loading model is used to simulate OA in vitro. Our findings suggest that shear loads may be more damaging to surrounding tissue than a normal load alone. In the shear specimens, the adjacent tissue showed reduced levels of appropriate collagen expression, and increased expression of degradative enzymes and an apoptosis related gene. This suggests chondrocytes in the shear tissue responded with attempts at repairs but with an ineffective response.
The authors would like to thank Christian Maltecca, Ph.D. for his assistance in developing the mixed model used for analysis of the results.
Funding for this work was provided by National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS) / National Institute of Health Grant 1 R03 AR054557-01A1. The funding agency was not involved in any aspects of the study design, including data collection, interpretation of data or manuscript preparation.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
RSM conducted the experiment and drafted the manuscript. PLM and MSA developed the experimental design, assisted in conducting the experiment, and MSA helped with drafting the manuscript. ATO assisted with sample preparation and results analysis. All authors (with the exception of PLM – deceased) read and approved the final manuscript.
Ethics approval and consent to participate
The samples used in this study were obtained from a local meat processing plant, therefore the animals were slaughtered and introduced into the food supply regardless of this experiment, allowing knowledge to be gained without having to sacrifice additional animals and not requiring IACUC approval. The ethics committee deemed no oversight was required as the tissue was a by-product of normal slaughterhouse operations and no animals were sacrificed for the purpose of this study. Because no oversight was deemed necessary, we did not require written consent, and verbal consent was provided by Andrea Neese of Neese Sausage to use the obtained animal tissue for research purposes.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Arden N, Nevitt MC. Osteoarthritis: epidemiology. Best Pract Res Clin Rheumatol. 2006;20:3–25.View ArticleGoogle Scholar
- Rai MF, Duan X, Quirk JD, Holguin N, Schmidt EJ, Chinzei N, Silva MJ, Sandell LJ. Post-traumatic osteoarthritis in mice following mechanical injury to the synovial joint. Sci Rep. 2017;7:45223.View ArticleGoogle Scholar
- Goldring MB, Culley KL, Otero M. Pathogenesis of osteoarthritis in General. Cham, Switzerland: Springer; 2017.Google Scholar
- Buckwalter JA, Lane NE. Athletics and osteoarthritis. Am J Sports Med. 1997;25:873–81.View ArticleGoogle Scholar
- Ficklin T, Thomas G, Barthel JC, Asanbaeva A, Thonar EJ, Masuda K, Chen AC, Sah RL, Davol A, Klisch SM. Articular cartilage mechanical and biochemical property relations before and after in vitro growth. J Biomech. 2007;40:3607–14.View ArticleGoogle Scholar
- Wilson W, van Burken C, van Donkelaar C, Buma P, van Rietbergen B, Huiskes R. Causes of mechanically induced collagen damage in articular cartilage. J Orthop Res. 2006;24:220–8.View ArticleGoogle Scholar
- Ashwell MS, O'Nan AT, Gonda MG, Mente PL. Gene expression profiling of chondrocytes from a porcine impact injury model. Osteoarthr Cartil. 2008;16:936–46.View ArticleGoogle Scholar
- Aigner T, Dudhia J. Genomics of osteoarthritis. Curr Opin Rheumatol. 2003;15:634–40.View ArticleGoogle Scholar
- Salzmann GM, Nuernberger B, Schmitz P, Anton M, Stoddart MJ, Grad S, Milz S, Tischer T, Vogt S, Gansbacher B. Physicobiochemical synergism through gene therapy and functional tissue engineering for in vitro chondrogenesis. Tissue Eng Part A. 2009;15:2513–24.View ArticleGoogle Scholar
- Aigner T, Zien A, Hanisch D, Zimmer R. Gene expression in chondrocytes assessed with use of microarrays. J Bone Jt Surg. 2003;85:117–23.View ArticleGoogle Scholar
- Kapitanov GI, Ayati BP, Martin JA. Modeling the effect of blunt impact on mitochondrial function in cartilage: implications for development of osteoarthritis. PeerJ. 2017;5:e3468.View ArticleGoogle Scholar
- Maerz T, Newton MD, Kurdziel MD, Altman P, Anderson K, Matthew HW, Baker KC. Articular cartilage degeneration following anterior cruciate ligament injury: a comparison of surgical transection and noninvasive rupture as preclinical models of post-traumatic osteoarthritis. Osteoarthr Cartil. 2016;24(11):1918–27.View ArticleGoogle Scholar
- Wiegant K, Intema F, Van Roermund PM, Barten-van Rijbroek AD, Doornebal A, Hazewinkel HA, Lafeber FP, Mastbergen SC. Evidence of cartilage repair by joint distraction in a canine model of osteoarthritis. Arthritis Rheumatol. 2015;67(2):465–74.View ArticleGoogle Scholar
- Martin JA, Anderson DD, Goetz JE, Fredericks D, Pedersen DR, Ayati BP, Marsh JL, Buckwalter JA. Complementary models of post-traumatic osteoarthritis reveal cellular responses to contact stress that damage articular cartilage and present targets for intervention. J Orthop Res. 2017;35(3):515.View ArticleGoogle Scholar
- Zhao Y, Liu X, Mo D, Chen Q, Chen Y. Analysis of reasons for sow culling and seasonal effects on reproductive disorders in southern China. Anim Reprod Sci. 2015;159:191–7.View ArticleGoogle Scholar
- Crenshaw TD. Arthritis or OCD-identification and prevention. Adv Pork Prod. 2006;17:199.Google Scholar
- Kirk RK, Svensmark B, Ellegaard LP, Jensen HE. Locomotive disorders associated with sow mortality in Danish pig herds. Transbound Emerg Dis. 2005;52:423–8.Google Scholar
- Stavrakakis S, Guy JH, Syranidis I, Johnson GR, Edwards SA. Pre-clinical and clinical walking kinematics in female breeding pigs with lameness: a nested case-control cohort study. Vet J. 2015;205:38–43.View ArticleGoogle Scholar
- Muirhead MR, Alexander TJL. Chapter 9: Managing health in the weaner, grower and finishing periods. In: Carr J, editor. Managing Pig Heath: A Reference for the Farm. 2nd ed; 2013. p. 388.Google Scholar
- Lee KB, Hui JH, Song IC, Ardany L, Lee EH. Injectable mesenchymal stem cell therapy for large cartilage defects—a porcine model. Stem Cells. 2007;25(11):2964–71.View ArticleGoogle Scholar
- Gelse K, Mühle C, Franke O, Park J, Jehle M, Durst K, Göken M, Hennig F, Mark KV, Schneider H. Cell-based resurfacing of large cartilage defects: long-term evaluation of grafts from autologous transgene-activated periosteal cells in a porcine model of osteoarthritis. Arthritis & Rheumatism: official. J Am Coll Rheumatol. 2008;58(2):475–88.Google Scholar
- Goldring MB. The role of cytokines as inflammatory mediators in osteoarthritis: lessons from animal models. Connect Tissue Res. 1999;40(1):1–1.View ArticleGoogle Scholar
- Sieker JT, Proffen BL, Waller KA, Chin KE, Karamchedu NP, Akelman MR, Perrone GS, Kiapour AM, Konrad J, Murray MM, Fleming BC. Transcriptional profiling of articular cartilage in a porcine model of early post-traumatic osteoarthritis. J Orthop Res. 2018;36(1):318–29.PubMedGoogle Scholar
- McCulloch RS, Ashwell MS, Maltecca C, O'Nan AT, Mente PL. Progression of gene expression changes following a mechanical injury to articular cartilage as a model of early stage osteoarthritis. Arthritis. 2014:371426.Google Scholar
- McCulloch RS, Ashwell MS, T O’Nan A, Mente PL. Identification of stable normalization genes for quantitative real-time PCR in porcine articular cartilage. J Anim Sci Biotechnol 2012;3:36.Google Scholar
- Steibel JP, Poletto R, Coussens PM, Rosa GJ. A powerful and flexible linear mixed model framework for the analysis of relative quantification RT-PCR data. Genomics. 2009;94:146–52.View ArticleGoogle Scholar
- Reiner A, Yekutieli D, Benjamini Y. Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003;19:368–75.View ArticleGoogle Scholar
- Xu Y, Barter MJ, Swan DC, Rankin KS, Rowan AD, Santibanez-Koref M, Loughlin J, Young DA. Identification of the pathogenic pathways in osteoarthritic hip cartilage: commonality and discord between hip and knee OA. Osteoarthr Cartil. 2012;20:1029–38.View ArticleGoogle Scholar
- Peffers MJ, Liu X, Clegg PD. Transcriptomic signatures in cartilage ageing. Arthritis Res Ther. 2013;15:R98.View ArticleGoogle Scholar
- Fitzgerald JB, Jin M, Dean D, Wood DJ, Zheng MH, Grodzinsky AJ. Mechanical compression of cartilage explants induces multiple time-dependent gene expression patterns and involves intracellular calcium and cyclic AMP. J Biol Chem. 2004;279:19502–11.View ArticleGoogle Scholar
- Van der Kraan PM, Van den Berg WB. Chondrocyte hypertrophy and osteoarthritis: role in initiation and progression of cartilage degeneration? Osteoarthr Cartil. 2012;20:223–32.View ArticleGoogle Scholar
- Swingler TE, Waters JG, Davidson RK, Pennington CJ, Puente XS, Darrah C, Cooper A, Donell ST, Guile GR, Wang W. Degradome expression profiling in human articular cartilage. Arthritis Res Ther. 2009;11:R96.View ArticleGoogle Scholar
- Kim HA, Lee YJ, Seong S, Choe KW, Song YW. Apoptotic chondrocyte death in human osteoarthritis. J Rheumatol. 2000;27:455–62.PubMedGoogle Scholar