Biochemical and genomic identification of novel biomarkers in progressive sarcoidosis: HBEGF, eNAMPT, and ANG-2
Nancy G. Casanova1 
Vivian Reyes-Hernon1 Taylor Gregory1 Belinda Sun2 Tadeo Bermudez1 Matthew K. Hufford1 Radu C. Oita1 Sara M. Camp1 
Gabriela Hernandez-Molina3 Jorge Rojas Serrano4 Xiaoguang Sun1 Jocelyn Fimbres3 Mehdi Mirsaeidi5 Saad Sammani1 Christian Bime1 
Joe G. N. Garcia1*- 1Department of Medicine, University of Arizona Health Sciences, Tucson, AZ, United States
- 2Department of Pathology, University of Arizona Health Sciences, Tucson, AZ, United States
- 3Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, México City, Mexico
- 4Instituto Nacional de Enfermedades Respiratorias, México City, Mexico
- 5Department of Medicine, College of Medicine, University of Florida, Jacksonville, FL, United States
Background: Progressive pulmonary fibrosis is a serious complication in subjects with sarcoidosis. The absence of reliable, non-invasive biomarkers that detect early progression exacerbates the difficulty in predicting sarcoidosis severity. To potentially address this unmet need, we evaluated a panel of markers for an association with sarcoidosis progression (HBEGF, NAMPT, IL1-RA, IL-6, IL-8, ANG-2). This panel encompasses proteins related to inflammation, vascular injury, cell proliferation, and fibroblast mitogenesis processes.
Methods: Plasma biomarker levels and biomarker protein expression in lung and lymph nodes tissues (immunohistochemical studies) from sarcoidosis subjects with limited disease and progressive (complicated) sarcoidosis were performed. Gene expression of the protein-coding genes included in this panel was analyzed using RNAseq in sarcoidosis granulomatous tissues from lung and lymph nodes.
Results: Except for IL-8, plasma levels of each biomarker—eNAMPT, IL-1RA, IL-6, ANG-2, and HBEGF—were significantly elevated in sarcoidosis subjects compared to controls. In addition, plasma levels of HBEGF were elevated in complicated sarcoidosis, while eNAMPT and ANG-2 were observed to serve as markers of lung fibrosis in a subgroup of complicated sarcoidosis. Genomic studies corroborated HBEGF and NAMPT among the top dysregulated genes and identified cytokine-related and fibrotic pathways in lung granulomatous tissues from sarcoidosis.
Conclusion: These findings suggest HBEGF, eNAMPT, and ANG-2 may serve as potential novel indicators of the clinical severity of sarcoidosis disease.
Introduction
Sarcoidosis is a multisystemic disease with 20% of afflicted patients estimated to develop progressive lung fibrosis (1). Despite attempts at utilizing genomic (2, 3) and genetic modalities (4–6), the development of progressive pulmonary fibrosis, i.e., complicated sarcoidosis, currently evades prediction due to the absence of reliable, non-invasive biomarkers of early progression. In prior sarcoidosis gene expression profiling studies in peripheral blood mononuclear cells (PBMCs) designed to address this important need, we showed dysregulated candidate genes associated with T cell receptor, cytokine-cytokine, and Jak-STAT signaling in complicated progressive phenotypes (2, 3). Similar studies in pulmonary fibrosis, analyzed peripheral blood gene expression, identifying transcripts associated with disease severity (7); including the activation of STAT3 by IL-6 in progressive fibrosis (8). Interrogation of specific up-regulated genes that link to sarcoidosis severity identified the HBEGF gene as differentially expressed in sarcoidosis compared to healthy controls and accurately discriminated subjects with limited sarcoidosis from subjects with progressive lung sarcoidosis (2), HBEGF encodes HBEGF (heparin-binding epidermal growth factor), a member of the EGF growth factor family expressed in the bronchial epithelium, smooth muscle, CD4 + T cells, and cardiac muscle. HBEGF signaling is a critical contributor to repair and regeneration processes that restore tissue homeostasis following injury, promoting wound healing, angiogenesis, and adipogenesis (9). The finding of increased HBEGF expression linked to disease severity suggested a potential utility of HBEGF as a plasma-based sarcoidosis biomarker.
Identical genomic approaches in subjects with ARDS and pulmonary hypertension (10–13) identified extracellular nicotinamide phosphoribosyl transferase (eNAMPT) as a damage-associated molecular pattern protein (DAMP) and master regulator of innate immunity pathways via ligation of the Toll-Like receptor 4 (TLR4) and subsequent NF-κB activation (14), processes directly implicated involvement in organ fibrosis (15, 16). Angiopoietin2 (ANG-2) is an endothelial-derived protein that increases junction instability (17) and has been implicated in directly controlling multiple inflammatory-related signals. Elevated ANG-2, along with circulating levels of eNAMPT, Interleukin-6 (IL-6), and Interleukin-1 receptor agonist (IL-1RA) were predictors of 28-day mortality in ARDS subjects (10, 18). In sarcoidosis, a IL-1a SNP was significantly overrepresented in sarcoidosis subjects (compared to controls) (19). IL-6 has been postulated as an important proinflammatory cytokine involved in Th1-mediated immunopathobiology of sarcoidosis (20).
Based upon these genomic and epigenomic studies identifying potentially novel markers in sarcoidosis and other inflammatory lung disorders (19–25), the current study selected HBEGF, NAMPT, IL-1RA, IL-6, IL-8, Ang-2 as a panel of plasma biomarkers postulated as potentially useful for identifying subjects with progressive sarcoidosis and lung fibrosis. In addition to measurements of this panel of plasma biomarkers in 127 sarcoidosis subjects and 82 healthy controls, we performed immunohistochemistry (IHC) studies and RNA sequencing to validate the potential contribution for these biomarkers to serve as diagnostic and prognostic biomarkers in subjects with sarcoidosis. These complementary biochemical and genomic approaches indicate that HBEGF, eNAMPT, and ANG-2 may serve as potentially novel indicators of the clinical severity of sarcoidosis.
Materials and methods
Cohort samples and demographics
A total of 127 plasma samples from sarcoidosis and 82 healthy controls were analyzed. The clinical characteristics and demographic information are presented in Table 1. Subjects with sarcoidosis were stratified for complicated phenotypes (n = 62). Complicated sarcoidosis was assigned to those exhibiting lung involvement with parenchymal lung disease by CT scan or radiographic Scadding stages III and IV; Forced vital capacity (FVC) < 50%; and cardiac or neurological involvement in addition to lung involvement, as previously described (26). We also identified within the complicated group a subgroup of sarcoidosis with pulmonary fibrosis (n = 19) and analyzed this subphenotype independently. For those with missing data for chest involvement, we used the ICD-9 and ICD10 codes to define the complicated status. Written informed consent was obtained from all participating subjects with protocols approved by the Institutional Review Board (IRB) at each participating institution. University of Arizona (IRB #1509097312R001), University of Miami (IRB# 20150612). De-identified samples were also obtained from a sarcoidosis cohort from the Instituto Nacional de Ciencias Medicas y Nutricion in Mexico City (Ref.#1711).
Table 1. Patient characteristics.
Plasma biomarker detection
Plasma was obtained from venous blood collected in EDTA tubes. Following centrifugation (2,500 × g for 15 min), samples were stored at ≤ –70°C until biomarker measurement was conducted. 150 μl of plasma was used to quantify HBEGF levels using Quantikine® ELISA Human HBEGF immunoassay (R&D Systems®), quantitative sandwich enzyme immunoassay technique. Plasma levels of eNAMPT, IL-6, IL-8, Ang-2, and IL-1RA, were measured in 50 μl of plasma utilizing an electrochemiluminescent multiplex immunoassay predesigned panel from MesoScale (Meso Scale Discovery, MSD®) (27–29) (see details in Supplementary materials and methods). Supplementary Table 1 shows the detection limits for each biomarker on every assay.
Immunohistochemistry analyses
Formalin-fixed, paraffin-embedded lung tissues were obtained from biopsy or autopsy specimens of individuals with sarcoidosis (n = 11) and healthy controls (n = 5) in compliance with their respective IRB. Lung (n = 8) and lymph node (n = 8) tissues fixed in 10% formalin underwent H&E and IHC staining for three different biomarkers (HBEGF, ANG-2, and NAMPT) using anti-Human HBEGF antibody (R&D Systems, AF-259-NA), ANG-2 recombinant rabbit mAb (ThermoFisher, JM71-34), and rabbit anti-human NAMPT pAb (Bethyl Laboratories, Montgomery, TX, USA) as previously described (27). Histopathology images, graded by a pathologist blinded to study groups, were selected for quantification analysis of H&E, HBEGF, ANG-2, and NAMPT staining area by Image J software version 1.53 h (30).
Gene expression
Previously generated RNA sequencing datasets for sarcoidosis granulomatous lung and lymph node tissues and healthy tissue were utilized to assess the transcriptomic profile of targeted genes, GO accession number: GSE157671. Gene expression was assayed using an Illumina HiSeq 2000 HTG EdgeSeq™ Oncology-biomarker panel. The differential expression analysis was performed using the Limma package (31). Transcripts with a fold change (FC) > 2 and a q-value < 0.01 were deemed as differentially expressed. Detailed methods and analysis have been previously reported (32).
Statistical analysis
Numerical variables were analyzed using non-parametric Mann–Whitney test and Kruskal–Wallis test when parametric assumptions were not met. Group comparisons were performed using Spearman correlation. Sensitivity and specificity of the variables were calculated to evaluate the area under the curve (AUC) of the Receiver Operator characteristics (ROC) in each of the plasma markers. Cut-off values were selected using J statistics (Youden index). A significance of 0.05% (p < 0.05) was required to consider statistical significance. All analyses were performed with Stata v.17 (StataCorp, TX), and Graphpad Prism v. 8.0 software (San Diego, CA).
Results
Patient characteristics
The patient characteristics in the sarcoidosis (n = 127) and healthy control groups (n = 82) are presented in Table 1. There was no significant difference in age between groups (p = 0.37). The sarcoidosis group exhibited a higher proportion of females and African Americans compared to controls (p < 0.01). Radiographic or chest CT scan data indicating lung involvement stage was available on 71 of the 127 sarcoidosis subjects, with pulmonary fibrosis reported in 19 (15%) of total sarcoidosis subjects. 62 subjects (49%) of the sarcoidosis cohort were classified as complicated sarcoidosis; documented parenchymal lung disease by CT scan or radiographic Scadding stages III and IV; FVC < 50%; and/or cardiac or neurological involvement. The characteristics of the complicated and uncomplicated groups are presented in Table 2.
Table 2. Characteristics of the phenotypic groups.
Circulating biomarkers in sarcoidosis
Plasma levels of eNAMPT, IL-6, ANG-2, IL-1RA, and HBEGF, but not IL-8 levels, were significantly elevated in sarcoidosis subjects compared to healthy controls (Table 3 and Figure 1). The diagnostic accuracy of the plasma biomarkers was assessed via estimates of sensitivity and specificity in addition to 95% confidence intervals and AUCs. IL-1RA exhibited the highest discriminatory accuracy between sarcoidosis and healthy controls (AUC 0.93, 95% CI of 0.90–0.95), followed by ANG-2 (AUC 0.85, 95% CI of 0.79–0.90). eNAMPT, IL-6, and HBEGF exhibited similar AUC values (0.74–0.78) (Figure 2). Spearman correlations to assess inter-marker relationships showed a strong positive correlation between IL-1RA and eNAMPT plasma levels (r = 0.42, p < 0.00001) and a moderate correlation between IL-1RA and IL-6 levels (r = 0.32, p < 0.0002). A predictive model with regression analysis was applied to examine the correlation between eNAMPT and IL-1RA and showed IL-1RA values to significantly predict eNAMPT levels (F2,124 = 14.15, prob F < 0.001). This significance was not observed for IL-6 levels after adjusting for regression standard error (Supplementary Figure 1).
Table 3. Biomarker plasma levels in sarcoidosis.
Figure 1. Biomarker plasma levels in sarcoidosis and healthy controls. Plasma levels medians are plotted. Compared to healthy controls, each biomarker exhibited higher values in plasma of sarcoidosis patients. Except for IL-8 all of plasma biomarkers were significantly higher in sarcoidosis. eNAMPT (2.1 vs. 1.2 ng/ml), IL-6 (1.2 vs. 0.6 pg/ml), ANG-2 (2.7 vs. 1.0 pg/ml), HBEGF (5.9 vs. 1 pg/ml), and IL-1 RA (227 vs. 82.3 pg/ml). Mann–Whitney test was calculated for each comparison; ****p-values < 0.0001 and *p < 0.05.
Figure 2. Circulating biomarkers in complicated sarcoidosis and pulmonary fibrosis. (A) The plasma levels of HBEGF were significantly elevated (**p < 0.003) in a group of 62 complicated sarcoidosis (parenchymal lung disease by CT scan, radiographic Scadding stages III and IV; Forced vital capacity (FVC) < 50%; cardiac or neurological involvement in addition to lung involvement) compared to non-complicated sarcoidosis. (B) ANG-2 plasma levels were significantly (p < 0.017) elevated in 19 subjects with sarcoidosis and pulmonary fibrosis (confirmed fibrotic changes on high-resolution chest CT scan or radiographic stage IV) compared to those with non-fibrosis.
Circulating biomarkers in complicated sarcoidosis including pulmonary fibrosis
We next assessed the capacity for putative sarcoidosis biomarkers to accurately reflect sarcoidosis severity and compared biomarkers levels in 62 subjects with complicated sarcoidosis. In this group, we included subjects with the following; parenchymal lung disease by CT scan, radiographic Scadding stages III and IV; FVC < 50%; cardiac or neurological involvement in addition to lung involvement. This group was compared to biomarkers values in the non-complicated sarcoidosis group (63 subjects). These analyses identified HBEGF as the sole biomarker to significantly differentiate these two groups, with higher mean HBEGF levels in the complicated group, 47.8 pg/ml ± 109.6 vs. 21.4 pg/ml ± 50.5 (AUC 0.70, 95% CI 0.56–0.86, p = 0.007) (Figure 2). We next compared biomarker levels in a subgroup of 19 sarcoidosis subjects with pulmonary fibrosis, with confirmed fibrotic changes on high-resolution chest CT scan or radiographic stage IV (Table 1) compared to 17 sarcoidosis subjects without parenchymal lung involvement (Scadding stages II and I) (Figure 2). ANG-2 was the sole marker significantly distinguishing the two groups, (3.8 ± 2.9 ng/ml in fibrosis subjects vs. 2.2 ± 0.8 ng/ml in sarcoidosis subjects without lung involvement, p = 0.017), with significant discriminatory power (AUC 0.73, 95% CI 0.57–0.9, p = 0.02). We next assessed the ANG-2 cut-off values predictive of pulmonary fibrosis and identified a range of 2.4 ng/ml and 2.8 ng/ml of ANG-2 in plasma to produce a sensitivity of 72 and 66% and a specificity of 65 and 75%, respectively (Figure 3). ANG-2 levels were not significantly different in between sarcoid subjects with pulmonary involvement (Stage III and IV) and subjects without pulmonary involvement (Stage I and Stage II). Summary statistics of all the biomarker levels in controls and sarcoidosis stratified by comparisons groups according to clinical phenotypes are presented in Table 3.
Figure 3. The receiver-operator characteristics curves (ROC) assessing the diagnostic accuracy power. (A–E) Sarcoidosis vs. controls. eNAMPT (A), IL-6 (B), IL-1RA (C), ANG-2 (D), and HBEGF (E). IL-1RA showed the highest discriminatory accuracy between sarcoidosis and healthy controls (AUC of 0.93, 95% CI of 0.90–0.95), followed by ANG-2 (AUC 0.85, 95% CI of 0.79–0.90). (F,G) Progressive sarcoidosis. ANG-2 diagnostic accuracy performed well to distinguish sarcoidosis pulmonary fibrosis vs. sarcoidosis without pulmonary involvement (Stages I and II) (F) and HBEGF showed a good discriminatory power to differentiate complicated sarcoidosis.
Sarcoidosis lung and lymph node immunohistochemistry
Based upon performance in the immunoassays and association with progression, HBEGF and ANG-2 were next selected for immunohistochemistry staining with Image J quantitation in tissue sections from sarcoidosis lung and lymph node tissue specimens from different donors (Figure 4). We also assessed NAMPT tissue expression due to its strong positive correlation with IL-1RA levels and known linkage to tissue fibrosis (15, 16). HBEGF expression in sarcoidosis lymph nodes was nine times greater than normal lymph node tissues (p = 0.000038) and a three-fold increase in expression in sarcoidosis lung compared to normal lung tissues (p = 0.00004). Similarly, ANG-2 expression in sarcoidosis lymph node tissues (p = 0.0001) and lung tissues (p = 0.0001) was significantly higher compared to healthy tissue. Finally, NAMPT expression in sarcoidosis lymph node (p = 0.0009) and lung tissues (p = 0.005) was significantly higher compared to healthy tissue (Figure 4). These results confirm the elevated expression of the three HBEGF, ANG-2, and eNAMPT plasma biomarkers in sarcoidosis granulomatous tissues.
Figure 4. (A–F) Representative immunohistochemical staining of lymph nodes (A,C,E) and lung tissues (B,D,F) from sarcoidosis compared healthy controls 40 × magnification. ANG-2 (A,B,G), HBEGF (C,D,H), NAMPT (E,F,I). (G–I) Bar graphs (G) illustrate quantification of mean proportions of ANG-2, HBEGF, and NAMPT *p < 0.05. Sarcoidosis in dark gray and normal tissue in light gray. ANG-2 expression in lung tissue and lymph node in sarcoidosis was 32.2 and 33.7%, respectively, in contrast to the expression observed in normal lung and lymph nodes was 9.8 and 11.3%, respectively (p = 1.2 e-04). HBEGF expression in sarcoidosis lymph nodes was 44.8% vs. 4.8% in normal lymph nodes (p = 3.8 e-05); in normal lungs was 18.3% contrasting with 60.4% in lung of sarcoidosis (p = 4.0 e-05). NAMPT mean proportions in lymph nodes were 12.7% in normal tissue and 52% in sarcoidosis (p = 4.9 e-04); while in normal lungs was 23.3% vs. 53% in sarcoidosis (p = 9.2 e-040. in lymph nodes).
Gene expression in sarcoidosis lung tissue
We next interrogated the expression of genes encoding the six plasma biomarkers (HBEGF, IL-1RA, NAMPT, IL6, ANG2, and IL8) using our previously reported publicly available gene expression dataset (GSE157671) generated by NextGen sequencing of micro-dissected sarcoidosis granulomas within lung tissues (6) and mediastinal lymph nodes (11, 32). Comparisons of biomarker panel gene expression failed to identify significant differential expression of the targeted genes between sarcoidosis lymph node granulomas and healthy lymph node tissues. In contrast, HBEGF and NAMPT were among the 73 DEGs in sarcoidosis lung granulomas (Supplementary Table 2) with significant dysregulation in sarcoidosis lung granulomas compared to lung tissue from controls (FRD 10%, FC1.5). Enrichment analysis of the DEGs and protein–protein interaction analysis using String database (33, 34) demonstrated an association of HBEGF and NAMPT together with proteins involved in inflammatory and lung fibrosis pathways (Figure 5). These included MMP9, a matrix metalloproteinase strongly associated with lung injury and fibrosis (35), and NOTCH4, an important regulator of inflammation and lung remodeling previously implicated in sarcoidosis (4, 32, 36) (Figure 5). A total of 33 KEGG pathways were identified with the top significant dysregulated pathways (q-value < 0.03) comprised of cytokine-cytokine receptor interaction, MAPK signaling, IL-17, Jak-STAT signaling, chemokine signaling, NF-kappa B, and TNF signaling pathways (Figure 5). A complete list of all the pathways is presented in Supplementary Table 3. These results are highly consistent with recent studies strongly implicating NAMPT-MMP9 protein interactions in influencing the development and severity of murine radiation-induced lung fibrosis (15, 37).
Figure 5. Enrichment analysis of DEGs in lung tissues from subjects with sarcoidosis compared to healthy lung. Functional protein network interaction showing eNAMPT and HBEGF relationship with regulators of lung fibrosis (MMP9), inflammation, and lung remodeling and NOTCH4, a T-cell activity and angiogenesis regulator. Each node represents a coded protein. The edges represent protein–protein associations (line color indicates type of interaction: green—gene neighborhood, red—gene fusion, blue—co-occurrence, and black—co-expression). Pathway analysis and first shell interactors members included are indicated by the color code. A complete list of dysregulated pathways is presented in Supplementary Table 2.
Discussion
This study presents a multi-dimensional exploration to validate potential inflammatory and lung remodeling biomarkers with tri-level assessment of circulating plasma levels, tissue protein expression, and lung mRNA transcript expression. With the exception of IL-8, each biomarker in the sarcoidosis panel (eNAMPT, IL-6, ANG-2 IL-1RA, and HBEGF) accurately discriminated between healthy controls and sarcoidosis subjects, with an optimal threshold > 74%. Our plasma eNAMPT results strongly contrast with a previous report indicating no difference between sarcoidosis and healthy controls (38), possibly reflecting differences in studied cohorts as well as biomarker analytic methodologies (commercial ELISA vs. MSD). Each panel biomarker exhibited unique biomarker properties with IL-1RA and ANG-2 exhibiting the most diagnostic accuracy in distinguishing sarcoidosis from healthy subjects. We also identified a strong positive correlation between plasma levels of eNAMPT, IL-1RA, and IL-6 levels, a novel finding in sarcoidosis, previously reported in juvenile idiopathic arthritis (39). The synergistic effect of IL-1 and IL-6 was previously demonstrated to be a necessary factor for T-cell activation in sarcoidosis (22).
ANG-2 and HBEGF proved to be useful biomarkers in assessing disease severity with ANG-2 differentiating sarcoidosis pulmonary fibrosis (stage IV) from sarcoidosis without pulmonary involvement (stages I/II). We noted that when Stage III was incorporated in the analysis, ANG-2 levels were no longer discriminatory, an observation potentially related to the insensitivity of CXRs for assessment of early fibrotic changes of the parenchyma compared to HRCT or PET scans (40). HBEGF plasma levels distinguished a more heterogenous group of complicated sarcoidosis based on the presence of progressive parenchymal lung disease (defined by CT scan, radiographic stages III and IV and/or PFT with FVC < 50%) or vital cardiac or neurological involvement. Peripheral blood concentrations of ANG-2 have been reported to be significantly associated to fibrosis in liver (41) and in idiopathic interstitial pneumonias (23). Furthermore, ANG-2 and HBEGF crosstalk signaling pathways are recognized to upregulate fibronectin synthesis and release of metalloproteinases and collagenases in mesangial cells glomerulosclerosis (16). Our results support this mechanism to be relevant in sarcoidosis pathogenesis as well.
Quantification of biomarker tissue expression by IHC was largely concordant with plasma biomarker results, with elevated expression of NAMPT, HBEGF, and ANG-2 in lung and lymph node tissues from subjects with sarcoidosis. The role of these proteins in the development of sarcoidosis pathobiology is poorly understood. HBEGF signaling responses participate in tissue injury/repair and regeneration, wound healing, and angiogenesis (42). HBEGF serum levels were positively correlated with pulmonary ground-glass score indicative of severe fibrosis (43). Our enrichment analysis of DEGs in lung tissues from subjects with sarcoidosis also revealed the network of interaction of HBEGF and NAMPT with MMP-9, a well-known matrix metalloproteinase target in IPF that is strongly associated with lung injury and fibrosis (35) and associated in progressive sarcoidosis with pulmonary infiltrates or fibrosis (44). MMP9 expression was increased more than 30-fold in sarcoidosis tissues (Supplementary Table 2). These results are highly consistent with our recent studies strongly implicating NAMPT-MMP9 protein interactions in influencing the development and severity of murine radiation-induced lung fibrosis (37). HBEGF interacts with SOCS3, a suppressor of cytokine signaling with an ability to inhibit Jak/STAT signaling, with elevated expression in lymph nodes sarcoidosis (45), we previously identified dysregulated HBEGF gene expression in PBMCs sarcoidosis (2). Whether HBEGF is a protective repair response to tissue injury or a direct participant in developing granulomas in sarcoidosis is still unknown, the current study confirms increased HBEGF expression in blood and in lung tissues of sarcoidosis.
Our results also identified other fibrogenesis-associated genes as DEGs in sarcoidosis lungs including AXL, involved in migration, aggregation, and anti-inflammation through inhibition of Toll-like receptors (TLRs). AXL targeting reduces fibrosis development in experimental renal fibrosis and in human intestinal organoid models (46, 47). Similarly, epidermal growth factor receptor (EGRF) in animal models to study the development of SARS-CoV-induced fibrosis, evidence that pulmonary fibrosis is caused by a hyperactive host response to lung injury mediated by EGFR signaling (48).
The protein interactive network identified an important correlation between HBEGF and NOTCH4, a key regulator of T cell activity and the branching angiogenesis associated with sarcoidosis (4, 32, 36). A GWAS study (4) identified NOTCH4 SNPs associated with sarcoidosis severity in African Americans. The most dysregulated DEG identified in sarcoidosis tissues was NPR1, whose expression was markedly suppressed (> 7 fold reduced expression). Chemokine (C-X-C motif) ligand 2 (CXCL2) also called macrophage inflammatory protein (Figure 5), CXCL2 affects neutrophil recruitment and activation (49) and promotes airway smooth muscle cell migration in asthma induced by IL-17 (50). Expression of chemokine (C–C motif) ligand 19 (CCL19), was markedly increased in sarcoidosis tissues, CCL19 is a critical regulator of the induction of T cell activation. Another highly downregulated gene was NR4A1, a key general regulator in the induction of T cell dysfunction (51). Loss of NR4A1 exacerbates organ fibrosis by dysregulating TGF-β pathway (52, 53) a mechanism potentially related to sarcoidosis lung progressive disease.
The identification of eNAMPT as a sarcoidosis biomarker and candidate gene is highly novel. Secreted eNAMPT, a DAMP and ligand for TLR4 (14), is a therapeutic target (54–56). eNAMPT-triggered inflammation dysregulates signaling pathway that disrupts fibrosis/remodeling resolution. The identification of eNAMPT, as a target and biomarker in sarcoidosis is consistent with the top significant dysregulated pathways identified comprised of cytokine-cytokine receptor interaction, MAPK signaling, IL-17, Jak-STAT signaling, chemokine signaling, NF-kappa B, and TNF signaling pathways (Figure 5). These results also support eNAMPT role in sarcoidosis pathogenesis potentially related to stimulation of CD14 + monocytes, promoting leukocytes and fibroblasts motility through the activation of transcription factors, IL-1RA, TNFα, and IL-6 (22, 25). Together, these data reflect the utility of eNAMPT, ANG-2, and HBEGF as diagnostic blood/tissue markers of sarcoidosis severity and as novel therapeutic targets.
There are limitations to our study including the underpowered nature of the studies performed. Our inability to discern significant IHC differences in marker expression between lung and lymph nodes and differences based on clinical phenotype and the effect of treatment in our results are important limitations. We acknowledge that the biomarkers studied here cannot yet be considered as prognostic biomarkers in newly diagnosed sarcoidosis as the study is not longitudinal and the disease duration was largely unknown in our cohorts. Further investigations are needed to corroborate the gene expression of the targeted markers in an increased number of samples and correlate the plasma/tissue expression and the transcriptome profile with fibrosis progression in sarcoidosis and in other interstitial lung diseases and longitudinal assessment of stability of phenotyping over time.
In summary, the current study has demonstrated the diagnostic accuracy of plasma levels of eNAMPT, IL-6, ANG-2, IL-1RA, and HBEGF as a potential panel of biomarkers for sarcoidosis diagnosis and HBEGF, Enampt, and ANG-2 as markers of sarcoidosis progression and lung fibrosis in subjects with complicated sarcoidosis. Genomic studies corroborated plasma protein dysregulation and identified cytokine-related KEGG pathways in lung granulomatous tissues from sarcoidosis subjects and involvement in inflammatory and lung fibrosis pathways. Together these findings have important implications in the implementation of clinical markers to assess sarcoidosis disease severity and suggest HBEGF, eNAMPT, and ANG-2 may serve as potentially novel indicators of the clinical severity of sarcoidosis disease.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE157671.
Ethics statement
The studies involving human participants were reviewed and approved by the University of Arizona (IRB #1509097312R001) and University of Miami (IRB# 20150612). De-identified samples were also obtained from a sarcoidosis cohort from the Instituto Nacional de Ciencias Medicas y Nutrición Salvador Zubirán in Mexico City (Ref. #1711). The patients/participants provided their written informed consent to participate in this study.
Author contributions
NC, CB, JG, MM, and BS: conception and design of the work, the analysis and interpretation of data for the work, the drafting and revision of the manuscript, and approval of the final version to be published. VR-H, TB, TG, GH-M, JS, BS, SS, JF, and MM: collection and analysis of data, revision of the manuscript, and approval of the final version to be published. VR-H, TB, MH, TG, SS, BS, JF, RO, and SC: data collection and assistance with processing and manuscript production and revision. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the R42-HL145930, R42 HL152888, and R01HL141387.
Acknowledgments
We thank all patients for participating in sarcoidosis research and for contributing samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.1012827/full#supplementary-material
Supplementary Figure 1 | Predictive model with regression analysis was used to assess the relationship between our plasma markers. Scatterplot shows the only significant correlation present, it was between plasma levels of eNAMPT and IL-1RA in sarcoidosis. Y-axis shows IL-1RA (pg/mL) and X-axis NAMPT (fitted values). Prediction regression model indicates that IL-1RA values significantly predict plasma eNAMPT levels, with an R2 0.19 (F2,124 = 14.15, prob F < 0.00).
Supplementary Table 1 | The upper limit of detection (ULOD) and lower limit of detection (LLOD) and background concentration of each plasma biomarker.
Supplementary Table 2 | The 73 differentially expressed genes in sarcoidosis lung granulomatous tissues compared to healthy lung tissue. Listed by p-value.
Supplementary Table 3 | KEGG Pathways identified from Differentially Expressed Genes (DEGs) in lung granulomatous tissues from sarcoidosis subjects.
References
1. Judson MA. The clinical features of sarcoidosis: a comprehensive review. Clin Rev Allergy Immunol. (2015) 49:63–78. doi: 10.1007/s12016-014-8450-y
2. Zhou T, Zhang W, Sweiss NJ, Chen ES, Moller DR, Knox KS, et al. Peripheral blood gene expression as a novel genomic biomarker in complicated sarcoidosis. PLoS One. (2012) 7:e44818. doi: 10.1371/journal.pone.0044818
3. Su R, Li MM, Bhakta NR, Solberg OD, Darnell EP, Ramstein J, et al. Longitudinal analysis of sarcoidosis blood transcriptomic signatures and disease outcomes. Eur Respir J. (2014) 44:985–93. doi: 10.1183/09031936.00039714
4. Adrianto I, Lin CP, Hale JJ, Levin AM, Datta I, Parker R, et al. Genome-wide association study of African and European Americans implicates multiple shared and ethnic specific loci in sarcoidosis susceptibility. PLoS One. (2012) 7:e43907. doi: 10.1371/journal.pone.0043907
5. Moller DR, Koth LL, Maier LA, Morris A, Drake W, Rossman M, et al. Rationale and design of the genomic research in alpha-1 antitrypsin deficiency and sarcoidosis (GRADS) study. Sarcoidosis Protocol. Ann Am Thorac Soc. (2015) 12:1561–71. doi: 10.1513/AnnalsATS.201503-172OT
6. Rossman MD, Kreider ME. Lesson learned from ACCESS (A case controlled etiologic study of sarcoidosis). Proc Am Thorac Soc. (2007) 4:453–6. doi: 10.1513/pats.200607-138MS
7. Yang IV, Luna LG, Cotter J, Talbert J, Leach SM, Kidd R, et al. The peripheral blood transcriptome identifies the presence and extent of disease in idiopathic pulmonary fibrosis. PLoS One. (2012) 7:e37708. doi: 10.1371/journal.pone.0037708
8. O’Reilly S, Ciechomska M, Cant R, van Laar JM. Interleukin-6 (IL-6) trans signaling drives a STAT3-dependent pathway that leads to hyperactive transforming growth factor-beta (TGF-beta) signaling promoting SMAD3 activation and fibrosis via Gremlin protein. J Biol Chem. (2014) 289:9952–60. doi: 10.1074/jbc.M113.545822
9. Iwamoto R, Mine N, Kawaguchi T, Minami S, Saeki K, Mekada E. HB-EGF function in cardiac valve development requires interaction with heparan sulfate proteoglycans. Development. (2010) 137:2205–14. doi: 10.1242/dev.048926
10. Bime C, Casanova N, Oita RC, Ndukum J, Lynn H, Camp SM, et al. Development of a biomarker mortality risk model in acute respiratory distress syndrome. Crit Care. (2019) 23:410. doi: 10.1186/s13054-019-2697-x
11. Lynn H, Sun X, Casanova N, Gonzales-Garay M, Bime C, Garcia JGN. Genomic and genetic approaches to deciphering acute respiratory distress syndrome risk and mortality. Antioxid Redox Signal. (2019) 31:1027–52. doi: 10.1089/ars.2018.7701
12. Casanova N, Zhou T, Gonzalez-Garay M, Lussier Y, Sweiss N, Ma SF, et al. MicroRNA and protein-coding gene expression analysis in idiopathic pulmonary fibrosis yields novel biomarker signatures predicting survival. Transl Res. (2020) 228:1–12. doi: 10.1016/j.trsl.2020.07.009
13. Singla S, Zhou T, Javaid K, Abbasi T, Casanova N, Zhang W, et al. Expression profiling elucidates a molecular gene signature for pulmonary hypertension in sarcoidosis. Pulm Circ. (2016) 6:465–71. doi: 10.1086/688316
14. Camp SM, Ceco E, Evenoski CL, Danilov SM, Zhou T, Chiang ET, et al. Unique toll-like receptor 4 activation by NAMPT/PBEF induces NFkappaB signaling and inflammatory lung injury. Sci Rep. (2015) 5:13135. doi: 10.1038/srep13135
15. Garcia AN, Casanova NG, Kempf CL, Bermudez T, Valera DG, Song JH, et al. eNAMPT is a novel DAMP that contributes to the severity of radiation-induced lung fibrosis. Am J Respir Cell Mol Biol. (2022) 66:497–509. doi: 10.1165/rcmb.2021-0357OC
16. Pillai VB, Sundaresan NR, Kim G, Samant S, Moreno-Vinasco L, Garcia JG, et al. Nampt secreted from cardiomyocytes promotes development of cardiac hypertrophy and adverse ventricular remodeling. Am J Physiol Heart Circ Physiol. (2013) 304:H415–26. doi: 10.1152/ajpheart.00468.2012
17. Lucas R, Verin AD, Black SM, Catravas JD. Regulators of endothelial and epithelial barrier integrity and function in acute lung injury. Biochem Pharmacol. (2009) 77:1763–72. doi: 10.1016/j.bcp.2009.01.014
18. Kümpers P, Lukasz A, David S, Horn R, Hafer C, Faulhaber-Walter R, et al. Excess circulating angiopoietin-2 is a strong predictor of mortality in critically ill medical patients. Crit Care. (2008) 12:R147. doi: 10.1186/cc7130
19. Hutyrova B, Pantelidis P, Drabek J, Zurkova M, Kolek V, Lenhart K, et al. Interleukin-1 gene cluster polymorphisms in sarcoidosis and idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. (2002) 165:148–51. doi: 10.1164/ajrccm.165.2.2106004
20. Maver A, Medica I, Peterlin B. Search for sarcoidosis candidate genes by integration of data from genomic, transcriptomic and proteomic studies. Med Sci Monit. (2009) 15:SR22–8.
21. Bihl MP, Laule-Kilian K, Bubendorf L, Rutherford RM, Baty F, Kehren J, et al. Progressive pulmonary sarcoidosis–a fibroproliferative process potentially triggered by EGR-1 and IL-6. Sarcoidosis Vasc Diffuse Lung Dis. (2006) 23:38–50.
22. Homolka J, Muller-Quernheim J. Increased interleukin 6 production by bronchoalveolar lavage cells in patients with active sarcoidosis. Lung. (1993) 171:173–83. doi: 10.1007/BF00183946
23. Ando M, Miyazaki E, Abe T, Ehara C, Goto A, Masuda T, et al. Angiopoietin-2 expression in patients with an acute exacerbation of idiopathic interstitial pneumonias. Respir Med. (2016) 117:27–32. doi: 10.1016/j.rmed.2016.05.012
24. Fischer A, Ellinghaus D, Nutsua M, Hofmann S, Montgomery CG, Iannuzzi MC, et al. Identification of immune-relevant factors conferring sarcoidosis genetic risk. Am J Respir Crit Care Med. (2015) 192:727–36. doi: 10.1164/rccm.201503-0418OC
25. Pandey JP, Frederick M. TNF-alpha, IL1-beta, and immunoglobulin (GM and KM) gene polymorphisms in sarcoidosis. Hum Immunol. (2002) 63:485–91. doi: 10.1016/S0198-8859(02)00399-3
26. Zhou T, Casanova N, Pouladi N, Wang T, Lussier Y, Knox KS, et al. Identification of Jak-STAT signaling involvement in sarcoidosis severity via a novel microRNA-regulated peripheral blood mononuclear cell gene signature. Sci Rep. (2017) 7:4237. doi: 10.1038/s41598-017-04109-6
27. Quijada H, Bermudez T, Kempf CL, Valera DG, Garcia AN, Camp SM, et al. Endothelial eNAMPT amplifies pre-clinical acute lung injury: efficacy of an eNAMPT-neutralising monoclonal antibody. Eur Respir J. (2021) 57:2002536. doi: 10.1183/13993003.02536-2020
28. Bime C, Casanova NG, Nikolich-Zugich J, Knox KS, Camp SM, Garcia JGN. Strategies to DAMPen COVID-19-mediated lung and systemic inflammation and vascular injury. Transl Res. (2021) 232:37–48. doi: 10.1016/j.trsl.2020.12.008
29. Bastarache JA, Koyama T, Wickersham NE, Ware LB. Validation of a multiplex electrochemiluminescent immunoassay platform in human and mouse samples. J Immunol Methods. (2014) 408:13–23. doi: 10.1016/j.jim.2014.04.006
30. Arena ET, Rueden CT, Hiner MC, Wang S, Yuan M, Eliceiri KW. Quantitating the cell: turning images into numbers with ImageJ. Wiley Interdiscip Rev Dev Biol. (2017) 6:e260. doi: 10.1002/wdev.260
31. Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, et al. A comparison of background correction methods for two-colour microarrays. Bioinformatics. (2007) 23:2700–7. doi: 10.1093/bioinformatics/btm412
32. Casanova NG, Gonzalez-Garay ML, Sun B, Bime C, Sun X, Knox KS, et al. Differential transcriptomics in sarcoidosis lung and lymph node granulomas with comparisons to pathogen-specific granulomas. Respir Res. (2020) 21:321. doi: 10.1186/s12931-020-01537-3
33. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. (2017) 45:D362–8. doi: 10.1093/nar/gkw937
34. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. (2019) 47:D607–13. doi: 10.1093/nar/gky1131
35. Craig VJ, Zhang L, Hagood JS, Owen CA. Matrix metalloproteinases as therapeutic targets for idiopathic pulmonary fibrosis. Am J Respir Cell Mol Biol. (2015) 53:585–600. doi: 10.1165/rcmb.2015-0020TR
36. Rangel JR, Chung Y, Rosenbach M, Ziober AF, Zhang PJ, Schaffer A. Expression of Notch signaling components in cutaneous foreign body and sarcoidal granulomas and fusing macrophages. Am J Dermatopathol. (2014) 36:409–13. doi: 10.1097/DAD.0b013e3182a730cb
37. Garcia AN, Casanova NG, Valera DG, Sun X, Song JH, Kempf CL, et al. Involvement of eNAMPT/TLR4 signaling in murine radiation pneumonitis: protection by eNAMPT neutralization. Transl Res. (2021) 239:44–57. doi: 10.1016/j.trsl.2021.06.002
38. Tanriverdi E, Iliaz S, Cortuk M, Turan D, Chousein EGU, Gul S, et al. Evaluation of serum biomarkers in patients with sarcoidosis: can visfatin be a new biomarker for sarcoidosis? Turk Thorac J. (2020) 21:145–9.
39. Funk RS, Singh R, Pramann L, Gigliotti N, Islam S, Heruth DP, et al. Nicotinamide phosphoribosyltransferase attenuates methotrexate response in juvenile idiopathic arthritis and in vitro. Clin Transl Sci. (2016) 9:149–57. doi: 10.1111/cts.12399
40. Grenier P, Chevret S, Beigelman C, Brauner MW, Chastang C, Valeyre D. Chronic diffuse infiltrative lung disease: determination of the diagnostic value of clinical data, chest radiography, and CT and Bayesian analysis. Radiology. (1994) 191:383–90. doi: 10.1148/radiology.191.2.8153310
41. Hernandez-Bartolome A, Lopez-Rodriguez R, Borque MJ, Gonzalez-Moreno L, Real-Martinez Y, Garcia-Buey L, et al. Angiopoietin-2/angiopoietin-1 as non-invasive biomarker of cirrhosis in chronic hepatitis C. World J Gastroenterol. (2016) 22:9744–51. doi: 10.3748/wjg.v22.i44.9744
42. Mehta VB, Besner GE. HB-EGF promotes angiogenesis in endothelial cells via PI3-kinase and MAPK signaling pathways. Growth Factors. (2007) 25:253–63. doi: 10.1080/08977190701773070
43. Hirabayashi M, Asano Y, Yamashita T, Miura S, Nakamura K, Taniguchi T, et al. Possible pro-inflammatory role of heparin-binding epidermal growth factor-like growth factor in the active phase of systemic sclerosis. J Dermatol. (2018) 45:182–8. doi: 10.1111/1346-8138.14088
44. Henry MT, McMahon K, Mackarel AJ, Prikk K, Sorsa T, Maisi P, et al. Matrix metalloproteinases and tissue inhibitor of metalloproteinase-1 in sarcoidosis and IPF. Eur Respir J. (2002) 20:1220–7. doi: 10.1183/09031936.02.00022302
45. White GE, Cotterill A, Addley MR, Soilleux EJ, Greaves DR. Suppressor of cytokine signalling protein SOCS3 expression is increased at sites of acute and chronic inflammation. J Mol Histol. (2011) 42:137–51. doi: 10.1007/s10735-011-9317-7
46. Landolt L, Furriol J, Babickova J, Ahmed L, Eikrem Ø, Skogstrand T, et al. AXL targeting reduces fibrosis development in experimental unilateral ureteral obstruction. Physiol Rep. (2019) 7:e14091. doi: 10.14814/phy2.14091
47. Steiner CA, Rodansky ES, Johnson LA, Berinstein JA, Cushing KC, Huang S, et al. AXL is a potential target for the treatment of intestinal fibrosis. Inflamm Bowel Dis. (2021) 27:303–16. doi: 10.1093/ibd/izaa169
48. Venkataraman T, Frieman MB. The role of epidermal growth factor receptor (EGFR) signaling in SARS coronavirus-induced pulmonary fibrosis. Antiviral Res. (2017) 143:142–50. doi: 10.1016/j.antiviral.2017.03.022
49. Qin CC, Liu YN, Hu Y, Yang Y, Chen Z. Macrophage inflammatory protein-2 as mediator of inflammation in acute liver injury. World J Gastroenterol. (2017) 23:3043–52. doi: 10.3748/wjg.v23.i17.3043
50. Al-Alwan LA, Chang Y, Mogas A, Halayko AJ, Baglole CJ, Martin JG, et al. Differential roles of CXCL2 and CXCL3 and their receptors in regulating normal and asthmatic airway smooth muscle cell migration. J Immunol. (2013) 191:2731–41. doi: 10.4049/jimmunol.1203421
51. Liu X, Wang Y, Lu H, Li J, Yan X, Xiao M, et al. Genome-wide analysis identifies NR4A1 as a key mediator of T cell dysfunction. Nature. (2019) 567:525–9. doi: 10.1038/s41586-019-0979-8
52. Palumbo-Zerr K, Zerr P, Distler A, Fliehr J, Mancuso R, Huang J, et al. Orphan nuclear receptor NR4A1 regulates transforming growth factor-β signaling and fibrosis. Nat Med. (2015) 21:150–8. doi: 10.1038/nm.3777
53. Zeng X, Yue Z, Gao Y, Jiang G, Zeng F, Shao Y, et al. NR4A1 is Involved in Fibrogenesis in Ovarian Endometriosis. Cell Physiol Biochem. (2018) 46:1078–90. doi: 10.1159/000488838
54. Ahmed M, Zaghloul N, Zimmerman P, Casanova NG, Sun X, Song JH, et al. Endothelial eNAMPT drives EndMT and preclinical PH: rescue by an eNAMPT-neutralizing mAb. Pulm Circ. (2021) 11:20458940211059712. doi: 10.1177/20458940211059712
55. Bermudez T, Sammani S, Song JH, Hernon VR, Kempf CL, Garcia AN, et al. eNAMPT neutralization reduces preclinical ARDS severity via rectified NFkB and Akt/mTORC2 signaling. Sci Rep. (2022) 12:696. doi: 10.1038/s41598-021-04444-9
Keywords: biomarker, sarcoidosis, fibrosis, plasma, gene expression
Citation: Casanova NG, Reyes-Hernon V, Gregory T, Sun B, Bermudez T, Hufford MK, Oita RC, Camp SM, Hernandez-Molina G, Serrano JR, Sun X, Fimbres J, Mirsaeidi M, Sammani S, Bime C and Garcia JGN (2022) Biochemical and genomic identification of novel biomarkers in progressive sarcoidosis: HBEGF, eNAMPT, and ANG-2. Front. Med. 9:1012827. doi: 10.3389/fmed.2022.1012827
Received: 09 August 2022; Accepted: 07 October 2022;
Published: 25 October 2022.
Edited by:
John-David Aubert, Centre Hospitalier Universitaire Vaudois (CHUV), SwitzerlandReviewed by:
Camillo Ribi, Centre Hospitalier Universitaire Vaudois (CHUV), SwitzerlandNancy Lin, University of Colorado, United States
Copyright © 2022 Casanova, Reyes-Hernon, Gregory, Sun, Bermudez, Hufford, Oita, Camp, Hernandez-Molina, Serrano, Sun, Fimbres, Mirsaeidi, Sammani, Bime and Garcia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joe G. N. Garcia, skipgarcia@arizona.edu