Journal of Applied Physiology

Gene profiling of the rat medial collateral ligament during early healing using microarray analysis

Connie S. Chamberlain, Sabrina H. Brounts, David G. Sterken, Kevin I. Rolnick, Geoffrey S. Baer, Ray Vanderby


Ligament heals in a synchronized and complex series of events. The remodeling process may last months or years. Experimental evidence suggests the damaged ligament does not recover its normal functional properties. Specific mechanisms to prevent scar formation and to regenerate the original mechanical function remain elusive but likely involve regulation of creeping substitution. Creeping substitution creates a larger hypercellular, hypervascular, and disorganized granulation tissue mass that results in an inefficient and nonregenerative wound healing process for the ligament. Control of creeping substitution may limit the extent of this tissue compromise and reduce the time necessary for healing. The objective of this study is to better understand the mechanism behind scar formation by identifying the extracellular matrix factors and other unique genes of interest differentially expressed during rat ligament healing via microarray. For this study, rat medial collateral ligaments were either surgically transected or left intact. Ligaments were collected at day 3 or 7 postinjury and used for microarray, quantitative PCR, and/or immunohistochemistry. Results were compared with the normal intact ligament. We demonstrate that early ligament healing is characterized by the modulation of several inflammatory and extracellular matrix factors during the first week of injury. Specifically, a number of matrix metalloproteinases and collagens are differentially and significantly expressed during early ligament healing. Additionally, we demonstrate the modulation of three novel genes, periostin, collagen-triple helix repeat containing-1, and serine protease 35 in our ligament healing model. Together, control of granulation tissue creeping substitution and subsequent downstream scar formation is likely to involve these factors.

  • collagen
  • matrix metalloproteinase
  • periostin

ligament healing involves a complex, coordinated series of events that results in a ligament that is more scar-like in character than the native tissue. The repair process may extend from months to years. Experimental evidence suggests the injured ligament incompletely recovers its original mechanical properties (25, 28). In an effort to develop an efficient repair process, researchers have tested a number of treatments including tissue engineering approaches, nonsteroidal anti-inflammatory drugs, and ultrasonic or electrical stimulation (28, 35). All show promise, but an incomplete understanding of the healing process makes optimal treatment regimes elusive. With tissue engineering and regenerative medicine, it is essential to understand the normal healing process in a ligament, thereby providing a basis to formulate and evaluate innovative new treatments.

Previous work from our laboratory (8) established a highly active cellular and vascular scenario within the first week after ligament injury. Specifically, macrophage infiltration crested between days 3–5 while blood vessel accumulation peaked between days 7–11 postinjury (Fig. 1). Creeping substitution of damaged tissue into normal uninjured ligament, resulting in the expansion of the original wound size, was also evident during this time (8). The reduction of creeping substitution during inflammation and early proliferation may limit downstream scar formation, but the genetic players remain largely undefined.

Fig. 1.

Typical angiogenic, granulation tissue and macrophage response of a wounded medial collateral ligament (MCL) at 3, 5, and 7 days postinjury. Microangiography of the healing MCL demonstrates low blood vessel infiltration 3 days postinjury (top). By day 5 and 7, blood vessels are prominent. Granulation tissue forms by day 3 and continues to expand within the ligament beyond day 7 (i.e., creeping substitution; middle). Macrophages also appear early during healing and peak around day 5 and 7 (bottom).

Matrix metalloproteinases (MMP) regulate normal tissue remodeling of the extracellular matrix (ECM) by degrading collagen and regulating cell recruitment into tissues. Neutrophils, macrophages, and fibroblasts release various MMPs, modulating inflammation and tissue degradation. During pathologic inflammatory conditions, excessive proteolytic activity is produced to facilitate tissue injury; rheumatoid arthritis and osteoarthritis are associated with increased MMP activity and subsequent collagen and proteoglycan degradation. Control of creeping substitution and subsequent scar formation may therefore modify collagen degradation by the MMPs. Therefore, the objective of this study is to better delineate the genetic factors involved during early ligament healing via microarray analysis. Specifically, this study focuses on identification of MMPs and collagens during early healing.


Animal model.

This study was approved by the University of Wisconsin Institutional Animal Use and Care Committee. Forty-two skeletally mature Wistar rats were used as an animal model for ligament healing. Animals were randomly placed in one of three groups (n = 14 animals/group): intact control, day 3 postinjury, or day 7 postinjury. Rats were anesthetized via isofluorane. A surgically transected, rather than torn, medial collateral ligament (MCL) was used as an experimental model to create a uniform defect for healing. Rats were subjected to bilateral MCL transection using sterile techniques (n = 37 animals). A 1-cm skin incision was made over the medial aspect of each stifle, and the subcutaneous tissue was dissected to expose the sartorius muscle and underlying MCL. The midpoint of the MCL was completely transected. Severed ends were not sutured and allowed to retract. The muscular, subcutaneous, and subdermal tissue layers were each closed with 4–0 Dexon suture. The remaining 14 animals underwent no transection and served as age-matched intact controls. All animals were allowed unrestricted cage movement immediately after surgery. At day 3 or 7 postinjury (n = 14 animals/day), the MCLs were carefully dissected, measured, weighed, and immediately placed in liquid nitrogen. MCLs from each time point were divided into three pools (3 MCLs/pool) and used for microarray analysis. The contralateral MCLs were used for quantitative PCR (qPCR). MCLs from the remaining five animals per day were snap frozen in optimal cutting temperature media and used for immunohistochemistry.

RNA isolation for microarray and qPCR analysis.

Total RNA was isolated from the MCLs using methods similar to Reno et al. (56), combining the TRIzol (Invitrogen, Carlsbad, CA) method with column fractionalization steps of the RNeasy total RNA kit (Qiagen, Valencia, CA). Tissue was homogenized using a PowerGen 500 homogenizer (Fisher Scientific, Pittsburgh, PA). Subsequent to RNA isolation and yield and purity of RNA were quantified by nanodrop spectrophotometric measurement at 260 nm (Nanodrop Technologies, Wilmington, DE).

Purification and preparation of RNA for affymetrix GeneChip analysis.

Nucleic acid qc, labeling, and hybridizations were performed by the Gene Expression Center, a core research facility within the University of Wisconsin Biotechnology Center. The quality and quantity of each RNA sample were analyzed on the Agilent 2100 BioAnalyzer using the RNA6000 NanoChip and the NanoDrop 1000 Spectrophotometer (Thermo Scientific, Wilmington, DE), respectively. Total RNA was converted to double-stranded cDNA using an oligo dT primer containing the T7 RNA polymerase promoter (5′-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-T24-3′). Briefly, first-strand synthesis was performed by incubating 3–7 μg of total RNA with 100 pmol of T7-d(T)24 primer (University of Wisconsin-Madison DNA Synthesis Facility), 1× first-stand buffer (Invitrogen), 10 mM DTT (Invitrogen), 200 U SuperScript II RT (Invitrogen), and 500 μM dNTPs (GE Healthcare Life Sciences) for 1 h at 37°C. Second-strand synthesis was accomplished by incubation with 1× second-strand buffer, 200 μM dNTPs, 10 U DNA ligase (Invitrogen), 40 U DNA polymerase I (Invitrogen), and 2 U RNaseH (Invitrogen) for 2 h at 16°C followed by a 5-min incubation with 10 U T4 polymerase (Invitrogen). Double-stranded cDNA was purified using the GeneChip Sample Cleanup Module kit (Affymetrix, Santa Clara, CA). A biotin-labeled in vitro transcription reaction was performed using the cDNA template and the Affymetrix 3′-amplification IVT labeling kit. Briefly, double-stranded cDNA was incubated with 1× IVT labeling buffer, 1× labeling NTP mix, and 1× T7 labeling enzyme mix for 16 h at 37°C to yield cRNA. cRNA was purified using the GeneChip Sample Cleanup Module. cRNA yield was determined by absorbance at 260 nM. The cRNA was fragmented at 0.5 μg/μl final concentration in 1× fragmentation buffer (40 mM Tris-acetate pH 8.1, 100 mM potassium acetate, and 30 mM magnesium acetate) for 32 min at 94°C. The size range of the cRNA before (0.1 kb and longer) and after (35–200 base fragments) was checked by agarose gel electrophoresis.

Hybridization and washing.

The Affymetrix Rat GeneChips (Affymetrix) contains >31,000 probe sets representing 28,700 rat genes. A total of nine chips were hybridized using three biological triplicates (containing 3 ligaments per sample) at three time points (intact and days 3 and 7). Microarray hybridization was performed by the University of Wisconsin Biotechnology Gene Expression Center (Madison, WI). GeneChips (Affymetrix) were hybridized with biotin labeled cRNA in 300 μl in the presence of 1× hybridization buffer [100 mM MES, 1 M NaCl (Ambion, Austin, TX), 20 mM EDTA (Ambion), and 0.01% Tween-20 (Pierce Chemical)], 50 pM control oligonucleotide B2 (Affymetrix), 0.1 mg/ml herring sperm DNA (Promega, Madison, WI), 0.5 mg/ml acetylated BSA (Invitrogen), and 1× eukaryotic hybridization controls (bioB, bioC, bioD, and cre at 1.5, 5, 25, and 100 pM, respectively; Affymetrix) for 16 h at 45°C on a rotisserie at 60 rpm. Before application to the GeneChip, samples were heated at 95°C for 5 min, followed by incubation at 45°C for 5 min and spun at 14,000 g for 5 min. Following hybridization, the labeled samples were removed from the GeneChip, stored in the appropriate vial at −20°C, and immediately filled with nonstringent buffer A that contains 6× SSPE (Ambion), 0.01% Tween 20. All GeneChips were postprocessed using the automated Affymetrix GeneChip Fluidics Station 450.

The postprocessing protocol of the Rat 2.0 Genome GeneChip included the following: one wash of 10 cycles with 2 mixes/cycle of nonstringent buffer A at 25°C; one wash of 4 cycles comprised of 15 mixes/cycle with stringent buffer B (100 mM MES, 0.1 M NaCl, and 0.01% Tween 20) at 50°C; staining of the probe array for 10 min at 25°C in streptavidin-phycoerythrin (SAPE) solution [1× MES stain buffer (100 mM MES, 1 M NaCl, 0.05% Tween 20), 2 mg/ml acetylated BSA, and 10 μg/ml SAPE (Molecular Probes)]; poststaining consisting of 10 wash cycles of four mixes/cycle with nonstringent buffer A at 25°C; a second staining of the probe array for 10 min in antibody solution [1× MES stain buffer, 2 mg/ml acetylated BSA, 0.1 mg/ml Normal Goat IgG (Sigma-Aldrich, St. Louis, MO), 3 μg/ml biotinylated antibody (Vector Laboratories)]; a third stain of the probe array for 10 min in SAPE solution at 25°C; and a final wash of 15 cycles of 4 mixes/cycle with nonstringent buffer A at 30°C.

Probe array scan and data acquisition.

To quantify the fluorescent signal from each feature on the GeneChip, all GeneChips were scanned at a wavelength of 570 nm, using the GC3000 G7 scanner. Fluorescent signals corresponding to hybridization intensities were analyzed with the Affymetrix GCOS software.

Microarray analysis.

Quantification of GeneChips was performed by Affymetrix GeneChip Operating Software (AGCC version 2.0). Microarray data were log transformed [log (base 2)], and the log ratio between the intact ligament with day 3 or day 7 was calculated and expressed as fold change. A gene was considered up or downregulated if the day 3 or day 7 samples were at least twofold different than the intact controls. Additionally, data were analyzed using the two-sample Welch t-test to determine significance per gene. To minimize the false discovery rate, the Storey q-value method was performed, and significance was based on q < 0.01. Any q value ≤0.01 per gene was considered significant. Paired t-tests were used to evaluate differences between intact and day 3, intact and day 7, or day 3 and day 7. A critical value of 0.05 was considered as the criterion to select a significant fold change in gene expression between days. Data are presented as fold change of the means of day 3 and day 7 from the intact ligament values (fold change ± SD or fold change ± SE). To validate the expression patterns of microarray analysis with qPCR, the qPCR data was first normalized to the housekeeping gene β-actin. A fold change was then determined between the intact and day 3 samples or the intact and day 7 samples. Paired t-tests were performed to determine qPCR significance between days.

Functional annotation.

Microarray data were functionally grouped using the Database for Annotation, Visualization, and Integrated Discovery (D.A.V.I.D.;, integrating functional genomic annotations with intuitive graphical summaries. With the use of the Affymetrix probes as identifiers, the samples were subjected to functional annotation clustering based on function and gene ontology. The categories of interest were further based on the enrichment score of the gene group. The enrichment score, determined by the minus log transformation of the geometric mean P values from all enriched annotation terms, ranks the biological significance of a gene cluster based on the scores of all the annotation terms.

RT and qPCR.

RNA was DNAase treated and reverse transcribed into cDNA using the Superscript III first strand kit (Invitrogen). RNA was added to 50 ng/μl random hexamer, annealing buffer, and RNAse/DNase free water, followed by 2× first-strand reaction mix and RNase enzyme mix. Reverse transcribed samples were diluted 2.5 ng/μl.

To validate microarray genes of interest, quantitative PCR was performed on ligament samples. Primers for collagen type I (Col1a1, Col1a2), type III (Col3a1), type V (Col4a3), type VI (Col6a1, Col6a3), type XI (Col11a1), type XXVII (Col27a1), collagen triple helix repeat containing-1 (Ctrhc1), periostin (Postn), serine protease 35 (Prss35), Mmp3, Mmp11, Mmp12, Mmp13, Mmp14, Mmp19, Timp1, and β-actin primers were utilized (Table 1). Primer sets were either obtained from published reports (11, 18, 26, 33, 36, 64, 66, 74, 77) or designed from sequences available from Genbank from Blast (Basic Local Alignment Search Tool;

View this table:
Table 1.

List of primer sequences for quantitative PCR, the annealing temperature for each primer, and where the primers were obtained

qPCR was performed using a Bio-Rad thermocycler (Bio-Rad, Hercules, CA). All reactions were carried out using Bio-Rad Sybr Green Supermix (Bio-Rad). Assessment of genes was performed by including negative controls in each qPCR experiment, examining PCR melt curves to ensure specificity, and verifying the primer efficiency was within 90–110% based on the standard curve. The cDNA of the target gene was also normalized to the housekeeping gene β-actin.


Immunohistochemistry was performed on frozen sections using mouse monoclonal or rabbit polyclonal antibodies. Sections were fixed with acetone, exposed to 3% hydrogen peroxide to eliminate endogenous peroxidase activity, and blocked using Background Buster (Innovex Biosciences, Richmond, CA). Sections were then incubated for 2 h with rabbit polyclonal periostin (1:1,000; Biovendor, Candler, NC) and mouse monoclonal CTHRC1 (1:500; ab72527; Abcam, Cambridge, MA). After primary antibody exposure, sections were treated with biotin and streptavidin conjugated to horseradish peroxidase using the Stat Q staining kit (Innovex Biosciences). The bound antibody complex was then visualized using DAB. Stained sections were dehydrated, cleared, coverslipped, and viewed under light microscopy. Negative controls omitting the primary antibody were included with each experiment. Positive controls of gut or spleen were also included.


Microarray results.

Microarray analysis was conducted to identify differences in the ligament transcriptome between day 3 and day 7 postinjury with the uninjured ligament. A total of 31,000 probe sequences were analyzed. Overall, 2,652 genes from day 3 and/or day 7 were at least twofold greater than the intact controls and achieved a q value of 0.01 or less (Table 2). Of those genes, 2,046 genes were significantly regulated at day 7. Of the day 7 significant targets, 1,161 were downregulated and 883 were upregulated. A total of 1,991 genes were different at day 3 from the intact ligament, with 1,050 upregulated targets and 941 downregulated. Those day 3 and/or day 7 genes significantly different from the intact were then subjected to D.A.V.I.D. and clustered via functionality, resulting in 238 groups. Table 3 lists functional clusters of interest. Seventeen clusters of interest were identified with an enrichment score ranging from 3.15 to 0.15. Functional groups include genes involved in ECM control, metalloenzyme/peptidase/metallopeptidase/proteolysis regulation, immune cell activation/proliferation/differentiation/regulation, proteoglycan/glycoprotein binding, regeneration/limb morphogenesis/wound response, endocytosis/phagocytosis, angiogenesis, cell chemotaxis, and innate response regulation. To further identify relevant gene clusters based on day of injury, genes were selected only if they were significant on one day but not the other and subjected to D.A.V.I.D. analysis. A total of 600 and 655 genes were exclusively regulated at day 3 and day 7, respectively. The remaining 1,391 genes were modulated at both day 3 and 7. The top five functional groups from each day are listed in Tables 4 and 5. Functional annotation grouping via D.A.V.I.D. indicated that many genes are involved in immune regulation at day 3 (Table 4) whereas the day 7 exclusive genes are primarily involved in ECM regulation and angiogenesis (Table 5). Of the 600 genes modulated at day 3, 434 were upregulated and 166 were downregulated compared with the uninjured ligament (Fig. 1). In comparison, 376 genes were upregulated and 279 genes were downregulated in the day 7 group (Fig. 2). Genes of interest are listed along the Venn diagrams and include inflammatory factors, growth factors, and ECM components, including collagens and MMPs.

View this table:
Table 2.

Number of genes significantly regulated during early wound healing

View this table:
Table 3.

Functional Annotation Clusters of interest obtained from D.A.V.I.D. analysis

View this table:
Table 4.

Functional annotation clustering of the top 5 day 3-specific genes obtained from D.A.V.I.D. analysis

View this table:
Table 5.

Functional annotation clustering of the top 5 day 7-specific genes obtained from D.A.V.I.D. analysis

Fig. 2.

Venn diagrams of specific extracellular matrix (ECM)-related genes up- or downregulated at days 3 and/or day 7 postinjury. Genes of interest involved in specific immune response or ECM regulation at day 3 and 7, respectively, are listed. Genes modulated at both day 3 and 7 are also listed and primarily include factors involved in ECM regulation. Factors analyzed further via quantitative PCR are indicated in bold.

Collagen expression.

Twenty-one collagen isoforms were detected within the healing ligament via microarray. Ten of the 21 isoforms were significantly changed after injury (Fig. 3). Types I and III collagens did not experience a significant fold change but were the most abundant collagens identified within the ligament regardless of injury. Types IVa3bp, Va1, Va2, VIa2, VIIa1, XVa1, XVIa1, and XVIIIa1 collagens were also constitutively expressed and did not exhibit a significant fold change. In contrast, type IVa5, Va3, VIa3, XIa1, XIa2, XIIa1, XIIIa1, XXIIa1, XXIVa1, and XXVIIa1 were all differentially expressed during healing. At both day 3 and 7, types IVa5 and XIa2 were downregulated while type Va3 was upregulated. The other differentially expressed collagens were only regulated on day 3 (collagen types XIa1 and XXIVa1) or day 7 (collagen types VIa3, XXVIIa1, XIIa1, XIIIa1, and XXIIa1).

Fig. 3.

Heatmap showing changes in genes of interest. Data indicate fold change of day 3 or day 7 genes from the intact normal ligament. Green and red bars indicate an upregulation and downregulation of genes, respectively.

MMP expression.

Ten MMPs were detected in the ligament via microarray. Mmp2 was the primary MMP expressed within the ligament albeit in a constitutive manner (Fig. 3). Mmp3 was the second highest MMP expressed in the intact ligament, followed by Mmp14, 16, 13, 12, 11, 28, 19, and 7. Mmp7, 16, and 28 did not exhibit a twofold difference after injury. In contrast, Mmp3, 11, 12, 13, 14, and 19 all exhibited at least a twofold change between groups. Additionally, Timp1 and Timp3 were significantly up and downregulated at day 3 and day 7 postinjury, respectively.

Validation of selected genes.

Several of the ECM and collagen genes of interest identified by microarray analysis were further validated via qPCR (Fig. 4). The genes at day 3 and day 7 were first normalized to the housekeeping gene β-actin. The day 3 and day 7 values were then compared (via fold change) to the intact ligament. Similar to the microarray results, type I collagen was not changed over day 3 and day 7 compared with the intact ligament. Type III collagen increased at day 7 compared with day 3 and the intact ligament, but the fold change was under two. Type VI alpha 1 and XI alpha 1 were also increased with injury, but the changes were not over twofold. In contrast, collagen type V alpha 3, VI alpha 1, and XXVII alpha 1 isoforms were significantly increased over the intact control. Expression of type VI alpha 1 and XXVII alpha 1 was highest at day 7 whereas type V alpha 3 was highest at day 3 and then decreased.

Fig. 4.

qPCR results of the genes of interest based on microarray results. A: qPCR analysis of genes found significantly regulated during microarray analysis or typically found in the ligament (as in the case for types I and III collagen). B: Mmps most significantly differentiated (Mmp11 and 13). C: other Mmps and Timp qPCR results. Data are expressed as fold change ± SE. Data were normalized with the housekeeping gene β-actin and compared with the normal intact ligament (and expressed as fold change).

Similar to the microarray results, Mmp3, 11, 12, 13, 14, and 19 all exhibited a twofold increase in gene expression after injury. Mmp11 and 13 demonstrated the greatest change after injury. All Mmps tested were greatest in expression at day 7 excluding Mmp13 and Mmp3, which decreased. Timp1 was also highest at day 7.

Cthrc1, Prss35, and Postn expression.

Microarray analysis detected three novel genes, Postn, Cthrc1, and Prss35, which were significantly upregulated after injury. Periostin (Postn), collagen triple helix repeat containing-1 (Cthrc1), and serine protease 35 (Prss35) were all significantly increased 7 days postinjury (Fig. 4). qPCR analysis verified these increases. Protein localization of periostin and CTHRC1 was further analyzed via immunohistochemistry. Few regions were positive for either factor in the intact ligament (Fig. 5, A and D). Injury stimulated an increase in both proteins at day 3 (Fig. 5, B and E). Levels were further increased at day 7 (Fig. 5, C and F).

Fig. 5.

Periostin (Postn), collagen triple helix repeat containing 1(Cthrc1), and serine protease 35 (Prss35) expression/localization during healing. qPCR and microarray comparison of Postn (top left), Cthrc1 (middle right), and Prss35 (bottom left) at 3 and 7 days postinjury in relation to the normal intact MCL. Immunohistochemistry (right) of Postn (A, B, and C) and CTHRC1 (D, E, and F) of the intact (A and D), day 3 (B and E), and day 7 (C and F) MCL. Staining for both periostin and CTHRC1 were low in the intact MCL (A and D, respectively). After injury, both periostin and CTHRC1 increased. qPCR and microarray data are expressed as fold change ± SD. qPCR data were normalized with the housekeeping gene β-actin and compared with the normal intact ligament (and expressed as fold change). Microarray data are compared with the normal intact ligament (and expressed as fold change).


The aim of this study was to profile transcriptional changes and identify candidate genes to regulate inflammation, ECM turnover, and scar formation during ligament healing. Our previous research (8) demonstrated significant vascular and cellular upregulation within the first week of ligament healing. These results indicated a significant increase in macrophages on days 1–5 and angiogenesis on days 7–11 postinjury. Additionally, the original wound size, based on size of granulation tissue, became increasingly larger (i.e., creeping substitution) at these times. We believe these overlapping immunogenic, cellular, angiogenic, and wound size changes control downstream scar formation. Based on the previous study, days 3 and 7 postinjured ligaments were chosen as time points to determine genetic changes from the intact ligament. Our current results identified a number of Mmps, Timps, and collagens modulated during healing. Specifically, minor collagens such as types XI, XXII, XXIV, and XXVII were identified and differentially expressed within the healing MCL. Unique genes, such as collagen triple helix repeat containing 1, periostin, and serine protease 35 were also identified in our healing model. Of the chosen factors, collagen types Va3, Via2, XIVa1, Cthrc1, Mmp3, Mmp11, Mmp12, Mmp13, Mmp14, and Mmp19, Timp3, Postn, and Prss35 were highly differentiated from the intact ligament, exhibiting a fourfold or greater change. To our knowledge, this is the first study to comprehensively identify and characterize genes involved during early MCL healing. Based on D.A.V.I.D. results, day 3 healing primarily includes genes associated with immune regulation. Our previous study (8) defined the spatial and temporal factors of normal ligament healing and demonstrated a strong upregulation of macrophages, elucidating the enhanced day 3 inflammatory factors in the current study. D.A.V.I.D. analysis further associated the day 7-specific genes to cell adhesion, ECM, and angiogenic functions that coincide with normal granulation tissue formation involving the upregulation of blood vessels, fibroblasts, and collagen turnover. Additional functional group analysis identified MMPs and collagens specific to the ligament. Similar inflammatory factors, MMPs, and collagens were also modulated in the early healing Achilles tendon (17). Recognition of the ligament MMPs and collagens provides a novel report of ligament injury-induced gene expression and sheds light on the possible players involved in the typical functional compromise (creeping substitution and scar formation) associated with normal ligament healing.

Two groups of collagens encoded by 42 genes to generate 29 types currently exist, including the fibril-forming (collagen types I, II, III, V, XI, XXIV, and XXVII) collagens or the nonfibril-forming collagens that are further subgrouped into the fibril associated with interrupted triple helices (FACIT; collagen types IX, XII, XIV, XIX, XX, and XXI), network-forming (collagen types VIII and X), microfibrillar (collagen type VI), anchoring fibril (collagen type VII), transmembrane (collagen types XIII and XVII), basement membrane (collagen type IV), multiplexin (collagen types XV, XVI, and XVIII), or other (type XXII, XXIII, XXV, and XXVI) collagens. In the current study, the fibril-forming (collagen type V, XI, and XXVII), FACIT (type XII collagen), microfibrillar (type VI collagen), transmembrane (type XIII collagen), and basement membrane (type IV collagen) collagens were differentially expressed after injury. The fibril-forming types I and III were the predominant collagens present but were expressed in a constitutive manner.

The fibrillar collagens were the most abundantly expressed in the current MCL model. Types I, III, and V collagens form heterotypic collagen fibrils, the fundamental functional unit of the ligament providing structural support. Type I collagen is the primary ligamentous collagen, consisting of two α1-chains and one α2-chain synthesized in a 2:1 ratio. Type I collagen was constitutively expressed within our healing model, regardless of injury. These results compare to our recent spatial and temporal study describing type I procollagen protein levels in the normal healing ligament (unpublished observations). Type I procollagen did not change the first week of healing but steadily decreased thereafter. These combined results support that type I collagen levels of the intact ligament are not different than the healing ligament within the first week of injury.

Fibrillar type III collagen intersperses within the type I collagen matrix. Type III collagen is commonly used as a marker for scar formation in the ligament. The microarray results showed no significant increase in type III collagen after injury, but the qPCR data indicated a slight increase at day 7. Previous results from our normal healing experiment localizing type III collagen via immunohistochemistry likewise demonstrated no change at day 3 but evidenced an increase at day 7 (unpublished observations). As scar formation progressed, type I collagen decreased and type III collagen increased (unpublished observations).

In contrast to the type I and III collagen microarray results, type V collagen was significantly upregulated at days 3 and 7 postinjury. Collagen type V is a widely distributed, low abundance fibrillar collagen believed to stimulate fibril nucleation and regulate fibril size and organization. Type V collagen mutations are linked to Ehlers-Danlos syndrome, a connective tissue disorder characterized by hyper-expandable skin (61). Levels also increase with hypertrophic scar formation, aging patellar tendons, and degenerating tendons (9, 12). Type V collagen abnormalities are accompanied by unusually large type I collagen diameters and increased fibril disorganization (2). More recently, type V collagen has been identified as an autoantigen of chronic inflammation/fibroproliferation processes and is believed to serve as a candidate marker for tendon pathology (7, 38). Based on these reports and the apparent importance in regulating type I collagen organization, type V collagen may play a critical role in scar formation as well as replacement fibrosis. Type V collagen may also play a fundamental role as an identifier of normal ligament injury.

To our knowledge, this is the first report to identify the fibril-forming types XI, XXIV, and XXVII collagens as well as type XXII collagen within the healing ligament. Type XI collagen is primarily known for its localization in cartilaginous tissues. However, in the current study, we verified the presence and the reduced expression of type XI collagen in the healing ligament. Type XI collagen is assembled as a heterotrimer consisting of α1 (XI)-, α2 (XI)-, and α3 (XI)-chains. Type XI collagen mutants develop Stickler syndrome, a connective tissue disorder resulting in beaded bundles of irregular collagen diameter causing premature arthritis and craniofacial disorders (1, 6, 57). Although its function is largely unknown, these results suggest that type XI collagen functions to regulate fibril diameters. Even less is known about types XXII and XXIV, but both participate in controlling bone and cartilage physiological processes such as collagen fibrillogenesis. Type XXIV is found primarily in the developing eye and skeleton and, to a lesser degree, skin and tendon during periods of accentuated type I collagen fibrillogenesis (20). Unlike types XI, XXII, and XXIV, type XXVII collagen may function in cartilage mineralization, cellular apoptosis, blood vessel invasion, matrix degradation, and type I collagen endochondral bone formation (15). Type XXVIIa1 collagen is primarily, but transiently, expressed in developing cartilage and dermis, the eye, and some epithelial layers during a narrow window of development in several tissues (48). The temporary expression of collagen type XXVII may function in cell signaling or as a supporting environment for invading blood vessels rather than structural stability. Type XXVII collagen upregulation in the current model may also support the transient blood vessel infiltration and matrix turnover during ligament healing. Taken together, these minor collagens may play a role in collagen fibrillogenesis and/or structural support to the MCL.

The fibrillar collagens provide structural support for the matrix, but the nonfibrillar collagens regulate anchoring and organization of the ECM. Type IV collagen is the primary component of basement membranes, forming a collagen, proteoglycan, glycoprotein meshwork. Type IV collagen forms the basement membranes around blood vessels. Angiogenesis increases substantially after injury, peaking around days 7–11 postinjury (8) but collagen type IV expression decreased in the current healing model. Type IV collagen is not required for basement membrane formation but rather for maintenance by mediating cell adhesion to the basement membrane. Injury disrupts normal ligament matrix, and the initial reduction in type IV collagen may permit cell migration to the damaged region, followed by reparative and replacement fibrosis via fibrillar collagen deposition and the eventual deposition of type IV collagen.

The microfibrillar type VI collagen contributes to a network of beaded filaments interwoven with other collagen fibrils in all connective tissues excluding bone. Type VI collagen acts as an adhesive between ECM components and the basement membrane (10). More recently, type VI has been reported to induce myofibroblast differentiation (40). Dermal wound healing studies reported type VI collagen from day 4 to 7 postinjury (59). The current results demonstrated a similar increase in type VI expression at day 7 postinjury, which coincides with myofibroblast upregulation in the healing ligament (unpublished observations). The later increase in type VI collagen may thereby maintain the number of myofibroblasts, which may contribute to creeping substitution and scar formation.

Only one FACIT collagen (type XII collagen) was differentially expressed in our wound healing model. FACIT collagens do not form supramolecular aggregates but link fibrils to other ECM proteins and bind collagen fibrils to regulate fibril growth, spacing, and diameter. Type XII collagen colocalizes with type I collagen and binds other ECM proteins, such as tenascin, to render the tissue more resistant to mechanical stress (42). Type XII collagen is also associated with enhancing the fluidity of collagen fibrils, presumably by reducing the friction between sliding fibrils (42). Additionally, type XII collagen mutants are associated with augmented female anterior cruciate ligament injuries (49). Based on these results, the day 7 increase in type XII collagen likewise provides a structural support between ECM and fibrillar collagens within the MCL.

The transmembrane type XIII collagen exists both as a membrane-bound and soluble molecule. Although its exact role is unknown, it plays an important role in certain adhesive interactions necessary for normal development. Research with type XIII collagen mutants mice developed lymphomas, abnormal basement membranes, and enhanced inflammatory gene expression (68). The reduced expression of type XIII collagen at day 7 in the current study likely reflects the disruption of collagen organization and hence cell-cell adhesion after injury.

Collagen triple helix repeat containing-1 (Cthrc1) is a gene product able to control collagen deposition and modulate transforming growth factor-β (TGF-β; Refs. 23, 51). TGF-β modulates transdifferentiation of fibroblasts to myofibroblasts by promoting expression of smooth muscle actin and subsequent scar formation (29, 34, 46). Cthrc1 overexpression stimulates faster cell migration to the wound and reduces TGF-β signaling, resulting in decreased wound formation, smooth muscle cell proliferation, and collagen deposition (24). Overexpression of Cthrc1 also reduces type Ia1, a2, and type IIIa1 collagen mRNA (24, 51). Based on these results, an increase of CTHRC1 may limit scar formation in the injured ligament by modulating TGF-β signaling. This concept, however, requires additional ligament healing studies using overexpressing and knockout Cthrc1 animals.

The MMPs and tissue inhibitors of MMP (TIMPs) are key enzymes involved in ECM remodeling during ligament healing. MMPs degrade numerous substances in the ECM, including collagen and fibrin, to expedite cell migration, deposit new ECM, and develop new tissue. MMPs also regulate the ECM structure by activating and releasing cytokines and growth factors from the ECM, sending cues to influence cell behavior and action. Four subsets of MMPs exist, including the collagenases, gelatinases, stromelysins, and membrane-type MMPS. In the present study, the collagenase (Mmp13), the stromelysins (Mmp3, Mmp12, and Mmp11), and membrane-bound (Mmp14 and Mmp19) collagen were expressed differentially in the healing ligament. The gelatinase Mmp2 was also expressed strongly in the ligament but in a constitutive manner.

Collagenases are proteases capable of degrading native triple helix fibrillar collagens types I, II, and III. In the current study, only the collagenase Mmp13 significantly increased in expression after injury. This Mmp was also the most significantly influenced by injury. Previous reports (21, 54, 55, 63, 71) demonstrated increased expression of MMP-13 in areas of excessive collagen destruction including osteoarthritis, rheumatoid arthritis, and periodontal inflammation. Message levels also significantly increased in the tendon injury model (5, 65). Cutaneous wound healing in Mmp13 knockout mice resulted in delayed wound closure and reepithelialization, decreased vascularization in granulation tissue, and decreased wound contraction and myofibroblast formation (14). Given these results, MMP-13 appears as the primary MMP involved in fibrillar collagen turnover, cell migration, angiogenesis, and MCL remodeling. Although not tested, MMP-13 may also contribute to the decrease in fibrillar collagens types XI, XXIV, and XXVII. Control of MMP-13-induced collagen degradation may therefore control creeping substitution and subsequent scar formation.

Stromelysins contribute a broad role in tissue development, remodeling, and wound healing. This class is also associated with collagen degradation of types IV, V, IX, and X. The stromelysin MMP-11 is expressed in normal and pathological remodeling processes including inflammation, wound healing, postpartum involution, and tumor proliferation (13, 58, 60, 67, 75). MMP-11 is primarily expressed in fibroblasts and does not degrade the typical ECM substrates (including collagens, gelatins, laminin, and fibronectin) but is involved in tissue remodeling, antiapoptotic events, and inflammation to favor cell survival and tumor progression. Collectively, these results suggest the significant increase in MMP-11 by the day 7 healing ligament may control inflammation and collagen degradation, before angiogenesis and after the inflammatory response. The role of MMP-11 in tumor progression suggests its possible role in modulating ligament creeping substitution.

The current ligament healing model demonstrates an increase in stromelysin MMP-12 starting at day 3 and continuing to day 7. MMP-12 is primarily produced by macrophages and is essential in stimulating macrophage migration. It also functions to digest elastin and other ECM molecules. Expression of MMP-12 is upregulated in several lung diseases (22, 27, 41, 76). A study (32) examining Mmp12 knockout mice reported an increase in Mmp2 and IL-13/IL-4R-induced Mmp13 expression after Mmp12 deficiency. Based on these results, control of MMP-12 may control phenotype-specific macrophage infiltration into the healing ligament and thereby modulate creeping substitution and subsequent scar formation.

MMP-3 degrades a wide range of ECM proteins including types IV, V, XI, and X collagens, proteoglycans, fibronectin, and laminin. Studies (47) with MMP-3 deficient animals suggest this protease functions to detach cells from the underlying basement membrane to initiate cell migration. Additional MMP-3 knockout studies (73) using the lung injury model resulted in reduced lung injury by inhibiting leukocyte recruitment. These studies suggest that the transient increase of MMP-3 in the day 3 ligament recruits cells to the wound site, and promotes inflammation.

Unlike other MMPs, membrane type-MMPs are not secreted into the ECM. Rather, they reside on cell membranes and function by binding to and activating other MMPs. MMP-14 is associated with stimulating MMP-2 and MMP-9 (70). In the current study, Mmp9 was not identified in the healing ligament, whereas Mmp2 was expressed. The increased Mmp14 expression observed at days 3 and 7 within the healing ligament could therefore modulate MMP-2 activity within the healing ligament model. Additionally, MMP-14 is collagenolytic towards type I, II, and III collagens (62, 70) and may likewise mediate collagen breakdown.

Membrane bound MMP-19 functions to degrade ECM molecules including the basement membrane and type IV collagen. This enzyme does not cleave type I collagen. MMP-19 has been localized to vascular smooth muscle cells, endothelial cells, fibroblasts of granulation tissue, monocytes, and macrophages and is upregulated during inflammatory conditions such as arthritis and multiple sclerosis (3, 4, 31, 53). MMP-19-deficient mice demonstrated enhanced tumor invasion, reduced CD8+ T-cell activation and distribution, and a diminished inflammatory reaction determined by lowered cytokine expression, decreased influx of inflammatory cells, lessened cell proliferation, and increased angiogenesis (3, 19). Additional reports described the downregulation of MMP-19 in invasive carcinomas, suggesting its role in tissue homeostasis. Based on these reports, the increase in Mmp19 (likely released by macrophages and fibroblasts) may recruit inflammatory cells and modulate degradation of ECM during inflammation.

TIMPs are natural inhibitors of MMPs. Currently, four members of the TIMP family have been identified and include TIMP-1, TIMP-2, TIMP-3, and TIMP-4. Timp1 was identified in the current healing model and inhibits most MMP activity excluding MMP-14 and MMP-2. Timp3 was also differentially regulated during ligament healing. TIMP-3 is bound to ECM and is not a freely soluble protein but is able to inhibit other molecules outside of MMPs including a disintegrin and metalloproteinase domain (ADAM) and ADAMS with thrombospondin domains (ADAM-TS). Timp2 and Timp4 were not identified in the current healing model.

The interaction between collagens and MMPs during ligament healing is not fully understood but plays a significant role in matrix turnover. In conjunction with previous reports and the current study, types I and III collagens are eventually degraded by collagenases that breakdown the helical structure. In the ligament, MMP-13 primarily controls fibrillar collagen turnover, but membrane-bound MMP-14 may also contribute. Collagen types IV and V are not targeted by collagenases but by gelatinases or stromelysins (30, 69). Gelatinases were not differentially regulated in the current study but the stromelysins were. Thus type IV and V collagens appear to be modulated by MMP-3 and -12 in the ligament model. MMP-3 and -12 also stimulate inflammatory cell migration and detachment while degrading collagen. Type VI collagen is not subject to the same breakdown by the MMPs and is degraded by serine proteases produced primarily by neutrophils and mast cells during inflammation. This suggests type VI collagen levels remain in check in the presence of inflammatory cells but increase as cell number decreases. Later, fibroblast-stimulated MMP-11 may ameliorate the macrophage MMP-induced degradation and initiate collagen deposition. Further study is required to elucidate the functional interactions among the MMPs, TIMPs, CTHRC1, and collagens.

Periostin is a gene of interest based on microarray data and current literature. Prior studies of periodontal ligament (16, 16), remodeling bone (16, 16), carcinomas (50), and myocardial infarcted heart (39, 43, 43) reported an increase in periostin and subsequent increase in deposition and attachment of collagen. Periostin knockout mice exhibit reduced stiffness and diminished collagen formation after myocardial infarction, suggesting improved cardiac function (16, 44, 45). Although these results indicate that periostin participates in collagen synthesis during wound healing, successful ligament healing also requires the spatial and temporal aspects of collagen alignment and production of the appropriate type of collagen. Previous research from our laboratory demonstrates a significant increase in type III collagen and a concomitant decrease in type I procollagen during the first week of healing. These immunohistochemistry results suggest that periostin colocalizes with type III collagen of the healing ligament. Current results indicate that periostin increases during the first week of healing to stimulate fibroblast proliferation and regulate collagen synthesis. Periostin knockout studies are currently underway to determine its role during wound healing.

This study is the first to report the presence of Prss35 within the healing ligament. A paucity of information is known about this extracellular serine protease, but its presence has been reported in the mouse ovary (37, 72). The reported spatial and temporal regulation of PRSS 35 within the mouse ovary suggests it may be involved in tissue remodeling that occurs with ovulation (37, 72). The significant upregulation of Prss35 within the healing MCL at day 7 also suggests its possible contribution to ligament remodeling. More research needs performed to determine the role of this protease.

When interpreting results, several things must be borne in mind. First, the study investigates healing in a rat MCL healing after surgical transection. Second, within this model, protein synthesis and functional expression may not directly correlate with the mRNA levels reported in the study. Third, periostin and CTHRC-1 immunohistochemistry was also performed to verify protein expression of these novel factors. However, many additional factors reported herein would benefit from further protein analysis. In addition, microarray results were not always consistent with qPCR results. Several causes may lead to such discrepancies. For instance, specific gene expression could be defined as absent via microarray but may be detected by the more sensitive qPCR technology. Further, the intrinsic variability in the affinity of labeled samples to the microarray probes could influence the fold change in gene activity between experimental groups. In our study, the entire ligament, rather than the injured region within the ligament, was assayed. The inclusion of the total ligament diluted expression particular genes regulated specifically within the healing region. Lastly, microarray results are influenced by array production, probe labeling, hybridization conditions, and image analysis (52), which may also explain any discrepancies with the qPCR analysis.

Our experiments focused on the first week of ligament healing during inflammation and proliferation, rather than the subsequent remodeling. We believe scar formation is a downstream result of various immune cell-mediated responses to injury. If so, identifying and modulating factors involved in creeping substitution, macrophage infiltration, and angiogenesis would limit scar formation and reduce healing time. A similarly robust analysis of later time points during healing would also supply a wealth of information on the formation and remodeling of the neo-ligament.

In summary, differences in the genetic expression of MMPs, collagens, and the novel genes Postn, Cthrc1, and Prss35 were observed in the early healing ligament. To our knowledge, these results provide the first comprehensive assessment of ECM expression in the MCL and will lead to a better understanding of the molecular mechanisms of creeping substitution and scar formation. Identification of these factors during early healing is novel, but their specific functional roles and interactions require further elucidation.


Financial support was provided by National Institute of Arthritis and Musculoskeletal and Skin Diseases Grants AR049266 and AR059916.


No conflicts of interest, financial or otherwise, are declared by the authors.


We acknowledge the petit autem pro rat sacrificium ultimum, the technical and computer support of Sandra Splinter BonDurant and Jean-Yves Sgro, heatmap development of Kayt Frisch, and statistical work of Alejandro Munoz del Rio.


  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
  46. 46.
  47. 47.
  48. 48.
  49. 49.
  50. 50.
  51. 51.
  52. 52.
  53. 53.
  54. 54.
  55. 55.
  56. 56.
  57. 57.
  58. 58.
  59. 59.
  60. 60.
  61. 61.
  62. 62.
  63. 63.
  64. 64.
  65. 65.
  66. 66.
  67. 67.
  68. 68.
  69. 69.
  70. 70.
  71. 71.
  72. 72.
  73. 73.
  74. 74.
  75. 75.
  76. 76.
  77. 77.
View Abstract