Next Article in Journal
Assessing the Thiamine Diphosphate Dependent Pyruvate Dehydrogenase E1 Subunit for Carboligation Reactions with Aliphatic Ketoacids
Next Article in Special Issue
The Aryl Hydrocarbon Receptor in Asthma: Friend or Foe?
Previous Article in Journal
Whole Genome Sequencing and Tn5-Insertion Mutagenesis of Pseudomonas taiwanensis CMS to Probe Its Antagonistic Activity Against Rice Bacterial Blight Disease
Previous Article in Special Issue
Dendritic Cells: Critical Regulators of Allergic Asthma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 22
10.3390/ijms21228640
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Whole Transcriptome Analysis of Myeloid Dendritic Cells Reveals Distinct Genetic Regulation in Patients with Allergies

by
Kijeong Lee
1,†,
Mi-Ryung Han
2,†,
Ji Woo Yeon
1,
Byoungjae Kim
1 and
Tae Hoon Kim
1,*
1
Department of Otorhinolaryngology-Head & Neck Surgery, College of Medicine, Korea University, Seoul 02841, Korea
2
Division of Life Sciences, College of Life Sciences and Bioengineering, Incheon National University, Incheon 22012, Korea
*
Author to whom correspondence should be addressed.
†
Co-first author: These authors contributed equally to this work.
Int. J. Mol. Sci. 2020, 21(22), 8640; https://doi.org/10.3390/ijms21228640
Submission received: 15 October 2020 / Revised: 9 November 2020 / Accepted: 11 November 2020 / Published: 16 November 2020
(This article belongs to the Special Issue Molecular Mechanisms of Asthma)
Browse Figures
Versions Notes

Abstract

:
Dendritic cells (DCs) play critical roles in atopic diseases, orchestrating both innate and adaptive immune systems. Nevertheless, limited information is available regarding the mechanism through which DCs induce hyperresponsiveness in patients with allergies. This study aims to reveal novel genetic alterations and future therapeutic target molecules in the DCs from patients with allergies using whole transcriptome sequencing. Transcriptome sequencing of human BDCA-3+/CD11c+ DCs sorted from peripheral blood monocytes obtained from six patients with allergies and four healthy controls was conducted. Gene expression profile data were analyzed, and an ingenuity pathway analysis was performed. A total of 1638 differentially expressed genes were identified at p-values < 0.05, with 11 genes showing a log2-fold change ≥1.5. The top gene network was associated with cell death/survival and organismal injury/abnormality. In validation experiments, amphiregulin (AREG) showed consistent results with transcriptome sequencing data, with increased mRNA expression in THP-1-derived DCs after Der p 1 stimulation and higher protein expression in myeloid DCs obtained from patients with allergies. This study suggests an alteration in the expression of DCs in patients with allergies, proposing related altered functions and intracellular mechanisms. Notably, AREG might play a crucial role in DCs by inducing the Th2 immune response.

1. Introduction

In atopic diseases, such as asthma and allergic rhinitis, CD4+ Th2-driven immune responses play a crucial role [1]. Dendritic cells (DCs) are important antigen-presenting cells that initiate allergic inflammation by orchestrating both the innate and adaptive immune systems [2,3]. In humans, heterogeneous subtypes of DCs include conventional CD11c+ myeloid DCs (mDCs) and CD123+ plasmacytoid DCs (pDCs). mDCs effectively present antigens to CD4+ and CD8+ T cells and induce high levels of pro-inflammatory chemokines, whereas pDCs are specialized in controlling overall airway inflammation by inducing regulatory T cells [4,5].
Immature DCs are mainly distributed in tissues where they are exposed to pathogens or allergens [6]. After being stimulated by various microbes, DCs migrate toward the lymphoid organs, undergo maturation, and acquire the ability to activate naïve T cells into effector T cells [7,8]. The differentiation of naïve T cells into Th1/Th17 cells is related to cell-mediated immunity, which is dominant in autoimmune inflammatory diseases, whereas Th2-skewed immune imbalance leads to the dysregulation of antibody-mediated immunity in allergic diseases, including asthma, atopic dermatitis, and allergic rhinitis [9,10,11]. Several studies have identified molecular mechanisms in DCs that elicit Th1, Th17, and Treg immune responses through pattern recognition receptors. Liu et al. revealed that the SIRT1-HIF1α signaling pathway modulates the differentiation of naïve T cells to either Th1 cells or Treg cells [12]. Huang et al. proposed the importance of MKP-1 signaling in DCs for the induction of Th1 and Th17 cell differentiation and p38α signaling for Th17 differentiation [13,14]. Furthermore, Suryawanshi et al. recently reported that canonical Wnt signaling in DCs regulates Th1 and Th17 cell differentiation [15]. Th2-mediated induction of naïve T cells has not been extensively investigated, and functional genomes or their interactions in DCs that lead to Th2-driven immune responses remain unknown.
Compared to other hybridization-based methods, such as microarray and real-time polymerase chain reaction (PCR), high-throughput RNA sequencing methods are advantageous because they allow the global analysis of entire transcriptomes, do not require previously defined probes, and are more sensitive [16,17]. Furthermore, cell-type-specific sequencing enables the discovery of novel transcriptomes and helps to identify previously unknown intracellular mechanisms of specific cells, particularly immune cells, including B cells, T cells, and macrophages, in cancers or autoimmune diseases [18,19,20]. There have been attempts to identify genetic biomarkers for allergies using next-generation sequencing (NGS) data obtained from the analysis of biopsy samples of patients with asthma and atopic dermatitis [21,22]. In addition, Sheu et al. analyzed microRNA expression and identified potential interactions with mRNA in asthmatic bronchial epithelial cells, and Seumois et al. conducted a transcriptomic analysis of Th2 cells to elucidate the molecular mechanisms underlying the pathogenesis of asthma and allergic rhinitis [23,24]. Nevertheless, genetic differences in antigen-presenting cells, such as DCs, between patients with allergies and controls have not been identified using NGS.
Therefore, in the present study, we conducted whole transcriptome sequencing of human peripheral blood mDCs from healthy participants and participants with allergies using NGS and a bioinformatics approach. Our overall aim was to reveal novel target molecules and identify genomic networks regulating DCs in patients with allergies using whole transcriptome sequencing. We believe that the results of the present study will provide perspectives for developing new therapeutic strategies to treat allergic diseases.

2. Results

2.1. mDCs from Patients with Allergies Showed 11 Differentially Expressed Genes (DEGs) Compared with Healthy Controls

We performed differential gene expression analysis of RNA sequencing data. A heat map of gene expression values was constructed (Figure 1A), and a volcano plot was generated to visualize the identified genes (Figure 1B) based on log2-fold change (Log2FC). We identified 1658 genes (771 upregulated and 887 downregulated genes) that showed differential expression with p-values < 0.05 in mDCs from patients with allergies versus healthy controls (Figure S1, Table S1). Among these genes, five upregulated and six downregulated genes in mDCs from patients with allergies showed meaningful results with p-values < 0.05 and |Log2FC| ≥ 1.5 (Table 1). The gene TPST1 showed the largest Log2FC value of −2.7 with a p-value = 0.026 in DCs from patients with allergies.

2.2. Altered Biological Functions and Molecular Networks Were Identified in mDCs from Patients with Allergies

UPTOHERE A DEG dataset was used for performing ingenuity pathway analysis (IPA, http://www.ingenuity.com) to examine the functional and molecular pathway enrichment of DEGs based on Fisher’s exact test. Canonical pathway and network analyses were performed based on functional annotations of the DEGs and corresponding bibliographic data to determine significantly altered signaling pathways and regulatory networks. The Ingenuity s Knowledge Base was used as the reference set. The activation score (z-score) was used to infer the activation state of transcriptional regulators (z-score < −2 and >2).
Enriched functional annotations of expression profiles were obtained following IPA with an input of the top 100 DEGs. The top 100 genes were yielded in the order of |Log2FC| among the genes with p-values < 0.05 (Table S2). Table 2 shows the most altered biological functions of the identified DEGs at the cellular level, including signal transduction, cell death, antigen presentation, tyrosine sulfation of proteins, and cell migration. As for canonical pathway analysis, pathways associated with innate immunity and antigen presentation, including the following, were identified: communication between innate and adaptive immune cells, the role of pattern recognition receptors in the recognition of bacteria and viruses, phagosome formation, and toll-like receptor signaling. Table 3 shows canonical pathway analysis results with associated molecules and ratio of the number of pathway-related molecules to the total number of molecules. These pathways possessed core signaling molecules among overlapping canonical pathways, making connections with pathways related to Th1/Th2 activation and prothrombin activation (Figure S2).
A total of 8 related networks were found under the IPA network based on functional annotations. The network with the highest score (score of 36) consisting of 20 focus molecules was associated with cell death/survival, cancer, and organismal injury/abnormality (Figure 2A). Extracellular signal-regulated kinase (ERK) 1/2 and protein kinase B (Akt) were the key nodes of this network. Other genes, such as amphiregulin (AREG), colony-stimulating factor 2 (CSF2), and the cluster of differentiation 1 (CD1c), were also involved in this network. The second-ranking network with a score of 18 included 12 focus molecules associated with cardiovascular system development/function, cell-to-cell signaling/interaction, and cellular movement (Figure 2B). The second network was composed of cytokines, chemokine receptors, and transmembrane receptors. Of the 11 DEGs identified, tyrosylprotein sulfotransferase 1 (TPST1), filamin B (FLNB) and thrombomodulin (THBD) were in the second network.

2.3. Validation of Sequencing Data Revealed Increased Expression of AREG in Der p 1-Stimulated THP-1-Derived DCs and DCs Sorted from Peripheral Blood Monocytes (PBMCs) of Patients with Alergies

Of the 11 discovered DEGs, we conducted experiments to verify three genes that appear to have a role in the top 2 functional networks: AREG, TPST1, FLNB.
The level of mRNA expression for each gene in THP-1-derived DCs was evaluated after Der p 1 stimulation. The expression of AREG was increased by up to 3.26-fold after 2 h of Der p 1 stimulation (p = 0.0034), and it gradually decreased 24 h after stimulation. (Figure 3A). Up to 12 h (1.87-fold, p-value = 0.0331) of stimulation, AREG expression was statistically higher than that of the control. The mRNA level of TPST-1 was decreased by 0.76-fold following 2 h of Der p 1 stimulation (p-value = 0.0238); next, it gradually increased up to 12 h of stimulation followed by decreased expression at 24 h (0.82-fold, p-value = 0.0538). In addition, FLNB mRNA expression increased immediately after Der p 1 stimulation; however, it gradually decreased, showing significantly decreased expression at 24 h of stimulation compared to the control (0.74-fold, p-value = 0.0011).
Western blot analysis was conducted to validate the expression of 3 genes in DCs sorted from PBMCs of patients with allergies and healthy participants. The protein levels of AREG was higher in allergic DCs than in control DCs, showing consistent results with the transcriptome sequencing results (Figure 3D). TPST-1 and FLNB expression was higher in allergic DCs than in transcriptome sequencing results.

3. Discussion

In the present study, whole transcriptomes of peripheral blood mDCs were sequenced successfully, and DEGs between patients with allergies and normal controls were determined. In addition, enriched functional annotation and network identification using IPA showed that cellular functions associated with cell survival, cell-to-cell signaling/interaction, and cellular movement are key mechanisms for immature DCs to induce allergic reactions. Among the selected DEGs, AREG was validated through experiments in THP-1-derived DCs stimulated with an allergen and western blot analysis in DCs extracted from allergy patients.
Among the 1658 DEGs with p-values < 0.05, five upregulated and six downregulated genes were identified with |Log2FC| ≥ 1.5. Most of these genes have not been reported to be associated with allergic DCs to date, except for some such as THBD and CD1c. THBD is known as blood dendritic cell antigen 3 (BDCA3) or CD141, antibodies against which were used for sorting mDCs in this study. Although all samples sequenced in this experiment were THBD (BDCA3)-positive DCs, the expression of THBD was significantly higher in DCs from patients with allergies than control DCs. This result was consistent with the results from a previous microarray-based study, which showed that THBD expression on house dust mite (HDM)-stimulated DCs from individuals with allergies was higher than that in the controls [25]. In this study, the authors introduced THBD as the gene with highest degree of differential expression among 156 DEGs identified from microarray. In addition, the upregulation of CD1c expression on DCs from patients with allergies observed in the current study was supported by previous studies, including a study, which revealed that pattern recognition receptor expression profiles of CD1c+ DCs are more prone to interact with various allergens than those of other DC subtypes. Another study reported that CD1c+ DCs on human nasal mucosa induce CRR7-dependent migration and Th2 immune response after being stimulated by thymic stromal lymphopoietin [26,27]. Furthermore, it has been previously identified that Jagged 1 (JAG1) expression on DCs is responsible for Th2 differentiation by interacting with Notch on CD4+ T cells, and our result was consistent with the finding that JAG1 was upregulated in DCs from patients with allergies, with statistical significance unless the Log2FC value was below 1.5 [28].
IPA with an input of the top 100 DEGs revealed that genetic changes in DCs from patients with allergies might be related to pathways involving innate immunity, antigen presentation, and Th1/Th2 activation pathways. Enriched functional annotation and network analysis yielded major regulatory functions of the top 100 DEGs and their interactions at the cellular level. Functions associated with signal transduction, cell death/viability, and cellular migration are essential in the affected biological functions, and colony-stimulating factor 2 (CSF2) is associated with a majority of these functions. Among the 11 identified DEGs, AREG, FLNB, THBD, HTRA1 (HtrA serine peptidase 1), and TPST1 were related to the top functional annotations, suggesting that the products of these genes could play critical roles in the induction of allergic reactions. Furthermore, IPA upstream regulator analysis, which predicts the causal regulator of genetic expression changes, showed that IL-13 is an activated key regulator of target molecules, including C3AR1, CD1C, DUSP10, F13A1, HTRA1, PLCXD1, ST8SIA4, and TLR1 (z-score 2.828, overlap p-value < 0.0001) (data not shown). Belinghausen et al. reported that IL-13 might be a critical factor for inducing a Th2 immune response by increasing the phosphorylation of STAT6 in immature DCs rather than in mature DCs or T cells [29]. Our current findings are consistent with this finding, and we suggest candidates for possible intracellular processes downstream of IL-13 to induce allergic reactions.
The top network with the highest score yielded under IPA data was associated with cell death/survival, cancer and organismal injury/abnormality consisting of ERK1/2, Akt, and CSF2 as key nodes. Although cell death/survival might be a pathway that could be affected by experimental factors, we believed that this pathway could play a role in inducing Th2 inflammation in that the life span and maturation of dendritic cells are known to be two important factors in regulating adoptive immunity [30]. ERK1/2 is a well-known intracellular signaling molecule in DCs as it stabilizes the transcription factor c-Fos, thus promoting IL-10 production, suppressing IL-12 (p70) production, and finally inducing T cell response biased toward Th2 differentiation [31,32]. CSF2, also known as granulocyte-macrophage colony-stimulating factor (GM-CSF), is reported to be an epithelial cell-derived cytokine that induces DC maturation and increases the expression of co-stimulatory molecules, thus promoting Th2 inflammation [33,34]. In addition, Acciani et al. reported that the activation of AREG and epidermal growth factor receptor (EGFR), which leads to GM-CSF expression, is an essential signaling pathway in HDM-induced allergic airway epithelial cells [35]. AREG-EGFR signaling is suggested to play a key role in type 2 immune responses, as it is expressed on epithelial and immune cells, including CD4+ T cells and innate lymphoid cells type 2 [36]. AREG expression on DCs has been reported to trigger tumorigenesis in response to ATP and promote pulmonary fibrosis in CD11c+ bone marrow DCs; however, previous studies have not addressed its role in allergic responses [37,38]. From the network analysis yield from IPA data in the present study, the top network consisted of AREG showing interactions with CSF2 and ERK1/2, suggesting that intracellular downstream signals of AREG in DCs could also play a key role in inducing Th2 immune responses. Our validation study showed increased AREG mRNA expression in Der p 1-treated THP-1-derived DCs and higher protein expression in DCs harvested from patients with allergies compared to the control. These findings also support the IPA data-based analysis results.
The second most significant network with the functions of cell-to-cell signaling/interaction and cellular movements consisted of various receptors associated with antigen presentation and chemokines, which might interact with extracellular IL-4 and TNF family members. The majority of the subunits constituting the HLA-DR group (HLA-DR alpha, beta1, and beta5) and representing the MHC class II receptor (signal 1), tended to be upregulated in our analysis; however, this upregulation was not statistically significant. Additionally, chemokines (signal 3) in the network, such as CCR5, were downregulated with p-value < 0.05 in the DCs of patients with allergies. CCR5 expression on DCs and T cells has been reported to be responsible for inducing a Th1 immune response by secreting cytokines, such as IL-12, thus supporting our result [39,40]. TPST1, one of the 6 downregulated genes in our data, was included in this network. TPST1, which catalyzes post-transcriptional tyrosine sulfation between various membrane proteins, is reported to act on chemokine receptors such as CCR5 and CXCR4. These previous studies support our data, indicating that the downregulation of TPST1 might reduce the interaction of CCR5-CCL3, thus triggering an immune response against Th1 [12,41,42]. Moreover, FLNB, which regulates cellular migration by binding to β-integrin, was involved in this network [43,44]. Previous studies have reported the role of β–integrin in restricting DC maturation and migration to lymph nodes from tissues, supporting our result of the potential function of FLNB [45,46]. Furthermore, another study revealed that the downstream regulation of FLNB suppresses ERK1/2 signaling, suggesting another functional mechanism for FLNB downregulation in DCs from individuals with allergies [47]. From the validation experiment of mRNA expression after allergen stimulation in THP-1-derived DCs, the expression of TPST-1 and FLNB varied over stimulation time; nevertheless, the expression of both genes showed a decreasing tendency at 24 h compared to that before allergen stimulation. The western blotting results of these genes in our study were contrary to those of the transcriptome analysis, and we presumed that various factors, including translation efficiency and transcription-translation feedback, could be associated. Therefore, further research is required.
Notably, there are a few limitations of the present study. First, the number of samples in each group was small and showed heterogeneity in terms of age and sex. However, to characterize the allergic group and the controls, we thoroughly excluded samples with clinical histories that could affect the study’s outcome, such as having any other significant medical histories or receiving medication. The allergic and control groups showed a clear difference in objective allergy tests, the proportion of blood eosinophils, and nasal endoscopic findings. Furthermore, to overcome the analysis results based on unadjusted p-value via validation of the transcriptome analysis results, experimental studies using cell line-derived DCs and DCs harvested from patients were conducted. Second, the DCs analyzed in our study were not derived from respiratory tissues. We attempted to isolate DCs from the nasal mucosa and extract their RNA for whole transcriptome sequencing. However, because of the difficulty in extracting a sufficient amount of DCs for transcriptome sequencing from nasal tissue, we decided to perform the analysis in mDCs sorted from peripheral blood. Nevertheless, the strength of our experiments is that mDCs of PBMCs were sorted through flow cytometry using surface antibodies, and this is the first attempt to compare whole transcriptome profiles of DCs between patients with allergies and controls.

4. Materials and Methods

4.1. Subjects

Samples from a total of 10 East Asians who visited Korea University Anam Hospital, including patients with allergic diseases (n = 6) and healthy controls (n = 4), were analyzed (Table 4). The exclusion criteria were as follows: history of smoking, recurrent or recent exacerbation of upper respiratory infection within 6 weeks preceding the study, ongoing medical treatment, and history of other systemic medical conditions. All participants underwent either a skin-prick test (SPT) or multiple-allergen simultaneous test (MAST) for inhalant allergens as well as serum total IgE determination by enzyme-linked immunosorbent assay (CAP system). All patients with allergies showed positive MAST or SPT results to HDM, positive reactions to the allergen provocation test, and nasal eosinophilia and active upper airway allergic symptoms consisting of watery rhinorrhea, itching, nasal obstruction, or sneezing. Participants with allergies who had a history of receiving allergic medication, including antihistamines, leukotriene modifiers, and corticosteroids, within 8 weeks preceding the study were excluded. Participants with negative results in the MAST and no history of allergic symptoms were classified into the regular group.
Physical examination was performed, including anterior rhinoscopy and rigid endoscopy, and inferior turbinate findings were graded according to the guidelines provided by Camacho et al. [48]. Additionally, the proportion of eosinophils and mononuclear cells among inflammatory cells in peripheral blood was calculated. The clinical data of the participants are summarized in Table 4. Comparisons of clinical characteristics between the healthy controls and patients with allergies were analyzed using Mann–Whitney U-tests. A p-value < 0.05 was considered statistically significant.
Before obtaining peripheral blood samples, the protocol was approved by the Institutional Review Board for Human Studies at Korea University Hospital, and all enrolled participants provided written informed consent.

4.2. Preparation of Human PBMCs (hPBMCs) from Blood and Flow Cytometry

Peripheral blood from all participants was collected into a mononuclear cell preparation tube containing sodium heparin and Ficoll Hypaque solution (BD Vacutainer CPT, BD Biosciences, Franklin Lakes, NJ, USA). hPBMCs were enriched between plasma and the Ficoll layer by centrifugation at 1500× g for 20 min at 25 °C. The hPBMCs were carefully transferred to a new tube and washed with 10 volumes of PBS, followed by centrifugation at 300× g for 15 min at 25 °C. Precipitated hPBMCs were rewashed with PBS by centrifugation at 300× g for 15 min at 25 °C and resuspended in staining buffer (2% FBS and 1 mM EDTA in PBS).
hPBMCs were washed with 2 mL of staining buffer by centrifugation at 300× g for 5 min. Washed hPBMCs were counted, and then the Fc receptors were blocked with 1 μg of pre-immune IgG per 1 × 105 cells. The Fc receptor-blocked cells were stained with a multi-color flow cytometry kit (R&D Systems, Inc., Minneapolis, MN, USA) for human myeloid DCs. The antibodies used for staining myeloid DCs were BDCA-3/CD141-PE Mouse IgG1 (Clone 501733), BDCA-1/CD1c-APC Goat IgG, CD11c-CFS Mouse IgG1 (Clone ICRF 3.9), and CD16/Fcγ RIII-PerCP Mouse IgG2A (Clone 245536). The cells were incubated with each antibody or corresponding isotype control antibody for 45 min at room temperature in the dark. After washing with staining buffer, stained cells were sorted on a BD FACSAria II (BD Biosciences) and analyzed on a FACSDiva (BD Biosciences). By measuring forward and side light scattering, viable cells were detected, and cell debris was excluded from the analysis. DC subsets stained with anti-BDCA-3 and CD11c antibodies were gated and sorted. The sorted DCs were lysed with QIAzol lysis reagent (QIAGEN, Hilden, Germany) and subjected to mRNA sequencing.

4.3. RNA Sequencing, Differentially Expression Analysis and Pathway Analysis

The preparation of whole transcriptome libraries and sequencing were conducted by Macrogen Inc. (Seoul, Korea). The RNAs from sorted DCs were isolated and underwent quality control. Quant-IT RiboGreen (Invitrogen, Carlsbad, CA, USA) was used to measure the total RNA concentration and the TapeStation RNA screentape (Agilent Technologies, Santa Clara, CA, USA) was used to verify the RNA. The RNA library construction consisted of only high-quality RNA samples (1 μg) with RNA integrity number (RIN) over 7.0. After removing ribosomal RNAs, RNA fragments underwent reverse transcription using random primers, and 100 nt paired-end sequencing was conducted by Illumina HiSeq2500. The library was quantified using qPCR according to the qPCR Quantification Protocol Guide (KAPA BIOSYSTEMTS, Wilmington, MA, USA) and verified using Agilent Technologies 2100 Bioanalyzer (Agilent Technologies).
mRNA-seq libraries were prepared using a paired-end mRNA sequencing sample preparation kit (TruSeq RNA Access Library Kit) following the manufacturer’s protocols. A flow cell containing unique clusters, using a unique bridge amplification reaction, was loaded into the Illumina HiSeq 2500 for automated extension and imaging cycles. Low-quality and adapter sequences from produced paired-end reads were trimmed using Trimmomatic (version 0.36) [49]. To extract RNA sequence variants, we employed the Genome Analysis Toolkit (GATK) [50] Best Practices workflow for SNP and InDel calling on RNA sequencing data. Briefly, the STAR 2-pass method was used to align raw data to the human reference genome (hg19) using STAR aligner [51]. SAM files produced in the above step were processed by Picard for adding read group information, sorting, marking duplicates, and indexing. Variants were called and filtered using the GATK HaplotypeCaller and VariantFiltration tools, respectively. All variants were annotated using ANNOVAR [52]. Enriched function annotation and pathway analysis were performed using IPA software (Ingenuity Systems Inc., Redwood City, CA, USA).
Differential expression analysis of individual genes was carried out using the “new tuxedo” protocol [53]. Trimmed RNA sequence paired-end reads were aligned to the human reference genome (hg19) using HISAT2 (version 2.0.5) [54]. StringTie (version 1.3.3b) was used to produce reconstructions of genes and estimate expression levels [55]. Ballgown was used to calculate differential gene expression from RNA sequencing data [56]. In detail, fragments per kilo base of transcript per million base pairs sequenced was used to estimate the level of gene expression, and the p-value for differential expression was extracted using a parametric F-test comparing nested linear models. The log2 ratio of fold change (log2FC) of the gene expression value between patients with allergies and controls was calculated using the Ballgown “stattest” function. Differential gene expression results were visualized using a volcano plot and heap map using R packages “calibrate (version 1.7.7)” and “ggplot2 (version 3.3.2)”, respectively.

4.4. THP-1-Derived Dendritic Cell Culture

To investigate changes in the expression of selected genes from transcriptome sequencing analysis in allergen-treated dendritic cells, a human monocytic cell line, THP-1, was used. THP-1 cells were incubated with 10% fetal bovine serum (FBS), penicillin/streptomycin, 50 µM β-mercaptoethanol (Sigma-Aldrich, Saint Louis, MO, USA), and 2 mM glutamine in RPMI 1640 (WELGENE, Gyeongsan, Korea) at 37 °C. To prepare DCs, THP-1 cells were incubated with IL-4 (100 ng/mL) (Pepro Tech, Rocky Hill, NJ, USA) and GM-CSF (100 ng/mL) (Pepro Tech). At day 2, the cells were replenished with IL-4 and GM-CSF, followed by incubation for 2 days to differentiate into immature DCs. At day 6, Der p 1 (2 μg/mL, Indoor Biotechnologies, Charlottesville, VA, USA) was added to the media for the maturation of DCs, and cells were cultured for up to 6 days. DCs were used to validate selected gene expression changes after Der p 1 stimulation for 2, 4, 8, 12, and 24 h.

4.5. Real-Time PCR

THP-1-derived dendritic cells were plated on 12-well plates with g/mL Der p 1 (Indoor biotechnology, Charlottesville, VA, USA). RNA was extracted from the control and Der p 1 treated cells after 2, 4, 8, 12, and 24 h of incubation. Detached THP-1- derived dendritic cells were lysed with QIAzol (Qiagen, Valencia, CA, USA) to extract total RNA. Sequentially, the lysates were treated with chloroform, ethanol, and isopropanol for the purification of mRNA. cDNA synthesis from extracted RNA was conducted using the cDNA synthesis master mix (GenDEPOT, Katy, TX, USA).
Gene expression in THP-1 derived dendritic cells was measured by quantitative real-time PCR. Amplification and quantification of cDNA were performed using SYBR Green master mix (Qiagen, Valencia, CA, USA) with specific primers (Table S3). PCRs were conducted using a real-time thermal cycler system (TP850) (Takara, Shiga, Japan) with 50 cycles of a two-step reaction: (1) denaturation at 95 °C for 15 s, and (2) annealing extension at 60 °C for 45 s. Data analysis was conducted using the ΔCt method. To determine the difference in gene expression among Der p 1 treated THP-1 derived dendritic cells, a one-way analysis of variance was used.

4.6. Western Blotting

Validation of selected genes (AREG, TPST-1, and FLNB) from RNA sequencing analysis at the protein level was conducted by western blotting. Protein extracts were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to nitrocellulose membranes. Membranes were incubated with anti-AREG (Santa Cruz Biotechnology, Dallas, TX, USA), -TPST-1 (Abcam, Cambridge, UK), and -FLNB antibodies (SantaCruz Biotechnology, Dallas, TX, USA) for the target molecule anti-GAPDH antibody (SantaCruz Biotechnology, Dallas, TX, USA) as a reference. Protein visualization was conducted using a ChemiDoc imaging system (Bio-Rad, Hercules, CA, USA).

5. Conclusions

In the present study, whole transcriptome profiling of immature DCs from patients with allergies and normal controls revealed 11 DEGs between the two groups and their expected functions regulating cell viability/activation, cell migration, and cytokine production through intracellular signaling or DC-T cell interaction. Among the several novel transcriptomes identified, AREG was verified through validation, and is thus proposed as a possible therapeutic target for allergic diseases. Nonetheless, further experiments with functional studies are needed to clarify the effects of this gene on the induction of Th2 immune responses.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/21/22/8640/s1.

Author Contributions

Conceptualization, T.H.K. and K.L.; methodology, B.K. and J.W.Y.; software, M.-R.H.; validation, K.L., J.W.Y. and B.K.; formal analysis, M.-R.H.; investigation, J.W.Y; resources, T.H.K; data curation, K.L.; writing—original draft preparation, K.L., and M.-R.H.; writing—review and editing, T.H.K. and B.K.; visualization, J.W.Y.; supervision, and T.H.K.; project administration, T.H.K.; funding acquisition, T.H.K. and M.-R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Basic Science Research Program, National Research Foundation of Korea, funded by the Ministry of Science and Technology and the Ministry of Science, ICT & Future Planning (2017R1A2B2003575, NRF-2020R1A2C1006398), the Ministry of Science and ICT (2020R1C1C1012288), Korea, under the ICT Creative Consilience program(IITP-2020-0018190011001) supervised by the IITP (Institute for Information & Communications Technology Planning & Evaluation, and the Korea Health Technology R&D Project (HI17C0387), Korea Health Industry Development Institute (KHIDI), and the Ministry of Health & Welfare. This research was also supported by a Korea University grant and a grant from Korea University Medical Center and Anam Hospital, Seoul, Korea.

Acknowledgments

We would like to acknowledge the support of Byung-Joo Ham (Department of Psychiatry, Korea University Anam Hospital, Korea University College of Medicine, Seoul, Korea) for enriched functional annotation and network identification using ingenuity pathway analysis.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DCDendritic cell
mDCMyeloid dendritic cell
pDCPlasmacytoid dendritic cell
PCRPolymerase chain reaction
NGSNext-generation sequencing
DEGDifferentially expressed genes
IPAIngenuity pathway analysis
Log2FCLog2-fold change
HDMHouse dust mite
PBMCPeripheral blood monocyte
hPBMCHuman peripheral blood monocyte
SPTSkin prick test
MASTMultiple-allergen simultaneous test

References

  1. Froidure, A.; Shen, C.; Pilette, C. Dendritic cells revisited in human allergic rhinitis and asthma. Allergy 2016, 71, 137–148. [Google Scholar] [CrossRef] [PubMed]
  2. Lambrecht, B.N.; Hammad, H. Biology of Lung Dendritic Cells at the Origin of Asthma. Immunity 2009, 31, 412–424. [Google Scholar] [CrossRef] [PubMed]
  3. Lambrecht, B.N.; De Veerman, M.; Coyle, A.J.; Gutierrez-Ramos, J.-C.; Thielemans, K.; Pauwels, R.A. Myeloid dendritic cells induce Th2 responses to inhaled antigen, leading to eosinophilic airway inflammation. J. Clin. Investig. 2000, 106, 551–559. [Google Scholar] [CrossRef] [PubMed]
  4. Rissoan, M.C.; Soumelis, V.; Kadowaki, N.; Grouard, G.; Briere, F.; de Waal Malefyt, R.; Liu, Y.J. Reciprocal Control of T Helper Cell and Dendritic Cell Differentiation. Science 1999, 283, 1183–1186. [Google Scholar] [CrossRef]
  5. Shortman, K.; Liu, Y.-J. Mouse and human dendritic cell subtypes. Nat. Rev. Immunol. 2002, 2, 151–161. [Google Scholar] [CrossRef]
  6. Rescigno, M.; Martino, M.; Sutherland, C.L.; Gold, M.R.; Ricciardi-Castagnoli, P.; Williams, N.S.; Moore, T.A.; Schatzle, J.D.; Puzanov, I.J.; Sivakumar, P.; et al. Dendritic Cell Survival and Maturation Are Regulated by Different Signaling Pathways. J. Exp. Med. 1998, 188, 2175–2180. [Google Scholar] [CrossRef]
  7. Na, H.; Cho, M.; Chung, Y. Regulation of Th2 Cell Immunity by Dendritic Cells. Immune Netw. 2016, 16, 1–12. [Google Scholar] [CrossRef] [Green Version]
  8. Macdonald, A.S.; Maizels, R.M. Alarming dendritic cells for Th2 induction. J. Exp. Med. 2008, 205, 13–17. [Google Scholar] [CrossRef] [Green Version]
  9. Dardalhon, V.; Korn, T.; Kuchroo, V.K.; Anderson, A.C. Role of Th1 and Th17 cells in organ-specific autoimmunity. J. Autoimmun. 2008, 31, 252–256. [Google Scholar] [CrossRef] [Green Version]
  10. Endo, Y.; Hirahara, K.; Yagi, R.; Tumes, D.J.; Nakayama, T. Pathogenic memory type Th2 cells in allergic inflammation. Trends Immunol. 2014, 35, 69–78. [Google Scholar] [CrossRef]
  11. Maggi, E. The TH1/TH2 paradigm in allergy. Immunotechnology 1998, 3, 233–244. [Google Scholar] [CrossRef]
  12. Liu, G.; Bi, Y.; Xue, L.; Zhang, Y.; Yang, H.; Chen, X.; Lu, Y.; Zhang, Z.; Liu, H.; Wang, X.; et al. Dendritic cell SIRT1-HIF1alpha axis programs the differentiation of CD4+ T cells through IL-12 and TGF-beta1. Proc. Natl. Acad. Sci. USA 2015, 112, E957–E965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Huang, G.; Wang, Y.; Vogel, P.; Kanneganti, T.-D.; Otsu, K.; Chi, H. Signaling via the kinase p38α programs dendritic cells to drive TH17 differentiation and autoimmune inflammation. Nat. Immunol. 2012, 13, 152–161. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, G.; Wang, Y.; Shi, L.Z.; Kanneganti, T.-D.; Chi, H. Signaling by the Phosphatase MKP-1 in Dendritic Cells Imprints Distinct Effector and Regulatory T Cell Fates. Immunity 2011, 35, 45–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Suryawanshi, A.; Manoharan, I.; Hong, Y.; Swafford, D.; Majumdar, T.; Taketo, M.M.; Manicassamy, B.; Koni, P.A.; Thangaraju, M.; Sun, Z.; et al. Canonical wnt signaling in dendritic cells regulates Th1/Th17 responses and suppresses autoimmune neuroinflammation. J. Immunol. 2015, 194, 3295–3304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  17. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef]
  18. Mahendra, A.; Yang, X.; Abnouf, S.; Adolacion, J.R.; Park, D.; Soomro, S.; Roszik, J.; Coarfa, C.; Romain, G.; Wanzeck, K.; et al. Beyond Autoantibodies: Biologic Roles of Human Autoreactive B Cells in Rheumatoid Arthritis Revealed by RNA-Sequencing. Arthritis Rheumatol. 2019, 71, 529–541. [Google Scholar] [CrossRef]
  19. Sumida, T.; Lincoln, M.R.; Ukeje, C.M.; Rodriguez, D.M.; Akazawa, H.; Noda, T.; Naito, A.T.; Komuro, I.; Dominguez-Villar, M.; Hafler, D.A. Activated beta-catenin in Foxp3(+) regulatory T cells links inflammatory environments to autoimmunity. Nat. Immunol. 2018, 19, 1391–1402. [Google Scholar] [CrossRef]
  20. Lyons, Y.A.; Wu, S.Y.; Overwijk, W.W.; Baggerly, K.A.; Sood, A.K. Immune cell profiling in cancer: Molecular approaches to cell-specific identification. NPJ Precis. Oncol. 2017, 1, 1–8. [Google Scholar] [CrossRef] [Green Version]
  21. Yick, C.Y.; Zwinderman, A.H.; Kunst, P.W.; Grunberg, K.; Mauad, T.; Dijkhuis, A.; Bel, E.H.; Baas, F.; Lutter, R.; Sterk, P.J. Transcriptome sequencing (RNA-Seq) of human endobronchial biopsies: Asthmaversuscontrols. Eur. Respir. J. 2013, 42, 662–670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Suarez-Farinas, M.; Ungar, B.; Da Rosa, J.C.; Ewald, D.A.; Rozenblit, M.; Gonzalez, J.; Xu, H.; Zheng, X.; Peng, X.; Estrada, Y.D.; et al. RNA sequencing atopic dermatitis transcriptome profiling provides insights into novel disease mechanisms with potential therapeutic implications. J. Allergy Clin. Immunol. 2015, 135, 1218–1227. [Google Scholar] [CrossRef] [PubMed]
  23. Seumois, G.; Zapardiel-Gonzalo, J.; White, B.; Singh, D.; Schulten, V.; Dillon, M.; Hinz, D.; Broide, D.H.; Sette, A.; Peters, B.; et al. Transcriptional Profiling of Th2 Cells Identifies Pathogenic Features Associated with Asthma. J. Immunol. 2016, 197, 655–664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sheu, C.-C.; Tsai, M.-J.; Chen, F.-W.; Chang, K.-F.; Chang, W.-A.; Chong, I.-W.; Kuo, P.-L.; Hsu, Y.-L. Identification of novel genetic regulations associated with airway epithelial homeostasis using next-generation sequencing data and bioinformatics approaches. Oncotarget 2017, 8, 82674–82688. [Google Scholar] [CrossRef] [Green Version]
  25. Yerkovich, S.T.; Roponen, M.; Smith, M.E.; McKenna, K.; Bosco, A.; Subrata, L.S.; Mamessier, E.; Wikström, M.E.; Le Souef, P.; Sly, P.D.; et al. Allergen-enhanced thrombomodulin (blood dendritic cell antigen 3, CD141) expression on dendritic cells is associated with a TH2-skewed immune response. J. Allergy Clin. Immunol. 2009, 123, 209–216.e4. [Google Scholar] [CrossRef]
  26. Lundberg, K.; Rydnert, F.; Greiff, L.; Lindstedt, M. Human blood dendritic cell subsets exhibit discriminative pattern recognition receptor profiles. Immunology 2014, 142, 279–288. [Google Scholar] [CrossRef] [Green Version]
  27. Melum, G.R.; Farkas, L.; Scheel, C.; Van Dieren, B.; Gran, E.; Liu, Y.-J.; Johansen, F.-E.; Jahnsen, F.L.; Baekkevold, E.S. A thymic stromal lymphopoietin–responsive dendritic cell subset mediates allergic responses in the upper airway mucosa. J. Allergy Clin. Immunol. 2014, 134, 613–621.e7. [Google Scholar] [CrossRef] [Green Version]
  28. Okamoto, M.; Matsuda, H.; Joetham, A.; Lucas, J.J.; Domenico, J.; Yasutomo, K.; Takeda, K.; Gelfand, E.W. Jagged1 on Dendritic Cells and Notch on CD4+ T Cells Initiate Lung Allergic Responsiveness by Inducing IL-4 Production. J. Immunol. 2009, 183, 2995–3003. [Google Scholar] [CrossRef] [Green Version]
  29. Bellinghausen, I.; Brand, P.; Bottcher, I.; Klostermann, B.; Knop, J.; Saloga, J. Production of interleukin-13 by human dendritic cells after stimulation with protein allergens is a key factor for induction of T helper 2 cytokines and is associated with activation of signal transducer and activator of transcription-6. Immunology 2003, 108, 167–176. [Google Scholar] [CrossRef]
  30. Maney, N.J.; Reynolds, G.; Krippner-Heidenreich, A.; Hilkens, C.M.U. Dendritic cell maturation and survival are differentially regulated by TNFR1 and TNFR2. J. Immunol. 2014, 193, 4914–4923. [Google Scholar] [CrossRef] [Green Version]
  31. Dillon, S.; Agrawal, A.; Van Dyke, T.; Landreth, G.; McCauley, L.; Koh, A.; Maliszewski, C.; Akira, S.; Pulendran, B. A Toll-Like Receptor 2 Ligand Stimulates Th2 Responses In Vivo, via Induction of Extracellular Signal-Regulated Kinase Mitogen-Activated Protein Kinase and c-Fos in Dendritic Cells. J. Immunol. 2004, 172, 4733–4743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Nakahara, T.; Moroi, Y.; Uchi, H.; Furue, M. Differential role of MAPK signaling in human dendritic cell maturation and Th1/Th2 engagement. J. Dermatol. Sci. 2006, 42, 1–11. [Google Scholar] [CrossRef] [PubMed]
  33. Gill, M.A. The role of dendritic cells in asthma. J. Allergy Clin. Immunol. 2012, 129, 889–901. [Google Scholar] [CrossRef] [PubMed]
  34. Deckers, J.; De Bosscher, K.; Lambrecht, B.N.; Hammad, H. Interplay between barrier epithelial cells and dendritic cells in allergic sensitization through the lung and the skin. Immunol. Rev. 2017, 278, 131–144. [Google Scholar] [CrossRef] [PubMed]
  35. Acciani, T.H.; Suzuki, T.; Trapnell, B.C.; Le Cras, T.D. Epidermal growth factor receptor signalling regulates granulocyte-macrophage colony-stimulating factor production by airway epithelial cells and established allergic airway disease. Clin. Exp. Allergy 2016, 46, 317–328. [Google Scholar] [CrossRef] [Green Version]
  36. Zaiss, D.M.W.; Gause, W.C.; Osborne, L.C.; Artis, D. Emerging Functions of Amphiregulin in Orchestrating Immunity, Inflammation, and Tissue Repair. Immunity 2015, 42, 216–226. [Google Scholar] [CrossRef] [Green Version]
  37. Ding, L.; Liu, T.; Wu, Z.; Hu, B.; Nakashima, T.; Ullenbruch, M.; De Los Santos, F.G.; Phan, S.H. Bone Marrow CD11c+ Cell-Derived Amphiregulin Promotes Pulmonary Fibrosis. J. Immunol. 2016, 197, 303–312. [Google Scholar] [CrossRef] [Green Version]
  38. Bles, N.; Di Pietrantonio, L.; Boeynaems, J.-M.; Communi, D. ATP confers tumorigenic properties to dendritic cells by inducing amphiregulin secretion. Blood 2010, 116, 3219–3226. [Google Scholar] [CrossRef]
  39. Luther, S.A.; Cyster, J.G. Chemokines as regulators of T cell differentiation. Nat. Immunol. 2001, 2, 102–107. [Google Scholar] [CrossRef]
  40. Aliberti, J.; e Sousa, C.R.; Schito, M.; Hieny, S.; Wells, T.; Huffnagle, G.B.; Sher, A. CCR5 provides a signal for microbial induced production of IL-12 by CD8α+ dendritic cells. Nat. Immunol. 2000, 1, 83–87. [Google Scholar] [CrossRef]
  41. Ludeman, J.P.; Stone, M.J. The structural role of receptor tyrosine sulfation in chemokine recognition. Br. J. Pharmacol. 2014, 171, 1167–1179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kehoe, J.W.; Bertozzi, C.R. Tyrosine sulfation: A modulator of extracellular protein–protein interactions. Chem. Biol. 2000, 7, R57–R61. [Google Scholar] [CrossRef] [Green Version]
  43. Calderwood, D.A.; Huttenlocher, A.; Kiosses, W.B.; Rose, D.M.; Woodside, D.G.; Schwartz, M.A.; Ginsberg, M.H. Increased filamin binding to beta-integrin cytoplasmic domains inhibits cell migration. Nat. Cell Biol. 2001, 3, 1060–1068. [Google Scholar] [CrossRef] [PubMed]
  44. Baldassarre, M.; Razinia, Z.; Burande, C.F.; Lamsoul, I.; Lutz, P.G.; Calderwood, D.A. Filamins Regulate Cell Spreading and Initiation of Cell Migration. PLoS ONE 2009, 4, e7830. [Google Scholar] [CrossRef] [Green Version]
  45. Schittenhelm, L.; Hilkens, C.M.; Morrison, V.L. Beta2 Integrins as Regulators of Dendritic Cell, Monocyte, and Macrophage Function. Front Immunol. 2017, 8, 1866. [Google Scholar] [CrossRef]
  46. Morrison, V.L.; James, M.J.; Grzes, K.; Cook, P.C.; Glass, D.G.; Savinko, T.; Lek, H.S.; Gawden-Bone, C.; Watts, C.; Millington, O.R.; et al. Loss of beta2-integrin-mediated cytoskeletal linkage reprogrammes dendritic cells to a mature migratory phenotype. Nat. Commun. 2014, 5, 5359. [Google Scholar] [CrossRef] [Green Version]
  47. Bandaru, S.; Zhou, A.X.; Rouhi, P.; Zhang, Y.; Bergo, M.O.; Cao, Y.; Akyurek, L.M. Targeting filamin B induces tumor growth and metastasis via enhanced activity of matrix metalloproteinase-9 and secretion of VEGF-A. Oncogenesis 2014, 3, e119. [Google Scholar] [CrossRef]
  48. Camacho, M.; Zaghi, S.; Certal, V.; Abdullatif, J.; Means, C.; Acevedo, J.; Liu, S.; Brietzke, S.E.; Kushida, C.A.; Capasso, R. Inferior Turbinate classification system, grades 1 to 4: Development and validation study. Laryngoscope 2015, 125, 296–302. [Google Scholar] [CrossRef]
  49. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  50. Van Der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  51. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  52. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  53. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. [Google Scholar] [CrossRef] [PubMed]
  54. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Frazee, A.C.; Pertea, G.; Jaffe, A.E.; Langmead, B.; Salzberg, S.L.; Leek, J.T. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 2015, 33, 243–246. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 08640 g001 550
Figure 1. Differentially expressed genes (DEGs) identified in dendritic cells (DCs) from patients with allergies versus controls: (A) Heat map constructed with five upregulated genes and six downregulated genes with p-values < 0.05 and |Log2FC| ≥ 1.5. Color key indicates the degree of differential gene expression in either direction (red color represents upregulation; green color represents downregulation); (B) Volcano plot was constructed by plotting the log2FC of individual genes on x-axis and the negative logarithm of their p-value to base 10 on y-axis. Colored circles refer to 1658 DEGs with p-values < 0.05, and among those, orange circles represent genes with |Log2FC| ≥ 1.5, while green circles refer to genes with |Log2FC| ≥ 2.0.
Figure 1. Differentially expressed genes (DEGs) identified in dendritic cells (DCs) from patients with allergies versus controls: (A) Heat map constructed with five upregulated genes and six downregulated genes with p-values < 0.05 and |Log2FC| ≥ 1.5. Color key indicates the degree of differential gene expression in either direction (red color represents upregulation; green color represents downregulation); (B) Volcano plot was constructed by plotting the log2FC of individual genes on x-axis and the negative logarithm of their p-value to base 10 on y-axis. Colored circles refer to 1658 DEGs with p-values < 0.05, and among those, orange circles represent genes with |Log2FC| ≥ 1.5, while green circles refer to genes with |Log2FC| ≥ 2.0.
Ijms 21 08640 g001
Ijms 21 08640 g002 550
Figure 2. Top two networks constructed by ingenuity pathway analysis (IPA) showing interactions between the top 100 differentially expressed genes (DEGs) in the dendritic cells (DCs) from patients with allergies. (A) The network with the highest score of 36, comprising amphiregulin (AREG), extracellular-signal-regulated kinase (ERK1/2), colony-stimulating factor (CSF), protein kinase B (Akt), was associated with cell death/survival, cancer, and organismal injury/abnormalities. (B) The second-highest scoring network with a score of 18, comprising filamin B (FLNB), tyrosylprotein sulfotransferase1 (TPST1), and C-X-C chemokine receptor type 4 (CXCR4), was associated with cell-to-cell signaling/interaction, cellular movement, and cardiovascular system development/function. The downregulated genes in allergic patients are colored in green, and the upregulated genes are colored in red. The gradient of the colored shape indicates DEG strength.
Figure 2. Top two networks constructed by ingenuity pathway analysis (IPA) showing interactions between the top 100 differentially expressed genes (DEGs) in the dendritic cells (DCs) from patients with allergies. (A) The network with the highest score of 36, comprising amphiregulin (AREG), extracellular-signal-regulated kinase (ERK1/2), colony-stimulating factor (CSF), protein kinase B (Akt), was associated with cell death/survival, cancer, and organismal injury/abnormalities. (B) The second-highest scoring network with a score of 18, comprising filamin B (FLNB), tyrosylprotein sulfotransferase1 (TPST1), and C-X-C chemokine receptor type 4 (CXCR4), was associated with cell-to-cell signaling/interaction, cellular movement, and cardiovascular system development/function. The downregulated genes in allergic patients are colored in green, and the upregulated genes are colored in red. The gradient of the colored shape indicates DEG strength.
Ijms 21 08640 g002
Ijms 21 08640 g003 550
Figure 3. Validation of selected differential expressed genes from transcriptome sequencing data. (A–C) Quantitative real-time polymerase chain reaction of Der p 1-stimulated THP-1 derived dendritic cells. (D) Western blotting of human myeloid dendritic cells sorted from peripheral blood monocytes of healthy controls and patients with allergies. * p < 0.05.
Figure 3. Validation of selected differential expressed genes from transcriptome sequencing data. (A–C) Quantitative real-time polymerase chain reaction of Der p 1-stimulated THP-1 derived dendritic cells. (D) Western blotting of human myeloid dendritic cells sorted from peripheral blood monocytes of healthy controls and patients with allergies. * p < 0.05.
Ijms 21 08640 g003
Table
Table 1. Eleven upregulated and downregulated genes with p-values < 0.05 and |Log2FC| ≥ 1.5 in myeloid dendritic cells from patients with allergies compared with controls, as determined by RNA sequencing data.
Table 1. Eleven upregulated and downregulated genes with p-values < 0.05 and |Log2FC| ≥ 1.5 in myeloid dendritic cells from patients with allergies compared with controls, as determined by RNA sequencing data.
Gene SymbolLog2FCp-ValueGene Name
Upregulated in patients with allergies
AREG2.10.0160Amphiregulin
THBD1.70.0463Thrombomodulin
CD1C1.50.0396Cluster of differentiation 1
RBM121.50.0353RNA-binding protein 12
SLC7A5P11.50.0349Solute carrier family 7 member 5 pseudogene 1
Downregulated in patients with allergies
TPST1−2.70.0258Tyrosylprotein Sulfotransferase 1
C3AR1−2.20.0030Complement Component 3a Receptor 1
HTRA1−1.90.0448High-Temperature Requirement A Serine Peptidase 1
FLNB−1.60.0198Filamin B
FUCA1−1.50.0127Alpha-L-Fucosidase 1
CCL3L3−1.50.0467Chemokine ligand 3-like 3
Table
Table 2. Enriched functional annotation of top 100 differentially expressed genes in dendritic cells from patients with allergies.
Table 2. Enriched functional annotation of top 100 differentially expressed genes in dendritic cells from patients with allergies.
Functional Annotationp-ValueRelated Gene SymbolNumber of Molecules
Signal transduction0.0001AREG, BRD8, CCR5, CD79B, CLEC5A, CREM, FLNB, MX1, NFATC1SIK1/SIK1B, TLR1, TLR612
Cell death0.0004AREG, C3AR1, CCR5, CREM, CSF2, GADD45A, ITGAV, MAP4K1, MX1, NFATC1, NOTCH3, THBD, SERPINE1, SIK1/SIK1B, TFRC, TLR1, TLR617
Activation of cells0.0005AREG, CD1C, CSF2, CTNND1, THBD, TLR1, TLR67
Antigen presentation0.0016CD1C, CSF22
Cell viability0.0028CSF2, MX12
Cell spreading0.0052FLNB, ITGAV, SERPINE13
Tyrosine sulfation of protein0.0069TPST11
Secretion of molecule0.0070ABCG1, C3AR1, CSF2, STX114
Migration of cells0.0110AREG, C8orf44−SGK3/SGK3, CCR5, CSF2, FLNB, HTRA1, ITGAV, NFATC1, SERPINE1, THBD10
Anchoring of cytoskeleton0.0172FLNB1
Differentiation of monocyte-derived dendritic cells0.0273CSF21
Synthesis of DNA0.0277AREG, CSF2, NFATC13
Table
Table 3. Canonical ingenuity pathway analysis of the top 100 differentially expressed genes in dendritic cells from patients with allergies.
Table 3. Canonical ingenuity pathway analysis of the top 100 differentially expressed genes in dendritic cells from patients with allergies.
Ingenuity Canonical Pathways−Log(p-Value)RatioMolecules
Communication between Innate and Adaptive Immune Cells4.900.0575CD79B, CCL3L3, TLR6, TLR1, CSF2
Altered T Cell and B Cell Signaling in Rheumatoid Arthritis3.710.0482CD79B, TLR6, TLR1, CSF2
Coagulation System3.620.0857F13A1, SERPINE1, THBD
SAPK/JNK Signaling3.260.0367GADD45A, DUSP10, MAP4K1, NFATC1
Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses2.960.0305TLR6, TLR1, CSF2, C3AR1
Extrinsic Prothrombin Activation Pathway2.860.1250F13A1, THBD
TREM1 Signaling2.740.0429TLR6, TLR1, CSF2
Caveolar-mediated Endocytosis Signaling2.720.0423ITGAE, FLNB, ITGAV
VDR/RXR Activation2.620.0390GADD45A, THBD, CSF2
Interferon Signaling2.160.0556MX1, ISG15
Intrinsic Prothrombin Activation Pathway2.050.0488F13A1, THBD
Phagosome Formation2.050.0244TLR6, TLR1, C3AR1
Th1 Pathway1.990.0231CCR5, NOTCH3, NFATC1
S-adenosyl-L-methionine Biosynthesis1.990.3330MAT2A
Epithelial Adherens Junction Signaling1.840.0204NOTCH3, TUBB2A, CTNND1
Remodeling of Epithelial Adherens Junctions1.660.0303TUBB2A, CTNND1
Role of Macrophages, Fibroblasts and Endothelial Cells in Rheumatoid Arthritis1.650.0129TLR6, TLR1, CSF2, NFATC1
Th1 and Th2 Activation Pathway1.610.0167CCR5, NOTCH3, NFATC1
Toll-like Receptor Signaling1.560.0267TLR6, TLR1
Glucocorticoid Receptor Signaling1.510.0117POLR2J2/POLR2J3, SERPINE1, CSF2, NFATC1
Pathogenesis of Multiple Sclerosis1.510.1110CCR5
Leucine Degradation I1.510.1110BCAT1
Hematopoiesis from Multipotent Stem Cells1.390.0833CSF2
Fcγ Receptor-mediated Phagocytosis in Macrophages and Monocytes1.390.0215PLD4, CSF2
NF-κB Activation by Viruses1.390.0215CCR5, ITGAV
TGF-β Signaling1.390.0215MAP4K1, SERPINE1
Systemic Lupus Erythematosus Signaling1.380.0135CD79B, CREM, NFATC1
Choline Biosynthesis III1.360.0769PLD4
Protein Kinase A Signaling1.360.0104FLNB, DUSP10, CREM, NFATC1
Isoleucine Degradation I1.330.0714BCAT1
Phospholipase C Signaling1.310.0127PLD4, CD79B, NFATC1
Table
Table 4. Characteristics of patients with allergies and healthy controls participating in the study.
Table 4. Characteristics of patients with allergies and healthy controls participating in the study.
CharacteristicsPatients with AllergiesHealthy Controls
No. of patients64
No. of female/male1:53:1
Age(y), mean range33.4(18–52)55.3(47–66)
Serum total IgE (IU/mL) 1537.3 ± 267.535.7 ± 12.6
Serum eosinophil (%) 14.4 ± 1.11.3 ± 0.6
Serum mononuclear cells (%)7.2 ± 1.46.7 ± 0.9
Endoscopic turbinate grade 13.7 ± 0.21
Duration of upper airway hyperresponsiveness(y)12.3 ± 1.9-
Specific IgEs with positive results (%)
D. pteronyssinus1000
D. farinae800
House dust800
Pollen500
Animal500
Other500
1 The values in patients with allergies were significantly higher than those in healthy controls.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, K.; Han, M.-R.; Yeon, J.W.; Kim, B.; Kim, T.H. Whole Transcriptome Analysis of Myeloid Dendritic Cells Reveals Distinct Genetic Regulation in Patients with Allergies. Int. J. Mol. Sci. 2020, 21, 8640. https://doi.org/10.3390/ijms21228640

AMA Style

Lee K, Han M-R, Yeon JW, Kim B, Kim TH. Whole Transcriptome Analysis of Myeloid Dendritic Cells Reveals Distinct Genetic Regulation in Patients with Allergies. International Journal of Molecular Sciences. 2020; 21(22):8640. https://doi.org/10.3390/ijms21228640

Chicago/Turabian Style

Lee, Kijeong, Mi-Ryung Han, Ji Woo Yeon, Byoungjae Kim, and Tae Hoon Kim. 2020. "Whole Transcriptome Analysis of Myeloid Dendritic Cells Reveals Distinct Genetic Regulation in Patients with Allergies" International Journal of Molecular Sciences 21, no. 22: 8640. https://doi.org/10.3390/ijms21228640

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

Article Metrics

Back to TopTop