ANXA5: A Key Regulator of Immune Cell Infiltration in Hepatocellular Carcinoma.

BACKGROUND Hepatocellular carcinoma (HCC) poses a significant threat to human life and is the most prevalent form of liver cancer. The intricate interplay between apoptosis, a common form of programmed cell death, and its role in immune regulation stands as a crucial mechanism influencing tumor metastasis. MATERIAL AND METHODS Utilizing HCC samples from the TCGA database and 61 anoikis-related genes (ARGs) sourced from GeneCards, we analyzed the relationship between ARGs and immune cell infiltration in HCC. Subsequently, we identified long non-coding RNAs (lncRNAs) associated with ARGs, using the least absolute shrinkage and selection operator (LASSO) regression analysis to construct a robust prognostic model. The predictive capabilities of the model were then validated through examination in a single-cell dataset. RESULTS Our constructed prognostic model, derived from lncRNAs linked to ARGs, comprised 11 significant lncRNAs: NRAV, MCM3AP-AS1, OTUD6B-AS1, AC026356.1, AC009133.1, DDX11-AS1, AC108463.2, MIR4435-2HG, WARS2-AS1, LINC01094, and HCG18. The risk score assigned to HCC samples demonstrated associations with immune indicators and the infiltration of immune cells. Further, we identified Annexin A5 (ANXA5) as the pivotal gene among ARGs, with it exerting a prominent role in regulating the lncRNA gene signature. Our validation in a single-cell database elucidated the involvement of ANXA5 in immune cell infiltration, specifically in the regulation of mononuclear cells. CONCLUSIONS This study delves into the intricate correlation between ARGs and immune cell infiltration in HCC, culminating in the development of a novel prognostic model reliant on 11 ARGs-associated lncRNAs. Furthermore, our findings highlight ANXA5 as a promising target for immune regulation in HCC, offering new perspectives for immune therapy in the context of HCC.

Identification of shared gene signatures and pathways for diagnosing osteoporosis with sarcopenia through integrated bioinformatics analysis and machine learning.

BACKGROUND: Prior studies have suggested a potential relationship between osteoporosis and sarcopenia, both of which can present symptoms of compromised mobility. Additionally, fractures among the elderly are often considered a common outcome of both conditions. There is a strong correlation between fractures in the elderly population, decreased muscle mass, weakened muscle strength, heightened risk of falls, and diminished bone density. This study aimed to pinpoint crucial diagnostic candidate genes for osteoporosis patients with concomitant sarcopenia. METHODS: Two osteoporosis datasets and one sarcopenia dataset were obtained from the Gene Expression Omnibus (GEO). Differential expression genes (DEGs) and module genes were identified using Limma and Weighted Gene Co-expression Network Analysis (WGCNA), followed by functional enrichment analysis, construction of protein-protein interaction (PPI) networks, and application of a machine learning algorithm (least absolute shrinkage and selection operator (LASSO) regression) to determine candidate hub genes for diagnosing osteoporosis combined with sarcopenia. Receiver operating characteristic (ROC) curves and column line plots were generated. RESULTS: The merged osteoporosis dataset comprised 2067 DEGs, with 424 module genes filtered in sarcopenia. The intersection of DEGs between osteoporosis and sarcopenia module genes consisted of 60 genes, primarily enriched in viral infection. Through construction of the PPI network, 30 node genes were filtered, and after machine learning, 7 candidate hub genes were selected for column line plot construction and diagnostic value assessment. Both the column line plots and all 7 candidate hub genes exhibited high diagnostic value (area under the curve ranging from 1.00 to 0.93). CONCLUSION: We identified 7 candidate hub genes (PDP1, ALS2CL, VLDLR, PLEKHA6, PPP1CB, MOSPD2, METTL9) and constructed column line plots for osteoporosis combined with sarcopenia. This study provides reference for potential peripheral blood diagnostic candidate genes for sarcopenia in osteoporosis patients.

Sample-Wise and Gene-Wise Comparisons Confirm a Greater Similarity of RNA and Protein Expression Data at the Level of Molecular Pathways and Suggest an Approach for the Data Quality Check in High-Throughput Expression Databases.

Identification of genes and molecular pathways with congruent profiles in the proteomic and transcriptomic datasets may result in the discovery of promising transcriptomic biomarkers that would be more relevant to phenotypic changes. In this study, we conducted comparative analysis of 943 paired RNA and proteomic profiles obtained for the same samples of seven human cancer types from The Cancer Genome Atlas (TCGA) and NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) [two major open human cancer proteomic and transcriptomic databases] that included 15,112 protein-coding genes and 1611 molecular pathways. Overall, our findings demonstrated statistically significant improvement of the congruence between RNA and proteomic profiles when performing analysis at the level of molecular pathways rather than at the level of individual gene products. Transition to the molecular pathway level of data analysis increased the correlation to 0.19-0.57 (Pearson) and 0.14-057 (Spearman), or 2-3-fold for some cancer types. Evaluating the gain of the correlation upon transition to the data analysis the pathway level can be used to refine the omics data by identifying outliers that can be excluded from the comparison of RNA and proteomic profiles. We suggest using sample- and gene-wise correlations for individual genes and molecular pathways as a measure of quality of RNA/protein paired molecular data. We also provide a database of human genes, molecular pathways, and samples related to the correlation between RNA and protein products to facilitate an exploration of new cancer transcriptomic biomarkers and molecular mechanisms at different levels of human gene expression.

Serum 25-hydroxyvitamin D levels and dermatomyositis: A 2-sample mendelian randomization study.

BACKGROUND: Previous studies have reported low serum 25-hydroxyvitamin D [25(OH)D] levels in dermatomyositis (DM) patients, but the exact causal relationship between them remains elusive. Our aim is to confirm the causal relationship between 25(OH)D and DM risk through a Mendelian randomization study. METHODS: Retrieve genome-wide association study (GWAS) data on 25(OH)D (n = 441 291) and DM (n cases = 201, n controls = 172 834) from the GWAS database (https://gwas.mrcieu.ac.uk/). Select single-nucleotide polymorphisms (SNPs) strongly correlated with 25(OH)D as instrumental variables (IVs). The primary analytical approach involves the use of the inverse-variance weighted method (IVW), supplemented by MR-Egger regression and weighted median methods to enhance the reliability of the results. Heterogeneity and sensitivity analyses were conducted using Cochran's Q and leave-one-out approaches, respectively. RESULTS: The IVW analysis confirmed a positive causal relationship between genetic variation in 25(OH)D levels and DM (OR = 2.36, 95% CI = 1.01-5.52, p = .048). Although not statistically significant (all p > .05), the other methods also suggested a protective effect of 25(OH)D on DM. Based on MR-Egger intercepts and Cochran's Q analysis, the selected SNPs showed no horizontal pleiotropy and heterogeneity. Sensitivity analysis demonstrated the robustness of the results against individual SNPs. CONCLUSION: We provide the first evidence of a causal relationship between 25(OH)D levels and DM. Our findings support the importance of measuring serum 25(OH)D levels and considering vitamin D supplementation in clinical practice for patients with DM.

An atlas and database of neuropeptide gene expression in the adult zebrafish forebrain.

Zebrafish is a useful model organism in neuroscience; however, its gene expression atlas in the adult brain is not well developed. In the present study, we examined the expression of 38 neuropeptides, comparing with GABAergic and glutamatergic neuron marker genes in the adult zebrafish brain by comprehensive in situ hybridization. The results are summarized as an expression atlas in 19 coronal planes of the forebrain. Furthermore, the scanned data of all brain sections were made publicly available in the Adult Zebrafish Brain Gene Expression Database (https://ssbd.riken.jp/azebex/). Based on these data, we performed detailed comparative neuroanatomical analyses of the hypothalamus and found that several regions previously described as one nucleus in the reference zebrafish brain atlas contain two or more subregions with significantly different neuropeptide/neurotransmitter expression profiles. Subsequently, we compared the expression data in zebrafish telencephalon and hypothalamus obtained in this study with those in mice, by performing a cluster analysis. As a result, several nuclei in zebrafish and mice were clustered in close vicinity. The present expression atlas, database, and anatomical findings will contribute to future neuroscience research using zebrafish.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research.

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

Bioinformatics analysis of signature genes related to cell death in keratoconus.

Keratoconus is corneal disease in which the progression of conical dilation of cornea leads to reduced visual acuity and even corneal perforation. However, the etiology mechanism of keratoconus is still unclear. This study aims to identify the signature genes related to cell death in keratoconus and examine the function of these genes. A dataset of keratoconus from the GEO database was analysed to identify the differentially expressed genes (DEGs). A total of 3558 DEGs were screened from GSE151631. The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that they mainly involved in response to hypoxia, cell-cell adhesion, and IL-17 signaling pathway. Then, the cell death-related genes datasets were intersected with the above 3558 DEGs to obtain 70 ferroptosis-related DEGs (FDEGs), 32 autophagy-related DEGs (ADEGs), six pyroptosis-related DEGs (PDEGs), four disulfidptosis-related DEGs (DDEGs), and one cuproptosis-related DEGs (CDEGs). After using Least absolute shrinkage and selection operator (LASSO), Random Forest analysis, and receiver operating characteristic (ROC) curve analysis, one ferroptosis-related gene (TNFAIP3) and five autophagy-related genes (CDKN1A, HSPA5, MAPK8IP1, PPP1R15A, and VEGFA) were screened out. The expressions of the above six genes were significantly decreased in keratoconus and the area under the curve (AUC) values of these genes was 0.944, 0.893, 0.797, 0.726, 0.882 and 0.779 respectively. GSEA analysis showed that the above six genes mainly play an important role in allograft rejection, asthma, and circadian rhythm etc. In conclusion, the results of this study suggested that focusing on these genes and autoimmune diseases will be a beneficial perspective for the keratoconus etiology research.

A Comprehensive Prognostic Model for Colon Adenocarcinoma Depending on Nuclear-Mitochondrial-Related Genes.

Background: Colon adenocarcinoma (COAD) has increasing incidence and is one of the most common malignant tumors. The mitochondria involved in cell energy metabolism, oxygen free radical generation, and cell apoptosis play important roles in tumorigenesis and progression. The relationship between mitochondrial genes and COAD remains largely unknown. Methods: COAD data including 512 samples were set out from the UCSC Xena database. The nuclear mitochondrial-related genes (NMRGs)-related risk prognostic model and prognostic nomogram were constructed, and NMRGs-related gene mutation and the immune environment were analyzed using bioinformatics methods. Then, a liver metastasis model of colorectal cancer was constructed and protein expression was detected using Western blot assay. Results: A prognostic model for COAD was constructed. Comparing the prognostic model dataset and the validation dataset showed considerable correlation in both risk grouping and prognosis. Based on the risk score (RS) model, the samples of the prognostic dataset were divided into high risk group and low risk group. Moreover, pathologic N and T stage and tumor recurrence in the two risk groups were significantly different. The four prognostic factors, including age and pathologic T stage in the nomogram survival model also showed excellent predictive performance. An optimal combination of nine differentially expressed NMRGs was finally obtained, including LARS2, PARS2, ETHE1, LRPPRC, TMEM70, AARS2, ACAD9, VARS2, and ATP8A2. The high-RS group had more inflamed immune features, including T and CD4+ memory cell activation. Besides, mitochondria-associated LRPPRC and LARS2 expression levels were increased in vivo xenograft construction and liver metastases assays. Conclusion: This study established a comprehensive prognostic model for COAD, incorporating nine genes associated with nuclear-mitochondrial functions. This model demonstrates superior predictive performance across four prognostic factors: age, pathological T stage, tumor recurrence, and overall prognosis. It is anticipated to be an effective model for enhancing the prognosis and treatment of COAD.

RicePilaf: a post-GWAS/QTL dashboard to integrate pangenomic, coexpression, regulatory, epigenomic, ontology, pathway, and text-mining information to provide functional insights into rice QTLs and GWAS loci.

BACKGROUND: As the number of genome-wide association study (GWAS) and quantitative trait locus (QTL) mappings in rice continues to grow, so does the already long list of genomic loci associated with important agronomic traits. Typically, loci implicated by GWAS/QTL analysis contain tens to hundreds to thousands of single-nucleotide polmorphisms (SNPs)/genes, not all of which are causal and many of which are in noncoding regions. Unraveling the biological mechanisms that tie the GWAS regions and QTLs to the trait of interest is challenging, especially since it requires collating functional genomics information about the loci from multiple, disparate data sources. RESULTS: We present RicePilaf, a web app for post-GWAS/QTL analysis, that performs a slew of novel bioinformatics analyses to cross-reference GWAS results and QTL mappings with a host of publicly available rice databases. In particular, it integrates (i) pangenomic information from high-quality genome builds of multiple rice varieties, (ii) coexpression information from genome-scale coexpression networks, (iii) ontology and pathway information, (iv) regulatory information from rice transcription factor databases, (v) epigenomic information from multiple high-throughput epigenetic experiments, and (vi) text-mining information extracted from scientific abstracts linking genes and traits. We demonstrate the utility of RicePilaf by applying it to analyze GWAS peaks of preharvest sprouting and genes underlying yield-under-drought QTLs. CONCLUSIONS: RicePilaf enables rice scientists and breeders to shed functional light on their GWAS regions and QTLs, and it provides them with a means to prioritize SNPs/genes for further experiments. The source code, a Docker image, and a demo version of RicePilaf are publicly available at https://github.com/bioinfodlsu/rice-pilaf.

Identification and validation of potential diagnostic signature and immune cell infiltration for HIRI based on cuproptosis-related genes through bioinformatics analysis and machine learning.

Background and aims: Cuproptosis has emerged as a significant contributor in the progression of various diseases. This study aimed to assess the potential impact of cuproptosis-related genes (CRGs) on the development of hepatic ischemia and reperfusion injury (HIRI). Methods: The datasets related to HIRI were sourced from the Gene Expression Omnibus database. The comparative analysis of differential gene expression involving CRGs was performed between HIRI and normal liver samples. Correlation analysis, function enrichment analyses, and protein-protein interactions were employed to understand the interactions and roles of these genes. Machine learning techniques were used to identify hub genes. Additionally, differences in immune cell infiltration between HIRI patients and controls were analyzed. Quantitative real-time PCR and western blotting were used to verify the expression of the hub genes. Results: Seventy-five HIRI and 80 control samples from three databases were included in the bioinformatics analysis. Three hub CRGs (NLRP3, ATP7B and NFE2L2) were identified using three machine learning models. Diagnostic accuracy was assessed using a receiver operating characteristic (ROC) curve for the hub genes, which yielded an area under the ROC curve (AUC) of 0.832. Remarkably, in the validation datasets GSE15480 and GSE228782, the three hub genes had AUC reached 0.904. Additional analyses, including nomograms, decision curves, and calibration curves, supported their predictive power for diagnosis. Enrichment analyses indicated the involvement of these genes in multiple pathways associated with HIRI progression. Comparative assessments using CIBERSORT and gene set enrichment analysis suggested elevated expression of these hub genes in activated dendritic cells, neutrophils, activated CD4 memory T cells, and activated mast cells in HIRI samples versus controls. A ceRNA network underscored a complex regulatory interplay among genes. The genes mRNA and protein levels were also verified in HIRI-affected mouse liver tissues. Conclusion: Our findings have provided a comprehensive understanding of the association between cuproptosis and HIRI, establishing a promising diagnostic pattern and identifying latent therapeutic targets for HIRI treatment. Additionally, our study offers novel insights to delve deeper into the underlying mechanisms of HIRI.

PESSA: A web tool for pathway enrichment score-based survival analysis in cancer.

The activation levels of biologically significant gene sets are emerging tumor molecular markers and play an irreplaceable role in the tumor research field; however, web-based tools for prognostic analyses using it as a tumor molecular marker remain scarce. We developed a web-based tool PESSA for survival analysis using gene set activation levels. All data analyses were implemented via R. Activation levels of The Molecular Signatures Database (MSigDB) gene sets were assessed using the single sample gene set enrichment analysis (ssGSEA) method based on data from the Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA), The European Genome-phenome Archive (EGA) and supplementary tables of articles. PESSA was used to perform median and optimal cut-off dichotomous grouping of ssGSEA scores for each dataset, relying on the survival and survminer packages for survival analysis and visualisation. PESSA is an open-access web tool for visualizing the results of tumor prognostic analyses using gene set activation levels. A total of 238 datasets from the GEO, TCGA, EGA, and supplementary tables of articles; covering 51 cancer types and 13 survival outcome types; and 13,434 tumor-related gene sets are obtained from MSigDB for pre-grouping. Users can obtain the results, including Kaplan-Meier analyses based on the median and optimal cut-off values and accompanying visualization plots and the Cox regression analyses of dichotomous and continuous variables, by selecting the gene set markers of interest. PESSA (https://smuonco.shinyapps.io/PESSA/ OR http://robinl-lab.com/PESSA) is a large-scale web-based tumor survival analysis tool covering a large amount of data that creatively uses predefined gene set activation levels as molecular markers of tumors.

Phenome-wide identification of therapeutic genetic targets, leveraging knowledge graphs, graph neural networks, and UK Biobank data.

The ongoing expansion of human genomic datasets propels therapeutic target identification; however, extracting gene-disease associations from gene annotations remains challenging. Here, we introduce Mantis-ML 2.0, a framework integrating AstraZeneca's Biological Insights Knowledge Graph and numerous tabular datasets, to assess gene-disease probabilities throughout the phenome. We use graph neural networks, capturing the graph's holistic structure, and train them on hundreds of balanced datasets via a robust semi-supervised learning framework to provide gene-disease probabilities across the human exome. Mantis-ML 2.0 incorporates natural language processing to automate disease-relevant feature selection for thousands of diseases. The enhanced models demonstrate a 6.9% average classification power boost, achieving a median receiver operating characteristic (ROC) area under curve (AUC) score of 0.90 across 5220 diseases from Human Phenotype Ontology, OpenTargets, and Genomics England. Notably, Mantis-ML 2.0 prioritizes associations from an independent UK Biobank phenome-wide association study (PheWAS), providing a stronger form of triaging and mitigating against underpowered PheWAS associations. Results are exposed through an interactive web resource.

Machine learning and weighted gene co-expression network analysis identify a three-gene signature to diagnose rheumatoid arthritis.

Background: Rheumatoid arthritis (RA) is a systemic immune-related disease characterized by synovial inflammation and destruction of joint cartilage. The pathogenesis of RA remains unclear, and diagnostic markers with high sensitivity and specificity are needed urgently. This study aims to identify potential biomarkers in the synovium for diagnosing RA and to investigate their association with immune infiltration. Methods: We downloaded four datasets containing 51 RA and 36 healthy synovium samples from the Gene Expression Omnibus database. Differentially expressed genes were identified using R. Then, various enrichment analyses were conducted. Subsequently, weighted gene co-expression network analysis (WGCNA), random forest (RF), support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage and selection operator (LASSO) were used to identify the hub genes for RA diagnosis. Receiver operating characteristic curves and nomogram models were used to validate the specificity and sensitivity of hub genes. Additionally, we analyzed the infiltration levels of 28 immune cells in the expression profile and their relationship with the hub genes using single-sample gene set enrichment analysis. Results: Three hub genes, namely, ribonucleotide reductase regulatory subunit M2 (RRM2), DLG-associated protein 5 (DLGAP5), and kinesin family member 11 (KIF11), were identified through WGCNA, LASSO, SVM-RFE, and RF algorithms. These hub genes correlated strongly with T cells, natural killer cells, and macrophage cells as indicated by immune cell infiltration analysis. Conclusion: RRM2, DLGAP5, and KIF11 could serve as potential diagnostic indicators and treatment targets for RA. The infiltration of immune cells offers additional insights into the underlying mechanisms involved in the progression of RA.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Exploring the common mechanism of vascular dementia and inflammatory bowel disease: a bioinformatics-based study.

Objective: Emerging evidence has shown that gut diseases can regulate the development and function of the immune, metabolic, and nervous systems through dynamic bidirectional communication on the brain-gut axis. However, the specific mechanism of intestinal diseases and vascular dementia (VD) remains unclear. We designed this study especially, to further clarify the connection between VD and inflammatory bowel disease (IBD) from bioinformatics analyses. Methods: We downloaded Gene expression profiles for VD (GSE122063) and IBD (GSE47908, GSE179285) from the Gene Expression Omnibus (GEO) database. Then individual Gene Set Enrichment Analysis (GSEA) was used to confirm the connection between the two diseases respectively. The common differentially expressed genes (coDEGs) were identified, and the STRING database together with Cytoscape software were used to construct protein-protein interaction (PPI) network and core functional modules. We identified the hub genes by using the Cytohubba plugin. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to identify pathways of coDEGs and hub genes. Subsequently, receiver operating characteristic (ROC) analysis was used to identify the diagnostic ability of these hub genes, and a training dataset was used to verify the expression levels of the hub genes. An alternative single-sample gene set enrichment (ssGSEA) algorithm was used to analyze immune cell infiltration between coDEGs and immune cells. Finally, the correlation between hub genes and immune cells was analyzed. Results: We screened 167 coDEGs. The main articles of coDEGs enrichment analysis focused on immune function. 8 shared hub genes were identified, including PTPRC, ITGB2, CYBB, IL1B, TLR2, CASP1, IL10RA, and BTK. The functional categories of hub genes enrichment analysis were mainly involved in the regulation of immune function and neuroinflammatory response. Compared to the healthy controls, abnormal infiltration of immune cells was found in VD and IBD. We also found the correlation between 8 shared hub genes and immune cells. Conclusions: This study suggests that IBD may be a new risk factor for VD. The 8 hub genes may predict the IBD complicated with VD. Immune-related coDEGS may be related to their association, which requires further research to prove.

Identification of potential biomarkers of leprosy: A study based on GEO datasets.

Leprosy has a high rate of cripplehood and lacks available early effective diagnosis methods for prevention and treatment, thus novel effective molecule markers are urgently required. In this study, we conducted bioinformatics analysis with leprosy and normal samples acquired from the GEO database(GSE84893, GSE74481, GSE17763, GSE16844 and GSE443). Through WGCNA analysis, 85 hub genes were screened(GS > 0.7 and MM > 0.8). Through DEG analysis, 82 up-regulated and 3 down-regulated genes were screened(|Log2FC| > 3 and FDR < 0.05). Then 49 intersection genes were considered as crucial and subjected to GO annotation, KEGG pathway and PPI analysis to determine the biological significance in the pathogenesis of leprosy. Finally, we identified a gene-pathway network, suggesting ITK, CD48, IL2RG, CCR5, FGR, JAK3, STAT1, LCK, PTPRC, CXCR4 can be used as biomarkers and these genes are active in 6 immune system pathways, including Chemokine signaling pathway, Th1 and Th2 cell differentiation, Th17 cell differentiation, T cell receptor signaling pathway, Natural killer cell mediated cytotoxicity and Leukocyte transendothelial migration. We identified 10 crucial gene markers and related important pathways that acted as essential components in the etiology of leprosy. Our study provides potential targets for diagnostic biomarkers and therapy of leprosy.

Unlocking the genetic tapestry of autoimmune diseases: Unveiling common genes across multiple conditions.

OBJECTIVES: This study aimed to unravel the complexities of autoimmune diseases by conducting a comprehensive analysis of gene expression data across 10 conditions, including systemic lupus erythematosus (SLE), psoriasis, Sjögren's syndrome, sclerosis, immune-associated diseases, osteoarthritis, cystic fibrosis, inflammatory bowel disease (IBD), type 1 diabetes, and Guillain-Barré syndrome. METHODS: Gene expression profiles were rigorously examined to identify both upregulated and downregulated genes specific to each autoimmune disease. The study employed visual representation techniques such as heatmaps, volcano plots, and contour-MA plots to provide an intuitive understanding of the complex gene expression patterns in these conditions. RESULTS: Distinct gene expression profiles for each autoimmune condition were uncovered, with psoriasis and osteoarthritis standing out due to a multitude of both upregulated and downregulated genes, indicating intricate molecular interplays in these disorders. Notably, common upregulated and downregulated genes were identified across various autoimmune conditions, with genes like SELENBP1, MMP9, BNC1, and COL1A1 emerging as pivotal players. CONCLUSION: This research contributes valuable insights into the molecular signatures of autoimmune diseases, highlighting the unique gene expression patterns characterizing each condition. The identification of common genes shared among different autoimmune conditions, and their potential role in mitigating the risk of rare diseases in patients with more prevalent conditions, underscores the growing significance of genetics in healthcare and the promising future of personalized medicine.

OMD Curation Toolkit: a workflow for in-house curation of public omics datasets.

BACKGROUND: Major advances in sequencing technologies and the sharing of data and metadata in science have resulted in a wealth of publicly available datasets. However, working with and especially curating public omics datasets remains challenging despite these efforts. While a growing number of initiatives aim to re-use previous results, these present limitations that often lead to the need for further in-house curation and processing. RESULTS: Here, we present the Omics Dataset Curation Toolkit (OMD Curation Toolkit), a python3 package designed to accompany and guide the researcher during the curation process of metadata and fastq files of public omics datasets. This workflow provides a standardized framework with multiple capabilities (collection, control check, treatment and integration) to facilitate the arduous task of curating public sequencing data projects. While centered on the European Nucleotide Archive (ENA), the majority of the provided tools are generic and can be used to curate datasets from different sources. CONCLUSIONS: Thus, it offers valuable tools for the in-house curation previously needed to re-use public omics data. Due to its workflow structure and capabilities, it can be easily used and benefit investigators in developing novel omics meta-analyses based on sequencing data.

Identification of Breast Cancer Subtypes Based on Endoplasmic Reticulum Stress-Related Genes and Analysis of Prognosis and Immune Microenvironment in Breast Cancer Patients.

Introduction: Endoplasmic reticulum stress (ERS) was a response to the accumulation of unfolded proteins and plays a crucial role in the development of tumors, including processes such as tumor cell invasion, metastasis, and immune evasion. However, the specific regulatory mechanisms of ERS in breast cancer (BC) remain unclear. Methods: In this study, we analyzed RNA sequencing data from The Cancer Genome Atlas (TCGA) for breast cancer and identified 8 core genes associated with ERS: ELOVL2, IFNG, MAP2K6, MZB1, PCSK6, PCSK9, IGF2BP1, and POP1. We evaluated their individual expression, independent diagnostic, and prognostic values in breast cancer patients. A multifactorial Cox analysis established a risk prognostic model, validated with an external dataset. Additionally, we conducted a comprehensive assessment of immune infiltration and drug sensitivity for these genes. Results: The results indicate that these eight core genes play a crucial role in regulating the immune microenvironment of breast cancer (BRCA) patients. Meanwhile, an independent diagnostic model based on the expression of these eight genes shows limited independent diagnostic value, and its independent prognostic value is unsatisfactory, with the time ROC AUC values generally below 0.5. According to the results of logistic regression neural networks and risk prognosis models, when these eight genes interact synergistically, they can serve as excellent biomarkers for the diagnosis and prognosis of breast cancer patients. Furthermore, the research findings have been confirmed through qPCR experiments and validation. Conclusion: In conclusion, we explored the mechanisms of ERS in BRCA patients and identified 8 outstanding biomolecular diagnostic markers and prognostic indicators. The research results were double-validated using the GEO database and qPCR.