Journal of Diabetes Research

Journal of Diabetes Research / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 214353 |

Michael E. Munson, James S. Wrobel, Crystal M. Holmes, David A. Hanauer, "Data Mining for Identifying Novel Associations and Temporal Relationships with Charcot Foot", Journal of Diabetes Research, vol. 2014, Article ID 214353, 13 pages, 2014.

Data Mining for Identifying Novel Associations and Temporal Relationships with Charcot Foot

Academic Editor: Mark A. Yorek
Received31 Jan 2014
Revised28 Mar 2014
Accepted28 Mar 2014
Published27 Apr 2014


Introduction. Charcot foot is a rare and devastating complication of diabetes. While some risk factors are known, debate continues regarding etiology. Elucidating other associated disorders and their temporal occurrence could lead to a better understanding of its pathogenesis. We applied a large data mining approach to Charcot foot for elucidating novel associations. Methods. We conducted an association analysis using ICD-9 diagnosis codes for every patient in our health system ( million with 41.2 million time-stamped ICD-9 codes). For the current analysis, we focused on the 388 patients with Charcot foot (ICD-9 713.5). Results. We found 710 associations, 676 (95.2%) of which had a P value for the association less than and 603 (84.9%) of which had an odds ratio > 5.0. There were 111 (15.6%) associations with a significant temporal relationship . The three novel associations with the strongest temporal component were cardiac dysrhythmia, pulmonary eosinophilia, and volume depletion disorder. Conclusion. We identified novel associations with Charcot foot in the context of pathogenesis models that include neurotrophic, neurovascular, and microtraumatic factors mediated through inflammatory cytokines. Future work should focus on confirmatory analyses. These novel areas of investigation could lead to prevention or earlier diagnosis.

1. Introduction

Charcot foot is a rare disease that often results in significant complications including the progressive destruction of bones and joints, loss of mobility, infections, amputations, and even death. There are approximately 40,000 new cases of Charcot foot each year. While no population-based studies of patients with Charcot foot have been described, single center prevalence estimates range from 0.1% to 0.9% [13].

Charcot foot most often occurs in the setting of diabetes, and an increased incidence of Charcot foot has been associated with the following: obesity, neuropathy, increased age, diabetes > 6 years, elevated hemoglobin A1c, renal failure, iron deficiency anemia, osteoporosis, and rheumatoid arthritis [4, 5]. Further, Matricali and colleagues noted a high rate of Charcot foot in simultaneous pancreas-kidney transplant patients [6], a treatment that is often used in patients with diabetes and end-stage renal failure.

Despite the well-known association of Charcot foot with diabetes and related disorders, the pathogenesis is still unknown. There are multiple theories describing the pathogenesis of Charcot foot, although the inciting mechanism is generally thought to be caused by trauma in a susceptible person with sensorimotor and/or autonomic neuropathy leading to inflammation, bone destruction, subluxation, and disordered repair in later stages.

In order to develop a better understanding of Charcot foot, we sought to use a novel approach to more extensively explore all possible diagnoses that may be associated with Charcot foot.

The increasing adoption of electronic health records (EHRs) presents the opportunity to mine vast amounts of data generated by many physicians over time in order to uncover novel associations. New computational approaches, often collectively referred to as “data mining” [79], have been developed that can explore large populations of patients to detect patterns or relationships that might otherwise go undetected. This may be especially true when diseases are rare and the patients are seen by a wide variety of providers and specialists, all of whom may treat aspects of the disease but may not observe some of the larger-scale patterns that might exist. Patients with Charcot foot, for example, may be seen by multiple providers including primary care physicians, podiatrists, endocrinologists, neurologists, orthopedic surgeons, and others.

Here we describe an exploratory analysis using a data mining approach applied to patients from a large, integrated health system [10, 11]. The focus of this study is not on the data mining methodology itself, but rather on the results obtained using data mining for patients who had a diagnosis of Charcot foot. We present findings that include well-known and well-established associations, which helps to confirm the validity of the approach, as well as other poorly known or even unknown associations. The technique is also able to determine the temporal relationships between Charcot foot and its associated disorders, and these relationships could have implications for explaining disease mechanisms and progression. Additionally, some of the potentially novel associations could be used to help form hypotheses about the etiology and pathogenesis of Charcot foot that could lead to future investigations.

2. Methods

2.1. Clinical Setting

This study was conducted at the University of Michigan Health System (UMHS), a large, integrated, academic, tertiary care facility with three hospitals and nearly 1,000 inpatient beds. At UMHS there are nearly 50,000 surgeries and over 45,000 inpatient visits annually. Additionally, there are almost 2 million outpatient visits per year across more than 120 outpatient clinical locations scattered primarily throughout Southeastern Michigan. There are six specialty centers including a Comprehensive Diabetes Center which houses a high volume podiatry clinic with three board certified podiatrists who specialize in the management of diabetic foot complications and amputation prevention. This study was considered exempt by our medical school’s institutional review board because the dataset contained no identifiers.

2.2. Analytic Approach

We conducted an association analysis using international classification of disease version 9 (ICD-9) diagnosis codes among every patient in our health system. These codes had been assigned for administrative (i.e., billing/reimbursement) purposes over an approximately twenty-year period. A detailed description of the methodological approach, as well as its limitations, has been described elsewhere [10, 11]. Briefly, for each pair of ICD-9 codes (i.e., diagnoses) we calculated the chi-square statistic and its associated value to determine if the two codes were significantly associated with one another. Further, we used the binomial test to determine if there was a significant temporal relationship between each pair of codes. That is, we tested if one code preceded the other code in a nonrandom manner. Such a temporal relationship can suggest, but in no way establishes, potential causality. To determine the potential role diabetes had in the associations, we also compared the frequency of diagnoses associated with Charcot foot in patients with diabetes (defined as codes 250.x, where x is one of many subcodes for diabetes) as well as those without diabetes. This was done using the proportion test. All statistical analyses were conducted using R version 2.15.3. Data were visualized with network diagrams using Cytoscape version 2.6.3 [12].

In this analysis, values and statistical significance should be interpreted with several caveats in mind. The first is that no adjustments have been made for multiple hypothesis testing, which is often used in genomic analyses. With such an exploratory analysis, values can be interpreted as a “ranking” in which results can be compared, and with lower values suggesting more significant relationships, but we do not define any specific threshold for which a result is truly statistically significant. values at or near the “traditional” threshold of 0.05 should be interpreted with caution.

2.3. Expert Review

The complete dataset included approximately 1.6 million patients and 41.2 million time-stamped ICD-9 codes. For the current analysis, we focused only on the patients who had been assigned the ICD-9 code 713.5 (“arthropathy associated with a neurological disorder”) that has been used in other peer-reviewed publications to describe Charcot foot [5, 1316]. ICD-9 codes are most often used for administrative/billing purposes rather than directly for clinical care and thus can sometimes be inaccurately assigned [17]. It is not known how often the ICD-9 code 713.5 is inappropriately applied, but one study that included a small number of patients with the code 713.5 found that 92% of the cases did indeed have Charcot foot as verified by chart abstraction [16].

These associations were then independently reviewed by three Charcot foot domain experts. The associations were divided into three categories by the reviewers: (1) associations that are generally well-known, (2) associations that were potentially novel or poorly known, and (3) associations that provided minimal insight about the relationship because the codes either were nonspecific or likely represented a Charcot foot diagnosis that was initially misdiagnosed as a swollen foot. [18]. Examples of nonspecific codes were 079.99 (“viral infection, not otherwise specified”) and 693.8 (“nonspecific abnormal radiological finding of the breast”). Examples of codes that were likely misdiagnoses that were later found to be Charcot foot included 274.9 (“gout, not otherwise specified”) and 715.90 (“osteoarthrosis, unspecified site”). The decision to exclude probable misdiagnoses was supported by our observation that the misdiagnoses tended to precede the diagnosis of Charcot foot in most patients. For example, in the case of the two diagnoses described above, gout (274.9) preceded the diagnoses of Charcot foot 75% of the time and osteoarthritis (715.90) preceded Charcot foot 70% of the time. These temporal findings can be found in the Supplementary Table in the Supplementary Material (available online at

3. Results

Among the overall population of 1.6 million patients, approximately 61,400 had an ICD-9 code for diabetes mellitus (ICD-9 250.x). Additionally, there were 388 patients with an ICD-9 code for Charcot foot (ICD-9 713.5), resulting in a low incidence (0.6%) even among our diabetic population. Of these 388 Charcot foot patients, 282 (72.7%) also had one of the 250.x diabetes codes assigned to them. Among the 106 patients with Charcot foot who did not have a diabetes code, there was no prevailing comorbid condition shared by the patients. The top five diagnoses appearing across the 106 patients lacking a diabetes code were ICD-9 codes 729.5 (; pain in limb), 719.47 (; pain in joint, ankle/foot), 733.90 (; disorder, bone/cartilage), 724.2 (; lumbago), and 719.46 (; pain in joint, lower leg).

At a large, tertiary academic medical center such as ours, it is possible that patients with rare diseases are often referred for specialized care but receive a majority of their care elsewhere. As such, it is possible that patients with Charcot foot were referred to UMHS with that diagnosis being their initial diagnosis in many cases. However, we did not find that to be the case. Twenty (5.2%) of the patients in our cohort were given a diagnosis of Charcot foot on their first visit to our health system, while an additional 10 patients (30 overall; 7.7% of cohort) were given the diagnosis within the first 30 days of visiting our health system. However, there was a wide range between the time of first contact with our health system and a diagnoses of Charcot foot: the mean time interval was 6.6 years (95% CI 6.1–7.1 years), and the median time interval was 6.2 years. It is also possible that our patients were only given a diagnosis of Charcot foot and then sent back to their referring physician outside of our network. This could cause an underrepresentation of diagnoses for those patients in our database. However, only one patient among the 388 had a single diagnosis of Charcot foot and no other diagnoses. The mean number of distinct diagnoses based on ICD-9 codes for the patients in our cohort was 84.8 (95% CI 77.9–91.8), with a median of 68.5 distinct diagnoses per patient.

Our data mining approach revealed 710 pairs of associations between Charcot foot and other disorders, 676 (95.2%) of which had a value for the association less than 1.0 10−5 and 603 (84.9%) of which had an odds ratio > 5.0. There were 111 (15.6%) pairs that had a significant temporal relationship with a value of less than 1.0 10−3. A full list of all 710 associations with Charcot foot can be found in the Supplementary Table. Figure 1 provides a high-level overview of some of the diagnoses associated with Charcot foot or associated with other closely related diagnoses. This figure also shows the strongest temporal relationships, using arrows between two diagnoses to represent the temporality.

A subset of known associations that were substantiated by our method is summarized in Table 1. Among these known associations, type 2 diabetes mellitus (T2DM; ICD-9 250.00) had the strongest temporal component with 87.2% of the diagnoses preceding Charcot foot when the two diagnoses cooccurred in the same patient (temporal = 8.5 10−37). Many of the diagnoses presented in Table 1 are interrelated, primarily because diabetes is known to cooccur with numerous other problems including obesity, chronic kidney disease, and anemia. Further, some of the results presented in Table 1 may be misleading because they show only the frequency of single codes in the population. Thus, coding variability is not well captured in the tables. For example, code 356.9 (peripheral neuropathy) appeared in 41.5% of the patients with Charcot foot. However, when including additional codes that could represent a peripheral neuropathy (i.e., 250.60, 250.61, 250.62, 250.63, 355.8, 355.9, 356.0, 356.1, 356.9, 337.1, and 354.2), we found that these codes collectively covered 298 patients, or 76.8% of the Charcot foot population.

Disease categoryICD-9 codeICD-9 descriptionNumber (%) with both disordersAssociation odds ratioAssociation valueTemporal value

Endocrine, nutritional and metabolic diseases, and immunity disorders250.00*Type 2 diabetes mellitus243 (62.6)55.4<5 ×  
278.00Obesity52 (13.4)11.7

Nervous system356.9Peripheral neuropathy161 (41.5)101.2<5 ×  

Blood and blood-forming organs280.0Anemia due to chronic blood loss16 (4.1)9.9 0.80
285.9Anemia, not otherwise specified98 (25.3)14.9 0.085
285.29Anemia of other chronic diseases10 (2.6)43.1 1
280.9Iron deficiency anemia32 (8.2)14.5 0.60

Skin and subcutaneous tissue682.6Cellulitis/abscess of leg87 (22.4)41.9<5 ×  
682.7Cellulitis/abscess of foot98 (25.3)133.8<5 ×  

Genitourinary system593.9Renal & ureteral disorder, not otherwise specified119 (30.7)20.9<5 ×  
585Chronic kidney disease64 (16.5)29.4

Musculoskeletal system and connective tissue714.0Rheumatoid arthritis37 (9.5)13.5
733.00Osteoporosis 44 (11.3)12.1 0.29
733.90Bone/cartilage disorder, not otherwise specified233 (60.1)26.0<5 ×  

Injury and poisoning996.81Complication of kidney transplant44 (11.3)43.3<5 ×   3.9 × 10−4

External causes of injury and supplemental classificationV42.0Kidney transplant50 (12.9)42.1<5 ×  

This table only lists ICD-9 code 250.00, not the other diabetes codes that also appeared in this cohort of patients.

New or less well-known associations that were determined by expert review to be interesting in the context of current Charcot foot etiological theories, and that had a strong temporal component, are summarized in Table 2. The top three most significant associations with Charcot foot that have not been well described are related to circulatory system disorders (ICD-9 codes 414.00, 440.21, and 443.9). The top three novel/poorly known associations with the strongest temporal component were cardiac dysrhythmia (ICD-9 429.9), pulmonary eosinophilia (ICD-9 518.3), and volume depletion disorder (ICD-9 276.5). Two disorders related to lumbar disc disease were also associated with Charcot foot: (1) displacement of lumbar intervertebral disc without myelopathy (ICD-9 722.10; association value 3.48 10−14; odds ratio 5.0) and (2) degeneration of lumbar or lumbosacral intervertebral disc (ICD-9 722.52; association value 3.32 10−54; odds ratio 10.2). These two lumbar disc diagnoses generally preceded Charcot foot, but their levels of temporal significance were lower than our threshold for inclusion in Table 2 (temporal values 0.04 and 0.01, resp.).

Disease categoryICD-9 codeICD-9 descriptionNumber (%) with both disordersAssociation odds ratioAssociation valueTemporal value

Endocrine, nutritional and metabolic diseases, and immunity disorders244.9Hypothyroidism42 (10.8) 6.1
272.0Hypercholesterolemia82 (21.1)17.9
272.4Hyperlipidemia173 (44.6)17.4
276.3Alkalosis13 (3.4)14.0
276.5Volume depletion disorder84 (21.6)8.1

Sense organs366.16Senile nuclear cataract110 (28.4)13.3

Circulatory system401.1Benign essential hypertension113 (29.1)13.3
402.10Benign hypertensive disease without heart failure15 (3.9)6.6
411.1Intermediate coronary syndrome38 (9.8)13.9
413.9Angina pectoris41 (10.6)13.1
414.0Coronary atherosclerosis18 (4.6)7.2
414.00Coronary atherosclerosis, vessel type155 (39.9)18.6<5 ×  
414.01Coronary atherosclerosis, native artery97 (25.0)16.3
414.02Coronary atherosclerosis, autologous graft11 (2.8)15.7
414.9Chronic ischemic heart disease49 (12.6)14.7
424.0Mitral valve disorder32 (8.2)7.4
425.4Cardiomyopathy, primary 44 (11.3)11.6
427.9Cardiac dysrhythmia168 (43.3)7.6
428.0Congestive heart failure84 (21.6)11.0
428.9Heart failure, not otherwise specified38 (9.8)10.1
429.9Heart disease, not otherwise specified39 (10.1)10.5
440.21Atherosclerosis of native arteries of the extremities with intermittent claudication93 (24.0)48.5<5 ×  
443.9Peripheral vascular disease95 (24.5)38.5<5 ×  

Respiratory system465.9Acute upper respiratory infection93 (24.0)2.7
466.0Acute bronchitis75 (19.3)5.5
473.9Chronic sinusitis35 (9.0)4.4
477.9Allergic rhinitis40 (10.3)2.3
493.90Asthma (without status asthmaticus)46 (11.9)3.0
518.3Pulmonary eosinophilia125 (32.2)7.1

Digestive system530.81Esophageal reflux64 (16.5) 4.7
560.9Intestinal obstruction, not otherwise specified30 (7.7)6.5

Genitourinary system596.9Bladder disorder, not otherwise specified12 (3.1)7.3
607.84Impotence, organic origin33 (8.5)14.0

Skin and subcutaneous tissue692.9Contact dermatitis52 (13.4) 3.3

Musculoskeletal system and connective tissue724.2Lumbago117 (30.2) 6.7
729.1Myalgia/myositis, not otherwise specified55 (14.2)9.1
715.90Osteoarthrosis, not otherwise specified99 (25.5)17.3

Injury and poisoning860.0Traumatic pneumothorax27 (7.0)6.7
847.0Sprain/strain of neck21 (5.4)5.7

Figure 2 displays an association diagram that includes the new or poorly known associations with Charcot foot that are listed in Table 2. The specific diagnoses from the table are shown as red nodes. Some of the diagnoses associated with the red nodes are highly interrelated, which can be seen by the various “clusters” and the large number of edges (lines) connecting each of the related nodes to one another. Other red nodes are less connected. These clusters suggest that there are potentially multiple etiological pathways that could lead to Charcot foot.

To demonstrate the role diabetes may have had in the relationship between Charcot foot and the diagnoses reported in Tables 1 and 2, we also report on the frequency of these diagnoses among the population of Charcot foot patients where diabetes was either present (i.e., ICD codes 250.x) or absent. Table 3 shows the well-known associations from Table 1 stratified by diabetes, and Table 4 shows novel or poorly known associations from Table 2, also stratified by diabetes. Many of the associations (e.g., peripheral vascular disease) occurred more frequently in the diabetes cohort, whereas others were more evenly distributed among those with and without diabetes (e.g., asthma).

Disease categoryICD-9 codeICD-9 descriptionNumber (%) among those with diabetesNumber (%) among those without diabetes value for difference in proportions

Endocrine, nutritional and metabolic diseases, and immunity disorders278.00Obesity46 (16.3)6 (5.7)

Nervous system356.9Peripheral neuropathy143 (50.7)18 (17.0)

Blood and blood-forming organs280.0Anemia due to chronic blood loss14 (5.0)2 (1.9)
285.9Anemia, not otherwise specified88 (31.2)10 (9.4)
285.29Anemia of other chronic diseases8 (2.8)2 (1.9)
280.9Iron deficiency anemia29 (10.3)3 (2.8)

Skin and subcutaneous tissue682.6Cellulitis/abscess of leg80 (28.4)7 (6.6)
682.7Cellulitis/abscess of foot94 (33.3)4 (3.8)

Genitourinary system593.9Renal & ureteral disorder, not otherwise specified111 (39.4)8 (7.5)
585Chronic kidney disease63 (22.3)1 (0.9)

Musculoskeletal system and connective tissue714.0Rheumatoid arthritis28 (9.9)9 (8.5)
733.00Osteoporosis 36 (12.8)8 (7.5)
733.90Bone/cartilage disorder, not otherwise specified175 (62.1)58 (54.7)

Injury and poisoning996.81Complication of kidney transplant43 (15.2)1 (0.9)

External causes of injury and supplemental classificationV42.0Kidney transplant49 (17.4)1 (0.9)

Disease categoryICD-9 codeICD-9 descriptionNumber (%) among those with diabetesNumber (%) among those without diabetes value for difference in proportions

Endocrine, nutritional and metabolic diseases, and immunity disorders244.9Hypothyroidism34 (12.1)8 (7.5)0.28
272.0Hypercholesterolemia78 (27.7)4 (3.8)
272.4Hyperlipidemia163 (57.8)10 (9.4)<2.2 ×  
276.3Alkalosis11 (3.9)2 (1.9)0.51
276.5Volume depletion disorder75 (26.6)9 (8.5)

Sense organs366.16Senile nuclear cataract102 (36.2)8 (7.5)

Circulatory system401.1Benign essential hypertension103 (36.5)10 (9.4)
402.10Benign hypertensive disease without heart failure14 (5.0)1 (0.9)0.12
411.1Intermediate coronary syndrome38 (13.5)0 (0.0)
413.9Angina pectoris40 (14.2)1 (0.9)
414.0Coronary atherosclerosis18 (6.4)0 (0.0)0.017
414.00Coronary atherosclerosis, vessel type143 (50.7)12 (11.3)
414.01Coronary atherosclerosis, native artery89 (31.6)8 (7.5)
414.02Coronary atherosclerosis, autologous graft11 (3.9)0 (0.0)0.086
414.9Chronic ischemic heart disease49 (17.4)0 (0.0)
424.0Mitral valve disorder31 (11.0)1 (0.9)
425.4Cardiomyopathy, primary 42 (14.9)2 (1.9)
427.9Cardiac dysrhythmia143 (50.7)25 (23.6)
428.0Congestive heart failure78 (27.7)6 (5.7)
428.9Heart failure, not otherwise specified36 (12.8)2 (1.9)
429.9Heart disease, not otherwise specified38 (13.5)1 (0.9)
440.21Atherosclerosis of native arteries of the extremities with intermittent claudication89 (31.6)4 (3.8)
443.9Peripheral vascular disease91 (32.3)4 (3.8)

Respiratory system465.9Acute upper respiratory infection71 (25.2)22 (20.8)0.44
466.0Acute bronchitis60 (21.3)15 (14.2)0.15
473.9Chronic sinusitis26 (9.2)9 (8.5)0.98
477.9Allergic rhinitis30 (10.6)10 (9.4)0.87
493.90Asthma (without status asthmaticus)34 (12.1)12 (11.3)0.98
518.3Pulmonary eosinophilia104 (36.9)21 (19.8)

Digestive system530.81Esophageal reflux47 (16.7)17 (16.0)1
560.9Intestinal obstruction, not otherwise specified25 (8.9)5 (4.7)0.25

Genitourinary system596.9Bladder disorder, not otherwise specified9 (3.2)3 (2.8)1
607.84Impotence, organic origin29 (10.3)4 (3.8)0.065

Skin and subcutaneous tissue692.9Contact dermatitis 38 (13.5)14 (13.2)1

Musculoskeletal system and connective tissue724.2Lumbago83 (29.4)34 (32.1)0.70
729.1Myalgia/myositis, not otherwise specified45 (16.0)10 (9.4)0.14
715.90Osteoarthrosis, not otherwise specified74 (26.2)25 (23.6)0.69

Injury and poisoning860.0Traumatic pneumothorax25 (8.9)2 (1.9)0.029
847.0Sprain/strain of neck15 (5.3)6 (5.7)1

There were 111 association pairs that had a temporal value of less than 1.0 10−3, and nearly all of them were ones in which the other disorder significantly preceded Charcot foot, using our measure of temporality which was based on the binomial test. Only four diagnoses followed Charcot foot based on the binomial test for temporal statistical significance and these were (1) obstructive sleep apnea; ICD-9 327.23; association value 2.9 10−24; odds ratio 11.7; temporal value 4.9 10−4, (2) abnormal chest sounds; ICD-9 786.7; association value 3.2 10−19; odds ratio 10.3; temporal value 9.8 10−4, (3) phantom limb syndrome; ICD-9 353.6; association value < 5 10−324; odds ratio 137.5; temporal value 6.6 10−5, and (4) nonspecific symptoms involving the chest and respiratory system; ICD-9 786.9; association value 7.5 10−132; odds ratio 17.7; temporal value 9.0 10−6.

4. Discussion

Data mining has been increasingly applied to clinical data [19] and has revealed multiple associations that were previously unknown [2022]. For Charcot foot, there has been a paucity of literature that has used large databases to explore and control for comorbid conditions [5]. It is not surprising, then, that data mining techniques applied to a rare disorder such as Charcot foot might reveal new or poorly described associations. Indeed, a strength of our approach is the ability to identify associations that are not initially intuitive, without any preconceived bias. This is particularly important for Charcot foot where there is considerable debate over its etiology and pathogenesis [23, 24]. In fact, a recent 2011 international task force that was convened to summarize the Charcot foot evidence revealed few Charcot foot published studies describing its etiology and pathogenesis [25].

To better understand how Charcot foot might develop, it is helpful to discuss the findings from our study in the context of an etiological model proposed by Koeck et al. [24]. The mechanisms proposed in this model include (1) neurotrophic, (2) microtraumatic, and (3) neurovascular effects. The relationships between the model and the findings from our current study are speculative but could help shed light on a pathologic process that remains poorly understood.

Neurotrophic influences can involve local sensory loss and selective sympathetic denervation [24]. In line with this aspect of the model, we found lumbar disc disease to be associated with Charcot foot. Recent work in degenerative disc disease has implicated the role of neurotrophic factors [26]. We also found hypothyroidism to be associated with Charcot foot. Thyroid hormone deficiencies associated with brain-derived neurotrophic factor have been established [27]. Another association with Charcot foot was bladder disorder. Disruptions in afferent regulation of the bladder, as well as with brain-derived neurotrophic factor, have been implicated in bladder dysfunction [28, 29]. Overall, we found that 12 of 13 genitourinary conditions preceded the diagnosis of Charcot foot. As far as we are aware, each of the above associations is novel for Charcot foot.

Regardless of what model predominates, what is known about Charcot foot etiology is that it involves a stage of proinflammatory cytokine activity including elevated tissue necrosis factor (TNF) [30] and receptor activator nuclear factor K ligand (RANKL) [31]. This proinflammatory cytokine activity has also been described as a part of the etiology model [24, 31]. In our study, we found that coronary artery disease (CAD) demonstrated an odds ratio of 18.6 with Charcot foot. It is well established that CAD is a proinflammatory state [32]. Relatedly, we also found hyperlipidemia to be associated with Charcot foot. Treatment of CAD with 3-hydroxy-3-methylglutaryl-coenzyme A (HMGCoA) reductase inhibitors (i.e., “statins”) has also been described to have an anti-inflammatory effect [33, 34]. CAD has not been described in prior Charcot foot multivariate models [4, 5]. It could be that age and diabetes status are collinear with CAD using a multivariate approach.

Other associations that do not fit neatly into any model are also worth mentioning including the association between Charcot foot and alkalosis, which had an odds ratio of 11.2. Transient perturbations in parathyroid hormone, calcium, and phosphorus [35, 36] may explain why alkalosis preceded Charcot foot 100% of the time. We also found pulmonary eosinophilia to be associated with Charcot foot, and this relationship had the fourth smallest value for all the temporal associations we described. Treatment for pulmonary eosinophilia with inhaled [37] or oral steroids [38] and subsequent bone mineral density and bone metabolism changes could partly explain why this disorder was associated with Charcot foot. Furthermore, we found an association with esophageal reflux and Charcot foot. This association could be due to concurrent use of proton pump inhibitors that have been described to reduce bone mineral density [39]. As far as the authors are aware, the associations described above are also novel in Charcot foot.

While we describe several new associations with Charcot foot, our analysis also confirmed other known risk factors for Charcot foot. For example, we found associations of Charcot foot with obesity, peripheral neuropathy, decreased bone mineral density, and a history of pancreas and/or kidney transplant surgery [4]. Our findings also agree with Stuck and colleagues in finding associations with renal failure, rheumatoid arthritis, and anemia [5].

Our approach did have limitations. One of the main limitations is that we chose to discuss those new associations that fit our conceptual model. Our discussion was not exhaustive of the 111 significant temporal associations we described using our conservative acceptance of value thresholds. Thus, all of the associations are provided as a Supplementary Table. The data presented here can be used for hypothesis generation, not to make any firm conclusions about causality. The more interesting associations could potentially be followed up by a chart review of the clinical records to ascertain more details about the relationships, as has been done with other interesting associations [20].

It is also important to note that our analysis focused on ICD-9 codes, which have potential limitations in not only accuracy but also specificity. Some codes are nonspecific and therefore do not provide enough detail to fully understand what they represent. For example, in our population the code 733.90 (bone/cartilage disorder, not otherwise specified) was present in 60.1% of the patients with Charcot foot. From the code itself we are unable to determine if these represent a different condition from Charcot foot. Using the temporal data to provide additional insight, we can observe that in the majority of patients the bone/cartilage disorder preceded the diagnoses of Charcot foot, suggesting that 733.90 could have been an initial misdiagnosis of Charcot foot rather than a distinct clinical entity. We also found that in some cases groups of codes should be considered collectively rather than separately. This is the case for the diabetes codes as well as those potentially representing peripheral neuropathy. Additionally, some of the variation in the use of codes could represent differences in coding among disciplines, since the codes used in our analysis were assigned not just by podiatrists but by clinicians across the entire health system.

Another limitation of our study is that medications were not included in our analysis, and yet these could influence both the development and prevention of Charcot foot. Future work should include additional clinical features to develop a richer understanding of Charcot foot and its etiology. Further, it is possible that some patients with true Charcot foot continued to remain misdiagnosed (with gout or osteoarthropathy, e.g.) and thus never received the diagnosis code of 713.5. Future studies should examine such cases for potential misdiagnoses of Charcot foot, especially because earlier detection could provide substantial clinical benefit [40].

A primary purpose for presenting our results was to help generate hypotheses for further testing and confirmation and to better shed light on Charcot foot, for which the etiology remains poorly described. Future work should seek to confirm the validity of the associations found in our exploratory analysis. One next step would be to conduct a detailed chart review of the cases from the electronic health record to better understand the context in which the associations occurred and even to confirm the diagnosis of Charcot foot. Incorporating other data sources, such as from the Veterans Health Administration or the Centers for Medicare and Medicaid Services, could also provide additional opportunities for analysis of this rare disorder. It will also be important to understand in more detail how coding practices (and variability) could influence our findings.

5. Conclusion

We present here a novel application of an approach to mining large databases, such as those available through the electronic health record, for elucidating new associations with Charcot foot that might otherwise go unnoticed. We also uniquely describe the temporal relationships between these various diseases and Charcot foot. This approach has potential for generating new hypotheses regarding the cause of Charcot foot. At a minimum, we hope to create a heightened awareness of Charcot foot in patients with associated conditions to help identify these at-risk patients sooner. Also, if we can elucidate the pathway of Charcot foot, it would hopefully lead to better treatment modalities including drug development.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Supplementary Materials

The Supplementary Material contains a table in the Microsoft Excel format that lists all of the associations we found with Charcot foot. These include associations that were found but not discussed in the manuscript. Additional information about how to interpret the results in the supplementary table are also provided within the file.

  1. Supplementary Material


  1. L. A. Lavery, D. G. Armstrong, R. P. Wunderlich, J. Tredwell, and A. J. M. Boulton, “Diabetic foot syndrome: evaluating the prevalence and incidence of foot pathology in Mexican Americans and non-Hispanic whites from a diabetes disease management cohort,” Diabetes Care, vol. 26, no. 5, pp. 1435–1438, 2003. View at: Publisher Site | Google Scholar
  2. J. Fabrin, K. Larsen, and P. E. Holstein, “Long-term follow-up in diabetic charcot feet with spontaneous onset,” Diabetes Care, vol. 23, no. 6, pp. 796–800, 2000. View at: Google Scholar
  3. S. Sinha, C. S. Munichoodappa, and G. P. Kozak, “Neuro-arthropathy (Charcot joints) in diabetes mellitus (clinical study of 101 cases),” Medicine, vol. 51, no. 3, pp. 191–210, 1972. View at: Google Scholar
  4. R. G. Frykberg and R. Belczyk, “Epidemiology of the Charcot foot,” Clinics in Podiatric Medicine and Surgery, vol. 25, no. 1, pp. 17–28, 2008. View at: Publisher Site | Google Scholar
  5. R. M. Stuck, M. Sohn, E. Budiman-Mak, T. A. Lee, and K. B. Weiss, “Charcot arthropathy risk elevation in the obese diabetic population,” The American Journal of Medicine, vol. 121, no. 11, pp. 1008–1014, 2008. View at: Publisher Site | Google Scholar
  6. G. A. Matricali, B. Bammens, D. Kuypers, M. Flour, and C. Mathieu, “High rate of Charcot foot attacks early after simultaneous pancreas-kidney transplantation,” Transplantation, vol. 83, no. 2, pp. 245–246, 2007. View at: Publisher Site | Google Scholar
  7. R. G. Ramani and S. G. Jacob, “Data mining in clinical data sets: a review,” International Journal of Applied Information Systems, vol. 4, no. 6, pp. 15–26, 2012. View at: Google Scholar
  8. P. B. Jensen, L. J. Jensen, and S. Brunak, “Mining electronic health records: towards better research applications and clinical care,” Nature Reviews Genetics, vol. 13, no. 6, pp. 395–405, 2012. View at: Publisher Site | Google Scholar
  9. I. M. Mullins, M. S. Siadaty, J. Lyman et al., “Data mining and clinical data repositories: insights from a 667,000 patient data set,” Computers in Biology and Medicine, vol. 36, no. 12, pp. 1351–1377, 2006. View at: Publisher Site | Google Scholar
  10. D. A. Hanauer and N. Ramakrishnan, “Modeling temporal relationships in large scale clinical associations,” Journal of the American Medical Informatics Association, vol. 20, no. 2, pp. 332–341, 2013. View at: Publisher Site | Google Scholar
  11. D. A. Hanauer, D. R. Rhodes, and A. M. Chinnaiyan, “Exploring clinical associations using “-omics” based enrichment analyses,” PLoS ONE, vol. 4, no. 4, Article ID e5203, 2009. View at: Publisher Site | Google Scholar
  12. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site | Google Scholar
  13. L. N. McEwen, K. R. Ylitalo, W. H. Herman, and J. S. Wrobel, “Prevalence and risk factors for diabetes-related foot complications in Translating Research Into Action for Diabetes (TRIAD),” Journal of Diabetes and Its Complications, vol. 27, no. 6, pp. 588–592, 2013. View at: Publisher Site | Google Scholar
  14. M. S. Pinzur, “Benchmark analysis of diabetic patients with neuropathic (Charcot) foot deformity,” Foot and Ankle International, vol. 20, no. 9, pp. 564–567, 1999. View at: Google Scholar
  15. M. Sohn, T. A. Lee, R. M. Stuck, R. G. Frykberg, and E. Budiman-Mak, “Mortality risk of charcot arthropathy compared with that of diabetic foot ulcer and diabetes alone,” Diabetes Care, vol. 32, no. 5, pp. 816–821, 2009. View at: Publisher Site | Google Scholar
  16. M. W. Sohn, R. M. Stuck, M. Pinzur, T. A. Lee, and E. Budiman-Mak, “Lower-extremity amputation risk after Charcot arthropathy and diabetic foot ulcer,” Diabetes Care, vol. 33, no. 1, pp. 98–100, 2010. View at: Publisher Site | Google Scholar
  17. E. T. Rhodes, L. M. B. Laffel, T. V. Gonzalez, and D. S. Ludwig, “Accuracy of administrative coding for type 2 diabetes in children, adolescents, and young adults,” Diabetes Care, vol. 30, no. 1, pp. 141–143, 2007. View at: Publisher Site | Google Scholar
  18. G. Botek, M. A. Anderson, and R. Taylor, “Charcot neuroarthropathy: an often overlooked complication of diabetes,” Cleveland Clinic Journal of Medicine, vol. 77, no. 9, pp. 593–599, 2010. View at: Publisher Site | Google Scholar
  19. J. Iavindrasana, G. Cohen, A. Depeursinge, H. Müller, R. Meyer, and A. Geissbuhler, “Clinical data mining: a review,” Yearbook of Medical Informatics, pp. 121–133, 2009. View at: Google Scholar
  20. D. A. Hanauer, N. Ramakrishnan, and L. S. Seyfried, “Describing the relationship between cat bites and human depression using data from an electronic health record,” PLoS ONE, vol. 8, no. 8, Article ID e70585, 2013. View at: Publisher Site | Google Scholar
  21. N. P. Tatonetti, P. P. Ye, R. Daneshjou, and R. B. Altman, “Data-driven prediction of drug effects and interactions,” Science Translational Medicine, vol. 4, no. 125, Article ID 125ra31, 2012. View at: Publisher Site | Google Scholar
  22. A. Wright, E. S. Chen, and F. L. Maloney, “An automated technique for identifying associations between medications, laboratory results and problems,” Journal of Biomedical Informatics, vol. 43, no. 6, pp. 891–901, 2010. View at: Publisher Site | Google Scholar
  23. E. Chantelau and G. J. Onvlee, “Charcot foot in diabetes: farewell to the neurotrophic theory,” Hormone and Metabolic Research, vol. 38, no. 6, pp. 361–367, 2006. View at: Publisher Site | Google Scholar
  24. F. X. Koeck, V. Bobrik, A. Fassold, J. Grifka, S. Kessler, and R. H. Straub, “Marked loss of sympathetic nerve fibers in chronic charcot foot of diabetic origin compared to ankle joint osteoarthritis,” Journal of Orthopaedic Research, vol. 27, no. 6, pp. 736–741, 2009. View at: Publisher Site | Google Scholar
  25. L. C. Rogers, R. G. Frykberg, D. G. Armstrong et al., “The Charcot foot in diabetes,” Diabetes Care, vol. 34, no. 9, pp. 2123–2129, 2011. View at: Publisher Site | Google Scholar
  26. S. E. Navone, G. Marfia, L. Canzi et al., “Expression of neural and neurotrophic markers in nucleus pulposus cells isolated from degenerated intervertebral disc,” Journal of Orthopaedic Research, vol. 30, no. 9, pp. 1470–1477, 2012. View at: Publisher Site | Google Scholar
  27. M. E. Gilbert and S. M. Lasley, “Developmental thyroid hormone insufficiency and brain development: a role for brain-derived neurotrophic factor (BDNF)?” Neuroscience, vol. 239, pp. 253–270, 2013. View at: Publisher Site | Google Scholar
  28. T. Antunes-Lopes, R. Pinto, S. C. Barros et al., “Urinary neurotrophic factors in healthy individuals and patients with overactive bladder,” Journal of Urology, vol. 189, no. 1, pp. 359–365, 2013. View at: Publisher Site | Google Scholar
  29. W. C. de Groat and N. Yoshimura, “Afferent nerve regulation of bladder function in health and disease,” Handbook of Experimental Pharmacology, vol. 194, pp. 91–138, 2009. View at: Publisher Site | Google Scholar
  30. G. Mabilleau, N. Petrova, M. E. Edmonds, and A. Sabokbar, “Number of circulating CD 14-positive cells and the serum levels of TNF-α are raised in acute charcot foot,” Diabetes Care, vol. 34, no. 3, article e33, 2011. View at: Publisher Site | Google Scholar
  31. A. Ndip, A. Williams, E. B. Jude et al., “The RANKL/RANK/OPG signaling pathway mediates medial arterial calcification in diabetic charcot neuroarthropathy,” Diabetes, vol. 60, no. 8, pp. 2187–2196, 2011. View at: Publisher Site | Google Scholar
  32. S. Kaptoge, E. Di Angelantonio, G. Lowe et al., “C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis,” The Lancet, vol. 375, no. 9709, pp. 132–140, 2010. View at: Publisher Site | Google Scholar
  33. E. Mandosi, M. Fallarino, A. Gatti et al., “Atorvastatin downregulates monocyte CD36 expression, nuclear NFκB and TNFα levels in type 2 diabetes,” Journal of Atherosclerosis and Thrombosis, vol. 17, no. 6, pp. 539–545, 2010. View at: Google Scholar
  34. K. Yanagi, T. Monden, S. Ikeda, M. Matsumura, and K. Kasai, “A crossover study of rosuvastatin and pitavastatin in patients with type 2 diabetes,” Advances in Therapy, vol. 28, no. 2, pp. 160–171, 2011. View at: Publisher Site | Google Scholar
  35. T. Hamano, N. Fujii, I. Matsui et al., “Guideline-practice gap in the management of predialysis chronic kidney disease mineral bone disorder in Japan,” Therapeutic Apheresis and Dialysis, vol. 15, supplement 1, pp. 2–8, 2011. View at: Publisher Site | Google Scholar
  36. E. M. Omron, “Comparative quantitative acid-base analysis in coronary artery bypass, severe sepsis, and diabetic ketoacidosis,” Journal of Intensive Care Medicine, vol. 20, no. 6, pp. 317–326, 2005. View at: Google Scholar
  37. P. R. Ebeling, B. Erbas, J. L. Hopper, J. D. Wark, and A. R. Rubinfeld, “Bone mineral density and bone turnover in asthmatics treated with long-term inhaled or oral glucocorticoids,” Journal of Bone and Mineral Research, vol. 13, no. 8, pp. 1283–1289, 1998. View at: Publisher Site | Google Scholar
  38. F. L. Dente, E. Bacci, M. L. Bartoli et al., “Effects of oral prednisone on sputum eosinophils and cytokines in patients with severe refractory asthma,” Annals of Allergy, Asthma and Immunology, vol. 104, no. 6, pp. 464–470, 2010. View at: Publisher Site | Google Scholar
  39. L. E. Targownik, W. D. Leslie, K. S. Davison et al., “The relationship between proton pump inhibitor use and longitudinal change in bone mineral density: a population-based study [corrected] from the Canadian Multicentre Osteoporosis study (CaMos),” The American Journal of Gastroenterology, vol. 107, no. 9, pp. 1361–1369, 2012. View at: Publisher Site | Google Scholar
  40. D. K. Wukich, W. Sung, S. A. Wipf, and D. G. Armstrong, “The consequences of complacency: managing the effects of unrecognized Charcot feet,” Diabetic Medicine, vol. 28, no. 2, pp. 195–198, 2011. View at: Publisher Site | Google Scholar

Copyright © 2014 Michael E. Munson et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.