Fibroblasts from phenotypically normal palmar fascia exhibit molecular profiles highly similar to fibroblasts from active disease in Dupuytren's Contracture

Background Dupuytren's contracture (DC) is a fibroproliferative disorder characterized by the progressive development of a scar-like collagen-rich cord that affects the palmar fascia of the hand and leads to digital flexion contractures. DC is most commonly treated by surgical resection of the diseased tissue, but has a high reported recurrence rate ranging from 27% to 80%. We sought to determine if the transcriptomic profiles of fibroblasts derived from DC-affected palmar fascia, adjacent phenotypically normal palmar fascia, and non-DC palmar fascial tissues might provide mechanistic clues to understanding the puzzle of disease predisposition and recurrence in DC. Methods To achieve this, total RNA was obtained from fibroblasts derived from primary DC-affected palmar fascia, patient-matched unaffected palmar fascia, and palmar fascia from non-DC patients undergoing carpal tunnel release (6 patients in each group). These cells were grown on a type-1 collagen substrate (to better mimic their in vivo environments). Microarray analyses were subsequently performed using Illumina BeadChip arrays to compare the transcriptomic profiles of these three cell populations. Data were analyzed using Significance Analysis of Microarrays (SAM v3.02), hierarchical clustering, concordance mapping and Venn diagram. Results We found that the transcriptomic profiles of DC-disease fibroblasts and fibroblasts from unaffected fascia of DC patients exhibited a much greater overlap than fibroblasts derived from the palmar fascia of patients undergoing carpal tunnel release. Quantitative real time RT-PCR confirmed the differential expression of select genes validating the microarray data analyses. These data are consistent with the hypothesis that predisposition and recurrence in DC may stem, at least in part, from intrinsic similarities in the basal gene expression of diseased and phenotypically unaffected palmar fascia fibroblasts. These data also demonstrate that a collagen-rich environment differentially alters gene expression in these cells. In addition, Ingenuity pathway analysis of the specific biological pathways that differentiate DC-derived cells from carpal tunnel-derived cells has identified the potential involvement of microRNAs in this fibroproliferative disorder. Conclusions These data show that the transcriptomic profiles of DC-disease fibroblasts and fibroblasts from unaffected palmar fascia in DC patients are highly similar, and differ significantly from the transcriptomic profiles of fibroblasts from the palmar fascia of patients undergoing carpal tunnel release.


Background
Dupuytren's contracture (DC) is characterized by abnormal thickening of palmar fascia into collagen-rich cords that cause the fingers to bend and curl into a flexed and contracted state [1]. Although this disease can occur in both sexes, it is more common in men of Northern European descent [2][3][4] and typically presents in the 4th to 6th decade of life. DC has been reported to behave as a heritable genetic disorder, with evidence that it arises (in at least some cases) from an autosomal dominant gene on chromosome 16 with variable penetrance [3]. Lifestyle factors including smoking or heavy drinking [5,6], and heavy manual labor and hand trauma, have also been linked to development of DC [7,8], as have diabetes, epilepsy and hypercholesterolemia [9][10][11].
Treatment of DC remains problematic. A variety of non-surgical interventions, including injection of steroids [12] or gamma-interferon [13], use of creams based on vitamin E [14], dimethyl sulphoxides [15], and ultrasound therapy [16] etc. yield limited benefits. Recently, direct injection of clostridial collagenase has been evaluated with some promising results [17][18][19]. However, surgical excision of the involved contracted tissue remains the mainstay of therapy, supplemented with post-operative splinting and physical therapy [20,21]. Alternative therapeutic approaches remain desirable since surgery carries significant risks, including damage to the digital nerves and blood vessels, damage to the underlying flexor tendons, and wound healing failure with the possibility of skin necrosis.
A particularly vexing feature of DC is its propensity for recurrence despite the appearance of successful initial treatment. Many patients eventually require multiple surgeries with a cumulative risk of morbidity. However, it remains unclear what factors are responsible for recurrence of the disease. Apart from the possibility of a genetic predisposition, it has been hypothesized that undetected residual foci of incipient disease are present in the otherwise normal appearing and uninvolved palmar fascia, and that these cells represent sites of disease recurrence.
We have previously investigated the transcriptomic differences between fibroblasts derived from diseased DC cords versus fibroblasts from phenotypically normal palmar fascia in patients undergoing carpal tunnel (CT) release. These studies demonstrated intrinsic differences in gene expression between these cell populations that persisted even after propagation under cell culture. We have now extended these studies to include fibroblasts from macroscopically uninvolved (i.e. phenotypically normal) palmar fascia surgically removed from patients with DC. These fibroblasts, as well as fibroblasts from diseased DC cords and control carpal tunnel fibroblasts, were cultivated in cell culture on a type-1 collagen substrate to better approximate the in vivo collagen-enriched environment that these cells experience. The transcriptomic signatures of these three cell types were then compared to answer the question: do fibroblasts from phenotypically normal palmar fascia in DC more closely resemble their counterparts in phenotypically dissimilar DC cords, or cells from phenotypically similar carpal tunnel fascia?

Clinical specimens
Dupuytren's contracture (DC) cord samples and small samples of phenotypically normal palmar fascia tissue (PF) were surgically explanted at the Hand and Upper Limb Centre at St Joseph's Health Care (SJHC), London, ON, Canada. We also obtained phenotypically normal palmar fascia from patients undergoing carpal tunnel release (CT) in London, ON, Canada, and at the Allegheny General Hospital, Pittsburgh, PA. The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki. All specimens were collected under Institutional Review Board approval. Six patient samples of each tissue type were used to derive fibroblasts for use in these studies.

Primary cell culture
Primary cultures of fibroblasts were purified from the surgically resected DC cord and matching specimens of normal appearing palmar fascia (PF), and from specimens of normal palmar fascia of patients undergoing carpal tunnel surgery (CT) as previously described [22,23]. The cultures were maintained in α-MEM-medium supplemented with 10% fetal bovine serum (FBS, Invitrogen Corporation, Carlsbad, CA) and 1% antibiotic-antimycotic solution (Sigma-Aldrich, St Louis, MO). All cultures were harvested prior to the sixth passage, with no changes in cell morphology observed during the in vitro expansion protocol. For the present study we used 6 primary DC and patient-matched PF cell cultures as well as primary fibroblasts purified and expanded from 6 patients undergoing surgical carpal tunnel release. Primary cultures of fibroblasts were grown on collagen monolayers as described previously [24]. In brief, collagen was coated on 6-well tissue culture plates with each well containing 880 μl of rat tail (type-1) collagen and 200 μl of the neutralization solution (2 parts 0.34N NaOH and 3 parts 10x Waymouth media) to a final concentration of 1.9 mg/ml. Following collagen polymerization, primary cultures were grown in α-MEM-medium containing 10% FBS and 1% antibioticantimycotic solution and the cultures were left undisturbed at 37°C at 5% CO 2 for 72 hours.

Total RNA extraction
Each well containing cells grown on collagen was treated with 1 ml of 0.25 mg/ml of collagenase XI (Sigma-Aldrich) at 37°C with gentle rotation for 20 mins, until the cells detached from the underlying collagen substrate. The samples were collected and centrifuged for 4 minutes at 900 rpm to concentrate the cells into a pellet for RNA purification. Total RNA was extracted using the RNeasy Mini Kit (Qiagen Inc., Valencia, CA) according to the manufacturer's instructions. The concentration of the extracted RNA was quantified using a Nanodrop ND-1000 Spectrophotometer (NanoDrop, Wilmington, DE) (absorption ratio 260/280 > 1.8). The samples were then analyzed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) to ensure adequate sample quality and absence of RNA degradation (RIN value > 8.0).

Microarray assays
Total RNA (1.5 μg/sample) purified from fibroblasts derived from DC cord (n = 6), PF (n = 6) and CT (n = 6) specimens was subjected to analysis with Human WG-6 v3.0 Expression BeadChips (Illumina, SanDiego, CA), each of which contains six arrays on a single BeadChip. Each array is comprised of >48,000 probes derived from human genes in the NCBI RefSeq and UniGene databases. Each array on the BeadChip thus provides genome-wide transcriptional representation of well-characterized genes, gene candidates, and splice variants. A total of 18 separate samples were interrogated on the BeadChip arrays, i.e., there was no pooling of samples within groups. This allowed us to evaluate sample to sample variability in gene expression both within and across each group.
Total RNA was amplified using the Illumina ® Total-Prep™RNA Amplification Kit (Ambion, Austin, TX). Labeling and hybridization was performed according to the Illumina gene expression protocol, featuring a reverse transcription step to synthesize first strand cDNA followed by addition of second strand master mix containing DNA polymerase for second strand cDNA synthesis. This was followed by a single in vitro transcription (IVT) amplification step that incorporated biotin-labeled nucleotides. Subsequent steps included array hybridization, washing, blocking, and streptavidin-Cy3 staining followed by serial non-stringent washing steps to remove unbound conjugate. Following the final rinse, the bioarrays were dried by centrifugation and scanned. The Illumina bioarrays were read in an Illumina BeadArray Reader and the primary intensity data were obtained in standard file format, using Genome Studio software.

Microarray expression analysis
The data obtained from the three cell sources were analyzed for differences in expression as previously described using the multiclass analysis tool of the Significance Analysis of Microarrays program (SAM version 3.02) [23,25]. We then performed unsupervised hierarchical clustering, concordance mapping and created Venn diagrams using the Partek Genomics Suite version 6.4 (St. Louis, MO) to evaluate similarities and differences in the gene expression patterns among the fibroblasts from each of the three cell sources.

Pathway analysis
Statistically significant and differentially expressed genes determined by SAM were subjected to systematic network analysis to determine the primary biological processes and pathways associated with each group using "Ingenuity Pathways Analysis" (IPA, ver.5.0, Ingenuity ® Systems, http://www.ingenuity.com, Mountain View, CA). Only genes that were significantly increased or decreased were included as a defined parameter for the core analysis. Using information stored in the Ingenuity Pathways Knowledge Base (IPKB), genes that mapped to genetic networks were ranked by a score based on the number of eligible molecules contained in the network, which also defines the probability that each network can be achieved by chance alone. A score of 3 is considered threshold of significance which means that there is a 1 in 1000 chance that genes are integrated in the network by pure chance. Canonical pathways were also identified from the IPA library. The significance of the association between the data set and the canonical pathway was measured in two ways within the Ingenuity program: i) as a ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to the canonical pathway was calculated, and ii) by the Fischer's exact test which was used to calculate a p-value determining the probability that the association between dataset genes and the canonical pathway is significant.

Quantitative real time RT-PCR
The protocol for reverse transcription reactions and real time PCR was followed as previously described [23]. Ninety ng of total RNA and 150 ng of random primers (Invitrogen Corporation) were used in 30 μL reactions and real time RT-PCR was performed with triplicate sample reads. The 18 samples were prepared separately such that gene expression analysis was performed using individual samples without pooling. Real time PCR amplification and detection of template was carried out using Applied Biosystems transcript-specific assays including: angiopoietin-like 7 (ANGPTL7) -Hs00221727_m1, Laminin alpha 5(LAMA5) -Hs00966585_m1, Shroom 2 -Hs01113 636_m1, and GAPDH -Hs02758991_g1. Using the comparative critical cycle (Ct) method and using GAPDH as the endogenous control, the expression levels of the target genes were normalized using a 95% confidence interval. The relative expression of the 18 samples (6 for each subset) was averaged and statistical analysis for significance was performed using a Student's t-test. Results shown are from two independent experiments performed in triplicate.

Results
The raw intensity values representing expression of individual transcripts were corrected by subtraction of the background intensity for each array. All values that were at or below background concomitantly in all 3 sample groups were removed yielding 27214 transcripts for further analysis. Box plot analysis of the distribution of these intensity values across the three groups demonstrated comparable dispersion of the individual data sets regardless of patient source, eliminating the need for log normalization or secondary smoothing of the raw intensity data prior to statistical evaluation ( Figure 1).
Transcript intensity values were analyzed for differences in individual gene expression levels using the multiclass comparison of the SAM program (Tusher et al., 2001). The overall raw data distribution was not modified prior to analysis except for median centering of the array data within the SAM program. Results were considered significant at a false detection rate (FDR) of 5% (q ≤ 0.05) and the multiclass comparison revealed 959 transcripts differentially expressed among the 3 classes (delta = 0.47 and FDR = 4.89%). Removal of cDNA clones, open reading frame transcripts, redundant probes and transcripts for hypothetical proteins among the 959 significant transcripts revealed 724 distinct, functionally annotated gene products that were differentially expressed among the 3 groups (Additional file 1).
Post hoc comparisons were performed on these 959 transcripts against the Student's t-distribution using the Partek Software for paired (Dupuytren samples: DC and PF) and the unpaired samples obtained during carpal tunnel release (CT). Three separate univariate comparisons resulted in 894 transcripts detected as significantly different in DC vs CT-derived fibroblasts, 816 significantly different transcripts between PF and CT-derived fibroblasts, and 308 differentially expressed transcripts between DC and PFderived fibroblasts. These lists are depicted in a Venn Diagram ( Figure 2) that illustrates the significant differences specific to each set of individual comparisons (within a circle) versus those shared in more than one of the statistical comparisons (shared between circles). For example, the comparison of DC vs CT-derived fibroblasts yielded 21 transcripts in the non-overlapping area (blue) that were different between these cell types but were not significantly different in the other two comparisons. In contrast, the overlapping region between the PF versus DC circle (green) and the DC vs CT circle (blue) contained 120 transcripts. These 120 transcripts were significantly different in comparisons of DC to either PF or CT-derived fibroblasts but were not different between PF and CT-derived fibroblasts (orange). There were 131 transcripts in the center of the Venn Diagram that overlapped among all three circles, indicating that these transcripts were significantly different in all three comparisons. The most interesting observation, however, is that 622 transcripts (pink) were significantly different in comparisons of CT-to both PF-and DC-derived fibroblasts, but were similarly expressed in DC and PF cells.
Unsupervised hierarchical clustering was performed using the Partek program on log base 2 intensity values of the differentially expressed transcripts identified by the SAM program. Agglomerative clustering was performed to detect dissimilarity based on Euclidean distance with clusters linked based on average values in order to determine if major variations existed in gene expression profiles among the 3 groups [26]. The program was set to delineate the maximum number of clusters based on individual sample comparisons. The primary branch of the resulting dendrogram separated the transcripts of the PF and DC-derived fibroblasts from those of the CTderived fibroblasts, indicating that significant differences in expression were based on patient origin ( Figure 3). However, the DC-and PF-derived fibroblasts also formed clearly distinct clusters separate from each other at the next level of branching. Thus, all three cell groups exhibited distinctly different expression signatures with the greatest difference between the unpaired and the paired (DC-derived) samples.
Concordance analysis was performed to delineate similarities in expression profiles among samples. Concordance was calculated using the initial dataset of 959 significant transcripts (log base 2) with similarity between samples based on Euclidean distance ( Figure 4). Relative similarity is depicted by a color gradient where blue represents the highest similarity diminishing to the most different values represented in red color. The light blue color at the blue-red interface was indicative of sample identity. Based on this analysis, expression of PF- [1][2][3][4][5][6] and DC-derived fibroblasts [7][8][9][10][11][12] exhibited the highest similarity (blue) while both were markedly different from CT-derived fibroblast [13][14][15][16][17][18] profiles, consistent with the results obtained via hierarchical clustering.

Identification of biological pathways
Biological interactions among genes differentially expressed between Dupuytren's-derived and CT-derived fibroblasts were identified using Ingenuity ® System Analysis. Two distinct gene sets of interest were analyzed: 1) genes concordantly expressed between PF-and CT-derived fibroblasts but significantly different in DC-fibroblasts (120 genes); and 2) genes concordantly expressed between DC-and   The program identified 19 and 172 genes respectively that were eligible for function/pathway analysis, comprising 2 and 12 networks respectively with scores of 3 or greater (Tables 1 and 2). The 19 eligible genes identified from the first gene set represented functional categories such as cell death, cell cycle, cellular growth and proliferation-all understandably relevant to a fibroproliferative process like Dupuytren's (Table 1).
In contrast, the 172 genes and 12 networks identified from the second gene set included a much broader array of functional categories, ranging from DNA replication, recombination and repair, to cellular development and cell morphology, to cell-to-cell signaling and interaction. The genes identified were also implicated in a host of other functions including drug metabolism, endocrine system development and function, tumor morphology, molecular transport and even lipid metabolism ( Table 2). What role each of these may play in the progression or recurrence of Dupuytren's pathology is as yet unclear.
An interesting feature of the Ingenuity analysis is that the data suggest the involvement of microRNAs (miRNAs) in the dysregulated processes leading to DC pathology. In examining the 120 genes concordantly expressed between PF-and CT-derived fibroblasts but significantly different in DC-fibroblasts, Ingenuity analysis identified two networks which showed direct interaction of three miRNAs, namely miRNA106A, miRNA17 and miRNA20B (Table 1). Micro-RNA involvement is also suggested by analysis of the genes concordantly expressed between DC-and PF-fibroblasts but significantly different in CT. Three microRNAs are putatively implicated: miRNA202, miRNA122, miRNA214 ( Table 2).

Validation of the differential expression of select gene products by quantitative RT-PCR
Quantitative real time RT-PCR was performed on the same RNA extracted for the microarray analysis on select genes to validate the microarray results (shown in Figure 5). The genes selected for such confirmation have not previously been reported to be differentially expressed in CT-, PF-or DC-derived fibroblasts. Angiopoietin-like7 mRNA (ANGPTL7) was significantly decreased in DC-derived fibroblasts compared to both PFand CT-derived cells, amongst which no significant difference was seen (Figure 6a). LAMA5 mRNA was similarly least in DC-derived cells, but was also significantly less in PF-compared to DC-derived cells (Figure 6b). In contrast, Shroom2 message was dramatically increased in DCderived fibroblasts compared to both CT-and PF-derived fibroblasts ( Figure 6c). As can be seen by comparing the gene expression patterns in Figures 5 and 6, in each case quantitative RT-PCR findings closely matched the results from microarray analysis.

Discussion
We previously reported significant differences in gene expression that distinguish fibroblasts obtained from Dupuytren's Contracture cords versus fibroblasts obtained from the normal fascia of patients undergoing surgical therapy for carpal tunnel syndrome [23]. The culture conditions in this study were modified from our previous report through the use of a type-1 collagen substrate as opposed to an uncoated standard plastic surface. Although the morphological phenotype of the fibroblasts appeared unaffected, many more differentially expressed genes (894 total) were detected in DC versus CT fibroblasts compared to the previous study (40 differentially expressed genes by Illumina BeadArray). Our results can be interpreted to indicate that palmar fascia fibroblasts are sensitive to their culture substrate at a molecular level despite minimal changes in morphological appearance and behavior. These findings can also be interpreted in light of our previous studies indicating that total cellular ß catenin levels in DC cells are altered by interactions with a type-1 collagen substrate [24]. Consistent with this, Ingenuity analysis on   Table 2 Ingenuity analysis of genes concordantly expressed between DC-and PF-fibroblasts but significantly different in CT (Continued) Table 2 Ingenuity analysis of genes concordantly expressed between DC-and PF-fibroblasts but significantly different in CT (Continued)   11 ASB1, BEGAIN, C1ORF9, C8ORF45, CCDC85B,  CTSD, DEFB103A, Eotaxin, EREG, ganglioside  GD1a, ganglioside GD1b, ganglioside GM1,  ganglioside GM2, ganglioside GM3,  ganglioside GT1, GHR, IL32, KIAA0408, LHX4 genes concordantly expressed between DC-and PF-but significantly different in CT-derived cells suggests the involvement of molecules in the Wnt/β-catenin signaling pathway. The gene encoding ß catenin, CTNNB1, was not amongst the dysregulated genes identified in this analysis, consistent with the primary regulation of cellular ß catenin levels through post-translational mechanisms rather than transcriptional induction [27]. Viewed in combination, our findings suggest that the sensitivity of fibroblasts derived from the palmar fascia of DC patients to their immediate environment is evident at both transcriptomic and proteomic levels. These findings may have profound implications for in vitro culture models designed to replicate DC in vivo and highlight the potentially important effects of culture substrates on cellular phenotypes.
There are at least two non-exclusive, potential explanations for the marked differences in gene expression between primary palmar fascia fibroblasts grown on plastic tissue culture trays and those grown on a type-1 collagen substrate. The collagen concentration used for these analyses was 1.9 mg/ml, yielding a relatively soft substrate that can be readily contracted by differentiated myofibroblasts in fibroblast populated collagen lattice assays [24]. A stiff and non-deformable substrate, such as tissue culture plastic, may promote myofibroblast differentiation in a similar manner to that induced by increased tissue density [28]. This differentiation may be due, at least in part, to the activation of latent TGF-ß1 secreted by these cells into their extra-cellular matrix [29]. Under these culture conditions, the differences in gene expression in primary fibroblasts from different tissues may be masked by the changes in gene expression induced during myofibroblast differentiation. In contrast, a soft and readily deformable type-1 collagen substrate may have prevented substrate induced myofibroblast differentiation and allowed for more sensitive detection of the original changes in gene expression that more closely reflect the tissues of origin. Additionally, collagens have been reported to induce a "proliferation permissive" signal through a beta1 integrinmediated PI3 Kinase/Akt pathway in fibroblasts in other systems [30]. Interactions with type-1 collagen may, therefore, activate gene transcription through this, and possibly other, pathways in CT, PF and DC cells and induce the differential expression of genes that are not activated in the absence of collagen. Thus, culturing primary fibroblasts on a type-1 collagen substrate could induce any or all of these effects and explain the marked increase in differentially expressed genes identified in our analyses.
It is notable that each different cell group (DC vs. PF vs. CT) displays a remarkable internal consistency. Six unrelated individuals are compared in each cell group and were found to have highly similar transcriptomic profiles, as evidenced by both hierarchical clustering and concordance mapping. This lends confidence to interpretation of our data as a meaningful representation of the actual differences between the tissue and disease states examined. The most compelling observation is that palmar fascia fibroblasts from phenotypically normal tissue in DC patients closely resemble (but, importantly, are not transcriptomically identical to) cells derived from clearly phenotypically abnormal tissue. As a corollary, both fibroblasts from diseased DC cords and uninvolved DC fascia are markedly transcriptomically dissimilar to fibroblasts from carpal tunnel fascia. This observation is strongly suggestive that a major contributor to the tendency for DC recurrence may be a genetic predisposition towards disease of phenotypically normal palmar fascia in DC patients. Such an interpretation would be consistent with previous reports that DC has a genetic basis {eg. Hu et al., 2005 [3]}. Our findings may also reflect on recent reports suggesting that DC cells may be derived from mesenchymal stem cells (MSCs) in perinodular fat and skin [31]. While our analyses do not rule out the possibility that MSCs from these tissues contribute to (or perhaps initiate) disease progression or recurrence, the marked similarity between the transcriptomes of DC and PF cells and their relative dissimilarity to CT cells strongly implies that the vast majority of cells in DC cord tissue are likely to be derived from PF cells.
These observations make it possible to consider the progression and recurrence of DC in a model based on principles first elucidated in tumor biology, the so-called "two-hit hypothesis" [32] first described by Nordling in 1953 [33]. In such a model we hypothesize that an inherited allelic defect (or possibly a group of allelic defects) constituting a first "hit" alters palmar fascial physiology in DC patients, but not sufficiently so as to lead to the DC nodule/cord phenotype. Progression to frank and active disease would require a second somatic mutation (the second "hit") in the remaining allele in the affected tissue type. Such a model would account for the observation that DC is a heritable disorder, usually presenting later in life (with correspondingly greater chances of having accrued the second requisite mutation), but with widely varying penetrance. Given that both PF and DC cells display hundreds of genes differentially expressed at the mRNA level compared to control CT cells, the putative etiological locus may be involved in transcriptional regulation. Of course, it remains possible that DC is a heterogeneous disorder with a similar phenotype arising from multiple disparate underlying genetic loci, nor does such a model eliminate the possibility of de novo sporadic cases of DC.
These findings also have potential implications for in vitro studies designed to dissect the molecular mechanisms of DC. If patients with a genetic predisposition to develop DC have inherited PF cells that exhibit profound differences in gene expression to normal palmar fascia (CT) cells, then PF cells represent the most clinically relevant controls for testing treatments designed to prevent DC progression. In contrast, normal palmar fascia cells, such as CT cells, are valuable for analyses designed to identify the molecular characteristics that distinguish DC cells from normal cells. Thus, while identifying differences in signaling pathway activation between DC and CT cells [34] are useful analyses to characterize differences from normality, it is clearly essential that such pathways are also assessed in PF cells before they are considered potential therapeutic targets. A study by Pavelic et al., [35] employing a proteomic approach has identified several proteins interacting in various signaling pathways in DC tissues when compared to unaffected palmar fascia from DC patients. This identification of protein-protein interaction potentially opens new therapeutic targets to reduce the recurrence of DC. Optimally, comparisons of DC cells to both CT and PF cells are likely to provide different, but equally valuable, insights into DC progression and recurrence that will allow researchers to distinguish and identify molecular or other therapies that target DC, but not PF or CT, cells.
The remarkable internal consistency between samples in the microarray was also reflected in the findings obtained from real time RT-PCR assays used to validate the microarray results. In seeking potentially novel targets to limit the progression of DC, we chose to directly examine three select genes that have not been previously implicated in the progression or recurrence of DC. Human ANGPTL7 had been characterized as a potent target gene of the WNT/β-catenin signaling pathway and currently is a pharmacogenomics target in the fields of oncology and regenerative medicine [36]. Studies by Kuchtey et al., (2008) [37] and subsequent studies by Comes et al., (2010) [38] showed that ANGPTL7 has the potential to alter extracellular matrix formation, a finding that may also be relevant to DC. We find that in DC-derived fibroblasts ANGPTL7 mRNA is substantially decreased in comparison to CT-and PF-derived fibroblasts (Figure 5a and 6a). The functional importance of