Next Article in Journal
Transport of L-Arginine Related Cardiovascular Risk Markers
Next Article in Special Issue
The Acute Effect of Exercise on Arterial Stiffness in Healthy Subjects: A Meta-Analysis
Previous Article in Journal
Chronic Rhinosinusitis: MALDI-TOF Mass Spectrometry Microbiological Diagnosis and Electron Microscopy Analysis; Experience of the 2nd Otorhinolaryngology Clinic of Cluj-Napoca, Romania
Previous Article in Special Issue
Infrapopliteal Artery Occlusive Disease: An Overview of Vessel Preparation and Treatment Options
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 9
Issue 12
10.3390/jcm9123974
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fibrosis Distinguishes Critical Limb Ischemia Patients from Claudicants in a Transcriptomic and Histologic Analysis

by
Guangzhi Cong
1,2,
Xiangdong Cui
1,
Ricardo Ferrari
1,
Iraklis I. Pipinos
3,4,
George P. Casale
3,
Ansuman Chattopadhyay
4 and
Ulka Sachdev
1,*
1
Department of Surgery, University of Pittsburgh Medical Centre, Pittsburgh, PA 15217, USA
2
Department of Surgery, University of Nebraska at Medical Center, Omaha, NE 68198, USA
3
Department of Surgery and VA Research Service, VA Nebraska-Western Iowa Health Care System, Omaha, NE 68198, USA
4
Molecular Biology Information Service, Health Sciences Library System University of Pittsburgh, Pittsburgh, PA 15261, USA
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(12), 3974; https://doi.org/10.3390/jcm9123974
Submission received: 13 October 2020 / Revised: 24 November 2020 / Accepted: 26 November 2020 / Published: 8 December 2020
(This article belongs to the Special Issue Peripheral Artery Disease: From Diagnosis to Treatment)
Browse Figures
Versions Notes

Abstract

:
Most patients with critical limb ischemia (CLI) from peripheral arterial disease (PAD) do not have antecedent intermittent claudication (IC). We hypothesized that transcriptomic analysis would identify CLI-specific pathways, particularly in regards to fibrosis. Derivation cohort data from muscle biopsies in PAD and non-PAD (controls) was obtained from the Gene Expression Omnibus (GSE120642). Transcriptomic analysis indicated CLI patients (N = 16) had a unique gene expression profile, when compared with non-PAD controls (N = 15) and IC (N = 20). Ninety-eight genes differed between controls and IC, 2489 genes differed between CLI and controls, and 2783 genes differed between CLI and IC patients. Pathway enrichment analysis showed that pathways associated with TGFβ, collagen deposition, and VEGF signaling were enriched in CLI but not IC. Receiver operating curve (ROC) analysis of nine fibrosis core gene expression revealed the areas under the ROC (AUC) were all >0.75 for CLI. Furthermore, the fibrosis area (AUC = 0.81) and % fibrosis (AUC = 0.87) in validation cohort validated the fibrosis discrimination CLI from IC and control (all n = 12). In conclusion, transcriptomic analysis identified fibrosis pathways, including those involving TGFβ, as a novel gene expression feature for CLI but not IC. Fibrosis is an important characteristic of CLI, which we confirmed histologically, and may be a target for novel therapies in PAD.

1. Introduction

Peripheral artery disease (PAD) is a manifestation of systemic atherosclerosis with associated adverse cardiovascular mortality and risk of amputation [1]. According to the latest version of Heart Disease and Stroke Statistics reported by the American Heart Association, 6.5 million Americans aged ≥40 years (5.5%) are estimated to be afflicted with PAD [2]. Among them, 2 million individuals are currently suffering from critical limb ischemia (CLI) [3,4], which represents the most severe form of PAD. CLI is traditionally viewed as a progressive continuum from intermittent claudication (IC). However, most IC patients do not progress to CLI [5,6], which suggests that CLI may represent a unique feature of PAD with distinct etiology.
CLI is diagnosed based on clinical presentation of either tissue loss or rest pain, coupled with hemodynamic measurements of ankle bronchial index (ABI) or toe pressures. The anatomic location of vascular disease is confirmed with traditional or computed tomographic, with occasional use of magnetic resonance imaging (MRI) [7]. However, as stated by ‘’Critical Limb Ischemia: An Expert Statement’’, approximately 30% of patients with CLI have a near-normal or normal ABI (>0.90) [7], possibly due to vessel calcification [8,9]. This suggests there is significant heterogeneity of pathogenesis among PAD patients. The use of transcriptomic data provides promising workflows for personalized diagnosis and treatment of cardiovascular diseases in several landmark studies [10,11,12,13]. Because transcriptome-wide profiling enables us to detect novel, un-annotated genes, sequence-level alterations, and alternative splicing events, it offers greater information about tissue than microarrays. We, therefore, hypothesized that transcriptomic analysis will identify CLI-specific pathways, particularly in regards to fibrosis, that will distinguish critical limb ischemia patients from those with intermittent claudication.
To address this hypothesis, we performed an analysis of a publicly available and recently published RNA-sequencing database of PAD cohorts. Following this, we evaluated the diagnostic performance of fibrosis pathway gene signatures for distinguishing CLI from IC and validated the findings using fibrosis staining of skeletal muscle biopsies of PAD patients.

2. Experimental Section

2.1. Study Design

An observational, cross-sectional study was conducted to explore whether fibrosis identified by transcriptomic and histology assay can distinguish CLI patients from IC. We reasoned that RNA-sequencing could provide insight into a pathogenesis pathway and could potentially be a complementary diagnostic tool. To interpret transcriptional differences seen in patients with and without PAD, transcriptome profiles of PAD patients and non-PAD controls were obtained from the Gene Expression Omnibus (GEO; GSE120642) resulting from studies by Ryan et al. [14]. First, transcriptomic profiles that distinguished CLI patients from those with claudication and controls were identified. We then validated our findings using histologic techniques to confirm fibrosis (Figure 1).

2.2. Retrieval of Derivation Cohort Datasets of Patients from the GEO

A derivation cohort of PAD patients and control subjects without PAD with RNA-sequencing-based profiles were acquired from the public functional genomics data repository GEO and analyzed [14]. In this published cohort, gastrocnemius muscle was obtained from PAD patients and non-PAD controls for whole transcriptome sequencing. The diagnosis of IC and CLI was based on clinical characteristics and ankle-brachial index [14]. The diagnostic criteria and participant’s demographic and medical information were explicitly illustrated in the paper [14].

2.3. Bioinformatics Analysis

Paired FASTQ files storing sequencing data from gastrocnemius muscle was obtained [14]. RNA-seq data analysis was done using CLC Genomics Workbench v20.0 (Qiagen Digital Insights, Aarhus, Denmark)—a software suite that provides a variety of tools for next-generation sequence analysis. The RNA-seq analysis workflow, as described in the CLC Genomics manual, was followed [15]. Briefly, single-end RNA-seq reads, obtained in FASTQ format, were checked for sequence quality using the tool “QC for Sequencing Reads” available under the toolbox item “Prepare Sequencing Data.” Reads with Phred quality score > 20 were then aligned to the human reference genome GRCh38 assembly sequence using the default parameters provided by the “RNA-seq Analysis” tool available under the “RNA-Seq Analysis” toolbox item. Genes and transcripts annotations provided by Ensembl (release 92) were used. Around 94% of reads from each sample were mapped into the reference genome [15,16].
To explore transcriptomics features among subtype of PAD, we conducted cluster analysis and dimensional reduction analysis by principal component analysis (PCA), uniform manifold approximation and projection (UMAP), and hierarchical cluster analysis. Empirical analysis of differentially expressed reads was conducted to identify differentially expressed genes (DEGs) with default parameters. DEGs were filtered according to pairwise comparisons for biological and statistical significance (fold change (FC) > 2, false discovery rate (FDR) < 0.05) [16].

2.4. Pathway Enrichment Analysis

To explore the functional enrichment of biological pathways and upstream molecules that are enriched in DEGs, pathway enrichment analysis and visualization was performed using IPA (QIAGEN Inc., Redwood City, CA, USA). Briefly, IPA regulation z-score algorithm will be used to determine the activation state of canonical pathways and upstream regulators using the general molecular network implemented in the ingenuity pathway analysis (IPA) knowledge base. To enhance the stringency of our analysis, only z-score >2 or <−2 (two folds) is considered significant as positive or negative [17].

2.5. Fibrosis Core Gene Expression Analysis

Transforming Growth Factor-β, collagen deposition, and vascular endothelial growth factor (VEGF) signaling pathway are the main contributors to the formation of fibrosis [18]. Therefore, the fibrosis core gene expression in skeletal muscle transcriptomics, including TGFβ1, TGFβ2, TGFβ3, COL1A1, COL1A2, VEGFA, VEGFB, VEGFC, and VEGFD were counted specifically as total counts and Transcripts Per Kilobase Million (TPM). The calculation and visualization were run by CLC Genomics Workbench 20.0 (Qiagen Digital Insight, Aarhaus, Denmark) and Partek® Flow® software, version 7.0 Copyright©; 2020 Partek Inc., St. Louis, MO, USA.

2.6. Validation Cohort Data

Validation cohort sets of gastrocnemius muscle biopsies slides were obtained from the University of Nebraska Medical Center (UNMC). The diagnosis of IC and CLI was also based on clinical characteristics and ankle-brachial index. Biopsy techniques and tissue processing were performed as described previously using a 6 mm Bergstrom needle under anesthesia [19]. Specimens were paraffin-fixed and sequentially sectioned. Masson’s trichrome staining was performed to identify areas of fibrosis.

2.7. Ethical Compliance

The study complied with all relevant ethical regulations consistent with the policies of the Omaha VA Medical Center, the University of Nebraska Medical Center, and the University of Pittsburgh Medical Center. The study guiding the procurement of muscle tissue was approved by the IRB of the Omaha VA Medical Center in accordance with the Declaration of Helsinki (ID# 00086, Protocol #01), and all patients provided written informed consent.

2.8. Fibrosis Evaluation in Validation Cohort

The fibrosis levels of gastrocnemius muscle biopsies were determined by Masson’s trichrome staining [20]. Photographs were taken with Olympus Provus-1 fluorescence microscope using a 20X objective. To quantify the extent of fibrosis, stained sections were analyzed using Image J analysis program [21] after color separation of the figures [22]. Percent and area of fibrosis were determined based on the number of pixels per image that were positive for fibrosis-specific staining (blue) in each section, represented as a percentage of the total area. Images were analyzed in a blinded fashion.

2.9. ROC Analysis Both in Derivation and Validation Cohort

To determine the specificity and sensitivity of the fibrosis pathway, receiver operating curve (ROC) analysis of the expression of nine core fibrosis genes in a derivation cohort was performed. Then, ROC analysis of percentage and area of fibrosis was performed in the validation. Specificity was designated on the x-axis, and sensitivity was designated on the y-axis. An area under the curve (AUC) = 1.0 suggests an ideal test, while an AUC < 0.5 indicates a non-diagnostic test.

2.10. Statistics for Clinical Data

Categorical variables were compared between non-PAD controls, IC, and CLI group using Χ2 test. Continuous variables were compared between HA, IC, and CLI group using single one-way analysis of variance (ANOVA) followed by Bonferroni t-test for data involving more than two groups. All p values reported are considered to be descriptive, and no adjustments for multiple comparisons were made. Statistical software package R-3.4.3 [23], The R Foundation) and Empower-Stats (X&Y Solutions, Inc., Boston, MA, USA) were used for the demographic analyses.

3. Results

3.1. Subjects’ Age and ABI Derivation and Validation Cohort

Table 1 demonstrates clinical characteristics of the PAD and control subjects from the derivation and validation cohorts.
The results of derivation cohort cited to published paper in JCI insight [14].

3.2. CLI Patients Have a Unique Gene Expression Profile Compared to Claudicants and Non-PAD Controls Subjects

To explore the transcriptomic features in these two different clinical manifestations of PAD, we conducted a dimensionality reduction by PCA, UMAP, and hierarchical cluster analysis. As shown in Figure 2A,B, PCA and UMAP revealed that CLI samples clustered in a distinct pattern from both control and IC patients. Meanwhile, control and IC patients’ transcriptomic profiles overlapped significantly. By hierarchical cluster analysis and differential expression analysis (DEA), we observed that 2783 genes out of 29,191 total genes (2012 up-regulated, and 771 down-regulated) were differentially expressed in CLI when compared to IC (Figure 2C,D), and 2489 genes out of 29,191 total genes (1765 up-regulated, and 724 down-regulated) were differentially expressed in CLI when compared to control (Figure 2D). However, there are only 98 genes out of 29,191 total genes (60 up-regulated, and 38 down-regulated) that were differentially expressed in IC when compared to control (Figure 2C,D).

3.3. Differentially Expressed Genes (DEG) and Canonical Pathway Analysis Distinguishes CLI Patients

A nonbiased analysis of the canonical pathways enriched by the differentially expressed genes revealed that oxidative phosphorylation, mitochondrial dysfunction, sirtuin signaling pathway, TCA cycle II, and fibrosis signaling pathways were distinct in CLI when compared with IC and controls (Figure 3A,B,D,E). Among these five canonical pathways, the fibrosis pathway and sirtuin signaling were activated, while the other pathways were inhibited. By overlaying these pathways, we found that oxidative phosphorylation, mitochondrial dysfunction, sirtuin signaling pathway, and TCA cycle II were connected to each other and shared common genes. On the contrary, the fibrosis pathway did not share common genes with the other pathways (Figure 3C,F). Therefore, we selected fibrosis to determine whether it was a candidate pathway differentiating CLI from IC and controls. Furthermore, there is only activation of actin cytoskeleton signaling with a small proportion of genes involved as showed in Figure 3G. Canonical pathway heatmap of three comparisons confirms that fibrosis pathway was activated in CLI but not IC when compared with non-PAD patients (Figure 3H).
To illustrate the potential differences in the canonical fibrosis pathway in CLI patients, we plotted the core molecular networks that may constitute the fibrosis pathway by IPA software, which compares Genbank Accession number/expression information with a proprietary protein-interaction database to establish the probability of a given signaling. When comparing CLI samples to controls, gene clusters that participate in TGFβ signaling, collagen deposition, VEGF pathways, and DNA synthesis were either widely activated or predicted to be activated (Figure 4A). Similarly, TGFβ signaling, collagen deposition, and VEGF pathways were activated in CLI compared to IC (Figure 4B). In contrast, these pathways did not differ when IC was compared to control patients (Figure 4C).

3.4. Fibrosis Pathway Core Genes Expression in Transcriptomics and ROC in the Derivation Cohort

TGFβ signaling, collagen synthesis, and deposition and VEGF signaling are key components of healing and fibrosis and likely coordinate with each other in the formation of skeletal muscle fibrosis [24]. The gene expression of core genes of TGFβ signaling (TGFβ1, TGFβ2, and TGFβ3), collagen synthesis and deposition (COL1A1 and COL1A2), and VEGF signaling (VEGF A, VEGF B, VEGF C, and VEGF D) were therefore calculated and expressed as total counts and TPM by CLC Genomics Workbench 20.0. Most of the core genes were activated when comparing CLI to IC and controls (Figure 5A). We then evaluated the diagnostic performance of these nine core genes to discriminate CLI from IC and controls using AUC. As shown in Figure 5B, the TGFβ signaling (TGFβ1, TGFβ2, and TGFβ3), collagen synthesis (COL1A1 and COL1A2) and VEGF signaling (VEGFA, VEGFB, EGFC, and VEGFD) yielded good discriminatory power (all AUC > 0.75).

3.5. Fibrosis Quantification and ROC in the Validation Cohort

To validate whether fibrosis was a prominent feature of CLI, we used Masson’s trichrome staining to detect fibrosis within muscle samples from CLI, IC, and control patients without PAD. Then, we accessed the area and percentage of fibrosis in patient’s biopsies slides and explored their discriminatory power for CLI. As depicted in Figure 6A,B, the fibrosis area and percentage were significantly increased in CLI compared to IC and control (all p < 0.05). Furthermore, the AUC of the fibrosis area and percentage were over 75%, which indicated that the discriminatory power of fibrosis is validated histologically (Figure 6C).

4. Discussion

CLI uncommonly develops as a progression of IC [5,6]. Instead, patients tend to present either with CLI or IC, suggesting that CLI has distinct pathogenesis from IC despite a similar ischemic etiology [24,25]. We used a transcriptomics approach to determine some of the identifying features of CLI and compared them to IC and normally perfused muscle. Previous PAD transcriptomic studies have not specifically defined gene expression levels and pathway analysis in relation to IC or CLI states. Sample sizes were also limited [26,27,28]. In this cross-sectional study, we found that (1) CLI patients present with a unique transcriptomic profile that is different from IC patients and non-PAD control; (2) fibrosis pathways that highlight TGFβ collagen deposition and VEGF distinguish CLI from IC and non-PAD control using transcriptomic data and IPA; and (3) increased fibrosis could be identified in CLI patients, validating the data we identified from the transcriptomic analysis (Figure 7). Our results provide a systems-level, data-driven determination of a complex network of molecular responses underpinning PAD, which was confirmed histologically. The data also provide important granularity for multiple distinct biological pathways, the relationships among these pathways, and their potential application in terms of the precision diagnosis and treatment of PAD.
Our findings are consistent with several proposed mechanisms of PAD progression, and confirm their importance in a relevant human study [29,30,31]. For example, TGFβ produced by vascular smooth muscle cells predicts fibrosis in the gastrocnemius of patients with PAD [30]. However, our study also determined other areas of difference between the groups that may be relevant. For example, our results suggest that there is inhibition of mitochondrial function in CLI. In contrast to claudicating and non-PAD control subjects, the pathways governing energy metabolism, including oxidative phosphorylation, mitochondrial dysfunction, sirtuin signaling, and TCA cycle II, were significantly suppressed. This is consistent with results in the study by Ryan et al. showing that skeletal muscle from CLI patients harbor a unique muscle mitochondriopathy when compared with IC patients [14]. Indeed, the RNA-sequencing data from the referenced study was made publicly available and was analyzed for this report. Muscle mitochondrial dysfunction has been demonstrated to be a key event in PAD [32]. Decreased oxidative capacity due to mitochondrial respiratory chain impairment is associated with an increased release of reactive oxygen species, leading to a reduction of calcium retention, and apoptosis [32]. However, detailing the intricacy of mitochondrial dysfunction is time-consuming and expensive, which may limit its utility in precision diagnosis and treatment [33,34].
Most notably, we identified fibrosis signaling as a critical event in CLI both by transcriptomic and histologic analysis. Fibrosis is the abnormal deposition of extracellular matrix, which is largely driven by activation of TGFβ signaling, collagen deposition, and VEGF signaling [18]. As shown in the fibrosis pathway enrichment analysis, these principal fibrosis component pathways were broadly activated. Ha et al. [30] reported that TGFβ signal was significantly activated in CLI patients and was generated in micro-vessels in response to tissue hypoxia. Our imaging was consistent with this study, and fibrosis in our images was often seen around vessels. Our observation may indicate a mechanism central to the association between ischemia and PAD progression [35] that relates to TGFβ activity; collagen deposition; and, ultimately, limb dysfunction.
Our results support the promising application of fibrosis signaling to distinguish CLI from IC. ABI is often the first-line noninvasive measure of limb ischemia in the clinic, but may be inaccurate if the blood vessels are very calcified [36]. Thus, there is a notable heterogeneity of ABI in patients with CLI [37]. Our results suggest that fibrosis may help distinguish the severity of PAD in a clinical setting, or can help guide future studies to mitigate disease. Fibrosis can be detected noninvasively using multiple imaging techniques including MRI [38]. We present that transcriptomic data may be able to identify activation of fibrotic pathways in patients with PAD. Transcriptomics has shown some promise in differentiating myopathies and muscular dystrophies in other clinical settings as described by Cummings [39]. Furthermore, transcriptomics of cardiac biopsies disclosed differences in patients with or without heart failure with preserved ejection fraction [40]. These studies consistently portrayed the potential of transcriptomic analysis to identify severity and pathogenesis in a variety of disease states, although its use in PAD is not optimized at this time.

Strengths and Limitations

A major strength of our study is that all the data were derived from skeletal muscle biopsies of PAD and non-PAD patients instead of animal models, which facilitates translatability. By using transcriptome network analysis, we were able to investigate our aims in an unbiased manner, leading to the discovery of new findings while also confirming and expanding upon anticipated results. By using a bioinformatics network analysis validated histologically, we were able to distill a large dataset into a coherent, statistically rigorous, and biologically meaningful finding. However, our study also has important limitations. First, the study had a cross-sectional design, and specific interventions could not be studied. In the future, the application of transcriptomics can be performed before and after treatments including revascularization. Secondly, the study population we utilized is considered small from a clinical perspective with limited demographic and commodities information, which limits its generalization. Future clinical trials with a bigger sample size are needed to determine the chronology of budget constraints that usually exist and can significantly affect the trade-off decision between the sample size and sequencing depth [41]. Thus, reliance on GEO datasets, which are publicly available, allow for re-analysis of existing data to define potentially new avenues of further investigation. The biopsies that underwent staining for fibrosis were not obtained from the patients in whom the genetic analysis was performed, but rather were performed on a “validation” cohort of sorts. While specific staining for TGFb and other proteins highlighted by genetic analysis would have been ideal, our ability to do so was limited by the COVID-19 pandemic.

5. Conclusions

In summary, our results support that fibrosis identified by transcriptomic and histologic analysis enables us to distinguish CLI from IC. Fibrosis is an important consequence of CLI, which we confirmed histologically, and may be a target for novel therapies in PAD. This study demonstrates the power of RNA-sequencing to determine features that stratify risk and suggests potential therapy. This study further supports a transcriptomic approach to diseases along a spectrum such as PAD.

Author Contributions

Conceptualization, G.C. and U.S.; Fibrosis staining, X.C.; software, G.C. and A.C.; validation, G.C., R.F. and U.S.; writing and review, G.C. and U.S.; Patient recruitment and sample preparation, I.I.P. and G.P.C. All authors have read and agreed to the published version of the manuscript.

Funding

Research reported in this publication was supported by the National Heart, Lung, and Blood Institute of the National Institutes of Health under Award Number R01HL136556 to US and R01AG049868 to IP. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Acknowledgments

We thank Katherine Helfrich for fibrosis image analysis performed at Center of Biologic Institute of University of Pittsburgh. Guangzhi Cong permanent address is at the Heart Centre of General Hospital of Ningxia Medical University, Yinchuan, Ningxia, 750004. CN. He was a visiting scholar at UPMC from October 2019 to November 2020, where the work was completed.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Data Availability Statement

All metadata from the study cohort have been deposited to Dryad Dataset (Cong, Guangzhi et al. (2020). Fibrosis distinguishes critical limb ischemia patients from claudicants in a transcriptomic and histologic analysis, Dryad, Dataset, [42]).

References

  1. Hiatt, W.R.; Fowkes, F.G.; Heizer, G.; Berger, J.S.; Baumgartner, I.; Held, P.; Katona, B.G.; Mahaffey, K.W.; Norgren, L.; Jones, W.S.; et al. Ticagrelor versus Clopidogrel in Symptomatic Peripheral Artery Disease. N. Engl. J. Med. 2017, 376, 32–40. [Google Scholar] [CrossRef] [Green Version]
  2. Benjamin, E.J.; Muntner, P.; Alonso, A.; Bittencourt, M.S.; Callaway, C.W.; Carson, A.P.; Chamberlain, A.M.; Chang, A.R.; Cheng, S.; Das, S.R.; et al. Heart Disease and Stroke Statistics-2019 Update A Report from the American Heart Association. Circulation 2019, 139, E56–E528. [Google Scholar] [CrossRef]
  3. Davies, M.G. Critical limb ischemia: Introduction. Methodist Debakey Cardiovasc. J. 2012, 8, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Agarwal, S.; Sud, K.; Shishehbor, M.H. Nationwide Trends of Hospital Admission and Outcomes among Critical Limb Ischemia Patients: From 2003–2011. J. Am. Coll. Cardiol. 2016, 67, 1901–1913. [Google Scholar] [CrossRef] [PubMed]
  5. Aquino, R.; Johnnides, C.; Makaroun, M.; Whittle, J.C.; Muluk, V.S.; Kelley, M.E.; Muluk, S.C. Natural history of claudication: Long-term serial follow-up study of 1244 claudicants. J. Vasc. Surg. 2001, 34, 962–970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Matzke, S.; Lepantalo, M. Claudication does not always precede critical leg ischemia. Vasc. Med. 2001, 6, 77–80. [Google Scholar] [CrossRef] [PubMed]
  7. Shishehbor, M.H.; White, C.J.; Gray, B.H.; Menard, M.T.; Lookstein, R.; Rosenfield, K.; Jaff, M.R. Critical Limb Ischemia: An Expert Statement. J. Am. Coll. Cardiol. 2016, 68, 2002–2015. [Google Scholar] [CrossRef]
  8. Bunte, M.C.; Jacob, J.; Nudelman, B.; Shishehbor, M.H. Validation of the relationship between ankle-brachial and toe-brachial indices and infragenicular arterial patency in critical limb ischemia. Vasc. Med. 2015, 20, 23–29. [Google Scholar] [CrossRef]
  9. Shishehbor, M.H.; Hammad, T.A. Treatment of Infrapopliteal Disease in Critical Limb Ischemia: Beyond Angioplasty. Circ. Cardiovasc. Interv. 2016, 9, e003882. [Google Scholar] [CrossRef] [Green Version]
  10. Schiano, C.; Costa, V.; Aprile, M.; Grimaldi, V.; Maiello, C.; Esposito, R.; Soricelli, A.; Colantuoni, V.; Donatelli, F.; Ciccodicola, A.; et al. Heart failure: Pilot transcriptomic analysis of cardiac tissue by RNA-sequencing. Cardiol. J. 2017, 24, 539–553. [Google Scholar] [CrossRef] [Green Version]
  11. Matkovich, S.J. Multiomic approaches to delineate the pathogenesis of cardiac disease. Curr. Opin. Cardiol. 2019, 34. [Google Scholar] [CrossRef] [PubMed]
  12. Luo, Y.L.; Huang, L.J.; Luo, W.J.; Ye, S.; Hu, Q.H. Genomic analysis of lncRNA and mRNA profiles in circulating exosomes of patients with rheumatic heart disease. Biol. Open 2019, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lam, C.K.; Tian, L.; Belbachir, N.; Wnorowski, A.; Shrestha, R.; Ma, N.; Kitani, T.; Rhee, J.W.; Wu, J.C. Identifying the Transcriptome Signatures of Calcium Channel Blockers in Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes. Circ. Res. 2019, 125, 212–222. [Google Scholar] [CrossRef]
  14. Ryan, T.E.; Yamaguchi, D.J.; Schmidt, C.A.; Zeczycki, T.N.; Shaikh, S.R.; Brophy, P.; Green, T.D.; Tarpey, M.D.; Karnekar, R.; Goldberg, E.J.; et al. Extensive skeletal muscle cell mitochondriopathy distinguishes critical limb ischemia patients from claudicants. JCI Insight 2018, 3. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, C.H.; Di, Y.P. Analysis of RNA Sequencing Data Using CLC Genomics Workbench. Methods Mol. Biol. 2020, 2102, 61–113. [Google Scholar] [CrossRef]
  16. Conesa, A.; Madrigal, P.; Tarazona, S.; Gomez-Cabrero, D.; Cervera, A.; McPherson, A.; Szczesniak, M.W.; Gaffney, D.J.; Elo, L.L.; Zhang, X.; et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016, 17, 13. [Google Scholar] [CrossRef] [Green Version]
  17. Kramer, A.; Green, J.; Pollard, J., Jr.; Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 2014, 30, 523–530. [Google Scholar] [CrossRef]
  18. Zhao, X.; Kwan, J.Y.Y.; Yip, K.; Liu, P.P.; Liu, F.F. Targeting metabolic dysregulation for fibrosis therapy. Nat. Rev. Drug Discov. 2020, 19, 57–75. [Google Scholar] [CrossRef]
  19. Xu, J.; Cui, X.; Li, J.; Koutakis, P.; Pipinos, I.; Tzeng, E.; Chen, A.; Sachdev, U. Chloroquine improves the response to ischemic muscle injury and increases HMGB1 after arterial ligation. J. Vasc. Surg. 2018, 67, 910–921. [Google Scholar] [CrossRef] [Green Version]
  20. Schipke, J.; Brandenberger, C.; Rajces, A.; Manninger, M.; Alogna, A.; Post, H.; Muhlfeld, C. Assessment of cardiac fibrosis: A morphometric method comparison for collagen quantification. J. Appl. Physiol. 2017, 122, 1019–1030. [Google Scholar] [CrossRef]
  21. Available online: https://imagej.nih.gov/ij/ (accessed on 1 October 2020).
  22. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Available online: http://www.R-project.org (accessed on 1 October 2020).
  24. Mahdy, M.A.A. Skeletal muscle fibrosis: An overview. Cell Tissue Res. 2019, 375, 575–588. [Google Scholar] [CrossRef] [PubMed]
  25. Mizzi, A.; Cassar, K.; Bowen, C.; Formosa, C. The progression rate of peripheral arterial disease in patients with intermittent claudication: A systematic review. J. Foot Ankle Res. 2019, 12, 40. [Google Scholar] [CrossRef] [PubMed]
  26. Stather, P.W.; Sylvius, N.; Wild, J.B.; Choke, E.; Sayers, R.D.; Bown, M.J. Differential microRNA expression profiles in peripheral arterial disease. Circ. Cardiovasc. Genet. 2013, 6, 490–497. [Google Scholar] [CrossRef] [Green Version]
  27. Cooke, J.P.; Wilson, A.M. Biomarkers of peripheral arterial disease. J. Am. Coll. Cardiol. 2010, 55, 2017–2023. [Google Scholar] [CrossRef] [Green Version]
  28. Baczynska, D.; Michalowska, D.; Barc, P.; Skora, J.; Karczewski, M.; Sadakierska-Chudy, A. The expression profile of angiogenic genes in critical limb ischemia popliteal arteries. J. Physiol. Pharm. 2016, 67, 353–362. [Google Scholar]
  29. Hamburg, N.M.; Creager, M.A. Pathophysiology of Intermittent Claudication in Peripheral Artery Disease. Circ. J. 2017, 81, 281–289. [Google Scholar] [CrossRef] [Green Version]
  30. Ha, D.M.; Carpenter, L.C.; Koutakis, P.; Swanson, S.A.; Zhu, Z.; Hanna, M.; DeSpiegelaere, H.K.; Pipinos, I.I.; Casale, G.P. Transforming growth factor-beta 1 produced by vascular smooth muscle cells predicts fibrosis in the gastrocnemius of patients with peripheral artery disease. J. Transl. Med. 2016, 14, 39. [Google Scholar] [CrossRef] [Green Version]
  31. Casanegra, A.I.; Stoner, J.A.; Tafur, A.J.; Pereira, H.A.; Rathbun, S.W.; Gardner, A.W. Differences in galectin-3, a biomarker of fibrosis, between participants with peripheral artery disease and participants with normal ankle-brachial index. Vasc. Med. 2016, 21, 437–444. [Google Scholar] [CrossRef] [Green Version]
  32. Pizzimenti, M.; Riou, M.; Charles, A.L.; Talha, S.; Meyer, A.; Andres, E.; Chakfé, N.; Lejay, A.; Geny, B. The Rise of Mitochondria in Peripheral Arterial Disease Physiopathology: Experimental and Clinical Data. J. Clin. Med. 2019, 8, 2125. [Google Scholar] [CrossRef] [Green Version]
  33. Rodenburg, R.J. Biochemical diagnosis of mitochondrial disorders. J. Inherit. Metab. Dis. 2011, 34, 283–292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Iannetti, E.F.; Prigione, A.; Smeitink, J.A.M.; Koopman, W.J.H.; Beyrath, J.; Renkema, H. Live-Imaging Readouts and Cell Models for Phenotypic Profiling of Mitochondrial Function. Front. Genet. 2019, 10, 131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Conte, M.S.; Bradbury, A.W.; Kolh, P.; White, J.V.; Dick, F.; Fitridge, R.; Mills, J.L.; Ricco, J.B.; Suresh, K.R.; Murad, M.H.; et al. Global vascular guidelines on the management of chronic limb-threatening ischemia. J. Vasc. Surg. 2019, 69, 3S–125S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Herraiz-Adillo, A.; Martinez-Vizcaino, V.; Cavero-Redondo, I.; Alvarez-Bueno, C.; Garrido-Miguel, M.; Notario-Pacheco, B. Diagnostic Accuracy Study of an Oscillometric Ankle-Brachial Index in Peripheral Arterial Disease: The Influence of Oscillometric Errors and Calcified Legs. PLoS ONE 2016, 11, e0167408. [Google Scholar] [CrossRef]
  37. Sukul, D.; Grey, S.F.; Henke, P.K.; Gurm, H.S.; Grossman, P.M. Heterogeneity of Ankle-Brachial Indices in Patients Undergoing Revascularization for Critical Limb Ischemia. JACC Cardiovasc. Interv. 2017, 10, 2307–2316. [Google Scholar] [CrossRef]
  38. Baues, M.; Dasgupta, A.; Ehling, J.; Prakash, J.; Boor, P.; Tacke, F.; Kiessling, F.; Lammers, T. Fibrosis imaging: Current concepts and future directions. Adv. Drug Deliv. Rev. 2017, 121, 9–26. [Google Scholar] [CrossRef]
  39. Cummings, B.B.; Marshall, J.L.; Tukiainen, T.; Lek, M.; Donkervoort, S.; Foley, A.R.; Bolduc, V.; Waddell, L.B.; Sandaradura, S.A.; O’Grady, G.L.; et al. Improving genetic diagnosis in Mendelian disease with transcriptome sequencing. Sci. Transl. Med. 2017, 9. [Google Scholar] [CrossRef] [Green Version]
  40. Das, S.; Frisk, C.; Eriksson, M.J.; Walentinsson, A.; Corbascio, M.; Hage, C.; Kumar, C.; Asp, M.; Lundeberg, J.; Maret, E.; et al. Transcriptomics of cardiac biopsies reveals differences in patients with or without diagnostic parameters for heart failure with preserved ejection fraction. Sci. Rep. 2019, 9, 3179. [Google Scholar] [CrossRef] [Green Version]
  41. Ching, T.; Huang, S.; Garmire, L.X. Power analysis and sample size estimation for RNA-Seq differential expression. RNA 2014, 20, 1684–1696. [Google Scholar] [CrossRef] [Green Version]
  42. Available online: https://doi.org/10.5061/dryad.rr4xgxd67 (accessed on 1 October 2020).
Jcm 09 03974 g001 550
Figure 1. Flow Diagram of Study design and data processing. RNA sequence data of 51 skeleton muscle biopsies in peripheral arterial disease (PAD) (16 of critical limb ischemia (CLI) and 20 of intermittent claudication (IC)) and non-PAD controls (N = 15) were retrieved from the Gene Expression Omnibus (GEO) database (GSE120642). After principal component analysis, Uniform Manifold Approximation and Projection (UMAP), differential expression analysis, hierarchical cluster analysis, and pathway enrichment analysis were completed using Ingenuity Pathway Analysis (http://www.ingenuity.com). The ability of fibrosis modular gene signatures to discriminate CLI from IC was assessed by ANOVA and receiver operating curve (ROC). For validation, we evaluated fibrosis in 36 muscle biopsies from non-PAD controls, IC, and CLI patients using Masson’s trichrome staining. After the measurement of fibrosis area and percentage, we conducted ROC analysis. PCA, principal component analysis.
Figure 1. Flow Diagram of Study design and data processing. RNA sequence data of 51 skeleton muscle biopsies in peripheral arterial disease (PAD) (16 of critical limb ischemia (CLI) and 20 of intermittent claudication (IC)) and non-PAD controls (N = 15) were retrieved from the Gene Expression Omnibus (GEO) database (GSE120642). After principal component analysis, Uniform Manifold Approximation and Projection (UMAP), differential expression analysis, hierarchical cluster analysis, and pathway enrichment analysis were completed using Ingenuity Pathway Analysis (http://www.ingenuity.com). The ability of fibrosis modular gene signatures to discriminate CLI from IC was assessed by ANOVA and receiver operating curve (ROC). For validation, we evaluated fibrosis in 36 muscle biopsies from non-PAD controls, IC, and CLI patients using Masson’s trichrome staining. After the measurement of fibrosis area and percentage, we conducted ROC analysis. PCA, principal component analysis.
Jcm 09 03974 g001
Jcm 09 03974 g002 550
Figure 2. Unique gene expression profile presenting in CLI patients. (A) Principal component analysis (PCA) of the samples. Each dot represents one subject. (B) (UMAP) of the samples. Each dot represents one subject. (C) Heatmap of all differentially expressed genes among groups. Samples are on columns, genes are on columns, and the heat map is based on standardized gene expression values. (D) Volcano plot representing differential gene expression among groups. Each dot on the plot is a single gene. Horizontal axis: fold change; vertical axis: false discovery rate (FDR) p-value (in log10 scale) by Wald test. Color coding is based on the fold change. Thick vertical lines highlight fold changes of −2 and +2, while a thick horizontal line represents a FDR p-value of 0.05. N = 15 for non-PAD controls, N = 20 for IC, N = 16 for CLI. PC1, principal component; CLI, critical limb ischemia; IC, intermittent claudication.
Figure 2. Unique gene expression profile presenting in CLI patients. (A) Principal component analysis (PCA) of the samples. Each dot represents one subject. (B) (UMAP) of the samples. Each dot represents one subject. (C) Heatmap of all differentially expressed genes among groups. Samples are on columns, genes are on columns, and the heat map is based on standardized gene expression values. (D) Volcano plot representing differential gene expression among groups. Each dot on the plot is a single gene. Horizontal axis: fold change; vertical axis: false discovery rate (FDR) p-value (in log10 scale) by Wald test. Color coding is based on the fold change. Thick vertical lines highlight fold changes of −2 and +2, while a thick horizontal line represents a FDR p-value of 0.05. N = 15 for non-PAD controls, N = 20 for IC, N = 16 for CLI. PC1, principal component; CLI, critical limb ischemia; IC, intermittent claudication.
Jcm 09 03974 g002
Jcm 09 03974 g003 550
Figure 3. IPA of genes differentially expressed between CLI and IC and Con. Differentially expressed genes calculated using DeSeq2 (FC > 2, FDR < 0.05) were merged and submitted for IPA analysis. (A,D) A selection of highly significant functions and their activation z-scores are shown in the bar graphs (orange for positive Z score, blue for negative Z score, and gray for no activation pattern available). The x-axis corresponds to the –log of the p-value (Fisher’s exact test); (B,E) a selection of highly significant functions and their activation percentage are shown in the stack bar graphs (red for activation, green for inhibition. Number for each bar for the total molecular in this pathway.) (C,F) Overlaying pathways. (C) CLI vs. Con; (F) CLI vs. IC; (A–C) for CLI vs. Con; (D–F) for CLI vs. IC. (G) Bar and stacked bar graph of first five function and their activation when IC vs. Con. (H) Canonical pathway heatmap of three comparisons. 3.4 Fibrosis pathways are enriched in CLI patients compared to IC and controls. TCA, Tricarboxylic Acid Cycle; FC, Fold Change; ILK, Integrin Linked Kinase; IPA, Ingenuity Pathway Analysis.
Figure 3. IPA of genes differentially expressed between CLI and IC and Con. Differentially expressed genes calculated using DeSeq2 (FC > 2, FDR < 0.05) were merged and submitted for IPA analysis. (A,D) A selection of highly significant functions and their activation z-scores are shown in the bar graphs (orange for positive Z score, blue for negative Z score, and gray for no activation pattern available). The x-axis corresponds to the –log of the p-value (Fisher’s exact test); (B,E) a selection of highly significant functions and their activation percentage are shown in the stack bar graphs (red for activation, green for inhibition. Number for each bar for the total molecular in this pathway.) (C,F) Overlaying pathways. (C) CLI vs. Con; (F) CLI vs. IC; (A–C) for CLI vs. Con; (D–F) for CLI vs. IC. (G) Bar and stacked bar graph of first five function and their activation when IC vs. Con. (H) Canonical pathway heatmap of three comparisons. 3.4 Fibrosis pathways are enriched in CLI patients compared to IC and controls. TCA, Tricarboxylic Acid Cycle; FC, Fold Change; ILK, Integrin Linked Kinase; IPA, Ingenuity Pathway Analysis.
Jcm 09 03974 g003
Jcm 09 03974 g004 550
Figure 4. Network overview of modular expression patterns of fibrosis pathway. (A) CLI vs. Con; (B) CLI vs. IC; (C) IC vs Con. Ingenuity Pathway Analysis of the differentially regulated genes when CLI vs. control and CLI vs IC. The network is displayed graphically as nodes (genes/gene products) and edges (biological relationship between nodes). The node color intensity indicates the fold change expression of genes, with red representing activation, green representing down-regulation of genes, orange representing prediction activation, and blue representing prediction inhibition; the lines indicate the type of interaction.
Figure 4. Network overview of modular expression patterns of fibrosis pathway. (A) CLI vs. Con; (B) CLI vs. IC; (C) IC vs Con. Ingenuity Pathway Analysis of the differentially regulated genes when CLI vs. control and CLI vs IC. The network is displayed graphically as nodes (genes/gene products) and edges (biological relationship between nodes). The node color intensity indicates the fold change expression of genes, with red representing activation, green representing down-regulation of genes, orange representing prediction activation, and blue representing prediction inhibition; the lines indicate the type of interaction.
Jcm 09 03974 g004
Jcm 09 03974 g005 550
Figure 5. Fibrosis pathway core genes distinguish CLI from IC and control. (A) Representative gene expression of fibrosis pathway by transcriptomic by box and whiskers, violins, and points plot; each dot is a sample. The plot title is the gene symbol of the selected gene. Expression levels are log (total counts) on the vertical axis. * p < 0.05, compared with non-PAD controls by wald test, # p < 0.05, compared with IC by wald test. (B) Receiver operating curve (ROC) (spell out ROC) for the diagnosis of CLI of TGF β signaling, collagen deposition signaling, and vascular endothelial growth factor (VEGF) signaling. All area under curve (AUC) > 0.75. N = 15 for non-PAD control, N = 20 for IC, and N = 16 for CLI.
Figure 5. Fibrosis pathway core genes distinguish CLI from IC and control. (A) Representative gene expression of fibrosis pathway by transcriptomic by box and whiskers, violins, and points plot; each dot is a sample. The plot title is the gene symbol of the selected gene. Expression levels are log (total counts) on the vertical axis. * p < 0.05, compared with non-PAD controls by wald test, # p < 0.05, compared with IC by wald test. (B) Receiver operating curve (ROC) (spell out ROC) for the diagnosis of CLI of TGF β signaling, collagen deposition signaling, and vascular endothelial growth factor (VEGF) signaling. All area under curve (AUC) > 0.75. N = 15 for non-PAD control, N = 20 for IC, and N = 16 for CLI.
Jcm 09 03974 g005
Jcm 09 03974 g006 550
Figure 6. Validation of fibrosis pathway to distinguish CLI from IC and control. Masson trichrome staining of gastrocnemius muscle biopsies. (A) Representative fibrosis staining image (×20 magnification). (B) Assessment of fibrosis by area and percentage. Each dot is a sample. Comparison between HA, IC, and CLI group using single one-way analysis of variance (ANOVA) followed by Bonferroni t-test; ** p < 0.01, *** p < 0.001, **** p < 0.0001. (C) ROC for the diagnosis of CLI. Fibrosis area (AUC area = 0.81) and fibrosis percentage (AUC area = 0.87). N = 12 for non-PAD control, IC, and CLI. NS, not significant.
Figure 6. Validation of fibrosis pathway to distinguish CLI from IC and control. Masson trichrome staining of gastrocnemius muscle biopsies. (A) Representative fibrosis staining image (×20 magnification). (B) Assessment of fibrosis by area and percentage. Each dot is a sample. Comparison between HA, IC, and CLI group using single one-way analysis of variance (ANOVA) followed by Bonferroni t-test; ** p < 0.01, *** p < 0.001, **** p < 0.0001. (C) ROC for the diagnosis of CLI. Fibrosis area (AUC area = 0.81) and fibrosis percentage (AUC area = 0.87). N = 12 for non-PAD control, IC, and CLI. NS, not significant.
Jcm 09 03974 g006
Jcm 09 03974 g007 550
Figure 7. Derivation cohort included 15 non-PAD controls, 20 IC, and 16 CLI patients. Validation cohort included 12 non-PAD controls, 12 IC, and 12 CLI patients. Both fibrosis pathway identified by transcriptomics and fibrosis quantification by Masson trichrome of skeletal muscle demonstrated an area under the curve (AUC) over 0.75. This indicates that fibrosis distinguishes CLI patients from IC in a transcriptomic and histologic analysis.
Figure 7. Derivation cohort included 15 non-PAD controls, 20 IC, and 16 CLI patients. Validation cohort included 12 non-PAD controls, 12 IC, and 12 CLI patients. Both fibrosis pathway identified by transcriptomics and fibrosis quantification by Masson trichrome of skeletal muscle demonstrated an area under the curve (AUC) over 0.75. This indicates that fibrosis distinguishes CLI patients from IC in a transcriptomic and histologic analysis.
Jcm 09 03974 g007
Table
Table 1. Patient characteristics.
Table 1. Patient characteristics.
Non-PAD ControlsIntermittent ClaudicantCritical Limb Ischemiap Value
Derivation cohortn = 32n = 27n = 19
Age (Mean ± SD)61(7.3)61(7.5)64(10)0.18
ABI (Mean ± SD) 0.65(0.21)0.35(0.30)0.0007
Validation cohortn = 12n = 12n = 12
Age (Mean ± SD)66.75(4.93)68.55(8.68)63.58(5.05)0.01
ABI (Mean ± SD)1.10(0.10)0.62(0.13)0.31(0.24)0.0001
PAD, peripheral artery disease; SD, standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cong, G.; Cui, X.; Ferrari, R.; Pipinos, I.I.; Casale, G.P.; Chattopadhyay, A.; Sachdev, U. Fibrosis Distinguishes Critical Limb Ischemia Patients from Claudicants in a Transcriptomic and Histologic Analysis. J. Clin. Med. 2020, 9, 3974. https://doi.org/10.3390/jcm9123974

AMA Style

Cong G, Cui X, Ferrari R, Pipinos II, Casale GP, Chattopadhyay A, Sachdev U. Fibrosis Distinguishes Critical Limb Ischemia Patients from Claudicants in a Transcriptomic and Histologic Analysis. Journal of Clinical Medicine. 2020; 9(12):3974. https://doi.org/10.3390/jcm9123974

Chicago/Turabian Style

Cong, Guangzhi, Xiangdong Cui, Ricardo Ferrari, Iraklis I. Pipinos, George P. Casale, Ansuman Chattopadhyay, and Ulka Sachdev. 2020. "Fibrosis Distinguishes Critical Limb Ischemia Patients from Claudicants in a Transcriptomic and Histologic Analysis" Journal of Clinical Medicine 9, no. 12: 3974. https://doi.org/10.3390/jcm9123974

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop