Differential long non-coding RNA expression profile and function analysis in primary Sjogren's syndrome

Associations of placental lncRNA expression with maternal pre-pregnancy BMI and infant birthweight in two birth cohorts

Pre-pregnancy obesity (ppOB) is linked to pregnancy complications and abnormal fetal growth through placental mechanisms, and long non-coding RNAs (lncRNAs) may play an epigenetic role in these processes. We investigated overall and sex-specific associations of pre-pregnancy body mass index (ppBMI), ppOB, and birthweight with placental lncRNA transcripts in two birth cohorts. Study participants were mother-child dyads recruited to the CANDLE (Memphis, TN)(n = 725) and GAPPS (Seattle and Yakima, WA)(n = 159) cohorts. Maternal ppBMI was assessed at enrollment using interviewer-administered questionnaires. LncRNAs (1,077 and 1,033 for CANDLE and GAPPS, respectively) were sequenced from placental samples collected at birth. Placental lncRNA was regressed on ppBMI, ppOB (ppBMI ≥30kg/m2), or continuous birthweight in cohort-specific weighted linear models controlling for a priori-specified confounders and experimental variables. Potential effect modification by infant-sex was examined in sex-stratified analyses and models including BMI-infant-sex interaction terms. No lncRNA transcripts were significantly associated with ppBMI, ppOB, or birthweight in primary models. Among male infants in CANDLE, expression of three lncRNA transcripts (ERVH48-1, AC139099.1, CEBPA-DT) was associated with ppBMI and one transcript (AC104083.1) with birthweight. In GAPPS, ppBMI was associated with two lncRNA transcripts (AP000879.1 and AL365203.2) among males, and birthweight was associated with 17 lncRNA transcripts (including LINC02709, KANSL1-AS1, DANCR, EPB41L4A-AS1, and GABPB1-AS1) among females. No BMI-infant-sex interactions were observed. Though many of these potential associations are for uncharacterized transcripts, several identified lncRNAs (e.g., ERVH48-1 and CEBPA-DT) have been linked to pathways controlling cancer or placental growth, trophoblast differentiation, and gene expression. These associations warrant validation in future studies.

IFI27, a potential candidate molecular marker for primary Sjogren's syndrome

OBJECTIVE

The etiology of primary Sjogren's syndrome (pSS) is complex and not completely clear. This study was to identify key genes in pSS based on Gene Expression Omnibus (GEO).

METHODS

We downloaded the GSE40568, GSE80805, GSE127952, and GSE164885 mRNA expression profiles from GEO. Differentially expressed genes (DEGs) analyses were carried out by using the online analysis tool GEO2R and R. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to elucidate the biological processes, molecular function, cellular component, and KEGG signaling pathways for the DEGs in salivary glands (SGs) and peripheral blood mononuclear cells (PBMCs). Genes co-expressed were found in PBMCs and SGs of pSS patients. RT-qPCR was performed for validation. Finally, clinical correlation analysis and receiver operator characteristic (ROC) curve analysis were performed.

RESULTS

A total of thirty-nine up-regulated and one down-regulated genes were identified in pSS SGs. GO and KEGG pathway revealed that these DEGs were related to response to virus, and type I interferon signaling pathway. It was verified that fourteen genes were up-regulated in the SGs of pSS by RT-qPCR. Twenty up-regulated genes were identified in pSS patients PBMCs. Two genes were up-regulated in SGs and PBMCs of pSS patients, including IFI27 and IFI44L. The mRNA level of IFI27 was positively correlated with the disease activity of pSS patients. Furthermore, ROC analyses proved IFI27 may have diagnostic value for pSS.

CONCLUSION

IFI27 might serve as a potential biomarker for the early diagnosis and therapy of pSS. Key Points • IFI27 expression was significantly increased in PBMCs and SGs of pSS patients. • IFI27 was positively correlated with disease activity in pSS patients. • IFI27 might have a good diagnostic value for pSS.

Decoding and reconstructing disease relations between dry eye and depression: a multimodal investigation comprising meta-analysis, genetic pathways and Mendelian randomization

INTRODUCTION

The clinical presentations of dry eye disease (DED) and depression (DEP) often comanifest. However, the robustness and the mechanisms underlying this association were undetermined.

OBJECTIVES

To this end, we set up a three-segment study that employed multimodality results (meta-analysis, genome-wide association study [GWAS] and Mendelian randomization [MR]) to elucidate the association, common pathways and causality between DED and DEP.

METHODS

A meta-analysis comprising 26 case-control studies was first conducted to confirm the DED-DEP association. Next, we performed a linkage disequilibrium (LD)-adjusted GWAS and targeted phenotype association study (PheWAS) in East Asian TW Biobank (TWB) and European UK Biobank (UKB) populations. Single-nucleotide polymorphisms (SNPs) were further screened for molecular interactions and common pathways at the functional gene level. To further elucidate the activated pathways in DED and DEP, a systemic transcriptome review was conducted on RNA sequencing samples from the Gene Expression Omnibus. Finally, 48 MR experiments were implemented to examine the bidirectional causation between DED and DEP.

RESULTS

Our meta-analysis showed that DED patients are associated with an increased DEP prevalence (OR = 1.83), while DEP patients have a concurrent higher risk of DED (OR = 2.34). Notably, cross-disease GWAS analysis revealed that similar genetic architecture (rG = 0.19) and pleiotropic functional genes contributed to phenotypes in both diseases. Through protein-protein interaction and ontology convergence, we summarized the pleiotropic functional genes under the ontology of immune activation, which was further validated by a transcriptome systemic review. Importantly, the inverse variance-weighted (IVW)-MR experiments in both TWB and UKB populations (p value <0.001) supported the bidirectional exposure-outcome causation for DED-to-DEP and DEP-to-DED. Despite stringent LD-corrected instrumental variable re-selection, the bidirectional causation between DED and DEP remained.

CONCLUSION

With the multi-modal evidence combined, we consolidated the association and causation between DED and DEP.

Sjogren's syndrome: Everything you always wanted to know about genetic and epigenetic factors

Primary Sjögren's syndrome (pSS) is a chronic, systemic autoimmune disease characterized by a wide spectrum of glandular and extra-glandular features. Genetic and epigenetic factors play an important role in the disease susceptibility and phenotype. There are a multitude of genes that have been identified as implicated in the pathogenesis of pSS, both in HLA and extra-HLA regions with a strong contribution given by genes in interferon signalling pathways. Among the HLA alleles, the most consistent associations have been found with DR2 and DR3 alleles at the DRB1 locus. Moreover, several gene variants outside the MHC locus are in genes involved in NF-κB signalling, B- and T-cell function and methylation processes possibly responsible for lymphomagenesis. There is still a lack of knowledge on precise genetic patterns and prediction models of diseases, and data on pharmacogenetics is scarce. A comprehensive summary of the common genetic factors and an extensive analysis of novel epigenetic aspects is provided, together with a view on the relationships between novel therapeutic agents for pSS and genetic targets in signalling pathways, aiming at improving tailored treatment strategies in the view of a more personalized medicine.

LncRNA SNHG1 serves as a biomarker for systemic lupus erythematosus and participates in the disease progression

LncRNAs play an important role in autoimmune diseases. The purpose of this study was to explore the role of lncRNA SNHG1 in systemic lupus erythematosus (SLE), and laid a theoretical foundation for the study of SLE. The basic clinical information of all subjects was first collected for statistical analysis, and SNHG1 expression in the serum of all subjects was detected by RT-qPCR. The value of SNHG1 in the diagnosis of SLE was assessed by ROC. The correlation between SNHG1 and each blood sample index was analyzed by Pearson correlation analysis. The role of SNHG1 in primary peripheral blood mononuclear cells (PBMCs) apoptosis was explored. SNHG1 expression is relatively upregulated in patients with SLE compared to healthy people. SNHG1 expression was positively correlated with SLEDAI score, IgG, CRP, and ESR, and negatively correlated with C3 and C4. ROC indicated that SNHG1 has the potential to assist in the diagnosis of SLE. PBMCs apoptosis in SLE was higher than that in control group, the knockdown and overexpression of SNHG1 could correspondingly inhibit and promote PBMCs apoptosis. SNHG1 has the potential to be a diagnosis marker for SLE and may be involved in regulating PBMCs apoptosis.

Long Non-Coding RNAs in Sjögren's Disease

Sjögren's disease (SjD) is a heterogeneous autoimmune disease characterized by severe dryness of mucosal surfaces, particularly the mouth and eyes; fatigue; and chronic pain. Chronic inflammation of the salivary and lacrimal glands, auto-antibody formation, and extra-glandular manifestations occur in subsets of patients with SjD. An aberrant expression of long, non-coding RNAs (lncRNAs) has been described in many autoimmune diseases, including SjD. Here, we review the current literature on lncRNAs in SjD and their role in regulating X chromosome inactivation, immune modulatory functions, and their potential as biomarkers.

CXCL9 may serve as a potential biomarker for primary Sjögren's syndrome with extra-glandular manifestations

BACKGROUND

Primary Sjögren's syndrome (pSS) is an autoimmune condition that causes harm to exocrine glands and also has extra-glandular manifestations (EGM). pSS patients with EGM have a worse prognosis than those with only sicca symptoms. Previous studies have shown that the minor salivary glands (MSG) of pSS patients exhibit a unique profile of cytokines and chemokines compared to healthy controls. However, there is a lack of research comparing pSS with EGM (pSS-EGM) and pSS without EGM (pSS-non-EGM). This study aims to explore potential biomarkers associated with pSS, particularly pSS with EGM.

METHODS

By utilizing RNA sequencing, we conducted an analysis on the gene expression profiles of MSG in 63 patients diagnosed with pSS, as well as 12 non-pSS individuals. Furthermore, we also investigated the MSG of pSS patients, both with and without EGM. Through bioinformatics analysis, we identified genes with differential expression (DEGs) and determined the core hub genes using PPI network. We then analyzed the top 20 DEGs and their correlation with the patients' clinical characteristics, and validated our findings using peripheral blood plasma.

RESULTS

A total of 725 differentially expressed genes (DEGs) were identified in the comparison between pSS and non-pSS groups, and 727 DEGs were observed between pSS-EGM and pSS-non-EGM. It is noteworthy that the expression levels of CXCL9 were higher in both pSS patients and pSS-EGM when compared to the control group. Taking into consideration the significance of the top 20 DEGs in relation to clinical parameters and the central hub genes, we ultimately chose CXCL9. In comparison to the non-pSS group, pSS patients exhibited notably greater expression of the CXCL9 gene in the MSG, as well as higher levels of CXCL9 protein in their plasma (p < 0.001). Furthermore, the expression of the CXCL9 gene and levels of CXCL9 protein were notably higher in pSS patients accompanied by EGM and those with SSA antibodies. Additionally, a correlation was found between the expression of the CXCL9 gene and the EULAR Sjogren's Syndrome Disease Activity Index (ESSDAI), as well as with immunoglobulin G (IgG) levels and erythrocyte sedimentation rate (ESR). Meanwhile, the protein levels of CXCL9 were found to be correlated with IgG levels and ESSDAI.

CONCLUSION

CXCL9 proves to be a valuable biomarker in pSS, specifically due to its strong ability to differentiate between pSS patients with EGM and those without EGM. There is a significant correlation between CXCL9 and various clinical parameters both at the gene and protein level. Therefore, CXCL9 could be a potential target for future treatment of pSS.

Application of omics in Sjögren's syndrome

OBJECTIVE

The pathogenesis, diagnosis, and treatment of Sjögren's syndrome (SS) face many challenges, and there is an urgent need to develop new technologies to improve our understanding of SS.

METHODS

By searching the literature published domestically and internationally in the past 20 years, this artical reviewed the research of various omics techniques in SS.

RESULTS

Omics technology provided valuable insights into the pathogenesis, early diagnosis, condition and efficacy evaluation of SS. It is helpful to reveal the pathogenesis of the disease and explore new treatment schemes, which will open a new era for the study of SS.

CONCLUSION

At present, omics research has made some gratifying achievements, but there are still many uncertainties. Therefore, in the future, we should improve research techniques, standardize the collection of samples, and adopt a combination of multi-omics techniques to jointly study the pathogenesis of SS and provide new schemes for its treatment.

Identification and functional prediction of long non-coding RNA and mRNA related to connective tissue disease-associated interstitial lung diseases

Objective

Recently, the role of long non-coding RNA (lncRNA) in rheumatic immune diseases has attracted widespread attention. However, knowledge of lncRNA in connective tissue disease-associated interstitial lung disease (CTD-ILD) is limited. This study explored the expression profile and possible mechanisms of lncRNA and mRNA in peripheral blood mononuclear cells (PBMCs) of CTD-ILD patients, especially systemic sclerosis (SSc)-ILD and rheumatoid arthritis (RA)-ILD.

Methods

LncRNA microarray analysis identified 240 diferentially expressed lncRNAs and 218 diferentially expressed mRNA in the CTD-ILD group and the connective tissue disease without associated interstitial lung disease (CTD-NILD) group. The bioinformatics analysis of diferential genes has identified several important biological processes and signal pathways, including nuclear factor kappa B (NF-kappa B) signaling pathway, interleukin 17 (IL-17) signaling pathway, B cell receptor signaling pathway. Relative expression levels of five diferentially expressed lncRNAs and one mRNA in 120 SSc and RA patients with or without ILD were detected by quantitative reverse-transcription (PCR).

Results

The ENST00000604692 expression level was significantly higher in the ILD than the without interstitial lung disease (NILD) group; T311354 and arginase-1 were significantly higher in SSc than RA group.

Conclusion

These data suggest that the specific profile of lncRNA in PBMCs of CTD-ILD patients and the potential signal pathways related to the pathogenesis of CTD-ILD, which may provide newfound insights for the diagnosis and treatment of CTD-ILD patients.

Current Views on Pathophysiology and Potential Therapeutic Targets in Sjögren's Syndrome: A Review from the Perspective of Viral Infections, Toll-like Receptors, and Long-Noncoding RNAs

Sjögren's syndrome (SS) is a rheumatic disease characterized by sicca and extraglandular symptoms, such as interstitial lung disease and renal tubular acidosis. SS potentially affects the prognosis of patients, especially in cases of complicated extraglandular symptoms; however, only symptomatic therapies against xerophthalmia and xerostomia are currently included in the practice guidelines as recommended therapies for SS. Considering that SS is presumed to be a multifactorial entity caused by genetic and environmental factors, a multidisciplinary approach is necessary to clarify the whole picture of its pathogenesis and to develop disease-specific therapies for SS. This review discusses past achievements and future prospects for pursuing the pathophysiology and therapeutic targets for SS, especially from the perspectives of viral infections, toll-like receptors (TLRs), long-noncoding RNAs (lncRNAs), and related signals. Based on the emerging roles of viral infections, TLRs, long-noncoding RNAs and related signals, antiviral therapy, hydroxychloroquine, and vitamin D may lower the risk of or mitigate SS. Janus-kinase (JAK) inhibitors are also potential novel therapeutic options for several rheumatic diseases involving the JAK-signal transducer and activator of transcription pathways, which are yet to be ascertained in a randomized controlled study targeting SS.

Major depression-related factor NEGR1 controls salivary secretion in mouse submandibular glands

Salivary gland cells, which secrete water in response to neuronal stimulation, are closely connected to other neurons. Transcriptomic studies show that salivary glands also express some proteins responsible for neuronal function. However, the physiological functions of these common neuro-exocrine factors in salivary glands are largely unknown. Here, we studied the function of Neuronal growth regulator 1 (NEGR1) in the salivary gland cells. NEGR1 was also expressed in mouse and human salivary glands. The structure of salivary glands of Negr1 knockout (KO) mice was normal. Negr1 KO mice showed tempered carbachol- or thapsigargin-induced intracellular Ca²⁺ increases and store-operated Ca²⁺ entry. Of interest, the activity of the large-conductance Ca²⁺-activated K⁺ channel (BK channel) was increased, whereas Ca²⁺-activated Cl^- channel ANO1 channel activity was not altered in Negr1 KO mice. Pilocarpine- and carbachol-induced salivation was decreased in Negr1 KO mice. These results suggest that NEGR1 influence salivary secretion though the muscarinic Ca²⁺ signaling.

Novel regulatory role of non-coding RNAs in ankylosing spondylitis

Ankylosing spondylitis (AS) is a type of arthritis that primarily affects the spine and involves disorders of the immune and skeletal systems. However, the exact pathogenesis of AS is not fully understood. Non-coding RNAs (ncRNAs), particularly, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and micro RNAs (miRNAs) and their interactions have been shown to influence many biological processes such as inflammatory responses, osteogenic differentiation and apoptosis, pyroptosis, and proliferation. In addition, ncRNAs reflect the disease activity of AS. In this review, we discuss the regulatory roles of ncRNAs in AS cell functions (inflammatory responses, cellular osteogenic differentiation and apoptosis, pyroptosis, and proliferation) and their potential applications in AS diagnosis and treatment. Understanding the role of ncRNAs in the pathogenesis of AS will lay the foundation for exploring potential new therapeutic approaches for AS.

Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses

OBJECTIVE

Sjögren's disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised.

METHODS

Whole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo-); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays.

RESULTS

Whole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene, RSAD2, in SjDRo+ (r≥0.65 or ≤-0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells.

CONCLUSION

LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Construction, bioinformatics analysis, and validation of competitive endogenous RNA networks in ulcerative colitis

Background: Ulcerative colitis (UC) is a common chronic disease of the digestive system. Recently, competitive endogenous RNAs (ceRNAs) have been increasingly used to reveal key mechanisms for the pathogenesis and treatment of UC. However, the role of ceRNA in UC pathogenesis has not been fully clarified. This study aimed to explore the mechanism of the lncRNA-miRNA-mRNA ceRNA network in UC and identify potential biomarkers and therapeutic targets.

Materials and Methods: An integrative analysis of mRNA, microRNA (miRNA), and long non-coding RNA (lncRNA) files downloaded from the Gene Expression Omnibus (GEO) was performed. Differentially expressed mRNA (DE-mRNAs), miRNA (DE-miRNAs), and lncRNA (DE-lncRNAs) were investigated between the normal and UC groups by the limma package. A weighted correlation network analysis (WGCNA) was used to identify the relative model for constructing the ceRNA network, and, concurrently, miRWalk and DIANA-LncBase databases were used for target prediction. Consecutively, the Gene Ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) pathway, and Reactome pathway enrichment analyses, protein-protein interaction (PPI) network, Cytohubba, and ClueGO were performed to identify hub genes. Additionally, we examined the immune infiltration characteristics of UC and the correlation between hub genes and immune cells using the immuCellAI database. Finally, the expression of potential biomarkers of ceRNA was validated via qRT-PCR in an experimental UC model induced by dextran sulfate sodium (DSS).

Result: The ceRNA network was constructed by combining four mRNAs, two miRNAs, and two lncRNAs, and the receiver operating characteristic (ROC) analysis showed that two mRNAs (CTLA4 and STAT1) had high diagnostic accuracy (area under the curve [AUC] > 0.9). Furthermore, CTLA4 up-regulation was positively correlated with the infiltration of immune cells. Finally, as a result of this DSS-induced experimental UC model, CTLA4, MIAT, and several associate genes expression were consistent with the results of previous bioinformatics analysis, which proved our hypothesis.

Conclusion: The investigation of the ceRNA network in this study could provide insight into UC pathogenesis. CTLA4, which has immune-related properties, can be a potential biomarker in UC, and MIAT/miR-422a/CTLA4 ceRNA networks may play important roles in UC.

Expression of LINC00847 in Peripheral Blood Mononuclear Cells of Children with Asthma and Its Prediction between Asthma Exacerbation and Remission

Objective. Asthma is defined as a heterogeneous disease that is usually characterized by chronic airway inflammation. Long noncoding RNAs play important roles in various biological processes including inflammation. To know more about the relationships between lncRNAs and asthma, we sought to the role of LINC00847 in peripheral blood mononuclear cells (PBMCs) of children with asthma exacerbation or asthma remission. Methods. Microarray analysis was performed on GSE143192 and GSE165934 datasets to screen differentially expressed lncRNAs (DElncRNAs) in human PBMCs between asthma patients and normal controls. LINC00847 was selected from DElncRNAs in human PBMCs between asthma patients and normal controls for further investigation. The expression levels of LINC00847 were quantified in PBMCs collected from 54 children with asthma exacerbation, 54 children with asthma remission, and 54 healthy children by real-time qPCR. The forced expiratory volume in the first second in percent predicted values (FEV1%), ratio of forced expiratory volume in 1 second to forced vital capacity (FEV1/FVC), and peak expiratory flow rate (PEF%) were tested for evaluation of lung function. The concentration of immunoglobulin E (IgE) and eosinophil count was examined. The serum levels of interleukin-4 (IL-4), interferon-γ (IFN-γ), and IL-17A were determined by the ELISA method. Results. The expression level of LINC00847 in PBMCs of asthma exacerbation children was remarkably higher than that in PBMCs of asthma remission children and healthy children ( $p<0.001$ ); the expression level of LINC00847 in PBMCs of asthma remission children was notably higher than that in PBMCs of healthy children ( $p<0.001$ ). Pearson correlation analysis revealed that the expression levels of LINC00847 in PBMCs of asthma children were negatively correlated with FEV1% (r = −0.489), FEV1/FVC (r = −0.436), PEF% (r = −0.626), and IFN-γ level (r = −0.614) of asthma children, but positively correlated with IgE concentration (r = 0.680), eosinophil count (r = 0.780), IL-4 (r = 0.524), and IL-17A (r = 0.622) levels. When LINC00847 expression was used to distinguish asthma exacerbation from asthma remission, a 0.871 AUC (95% CI: 0.805–0.936) was yielded with sensitivity of 79.63% and specificity of 77.78%. Conclusion. The study demonstrates that increased LINC00847 expression may be associated with the development and progression of asthma, possibly serving as a novel biomarker for predicting asthma exacerbation from asthma remission.