ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 09 February 2021

Sec. Molecular and Cellular Pathology

Volume 9 - 2021 | https://doi.org/10.3389/fcell.2021.636477

Long Non-coding RNAs and mRNAs Expression Profiles of Monocyte-Derived Dendritic Cells From PBMCs in AR

  • YZ

    Yumei Zhou 1†

  • XC

    Xuemei Chen 1†

  • YZ

    Yanfei Zheng 1†

  • RS

    Rongmin Shen 1

  • SS

    Shuxian Sun 1

  • FY

    Fei Yang 1

  • JM

    Jiayu Min 1

  • LB

    Lei Bao 1

  • YZ

    Yan Zhang 1

  • XZ

    Xiaoshan Zhao 2

  • JW

    Ji Wang 1*

  • QW

    Qi Wang 1*

  • 1. National Institute of TCM Constitution and Preventive Medicine, School of Chinese Medicine, Beijing University of Chinese Medicine, Beijing, China

  • 2. School of Traditional Chinese Medicine, Southern Medical University, Guangzhou, China

Abstract

Objective:

The objective of this study is to explore the long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) expression profiles of monocyte-derived dendritic cells (DCs) obtained from peripheral blood mononuclear cells (PBMCs). DCs are known to play a major role in the regulating function of allergic rhinitis (AR).

Methods:

PBMCs were separately isolated from the human peripheral blood of patients with AR and normal person (NP). The mixed lymphocyte reaction (MLR) assay was used to evaluate the function of DCs. Flow cytometry was used to determine the immune regulatory function of immature DCs (imDCs) and mature DCs (mDCs). lncRNAs and mRNAs in the NP group (DCs isolated from NP) and the test group (DCs isolated from patients with AR) were identified via chip technology and bioinformatic analyses. Moreover, bioinformatic analyses were employed to identify the related biological functions of monocyte-derived DCs and construct the functional networks of lncRNAs and mRNAs that are differentially expressed (DE) in imDCs and mDCs.

Results:

MLR was significantly higher in the mDCs group than that in the imDCs group. CD14 was highly expressed in imDCs, whereas HLA-DR, CD80, and CD86 were highly expressed in mDCs (p < 0.001). We identified 962 DE lncRNAs and 308 DE mRNAs in the imDCs of NP and patients with AR. Additionally, there were 601 DE lncRNAs and 168 DE mRNAs in the mDCs in the NP and test groups. Quantitative RT-qPCR was used to study the significant fold changes of lncRNAs and mRNAs. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found 16 significant regulated pathways in imDCs and 10 significant regulated pathways in mDCs, including the phagosome, cell adhesion signaling pathway, and inflammatory mediator regulation of TRP channels pathway.

Conclusion:

Our research studied the lncRNA and mRNA expression profiles of monocyte-derived DCs and demonstrated the functional networks that are involved in monocyte-derived DCs-mediated regulation in AR. These results provided possible molecular mechanisms of monocyte-derived DCs in the immunoregulating function and laid the foundation for the molecular therapeutic targets of AR.

Background

Allergic rhinitis (AR) is a very common allergic disease that affects 10–40% of the global population (Bousquet et al., 2008). Its remarkable prevalence and relapses put an extensive burden on its patients and the society. Furthermore, AR negatively impacts the quality of life of patients with AR. There had been a marked increase in the prevalence of AR during the past years (Wang et al., 2016). AR turns into asthma if it is not treated in time, and the adequate treatment of AR can alleviate the severity of asthma (Leynaert et al., 2004). Researchers have illustrated that mast cell infiltration, lymphocytes imbalance, and goblet cell hyperplasia are involved in the pathogenesis of AR (Ouyang et al., 2010; Poggi et al., 2012). AR, which is a type I allergic disorder that is mediated by IgE humoral immune response, is accompanied by an influx of eosinophils and T helper 2 cells that secrete pro-inflammatory cytokines, namely, IL-4, IL-5, and IL-13 (Wilson et al., 2005). Abnormal innate and adaptive immune responses play a major role in the pathogenesis of AR.

Dendritic cells (DCs), which are the most important antigen-presenting cells (APCs) that send signals to the T cells, mainly participate in the pathogenesis of many diseases with immunoregulatory mechanisms, such as AR. DCs link the innate and adaptive immune responses. The peripheral blood mononuclear cells (PBMCs) have a round nucleus (Delves, 2016). PBMCs include lymphocytes, monocytes, and DCs. In humans, the frequencies of these DCs vary among individuals. PBMCs are divided into various functional subtypes with respect to the specific cytokine expression profiles, surface markers, and the transcription factors. Phenotypic and functional assessments of PBMC research lay the foundation of the human immune system research; hence, the knowledge that population is represented in the peripheral blood and how they act with other immune cells is essential. Additionally, the results from human PBMC studies (Schiekofer et al., 2003; Tacconi et al., 2004; Chang et al., 2014) cannot be neglected. Therefore, it is important to know the progression of AR along with its expression profiles in PBMCs, especially DCs.

Long non-coding RNAs (lncRNAs), over 200 nt in length, is a type of RNA that does not a protein coding function (Ulitsky and Bartel, 2013). These RNAs have been regarded as indispensable epigenetic regulators and are probably involved in the cell’s biological behaviors (Kopp and Mendell, 2018). For example, they are involved in regulating the homeostasis of the immune system (Wang et al., 2014; Du et al., 2017). However, it is critical to find out whether lncRNA can immunoregulate DC in the progression of AR.

The combination of lncRNA–messenger RNA (mRNA) expression profiles and functional networks is adopted to analyze the DC-mediated regulation functions. These results improve our understanding of lncRNAs in the immunoregulatory function of monocyte-derived DCs and indicate the potential targets for the curative treatment of AR.

Materials and Methods

Subjects

Patients with AR visited doctors in the outpatient service in the Guo Yi Tang of Beijing University of Chinese Medicine. In this study, there were 24 subjects: 12 males and 12 females. They were divided into two groups: the AR group (patients with AR, 12 subjects: five males and seven females) and NP group (normal persons, 12 subjects: four males and eight females). With support/approval from the Ethics Committee of Beijing University of Chinese Medicine, this study was conducted while adhering to the principles of the Declaration of Helsinki. Patients in the AR group were positive for skin puncture test, including pollen, food, dust mites, paint, or molds as well as in specific IgE. Two weeks before study recruitment, these patients with AR received no topical or systemic corticosteroid therapy. We chose the study participants with no history of smoking or other immune system disorders, such as rheumatoid arthritis, systemic lupus erythematosus, and scleroderma.

Isolation of PBMCs and Generation of DCs

The whole blood samples obtained from the two groups were stored in vacuum tubes with heparin, and PBMCs were isolated from these samples by lymphocyte separation solution (Tianjin Haoyang Biological Manufacture Co., Ltd.). Mononuclear cells were seeded in 12-well plates with the RPMI 1640 medium that contains 10% heat-inactivated fetal calf serum (FCS, GIBCO, Germany) and 2 mM of L-glutamine (R10 medium, Sigma, St. Louis, MO, United States). After incubation at 37°C for 2 h, the non-adherent cells were removed and the adherent cells were cultured within the medium containing 100 ng/ml of rhGM-CSF and 100 ng/ml of rhIL-4 (R&D Systems, Minneapolis, MN, United States). On the sixth day, 50 ng/ml of TNF-α (R&D Systems, Minneapolis, MN, United States) was added into the samples; the method had the same protocol in the study of Andreia et al. (2005). Immature DCs (imDCs) were collected on the fifth day, and the mature DCs (mDCs) were collected on the seventh day.

Mixed Lymphocyte Reaction (MLR)

After being treated with 25 μg/mL of mitomycin at 37°C for 30 min, the DCs were placed at the concentrations of 2 × 108 cells per well at a quantity of 200 μl and incubated with non-adherent PBMCs obtained from the same healthy people. These samples were later stimulated by lymphocytes in the concentration proportions of 1:10, 1:50, and 1:100. Thereafter, the samples were mixed and incubated with non-adherent PBMC from the same healthy persons at the same concentrations in triplicate. The cells were treated with 10% fetal bovine serum. Then, 10 μl of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium bromide (MTT) solution (5 mg/ml, medium dilution, Sigma–Aldrich Chemical Co., St. Louis, MO, United States) was added to each well, and the cells were incubated for 72 h in the incubator. Then, 150 μl of DMSO was added followed by the addition of enzymes after 4 h. The absorbance was detected using a spectrophotometer at 570 nm.

DC Surface Marker Expression Analysis in imDCs and mDCs

The CD14 (PerCP-Cy 5.5, BD Biosciences, United States), HLA-DR (APC, BD Biosciences, United States), and isotype mouse IgG2a-PE (PE, BioLegend, United States) were added in the samples of imDCs. Moreover, CD86-APC, CD80-PE, and isotype mouse IgG1–FITC (BioLegend, United States) were added in the samples of mDCs. The cells were then suspended with precooled PBS, counted under a microscope, and centrifuged at 1000 g for 5 min. Data were acquired using a FACSCalibur cytometer (BD Biosciences, United States) and the ratios of CD14+, HLA-DR+, CD80+, and CD86+ DCs were determined.

RNA Extraction, Labeling, Chip Hybridization, and Scanning

The RNA extraction, labeling, chip hybridization, and scanning were all finished following the use of Agilent Human lncRNA–mRNA profiling chip (4∗180K, Design ID: 062918). All the chip results were detailed according to the processes of software operations. The screening criteria were to increase or decrease the fold change value ≥ 2.0 and p-value < 0.05.

Gene Ontology (GO) and Pathway Enrichment Analysis

The GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were used to study the differentially expressed (DE) mRNAs. The results of the target gene analysis of lncRNA and the mRNA expression results on the chip need to be correlated so that the upregulated and downregulated DE genes can be investigated. The corrected p-value < 0.05 by calculating the FDR and FDR < 0.05 was selected as the threshold.

LncRNA–mRNA-Weighted Co-expression Network

The correlation of the lncRNA–mRNA expression in the imDCs and mDCs was calculated. Then, the relationship pairs of LncRNA and mRNA based on the abovementioned criteria (p-value < 0.05, FDR < 0.05) were screened. The co-expression network of lncRNAs and mRNAs was constructed, and the co-expression network of lncRNAs and mRNAs was then established.

Real-Time PCR Verification of DE Genes

Quantitative RT-qPCR was used to investigate the different expressions of imDCs’ and mDCs’ genes between the patients with AR and NP. The total RNA was extracted from imDCs and DCs according to the kit for cells. In all, 20 μg of total RNA was converted into cDNA by using oligo (dT) and reverse transcriptase (Thermo, United States) to analyze the qPCR results. The thermal cycler conditions were set as follows: amplificated at 95°C for 10 min, 95°C for 15 s, 60°C for 60 s, 40 cycles of denaturation (15 s, 94°C), 15 s at 95°C, and a combined process of annealing and extension (1 min, 60°C). Supplementary Table S1 shows the primers for these genes.

Statistical Analysis

The SPSS 22.0 software was used for statistical analysis in this study. Results were expressed as mean ± standard deviation. The relationship of lncRNAs and mRNAs was determined by the Spearman’s correlation coefficient. p-values < 0.05 were considered significant values for this study.

Results

Allogeneic T-Cell Proliferation Experiment

In allergenic mixed lymphocyte reaction (MLR), the levels of T-cell proliferation were increased with the proportion of T cells. DCs in the patients with AR have a stronger stimulation ability than NP as shown in Figures 1A–C. When the ratio of DC cells to T cells was 1:50 and 1:100, the difference was significant (p < 0.001).

FIGURE 1

Immunophenotype of DCs in Patients With AR and NP

We evaluated the percentages of CD14, HLA-DR, CD80, and CD86 in the DCs. The immunophenotypic characteristics of DCs in patients with AR were compared with that of NP, as shown in Figure 2. The CD14 concentration of imDCs was lower in patients with AR than that in the NP group. The mean percentage of CD14+ DCs in the patients with AR was 0.84 ± 0.25% (median: 0.78%), and it was significantly lower than the value in the NP group (p = 0.0008), where the mean percentage of these cells was 0.012 ± 0.013% (median: 0.017%) (Figures 2A,B).

FIGURE 2

The mean proportion of HLA-DR+ mDCs in the patients with AR was 44.25 ± 8.64% (median: 44.8%) and was higher (p = 0.0006) than the NP group, where the mean percentage of these cells was 95.3 ± 3.84% (median: 99.0%) (Figures 2C,D).

The mean proportion of CD80+ mDCs in the patients with AR was 72.15 ± 7.64% (median: 67.78%) and was higher (p = 0.0003) than the NP group, where the mean percentage of these cells was 95.3 ± 3.84% (median: 96.92%) (Figures 2E,F).

The mean proportion of CD86+ mDCs in the patients with AR was 62.15 ± 7.64% (median: 64.69%) and was higher (p = 0.006) than the control group, where the mean percentage of these cells was 92.3 ± 5.84% (median: 95.12%) (Figures 2G,H).

Identification of DE lncRNAs and mRNAs in imDCs

In total, 308 DE mRNAs, including 175 upregulated mRNAs and 133 downregulated mRNAs, were found in the imDCs of patients with AR and NP. A clustergram (Figure 3A) and volcano plots (Figure 3C) are used to depict DE mRNAs. Table 1 shows 67 mRNAs with the largest fold changes. The list contains several genes, including HLA-C, MARCO, KIR2DS3, ITGAV, CD36, and IFNB1. Additionally, 168 DE mRNAs, including 77 upregulated mRNAs and 91 downregulated mRNAs, were found in the mDCs of patients with AR. Figures 4A,C show the clustergram and volcano plots of the DE mRNAs, respectively. Table 2 shows the 10 mRNAs with the greatest variation, and several genes, such as HLA-B, F11R, HLA-DQB1b, HLA-DQB1, and PTAFR, are also shown.

FIGURE 3

TABLE 1

GeneNameGenbank AccessionFC (abs)Regulation
Fatty acid 2-hydroxylaseNM_0243062.7454393Up
Uncharacterized LOC645984AK0954362.289749Up
Shisa homolog 9 (Xenopus laevis)NM_0011452052.3817623Down
NK2 homeobox 1NM_0033172.2473493Down
Sodium channel, voltage-gated, type VII, alphaNM_0029763.4931335Down
Nucleolar and spindle associated protein 1NM_0163592.0980568Up
Ventral anterior homeobox 1NM_1991312.798632Down
Chromosome 20 open reading frame 132NM_2136312.831788Up
Programmed cell death 1 ligand 2NM_0252393.0216887Down
Solute carrier family 29 (nucleoside transporters), member 1NM_0010781772.2991307Up
Very low density lipoprotein receptorNM_0010180562.4687624Down
Tropomyosin 2 (beta)NM_2136742.175945Down
Uncharacterized LOC100131129AK1271842.1942973Up
Peptidyl arginine deiminase, type IINM_0073653.0724247Up
CD300 molecule-like family member fNM_1390182.5870113Up
Programmed cell death 1 ligand 2NM_0252392.9215207Down
CD36 molecule (thrombospondin receptor)NM_0010015472.167322Down
Suppression of tumorigenicity 5NM_0054182.673464Up
T-cell acute lymphocytic leukemia 1NM_0031892.3056405Down
Fc receptor-like BNM_0010029012.3118186Up

The characteristics of mRNAs with the largest fold change in imDCs.

FIGURE 4

TABLE 2

mRNAGenelncRNAGeneCorrelation coefficientP-value
NM_018208ETNK2NON-HSAT008926ETNK20.99706213.24775E–10
NM_005891ACAT2NON-HSAG045301NON-HSAG0453010.99531732.09182E–09
NM_182985TRIM69NR_104175.1LOC400799−0.99140542.36259E–08
BX538082GPR17FR351114FR3511140.99064623.31169E–08
NM_002977SCN9ANON-HSAG029733NON-HSAG0297330.98820978.33527E–08
NM_014485HPGDSNON-HSAT097445HPGDS0.98776579.6583E–08
AK131565LOC100132368ENST00000588609LINC00906-0040.98638381.47941E–07
NM_005005NDUFB9NON-HSAT101913RP11-1113N2.40.98606541.62207E–07
XM_001719518LOC100128869NON-HSAG042749NON-HSAG042749−0.9857261.78529E–07
NM_032772ZNF503ENST00000438293RP11-88H9.2-0030.98533181.98986E–07

The top 10 co-expression of mRNA and lncRNA in imDCs.

In total, 962 DE lncRNAs, including 434 upregulated and 528 downregulated lncRNAs, were found in patients with AR and NP. The clustergram (Figure 3B) and volcano plots (Figure 3D) show DE lncRNAs. In total, 601 lncRNAs, including 200 upregulated and 401 downregulated lncRNAs, were found in the DE mDCs of the patients with AR. The clustergram in Figure 4B and the volcano plots in Figure 4D show the DE lncRNAs. Tables 1, 3 show the top 10 lncRNAs of imDCs and mDCs with the largest fold changes. The pathways, including interferon-gamma-mediated signaling pathway, membrane repolarization, and peptide antigen binding, that contribute to the phagocytosis function in imDCs and antigen-presenting function of mDCs were also identified.

TABLE 3

GeneNameAccession no.Fold changeRegulation
Armadillo repeat containing 9NM_0251392.1053998Up
Cat eye syndrome chromosome region, candidate 2NM_0314132.1667209Down
Semenogelin IINM_0030082.2934558Up
Polycystic kidney disease 1 like 1NM_1382953.8540776Up
Rho guanine nucleotide exchange factor (GEF) 35NM_0010037022.662969Up
Sushi domain containing 4NM_0179823.1600573Down
ST6 beta-galactosamide alpha-2,6-sialyltranferase 1NM_1732162.504244UP
Chromosome 2 open reading frame 71NM_0010298832.1819167Up
Dehydrogenase/reductase (SDR family) member 3NM_0047532.2156718Up
Sterile alpha motif and leucine zipper containing kinase AZKNM_0166532.2057335UP
Chromosome 5 open reading frame 62NM_0329472.0991595Up
Transmembrane protein ENSP00000343375BC0313042.0236611Down
Family with sequence similarity 101, member BNM_1827052.506473Up
Olfactory receptor, family 2, subfamily T, member 8NM_0010055222.5608928Down
Sorbin and SH3 domain containing 1NM_0010349542.0617623Down
Transmembrane protease, serine 6BC0390822.2889757Down
DEAD (Asp-Glu-Ala-Asp) box polypeptide 6NM_0043972.0167336Up
spondin 2, extracellular matrix proteinNM_0124452.7114112Up
HLA-DBQ1NM_0012439612.154Down

The characteristics of mRNAs with the largest fold change in mDCs.

Interaction, Co-expression Network Analysis of DE mRNAs in Patients With AR

Figure 5A shows the interactions of proteins that were coded by DE mRNAs in imDCs. Additionally, Figure 5B shows the co-expression network between DE lncRNAs and mRNAs. TRIM69 has the maximum target mRNAs including 58 DE mRNAs, and SIDT2 has the maximum co-expressed lncRNAs. Table 2 shows the top 10 co-expression pairs in imDCs. lncRNAs exert their biological function as ceRNAs (Wilfried et al., 2016).

FIGURE 5

We identified 268 target genes after analyzing the possible DE lncRNAs target genes in imDCs. Figure 5C shows the target genes with a combined score of more than 0.9. HLA-C was the target gene of AC108142.1-005, and CD36 was the target gene of FR264384. Moreover, IFNB1 was the target gene of MIR3150B-210. The Venn diagram analysis showed that 95 mRNAs were coincided between the 166 DE mRNAs and the 95 DE lncRNA target genes (Figure 5D). The 95 DE lncRNAs were all included in the 166 DE mRNAs.

Figure 6A shows the interaction proteins that were coded by DE mRNAs in mDCs. In this network, HLA-B, HLA-DQB1, HLA-DQB2, PTAFR, and F11R are important genes that interact with many other DE mRNAs. Furthermore, Figure 6B shows the co-expression network of DE lncRNAs and mRNAs. TRIM77P has the maximum target numbers including 39 DE mRNAs, in which FAM153A and ZNF396 have the maximum co-expressed lncRNAs. Table 4 shows the top 10 co-expression pairs of mDCs.

FIGURE 6

TABLE 4

mRNAGenelncRNAGeneCorrelation coefficientP-value
NM_012445SPON2NON-HSAG037252NON-HSAG0372520.99474993.30E–09
NM_138782FCHO2NON-HSAT102095FCHO20.99453683.87E–09
AK311167LOC100132352TCONS_l2_00028804linc-ANKRD20A4-20.98329513.34E–07
NM_003162STRNENST00000562064CTD-2015G9.2-0010.98259293.93E–07
NM_001034954SORBS1TCONS_00010265linc-FASTKD3-10.97945637.60E–07
NM_001463FRZBNR_038973.1FAM170B-AS10.97778941.04E–06
NM_001979EPHX2NON-HSAT021480AP000445.1−0.97698111.19E–06
AB018295FAM153ATCONS_00002078linc-LPHN2-20.97479581.71E–06
NM_007048BTN3A1NON-HSAT142393RP11-429P3.20.97373922.02E–06
NM_182487OLFML2AFR052374FR0523740.97238322.46E–06
NM_001146162TRIM77PENST00000541885RP11-439H13.2-0010.97232542.48E–06

The top 10 co-expression of mRNA and lncRNA in mDCs.

Additionally, we also investigated the possible presence of DE lncRNA target genes in mDCs. In our study, 99 target genes were identified in these DE lncRNAs. Figure 6C shows the target genes with the combined score of more than 0.9. In the figure, HLA-B has three target genes of lincRNAs, namely, DHRS3, FCHO2, and linc-PRR5-1. Figure 6D shows that 42 mRNAs coincided between the 95 DE mRNAs and the 42 DE lncRNA target genes in the Venn diagram analysis. Moreover, the 42 DE lncRNAs were included in the 95 DE mRNAs. It also includes some known inflammatory-related molecules.

Validation of DE mRNA and lncRNA Expression Levels by RT-qPCR

RT-qPCR was used to evaluate DE mRNAs and lncRNAs to verify our RNA chip results. We randomly detected three lncRNAs and 10 mRNAs. MARCO, KIR2DS3, F11R, HLA-B, HLA-C, NON-HSAG046717, and NON-HSAT089067 were upregulated, whereas ITGAV, CD36, IFNB1, PTAFR, HLA-DQB1, NON-HSAT 059748, NON-HSAT024276, and NON-HSAT098958 were downregulated. The results of RT-qPCR were consistent with those of RNA chip results, hence confirming that our chip data were reliable (Figure 7).

FIGURE 7

Discussion

In our research, we used rhGM-CSF, rhIL-4, and TNF-α to induce the mDCs from PBMCs. The generation of human monocyte-derived DCs from whole blood was recognized by all the scientists (Wilfried et al., 2016). We investigated phenotypic and functional features of DE DCs in vitro from patients with AR and the NPs. Besides that, we used RT-qPCR to confirm these findings. Our KEGG pathway analysis (Tables 5, 6) indicates that interferon-gamma-mediated signaling pathway, membrane repolarization, and peptide antigen binding pathways contribute to the phagocytosis function in imDCs, and the antigen-presenting function in mDCs contributes to the immunoregulatory function of DCs in AR.

TABLE 5

PathwayCountP-valueCorrected P-valueGene
Vitamin digestion and absorption30.0033774920.355867536PLB1/SLC19A1/SCARB1
Butanoate metabolism30.0052683820.355867536AACS/ACADS/ACAT2
Hypertrophic cardiomyopathy50.0057647260.355867536TPM2/TPM1/MYH7/ITGAV/IGF1
Dilated cardiomyopathy50.007550820.355867536TPM2/TPM1/MYH7/ITGAV/IGF1
Aldosterone synthesis and secretion50.0082233470.355867536CACNA1I/KCNJ5/PRKCB/SCARB1/POMC
GnRH secretion40.0088966880.355867536CACNA1I/KCNJ5/PRKCB/PIK3R3
Aldosterone-regulated sodium reabsorption30.0115020850.390782867PRKCB/PIK3R3/IGF1
Phagosome60.0130260960.390782867C1R/ITGAV/SCARB1/MARCO/CD36/HLA-C
Glioma40.0153035730.408095288PRKCB/PIK3R3/CDK6/IGF1
Fat digestion and absorption30.0172912780.414990664SCARB1/CD36/ACAT2
Valine, leucine, and isoleucine degradation30.0231520580.479306938AACS/ACADS/ACAT2
Cholesterol metabolism30.0257622580.479306938LRPAP1/SCARB1/CD36
Natural killer cell-mediated cytotoxicity50.0259624590.479306938IFNB1/PRKCB/PIK3R3/KIR2DS3/HLA-C
Terpenoid backbone biosynthesis20.0316535910.542632982MVD/ACAT2
Inflammatory mediator regulation of TRP channels40.0388488830.574014469PRKCB/PIK3R3/IGF1/P2RY2
Parathyroid hormone synthesis, secretion, and action40.0464829010.574014469MMP14/RUNX2/PRKCB/PDE4D

Pathways KEGG analysis of imDCs.

TABLE 6

PathwayCountP-valueCorrected P-valueGene
Neuroactive ligand–receptor interaction60.0089160460.526895235PTAFR, POMC, OPRM1, MCHR2, GNRH1, MLNR
Allograft rejection20.0173160970.526895235HLA-DQB1, HLA-B
Graft-versus-host disease20.0209251870.526895235HLA-DQB1, HLA-B
Type I diabetes mellitus20.0218724840.526895235HLA-DQB1, HLA-B
Other types of O-glycan biosynthesis20.0258355630.526895235ST6GAL1, LFNG
Arginine and proline metabolism20.0289847930.526895235AMD1, ALDH7A1
Autoimmune thyroid disease20.0322797320.526895235HLA-DQB1, HLA-B
Glutathione metabolism20.0368906580.526895235GSTM3, G6PD
Viral myocarditis20.040505540.526895235HLA-DQB1, HLA-B
Cell adhesion molecules30.0446521390.526895235HLA-DQB1, HLA-B, F11R
Retinol metabolism20.0507619920.544537736DHRS3, ALDH1A1

Pathways of mDCs with the largest significant difference in KEGG analysis.

Dendritic cells play a central role in allergic inflammation (Froidure et al., 2016). The latest in vitro techniques allow the in vitro differentiation of DCs (Thurner et al., 1999). Their ability to induce the proliferation of T cells in the MLR assay is commonly used for evaluating their functions (Cao et al., 2004). In MLR experiment, the stimulation index of mDCs was significantly higher than that of the imDCs. In our study, the expression of HLA-DR, CD80, and CD86 in patients with AR is indeed upregulated than that in NP, which is a sign of improved mDCs antigen presentation in patients with AR. This result is in accordance with a previous study by KleinJan et al. (2006). In the study, CD14 in the NP group was higher than that in patients with AR. CD14 is the marker of monocytes whose expression decreases gradually during DC differentiation from monocytes. In fact, CD80 and CD86 are important co-stimulating factors that affect the proliferation of T lymphocytes in the DCs (Duperrier et al., 2000; Ebner et al., 2001; Andreia et al., 2005; Wilfried et al., 2016). DCs are the most efficient APCs. It can present the antigens to the T cells for stimulating the adaptive immune response.

More number of studies have focused on exploring the mRNA expressed in DCs, but none have been conducted to reveal which and how mRNA affects DCs in patients with AR. To investigate the mechanism of DCs’ functions in patients with AR, the mRNA expression profile was determined and bioinformatics analysis was performed in our study. In total, 308 mRNAs were identified in the analysis. Among these DE genes, HLA-C, ITGAV, MARCO, CD36, IFNB1, and KIR2DS3 were found to be most significantly upregulated in the imDCs’ network. HLA-C plays an important role in promoting differential DC maturation (Raj et al., 2011). ITGAV is the expression of DC-specific transmembrane protein. MARCO promotes TLR activation, which validates a major role of MARCO in mounting an inflammatory response (Haydn et al., 2014). ImDCs play an important role in the phagocytosis of apoptotic cells, in particular, CD36 (Albert et al., 1998). Additionally, the GO analysis demonstrated specific molecular functions, for example, binding of transforming growth factor-beta, signaling pattern recognition receptor activity, and pattern recognition receptor activity, thereby indicating the critical role of these cytokines in the immunoregulatory functions of imDCs. The KEGG analysis identified 17 signaling pathways of the DE mRNAs, wherein interferon-gamma-mediated signaling pathway, membrane repolarization, and peptide antigen binding were the pathways with the significant differences that contributed to the phagocytosis function of imDCs (Figures 3E,F). This result is consistent with that of a previous study (Lambrecht, 2001). Although research works have been conducted to determine the functions of imDCs in AR, the role of IFNB1 and KIR2DS3 in imDCs was not known in AR (Tables 7, 8).

TABLE 7

TermDomainCountP-valueCorrected P-value
Apoptotic cell clearanceBiological process54.46014E–054.46014E–05
VasculogenesisBiological process60.0001645070.000164507
Response to fatty acidBiological process60.0002637380.000263738
Phospholipid metabolic processBiological process130.0003191480.000319148
Glycerophospholipid metabolic processBiological process110.0003239520.000323952
Oligodendrocyte differentiationBiological process60.0004544510.000454451
Actin filament-based movementBiological process70.000474940.00047494
Glycerolipid metabolic processBiological process120.0007484760.000748476
Negative regulation of tumor necrosis factor productionBiological process50.0007692090.000769209
Carbohydrate derivative transportBiological process50.0007692090.000769209
Synaptic vesicle membraneCellular component50.0026492070.002649207
Exocytic vesicle membraneCellular component50.0026492070.002649207
Integral component of synaptic vesicle membraneCellular component30.0038053050.003805305
Presynaptic cytosolCellular component20.0085538710.008553871
PresynapseCellular component110.0089233160.008923316
Cell leading edgeCellular component100.0097984180.009798418
Integral component of organelle membraneCellular component60.0098957440.009895744
Exocytic vesicleCellular component60.0098957440.009895744
Melanosome membraneCellular component20.009908250.00990825
ChitosomeCellular component20.009908250.00990825
Transforming growth factor beta bindingMolecular function47.69207E–057.69207E–05
1-Phosphatidylinositol bindingMolecular function30.0004521980.000452198
Lipoprotein particle receptor activityMolecular function30.0004521980.000452198
Signaling pattern recognition receptor activityMolecular function30.0008242250.000824225
Pattern recognition receptor activityMolecular function30.0009809960.000980996
Amyloid-beta bindingMolecular function50.0011677640.001167764
Lipase activityMolecular function50.0052149710.005214971
Intronic transcription regulatory region sequence-specific DNA bindingMolecular function20.0052360340.005236034
Intronic transcription regulatory region DNA bindingMolecular function20.0052360340.005236034
Lipoprotein particle bindingMolecular function30.0063134990.006313499

GO analysis of DE mRNA in imDCs.

TABLE 8

TermDomainCountP-valueCorrected P-value
Interferon-gamma-mediated signaling pathwayBiological process50.0001051450.000105145
Regulation of membrane repolarizationBiological process30.0008976650.000897665
SA node cell to atrial cardiac muscle cell communicationBiological process20.0011431060.001143106
Cellular response to interferon-gammaBiological process50.0011524830.001152483
Atrial cardiac muscle cell membrane repolarizationBiological process20.0013924640.001392464
Cerebral cortex cell migrationBiological process30.0014909030.001490903
Telencephalon developmentBiological process60.0016816840.001681684
Response to interferon-gammaBiological process50.0019812970.001981297
Membrane repolarizationBiological process30.0021579710.002157971
Regulation of macrophage cytokine productionBiological process20.0022809060.002280906
MHC protein complexCellular component30.0003143840.000314384
Integral component of lumenal side of endoplasmic reticulum membraneCellular component30.0004917640.000491764
Lumenal side of endoplasmic reticulum membraneCellular component30.0004917640.000491764
Side of membraneCellular component80.0006898060.000689806
Clathrin-coated endocytic vesicle membraneCellular component30.0012741690.001274169
Plasma membrane protein complexCellular component80.0029122820.002912282
ER to Golgi transport vesicle membraneCellular component30.0031928270.003192827
MHC class II receptor activityCellular component20.0032428540.003242854
Alpha-actinin bindingCellular component30.0038958230.003895823
Actinin bindingCellular component30.0060416890.006041689
Protein phosphatase 2A bindingMolecular function30.0003143840.000314384
Lipopolysaccharide bindingMolecular function30.0004917640.000491764
Peptide antigen bindingMolecular function30.0004917640.000491764
Oxidoreductase activity, acting on the aldehyde or oxo group of donors, NAD or NADP as acceptorMolecular function80.0006898060.000689806
Oxidoreductase activity, acting on the aldehyde or oxo group of donorsMolecular function30.0012741690.001274169
Protease bindingMolecular function80.0029122820.002912282
Growth factor receptor bindingMolecular function30.0031928270.003192827
MHC class II receptor activityMolecular function20.0032428540.003242854
Alpha-actinin bindingMolecular function30.0038958230.003895823
Actinin bindingMolecular function30.0060416890.006041689

GO analysis of DE mRNA in mDCs.

The imDCs migrate to the lymphoid organs that will be matured in the future. They present captured Ag to the naïve T cells (Banchereau and Steinman, 1998). Hence, the imDCs and mDCs had different functions in the presented Ag. In our study, we focused on the immature and mature stages of DCs. In the mDC’ mRNA analysis, 168 mRNAs were identified. Among these DE genes, HLA-B, HLA-DQB1, HLA-DQB2, PTAFR, and F11R were the most significantly upregulated genes in mDCs. HLA-B is the major histocompatibility complex (class I) antigens that present the processed antigens. HLA-DQB1 is the major histocompatibility complex (class II) antigens that have been identified as useful biomarkers of candidacy for effective allergy immunotherapy in patients with AR (Yanming et al., 2019). This result indicates that these cytokines affect the antigen-presenting process in the mDCs. We can determine from the KEGG analysis that there are 12 signaling pathways related to DE mRNAs. The cell adhesion molecules were the useful pathways that contribute to the antigen-presenting function in mDCs (Figures 4E,F), which is consistent with a previous study (Roche and Furuta, 2015). Although many genes have been identified to play an important role in mDCs, but the mechanisms by which HLA-DQB2, PTAFR, and F11R affect mDC in AR are unknown (Tables 9, 10).

TABLE 9

lncRNA_AccessionFC (abs)RegulationChromosomeStrandStartEndClassSize
TCONS_l2_000304386.1069527DownchrX−55714615644346lncRNA72885
NON-HSAT0169345.6181483Downchr10−127823936127843874lncRNA19938
TCONS_l2_000012745.147873Upchr1−6545088065451399lncRNA519
NON-HSAT0169334.9528494Downchr10−127779304127798357lncRNA19053
NON-HSAT0597484.615665Downchr18+6681706566832387lncRNA15322
NON-HSAG0539334.34199UpchrX−24840702527190lncRNA43120
NON-HSAG0297334.2480536Downchr2−167054881167055243lncRNA362
NON-HSAT0939334.0003886Downchr3+188985384189038493lncRNA53109
NON-HSAG0369573.8607297Downchr3+188985385189038493lncRNA53108
NON-HSAG0558553.8563106DownchrY−2103438721239448lncRNA205061

lncRNAs in imDCs with the largest fold change.

TABLE 10

lncRNA_AccessionFC (abs)RegulationChromosomeStrandStartEndClassSize
NON-HSAT0169334.767738Downchr10−127779304127798357lincRNA19053
NON-HSAT0052463.4715514Downchr1−113068497113084597lincRNA16100
NR_038346.13.311445Downchr7+7908227279100524non-coding RNA18252
NON-HSAG0389663.027513Downchr4−141364352141419531lincRNA55179
NON-HSAG0372523.011491Upchr4−11651711202750lincRNA37579
NON-HSAG0087002.9691923Upchr11−6555652265562174lincRNA5652
ENST000004508472.9003682Downchr1−248647546248648785antisense1239
ENST000005418852.8448925Downchr12−6490094664927418lincRNA26472
TCONS_000076882.8411803Upchr4+185427281185436808lincRNA9527

The characteristics of lncRNAs with the largest fold change in mDCs.

Long non-coding RNA is an important component in the mRNA expression profiles (Kopp and Mendell, 2018). Several studies have shown that the differentiation of DCs is closely related to lncRNAs (Pin et al., 2014; Majid et al., 2020). It is confirmed from our studies that lncRNA may play a key role in the development of DC. Our data showed that 172 lncRNAs of imDCs and 104 lncRNAs of mDCs were significantly expressed in patients with AR as compared to that in NP by at least twofold changes. This result helps in studying the AR-related global transcriptome.

Long non-coding RNAs do not have a protein coding function, but they can be used as a new modulator, such as cis- or trans-gene regulating expression, demethylation-promoting effect, and mRNA-processing control were the major mechanisms (Sone et al., 2007; Wahlestedt, 2013; Qiao et al., 2016). To analyze the functions of lncRNAs, mRNA–lncRNA was used to create a co-expression profile to predict the potential functions of the DE lncRNAs of patients with AR. We found that the DE lncRNAs were all included in the results of mRNA–lncRNA chip results.

We used RT-qPCR to validate the mRNA–lncRNA chip results with only one randomly selected transcript, and found that the RT-qPCR results are consistent with the chip results in better understanding the role of lncRNA in the pathogenesis of DC-mediated AR.

In summary, our study proved that several mRNAs and lncRNAs of DC affect a certain process of AR pathogenesis by regulating the target genes. This result points out the direction for future studies on determining and explaining the functions and mechanisms of AR-related mRNA and lncRNA, and provides new therapeutic targets for patients with AR.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author/s.

Ethics statement

The studies involving human participants were reviewed and approved by the medical and experimental animal ethics committee of Beijing University of Traditional Chinese Medicine. The patients/participants provided their written informed consent to participate in this study.

Author contributions

YZ, XC, and YZ carried out the experiments. JW, QW, YZ, XC, and YZ designed the study and edited the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (No. 81973715) and Beijing Natural Science Foundation of China (No. 7202110).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.636477/full#supplementary-material

References

  • 1

    AlbertM. L.PearceS. F.FranciscoL. M.SauterB.RoyP.SilversteinR. L.et al (1998). Immature dendritic cells phagocytose apoptotic cells via alphavbeta5 and CD36, and cross-present antigens to cytotoxic T lymphocytes.J Exp Med.1881359–1368. 10.1084/jem.188.7.1359

  • 2

    AndreiaR. N.LuisF. C. E.AnselmoL. B.LeiteK. R. M.BuzaidA. C.Câmara-LopesL. H.et al (2005). Dendritic cells derived from metastatic cancer patients vaccinated with allogeneic dendritic cell- autologous tumor cell hybrids express more CD86 and induce higher levels of interferon-gamma in mixed lymphocyte reactions.Cancer Immunol. Immunother.5461–66. 10.1007/s00262-004-0550-8

  • 3

    BanchereauJ.SteinmanR. M. (1998). Dendritic cells and the control of immunity.Nature392245–252. 10.1038/32588

  • 4

    BousquetJ.KhaltaevN.CruzAA.DenburgJ.FokkensW. J.TogiasA.et al (2008). Allergic rhinitis and its impact on asthma (ARIA) 2008 update (in collaboration with the World Health Organization, GA (2) LEN and Aller Gen).Allergy638–160.

  • 5

    CaoM.-D.ChenZ.-D.XingY. (2004). Gamma irradiation of human dendritic cells influences proliferation and cytokine profile of T cells in autologous mixed lym- phocyte reaction.Cell Biol. Int.28223–228. 10.1016/j.cellbi.2003.12.006

  • 6

    ChangC. C.LinC. C.HsiehW. L.LaiH. W.TsaiC. H.ChengY. W. (2014). MicroRNA expression profiling in PBMCs: a potential diagnostic biomarker of chronic hepatitis C.Dis. Markers2014:157.

  • 7

    DelvesP. J. (2016). Roitt’s essential immunology, 11th edn.Trans. R. Soc. Trop. Med. Hyg.96108–108.

  • 8

    DuM.YuanL.TanX.HuangD.WangX.ZhengZ.et al (2017). The LPS-inducible lncRNA Mirt2 is a negative regulator of inflammation.Nat. Commun.8:2049.

  • 9

    DuperrierK.EljaafariA.Dezutter-DambuyantC.BardinC.JacquetC.YonedaK.et al (2000). Distinct subsets of dendritic cells resembling dermal DCs can be generated in vitro from monocytes, in the presence of different serum supplements.J. Immunol. Methods238119–131. 10.1016/s0022-1759(00)00147-2

  • 10

    EbnerS.RatzingerG.KrosbacherB.SchmuthM.WeissA.ReiderD.et al (2001). Production of IL-12 by human monocyte-derived dendritic cells is optimal when the stimulus is given at the onset of maturation, and is further enhanced by IL-4.J. Immunol.166633–641. 10.4049/jimmunol.166.1.633

  • 11

    FroidureA.ShenC.PiletteC. (2016). Dendritic cells revisited in human allergic rhinitis and asthma.Allergy71137–148. 10.1111/all.12770

  • 12

    HaydnT. K.DunnL. K.GhoshS.NechamaM.KobzikL.ArredouaniM. S. (2014). The scavenger receptor MARCO modulates TLR- induced responses in dendritic cells.PLoS One9:e104148. 10.1371/journal.pone.0104148

  • 13

    KleinJanA.WillartM.vanRijtL. S.BraunstahlG. J.LemanK.JungS.et al (2006). An essential role for dendritic cells in human and expe- rimental allergic rhinitis.J. Allergy Clin. Immunol.11851117–1125. 10.1016/j.jaci.2006.05.030

  • 14

    KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs.Cell172393–407. 10.1016/j.cell.2018.01.011

  • 15

    LambrechtB. N. (2001). Allergen uptake and presentation by dendritic cells.Curr. Opin. Allergy Clin. Immunol.151–59. 10.1097/01.all.0000010985.57414.74

  • 16

    LeynaertB.NeukirchC.KonyS.GuénégouA.BousquetJ.AubierM.et al (2004). Association between asthma and rhinitis according to atopic sensitization in a population-based study.J. Allergy Clin. Immunol.11386–93. 10.1016/j.jaci.2003.10.010

  • 17

    MajidP. K.RakhshiN.AleaghaM. S. E.AbdariM.AlikhahA.SafarianG.et al (2020). Differential expression of STAT3 gene and its regulatory long non-coding RNAs, namely lnc-DC and THRIL, in two eastern Iranian ethnicities with multiple sclerosis.Neurol Sci.41561–568. 10.1007/s10072-019-04092-y

  • 18

    OuyangY.MiyataM.HatsushikaK.OhnumaY.KatohR.OgawaH.et al (2010). TGF-beta signaling may play a role in the development of goblet cell hyperplasia in a mouse model of allergic rhinitis.Allergol. Int.59313–319. 10.2332/allergolint.10-sc-0172

  • 19

    PinW.XueY.HanY.LinL.WuC.XuS.et al (2014). The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation.Science344310–313. 10.1126/science.1251456

  • 20

    PoggiA.CanevaliP.ContatoreM.CiprandiG. (2012). Higher frequencies of CD161+ circulating t lymphocytes in allergic rhinitis patients compared to healthy donors.Int. Arch. Allergy Immunol.158151–156. 10.1159/000330903

  • 21

    QiaoJ.YaoH.XiaY.ChuP.LiM.WuY.et al (2016). Long non-coding RNAs expression profiles in hepatocytes of mice after hematopoietic stem cell transplantation.IUBMB Life68232–241. 10.1002/iub.1479

  • 22

    RajH.KubalC.MooreJ.NeilD.CookM.BallS.et al (2011). KIR and HLA-C interactions promote differential dendritic cell maturation and is a major determinant of graft failure following kidney transplantation.PLoS One6:e23631. 10.1371/journal.pone.0023631

  • 23

    RocheP. A.FurutaK. (2015). The ins and outs of mhc class-II-mediated antigen processing and presentation.Nat. Rev. Immunol.15203–216. 10.1038/nri3818

  • 24

    SchiekoferS.AndrassyM.ChenJ.RudofskyG.SchneiderJ.WendtT.et al (2003). Acute hyperglycemia causes intracellular formation of CML and activation of ras, p42/44 MAPK, and nuclear factor kappaB in PBMCs.Diabetes52621–633. 10.2337/diabetes.52.3.621

  • 25

    SoneM.HayashiT.TaruiH.AgataK.TakeichiM.NakagawaS. (2007). The mRNA-like noncoding RNA Gomafu constitutes a novel nuclear domain in a subset of neurons.J. Cell Sci.1202498–2506. 10.1242/jcs.009357

  • 26

    TacconiS.PerriR.BalestrieriE.GrelliS.BernardiniS.AnnichiaricoR.et al (2004). Increased caspase activation in peripheral blood mononuclear cells of patients with Alzheimer’s disease.Exp. Neurol.190254–262. 10.1016/j.expneurol.2004.07.009

  • 27

    ThurnerB.RoderC.CieckmannD.HeuerM.KruseM.GlaserA.et al (1999). Generation of large numbers of fully mature and stable dendritic cells from leukapheresis products for clinical application.J. Immunol. Methods2231–15. 10.1016/s0022-1759(98)00208-7

  • 28

    UlitskyI.BartelD. P. (2013). lincRNAs: genomics, evolution, and mechanisms.Cell15426–46. 10.1016/j.cell.2013.06.020

  • 29

    WahlestedtC. (2013). Targeting long non-coding RNA to therapeutically upregulate gene expression.Nat. Rev. Drug Discov.12433–446. 10.1038/nrd4018

  • 30

    WangP.XueY.HanY.LinL.WuC.XuS.et al (2014). The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation.Science344310–313.

  • 31

    WangX. D.ZhengM.LouH. F.WangC. S.ZhangY.BoM. Y.et al (2016). An increased prevalence of self-reported allergic rhinitis in major Chinese cities from 2005 to 2011.Allergy711170–1180. 10.1111/all.12874

  • 32

    WilfriedP.Lass-FlörlC.WilflingsederD. (2016). Generation of human monocyte-derived dendritic cells from whole blood.J. Vis. Exp.118:54968.

  • 33

    WilsonM. S.TaylorM. D.BalicA.FinneyC. A.LambJ. R.MaizelsR. M. (2005). Suppression of allergic airway inflammation by helminth-induced regulatory t cells.J. Exp. Med.2021199–1212. 10.1084/jem.20042572

  • 34

    YanmingZ.ZhaoY.ZhangY.ZhangL. (2019). HLA-II genes are associated with outcomes of specific immunotherapy for allergic rhinitis.Int. Forum. Allergy Rhinol.91311–1317. 10.1002/alr.22384

Summary

Keywords

monocyte-derived dendritic cells, allergic rhinitis, long non-coding RNA, immunoregulation, mRNA

Citation

Zhou Y, Chen X, Zheng Y, Shen R, Sun S, Yang F, Min J, Bao L, Zhang Y, Zhao X, Wang J and Wang Q (2021) Long Non-coding RNAs and mRNAs Expression Profiles of Monocyte-Derived Dendritic Cells From PBMCs in AR. Front. Cell Dev. Biol. 9:636477. doi: 10.3389/fcell.2021.636477

Received

01 December 2020

Accepted

19 January 2021

Published

09 February 2021

Volume

9 - 2021

Edited by

Shanchun Guo, Xavier University of Louisiana, United States

Reviewed by

Peng Wang, Sun Yat-sen University, China; Mamunur Rashid, Columbia University, United States

Updates

Copyright

*Correspondence: Ji Wang, Qi Wang,

†These authors have contributed equally to this work

This article was submitted to Molecular Medicine, a section of the journal Frontiers in Cell and Developmental Biology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics