irGSEA: the integration of single-cell rank-based gene set enrichment analysis.

irGSEA is an R package designed to assess the outcomes of various gene set scoring methods when applied to single-cell RNA sequencing data. This package incorporates six distinct scoring methods that rely on the expression ranks of genes, emphasizing relative expression levels over absolute values. The implemented methods include AUCell, UCell, singscore, ssGSEA, JASMINE and Viper. Previous studies have demonstrated the robustness of these methods to variations in dataset size and composition, generating enrichment scores based solely on the relative gene expression of individual cells. By employing the robust rank aggregation algorithm, irGSEA amalgamates results from all six methods to ascertain the statistical significance of target gene sets across diverse scoring methods. The package prioritizes user-friendliness, allowing direct input of expression matrices or seamless interaction with Seurat objects. Furthermore, it facilitates a comprehensive visualization of results. The irGSEA package and its accompanying documentation are accessible on GitHub (https://github.com/chuiqin/irGSEA).

Benchmarking enrichment analysis methods with the disease pathway network.

Enrichment analysis (EA) is a common approach to gain functional insights from genome-scale experiments. As a consequence, a large number of EA methods have been developed, yet it is unclear from previous studies which method is the best for a given dataset. The main issues with previous benchmarks include the complexity of correctly assigning true pathways to a test dataset, and lack of generality of the evaluation metrics, for which the rank of a single target pathway is commonly used. We here provide a generalized EA benchmark and apply it to the most widely used EA methods, representing all four categories of current approaches. The benchmark employs a new set of 82 curated gene expression datasets from DNA microarray and RNA-Seq experiments for 26 diseases, of which only 13 are cancers. In order to address the shortcomings of the single target pathway approach and to enhance the sensitivity evaluation, we present the Disease Pathway Network, in which related Kyoto Encyclopedia of Genes and Genomes pathways are linked. We introduce a novel approach to evaluate pathway EA by combining sensitivity and specificity to provide a balanced evaluation of EA methods. This approach identifies Network Enrichment Analysis methods as the overall top performers compared with overlap-based methods. By using randomized gene expression datasets, we explore the null hypothesis bias of each method, revealing that most of them produce skewed P-values.

The glioblastoma multiforme tumor site promotes the commitment of tumor-infiltrating lymphocytes to the TH17 lineage in humans.

Although glioblastoma multiforme (GBM) is not an invariably cold tumor, checkpoint inhibition has largely failed in GBM. In order to investigate T cell-intrinsic properties that contribute to the resistance of GBM to endogenous or therapeutically enhanced adaptive immune responses, we sorted CD4+ and CD8+ T cells from the peripheral blood, normal-appearing brain tissue, and tumor bed of nine treatment-naive patients with GBM. Bulk RNA sequencing of highly pure T cell populations from these different compartments was used to obtain deep transcriptomes of tumor-infiltrating T cells (TILs). While the transcriptome of CD8+ TILs suggested that they were partly locked in a dysfunctional state, CD4+ TILs showed a robust commitment to the type 17 T helper cell (TH17) lineage, which was corroborated by flow cytometry in four additional GBM cases. Therefore, our study illustrates that the brain tumor environment in GBM might instruct TH17 commitment of infiltrating T helper cells. Whether these properties of CD4+ TILs facilitate a tumor-promoting milieu and thus could be a target for adjuvant anti-TH17 cell interventions needs to be further investigated.

vissE: a versatile tool to identify and visualise higher-order molecular phenotypes from functional enrichment analysis.

Functional analysis of high throughput experiments using pathway analysis is now ubiquitous. Though powerful, these methods often produce thousands of redundant results owing to knowledgebase redundancies upstream. This scale of results hinders extensive exploration by biologists and can lead to investigator biases due to previous knowledge and expectations. To address this issue, we present vissE, a flexible network-based analysis and visualisation tool that organises information into semantic categories and provides various visualisation modules to characterise them with respect to the underlying data, thus providing a comprehensive view of the biological system. We demonstrate vissE's versatility by applying it to three different technologies: bulk, single-cell and spatial transcriptomics. Applying vissE to a factor analysis of a breast cancer spatial transcriptomic data, we identified stromal phenotypes that support tumour dissemination. Its adaptability allows vissE to enhance all existing gene-set enrichment and pathway analysis workflows, empowering biologists during molecular discovery.

Identification of molecular targets and small drug candidates for Huntington's disease via bioinformatics and a network-based screening approach.

Huntington's disease (HD) is a gradually severe neurodegenerative ailment characterised by an increase of a specific trinucleotide repeat sequence (cytosine-adenine-guanine, CAG). It is passed down as a dominant characteristic that worsens over time, creating a significant risk. Despite being monogenetic, the underlying mechanisms as well as biomarkers remain poorly understood. Furthermore, early detection of HD is challenging, and the available diagnostic procedures have low precision and accuracy. The research was conducted to provide knowledge of the biomarkers, pathways and therapeutic targets involved in the molecular processes of HD using informatic based analysis and applying network-based systems biology approaches. The gene expression profile datasets GSE97100 and GSE74201 relevant to HD were studied. As a consequence, 46 differentially expressed genes (DEGs) were identified. 10 hub genes (TPM1, EIF2S3, CCN2, ACTN1, ACTG2, CCN1, CSRP1, EIF1AX, BEX2 and TCEAL5) were further differentiated in the protein-protein interaction (PPI) network. These hub genes were typically down-regulated. Additionally, DEGs-transcription factors (TFs) connections (e.g. GATA2, YY1 and FOXC1), DEG-microRNA (miRNA) interactions (e.g. hsa-miR-124-3p and has-miR-26b-5p) were also comprehensively forecast. Additionally, related gene ontology concepts (e.g. sequence-specific DNA binding and TF activity) connected to DEGs in HD were identified using gene set enrichment analysis (GSEA). Finally, in silico drug design was employed to find candidate drugs for the treatment HD, and while the possible modest therapeutic compounds (e.g. cortistatin A, 13,16-Epoxy-25-hydroxy-17-cheilanthen-19,25-olide, Hecogenin) against HD were expected. Consequently, the results from this study may give researchers useful resources for the experimental validation of Huntington's diagnosis and therapeutic approaches.

DhuFAP: a platform for gene functional analysis in Dendrobium huoshanense.

BACKGROUND: Dendrobium huoshanense, a traditional medicinal and food plant, has a rich history of use. Recently, its genome was decoded, offering valuable insights into gene function. However, there is no comprehensive gene functional analysis platform for D. huoshanense. RESULT: To address this, we created a platform for gene function analysis and comparison in D. huoshanense (DhuFAP). Using 69 RNA-seq samples, we constructed a gene co-expression network and annotated D. huoshanense genes by aligning sequences with public protein databases. Our platform contained tools like Blast, gene set enrichment analysis, heatmap analysis, sequence extraction, and JBrowse. Analysis revealed co-expression of transcription factors (C2H2, GRAS, NAC) with genes encoding key enzymes in alkaloid biosynthesis. We also showcased the reliability and applicability of our platform using Chalcone synthases (CHS). CONCLUSION: DhuFAP ( www.gzybioinformatics.cn/DhuFAP ) and its suite of tools represent an accessible and invaluable resource for researchers, enabling the exploration of functional information pertaining to D. huoshanense genes. This platform stands poised to facilitate significant biological discoveries in this domain.

Loss of chromosome Y in regulatory T cells.

BACKGROUND: Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. RESULTS: Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. CONCLUSIONS: Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to the increased risk for cancer observed among men with Y loss in leukocytes.

Effect of SNPs on Litter Size in Swine.

Although sows do not directly enter the market, they play an important role in piglet breeding on farms. They consume large amounts of feed, resulting in a significant environmental burden. Pig farms can increase their income and reduce environmental pollution by increasing the litter size (LS) of swine. PCR-RFLP/SSCP and GWAS are common methods to evaluate single-nucleotide polymorphisms (SNPs) in candidate genes. We conducted a systematic meta-analysis of the effect of SNPs on pig LS. We collected and analysed data published over the past 30 years using traditional and network meta-analyses. Trial sequential analysis (TSA) was used to analyse population data. Gene set enrichment analysis and protein-protein interaction network analysis were used to analyse the GWAS dataset. The results showed that the candidate genes were positively correlated with LS, and defects in PCR-RFLP/SSCP affected the reliability of candidate gene results. However, the genotypes with high and low LSs did not have a significant advantage. Current breeding and management practices for sows should consider increasing the LS while reducing lactation length and minimizing the sows' non-pregnancy period as much as possible.

Integrative in silico approaches to analyse microRNA-mediated responses in human diseases.

Advancements in sequencing technologies have facilitated omics level information generation for various diseases in human. High-throughput technologies have become a powerful tool to understand differential expression studies and transcriptional network analysis. An understanding of complex transcriptional networks in human diseases requires integration of datasets representing different RNA species including microRNA (miRNA) and messenger RNA (mRNA). This review emphasises on conceptual explanation of generalized workflow and methodologies to the miRNA mediated responses in human diseases by using different in silico analysis. Although, there have been many prior explorations in miRNA-mediated responses in human diseases, the advantages, limitations and overcoming the limitation through different statistical techniques have not yet been discussed. This review focuses on miRNAs as important gene regulators in human diseases, methodologies for miRNA-target gene prediction and data driven methods for enrichment and network analysis for miRnome-targetome interactions. Additionally, it proposes an integrative workflow to analyse structural components of networks obtained from high-throughput data. This review explains how to apply the existing methods to analyse miRNA-mediated responses in human diseases. It addresses unique characteristics of different analysis, its limitations and its statistical solutions influencing the choice of methods for the analysis through a workflow. Moreover, it provides an overview of promising common integrative approaches to comprehend miRNA-mediated gene regulatory events in biological processes in humans. The proposed methodologies and workflow shall help in the analysis of multi-source data to identify molecular signatures of various human diseases.

CXCL10 and its receptor in patients with chronic hepatitis B and their ability to predict HBeAg seroconversion during antiviral treatment with TDF.

The serum chemokine C-X-C motif ligand-10 (CXCL10) and its unique receptor (CXCR3) may predict the prognosis of patients with chronic hepatitis B (CHB) treated with tenofovir disoproxil fumarate (TDF). Nevertheless, there are few reports on the profile of CXCL10 and CXCR3 and their clinical application in HBeAg (+) CHB patients during TDF antiviral therapy. CXCL10 and CXCR3 were determined in 118 CHB patients naively treated with TDF for at least 96 weeks at baseline and at treatment weeks 12 and 24. In addition, gene set enrichment analysis was used to examine the associated dataset from Gene Expression Omnibus and explore the gene sets associated with HBeAg seroconversion (SC). The change of CXCL10 (ΔCXCL10, baseline to 48-week TDF treatment) and CXCR3 (ΔCXCR3) is closely related to the possibility of HBeAg SC of CHB patients under TDF treatment. Immunohistochemical analysis of CXCL10/CXCR3 protein in liver tissue shows that there is a significant difference between paired liver biopsy samples taken before and after 96 weeks of successful TDF treatment of CHB patients (11 pairs) but no significance for unsuccessful TDF treatment (14 pairs). Multivariate Cox analysis suggests that the ΔCXCL10 is an independent predictive indicator of HBeAg SC, and the area under the receiver operating characteristic curve of the ΔCXCL10 in CHB patients is 0.8867 (p < 0.0001). Our results suggest that a lower descending CXCL10 level is associated with an increased probability of HBeAg SC of CHB patients during TDF therapy. Moreover, liver tissue CXCL10 might be involved in the immunological process of HBeAg SC.

Identifying COVID-19 subtypes by single-sample gene set enrichemnt analysis and providing guidance for sensitive drug selection.

This study aimed at using single-sample gene set enrichment analysis scores to cluster naso/pharyngeal swab specimen samples from coronavirus disease 2019 (COVID-19) patients into two clusters. One cluster with higher fractions of immune cells and more active inflammatory-related pathways was called the Immunity-High (Immunity-H) group, and the other one was called the Immunity-Low group. We explored impacts of the method on COVID-19 treatment. First, given that the Immunity-H group was mainly enriched in inflammatory-related pathways and had higher fractions of inflammatory cells, the Immunity-H group may obtain more curative effects from anti-inflammatory treatment. Second, we searched some hot genes from the PubMed platform that had been studied by researchers and found these genes upregulated in the Immunity-H group, so we speculated the Immunity-H group and Immunity-Low group may have different curative effects from drugs targeting these genes. Finally, we screened out hub genes for the Immunity-H group and predicted potential drugs for these hub genes by a public data set (http://dgidb.genome.wustl.edu). These hub genes are significantly upregulated in the Immunity-H group and neutrophils so that the Immunity-H group may obtain different treatment results from potential drugs compared with the Immunity-Low group. Therefore, the cluster method may provide help in drug development and administration for COVID-19 patients.

LncRNA RPARP-AS1 promotes the progression of osteosarcoma cells through regulating lipid metabolism.

Osteosarcoma (OS) is a highly malignant tumor, and its dysregulated lipid metabolism is associated with tumorigenesis and unfavorable prognosis. Interestingly, long noncoding RNAs (lncRNAs) have emerged as pivotal regulators of lipid metabolism, exerting notable impacts on tumor proliferation. Nevertheless, the involvement of RPARP-AS1, a novel lipid metabolism-associated lncRNA, remains unexplored in the context of OS. This study aims to identify functionally relevant lncRNAs impacting OS proliferation and lipid metabolism and seeks to shed light on the upstream regulatory mechanisms governing lipogenic enzyme activity. Based on comprehensive bioinformatic analysis and the establishment of a risk model, we identified seven lncRNAs significantly associated with clinical characteristics and lipid metabolism-related genes in patients with OS. Among these, RPARP-AS1 was selected for in-depth investigation regarding its roles in OS proliferation and lipid metabolism. Experimental techniques including RT-qPCR, Western blot, cell viability assay, assessment, and quantification of free fatty acids (FFAs) and triglycerides (TGs) were utilized to elucidate the functional significance of RPARP-AS1 in OS cells and validate its effects on lipid metabolism. Manipulation of RPARP-AS1 expression via ectopic expression or siRNA-mediated knockdown led to alterations in epithelial-mesenchymal transition (EMT) and expression of apoptosis-associated proteins, thereby influencing OS cell proliferation and apoptosis. Mechanistically, RPARP-AS1 was found to augment the expression of key lipogenic enzymes (FABP4, MAGL, and SCD1) and potentially modulate the Akt/mTOR pathway, thereby contributing to lipid metabolism (involving alterations in FFA and TG levels) in OS cells. Collectively, our findings establish RPARP-AS1 as a novel oncogene in OS cells and suggest its role in fostering tumor growth through the enhancement of lipid metabolism.

Establishment and Evaluation of Exosomes-Related Gene Risk Model in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a challenging disease to evaluate in terms of prognosis, requiring close attention to the prognosis of HCC patients. Exosomes have been shown to play an important role in HCC development and have significant potential in managing HCC patient prognosis, as they are detectable in patients' blood. By using small extracellular vesicular RNA, liquid biopsies can reflect the underlying physiological and pathological status of the originating cells, providing a valuable assessment of human health. No study has explored the diagnostic value of mRNA expression changes in exosomes for liver cancer. The present study investigated establishing a risk prognosis model based on mRNA expression levels in exosomes from blood samples of liver cancer patients and evaluated its diagnostic and prognostic value, providing new targets for liver cancer detection. We obtained mRNA data from HCC patients and normal controls from the TCGA and exoRBase 2.0 databases and established a risk prognostic assessment model using exosomes-related risk genes selected through prognostic analysis and Lasso Cox analysis. The patients were divided into high-risk and low-risk groups based on median risk score values to validate the independence and evaluability of the risk score. The clinical value of the model was further analyzed using a nomograph model, and the efficacy of immunotherapy and cell-origin types of prognostic risk genes were further assessed in the high- and low-risk groups by immune checkpoint and single-cell sequencing. A total of 44 genes were found to be significantly associated with the prognosis of HCC patients. From this group, we selected six genes (CLEC3B, CYP2C9, GNA14, NQO1, NT5DC2, and S100A9) as exosomal risk genes and used them as a basis for the risk prognosis model. The clinical information of HCC patients from the TCGA and ICGC databases demonstrated that the risk prognostic score of the model established in this study was an independent prognostic factor with good robustness. When pathological stage and risk prognostic score were incorporated into the model to predict clinical outcomes, the nomograph model had the best clinical benefit. Furthermore, immune checkpoint assays and single-cell sequencing analysis suggested that exosomal risk genes were derived from different cell types and that immunotherapy in the high-risk groups could be beneficial. Our study demonstrated that the prognostic scoring model based on exosomal mRNA was highly effective. The six genes selected using the scoring model have been previously reported to be associated with the occurrence and development of liver cancer. However, this study is the first to confirm that these related genes existed in the blood exosomes, which could be used for liquid biopsy of patients with liver cancer, thereby avoiding the need for puncture diagnosis. This approach has a high value in clinical application. Through single-cell sequencing, we found that the six genes in the risk model originate from multiple cell types. This finding suggests that the exosomal characteristic molecules secreted by different types of cells in the microenvironment of liver cancer may serve as diagnostic markers.

Gestational Diabetes Mellitus-Induced Milk Fat Globule Membrane Protein Changes of Human Mature Milk Based on TMT Proteomic Analysis.

Breastfeeding by mothers with gestational diabetes mellitus (GDM) has been shown to reduce maternal insulin demands and diminish the risks of diabetes in infants, leading to improved long-term health outcomes. Milk fat globule membrane (MFGM) proteins play a crucial role in influencing the immunity and cognitive development of infants. Understanding the alterations in MFGM proteins in breastmilk from mothers with GDM is essential for enhancing their self-efficacy and increase breastfeeding rates. The objective of this study is to investigate and compare MFGM proteins in milk from mothers with GDM and without based on tandem mass tag (TMT) labeling and liquid chromatography tandem mass spectrometry (LC-MS) techniques. A total of 5402 proteins were identified, including 4 upregulated proteins and 24 downregulated proteins. These significantly altered proteins were found to be associated with human diseases, cellular processes, and metabolism pathways. Additionally, the oxidative phosphorylation pathway emerged as the predominant pathway through Gene Set Enrichment Analysis (GSEA) involving all genes.

RNAlysis: analyze your RNA sequencing data without writing a single line of code.

BACKGROUND: Among the major challenges in next-generation sequencing experiments are exploratory data analysis, interpreting trends, identifying potential targets/candidates, and visualizing the results clearly and intuitively. These hurdles are further heightened for researchers who are not experienced in writing computer code since most available analysis tools require programming skills. Even for proficient computational biologists, an efficient and replicable system is warranted to generate standardized results. RESULTS: We have developed RNAlysis, a modular Python-based analysis software for RNA sequencing data. RNAlysis allows users to build customized analysis pipelines suiting their specific research questions, going all the way from raw FASTQ files (adapter trimming, alignment, and feature counting), through exploratory data analysis and data visualization, clustering analysis, and gene set enrichment analysis. RNAlysis provides a friendly graphical user interface, allowing researchers to analyze data without writing code. We demonstrate the use of RNAlysis by analyzing RNA sequencing data from different studies using C. elegans nematodes. We note that the software applies equally to data obtained from any organism with an existing reference genome. CONCLUSIONS: RNAlysis is suitable for investigating various biological questions, allowing researchers to more accurately and reproducibly run comprehensive bioinformatic analyses. It functions as a gateway into RNA sequencing analysis for less computer-savvy researchers, but can also help experienced bioinformaticians make their analyses more robust and efficient, as it offers diverse tools, scalability, automation, and standardization between analyses.

Mitophagy and mitochondrion-related expression profiles in response to physiological and pathological hypoxia in the corneal epithelium.

To study the mitochondrial and cellular responses to physiological and pathological hypoxia, corneal epithelial cells were preconditioned under 21% O2, 8% O2 or 1% O2. The cell survival rate, mitochondrial fluorescence and mitophagy flux were quantified using flow cytometry. After RNA sequencing, gene set enrichment analysis (GSEA) was performed. When the oxygen level decreased from 21% to 8%, mitochondrial fluorescence decreased by 45% (p < 0.001), accompanied by an 80% increase in mitophagy flux (p < 0.001). When the oxygen level dropped to 1%, the cell survival rate and mitochondrial fluorescence decreased, while mitophagy flux further increased (each p < 0.001). Comparison of 1% O2 vs. 21% O2 revealed enrichment of the HYPOXIA hallmark. Most of the significantly enriched mitochondrion-related gene sets were involved in apoptosis. The corresponding foremost leading edge genes belonged to the BCL-2 family. Corneal epithelial cell fate decisions under hypoxia may involve noncanonical pathways of mitophagy.

MMP-3 Knockout Induces Global Transcriptional Changes and Reduces Cerebral Infarction in Both Male and Female Models of Ischemic Stroke.

Ischemic stroke followed by reperfusion (IR) leads to extensive cerebrovascular injury characterized by neuroinflammation and brain cell death. Inhibition of matrix metalloproteinase-3 (MMP-3) emerges as a promising therapeutic approach to mitigate IR-induced stroke injury. We employed middle cerebral artery occlusion with subsequent reperfusion (MCAO/R) to model ischemic stroke in adult mice. Specifically, we investigated the impact of MMP-3 knockout (KO) on stroke pathophysiology using RNA sequencing (RNA-seq) of stroke brains harvested 48 h post-MCAO. MMP-3 KO significantly reduced brain infarct size following stroke. Notably, RNA-seq analysis showed that MMP-3 KO altered expression of 333 genes (252 downregulated) in male stroke brains and 3768 genes (889 downregulated) in female stroke brains. Functional pathway analysis revealed that inflammation, integrin cell surface signaling, endothelial- and epithelial-mesenchymal transition (EndMT/EMT), and apoptosis gene signatures were decreased in MMP-3 KO stroke brains. Intriguingly, MMP-3 KO downregulated gene signatures more profoundly in females than in males, as indicated by greater negative enrichment scores. Our study underscores MMP-3 inhibition as a promising therapeutic strategy, impacting multiple cellular pathways following stroke.

Analysis of small EV proteomes reveals unique functional protein networks regulated by VAP-A.

Extracellular vesicles (EVs) influence cell phenotypes and functions via protein, nucleic acid, and lipid cargoes. EVs are heterogeneous, due to diverse biogenesis mechanisms that remain poorly understood. Our previous study revealed that the endoplasmic reticulum (ER) membrane contact site (MCS) linker protein vesicle associated protein associated protein A (VAP-A) drives biogenesis of a subset of RNA-enriched EVs. Here, we examine the protein content of VAP-A-regulated EVs. Using label-free proteomics, we identified down- and upregulated proteins in small EVs (SEVs) purified from VAP-A knockdown (KD) colon cancer cells. Gene set enrichment analysis (GSEA) of the data revealed protein classes that are differentially sorted to SEVs dependent on VAP-A. Search Tool for the Retrieval of Reciprocity Genes (STRING) protein-protein interaction network analysis of the RNA-binding protein (RBP) gene set identified several RNA functional machineries that are downregulated in VAP-A KD SEVs, including ribosome, spliceosome, mRNA surveillance, and RNA transport proteins. We also observed downregulation of other functionally interacting protein networks, including cadherin-binding, unfolded protein binding, and ATP-dependent proteins.

Empowering biologists to decode omics data: the Genekitr R package and web server.

BACKGROUND: A variety of high-throughput analyses, such as transcriptome, proteome, and metabolome analysis, have been developed, producing unprecedented amounts of omics data. These studies generate large gene lists, of which the biological significance shall be deeply understood. However, manually interpreting these lists is difficult, especially for non-bioinformatics-savvy scientists. RESULTS: We developed an R package and a corresponding web server-Genekitr, to assist biologists in exploring large gene sets. Genekitr comprises four modules: gene information retrieval, ID (identifier) conversion, enrichment analysis and publication-ready plotting. Currently, the information retrieval module can retrieve information on up to 23 attributes for genes of 317 organisms. The ID conversion module assists in ID-mapping of genes, probes, proteins, and aliases. The enrichment analysis module organizes 315 gene set libraries in different biological contexts by over-representation analysis and gene set enrichment analysis. The plotting module performs customizable and high-quality illustrations that can be used directly in presentations or publications. CONCLUSIONS: This web server tool will make bioinformatics more accessible to scientists who might not have programming expertise, allowing them to perform bioinformatics tasks without coding.

Ichthyosis linked to sphingosine 1-phosphate lyase insufficiency is due to aberrant sphingolipid and calcium regulation.

Sphingosine 1-phosphate lyase (SGPL1) insufficiency (SPLIS) is a syndrome which presents with adrenal insufficiency, steroid-resistant nephrotic syndrome, hypothyroidism, neurological disease, and ichthyosis. Where a skin phenotype is reported, 94% had abnormalities such as ichthyosis, acanthosis, and hyperpigmentation. To elucidate the disease mechanism and the role SGPL1 plays in the skin barrier we established clustered regularly interspaced short palindromic repeats-Cas9 SGPL1 KO and a lentiviral-induced SGPL1 overexpression (OE) in telomerase reverse-transcriptase immortalised human keratinocytes (N/TERT-1) and thereafter organotypic skin equivalents. Loss of SGPL1 caused an accumulation of S1P, sphingosine, and ceramides, while its overexpression caused a reduction of these species. RNAseq analysis showed perturbations in sphingolipid pathway genes, particularly in SGPL1_KO, and our gene set enrichment analysis revealed polar opposite differential gene expression between SGPL1_KO and _OE in keratinocyte differentiation and Ca2+ signaling genesets. SGPL1_KO upregulated differentiation markers, while SGPL1_OE upregulated basal and proliferative markers. The advanced differentiation of SGPL1_KO was confirmed by 3D organotypic models that also presented with a thickened and retained stratum corneum and a breakdown of E-cadherin junctions. We conclude that SPLIS associated ichthyosis is a multifaceted disease caused possibly by sphingolipid imbalance and excessive S1P signaling, leading to increased differentiation and an imbalance of the lipid lamellae throughout the epidermis.