Dupuytren's: a systems biology disease
© BioMed Central Ltd 2011
Published: 12 September 2011
Skip to main content
© BioMed Central Ltd 2011
Published: 12 September 2011
Dupuytren's disease (DD) is an ill-defined fibroproliferative disorder of the palm of the hands leading to digital contracture. DD commonly occurs in individuals of northern European extraction. Cellular components and processes associated with DD pathogenesis include altered gene and protein expression of cytokines, growth factors, adhesion molecules, and extracellular matrix components. Histology has shown increased but varying levels of particular types of collagen, myofibroblasts and myoglobin proteins in DD tissue. Free radicals and localised ischaemia have been suggested to trigger the proliferation of DD tissue. Although the existing available biological information on DD may contain potentially valuable (though largely uninterpreted) information, the precise aetiology of DD remains unknown. Systems biology combines mechanistic modelling with quantitative experimentation in studies of networks and better understanding of the interaction of multiple components in disease processes. Adopting systems biology may be the ideal approach for future research in order to improve understanding of complex diseases of multifactorial origin. In this review, we propose that DD is a disease of several networks rather than of a single gene, and show that this accounts for the experimental observations obtained to date from a variety of sources. We outline how DD may be investigated more effectively by employing a systems biology approach that considers the disease network as a whole rather than focusing on any specific single molecule.
Systems biology combines mechanistic modelling with quantitative experimentation in studies of networks [27–34] and aims at understanding how the interaction of multiple components within a cell, tissue, organ or indeed individual leads to much of biological function and obfuscates correlations with single genes. Systems-level approaches have begun to help comprehension of network control, (dys-)regulation, and function [35–38]. This has improved the understanding of certain disorders , and has provided new rationales for drug discovery [40–42]. The complex biology of DD may constitute an invitation to a systems level approach. In this review, we outline such an approach.
Clinical examples of fibrosis include renal interstitial fibrosis , scleroderma , sarcoidosis , idiopathic pulmonary fibrosis , retroperitoneal fibrosis  and DD . DD tissue shows increased deposition of collagen III relative to collagen I and increased levels of collagen hydroxylation and glycosylation . DD is thought to arise either from a defect in the wound repair process or from an abnormal response to wounding. The presence of immune cells and related phenomena in DD tissue suggests DD may be immune-related [50–53]. Cellularity (quantified as the cellular density) of the DD nodules (see below) is indicative of the activity of the disease . DD has been classified into three stages co-existing in the same specimen, that is, proliferative, involutional and residual, further subdivided into the essentially fibrous nodules, reactive tissue, and residual tissue. It contains two structurally distinct fibrotic elements: the nodule is a highly vascularised tissue containing many fibroblasts, with a high percentage being recognised as myofibroblasts due to their expression of the α-smooth muscle actin; and the cord is relatively avascular, acellular, and collagen-rich with few myofibroblasts. The nodule may develop into the cord as the disease progresses over time or the two structures represent independent stages of the disease. Macroscopically, neither the deep retinacular tissue that includes the transverse palmar ligament or fascia, also known as 'Skoog's fibres', nor the fibrous flexor tendon sheaths appear to be involved in DD. Other areas are affected macroscopically but at irregular depth and distribution, with the more superficial layers and ulnar side of the palm being affected most.
The specialised mesenchymal cells expressing smooth muscle α-actin may explain the contractility observed in DD [11, 54–56]; they resemble the myofibroblasts of granulation tissue thought to be responsible for contraction during wound healing. The Dupuytren myofibroblast synthesizes fibronectin, an extracellular glycoprotein that connects myofibroblasts and connects them to the extracellular stromal matrix through an integrin.
According to genome-wide gene expression profiles, fibroblasts come in various subtypes , perhaps due to 'topographic differentiation'; that is, distinct phenotypes persisted in vitro even when fibroblasts were isolated from the influence of other cell types . Chang et al.  did not evaluate diversity or cell heterogeneity in DD. All of those evaluated had the morphology of elongated, spindle-shaped cells. Fibroblast cultures were uniformly positive for a mesenchymal immunofluorescence marker, but negative for markers of epithelial, smooth muscle, endothelial, perineural, and histiocytic cells. Different passages of the same fibroblast culture clustered with each other, indicating that their in vitro phenotypes were stable. Several components implicated as modulators of transdifferentiation of DD fibroblasts into myofibroblasts have been reported [59–69]. Among the cytokines, transforming growth factor-β is thought to be a significant inducer of myofibroblast transdifferentiation because of its ability to up-regulate α-smooth muscle actin and collagen in fibroblasts, both in vivo and in vitro .
A study performed in a five generation Swedish family suggested that DD was inherited in an autosomal dominant pattern . Linkage analysis implicated a single region of approximately 6 cM between markers D16S419 and D16S3032 at a logarithm of the odds (LOD) score >1.5. Genotyping of four siblings affected by the disease but from another branch of the family together with the use of additional microsatellite markers produced a maximal LOD score of 3.2 (for D16S415), with four other markers producing LOD scores >1.5. When a disease is dominant, it is likely to be caused by a single allele of a single gene, and by the molecule it encodes. From this perspective, the above findings would suggest that DD is a single gene disease. To date, however, linkage to a single gene has not been reported at a LOD that is much more significant than the marginal value of 3 in this Swedish study and the penetrance in this study was incomplete. In addition, the disease develops at an advanced age, there are many more sporadic cases of DD, and there are few such families for which the genetic analysis has been performed. Indeed, other studies have shown association of the disease with other loci, including a positive association with HLA-DRB1*15 on chromosome 6 in Caucasians . A study of 20 British DD patients with a maternally transmitted inheritance pattern demonstrated a mutation within the mitochondrial genome (mitochondrial 16S ribosomal RNA region) in 90% of patients . The defective mitochondria generated abnormally high levels of free radicals and induced defects in apoptotic mechanisms.
A relation between localised ischaemia, reactive oxygen species (ROS) and DD was projected from a study in which palmar fascia from 10 DD individuals were subjected to 0 to 60 minutes of tourniquet ischaemia DD . The concentration of hypoxanthine was six-fold higher in Dupuytren's palmar fascia compared to the palmar fascia from ten suitable control patients (having carpal tunnel decompression), implying that measuring metabolites directly in tissue could help understand DD. Xanthine oxidase catalyzes the removal of hypoxanthine, generating super oxide free radicals and hydrogen peroxide as by-products that would damage the perivascular connective tissue, which fibroblasts would attempt to repair. Upon addition of free radicals to fibroblast cultures from DD palmar fascia, lower concentrations of ROS stimulated fibroblast proliferation, explaining the observed increase in collagen type III . The hypoxanthine was more abundant in nodular areas than in the tight fibrous cords. We propose that microvessel narrowing, leading to localised hypoxic conditions, may be one cause of DD, secondary to age, smoking and other environmental factors as discussed by Shih and Bayat .
The DD proteomics venture began in 2006 with the study of protein expression profiles in an attempt to identify potential disease biomarkers [77, 78]. In one study, two-dimensional gel electrophoresis was performed to extract proteins from diseased tissue (nodule and cord), the Skoog's fibres, and normal control tissues. MALDI-TOF-MS (matrix-assisted laser desorption ionization time-of-flight mass spectrometry) generated a peptide mass fingerprint that was used to search protein databases; however, the authors did not report names of identified proteomic changes in their abstract. SELDI-TOF-MS (surface enhanced laser desorption/ionization time-of-flight mass spectrometry) using Ciphergen's SELDI-TOF-MS Protein Biological System II (PBSII) ProteinChip reader  revealed several differentially expressed low molecular weight (<20 kDa) tissue proteins and identified three disease-associated protein features (4,600.8 Da, 10,254.5 Da, and 11,405.1 Da) that were elevated (5-, 12-, and 4-fold, respectively). Three potential low molecular weight protein markers (p4.6DC, p1ODC, p11.7DC) for DD were identified. An integrative proteomic-interactomic approach  coupling two-dimensional gel electrophoresis with mass spectrometry compared the proteomic profile of DD tissue with that of unaffected patient-matched palmar fasciae tissue and found several proteins correlated with DD. The findings were used to create a protein-protein interaction network (interactome) map on the basis of the proposed interactions in the Human Interactome Map (HiMAP)  and the Search Tool for the Retrieval of Interacting Proteins (STRING) ; however, several proteins were added to fill gaps in order to yield a complete network, which then involved extra- and intracellular signalling, oxidative stress, cytoskeletal changes, and alterations in cellular metabolism. In particular, ERBB-2 and insulin-like growth factor 1 receptor (IGF-1R) and the Akt signalling pathway emerged as novel components of pro-survival signalling in Dupuytren's fibroblasts. One should exercise care, however, not to over-interpret these results, as they are partly based on inference from other protein interaction data obtained in different contexts. In addition, increased activity of pathways need not involve increased protein levels  and increased pathway expression may be homeostatic rather than aetiological.
As more and more aspects associating with DD are revealed, we see less and less forest (understanding of the disease) for the trees (its many molecules) . Even if the disease were set in motion by a single genetic factor, its aetiology would involve many diverse processes such that DD will be co-determined by the many factors that regulate those processes. If indeed the networks governing differentiation of normal fibrocytes of the palm of the hand are perturbed irreversibly so that they differentiate into muscle-like tissue without the proper controllers of contraction and relaxation, then different sets of genetic perturbations could lead to DD. In this context DD may be much like cancer .
The dilemma is that although we now have an unprecedented set of methodologies for the identification and analysis of all the molecules in living cells, that methodology alone is not enough. We need something substantially more to understand how all those molecules interact to create functional networks. Seeing more molecules may not help our understanding; seeing the connections between them and more mechanism might.
Most diseases have multiple genes associated with them. Such diseases might be considered to be a group of single-gene diseases; that is, many different diseases each being caused by a different single gene lesion, but all with similar phenotypes . This would explain the association of multiple genes with the disease. In the case of a group of single gene diseases, no other faulty molecules should be important for that individual disease and no other gene changes (for example, polymorphisms) or conditions (for example, diet) should influence the disease severity. Notably, a single patient's transcriptome should then show only changes in the single molecular culprit and not in other factors controlling the network leading to the disease; and in the transcriptomes of different patients suffering from the same disease group, that single molecular culprit should be different. For DD this is not what is observed (see above). As illustrated in Figure 3b, in a systems biology disease the function that is compromised depends cooperatively on a number of pathways, the functioning of each of which again depends on many cooperating molecular factors. In systems biology diseases one would typically find multiple changes in the transcriptome or proteome of each patient, differing between individual patients but such that all have a very similar disabling effect on network function. Identifying a disease as a systems biology or network disease does not dispel molecules from its pathology: molecules are always involved. The issue is whether the change in networking of the molecules is crucial for the disease, that is, whether the disease is more a consequence of faulty networking than of an individual malfunctioning molecule.
DD has been identified as a disease inherited in an autosomal dominant pattern  (and see above). It was linked to a single 6cM region on chromosome 16. This would suggest that all DD patients should have a mutation in this part of their genome, and that transcriptomes of DD patients should be altered in terms of the level of the transcripts encoded by this part of the genome or in terms of the coding sequence of one of those transcripts. However, the dominance was incomplete, weak, and has only been observed in a single Swedish family. This suggests that the genes on chromosome 16 are only dominant when other genes in the genome are of certain allelic forms. Moreover, in many other cases the expression levels of many other mRNAs were changed, although it remains to be analyzed whether in those studies there was always a change in mRNAs from the 6 cM region on chromosome 16. In our own studies, DD nodule transcriptomes of individual patients have all exhibited multiple changes in mRNA levels, and although these changes overlapped, they were not identical between individuals. The proteome did not point to a single causative protein either. The functional studies pointed to myofibroblast enrichment, although not clearly as the sole cause, and neither was a causal relationship between a gene on chromosome 16 in the 6 cM region and differentiation of myofibroblasts established. This all shows that DD is not a single-gene disease and suggests that it is not just a group of pure single-gene diseases either. It is much more likely to be a systems biology disease.
Surgical intervention is still the current mainstay of treatment for DD, usually involving fasciotomy, fasciectomy or dermofasciectomy . A variety of non-operative techniques have been practiced but have failed to give long-lasting benefits. More recently, clinical investigation (phase II and III), including the use of clostridial collagenase injections, have shown encouraging results in some DD patients [89, 90], but long-term follow-up results are required before this can be advocated as standard procedure in place of surgery.
Since association of DD with a single gene is inconclusive, the present mainstream research paradigm may be unlikely to lead to a full understanding of it. The experimental data appear to be more consistent with DD being a systems biology disease. Therefore, a different approach to the disease should be considered; analysis, diagnosis and therapy should target pathways rather than genes or their products. The concept of a 'candidate gene' should be replaced with that of 'candidate pathway(s)'. Studies should be aimed at elucidating cause-effect chains, rather than disease correlations. From the experimental data, alterations in pathways should be inferred. Using transgenic and antisense approaches in cell lines, these pathway alterations should then be induced and the predicted development into a DD cellular phenotype tested. The pathways are expected to be integrals of gene expression, signalling and metabolic networks, as should be the approach and data analysis.
A hypothesis-driven systems biology would be based on a priori observations in human, in vitro or in vivo (linkage and expression studies, for example), or on knowledge of related diseases (such as plantar fibromatosis, peyronies, musculo-aponeurotic fibromatosis and even keloid disease). Inter-relationships would be sought between hypothesized underlying mechanisms governing these fibrotic disorders and physiological changes predicted based on molecular and environmental changes impacting on those mechanisms. This could then be extended to understand inter-versus intra-individual variability. Altering the networks using multiple molecular interventions in a tissue culture model system for DD would enable the hypotheses to be tested.
Such an approach should also help put into perspective existing inconclusive discoveries and maximize the utilization of data obtained from molecular approaches from molecular biology, perhaps in the sense of a live model repository . Indeed, omics experiments generate heterogeneous data and metadata torrents. Systems biology modelling would serve as a way to organize the data rationally and efficiently. Emerging standard languages, such as the Systems Biology Markup Language (SBML) , Cell Markup Language (CellML) , Cell System Markup Language (CSML) , Biological Pathway Exchange (BioPAX)  and Systems Biology Graphical Notation (SBGN) , and modelling tools such as COmplex PAthway SImulator (COPASI) , Cytoscape  and Pathway databases (for example, Ingenuity pathway analysis software) facilitate data representation and inter-operability from leading multidisciplinary research groups . Additionally, the Java Web Simulation (JWS) facility, which quality controls kinetic models and puts them into a live repository, enables through-web experimentation in silico for scientists naive with respect to modelling, although it is experiment-dependent . BioModels is a parallel model repository not built for in silico experimentation but focusing on annotation .
If we are right, and the systems biology research paradigm is adopted, the rewards should be substantial: no longer will the data collected in this field disappear into the diaspora of the experimental literature; they will be analyzed in terms of network models and, when informative, connected with proper data annotation . Once the network hypotheses are proven, they will underpin the development of new rational biomarker strategies, and become the starting point for potential therapeutic interventions and prophylaxis.
More concretely, even though myofibroblasts obtained from different stages of DD may exhibit features that could trigger contraction and uncontrollable growth, neither the diversity of these cells nor the extent or nature of their local specificity in situ with regard to differentiation have been examined systematically. Remodelling of vascular connective tissue should be of fundamental importance as DD progresses over time and the ability of that tissue to be remodelled could be an important factor in the development of the disease. Matrix remodelling and matrix turnover are controlled by a complex network of cell-cell and cell-matrix interactions . Mathematical modelling of these networks combined with targeted experimentation should help deduce the net outcome of the balance between proliferative and degradative processes.
Together, genomics, transcriptomics, proteomics, fluxomics, lipidomics, interactomics, glycomics and secretomics studies of biofluids within the DD system have the potential to improve our understanding of the disease immeasurably [104–106]; however, the 'together' concept requires an integration that is only plausible through systems approaches.
Recognition of DD as a systems biology disease should also affect therapy. Treatment need not target the reversal of a molecular event but rather return a network affected by a combination of molecular culprits (for example, resulting from single-nucleotide polymorphisms) to its normal state. This return may be attempted by interventions in the network that may otherwise be unrelated to the molecular culprits. In this sense the recognition that DD is a network disease should increase the number of treatment options. More specifically, one should look at a number of intracellular networks (metabolic, signal transducing, or otherwise) affecting the disease and then consider treatment with combinations of existing drugs (and behavioural advice) that target all those networks. The relative dosages should be optimized, allowing for the optimum to differ between patients.
logarithm of the odds
reactive oxygen species.
SR thanks BBSRC and EPSRC for her studentship via the Doctoral Training Centre, at the Manchester Centre for Integrative Systems Biology, UK. HW [http://www.systembiology.net/support/] and RG thank various funders but in particular the EPSRC and BBSRC for their support of the above two centres.