Free Access
Issue
Med Sci (Paris)
Volume 34, October 2018
Cancer biomarkers
Page(s) 121 - 133
DOI https://doi.org/10.1051/medsci/201834f121
Published online 07 November 2018

© 2018 médecine/sciences – Inserm

Introduction

Sarcoidosis is a multi-organ disease characterized by the presence of immune granulomas with uncertain etiology [1, 2]. The pulmonary involvements account for 90% to 95% of cases [3], and skin lesions are present in at least 20% of sarcoidosis cases [4, 5] and are the initial disease manifestation in nearly one third of these subjects [6]. The severity of sarcoidosis may range from an asymptomatic stage to a life-threatening condition. In a minority of cases, sarcoidosis may cause fibrosis of involved organs, leading to permanent organ damage that will not respond spontaneously or poorly to treatments [7- 9]. Sarcoidosis is rarely fatal with an age- and sex-adjusted mortality rate of less than 5/1,000,000 in the United States. However, it appears that the death rate from sarcoidosis is rising [10]. The majority of deaths result from pulmonary sarcoidosis where pulmonary fibrosis and/or pulmonary hypertension developed [11, 12], although it usually progresses slowly over decades [13].

Mounting research has focused on the pathologic mechanisms of sarcoidosis. It is suspected that exposure to one or more extrinsic antigens in a genetically susceptible individual leads to the overactivation of inflammatory pathways that promote the formation of sarcoidal granulomas. Studies have found an elevated risk of sarcoidosis in individuals who have been exposed to microbial agents [14, 15] and environmental substances [16, 17]. Disease susceptibility may be genetically determined, and several genes affecting the prevalence and course of sarcoidosis have been identified. Specifically, human leukocyte antigen (HLA) genes have been shown to influence the development and progression of sarcoidosis [18]. Cytokines like TNF-α, IL-12 and interferon gamma (IFN-γ) have been revealed to be involved in the formation of sarcoidosis [1, 19]. TGF-beta/Smad pathway is involved in the pathogenesis of pulmonary sarcoidosis [20]. The known mechanisms include inflammation, granuloma formation, genetic polymorphism etc., but the primary cause of sarcoidosis remains unclear.

Modern experimental technologies, such as microarray, have facilitated research in discovering novel pathogenetic mechanisms of sarcoidosis. Large amounts of information, specifically with regard to the microarray-based mRNA expression analysis of pathological tissues including lung, lymph nodes and blood cells of sarcoidosis patients, have been used to find potential pathogenic mediators and novel biomarkers [9, 21] such as Matrix metalloproteinase 12 (MMP12). Unsupervised bioinformatics analyses of gene expression data can identify highly regulated molecular pathways that are likely to promote abnormal granulomatous inflammation. Up-regulated hypoxia inducible factor-1A (HIF-1A) and vascular endothelial growth factor (VEGF) genes have been linked to acknowledged negative prognostics [22].

In this study, we utilized the gene expression profiles of different staged cutaneous sarcoidosis (shorten as CSC) and pulmonary sarcoidosis (abbreviated as PSC) as well as normal controls, to identify specific and mutual gene expression patterns as well as signaling pathways involved in disease occurrence and development and to search relevant biomarkers, reliable endpoints as well as potential efficient treatments.

Material and methods

Study datasets

The gene expression profiles of pulmonary sarcoidosis (PSC) were downloaded from the GEO [23] (http://www.ncbi.nlm.nih.gov/geo/). In GSE19976, there were eight patients with active but self-limiting pulmonary sarcoidosis and seven with active, progressive, fibrotic pulmonary disease [9]. The gene expression profile was based on the GPL6244 [HuGene-1_0-st] Affymetrix Human Gene 1.0 ST Array [transcript (gene) version]. GSE16538 was based on the GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 plus 2.0 Array including six active pulmonary sarcoidosis (n=6) and six normal lung anatomy [24]. In GSE32887, there were 15 samples of sarcoidosis with skin lesions, 11 unaffected skin from patients with cutaneous sarcoidosis (CSC) and 5 normal samples [21]. We marked the self-limiting pulmonary sarcoidosis as “SL_PSC”, the fibrotic pulmonary sarcoidosis as “F_PSC”, sarcoidosis with skin lesions as “TL_CSC” and unaffected skin from patients with CSC as “NL_CSC”.

Data analysis

With the “MetaDE” package of R project, we merged the GSE19976 and GSE16538, and finally acquired 8 samples of SL_PSC, 13 samples of F_PSC and 6 normal samples. To explore the pathogenesis of PSC, we considered the SL_PSC and F_PSC as the test groups and the normal samples as control groups. Through the “limma” package of R project, the genes with p-value <0.05 were considered as the differentially expressed genes (DEGs). We performed the function enrichment to find the significant pathways and biological functions with the “Database for Annotation, Visualization and Integrated Discovery” (DAVID) (http://david.abcc.ncifcrf.gov/). The functions with p-values for enrichment scores less than 0.05 were extracted.

Functional dynamic deviation scoring

Differentially expressed genes of the three groups were linked to different biological progresses. We used the Euclidean distance of quantitative method to calculate each significant GO term relative to the normal state of the skin, the dynamic functional variety under the stages of non-lesion and lesion sarcoidosis [25].

A ( P ) = 1 N i = 1 N ( X i - Y i ) 2 $$ A(P)=\frac{1}{N}\sum_{i=1}^N\sqrt{{({X}_i-{Y}_i)}^2} $$

A (P) represented deviation score of function P; N showed the number of differentially expressed genes in that function; Xi was the mean expression level of gene i in the non-lesion/lesion sarcoidosis; Yi was the mean expression level of gene i in normal skin. Through calculating the Euclidean distance and cumulative sum of all the genes in function P under the non-lesion/lesion sarcoidosis and normal skin, the degree of function P deviating from the normal state under the non-lesion/lesion sarcoidosis was figured out.

Feature selection

The intersections of disease-related genes and disease progression-related genes between PSC and CSC were applied as features, the gene expression levels were applied as feature values and both were combined with machine learning RFE algorithm for feature selection. Feature set was obtained through combination of a number of iterative features extracted in the process of selection. Feature set number increased gradually, thus making it possible to evaluate efficiency of the highest corresponding feature set, enabling feature selection.

Construction of a SVM classification model

The genes optimized by feature selection were applied as biomarkers to classify patients with PSC or CSC, using the classification of support vector machine (SVM) model. A five-fold cross-validation was performed. All samples were randomly divided into five groups, four of which being used as training set, the other one as test set. Each sample was then used as test set for prediction. A ROC curve was finally used to evaluate the effectiveness of the model.

Results

Analysis of mRNA expression profiles

We found 1,903 differentially expressed genes (DEGs) between CSC patients and normal control, and 1,498 DEGs between TL_CSC and NL_CSC patients. For the PSC, there were 1,058 DEGs between disease and normal samples, and 80 DEGs between F_PSC and SL_PSC samples. The DEGs are shown in Table 1.

Table 1.

The differentially expressed genes in different stages of sarcoidosis.

Twelve DEGs were significantly expressed at different stages of CSC ( Figure 1A) whereas 2 genes were found significantly deregulated at different stages of PSC ( Figure 1B). Between CSC and PSC, we found 97 common disease-related genes ( Figure 1C), and 13 common progression-related genes ( Figure 1D).

thumbnail Figure 1.

Numbers of DEGs at different stages of sarcoidosis. (A) DEGs of CSC; (B) DEGs of PSC; (C) relationships of DEGs between CSC and PSC; (D) relationships of process disease DEGs between CSC and PSC.

Functional analysis

Functional enrichment analysis was performed to better understand the role of these DEGs. Our data show that DEGs between F_PSC and normal tissue are significantly enriched in genes encoding proteins involved in cellular protein catabolic process, protein catabolic process, modification-dependent macromolecule catabolic process, sterol biosynthetic process, etc. ( Figure 2A). By contrast, DEGs between SL_PSC and normal tissue were significantly enriched in genes encoding molecules involved in chemotaxis, positive regulation of translation, angiogenesis, protein kinase cascade, etc. ( Figure 2C). Moreover, the DEGs in F_PSC and SL_PSC mainly regulate phosphorylation, protein kinase cascade, protein modification process, regulation of immune effector process etc. ( Figure 2B).

thumbnail Figure 2.

Biological process of DEGs regulated at different stages of sarcoidosis. (A) Functional enrichment of DEGs of F_ normal_PSC; (B) functional enrichment of DEGs of F_SL_PSC; (C) functional enrichment of DEGs of SL_normal_PSC; (D) functional enrichment of DEGs of NL_normal_CSC; (E) functional enrichment of DEGs of TL_normal_CSC; (F) functional enrichment of DEGs of TL_NLCSC. Vertical coordinates show the GO terms and horizontal ordinate the p-values.

Regarding to CSC, the DEGs between NL_CSC and normal tissue were significantly enriched in genes encoding DNA metabolic process, spermatogenesis, negative regulation of transcription etc. ( Figure 2D). By contrast, DEGs between TL_CSC and normal tissue were markedly enriched in immune cell-related functions including T cell activation, positive regulation of cell activation and of leukocyte activation etc. ( Figure 2E). Furthermore, DEGs in TL_CSC and NL_CSC mostly modulate T cell activation, positive regulation of immune response, positive regulation of cell activation etc. ( Figure 2F).

In addition, we also carried out a functional enrichment on the mutual disease-related genes and disease progression-related genes ( Figure 3). Results showed that the common disease-related genes between CSC and PSC are involved in the sterol biosynthetic process, steroid biosynthetic process and so on ( Figure 3A). The common progression-related genes mainly regulate the intracellular signaling cascade, phosphorylation, regulation of erythrocyte differentiation, protein modification process, protein kinase cascade ( Figure 3B).

thumbnail Figure 3.

Function analysis of (A) the common disease-related DEGs and of (B) common disease progression-related genes. Vertical coordinates represent the GO terms and horizontal ordinate was the p-values.

Dynamic functional changes

CSC

According to the deviation scores, 149 biological functions impacted were defined. These biological functions are differentially regulated during the transition from normal to non-invasive sarcoidosis and to skin lesion sarcoidosis (Supplement Table 1). The following biological processes changed markedly ( Figure 3A): immune response-activating signal transduction, immune response-regulating signal transduction, activation of immune response. The most significant differentially expressed genes are involved in apoptosis, T cell receptor signaling, B cell receptor signaling, viral myocarditis, antigen processing and presentation, Fcgamma R-mediated phagocytosis ( Figure 3B).

PSC

Changes in biological functions were also suggested when the gene expression between normal lung to minor lesion PSC and F_PSC was compared (Supplement Table 2). Biological processes such as the positive regulation of antigen receptor-mediated signaling, the regulation of antigen receptor-mediated signaling, the regulation of lymphocyte activation, the regulation of leukocyte activation, lymphocyte differentiation etc., were significantly impacted ( Figure 4A). Several signaling pathways are also remarkably changed, such as that involved in Primary immunodeficiency, in T cell receptor signaling pathway, in cytokine-cytokine receptor interaction, and in Natural Killer (NK) cell-mediated cytotoxicity ( Figure 4B).

thumbnail Figure 4.

The top ten dynamic functions in cutaneous sarcoidosis. (A) dynamic GO of CSC; (B) dynamic pathway of CSC.

Biomarker screening based on feature selection

Feature selection was performed on the intersections of disease-related genes and disease progression-related genes from PSC and CSC. Through iterative feature extraction, three biomarkers (CLDN7, NUDT11 and FAM73A) were obtained from the disease-related genes, and two biomarkers (TAP1 and TAP2) were acquired from the disease progression-related genes ( Figure 6AB). The selected biomarkers were used as features to construct SVM classifier for the supervised classification of patients with PSC or CSC. Combined with a five-fold cross-validation, a ROC curve was generated to evaluate the effectiveness of the model. Our results indicated that the classification model constructed on the basis of disease-related and disease progression-related biomarkers, optimized by feature selection, could effectively distinguish patients with PSC or CSC, with efficiencies of 87% and 92%, respectively ( Figure 6CD).

thumbnail Figure 5.

The top ten dynamic functions in pulmonary sarcoidosis. (A) dynamic GO of PSC; (B) dynamic pathway of PSC.

thumbnail Figure 6.

Biomarkers between pulmonary and cutaneous sarcoidosis. (A) feature selection for the common disease genes between CSC and PSC; (B) feature selection for the common processed disease genes between CSC and PSC; (C) ROC of SVM using disease-related biomarkers; (D) ROC of SVM using process-related biomarkers.

Discussion

Sarcoidosis is a multiple system granulomatous disease of unknown etiology [1, 2, 26]. The exact pathogenesis of sarcoidosis is still not fully understood and the link between the CSC and PSC remains unknown [27]. In the recent years, many researches have focused on the pathogenesis of CSC and PSC [26, 28]. Despite these efforts, it is still unclear how these two diseases correlate. Identifying the specific biomarkers closely associated with both CSC and PSC could largely contribute to early clinical diagnosis of the different types of sarcoidosis and to prevent CSC – PSC transition. Therefore, in this study, we have analyzed and compared the expression profiles of different staged CSC and PSC, identifying specific genes and the related impacted biological functions of the CSC and PSC, respectively. Moreover, the shared genes and functional pathways between CSC and PSC were also analyzed. 97 common disease-related genes and 13 mutual disease progression-related genes were identified. Further analysis showed that these common disease-related genes encode proteins mainly involved in sterol biosynthetic process and steroid biosynthetic process; by contrast, the mutual disease progression-related genes mostly regulate intracellular signaling cascade, phosphorylation, erythrocyte differentiation, protein modification process, and protein kinase cascade. We propose that these common genes and their related functions could represent the common features between CSC and PSC, and represent also potential markers for the diagnosis of the CSC-PSC transition. In addition, through feature selection, 3 disease-related biomarkers (CLDN7, NUDT11 and FAM73A) and 2 disease progression-related biomarkers (TAP1 and TAP2) were identified. Moreover, SVM classifier constructed on the basis of these features could effectively distinguish patients with PSC or CSC.

We considered that the shared genes are simultaneously involved in different types of sarcoidosis. Therefore, they are conservative genes significantly associated with sarcoidosis. Abnormal expression of these genes can lead to sarcoidosis. However, their expression is probably affected by other factors, such as distribution, genetic background [29- 31] and environmental stimulation [32], as well as other related genes [21, 33]. The intersection of DEGs embodies the universality of sarcoidosis while the specific genes defined in CSC and PSC represent the specificity of these two diseases. Because the specific genes are only differentially expressed in one type of sarcoidosis, these DEGs could represent clinically relevant diagnostic markers to distinguish different types of sarcoidosis. In addition, we also analyzed the differences in patients with sarcoidosis at different stages, identifying disease-related genes and progression-related genes. The abnormal expression of these genes fully displayed the dynamic nature of disease progression. Thus, utilizing these genes as diagnostic markers to evaluate a dynamic disease will contribute to identify the turning point of the disease, and could help preventing disease progression and improving the response to treatment.

During the functional analysis, we found that disease-related genes in both CSC and PSC are involved in a number of important biological functions, such as sterol biosynthesis and steroid biosynthetic process. During the occurrence of sarcoidosis, some important biological pathways are dysregulated. Correction of these abnormalities during the disease process could prevent disease progression and improve the efficacy of treatment. As the disease progresses, more downstream genes and biological functions are abnormally induced. Our results have identified common progression-related genes regulating intracellular signaling cascades, phosphorylation, regulation of erythrocyte differentiation, protein modification, and protein kinase cascade. Thus, the specific dysregulated functions probably reveal the biological mechanisms accompanying different types of sarcoidosis. This may help the identification of the etiology of the disease and the development of effective therapeutic strategies.

Data analyses have also demonstrated that several biological processes defined by DEGs contribute to the occurrence of CSC and/or to the development of the disease into cutaneous CSC. For example, T cell activation and the positive regulation of leukocyte activation may mark the occurrence of CSC. Moreover, T cell activation and a set of genes encoding proteins involved in the regulation of immune response were significantly differentially regulated by DEGs in TL_CSC and NL_CSC. It suggests a mechanism of development of the disease from non-lesion CSC to skin lesion CSC. Regarding to PSC, biological processes such as chemotaxis, positive regulation of translation, angiogenesis, protein kinase cascade, were remarkably differentially regulated by the DEGs in self-limiting PSC as compared to normal tissue, which may reveal the occurrence of PSC. DEGs in fibrotic PSC and normal tissue were enriched in genes encoding proteins involved in catabolic processes, in sterol biosynthetic process, showing the PSC progressed to a more serious stage. We also discovered that phosphorylation, protein kinase cascade, protein modification process, regulation of immune effector process, are dysregulated between F_PSC and SL_PSC, indicating that these biological processes play significant roles in the fibrotic progress of PSC. Consistent with our findings, previous researches have reported that persistent stimulation of the immune system plays a major role in granuloma formation and maintenance [26, 34]. Some immune-related cells such as macrophages, dendritic cells, lymphocytes and regulatory T cells, as well as immune-related factors are all closely involved in the cellular and molecular mechanisms of granuloma formation [8, 19, 35] and in the immune-pathogenesis of sarcoidosis [34]. Recent research indicated that immune-relevant factors such as IL-33 are important factors in sarcoidosis, making them potential biomarkers [36, 37]. Moreover, IFN-γ plays an important role in granuloma formation in this disease [37]. Although immunological events may determine the granuloma fate, other signaling pathways and factors also play important roles in pathogenesis and/or developmental process of sarcoidosis. Angiogenesis regulated by HIF-1A/VEGF/inhibitor of growth protein 4 (ING-4) axis may be crucial for the course and outcome of sarcoidosis [22]. Despite these advances, there are no clinically useful biomarkers that can assist the clinician in diagnosis, prognosis or assessment of treatment effects [17, 28, 38].

Disease progression is a dynamic process and we used a dynamic functional scoring algorithm to identify and analyze the differential functional pathways at different stages of disease progression. We found that functional pathways display a gradient dynamic change during disease progression. The functional pathways were mostly close to normal at the early stage of the disease, while the disease progression made the abnormalities of pathways increasing and finally greatly deviating from the normal levels. The correlations of function and disease stimuli were therefore ranked according to the extent of functional variation. The more obvious the dysfunction was, and the closer the correlation with the disease was. Thus, these dysregulated functions may be the key players in the mechanisms that drive the disease, and the genes encoding molecules involved in these functional pathways might represent therapeutic targets.

Furthermore, three disease-related biomarkers (CLDN7, NUDT11 and FAM73A) and two disease progression-related biomarkers (TAP1 and TAP2) were identified by feature selection. Claudin-7 has been previously reported to be expressed in sarcoidosis [39]. NUDT11 belongs to a subgroup of phosphohydrolases that preferentially attack diphosphoinositol polyphosphates [40]. Dysregulation of NUDT11 may contribute to the abnormal metabolism observed during the occurrence of sarcoidosis. However, there was no report on the role of FAM73A in sarcoidosis. The transporter associated with antigen processing (TAP) may be involved in determining sarcoidosis susceptibility [41, 42]. Previously, serum amyloid-A was reported to be potentially helpful for detecting the activity of sarcoidosis, and it has been also reported that soluble interleukin-2 receptor (sIL2-R) measurement could be useful in exploring the extra-pulmonary organ involvement [43]. Despite these, no single biomarker can be reliably used for a correct diagnosis or an exclusion of sarcoidosis [44,45]. Our results indicate that the use of these five biomarkers (CLDN7, NUDT11, FAM73A, TAP1, and TAP2) as features to construct SVM classifier, could effectively distinguish patients with PSC or CSC. Certainly, further clinical confirmations are needed.

To conclude, DEGs in pulmonary sarcoidosis and cutaneous sarcoidosis patients as well as in healthy individuals were analyzed to identify common disease-related genes and mutual disease progression-related genes which linked pulmonary sarcoidosis and cutaneous sarcoidosis together. DEGs and their biological functions are dynamically dysregulated. The biological functional pathways regulated by these DEGs may help defining the common mechanism shared by different types of sarcoidosis. They are encoding proteins likely involved in the pathogenesis of sarcoidosis and could serve as disease markers. Our study provides a novel insight into the common pathogenesis of sarcoidosis and could help developing new therapeutic strategies.

Disclosure statement

The authors disclose no conflict of interest.

References

  1. Valeyre D. et al. Sarcoidosis. Lancet 2014 ; 383 : 9923 1155–1167. [CrossRef] [PubMed] [Google Scholar]
  2. Iannuzzi M.C. Rybicki B.A. Teirstein A.S. Sarcoidosis. N Engl J Med 2007 ; 357 : 21 2153–2165. [CrossRef] [PubMed] [Google Scholar]
  3. Baughman R.P. et al. Clinical characteristics of patients in a case control study of sarcoidosis. Am J Respir Crit Care Med 2001 ; 164 : 10 Pt 1 1885–1889. [CrossRef] [PubMed] [Google Scholar]
  4. Mana J. et al. Cutaneous involvement in sarcoidosis. Relationship to systemic disease. Arch Dermatol 1997 ; 133 : 7 882–888. [Google Scholar]
  5. Lodha S. Sanchez M. Prystowsky S. Sarcoidosis of the skin: a review for the pulmonologist. Chest 2009 ; 136 : 2 583–596. [CrossRef] [PubMed] [Google Scholar]
  6. Costabel U. Guzman J. Baughman R.P. Systemic evaluation of a potential cutaneous sarcoidosis patient. Clin Dermatol 2007 ; 25 : 3 303–311. [CrossRef] [PubMed] [Google Scholar]
  7. Chappell A.G. Cheung W.Y. Hutchings H.A. Sarcoidosis: a long-term follow up study. Sarcoidosis Vasc Diffuse Lung Dis 2000 ; 17 : 2 167–173. [PubMed] [Google Scholar]
  8. Broos C.E. et al. Granuloma Formation in Pulmonary Sarcoidosis. Frontiers in Immunology 2013 ; 4 : 437. [CrossRef] [PubMed] [Google Scholar]
  9. Lockstone H.E. et al. Gene set analysis of lung samples provides insight into pathogenesis of progressive, fibrotic pulmonary sarcoidosis. Am J Respir Crit Care Med 2010 ; 181 : 12 1367–1375. [CrossRef] [PubMed] [Google Scholar]
  10. Swigris J.J. et al. Sarcoidosis-related mortality in the United States from 1988 to 2007. Am J Respir Crit Care Med 2011 ; 183 : 11 1524–1530. [CrossRef] [PubMed] [Google Scholar]
  11. Huang C.T. et al. Mortality in sarcoidosis. A changing pattern of the causes of death. Eur J Respir Dis 1981 ; 62 : 4 231–238. [PubMed] [Google Scholar]
  12. Arcasoy S.M. et al. Characteristics and outcomes of patients with sarcoidosis listed for lung transplantation. Chest 2001 ; 120 : 3 873–880. [CrossRef] [PubMed] [Google Scholar]
  13. Iwai K. et al. Racial difference in cardiac sarcoidosis incidence observed at autopsy. Sarcoidosis 1994 ; 11 : 1 26–31. [PubMed] [Google Scholar]
  14. Brownell I. et al. Evidence for mycobacteria in sarcoidosis. Am J Respir Cell Mol Biol 2011 ; 45 : 5 899–905. [CrossRef] [PubMed] [Google Scholar]
  15. Gupta D. et al. Molecular evidence for the role of mycobacteria in sarcoidosis: a meta-analysis. Eur Respir J 2007 ; 30 : 3 508–516. [CrossRef] [PubMed] [Google Scholar]
  16. Newman L.S. et al. A case control etiologic study of sarcoidosis: environmental and occupational risk factors. Am J Respir Crit Care Med 2004 ; 170 : 12 1324–1330. [CrossRef] [PubMed] [Google Scholar]
  17. Chen E.S. Moller D.R. Sarcoidosis–scientific progress and clinical challenges. Nat Rev Rheumatol 2011 ; 7 : 8 457–467. [CrossRef] [PubMed] [Google Scholar]
  18. Zhou Y. et al. Human leukocyte antigen-A, -B, and -DRB1 alleles and sarcoidosis in Chinese Han subjects. Hum Immunol 2011 ; 72 : 7 571–575. [CrossRef] [PubMed] [Google Scholar]
  19. Grunewald J. Eklund A. Role of CD4+ T cells in sarcoidosis. Proc Am Thorac Soc 2007 ; 4 : 5 461–464. [CrossRef] [Google Scholar]
  20. Piotrowski W.J. et al. TGF-beta and SMADs mRNA Expression in Pulmonary Sarcoidosis. Adv Exp Med Biol 2015 ; 852 : 59–69. [CrossRef] [PubMed] [Google Scholar]
  21. Judson, M.A., et al., Molecular profiling and gene expression analysis in cutaneous sarcoidosis: the role of interleukin-12, interleukin-23, and the T-helper 17 pathway. J Am Acad Dermatol, 2012. 66(6): p. 901–10, 910.e1-2. [CrossRef] [PubMed] [Google Scholar]
  22. Piotrowski, W.J., et al., Expression of HIF-1A/VEGF/ING-4 Axis in Pulmonary Sarcoidosis. Adv Exp Med Biol, 2015. [Google Scholar]
  23. Imai H. Kaira K. Minato K. Clinical significance of post-progression survival in lung cancer. Thorac Cancer 2017 ; 8 : 379–386. [CrossRef] [PubMed] [Google Scholar]
  24. Crouser E.D. et al. Gene expression profiling identifies MMP-12 and ADAMDEC1 as potential pathogenic mediators of pulmonary sarcoidosis. Am J Respir Crit Care Med 2009 ; 179 : 10 929–938. [CrossRef] [PubMed] [Google Scholar]
  25. Eickhoff S.B. et al. Analysis of neurotransmitter receptor distribution patterns in the cerebral cortex. Neuroimage 2007 ; 34 : 4 1317–1330. [CrossRef] [PubMed] [Google Scholar]
  26. Chen, E.S. and D.R. Moller, Etiologies of Sarcoidosis. Clin Rev Allergy Immunol, 2015. [Google Scholar]
  27. Ma Y. Gal A. Koss M.N. The pathology of pulmonary sarcoidosis: update. Semin Diagn Pathol 2007 ; 24 : 3 150–161. [CrossRef] [PubMed] [Google Scholar]
  28. Morgenthau A.S. Iannuzzi M.C. Recent advances in sarcoidosis. Chest 2011 ; 139 : 1 174–182. [CrossRef] [PubMed] [Google Scholar]
  29. Rybicki B.A. et al. Familial risk ratio of sarcoidosis in African-American sibs and parents. Am J Epidemiol 2001 ; 153 : 2 188–193. [CrossRef] [Google Scholar]
  30. Rybicki B.A. et al. Racial differences in sarcoidosis incidence: a 5-year study in a health maintenance organization. Am J Epidemiol 1997 ; 145 : 3 234–241. [CrossRef] [Google Scholar]
  31. Sverrild A. et al. Heredity in sarcoidosis: a registry-based twin study. Thorax 2008 ; 63 : 10 894–896. [CrossRef] [PubMed] [Google Scholar]
  32. Newman L.S. et al. A case control etiologic study of sarcoidosis: environmental and occupational risk factors. Am J Respir Crit Care Med 2004 ; 170 : 12 1324–1330. [CrossRef] [PubMed] [Google Scholar]
  33. Crouser E.D. et al. Differential expression of microRNA and predicted targets in pulmonary sarcoidosis. Biochem Biophys Res Commun 2012 ; 417 : 2 886–891. [CrossRef] [Google Scholar]
  34. Loke W.S.J. Herbert C. Thomas P.S. Sarcoidosis: Immunopathogenesis and Immunological Markers. International Journal of Chronic Diseases 2013 ; 2013 : 13. [Google Scholar]
  35. Ringkowski S. Thomas P.S. Herbert C. Interleukin-12 family cytokines and sarcoidosis. Front Pharmacol 2014 ; 5 : 233. [PubMed] [Google Scholar]
  36. Naumnik, W., et al., Interleukin-33 as a New Marker of Pulmonary Sarcoidosis. Adv Exp Med Biol, 2015. [Google Scholar]
  37. Fischer, A., et al., Identification of Immune-relevant Factors Conferring Sarcoidosis Genetic Risk. Am J Respir Crit Care Med, 2015. [Google Scholar]
  38. Mana J. et al. Granulomatous cutaneous sarcoidosis: diagnosis, relationship to systemic disease, prognosis and treatment. Sarcoidosis Vasc Diffuse Lung Dis 2013 ; 30 : 4 268–281. [PubMed] [Google Scholar]
  39. Kaarteenaho-Wiik R. Soini Y. Claudin-1, -2, -3, -4, -5, and -7 in Usual Interstitial Pneumonia and Sarcoidosis. Journal of Histochemistry and Cytochemistry 2009 ; 57 : 3 187–195. [CrossRef] [Google Scholar]
  40. Hidaka K. et al. An adjacent pair of human NUDT genes on chromosome X are preferentially expressed in testis and encode two new isoforms of diphosphoinositol polyphosphate phosphohydrolase. J Biol Chem 2002 ; 277 : 36 32730–32738. [CrossRef] [PubMed] [Google Scholar]
  41. Foley P.J. et al. Analysis of MHC encoded antigen-processing genes TAP1 and TAP2 polymorphisms in sarcoidosis. Am J Respir Crit Care Med 1999 ; 160 : 3 1009–1014. [CrossRef] [PubMed] [Google Scholar]
  42. Foley P. et al. HLA-DPB1 and TAP1 polymorphisms in sarcoidosis. Chest 1997 ; 111 : 6 Suppl 73S. [CrossRef] [PubMed] [Google Scholar]
  43. Gungor S. et al. Conventional markers in determination of activity of sarcoidosis. Int Immunopharmacol 2015 ; 25 : 1 174–179. [CrossRef] [PubMed] [Google Scholar]
  44. Wessendorf T.E. Bonella F. Costabel U. Diagnosis of Sarcoidosis. Clin Rev Allergy Immunol 2015 ; 49 : 1 54–62. [CrossRef] [PubMed] [Google Scholar]
  45. Li S. Sun X. Miao S. Liu J. Jiao W. Differential protein-coding gene and long noncoding RNA expression in smoking-related lung squamous cell carcinoma. Thorac Cancer 2017 ; 8 : 672–681. [CrossRef] [PubMed] [Google Scholar]

All Tables

Table 1.

The differentially expressed genes in different stages of sarcoidosis.

All Figures

thumbnail Figure 1.

Numbers of DEGs at different stages of sarcoidosis. (A) DEGs of CSC; (B) DEGs of PSC; (C) relationships of DEGs between CSC and PSC; (D) relationships of process disease DEGs between CSC and PSC.

In the text
thumbnail Figure 2.

Biological process of DEGs regulated at different stages of sarcoidosis. (A) Functional enrichment of DEGs of F_ normal_PSC; (B) functional enrichment of DEGs of F_SL_PSC; (C) functional enrichment of DEGs of SL_normal_PSC; (D) functional enrichment of DEGs of NL_normal_CSC; (E) functional enrichment of DEGs of TL_normal_CSC; (F) functional enrichment of DEGs of TL_NLCSC. Vertical coordinates show the GO terms and horizontal ordinate the p-values.

In the text
thumbnail Figure 3.

Function analysis of (A) the common disease-related DEGs and of (B) common disease progression-related genes. Vertical coordinates represent the GO terms and horizontal ordinate was the p-values.

In the text
thumbnail Figure 4.

The top ten dynamic functions in cutaneous sarcoidosis. (A) dynamic GO of CSC; (B) dynamic pathway of CSC.

In the text
thumbnail Figure 5.

The top ten dynamic functions in pulmonary sarcoidosis. (A) dynamic GO of PSC; (B) dynamic pathway of PSC.

In the text
thumbnail Figure 6.

Biomarkers between pulmonary and cutaneous sarcoidosis. (A) feature selection for the common disease genes between CSC and PSC; (B) feature selection for the common processed disease genes between CSC and PSC; (C) ROC of SVM using disease-related biomarkers; (D) ROC of SVM using process-related biomarkers.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.