- Open Access
Differential transcriptomics in sarcoidosis lung and lymph node granulomas with comparisons to pathogen-specific granulomas
Respiratory Research volume 21, Article number: 321 (2020)
Despite the availability of multi-“omics” strategies, insights into the etiology and pathogenesis of sarcoidosis have been elusive. This is partly due to the lack of reliable preclinical models and a paucity of validated biomarkers. As granulomas are a key feature of sarcoidosis, we speculate that direct genomic interrogation of sarcoid tissues, may lead to identification of dysregulated gene pathways or biomarker signatures.
To facilitate the development sarcoidosis genomic biomarkers by gene expression profiling of sarcoidosis granulomas in lung and lymph node tissues (most commonly affected organs) and comparison to infectious granulomas (coccidiodomycosis and tuberculosis).
Transcriptomic profiles of immune-related gene from micro-dissected sarcoidosis granulomas within lung and mediastinal lymph node tissues and compared to infectious granulomas from paraffin-embedded blocks. Differentially-expressed genes (DEGs) were profiled, compared among the three granulomatous diseases and analyzed for functional enrichment pathways.
Despite histologic similarities, DEGs and pathway enrichment markedly differed in sarcoidosis granulomas from lymph nodes and lung. Lymph nodes showed a clear immunological response, whereas a structural regenerative response was observed in lung. Sarcoidosis granuloma gene expression data corroborated previously reported genomic biomarkers (STAB1, HBEGF, and NOTCH4), excluded others and identified new genomic markers present in lung and lymph nodes, ADAMTS1, NPR1 and CXCL2. Comparisons between sarcoidosis and pathogen granulomas identified pathway divergences and commonalities at gene expression level.
These findings suggest the importance of tissue and disease-specificity evaluation when exploring sarcoidosis genomic markers. This relevant translational information in sarcoidosis and other two histopathological similar infections provides meaningful specific genomic-derived biomarkers for sarcoidosis diagnosis and prognosis.
Sarcoidosis is a multisystemic disease of unknown etiology characterized by the formation of granulomatous lesions, especially in lung tissues and thoracic lymph nodes [1, 2]. The clinical heterogeneity and unpredictable disease course in sarcoidosis represents a challenge in early diagnosis and prediction of progression. Despite significant progress in the understanding of genetic predisposition and role of immunity in sarcoidosis pathogenesis, predicting the clinical course of sarcoidosis in a given patient remains challenging. The majority of sarcoidosis patients experience spontaneous recovery, however, fully a third of subjects develop complicated sarcoidosis defined as multiorgan involvement, progressive fibrotic lung involvement and pulmonary function deterioration [3, 4]. Validated sarcoidosis biomarkers are desperately needed to distinguish patients who will spontaneously recover from patients who will worsen and develop severe manifestations of the disease. Thus, personalized medicine approaches are tasked with the capacity to predict the patient phenotype likely to develop progressive disease or that may eventually require lung transplantation. While there has been progress in identifying genetic variants that contribute to complicated sarcoidosis [5, 6], and we have previously utilized genomic expression profiling of peripheral mononuclear cells (PBMCs) to subphenotype patients with a variety of inflammatory disorders [7,8,9] including sarcoidosis , there remains an important and unmet need for molecular and specific genome-based biomarkers for diagnosis and prediction of disease severity in sarcoidosis.
The etiology of sarcoidosis is unknown, and the diagnosis remains one of exclusion . Specific infections, including a variety of Mycobacterium species and Propionibacterium acnes [11,12,13,14,15,16,17,18,19,20], and environmental exposures  produce granulomatous inflammatory lesions that mimic sarcoidosis and have been posited as potential etiologic agents. Exposure to mycobacterium tuberculosis (TB) and endemic fungi such coccidioidomycosis (CM), also known as Valley Fever, are major causes of granulomatous lung disease  and must be clinically excluded when sarcoidosis is suspected. Coccidioidomycosis caused by Coccidioides immitis and C. posadasii, a soil-dwelling fungi disease, are endemic in the arid southwestern USA and northern Mexico.
Granulomas are conglomerates of epithelioid and multinucleated giant cells surrounded by CD4+ and CD8+ T lymphocytes that result from the complex immunopathogenesis between host genetic factors, environmental or infectious triggers . Sarcoidosis array-based transcriptomic studies in blood, lung and lymph node, lacrimal gland, orbital tissue and skin identified tissue-specific differential gene expression [22,23,24,25], highlighting the importance of assessing different compartments in a multisystemic disease such as sarcoidosis. The majority of these transcriptome studies were focused on identifying gene signatures that differentiate controls from progressive forms of sarcoidosis. Additional studies in peripheral blood compared the gene expression of sarcoidosis and TB, identifying a significant overlap in gene expression associated with Type I and II IFN pathways . Approaches to quantify the expression of sarcoidosis granuloma-related genes using RT-PCR focused solely on cytokine-related genes associated with NFKB and STAT transcription factors .
As expression profiling of sarcoidosis granuloma using next generation sequencing expression has not been reported, one goal of the present study was to compare gene expression profiles in two sarcoidosis granulomatous tissues: lung and lymph node. The second goal of this study was to assess the ability of granuloma gene expression profiles to discriminate sarcoidosis from TB and CM. Our bioinformatic analyses included determination of unique pathways associated with sarcoidosis and overlapping pathways with CM and TB. Our results indicate that despite histologic similarities, sarcoidosis lung and lymph node granulomas exhibit distinct expression profiles. In addition, sarcoidosis pathway expression was significantly divergent compared to TB and CM. These findings suggest the importance of tissue-specific pathobiology considerations to be employed when exploring potential sarcoidosis genomic markers for potential diagnostic use and discriminate from other granulomatous diseases.
Materials and methods
Sample selection and design
We used microdissected granulomatous tissue from lung  and mediastinal lymph node  from 18 subjects with sarcoidosis, three with CM, four with TB and control tissue from lungs  and lymph nodes  from 6 healthy individuals. Patient and sample characteristics are described in Table 1. Tissue specimens were collected from clinically-indicated biopsies for diagnostic purposes. Archived, de-identified specimens were acquired from the Tissue Acquisition Shared Resource at the University of Arizona. The study was approved by the human subjects protection program IRB # 1509097312A001.
Gene expression profiling of selected granulomas
Formalin-fixed, paraffin-embedded (FFPE) microdissections from lung and lymph node granulomatous tissue and healthy tissue (controls) were assayed using Next Generation Sequencing-based gene expression by HTG EdgeSeq Oncology-biomarker (HTG Molecular Diagnostics, Inc.). This system uses targeted capture sequencing to quantitate RNA expression levels of gene targets in FFPE tissues. The panel included 2535 probes and 15 housekeeper genes for quantitative analysis of targeted mRNAs. A list of the complete biomarker panel genes and the biological pathways assayed are in Additional file 1: Table S1. This panel covers multiple immune-oncologic-related pathways previously published and identified as relevant in sarcoidosis including such as interferon, MAP kinase, NFКB, and JAK/STAT pathways [18,19,20]. Representative tissue sections from granulomatous tissue were stained with hematoxylin- and eosin (H&E). Sites for granuloma microdissection were selected by an expert pathologist (Additional file 2). Sections of ~ 1.5 mm2 thickness of a 5 μm FFPE tissue were obtained, 5–10 sections per subject were lysed followed by RNA extraction-free chemistry method using the HTG®’s assay kit. Library preparation for the HTG Edge processor for nuclease protection steps was used, followed by PCR tagging, library amplification, quantitation and normalization. Library was then sequenced on an Illumina MiSeq using 150-cycle V3 kits .
Analysis of biomarker transcriptomes and canonical pathways
We generated the transcriptome index using the salmon package  and a fasta file containing the sequence of the HTG EdgeSeq probes. Quantification was performed for each sample using the salmon quant command , the index file and the corresponding fastq files. The raw counts were loaded into R using the tximport package . The Differential Expression (DE) was calculated using Limma  and EdgeR  packages from the R v3.5.3 and Bioconductor v3.8 packages . A list with statistically-significant genes for each comparison was submitted for pathway enrichment unbiased comparison of the DE gene sets to the human ConsensusPathDB  website using the over-representation gene set analysis against the pathway databases and the gene ontology categories using a minimum overlap with input list of 2 and a p-value cutoff of 0.01.
Transcriptome dataset validation
We validated our transcriptome results in three independent cohorts of microarray datasets publicly available in the Gene Expression Omnibus (GEO), GSE16538 and GSE63548  and an unpublished microarray dataset, details in Additional file 3: Table S3. At the time this manuscript was written, no GEO data sets were available for humans infected with CM. The microarray datasets for TB and sarcoidosis were analyzed using R (v3.5.3) and Bioconductor (v3.8) . Oligo package (v1.46.0) was used for the processing of the data in combination with specific platform design library and the corresponding annotation library. For the Affymetrix array analyses, Oligo’s implementation of the Robust Multichip Average algorithm (RMA) for background correction, quantile normalization and summarization was used. For the Illumina arrays, the Illuminaio package  was used in combination with the corresponding bead Chip. The differential expression analysis was performed using the Limma package . After normalization gene expression was averaged and each tissue was compared against heathy tissue from the same origin. The inverse log2 (FC) generated by Limma was used to determine differences between comparisons, with the direction of change indicated by the sign. Transcripts with a fold change (FC) > 2 and a p value < 0.01 were deemed differentially expressed.
Microdissected lung and lymph node granulomas from 18 sarcoidosis subjects, four TB subjects, three CM subjects and six healthy controls were assessed. There were no significant differences in age (p < 0.1) in the three granulomatous disease groups (Table 1). The majority of sarcoidosis subjects were female and non-Hispanic whites while the demographic most frequent reported in the CM and TB groups was Hispanic. Clinical data related to organ involvement and radiographic Scadding stage was available for 70% of sarcoidosis subjects as detailed in Additional file 1: Table S1.
After excluding probes with multiple hits or partial hits, our gene expression data in granulomatous tissues included 2430 probes that uniquely mapped to single genes. We averaged gene expression levels after normalization and compared diseased groups versus healthy tissue controls from the same tissue origin, the fold change (FC) was obtained as the ratio between the averages. We identified 250 significantly dysregulated genes between all granulomatous-diseased tissues and healthy control tissues. Figure 1 represents the differentially expressed genes (DEGs) volcano plots, in each category; granulomas from CM exhibited the fewest number of differentiated transcripts, 34 in total (Fig. 1a), sarcoidosis granulomas from the lung (Fig. 1b) had more DEGs than those granulomas from lymph nodes (Fig. 1c), 88 vs 60 dysregulated transcripts respectively. Notably, in both sarcoidosis-affected tissues, the number of downregulated genes (FC < − 2) outpaced the number of upregulated genes (FC > 2), with downregulation more marked in lung granulomas (93% of the transcripts) than in lymph node granulomas (56% of the transcripts). Tuberculosis granulomas (Fig. 1d), showed the highest number of the differentiated transcripts, 140 in total. The heat maps comparing the gene expression in each group against healthy tissues are presented in Fig. 2. In addition, we screened the total of 250 DEGs among those four sets performing a Venn diagram analysis . Figure 2e shows the overlapping and unique transcripts in the different granulomas. Sarcoidosis granulomas exhibited 138 dysregulated transcripts in both lung (Fig. 2a) and lymph node (Fig. 2b), 87 transcripts were exclusive present in sarcoidosis, notoriously, 30% of the DEGs in lung granulomas were unique to sarcoidosis, and the overlapping transcripts dysregulated in both sarcoidosis tissues was only 4%. TB had 89 DEG that were exclusively dysregulated in TB granulomas. When compared TB to sarcoidosis 17% of the transcripts were dysregulated in both diseases, while only 4% of the transcripts were common between CM and lung sarcoidosis.
Comparisons of expression profiles in sarcoidosis granulomas
We identified that sarcoidosis granulomas from lung and lymph nodes share only ten DEGs. NR1H3 and CXCL9 were upregulated whereas the remainder (ADAMTS1, CXCL2, HSPB6, ITGA9, NPR1, NR4A1, CCL14, FABP4) showed downregulation in both lung and lymph nodes. Interrogation of the expression levels of these 10 sarcoidosis genes to TB and CM DEGs revealed CCL14, CXCL9, FABP4, NR1H3 were also significantly dysregulated in TB. Only six genes (ADAMTS1, HSPB6, NR4A1, CXCL2, NPR1, ITGA9) were exclusively dysregulated in both lung and lymph node in sarcoidosis granulomas but not in CM or TB (Fig. 3). These results are very consistent with the notion that granuloma gene expression is highly tissue- and disease-specific.
Comparison of DEGs in sarcoidosis vs CM and TB granulomas
Comparison of TB lymph node granulomas expression profiles to lymph node healthy controls identified 140 DEGs with a surprising number of DEGs common to sarcoidosis, 43 transcripts, 39 present in lymph granulomas. Unlike sarcoidosis where DEGs were predominantly downregulated, TB DEGs showed a balanced proportion of expression, 64 down-regulated and 76 upregulated. Notably, we observed the same direction in the transcript regulation, in all the 43-shared transcripts, except for one, OLR1, a low density lipoprotein receptor involved in Fas-induced apoptosis, this gene was upregulated in TB and down regulated in sarcoidosis, over-expression of OLR1 result in upregulation of NF-κB and its target pro-oncogenes .
Gene expression profiles from CM lung granulomas exhibited the fewest number of dysregulated genes, 34 DEGs with a balanced pattern of dysregulation, 18 DEGs downregulated, 16 DEGs upregulated. A total of 15 DEGs were exclusively dysregulated in CM granulomas. Eight upregulated transcripts were common to sarcoidosis lung granulomas (CSF3, SERTAD1, MMP9, CCL19, BCL3, AREG, PTGS, F3). Only three DEGs, SLAMF7, DC27, and CXCL13, were dysregulated in all three diseases (sarcoidosis, TB, CM). Expression of CXCL13, a chemokine and B-lymphocyte chemoattractant associated to calcium influx was downregulated in both sarcoidosis (lymph node) and TB granulomas but upregulated CM granulomas. A list of the top dysregulated genes by granuloma tissue origin is available in Additional file 5: Table S4.
Functional biologic pathway enrichment analysis
We first identified the pathways that were common to all three granulomatous diseases using the gene list with all the DEG (250), cytokine-cytokine receptor interaction, chemokine signaling, VEGFA-VEGFR2 and focal adhesion accounted for the most significant pathways with the most number of proportionally dysregulated genes involved, we then mapped the pathways associated to the genes that were exclusively present in sarcoidosis. Sarcoidosis pathways showed similar altered pathways to the rest of the granulomatous diseases with the exception of nuclear receptors meta-pathway (Table 2). Then we compared the differences in pathways between lung and lymph nodes in sarcoidosis, using the gene set of DEG for each tissue, differences were delineated by a clear immunological response, involving leucocyte migration and neutrophil chemotaxis in the lymph nodes while a structural regenerative response, characterized by cell migration and angiogenesis, was observed at the lung level (Fig. 4a). At the disease level, we identified pathways that were disease specific for each gene set as seen in Fig. 4b.
Independent validation of disease-specific granuloma transcriptome results
We further validated our results using independent cohorts of microarray expression in lymph node and lung tissues from sarcoidosis and tuberculosis (Additional file 5: Fig. S1). Validation analysis in these microarrays data sets confirmed the dysregulation of 90 genes (FC > 2 FDR < 0.012) of the DEG identified in our panel. Among those dysregulated transcripts, we validated 46 in TB and 44 in sarcoidosis, 14 of them differentially expressed in lung, 30 in lymph nodes and 4 genes dysregulated in both tissues ADAMST1, CXCL2, CXCL9, FABP4. The direction of the expression defined by fold change was confirmed in all the dysregulated genes in the microarray data.
Sarcoidosis is a heterogeneous disease lacking validated clinical markers that allow differentiation from other granulomatous diseases or that predict the severity or level of organ involvement. In this study, we unveil the transcriptome profiles of cells comprising granulomas from sarcoidosis, CM, and TB. Using the HTG EdgeSeq oncology biomarker panel system to sequence FFPE biopsied tissues collected at the time of the diagnosis, we avoid limitations such as low RNA yields associated to FFPE extractions. We generated information to further define the molecular mechanisms that occur at the compartment level in sarcoidosis and compared to diseases that share a common histopathological hallmark.
Sarcoidosis granulomas showed differences in gene pathway expression between compartments potentially associated to disease progression. Lymph node granulomas exhibited clear immunological responses, involving leucocyte migration and neutrophil chemotaxis corresponding to early clinical stages observed in sarcoidosis. A structural regenerative response was observed in lung granulomas, notoriously marked by cell migration and angiogenesis, that appeared to correspond to lung remodeling stages observed in complicated and chronic sarcoidosis.
Only ten DEGs were common to both sarcoidosis lung granulomas and lymph nodes granulomas. ADAMTS1, CXCL2, HSPB6, ITGA9, NPR1, NR4A1, CCL14, CXCL9, FABP4 and NR1H3. Interestingly NR1H3, a regulator of macrophage function and inflammation, and CXCL9 a chemokine ligand involved in T-cell trafficking, were upregulated in both granulomas independently of the tissue origin (FC > 2.6). However, CXCL9 was also upregulated in the TB granulomas. Downregulated genes included FABP4, a fatty acid-binding protein associated with activation of proinflammatory macrophages in breast cancer, obesity and leukemia progression [39, 40]. FABP4, was previously proposed as sarcoidosis candidate gene [22, 41]. We identified that it was also significantly downregulated in TB. CCL14, a chemokine ligand associated with changes in intracellular calcium concentration in monocytes. Hypercalcemia associated to granulomatous diseases is commonly reported [42,43,44], the potential pathogenic role of CXCL13 and CXCL14, both DEG in sarcoidosis and TB needs further investigation.
Despite similarities at the histological level, our analysis revealed sarcoidosis granulomas to be significantly divergent at the genomic level when compared to CM granulomas. TB and sarcoidosis granuloma profiles share stronger similarity at the transcriptional and pathway level. Interestingly, the number of DE genes in each comparison using the same thresholds varied significantly, indicating that the number of dysregulated genes does not correlate with the number of samples per comparison, or the size of the microdissection, suggesting that the granulomas are qualitatively distinct depending on the disease and tissue of origin. We identified only three shared dysregulated transcripts in CM, TB and sarcoidosis (lymph nodes); CXCL13, CD27 and SLAMF7, a finding that points to diverse pathogenic mechanisms in granulomatous diseases that result in a similar histopathological conglomeration of CD4+ T cells surrounded by CD8+ T cells, fibroblasts and B cells. CXCL13 is a chemokine with key role in B cell migration by regulation of Ca++ influx , while SLAMF7 is a regulator of T lymphocyte development and function such as lytic activity and a modulator of B cell activation and memory . Furthermore, CD27 regulates B-cell activation and immunoglobulin synthesis through the activation of NF-kappaB and MAPK8/JNK. Our microarray validation data corroborated SLAMF7 DE but only in sarcoidosis lymph node granulomas.
The activation of Th1 cells capable of producing interferon (IFN)-γ is important in the immunological pathogenesis of the granuloma formation . Interferon-inducible neutrophil-driven blood transcriptional signature was previously associated to TB and one of the two sarcoidosis profiles (weaker vs strong IFN-inducible profile), with a higher abundance and expression in tuberculosis . Our panel included probes for IFN related genes such as IFNGR1, IRF1, IFI27, IFIT2, IFNA2, IFNAR1, IFNB1, IFNG and IFN regulatory related genes; IRF2, IRF3, IRF4, IRF5, IRF7 and IRF8. Only IFI27 was significantly upregulated in TB (p 0.0000678 FC > 6) but didn’t achieve significance in Sarcoidosis. IFI27 is involved in type-I interferon-induced apoptosis characterized by a rapid release of cytochrome C, related pathways include innate immune system and interferon gamma signaling. The IFN immune related pathways that contribute to sarcoidosis etiopathogenesis are reported to be exacerbated by immunotherapy with IFN alpha [49,50,51]. Further investigation is needed to explain the pathogenic effect of the observed differences at the gene expression level.
This study also served to demonstrate the lack of specificity of prior suggested potential sarcoidosis genomic markers due to common dysregulation in fungal or mycobacterium granulomas. For example, CXCL9, a chemoattractant for lymphocytes previously proposed as a marker of sarcoidosis severity , was also over-expressed in TB and is also dysregulated in Beryllium disease . On the other hand, we confirmed downregulation of NOTCH4 expression to be exclusively present in lung sarcoidosis. NOTCH signaling is initiated after activation of toll-like receptors in macrophages and it has been previously reported that Notch receptor ligand Dll4 is expressed in cutaneous sarcoidosis granulomas . Supporting this association, SNPs in NOTCH4 have been reported as a sarcoidosis-associated locus in genome-wide associated study in African Americans .
We have previously generated a 20-sarcoidosis gene signature by microarray in peripheral blood from PBMCs . We now confirm three genes (SERTAD1, HBEGF, KLRB1) as DEGs by oncopanel profiling. KLRB1 expression was also downregulated in TB and SERTAD1 was significantly downregulated in CM. Therefore, these two genes appear to be part of common process in granulomatous diseases. In contrast, HBEGF, a heparin-binding epidermal growth factor, was down-regulated in sarcoidosis lung granulomas (FC − 2.36, p = 0.0016) and was not dysregulated in CM or TB. HBEGF is involved in epithelialization, wound contraction and angiogenesis with an active pro-inflammatory role in skin and lung alveolar regeneration , pneumonitis and early stages of systemic sclerosis . Importantly, within the sarcoidosis PBMC molecular signature, HBEGF expression was up-regulated in contrast to the observed down-regulation in granuloma tissues. This variation in gene expression across tissues may be explained by either heritability at cis loci, by trans regulation across tissues [53, 54] or by epigenetic mechanisms. Our previously conducted epigenetic studies in sarcoidosis identified IL6ST and STAB1 as part of a miRNA-derived gene signature for complicated sarcoidosis . In the current study, IL6ST was significantly dysregulated only in TB granulomas (FC − 1.42 and p 0.0008), while STAB1 gene was significantly downregulated in lymph nodes in sarcoidosis (FC − 1.64, p 0.0014). STAB1 codes for a transmembrane receptor that is expressed in endothelial cells and lymph nodes, with functions in angiogenesis, lymphocyte homing and cell adhesion which are key aspects in the chronicity of the granuloma.
Leveraging the power derived from publicly available microarray data sets from independent cohorts, we successfully validated ADAMST1, CXCL2 and FABP4, that were identified DEGs by the oncopanel. These three genomic markers were present exclusively in sarcoidosis in both lung and lymph node. NPR1, Natriuretic Peptide Receptor 1, is associated with cardiac hypertrophy and fibrosis  and has not been previously associated to sarcoidosis, was present in the lung and lymph node transcriptomic profiles and was also validated in microarray datasets but only in sarcoidosis granulomas from lymph node. ADAMTS1, gene encodes a disintegrin and metalloproteinase with thrombospondin motif associated with various inflammatory processes and a potentially important regulator of pathogenic granuloma formation in sarcoidosis in both compartments.
There are some limitations of the current study. One is related to RNA from FFPE tissues, despite the use of quantitative nuclease protection chemistry that enabled extraction-free quantitation of mRNA, RNA quantity and quality are limited in FFPE compared to cryopreserved tissues. Other limitations include the use of a sequencing panel and not whole genome sequencing, that limited number of gene probes and the expression levels of some genes was undetected. Our comparisons were against tissues and not immune cells, that and the difference in platforms (RNAseq vs microarrays) can explain the limited overlap to already published data sets. Last, the only gene expression data available in CM was microarray data in mice strains limiting the validation of our results in CM. Despite the observed differences at the pathway level in sarcoidosis granulomas at lung and lymph node that indicate disease progression, further investigations are needed to corroborate the expression of these genes in sarcoidosis heterogeneous subphenotypes since limited clinical data was available to establish these correlations.
Sarcoidosis is a complex systemic disease of unknown cause with clinical care hampered by the paucity of available tools for diagnosis and prognosis. Our study is the first to identify RNAseq-derived transcriptome differences in granulomas from lung and lymph nodes of sarcoidosis subjects compared to CM and TB. Our data suggest that immunological signaling related-pathways distinguishes sarcoidosis from other granulomatous diseases and fills a crucial gap in knowledge of gene expression within granulomas. This study corroborates previous genomic markers suggested for sarcoidosis such as STAB1, HBEGF, FABP4 and NOTCH4, with HBEGF and NOTCH4 only expressed in lung granulomas. We have identified new genomic markers for sarcoidosis, ADAMTS1, CXCL2, and NPR1. This relevant translational information in sarcoidosis and other two histopathological similar infections provides important information for the future development of meaningful specific genomic-derived biomarkers for sarcoidosis diagnosis and prognosis.
Databases generated from this study are available at GEO: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE157671
Differentially Expressed Genes
Gene Expression Omnibus
False Discovery Rate
Soto-Gomez N, Peters JI, Nambiar AM. Diagnosis and Management of Sarcoidosis. Am Fam Physician. 2016;93(10):840–8.
Judson MA. The clinical features of sarcoidosis: a comprehensive review. Clin Rev Allergy Immunol. 2015;49(1):63–78. https://doi.org/10.1007/s12016-014-8450-y.
Lockstone HE, Sanderson S, Kulakova N, Baban D, Leonard A, Kok WL, McGowan S, McMichael AJ, Ho LP. 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–75. https://doi.org/10.1164/rccm.200912-1855OC.
Chen ES, Moller DR. Sarcoidosis–scientific progress and clinical challenges. Nat Rev Rheumatol. 2011;7(8):457–67. https://doi.org/10.1038/nrrheum.2011.93.
Schupp JC, Vukmirovic M, Kaminski N, Prasse A. Transcriptome profiles in sarcoidosis and their potential role in disease prediction. Curr Opin Pulm Med. 2017;23(5):487–92. https://doi.org/10.1097/MCP.0000000000000403.
Casanova N, Zhou T, Knox KS, Garcia JG. Identifying novel biomarkers in sarcoidosis using genome-based approaches. Clin Chest Med. 2015;36(4):621–30. https://doi.org/10.1016/j.ccm.2015.08.005.
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(14):1027–52. https://doi.org/10.1089/ars.2018.7701.
Zhang W, Zhou T, Ma SF, Machado RF, Bhorade SM, Garcia JG. MicroRNAs implicated in dysregulation of gene expression following human lung transplantation. Transl Respir Med. 2013. https://doi.org/10.1186/2213-0802-1-12.
Huang Y, Ma SF, Vij R, Oldham JM, Herazo-Maya J, Broderick SM, Strek ME, White SR, Hogarth DK, Sandbo NK, Lussier YA, Gibson KF, Kaminski N, Garcia JG, Noth I. A functional genomic model for predicting prognosis in idiopathic pulmonary fibrosis. BMC Pulm Med. 2015;15:147. https://doi.org/10.1186/s12890-015-0142-8.
Zhou T, Zhang W, Sweiss NJ, Chen ES, Moller DR, Knox KS, Ma SF, Wade MS, Noth I, Machado RF, Garcia JG. Peripheral blood gene expression as a novel genomic biomarker in complicated sarcoidosis. PLoS ONE. 2012;7(9):e44818. https://doi.org/10.1371/journal.pone.0044818.
Brownell I, Ramirez-Valle F, Sanchez M, Prystowsky S. Evidence for mycobacteria in sarcoidosis. Am J Respir Cell Mol Biol. 2011;45(5):899–905. https://doi.org/10.1165/rcmb.2010-0433TR.
Drake WP, Dhason MS, Nadaf M, Shepherd BE, Vadivelu S, Hajizadeh R, Newman LS, Kalams SA. Cellular recognition of Mycobacterium tuberculosis ESAT-6 and KatG peptides in systemic sarcoidosis. Infect Immun. 2007;75(1):527–30. https://doi.org/10.1128/iai.00732-06.
Eishi Y, Suga M, Ishige I, Kobayashi D, Yamada T, Takemura T, Takizawa T, Koike M, Kudoh S, Costabel U, Guzman J, Rizzato G, Gambacorta M, du Bois R, Nicholson AG, Sharma OP, Ando M. Quantitative analysis of mycobacterial and propionibacterial DNA in lymph nodes of Japanese and European patients with sarcoidosis. J Clin Microbiol. 2002;40(1):198–204. https://doi.org/10.1128/jcm.40.1.198-204.2002.
Goldmann T, Zissel G, Sen Gupta R, Schlaak M, Vollmer E, Muller-Quernheim J. Formation of granulomas in the lungs of severe combined immunodeficient mice after infection with bacillus Calmette-Guerin. Am J Pathol. 2001;158(5):1890–1 (Epub 2001/05/05).
Goyal B, Sheikh JA, Agarwal R, Verma I. Levels of circulating immune complexes containing Mycobacterium Tuberculosis-specific antigens in pulmonary tuberculosis and sarcoidosis patients. Indian J Med Microbiol. 2017;35(2):290–2. https://doi.org/10.4103/0255-0857.209569 (Epub 2017/07/07).
Gupta D, Agarwal R, Aggarwal AN, Jindal SK. Molecular evidence for the role of mycobacteria in sarcoidosis: a meta-analysis. Eur Respir J. 2007;30(3):508–16. https://doi.org/10.1183/09031936.00002607 (Epub 2007/06/01).
Mortaz E, Adcock IM, Barnes PJ. Sarcoidosis: Role of non-tuberculosis mycobacteria and Mycobacterium tuberculosis. Int J Mycobacteriol. 2014;3(4):225–9. https://doi.org/10.1016/j.ijmyco.2014.10.008 (Epub 2016/01/21).
Oswald-Richter K, Sato H, Hajizadeh R, Shepherd BE, Sidney J, Sette A, Newman LS, Drake WP. Mycobacterial ESAT-6 and katG are recognized by sarcoidosis CD4+ T cells when presented by the American sarcoidosis susceptibility allele, DRB1*1101. J Clin Immunol. 2009;30(1):157–66. https://doi.org/10.1007/s10875-009-9311-y (Epub 2009/06/19).
Oswald-Richter KA, Culver DA, Hawkins C, Hajizadeh R, Abraham S, Shepherd BE, Jenkins CA, Judson MA, Drake WP. Cellular responses to mycobacterial antigens are present in bronchoalveolar lavage fluid used in the diagnosis of sarcoidosis. Infect Immun. 2009;77(9):3740–8. https://doi.org/10.1128/IAI.00142-09 (Epub 2009/07/15).
Ebe Y, Ikushima S, Yamaguchi T, Kohno K, Azuma A, Sato K, Ishige I, Usui Y, Takemura T, Eishi Y. Proliferative response of peripheral blood mononuclear cells and levels of antibody to recombinant protein from Propionibacterium acnes DNA expression library in Japanese patients with sarcoidosis. Sarcoidosis Vasc Diffuse Lung Dis. 2000;17(3):256–65.
Crowley LE, Herbert R, Moline JM, Wallenstein S, Shukla G, Schechter C, Skloot GS, Udasin I, Luft BJ, Harrison D, Shapiro M, Wong K, Sacks HS, Landrigan PJ, Teirstein AS. “Sarcoid like” granulomatous pulmonary disease in World Trade Center disaster responders. Am J Ind Med. 2011;54(3):175–84. https://doi.org/10.1002/ajim.20924 (Epub 2011/02/08).
Rosenbaum JT, Choi D, Wilson DJ, Grossniklaus HE, Harrington CA, Sibley CH, Dailey RA, Ng JD, Steele EA, Czyz CN, Foster JA, Tse D, Alabiad C, Dubovy S, Parekh P, Harris GJ, Kazim M, Patel P, White V, Dolman P, Korn BS, Kikkawa D, Edward DP, Alkatan H, Al-Hussain H, Yeatts RP, Selva D, Stauffer P, Planck SR. Parallel gene expression changes in sarcoidosis involving the lacrimal gland, orbital tissue, or blood parallel gene expression changes in sarcoidosis parallel gene expression changes in sarcoidosis. JAMA Ophthalmol. 2015;133(7):770–7. https://doi.org/10.1001/jamaophthalmol.2015.0726.
Crouser ED, Culver DA, Knox KS, Julian MW, Shao G, Abraham S, Liyanarachchi S, Macre JE, Wewers MD, Gavrilin MA, Ross P, Abbas A, Eng C. 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–38. https://doi.org/10.1164/rccm.200803-490OC (Epub 2009/02/12).
Koth L, Solberg O, Peng J, Bhakta N, Nguyen C, Woodruff P. Sarcoidosis blood transcriptome reflects lung inflammation and overlaps with tuberculosis. Am J Respir Crit Care Med. 2011;184(10):1153–63.
Sheng Y, Yang Y, Wu Y, Yang Q. Exploring the dynamic changes between pulmonary and cutaneous sarcoidosis based on gene expression. Med Sci (Paris). 2018;34:121–33.
Christophi GP, Caza T, Curtiss C, Gumber D, Massa PT, Landas SK. Gene expression profiles in granuloma tissue reveal novel diagnostic markers in sarcoidosis. Exp Mol Pathol. 2014;96(3):393–9. https://doi.org/10.1016/j.yexmp.2014.04.006.
Ravi RK, Walton K, Khosroheidari M. MiSeq: a next generation sequencing platform for genomic analysis. Methods Mol Biol. 2018;1706:223–32. https://doi.org/10.1007/978-1-4939-7471-9_12 (Epub 2018/02/10).
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9. https://doi.org/10.1038/nmeth.4197 (Epub 2017/03/06).
Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. 2015;4:1521. https://doi.org/10.12688/f1000research.7563.2 (Epub 2015/12/30).
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. doi: https://doi.org/10.1093/nar/gkv007. Epub 2015/01/20.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97. https://doi.org/10.1093/nar/gks042 (Epub 2012/01/31).
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):80. doi: https://doi.org/10.1186/gb-2004-5-10-r80. Epub 2004/09/15.
Kamburov A, Stelzl U, Lehrach H, Herwig R. The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 2013;41:D793-800. https://doi.org/10.1093/nar/gks1055 (Epub 2012/11/11).
Maji A, Misra R, Kumar Mondal A, Kumar D, Bajaj D, Singhal A, Arora G, Bhaduri A, Sajid A, Bhatia S, Singh S, Singh H, Rao V, Dash D, Baby Shalini E, Sarojini Michael J, Chaudhary A, Gokhale RS, Singh Y. Expression profiling of lymph nodes in tuberculosis patients reveal inflammatory milieu at site of infection. Sci Rep. 2015;5:15214. https://doi.org/10.1038/srep15214 (Epub 2015/10/15).
Smith ML, Baggerly KA, Bengtsson H, Ritchie ME, Hansen KD. illuminaio: an open source IDAT parsing tool for Illumina microarrays. F1000Res. 2013;2:264. https://doi.org/10.12688/f1000research.2-264.v1 (Epub 2013/01/01).
Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK. A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007;23(20):2700–7.
Oliveiros V. An interactive tool for comparing lists with Venn's diagrams 2015 [cited 2019].
Khaidakov M, Mitra S, Kang BY, Wang X, Kadlubar S, Novelli G, Raj V, Winters M, Carter WC, Mehta JL. Oxidized LDL receptor 1 (OLR1) as a possible link between obesity, dyslipidemia and cancer. PLoS ONE. 2011;6(5):e20277. doi: https://doi.org/10.1371/journal.pone.0020277. Epub 2011/06/04.
Yan F, Shen N, Pang JX, Zhang YW, Rao EY, Bode AM, Al-Kali A, Zhang DE, Litzow MR, Li B, Liu SJ. Fatty acid-binding protein FABP4 mechanistically links obesity with aggressive AML by enhancing aberrant DNA methylation in AML cells. Leukemia. 2017;31(6):1434–42. https://doi.org/10.1038/leu.2016.349 (Epub 2016/11/26).
Hao J, Zhang Y, Yan X, Yan F, Sun Y, Zeng J, Waigel S, Yin Y, Fraig MM, Egilmez NK, Suttles J, Kong M, Liu S, Cleary MP, Sauter E, Li B. Circulating adipose fatty acid binding protein is a new link underlying obesity-associated breast/mammary tumor development. Cell Metab. 2018;28(5):689–705. https://doi.org/10.1016/j.cmet.2018.07.006 (Epub 2018/08/14).
Maver A, Medica I, Peterlin B. Search for sarcoidosis candidate genes by integration of data from genomic, transcriptomic and proteomic studies. Med Sci Monitor. 2009;15(12):22–8 (Epub 2009/12/01).
Baughman RP, Janovcik J, Ray M, Sweiss N, Lower EE. Calcium and vitamin D metabolism in sarcoidosis. Sarcoidosis Vasc Diffuse Lung Dis. 2013;30(2):113–20.
Peces R. Hypercalcemia in granulomatous diseases. Nefrologia. 2000;20(5):393–5.
Sharma OP. Hypercalcemia in granulomatous disorders: a clinical review. Curr Opin Pulm Med. 2000;6(5):442–7. https://doi.org/10.1097/00063198-200009000-00010.
Rupprecht TA, Plate A, Adam M, Wick M, Kastenbauer S, Schmidt C, Klein M, Pfister HW, Koedel U. The chemokine CXCL13 is a key regulator of B cell recruitment to the cerebrospinal fluid in acute Lyme neuroborreliosis. J Neuroinflammation. 2009;6:42. https://doi.org/10.1186/1742-2094-6-42 (Epub 2010/01/01).
Detre C, Keszei M, Romero X, Tsokos GC, Terhorst C. SLAM family receptors and the SLAM-associated protein (SAP) modulate T cell functions. Semin Immunopathol. 2010;32(2):157–71. https://doi.org/10.1007/s00281-009-0193-0 (Epub 2010/02/11).
Takimoto H, Maruyama H, Shimada KI, Yakabe R, Yano I, Kumazawa Y. Interferon-gamma independent formation of pulmonary granuloma in mice by injections with trehalose dimycolate (cord factor), lipoarabinomannan and phosphatidylinositol mannosides isolated from Mycobacterium tuberculosis. Clin Exp Immunol. 2006;144(1):134–41. https://doi.org/10.1111/j.1365-2249.2006.03043.x (Epub 2006/03/18).
Bloom CI, Graham CM, Berry MP, Rozakeas F, Redford PS, Wang Y, Xu Z, Wilkinson KA, Wilkinson RJ, Kendrick Y, Devouassoux G, Ferry T, Miyara M, Bouvry D, Valeyre D, Gorochov G, Blankenship D, Saadatian M, Vanhems P, Beynon H, Vancheeswaran R, Wickremasinghe M, Chaussabel D, Banchereau J, Pascual V, Ho LP, Lipman M, O’Garra A. Transcriptional blood signatures distinguish pulmonary tuberculosis, pulmonary sarcoidosis, pneumonias and lung cancers. PLoS ONE. 2013;8(8):e70630. https://doi.org/10.1371/journal.pone.0070630.
Monteagudo B, Grueiro MC, Vilas-Sueiro A, Campo-Cerecedo F. Papular sarcoidosis of the knees following treatment with interferon alpha and ribavirin in a woman with hepatitis C. Actas Dermosifiliogr. 2019;110(10):871–3. https://doi.org/10.1016/j.ad.2018.04.016 (Epub 2019/08/08).
Flores-Robles BJ, Sanguesa-Gomez CP, Barbadillo Mateos C, Roustan-Gullon LG, Kovtun I. Sarcoidosis developing after treatment with interferon alpha in a patient with hepatitis C: An association every physician should know. Rev Gastroenterol Mex. 2017;82(2):185–6. https://doi.org/10.1016/j.rgmx.2016.01.005 (Epub 2016/05/11).
Trien R, Cooper CJ, Paez D, Colon E, Ajmal S, Salameh H. Interferon-alpha-induced sarcoidosis in a patient being treated for hepatitis C. Am J Case Rep. 2014;15:235–8. https://doi.org/10.12659/AJCR.890180 (Epub 2014/06/06).
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(5):409–13. https://doi.org/10.1097/DAD.0b013e3182a730cb (Epub 2014/01/08).
Price AL, Helgason A, Thorleifsson G, McCarroll SA, Kong A, Stefansson K. Single-tissue and cross-tissue heritability of gene expression via identity-by-descent in related or unrelated individuals. PLoS Genet. 2011;7(2):e1001317. doi: https://doi.org/10.1371/journal.pgen.1001317. Epub 2011/03/09.
Consortium GT, Laboratory DA, Coordinating Center -Analysis Working G, Statistical Methods groups-Analysis Working G, Enhancing Gg, Fund NIHC, Nih/Nci, Nih/Nhgri, Nih/Nimh, Nih/Nida, Biospecimen Collection Source Site N, Biospecimen Collection Source Site R, Biospecimen Core Resource V, Brain Bank Repository-University of Miami Brain Endowment B, Leidos Biomedical-Project M, Study E, Genome Browser Data I, Visualization EBI, Genome Browser Data I, Visualization-Ucsc Genomics Institute UoCSC, Lead a, Laboratory DA, Coordinating C, management NIHp, Biospecimen C, Pathology, e QTLmwg, Battle A, Brown CD, Engelhardt BE, Montgomery SB. Genetic effects on gene expression across human tissues. Nature. 2017;550(7675):204–13. doi: https://doi.org/10.1038/nature24277. Epub 2017/10/13.
Ellmers LJ, Scott NJ, Piuhola J, Maeda N, Smithies O, Frampton CM, Richards AM, Cameron VA. Npr1-regulated gene pathways contributing to cardiac hypertrophy and fibrosis. J Mol Endocrinol. 2007;38(1–2):245–57. https://doi.org/10.1677/jme.1.02138 (Epub 2007/02/13).
Special thanks to Betsy M. Dennison, Jocelyn Fimbres and Paulo Pena from TACMASR.
University of Arizona (CDA22) Nancy G. Casanova; Foundation for Sarcoidosis Research (Mini-grant) Joe G. N. Garcia.
Ethics approval and consent to participate
De-identified specimens were acquired from the Tissue Acquisition Shared Resource. The study was granted a waiver of consent and it was approved by the University of Arizona Human Subjects Protection Program IRB # 1509097312A001.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Oncopanel gene list. List of the complete biomarker panel genes and the biological pathways assayed.
Sarcoidosis clinical phenotype. Scadding stage based on radiographic imaging which consist of Stage I: hilar enlargement only; Stage II: hilar enlargement plus interstitial lung disease; Stage III: interstitial lung disease and Stage IV: lung fibrosis. Lung involvement and other organs affected based on clinical data. All subjects have at least one reported clinical symptom related to lung involvement (shortness of breath, cough). NA. Information not available.
Microarrays data sets used for validation of the DEG identified by the oncopanel.
Sarcoidosis, TB and CM top common dysregulated genes. Common DEGs in granulomas derived this oncopanel in lung sarcoidosis and CM and in granulomas derived from lymph nodes from TB and Sarcoidosis.
Microarray validation. Comparison of granuloma DEGs identified by the oncopanel to microarray datasets. Heat maps representing expression profiles of validated gene sets in A) Sarcoidosis Lung 14 transcripts validated, B) Sarcoidosis Lymph nodes, 30 transcripts validated and C). TB lymph nodes, 46 transcripts validated. Transcripts differentially regulated showed a fold change ±2 & FDR<0.01.
About this article
Cite this article
Casanova, N.G., Gonzalez-Garay, M.L., Sun, B. et al. Differential transcriptomics in sarcoidosis lung and lymph node granulomas with comparisons to pathogen-specific granulomas. Respir Res 21, 321 (2020). https://doi.org/10.1186/s12931-020-01537-3
- Valley fever
- Gene expression