Skip to main content
Download PDF

Long non-coding RNA expression profile in minor salivary gland of primary Sjögren’s syndrome

Arthritis Research & Therapy volume 18, Article number: 109 (2016) Cite this article

Abstract

Background

To examine the roles of long noncoding RNAs (lncRNAs) in the regulation of primary Sjögren’s syndrome (pSS) and reveal the expression profile of lncRNAs in labial salivary glands (LSGs) in pSS patients.

Method

The expression of 63,431 lncRNAs and 39,887 mRNAs were determined in the LSG of four pSS patients and four healthy controls using microarray experiments. Validation was performed in 30 pSS patients and 16 controls using real-time PCR. LncRNA-mRNA co-expression and gene-pathway networks were constructed using bioinformatics software.

Result

A total of 1243 lncRNAs (upregulated: 890, downregulated: 353) and 1457 mRNAs (upregulated: 1141, downregulated: 316) were differentially expressed in the LSGs of pSS patients (fold change >2, P <0.05). Eight of these lncRNAs were validated using real-time PCR. ENST00000420219.1 (3.13-fold), ENST00000455309.1 (2.51-fold), n336161 (2.45-fold), NR_002712 (2.41-fold), ENST00000546086.1 (1.94-fold), Lnc-UTS2D-1:1 (1.79-fold), n340599 (1.69-fold), and TCONS_l2_00014794 (1.28-fold) were significantly upregulated in pSS. There were strong correlations between these lncRNAs and β2 microglobulin, disease course, erythrocyte sedimentation rate (ESR), rheumatoid factor (RF), IgA, IgM, visual analogue scale (VAS) of parotid swelling and VAS of dry eyes. Computational analyses revealed that 28 of the differentially expressed (DE) mRNAs were associated with eight DE lncRNAs involved in chemokine signaling pathways, the nuclear factor-kappa B (NF-κB) signaling pathway, and tumor necrosis factor (TNF) signaling pathway.

Conclusions

Our study revealed the expression profile of lncRNAs in LSGs of pSS patients. Many novel lncRNA transcripts that play important roles in the pathogenesis of pSS were dysregulated in pSS. Therefore, this study will aid in the development of new diagnostic biomarkers and drug therapies.

Background

Primary Sjögren’s syndrome (pSS) is an autoimmune disease that is characterized by the dysfunction of exocrine glands, primarily the salivary and lachrymal glands, which results in dry mouth and eyes. The extensive infiltration of sensitized lymphocytes into target glands leads to the primary pathological manifestations of pSS. The etiology of pSS not clear, but previous studies have suggested that hereditary, hormonal, and environmental factors play crucial roles in the onset and progression of pSS [1]. Activation of the innate/adaptive immune system is the first line of defense against infections and damaged tissues. However, aberrantly activated inflammatory processes underlie autoimmune disease. Persistent perturbations of these inflammatory pathways are detrimental to the host and eventually result in disease conditions, including rheumatoid arthritis (RA), cardiovascular disease, and cancer [2–4]. Therefore, inflammatory signaling pathways require strict regulation at the transcriptional and posttranscriptional levels. Recent studies have demonstrated that noncoding RNAs, especially microRNAs, play crucial roles in the regulation of inflammatory signaling pathways [5, 6]. The roles of long noncoding RNAs (lncRNAs) as novel regulators of these pathways have emerged in more recent years. LncRNAs are a newly discovered class of regulatory molecules that are not translated into proteins and consist of more than 200 nucleotides. LncRNAs strongly affect a variety of biological processes in cells and organ systems. Several studies revealed the strong involvement of lncRNAs in the regulation of the immune response, including several pathways of innate immunity [7, 8]. LincRNA-Cox2 is one lncRNA that both negatively and positively regulates the expression of many important immune-related genes. This lncRNA mediates its repressive effects on interferon-stimulated genes via interactions with hnRNP-A/B and A2/B1 [9, 10]. Recent studies have also identified THRIL as a key regulator of tumor necrosis factor alpha (TNF-α) induction on Toll-like receptor (TLR) 1/2 signaling in human THP1 macrophages. Knockdown of THRIL downregulated the production of TNF-α mRNA [11]. Several other lncRNAs, such as NEAT [12, 13], lnc-IL7R [14], PACER [15], lnc-DC [16], IL1β-RBT46 [17] and AS-IL1α [18], control the innate immune response, immune cell development, and adaptive immunity. These studies indicate the crucial role of lncRNAs in the normal immune system. However, it is important to identify whether lncRNA dysfunction is involved in the pathogenesis of autoimmune diseases. Increasing evidence suggests that the dysregulation of lncRNAs plays an important role in autoimmune diseases, such as systemic lupus erythematosus (SLE), RA, type I diabetes mellitus (T1DM), and multiple sclerosis (MS) [19–24]. The present study analyzed the lncRNA expression profiles in labial salivary glands (LSGs) of pSS patients using lncRNA microarray to investigate the potential roles of lncRNAs in the pathogenesis of pSS. The results provide a new direction for the diagnosis and therapy of pSS.

Methods

Study subjects

Thirty patients diagnosed with pSS and 16 control subjects were recruited from the Department of Oral Surgery, Shanghai Ninth People’s Hospital, School of Medicine, Shanghai Jiao Tong University. All of the selected pSS patients fulfilled the American-European consensus group criteria for pSS. Specimens were collected during labial biopsy. No immunosuppressive treatment was administered to patients prior to diagnosis. Sixteen control subjects were diagnosed with labial gland mucocele and underwent mucocele excision. Normal labial glands adjacent to mucocele were collected during surgery. The absence of acute infections or systemic diseases was confirmed in all healthy donors. All specimens were immediately frozen in liquid nitrogen after resection and stored at −80 °C until RNA extraction. The discovery cohort was composed of eight patients mixed with pSS patients and control subjects for the screening of 63,431 lncRNA transcripts. The independent validation cohort consisted of 30 pSS patients and 16 normal samples. Table 1 summarizes the detailed demographic, clinical and laboratory characteristics of the 30 pSS patients. All healthy donors were females aged 25–65 years old. Consent was obtained from each participant prior to sample collection. The Ethics Committee, Faculty of Medicine, Shanghai Jiao Tong University approved this study.

Table 1 The detailed demographic, clinical, and laboratory characteristics of 30 pSS patients
Full size table

Microarray

Total RNA was extracted from eight samples (four pSS and four control subjects) using TRIzol reagent (Life Technologies, Carlsbad, CA, USA). The RNeasy Mini Kit (Qiagen, GmBH, Hilden, Germany) was used to purify the total RNA according to the manufacturer’s recommendation. Purified total RNA was quantified using a NanoDrop 1000 (Thermo Fisher Scientific, Waltham, MA, USA). The assessment of RNA integrity was determined using RNA LabChip™ kits and an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Only samples with 2100 RIN ≥7.0 and 28S/18S ≥0.7 were used.

Agilent SurePrint G3 microarray was used to investigate 63,431 lncRNAs and 39,887 mRNAs. Total RNA was amplified and labeled using a Low Input Quick Amp Labeling Kit, One-Color (Agilent Technologies, Santa Clara, CA, US). Labeled cRNA were purified using an RNeasy mini kit (Qiagen, GmBH, Hilden, Germany). The microarray hybridization was performed based on the manufacturer’s standard protocols (Agilent Technologies, Santa Clara, CA, US). Slides were washed in staining dishes and scanned. Raw data were normalized using a Quantile algorithm in Gene Spring Software 11.0 (Agilent Technologies, Santa Clara, CA, US).

LncRNA-mRNA co-expression network

An lncRNA-mRNA co-expression network was built based on the normalized signal intensity of differentially expressed lncRNAs and mRNAs to explore the dysregulation of lncRNAs in pSS patients. A co-expression network of control subjects and pSS patients was established. Differentially expressed lncRNAs and mRNAs that met the criteria (P value <0.05, fold change >2 or <0.5) were selected. Correlations between lncRNA and lncRNA, lncRNA and mRNA, mRNA and mRNA were investigated using Pearson’s correlations. Only strong correlations (P <0.001) were drawn in these renderings. The importance of a gene in this network is reflected by degree. A gene with a large degree indicated that it was at a central position in the network and it shared closer relationships with more genes. The last step investigated genes that exhibited different degrees in pSS patients and control subjects.

Screening of differentially expressed genes for validation

We screened differentially expressed (DE) lncRNAs for further validation using the following two approaches to validate the results of microarray experiments in an independent cohort and investigate the correlations between gene expression levels and clinical characteristics.

First, genes were evaluated based on the data revealed by microarray experiments using the following criteria: (a) the fold change of genes must be greater than fivefold compared to control subjects; (b) the signal value of the probes in each sample must be greater than seven; (c) the signal of the probes in each sample must be significantly different from the background signal; and (d) genes with lncRNA-mRNA repeated sequences and without information in databases were excluded.

Second, the degree of differentially expressed genes between pSS patients and control subjects was compared based on the co-expression network results. The top 30 genes with high different degrees and the qualified requirements for signal values were selected.

Real-time PCR

Total RNA was extracted as described above. cDNA was synthesized from 0.5 μg RNA using the iScript cDNA synthesis kit (Bio-Rad Laboratories, Hercules, CA, USA). Real-time PCR was performed using an ABI Power SYBR Green PCR Master Mix (ABI, Foster City, CA, USA) and 7900 HT Sequence Detection System (ABI, Foster City, CA, USA). Real-time PCR was performed with a 5-ng cDNA template using the 2× SYBR Green PCR buffer and 10-μmol PCR primers in a total volume of 10 μl. Reactions were performed in 384-well PCR microplates. Additional file 1: Table S1 lists the primer sequences.

The expression of each lncRNA was represented as fold changes using the △△Ct method to obtain quantitative results. Differences between groups were analyzed using a two-tailed Mann-Whitney U test or unpaired t test, based on the homogeneity of variance. Spearman’s test was used for correlation studies. A value of P <0.05 was considered significant.

Immunohistochemistry

Immunohistochemical staining to detect C-X-C chemokine receptor type 4 (CXCR4), CD19, CD21, Toll-like receptor 9 (TLR9) and intercellular cell adhesion molecule 1 (ICAM1) was performed on LSG biopsy sections as described previously [25]. CXCR4 (ab124824, Abcam, Cambridge, MA, USA), CD3 (ab699, Abcam, Cambridge, MA, USA), CD19 (ab134114, Abcam, Cambridge, MA, USA), CD21 (ab75985, Abcam Cambridge, MA, USA), ICAM1 (ab53013, Abcam, Cambridge, MA, USA), CD20 (ab78237, Abcam, Cambridge, MA, USA), and TLR9 (BA3861-1, Boster Biotechnology, Wuhan, China) antibodies were used in this experiment. Negative control staining was performed by replacing primary antibodies with PBS. Positive immunoreactivity appeared as a brown color. Double staining for CD3 and CD20 was used to analyze T/B cell segregation using the DouMaxvision™ double-stain system (KIT-9998, Maixin Biotechnology, Fuzhou, China). A scoring system was used to describe the results of double staining as described previously [26]. The images of grade 1 to 3 are shown in Additional file 2: Figure S3.

EBV-encoded RNA (EBER) in situ hybridization

In situ hybridization was performed in LSG of 30 pSS patients using Epstein-Barr Virus ISH Detection Kit (Triplex International Biosciences, Fuzhou, China), according to the manufacturer’s protocol. Patients with lymphepithelioma were used as positive controls. The images of EBER-positive cells in representative ectopic lymphoid structures are shown in Additional file 3: Figure S4.

Bioinformatics analysis

Aberrantly expressed lncRNAs and mRNAs with statistical significance were identified using Volcano Plot filtering. The threshold used to screen up- or downregulated RNAs was a fold change >2.0 (P <0.05). Hierarchical clustering was performed in Cluster 3.0, and heat maps were generated in Java Treeview. The DE mRNAs or DE lncRNAs-related mRNAs were analyzed using pathway annotation and gene ontology (GO) functional enrichment in the Cytoscape 3.3 software. The gene-pathway network was also constructed using Cytoscape 3.3.

Results

Overview of aberrantly expressed lncRNAs and mRNAs in pSS

The expression levels of mRNAs and lncRNAs in four LSGs of pSS patients and paired control samples were analyzed using gene expression microarrays. The expression signatures of lncRNAs and mRNAs were reviewed using scatter plot and hierarchical clustering analyses. The scatter plots revealed that many lncRNAs and mRNAs were differentially expressed between pSS patients and control subjects (Fig. 1a, b). The heat maps of DE lncRNAs or mRNAs indicated the high level of concordance in pSS or control samples (Fig. 1c, d). The above-listed data indicated that changes in lncRNAs and mRNAs in LSGs were associated with the pathogenesis of pSS. Of the 63,431 lncRNAs and 39,887 mRNAs in the microarray, 1243 lncRNAs and 1457 mRNAs were significantly differentially expressed in the labial glands of pSS patients compared to control subjects (fold change >2, P <0.05). A total of 890 lncRNAs and 1141 mRNAs were upregulated, and 353 lncRNAs and 316 mRNAs were downregulated. Table 2 lists the 20 most up- and downregulated DE lncRNAs. Table 3 lists the 20 most up- and downregulated DE mRNAs. All of the data from the microarray trials are stored in the GEO database with accession number GSE76013.

Fig. 1
figure 1

The expression profiling of lncRNA and mRNA in LSGs of pSS patients and healthy controls. a, b Scatter plot comparing global lncRNA or mRNA gene expression profiles between pSS and control subjects. Green lines and red line indicate twofold differences in either direction in lncRNA and mRNA expression. c, d Heat map showing hierarchical clustering of lncRNAs or mRNA with expression changes greater than twofold. Red and green colors represent up- and downregulated genes, respectively. lncRNA long noncoding RNA, SjS Sjögren’s syndrome

Full size image
Table 2 Top 20 significantly differential expressed lncRNAs between pSS patients and control subjects
Full size table
Table 3 Top 20 significantly differential expressed mRNAs between pSS patients and control subjects
Full size table

Real-time PCR: lncRNA expression profiles and correlations with clinical characteristics

Aberrantly expressed lncRNAs were identified in labial gland tissue samples from 30 pSS patients and 16 control subjects. The selection of representative lncRNAs for validation was based on the two approaches described in the Methods section. Nine lncRNAs were selected: NR_002712, n341833, lnc-UTS2D-1:1, TCONS_l2_00014794, n336161, ENST00000420219.1, ENST00000455309.1, n340599, and ENST00000546086.1. The results demonstrated that ENST00000420219.1 (3.13-fold), ENST00000455309.1 (2.51-fold), n336161 (2.45-fold), NR_002712 (2.41-fold), ENST00000546086.1 (1.94-fold), Lnc-UTS2D-1:1 (1.79-fold), n340599 (1.69-fold), and TCONS_l2_00014794 (1.28-fold) were significantly increased in pSS (Fig. 2). However, n341833 was not dysregulated in pSS patients. The next step was an analysis to determine whether any correlations existed between these DE lncRNA expression levels and the clinical characteristics. Figures 3a–d and 4a–b show that strong correlations were observed in eight pairs: (1) disease course with six lncRNA expression levels; (2) visual analogue scale (VAS) of parotid swelling with two lncRNAs; (3) VAS of dry eyes with one lncRNA; (4) β2 microglobulin with six lncRNAs; (5) erythrocyte sedimentation rate (ESR) with three lncRNAs; (6) rheumatoid factor (RF) with two lncRNAs; (7) immunoglobulin (Ig)A with two lncRNAs; and (8) IgM with one lncRNA. Four lncRNAs were significantly upregulated in SSB-positive patients compared to SSB-negative patients (Fig. 4c–d). No significant correlation existed between these DE lncRNAs and dry mouth, IgE, IgG, complement 3 (C3), complement 4 (C4), C-reactive protein (CRP) and grading of labial biopsy. Additional file 4: Table S2 shows the detailed correlation analysis results. A multivariate model was used to identify lncRNAs that correlated independently with disease characteristics. Additional file 5: Table S3 shows the results of this analysis.

Fig. 2
figure 2

The expression level of nine lncRNAs in LSG of pSS patients and control subjects were identified using real-time PCR. We confirmed that eight lncRNAs were significantly upregulated, except n341833. ENST00000420219.1 (3.13-fold, P <0.0001), ENST00000455309.1(2.51-fold, P <0.0001), n336161(2.45-fold, P = 0.0028), NR_002712(2.41-fold, P = 0.0005), ENST00000546086.1 (1.94-fold, P = 0.0047), Lnc-UTS2D-1:1(1.79-fold, P = 0.0428), n340599(1.69-fold, P = 0.0001), and TCONS_l2_00014794(1.28-fold, P = 0.0132) were significantly increased in pSS. pSS primary Sjögren’s syndrome

Full size image
Fig. 3
figure 3

Correlation between lncRNAs and clinical characteristics in pSS. a The expression level of ENST00000420219.1, ENST00000455309.1, n340599, n336161, ENST00000546086.1, and TCONS_l2_00014794 were significantly associated with the course of disease. b The expression level of ENST00000420219.1, ENST00000455309.1, NR_002712, n340599, ENST00000546086.1, and TCONS_l2_00014794 were strongly associated with β2 microglobulin. c The expression level of ENST00000455309.1, NR_002712, and n336161 were associated with ESR. d The expression levels of ENST00000455309.1 and ENST00000546086.1 were significantly associated with rheumatoid factor (RF). ESR erythrocyte sedimentation rate

Full size image
Fig. 4
figure 4

Correlation between DE lncRNAs and clinical characteristics in pSS. a The expression level of ENST00000420219.1 and NR_002712 were significantly correlated with VAS of parotid swelling, and the expression level of n336161 was associated with VAS of dry eyes. b The expression levels of ENST00000420219.1 and ENST00000455309.1 were significantly correlated with IgA, and the expression level of NR_002712 was associated with IgM. c The expression levels of ENST00000420219.1, NR_002712, ENST00000546086.1, and TCONS_l2_00014794 were significantly upregulated in SSB-positive patients compared with SSB-negative patients. Ig immunoglobulin, VAS visual analogue scale

Full size image

Functional prediction of DE mRNAs

GO and pathway analyses of the up- and downregulated genes in LSGs of pSS patients were performed to further examine the potential mechanism of pSS. The GO results indicated that the most significantly enriched cellular components of upregulated mRNAs in LSGs of pSS were the external side of the plasma membrane, immunological synapse, and MHC class II protein complex (Fig. 5a). The most significantly enriched molecular functions of upregulated mRNAs were peptide antigen binding, MHC class II receptor activity, and transmembrane signaling receptor activity (Fig. 5a). The most significantly enriched biological processes of upregulated mRNAs were immune response, inflammatory response, and cytokine-mediated signaling pathway (Fig. 5a). The most significantly enriched cellular components of downregulated mRNAs in LSGs of pSS were host cell nucleus, nucleus, and nucleolus (Fig. 5b). The most significantly enriched molecular functions of downregulated mRNAs were heparanase activity, monooxygenase activity, and mu-type opioid receptor binding (Fig. 5b). The most significantly enriched biological processes of downregulated mRNAs were transmembrane transport, mammary gland alveolus development, and potassium ion homeostasis (Fig. 5b). DE mRNAs were analyzed in the Kyoto Encyclopedia of Genes and Genomes (KEGG). The results revealed that the upregulated mRNAs in LSGs of pSS were significantly involved in graft-versus-host disease, cytokine-cytokine receptor interactions, and cell adhesion molecules (Fig. 5c). The downregulated mRNAs were significantly involved in gastric acid secretion, mineral absorption, and retinol metabolism (Fig. 5d).

Fig. 5
figure 5

Biological functions of DE mRNAs (P <0.05, fold changes >2). a The significant molecular function, biological process, and cellular component of upregulated mRNAs. b The significant molecular function, biological process, and cellular component of downregulated mRNAs. c Significantly enriched pathways of upregulated mRNAs. d Significantly enriched pathways of downregulated mRNAs

Full size image

Relational analyses of lncRNAs and mRNAs

Twenty-eight DE mRNAs from the co-expression network data strongly correlated with the eight validated lncRNAs (P <0.001, r >0.99). Table 4 shows the detailed information of these mRNAs. Figure 6 shows that gene-pathway network graph analyses revealed that the 28 DE mRNAs, including ICAM1, TLR9, TNF receptor-associated factor 1 (TRAF1), CXCR4, chemokine (C-C motif) ligand 20 (CCL20), and CD19 were likely involved in chemokine signaling pathways, the nuclear factor-kappa B (NF-Κb) signaling pathway, TNF signaling pathway, African trypanosomiasis, natural killer cell-mediated cytotoxicity, and Epstein-Barr virus infection [27–30], which are involved in the pathogenesis of pSS. Additional file 6: Figures S1 and Additional file 7: Figure S2 show that the results of immunohistochemistry in LSGs of pSS patients confirmed the overexpression of CD19, CXCR4, ICAM1, and TLR9 proteins. These results indicate the functional roles of DE lncRNAs in the progress of pSS. Since many of these lncRNA are supposed to regulate immune functions, the correlation between the eight upregulated lncRNAs and main histopathological data, such as the infection of Epstein-Barr virus, B/T cell segregation, and presence/absence of ectopic lymphoid structures, were analyzed. The results revealed no significant correlation exists (Additional file 8: Figure S5).

Table 4 Detailed information of DE lncRNAs linked DE mRNA
Full size table
Fig. 6
figure 6

Gene-pathway network graph of DE lncRNA-related mRNAs from Table 4. Circles and boxes represent DE lncRNAs-related genes and the corresponding pathways, respectively. The color of the pathway terms is defined by the P value of enrichment

Full size image

Discussion

LncRNAs were primarily investigated in genomic imprinting, cancers, and cell differentiation, but these molecules are emerging as important regulators of immune cell differentiation and the activation of innate immunity. However, the identification of lncRNA expression in autoimmune diseases is largely underexplored. Few studies reported that lncRNAs played a crucial role in autoimmune diseases, such as SLE, RA, T1DM and MS [31–36], but fewer studies have examined the expression profile of lncRNAs in pSS. The critical role of lncRNAs in the pathogenesis of pSS was demonstrated previously [37]. Therefore, differentially expressed lncRNAs in the LSGs of pSS patients were identified using microarray experiments.

The results demonstrated that upregulated mRNAs were far more numerous than downregulated mRNAs in pSS samples, which indicated the activation of many new biological processes or signaling pathways in pathological conditions. DE mRNAs were further analyzed in pSS patients using GO term enrichment and pathway enrichment analyses. The GO results indicated that the most significantly enriched cellular components of upregulated mRNAs in LSGs of pSS patients were the external side of plasma membrane, immunological synapse, and MHC class II protein complex. The most significantly enriched molecular functions of upregulated mRNAs were peptide antigen binding, MHC class II receptor activity, and transmembrane signaling receptor activity. These results are consistent with previous studies that demonstrated aberrantly expressed MHC class II and costimulation molecules in the epithelial cells of salivary glands in pSS patients [38, 39]. The abnormal expression of these molecules on the surface of salivary gland epithelial cells may favor the presentation of SSA and SSB epitopes to T cells and lead to autoantibody production [40]. The immune response, inflammatory response, and cytokine-mediated signaling pathway were the most significantly enriched biological processes of the upregulated mRNAs. The results suggest that an active autoimmune inflammatory response occurred in the epithelial cells of salivary glands in pSS patients. The results also demonstrated that several cellular components, molecular functions, and biological processes were inhibited. The pathway analysis results suggested that the upregulated mRNAs in LSGs of pSS patients were significantly involved in graft-versus-host diseases, cytokine-cytokine receptor interactions, and cell adhesion molecules. The results also suggest that deficiencies occurred in the labial gland epithelial cells of pSS patients and ended with the activation of autoimmune inflammation and the release of inflammatory cytokines by immune cells.

Nine of the aberrantly expressed lncRNAs, including NR_002712, n341833, lnc-UTS2D-1:1, TCONS_l2_00014794, n336161, ENST00000420219.1, ENST00000455309.1, n340599, and ENST00000546086.1, were further confirmed in the 30 pSS patients using real-time PCR. Eight lncRNAs were significantly upregulated, but n341833 was not upregulated. Recent prospective studies confirmed that serum β2 microglobulin levels were associated with the EULAR Sjögren’s syndrome disease activity index (ESSDAI) and EULAR SS patient-reported index (ESSPRI) in the pSS patients [41, 42]. Notably, the expression levels of NR_002712, ENST00000546086.1, TCONS_l2_00014794, n340599, ENST00000455309.1, and ENST00000420219.1 correlated with the β2 microglobulin levels, which suggest that the upregulated expression of these lncRNAs was associated with active disease states. These lncRNAs may be strongly involved in the progress of pSS, especially ENST00000455309.1, which significantly correlated with the disease course, ESR, RF, and IgA expression levels. These multiple correlations were also observed with other lncRNAs. The results further confirmed that these lncRNAs played critical roles in the pathophysiology of pSS, but further functional studies are needed to examine the potential mechanisms. The regulatory mechanisms and functional principles of lncRNAs were elucidated recently. LncRNAs that regulate the abundance of genomically neighboring or distal gene products are classified as cis-regulatory models or trans-regulatory models, respectively [43]. However, lncRNAs achieve regulatory functions via modularity, the collection of diverse combinations of proteins, and possible RNA or DNA interactions. This study investigated the potential targets of eight DE lncRNAs by using an lncRNA-mRNA co-expression network. This network revealed that 28 DE mRNAs were strongly associated with eight DE lncRNAs. A gene-pathway network of the eight DE lncRNA-related genes was established to explore the functional mechanism. The results revealed a significant involvement of these eight DE lncRNAs in several important signaling pathways that play crucial roles in the pathogenesis of pSS. However, the precise regulatory mechanism of lncRNAs requires further study. Animal models should be used in further research to help elucidate the biological processes of lncRNAs in pSS. Another limitation of this study was the lack of pSS mouse models. LncRNAs expression profiles should be detected in different animal models because of the variety of pSS mouse models.

Conclusions

This study revealed aberrant expression profiles of lncRNAs in LSGs of pSS patients for the first time. A total of 1243 lncRNAs and 1457 mRNAs were differentially expressed in the labial glands of pSS patients compared to the control subjects, and eight lncRNAs were further confirmed using real-time PCR in 30 pSS patients. Strong correlations between lncRNAs and clinical characteristics were observed. This study will lead to novel directions in pSS diagnosis and therapy.

Ethics approval and consent to participate

This study was approved by the Ethics Committee, Faculty of Medicine, Shanghai Jiao Tong University.

Abbreviations

C3:

complement 3

C4:

complement 4

CCL20:

chemokine (C-C motif) ligand 20

CRP:

C-reactive protein

CXCR4:

C-X-C chemokine receptor type 4

DE:

differentially expressed

EBER:

EBV-encoded RNA

ESR:

erythrocyte sedimentation rate

ESSDAI:

EULAR Sjögren’s syndrome disease activity index

ESSPRI:

EULAR SS patient-reported index

GO:

gene ontology

ICAM1:

intercellular cell adhesion molecule 1

Ig:

immunoglobulin

lncRNA:

long noncoding RNA

LSG:

labial salivary gland

MS:

multiple sclerosis

pSS:

primary Sjögren’s syndrome

RA:

rheumatoid arthritis

RF:

rheumatoid factor

SLE:

systemic lupus erythematosus

T1DM:

type I diabetes mellitus

TLR9:

Toll-like receptor 9

TNF-α:

tumor necrosis factor alpha

TRAF1:

TNF receptor-associated factor 1

VAS:

visual analogue scale

References

  1. Fox RI. Sjogren’s syndrome. Lancet. 2005;366(9482):321–31.

    Article  CAS  PubMed  Google Scholar 

  2. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. 2013;14(10):1014–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Gierut A, Perlman H, Pope RM. Innate immunity and rheumatoid arthritis. Rheum Dis Clin North Am. 2010;36(2):271–96.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Moore KJ, Sheedy FJ, Fisher EA. Macrophages in atherosclerosis: a dynamic balance. Nat Rev Immunol. 2013;13(10):709–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Li Y, Shi X. MicroRNAs in the regulation of TLR and RIG-I pathways. Cell Mol Immunol. 2013;10(1):65–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. O'Connell RM, Rao DS, Baltimore D. microRNA regulation of inflammatory responses. Annu Rev Immunol. 2012;30:295–312.

    Article  PubMed  Google Scholar 

  7. Fitzgerald KA, Caffrey DR. Long noncoding RNAs in innate and adaptive immunity. Curr Opin Immunol. 2014;26:140–6.

    Article  CAS  PubMed  Google Scholar 

  8. Heward JA, Lindsay MA. Long non-coding RNAs in the regulation of the immune response. Trends Immunol. 2014;35(9):408–19.

    Article  CAS  PubMed  Google Scholar 

  9. Sauvageau M, Goff LA, Lodato S, Bonev B, Groff AF, Gerhardinger C, et al. Multiple knockout mouse models reveal lincRNAs are required for life and brain development. eLife. 2013;2:e01749.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Carpenter S, Aiello D, Atianand MK, Ricci EP, Gandhi P, Hall LL, et al. A long noncoding RNA mediates both activation and repression of immune response genes. Science. 2013;341(6147):789–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Szmyrka-Kaczmarek M, Kosmaczewska A, Ciszak L, Szteblich A, Wiland P. Peripheral blood Th17/Treg imbalance in patients with low-active systemic lupus erythematosus. Postepy Hig Med Dosw (Online). 2014;68:893–8.

    Article  Google Scholar 

  12. Imamura K, Imamachi N, Akizuki G, Kumakura M, Kawaguchi A, Nagata K, et al. Long noncoding RNA NEAT1-dependent SFPQ relocation from promoter region to paraspeckle mediates IL8 expression upon immune stimuli. Mol Cell. 2014;53(3):393–406.

    Article  CAS  PubMed  Google Scholar 

  13. Hirose T, Virnicchi G, Tanigawa A, Naganuma T, Li R, Kimura H, et al. NEAT1 long noncoding RNA regulates transcription via protein sequestration within subnuclear bodies. Mol Biol Cell. 2014;25(1):169–83.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Cui H, Xie N, Tan Z, Banerjee S, Thannickal VJ, Abraham E, et al. The human long noncoding RNA lnc-IL7R regulates the inflammatory response. Eur J Immunol. 2014;44(7):2085–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Krawczyk M, Emerson BM. p50-associated COX-2 extragenic RNA (PACER) activates COX-2 gene expression by occluding repressive NF-kappaB complexes. eLife. 2014;3:e01776.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Wang P, Xue Y, Han Y, Lin L, Wu C, Xu S, et al. The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation. Science. 2014;344(6181):310–3.

    Article  CAS  PubMed  Google Scholar 

  17. IIott NE, Heward JA, Roux B, Tsitsiou E, Fenwick PS, Lenzi L, et al. Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes. Nat Commun. 2014;5:3979.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Chan J, Atianand M, Jiang Z, Carpenter S, Aiello D, Elling R, et al. Cutting edge: a natural antisense transcript, AS-IL1alpha, controls inducible transcription of the proinflammatory cytokine IL-1alpha. J Immunol. 2015;195(4):1359–63.

    Article  CAS  PubMed  Google Scholar 

  19. Cope AP, Feldmann M. Emerging approaches for the therapy of autoimmune and chronic inflammatory disease. Curr Opin Immunol. 2004;16(6):780–6.

    Article  CAS  PubMed  Google Scholar 

  20. Chatenoud L. Immune therapies of autoimmune diseases: are we approaching a real cure? Curr Opin Immunol. 2006;18(6):710–7.

    Article  CAS  PubMed  Google Scholar 

  21. Borchers AT, Leibushor N, Naguwa SM, Cheema GS, Shoenfeld Y, Gershwin ME. Lupus nephritis: a critical review. Autoimmun Rev. 2012;12(2):174–94.

    Article  PubMed  Google Scholar 

  22. Suarez-Gestal M, Calaza M, Endreffy E, Pullmann R, Ordi-Ros J, Sebastiani GD, et al. Replication of recently identified systemic lupus erythematosus genetic associations: a case–control study. Arthritis Res Ther. 2009;11(3):R69.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Yamanaka H. TNF as a target of inflammation in rheumatoid arthritis. Endocr Metab Immune Disord Drug Targets. 2015;15(2):129–34.

    Article  CAS  PubMed  Google Scholar 

  24. Gomez JA, Wapinski OL, Yang YW, Bureau JF, Gopinath S, Monack DM, et al. The NeST long ncRNA controls microbial susceptibility and epigenetic activation of the interferon-gamma locus. Cell. 2013;152(4):743–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Shi H, Yu CQ, Xie LS, Wang ZJ, Zhang P, Zheng LY. Activation of TLR9-dependent p38MAPK pathway in the pathogenesis of primary Sjogren’s syndrome in NOD/Ltj mouse. J Oral Pathol Med. 2014;43(10):785–91.

    Article  CAS  PubMed  Google Scholar 

  26. Barone F, Bombardieri M, Manzo A, Blades MC, Morgan PR, Challacombe SJ, et al. Association of CXCL13 and CCL21 expression with the progressive organization of lymphoid-like structures in Sjogren’s syndrome. Arthritis Rheum. 2005;52(6):1773–84.

    Article  CAS  PubMed  Google Scholar 

  27. Horvath S, Nazmul-Hossain AN, Pollard RP, Kroese FG, Vissink A, Kallenberg CG, et al. Systems analysis of primary Sjogren’s syndrome pathogenesis in salivary glands identifies shared pathways in human and a mouse model. Arthritis Res Ther. 2012;14(6):R238.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Greenwell-Wild T, Moutsopoulos NM, Gliozzi M, Kapsogeorgou E, Rangel Z, Munson PJ, et al. Chitinases in the salivary glands and circulation of patients with Sjogren’s syndrome: macrophage harbingers of disease severity. Arthritis Rheum. 2011;63(10):3103–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Hu S, Wang J, Meijer J, Ieong S, Xie Y, Yu T, et al. Salivary proteomic and genomic biomarkers for primary Sjogren’s syndrome. Arthritis Rheum. 2007;56(11):3588–600.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Lisi S, Sisto M, Lofrumento DD, D'Amore M. Sjogren’s syndrome autoantibodies provoke changes in gene expression profiles of inflammatory cytokines triggering a pathway involving TACE/NF-kappaB. Lab Invest. 2012;92(4):615–24.

    Article  CAS  PubMed  Google Scholar 

  31. Kino T, Hurt DE, Ichijo T, Nader N, Chrousos GP. Noncoding RNA gas5 is a growth arrest- and starvation-associated repressor of the glucocorticoid receptor. Sci Signal. 2010;3(107):ra8.

    PubMed  PubMed Central  Google Scholar 

  32. Haywood ME, Rose SJ, Horswell S, Lees MJ, Fu G, Walport MJ, et al. Overlapping BXSB congenic intervals, in combination with microarray gene expression, reveal novel lupus candidate genes. Genes Immun. 2006;7(3):250–63.

    Article  CAS  PubMed  Google Scholar 

  33. Koenders MI, van den Berg WB. Novel therapeutic targets in rheumatoid arthritis. Trends Pharmacol Sci. 2015;36(4):189–95.

    Article  CAS  PubMed  Google Scholar 

  34. Hakonarson H, Grant SF, Bradfield JP, Marchand L, Kim CE, Glessner JT, et al. A genome-wide association study identifies KIAA0350 as a type 1 diabetes gene. Nature. 2007;448(7153):591–4.

    Article  CAS  PubMed  Google Scholar 

  35. Wallace C, Smyth DJ, Maisuria-Armer M, Walker NM, Todd JA, Clayton DG. The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes. Nat Genet. 2010;42(1):68–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Burfoot RK, Jensen CJ, Field J, Stankovich J, Varney MD, Johnson LJ, et al. SNP mapping and candidate gene sequencing in the class I region of the HLA complex: searching for multiple sclerosis susceptibility genes in Tasmanians. Tissue Antigens. 2008;71(1):42–50.

    CAS  PubMed  Google Scholar 

  37. Sandhya P, Joshi K, Scaria V. Long noncoding RNAs could be potential key players in the pathophysiology of Sjogren’s syndrome. Int J Rheum Dis. 2015;18:898–905.

    Article  CAS  PubMed  Google Scholar 

  38. Jonsson R, Klareskog L, Backman K, Tarkowski A. Expression of HLA-D-locus (DP, DQ, DR)-coded antigens, beta 2-microglobulin, and the interleukin 2 receptor in Sjogren’s syndrome. Clin Immunol Immunopathol. 1987;45(2):235–43.

    Article  CAS  PubMed  Google Scholar 

  39. Manoussakis MN, Dimitriou ID, Kapsogeorgou EK, Xanthou G, Paikos S, Polihronis M, et al. Expression of B7 costimulatory molecules by salivary gland epithelial cells in patients with Sjogren’s syndrome. Arthritis Rheum. 1999;42(2):229–39.

    Article  CAS  PubMed  Google Scholar 

  40. Tsunawaki S, Nakamura S, Ohyama Y, Sasaki M, Ikebe-Hiroki A, Hiraki A, et al. Possible function of salivary gland epithelial cells as nonprofessional antigen-presenting cells in the development of Sjogren’s syndrome. J Rheumatol. 2002;29(9):1884–96.

    CAS  PubMed  Google Scholar 

  41. Pertovaara M, Korpela M. ESSPRI and other patient-reported indices in patients with primary Sjogren’s syndrome during 100 consecutive outpatient visits at one rheumatological clinic. Rheumatology (Oxford). 2014;53(5):927–31.

    Article  CAS  Google Scholar 

  42. Gottenberg JE, Seror R, Miceli-Richard C, Benessiano J, Devauchelle-Pensec V, Dieude P, et al. Serum levels of beta2-microglobulin and free light chains of immunoglobulins are associated with systemic disease activity in primary Sjogren’s syndrome. Data at enrollment in the prospective ASSESS cohort. PloS one. 2013;8(5):e59868.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154(1):26–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work is supported by Key Project of Shanghai Municipal Commission of Health and Family Planning (Project Number: 2014035), Science and Technology Commission of Shanghai Municipality (Project Number: 16ZR1419000).

Author information

Author notes

    Authors and Affiliations

    1. Department of Oral Surgery, Affiliated Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai Key Laboratory of Stomatology, Shanghai, China

      Huan Shi, Ningning Cao, Yiping Pu, Lisong Xie, Lingyan Zheng & Chuangqi Yu

    Authors
    1. Huan Shi
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Ningning Cao
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Yiping Pu
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Lisong Xie
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Lingyan Zheng
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Chuangqi Yu
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding authors

    Correspondence to Lingyan Zheng or Chuangqi Yu.

    Additional information

    Competing interests

    The authors declare that they have no competing interests.

    Authors’ contributions

    HS carried out the molecular genetic studies, participated in the sequence alignment, and drafted the manuscript. NC carried out the molecular genetic studies, the immunoperoxidase assay, and revised the manuscript. LX participated in the sequence alignment, collected tissue samples, and helped to draft the manuscript. YP carried out in situ hybridization experiment and collected tissue samples. LZ participated in the design of the study, performed the statistical analysis, and helped to revise the manuscript. CY conceived of the study, and participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.

    Additional files

    Additional file 1: Table S1.

    Primer sequences used in validation of lncRNAs. (DOCX 14 kb)

    Additional file 2: Figure S3.

    Characterization of the cellular infiltrate and histomorphologic grading of LSG. (A) Shown are representative examples of grade 1 (G1; <50 periductal lymphocytes); (B,C) grade 2 (G2; >50 periductal lymphocytes): nonsegregated and segregated aggregates (NS-G2 and S-G2, respectively); (D,F) grade 3 (G3; >50 periductal lymphocytes, with GC-like structures). (TIF 58923 kb)

    Additional file 3: Figure S4.

    EBER expression in LSG of pSS patients. (A) Significant number of EBER+ cells in representative ectopic lymphoid structure-positive LSG tissue. (B) Double staining for CD3 and CD20 in the same section of LSG tissue. (C) Staining for CD21 in the same section of LSG tissue. (TIF 53011 kb)

    Additional file 4: Table S2.

    Detailed correlation analysis results. (DOCX 20 kb)

    Additional file 5: Table S3.

    Multivariate model used to identify the lncRNAs that correlated independently with disease characteristics. (DOCX 14 kb)

    Additional file 6: Figure S1.

    Immunohistochemistry for CD19 and ICAM1 in LSG of pSS and healthy control. (A, E) Infiltrating lymphocytes in LSG of pSS were positively staining for CD19 (200× and 400× magnification); (B, F) glandular epithelial cells from control were negative (200× and 400× magnification); (C, G) infiltrating lymphocytes and adjacent glandular epithelial cells in LSG of pSS were positively staining for ICAM1 (200× and 400× magnification). (D, H) Glandular epithelial cells from control were negative (200× and 400× magnification). (TIF 28780 kb)

    Additional file 7: Figure S2.

    Immunohistochemistry for TLR9 and CXCR4 in LSG of pSS and healthy controls. (A, E) Infiltrating lymphocytes in LSG of pSS positively stained for TLR9 (200× and 400× magnification); (B, F) glandular epithelial cells from control were negative (200× and 400× magnification); (C, G) infiltrating lymphocytes and adjacent glandular epithelial cells in LSG of pSS were positively stained for CXCR4 (200× and 400× magnification). (D, H) Glandular epithelial cells from control were negative (200× and 400× magnification). (TIF 28419 kb)

    Additional file 8: Figure S5.

    Expression level analysis between the eight upregulated lncRNAs and main histopathological data. (A) Expression of eight lncRNAs in grade 1 (G1), nonsegregated grade 2 (NS-G2), segregated grade 2 (S-G2), and grade 3 (G3) foci in LSG from 30 patients with pSS. (B) Expression of eight lncRNAs in EBER-positive and -negative LSG from 30 patients with pSS. (TIF 2332 kb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Shi, H., Cao, N., Pu, Y. et al. Long non-coding RNA expression profile in minor salivary gland of primary Sjögren’s syndrome. Arthritis Res Ther 18, 109 (2016). https://doi.org/10.1186/s13075-016-1005-2

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s13075-016-1005-2

    Keywords

    Arthritis Research & Therapy

    ISSN: 1478-6362

    Contact us