Dataset Information

Integrated Bioinformatics and Validation Reveal Potential Biomarkers Associated With Progression of Primary Sjogren's Syndrome.

ABSTRACT:

Background

Primary Sjögren's syndrome (pSS) is a chronic systemic autoimmune disease of the exocrine glands characterized by specific pathological features. Previous studies have pointed out that salivary glands from pSS patients express a unique profile of cytokines, adhesion molecules, and chemokines compared to those from healthy controls. However, there is limited evidence supporting the utility of individual markers for different stages of pSS. This study aimed to explore potential biomarkers associated with pSS disease progression and analyze the associations between key genes and immune cells.

Methods

We combined our own RNA sequencing data with pSS datasets from the NCBI Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) via bioinformatics analysis. Salivary gland biopsies were collected from 14 pSS patients, 6 non-pSS patients, and 6 controls. Histochemical staining and scanning electron micrographs (SEM) were performed to macroscopically and microscopically characterize morphological features of labial salivary glands in different disease stages. Then, we performed quantitative PCR to validate hub genes. Finally, we analyzed correlations between selected hub genes and immune cells using the CIBERSORT algorithm.

Results

We identified twenty-eight DEGs that were upregulated in pSS patients compared to healthy controls. These were mainly involved in immune-related pathways and infection-related pathways. According to the morphological features of minor salivary glands, severe interlobular and periductal lymphocytic infiltrates, acinar atrophy and collagen in the interstitium, nuclear shrinkage, and microscopic organelle swelling were observed with pSS disease progression. Hub genes based on above twenty-eight DEGs, including MS4A1, CD19, TCL1A, CCL19, CXCL9, CD3G, and CD3D, were selected as potential biomarkers and verified by RT-PCR. Expression of these genes was correlated with T follicular helper cells, memory B cells and M1 macrophages.

Conclusion

Using transcriptome sequencing and bioinformatics analysis combined with our clinical data, we identified seven key genes that have potential value for evaluating pSS severity.

SUBMITTER: Li N

PROVIDER: S-EPMC8343000 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Globally, the most common form of arrhythmias is atrial fibrillation (AF), which causes severe morbidity, mortality, and socioeconomic burden. The application of machine learning algorithms in combination with weighted gene co-expression network analysis (WGCNA) can be used to screen genes, therefore, we aimed to screen for potential biomarkers associated with AF development using this integrated bioinformatics approach. Methods: On the basis of the AF endocardium gene expression profiles GSE79768 and GSE115574 from the Gene Expression Omnibus database, differentially expressed genes (DEGs) between AF and sinus rhythm samples were identified. DEGs enrichment analysis and transcription factor screening were then performed. Hub genes for AF were screened using WGCNA and machine learning algorithms, and the diagnostic accuracy was assessed by the receiver operating characteristic (ROC) curves. GSE41177 was used as the validation set for verification. Subsequently, we identified the specific signaling pathways in which the key biomarkers were involved, using gene set enrichment analysis and reverse prediction of mRNA-miRNA interaction pairs. Finally, we explored the associations between the hub genes and immune microenvironment and immune regulation. Results: Fifty-seven DEGs were identified, and the two hub genes, hypoxia inducible factor 1 subunit alpha inhibitor (HIF1AN) and mitochondrial inner membrane protein MPV17 (MPV17), were screened using WGCNA combined with machine learning algorithms. The areas under the receiver operating characteristic curves for MPV17 and HIF1AN validated that two genes predicted AF development, and the differential expression of the hub genes was verified in the external validation dataset. Enrichment analysis showed that MPV17 and HIF1AN affect mitochondrial dysfunction, oxidative stress, gap junctions, and other signaling pathway functions. Immune cell infiltration and immunomodulatory correlation analyses showed that MPV17 and HIF1AN are strongly correlated with the content of immune cells and significantly correlated with HLA expression. Conclusion: The identification of hub genes associated with AF using WGCNA combined with machine learning algorithms and their correlation with immune cells and immune gene expression can elucidate the molecular mechanisms underlying AF occurrence. This may further identify more accurate and effective biomarkers and therapeutic targets for the diagnosis and treatment of AF.

Project description:Background: Using machine learning based on metabolomics, this study aimed to construct an effective primary Sjogren's syndrome (pSS) diagnostics model and reveal the potential targets and biomarkers of pSS. Methods: From a total of 39 patients with pSS and 38 healthy controls (HCs), serum specimens were collected. The samples were analyzed by ultra-high-performance liquid chromatography coupled with high-resolution mass spectrometry. Three machine learning algorithms, including the least absolute shrinkage and selection operator (LASSO), random forest (RF), and extreme gradient boosting (XGBoost), were used to build the pSS diagnosis models. Afterward, four machine learning methods were used to reduce the dimensionality of the metabolomics data. Finally, metabolites with significant differences were screened and pathway analysis was conducted. Results: The area under the curve (AUC), sensitivity, and specificity of LASSO, RF and XGBoost test set all reached 1.00. Orthogonal partial least squares discriminant analysis was used to classify the metabolomics data. By combining the results of the univariate false discovery rate and the importance of the variable in projection, we identified 21 significantly different metabolites. Using these 21 metabolites for diagnostic modeling, the AUC, sensitivity, and specificity of LASSO, RF, and XGBoost all reached 1.00. Metabolic pathway analysis revealed that these 21 metabolites are highly correlated with amino acid and lipid metabolisms. On the basis of 21 metabolites, we screened the important variables in the models. Further, five common variables were obtained by intersecting the important variables of three models. Based on these five common variables, the AUC, sensitivity, and specificity of LASSO, RF, and XGBoost all reached 1.00.2-Hydroxypalmitic acid, L-carnitine and cyclic AMP were found to be potential targets and specific biomarkers for pSS. Conclusion: The combination of machine learning and metabolomics can accurately distinguish between patients with pSS and HCs. 2-Hydroxypalmitic acid, L-carnitine and cyclic AMP were potential targets and biomarkers for pSS.

Project description:Objectives The aim of this study was to investigate serum biomarkers linked to primary Sjögren's syndrome (pSS)-associated interstitial lung disease (ILD). Methods 69 pSS patients were consecutively enrolled and evaluated via quantitative ILD scoring based on high-resolution computed tomography (HRCT). Biomarkers of interest were assessed by multiplex enzyme-linked immunosorbent assays (ELISAs). Results Among consecutively enrolled patients with pSS, the presence of pSS–ILD was 50% based on the presence of radiographically defined interstitial lung abnormalities (ILA) meeting specified criteria for mild/moderate (ILA 2) or severe (ILA 3) disease. Age, immunoglobulin M (IgM), C-reactive protein (CRP), and serum levels of eotaxin/CCL11, Krebs von den Lungen-6 (KL-6), TNFα, and TGFα were significantly higher in the combined pSS–ILD group (ILA 2 + ILA 3) than in the pSS–no-ILD and pSS–indeterminate ILD groups (ILA 0 and ILA 1, respectively) in unadjusted analyses (p < 0.05 for all variables). A binary logistic regression model revealed that disease duration and KL-6 levels were associated with the presence of pSS–ILD (p < 0.05). Complementary least absolute shrinkage and selection operator (LASSO) modeling showed that age, KL-6, and TNF-α effectively differentiated pSS–ILD (ILA 2 + ILA3) from pSS without ILD (ILA 0 + ILA 1), with an area under the curve (AUC) of 0.883 (p value < 0.0001). Conclusions Patient age, disease duration, and serum levels of both KL-6 and TNFα were the most discriminating factors associated with the presence of ILD in our pSS patients. Higher levels of CRP, IgM, eotaxin, TGFα, and TNFα should also prompt the search for occult as well as clinically evident lung involvement based on statistically significant univariate associations with pSS–ILD. Clinical trial registration None. Supplementary Information The online version contains supplementary material available at 10.1186/s40001-022-00828-3.

Project description:BackgroundReduced salivation is considered a major clinical feature of most but not all cases of primary Sjögren's syndrome (pSS). Reduced saliva flow may lead to changes in the salivary microbiota. These changes have mainly been studied with culture that typically recovers only 65% of the bacteria present.ObjectiveThis study was to use high throughput sequencing, covering both cultivated and not-yet-cultivated bacteria, to assess the bacterial microbiota of whole saliva in pSS patients with normal salivation.MethodsBacteria of whole unstimulated saliva from nine pSS patients with normal salivation flow and from nine healthy controls were examined by high throughput sequencing of the hypervariable region V1V2 of 16S rRNA using the 454 GS Junior system. Raw sequence reads were subjected to a species-level, reference-based taxonomy assignment pipeline specially designed for studying the human oral microbial community. Each of the sequence reads was BLASTN-searched against a database consisting of reference sequences representing 1,156 oral and 12,013 non-oral species. Unassigned reads were then screened for high-quality non-chimeras and subjected to de novo species-level operational taxonomy unit (OTU) calling for potential novel species. Downstream analyses, including alpha and beta diversities, were analyzed using the Quantitative Insights into Microbial Ecology (QIIME) pipeline. To reveal significant differences between the microbiota of control saliva and Sjögren's saliva, a statistical method introduced in Metastats www.metastats.cbcb.umd.edu was used.ResultsSaliva of pSS patients with normal salivation had a significantly higher frequency of Firmicutes compared with controls (p=0.004). Two other major phyla, Synergistetes and Spirochaetes, were significantly depleted in pSS (p=0.001 for both). In addition, we saw a nearly 17% decrease in the number of genera in pSS (25 vs. 30). While Prevotella was almost equally abundant in both groups (25% in pSS and 22% in controls), about a twofold increase in pSS of Streptococcus (28% vs. 17%) and Veillonella (26% vs. 12%) was detected. Prevotella melaninogenica was the major species in controls (13%) while Veillonella atypica and the Veillonella parvula groups dominated in patient samples (14 and 14%). The scarcity in bacterial species in pSS compared with controls was also demonstrated by alpha and beta diversity analyses, as well as read abundance depicted in a phylogenetic tree.ConclusionsWhile Firmicutes was significantly higher in pSS patients than in controls, Synergistetes and Spirochaetes were significantly lower. The number of bacterial genera and species was also lower. These data showed that microbial dysbiosis is another key characteristic of pSS whole saliva which can occur independent of hyposalivation.

Project description:Epithelial ovarian cancer (EOC) is one of the malignancies in women, which has the highest mortality. However, the microlevel mechanism has not been discussed in detail. The expression profiles GSE27651, GSE38666, GSE40595, and GSE66957 including 188 tumor and 52 nontumor samples were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were filtered using R software, and we performed functional analysis using the clusterProfiler. Cytoscape software, the molecular complex detection plugin and database STRING analyzed DEGs to construct protein-protein interaction network. We identified 116 DEGs including 81 upregulated and 35 downregulated DEGs. Functional analysis revealed that they were significantly enriched in the extracellular region and biosynthesis of amino acids. We next identified four bioactive compounds (vorinostat, LY-294002,trichostatin A, and tanespimycin) based on ConnectivityMap. Then 114 nodes were obtained in protein-protein interaction. The three most relevant modules were detected. In addition, according to degree ≥ 10, 14 core genes including FOXM1, CXCR4, KPNA2, NANOG, UBE2C, KIF11, ZWINT, CDCA5, DLGAP5, KIF15, MCM2, MELK, SPP1, and TRIP13 were identified. Kaplan-Meier analysis, Oncomine, and Gene Expression Profiling Interactive Analysis showed that overexpression of FOXM1, SPP1, UBE2C, KIF11, ZWINT, CDCA5, UBE2C, and KIF15 was related to bad prognosis of EOC patients. CDCA5, FOXM1, KIF15, MCM2, and ZWINT were associated with stage. Receiver operating characteristic (ROC) curve showed that messenger RNA levels of these five genes exhibited better diagnostic efficiency for normal and tumor tissues. The Human Protein Atlas database was performed. The protein levels of these five genes were significantly higher in tumor tissues compared with normal tissues. Functional enrichment analysis suggested that all the hub genes played crucial roles in citrate cycle tricarboxylic acid cycle. Furthermore, the univariate and multivariate Cox proportional hazards regression showed that ZWINT was independent prognostic indictor among EOC patients. The genes and pathways discovered in the above studies may open a new direction for EOC treatment.