FLAME: A Web Tool for Functional and Literature Enrichment Analysis of Multiple Gene Lists

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

Oxidative stress governs mosquito innate immune signalling to reduce chikungunya virus infection in Aedes-derived cells

Arboviruses such as chikungunya, dengue and zika viruses cause debilitating diseases in humans. The principal vector species that transmits these viruses is the Aedes mosquito. Lack of substantial knowledge of the vector species hinders the advancement of strategies for controlling the spread of arboviruses. To supplement our information on mosquitoes' responses to virus infection, we utilized Aedes aegypti-derived Aag2 cells to study changes at the transcriptional level during infection with chikungunya virus (CHIKV). We observed that genes belonging to the redox pathway were significantly differentially regulated. Upon quantifying reactive oxygen species (ROS) in the cells during viral infection, we further discovered that ROS levels are considerably higher during the early hours of infection; however, as the infection progresses, an increase in antioxidant gene expression suppresses the oxidative stress in cells. Our study also suggests that ROS is a critical regulator of viral replication in cells and inhibits intracellular and extracellular viral replication by promoting the Rel2-mediated Imd immune signalling pathway. In conclusion, our study provides evidence for a regulatory role of oxidative stress in infected Aedes-derived cells.

A genomics perspective of personalized prevention and management of obesity

This review discusses the landscape of personalized prevention and management of obesity from a nutrigenetics perspective. Focusing on macronutrient tailoring, we discuss the impact of genetic variation on responses to carbohydrate, lipid, protein, and fiber consumption. Our bioinformatic analysis of genomic variants guiding macronutrient intake revealed enrichment of pathways associated with circadian rhythm, melatonin metabolism, cholesterol and lipoprotein remodeling and PPAR signaling as potential targets of macronutrients for the management of obesity in relevant genetic backgrounds. Notably, our data-based in silico predictions suggest the potential of repurposing the SYK inhibitor fostamatinib for obesity treatment in relevant genetic profiles. In addition to dietary considerations, we address genetic variations guiding lifestyle changes in weight management, including exercise and chrononutrition. Finally, we emphasize the need for a refined understanding and expanded research into the complex genetic landscape underlying obesity and its management.

Combined genomics to discover genes associated with tolerance to soil carbonate

Carbonate-rich soils limit plant performance and crop production. Previously, local adaptation to carbonated soils was detected in wild Arabidopsis thaliana accessions, allowing the selection of two demes with contrasting phenotypes: A1 (carbonate tolerant, c+) and T6 (carbonate sensitive, c-). Here, A1(c+) and T6(c - ) seedlings were grown hydroponically under control (pH 5.9) and bicarbonate conditions (10 mM NaHCO3 , pH 8.3) to obtain ionomic profiles and conduct transcriptomic analysis. In parallel, A1(c+) and T6(c - ) parental lines and their progeny were cultivated on carbonated soil to evaluate fitness and segregation patterns. To understand the genetic architecture beyond the contrasted phenotypes, a bulk segregant analysis sequencing (BSA-Seq) was performed. Transcriptomics revealed 208 root and 2503 leaf differentially expressed genes in A1(c+) versus T6(c - ) comparison under bicarbonate stress, mainly involved in iron, nitrogen and carbon metabolism, hormones and glycosylates biosynthesis. Based on A1(c+) and T6(c - ) genome contrasts and BSA-Seq analysis, 69 genes were associated with carbonate tolerance. Comparative analysis of genomics and transcriptomics discovered a final set of 18 genes involved in bicarbonate stress responses that may have relevant roles in soil carbonate tolerance.

MFGE-8 identified in fetal mesenchymal-stromal-cell-derived exosomes ameliorates acute hepatic failure pathology

Liver transplantation is the gold-standard therapy for acute hepatic failure (AHF) with limitations related to organ shortage and life-long immunosuppressive therapy. Cell therapy emerges as a promising alternative to transplantation. We have previously shown that IL-10 and Annexin-A1 released by amniotic fluid human mesenchymal stromal cells (AF-MSCs) and their hepatocyte progenitor-like (HPL) or hepatocyte-like (HPL) cells induce liver repair and downregulate systemic inflammation in a CCl4-AHF mouse model. Herein, we demonstrate that exosomes (EXO) derived from these cells improve liver phenotype in CCl4-induced mice and promote oval cell proliferation. LC-MS/MS proteomic analysis identified MEFG-8 in EXO cargo that facilitates rescue of AHF by suppressing PI3K signaling. Administration of recombinant MFGE-8 protein also reduced liver damage in CCl4-induced mice. Clinically, MEFG-8 expression was decreased in liver biopsies from AHF patients. Collectively, our study provides proof-of-concept for an innovative, cell-free, less immunogenic, and non-toxic alternative strategy for AHF.

SCALA: A complete solution for multimodal analysis of single-cell Next Generation Sequencing data

Analysis and interpretation of high-throughput transcriptional and chromatin accessibility data at single-cell (sc) resolution are still open challenges in the biomedical field. The existence of countless bioinformatics tools, for the different analytical steps, increases the complexity of data interpretation and the difficulty to derive biological insights. In this article, we present SCALA, a bioinformatics tool for analysis and visualization of single-cell RNA sequencing (scRNA-seq) and Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) datasets, enabling either independent or integrative analysis of the two modalities. SCALA combines standard types of analysis by integrating multiple software packages varying from quality control to the identification of distinct cell populations and cell states. Additional analysis options enable functional enrichment, cellular trajectory inference, ligand-receptor analysis, and regulatory network reconstruction. SCALA is fully parameterizable, presenting data in tabular format and producing publication-ready visualizations. The different available analysis modules can aid biomedical researchers in exploring, analyzing, and visualizing their data without any prior experience in coding. We demonstrate the functionality of SCALA through two use-cases related to TNF-driven arthritic mice, handling both scRNA-seq and scATAC-seq datasets. SCALA is developed in R, Shiny and JavaScript and is mainly available as a standalone version, while an online service of more limited capacity can be found at http://scala.pavlopouloslab.info or https://scala.fleming.gr.

Flame (v2.0): advanced integration and interpretation of functional enrichment results from multiple sources

Functional enrichment is the process of identifying implicated functional terms from a given input list of genes or proteins. In this article, we present Flame (v2.0), a web tool which offers a combinatorial approach through merging and visualizing results from widely used functional enrichment applications while also allowing various flexible input options. In this version, Flame utilizes the aGOtool, g: Profiler, WebGestalt, and Enrichr pipelines and presents their outputs separately or in combination following a visual analytics approach. For intuitive representations and easier interpretation, it uses interactive plots such as parameterizable networks, heatmaps, barcharts, and scatter plots. Users can also: (i) handle multiple protein/gene lists and analyse union and intersection sets simultaneously through interactive UpSet plots, (ii) automatically extract genes and proteins from free text through text-mining and Named Entity Recognition (NER) techniques, (iii) upload single nucleotide polymorphisms (SNPs) and extract their relative genes, or (iv) analyse multiple lists of differentially expressed proteins/genes after selecting them interactively from a parameterizable volcano plot. Compared to the previous version of 197 supported organisms, Flame (v2.0) currently allows enrichment for 14 436 organisms.

AVAILABILITY AND IMPLEMENTATION

Web Application: http://flame.pavlopouloslab.info. Code: https://github.com/PavlopoulosLab/Flame. Docker: https://hub.docker.com/r/pavlopouloslab/flame.

Dry side of the core: a meta-analysis addressing the original nature of the ABA signalosome at the onset of seed imbibition

The timing of seedling emergence is a major agricultural and ecological fitness trait, and seed germination is controlled by a complex molecular network including phytohormone signalling. One such phytohormone, abscisic acid (ABA), controls a large array of stress and developmental processes, and researchers have long known it plays a crucial role in repressing germination. Although the main molecular components of the ABA signalling pathway have now been identified, the molecular mechanisms through which ABA elicits specific responses in distinct organs is still enigmatic. To address the fundamental characteristics of ABA signalling during germination, we performed a meta-analysis focusing on the Arabidopsis dry seed proteome as a reflexion basis. We combined cutting-edge proteome studies, comparative functional analyses, and protein interaction information with genetic and physiological data to redefine the singular composition and operation of the ABA core signalosome from the onset of seed imbibition. In addition, we performed a literature survey to integrate peripheral regulators present in seeds that directly regulate core component function. Although this may only be the tip of the iceberg, this extended model of ABA signalling in seeds already depicts a highly flexible system able to integrate a multitude of information to fine-tune the progression of germination.

PriPath: identifying dysregulated pathways from differential gene expression via grouping, scoring, and modeling with an embedded feature selection approach

BACKGROUND

Cell homeostasis relies on the concerted actions of genes, and dysregulated genes can lead to diseases. In living organisms, genes or their products do not act alone but within networks. Subsets of these networks can be viewed as modules that provide specific functionality to an organism. The Kyoto encyclopedia of genes and genomes (KEGG) systematically analyzes gene functions, proteins, and molecules and combines them into pathways. Measurements of gene expression (e.g., RNA-seq data) can be mapped to KEGG pathways to determine which modules are affected or dysregulated in the disease. However, genes acting in multiple pathways and other inherent issues complicate such analyses. Many current approaches may only employ gene expression data and need to pay more attention to some of the existing knowledge stored in KEGG pathways for detecting dysregulated pathways. New methods that consider more precompiled information are required for a more holistic association between gene expression and diseases.

RESULTS

PriPath is a novel approach that transfers the generic process of grouping and scoring, followed by modeling to analyze gene expression with KEGG pathways. In PriPath, KEGG pathways are utilized as the grouping function as part of a machine learning algorithm for selecting the most significant KEGG pathways. A machine learning model is trained to differentiate between diseases and controls using those groups. We have tested PriPath on 13 gene expression datasets of various cancers and other diseases. Our proposed approach successfully assigned biologically and clinically relevant KEGG terms to the samples based on the differentially expressed genes. We have comparatively evaluated the performance of PriPath against other tools, which are similar in their merit. For each dataset, we manually confirmed the top results of PriPath in the literature and found that most predictions can be supported by previous experimental research.

CONCLUSIONS

PriPath can thus aid in determining dysregulated pathways, which applies to medical diagnostics. In the future, we aim to advance this approach so that it can perform patient stratification based on gene expression and identify druggable targets. Thereby, we cover two aspects of precision medicine.

Phospholamban R14del disease: The past, the present and the future

Arrhythmogenic cardiomyopathy affects significant number of patients worldwide and is characterized by life-threatening ventricular arrhythmias and sudden cardiac death. Mutations in multiple genes with diverse functions have been reported to date including phospholamban (PLN), a key regulator of sarcoplasmic reticulum (SR) Ca²⁺ homeostasis and cardiac contractility. The PLN-R14del variant in specific is recognized as the cause in an increasing number of patients worldwide, and extensive investigations have enabled rapid advances towards the delineation of PLN-R14del disease pathogenesis and discovery of an effective treatment. We provide a critical overview of current knowledge on PLN-R14del disease pathophysiology, including clinical, animal model, cellular and biochemical studies, as well as diverse therapeutic approaches that are being pursued. The milestones achieved in <20 years, since the discovery of the PLN R14del mutation (2006), serve as a paradigm of international scientific collaboration and patient involvement towards finding a cure.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

YMLA: A comparative platform to carry out functional enrichment analysis for multiple gene lists in yeast

Comparative analysis among multiple gene lists on their functional features is now a routine task due to the advancement of high-throughput experiments. Several enrichment analysis tools were developed in the past. However, these tools mainly focus on one gene list and contain only gene ontology or interaction features. What makes it worse, comparative investigation and customized feature set reanalysis are still unavailable. Therefore, we constructed the YMLA (Yeast Multiple List Analyzer) platform in this research. YMLA includes 39 yeast features and facilitates comparative analysis among multiple gene lists via tabular views, heatmaps, and network plots. Moreover, the customized feature set reanalysis function was implemented in YMLA to help form mechanism hypotheses based on a selected enriched feature subset. We demonstrated the biological applicability of YMLA via example lists consisting of genes with top/bottom translation efficiency values. The analysis results provided by YMLA reveal novel facts consistent with previous experiments. YMLA is available at https://cosbi7.ee.ncku.edu.tw/YMLA/.

Spatial Transcriptomic Analysis Reveals Associations between Genes and Cellular Topology in Breast and Prostate Cancers

BACKGROUND

Cancer is the leading cause of death worldwide with breast and prostate cancer the most common among women and men, respectively. Gene expression and image features are independently prognostic of patient survival; but until the advent of spatial transcriptomics (ST), it was not possible to determine how gene expression of cells was tied to their spatial relationships (i.e., topology).

METHODS

We identify topology-associated genes (TAGs) that correlate with 700 image topological features (ITFs) in breast and prostate cancer ST samples. Genes and image topological features are independently clustered and correlated with each other. Themes among genes correlated with ITFs are investigated by functional enrichment analysis.

RESULTS

Overall, topology-associated genes (TAG) corresponding to extracellular matrix (ECM) and Collagen Type I Trimer gene ontology terms are common to both prostate and breast cancer. In breast cancer specifically, we identify the ZAG-PIP Complex as a TAG. In prostate cancer, we identify distinct TAGs that are enriched for GI dysmotility and the IgA immunoglobulin complex. We identified TAGs in every ST slide regardless of cancer type.

CONCLUSIONS

These TAGs are enriched for ontology terms, illustrating the biological relevance to our image topology features and their potential utility in diagnostic and prognostic models.

Drug genetic associations with COVID-19 manifestations: a data mining and network biology approach

Available drugs have been used as an urgent attempt through clinical trials to minimize severe cases of hospitalizations with Coronavirus disease (COVID-19), however, there are limited data on common pharmacogenomics affecting concomitant medications response in patients with comorbidities. To identify the genomic determinants that influence COVID-19 susceptibility, we use a computational, statistical, and network biology approach to analyze relationships of ineffective concomitant medication with an adverse effect on patients. We statistically construct a pharmacogenetic/biomarker network with significant drug-gene interactions originating from gene-disease associations. Investigation of the predicted pharmacogenes encompassing the gene-disease-gene pharmacogenomics (PGx) network suggests that these genes could play a significant role in COVID-19 clinical manifestation due to their association with autoimmune, metabolic, neurological, cardiovascular, and degenerative disorders, some of which have been reported to be crucial comorbidities in a COVID-19 patient.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

Darling: A Web Application for Detecting Disease-Related Biomedical Entity Associations with Literature Mining

Finding, exploring and filtering frequent sentence-based associations between a disease and a biomedical entity, co-mentioned in disease-related PubMed literature, is a challenge, as the volume of publications increases. Darling is a web application, which utilizes Name Entity Recognition to identify human-related biomedical terms in PubMed articles, mentioned in OMIM, DisGeNET and Human Phenotype Ontology (HPO) disease records, and generates an interactive biomedical entity association network. Nodes in this network represent genes, proteins, chemicals, functions, tissues, diseases, environments and phenotypes. Users can search by identifiers, terms/entities or free text and explore the relevant abstracts in an annotated format.

PREGO: A Literature and Data-Mining Resource to Associate Microorganisms, Biological Processes, and Environment Types

To elucidate ecosystem functioning, it is fundamental to recognize what processes occur in which environments (where) and which microorganisms carry them out (who). Here, we present PREGO, a one-stop-shop knowledge base providing such associations. PREGO combines text mining and data integration techniques to mine such what-where-who associations from data and metadata scattered in the scientific literature and in public omics repositories. Microorganisms, biological processes, and environment types are identified and mapped to ontology terms from established community resources. Analyses of comentions in text and co-occurrences in metagenomics data/metadata are performed to extract associations and a level of confidence is assigned to each of them thanks to a scoring scheme. The PREGO knowledge base contains associations for 364,508 microbial taxa, 1090 environmental types, 15,091 biological processes, and 7971 molecular functions with a total of almost 58 million associations. These associations are available through a web portal, an Application Programming Interface (API), and bulk download. By exploring environments and/or processes associated with each other or with microbes, PREGO aims to assist researchers in design and interpretation of experiments and their results. To demonstrate PREGO’s capabilities, a thorough presentation of its web interface is given along with a meta-analysis of experimental results from a lagoon-sediment study of sulfur-cycle related microbes.

Biomolecule and Bioentity Interaction Databases in Systems Biology: A Comprehensive Review

Technological advances in high-throughput techniques have resulted in tremendous growth of complex biological datasets providing evidence regarding various biomolecular interactions. To cope with this data flood, computational approaches, web services, and databases have been implemented to deal with issues such as data integration, visualization, exploration, organization, scalability, and complexity. Nevertheless, as the number of such sets increases, it is becoming more and more difficult for an end user to know what the scope and focus of each repository is and how redundant the information between them is. Several repositories have a more general scope, while others focus on specialized aspects, such as specific organisms or biological systems. Unfortunately, many of these databases are self-contained or poorly documented and maintained. For a clearer view, in this article we provide a comprehensive categorization, comparison and evaluation of such repositories for different bioentity interaction types. We discuss most of the publicly available services based on their content, sources of information, data representation methods, user-friendliness, scope and interconnectivity, and we comment on their strengths and weaknesses. We aim for this review to reach a broad readership varying from biomedical beginners to experts and serve as a reference article in the field of Network Biology.