Gene Ontology Causal Activity Modeling (GO-CAM) moves beyond GO annotations to structured descriptions of biological functions and systems

Natural but threatening? (I) A systematic aquatic ecotoxicity evaluation of biopolymers and modified natural polymers

Synthetic polymers and modified natural polymers are subject to EC 2023/2055, whereas biopolymers can be considered as important unregulated substituents. However, as regulatory requirements for biopolymers and a regulatory accepted hazard evaluation concept for polymers in general is missing, potential ecotoxicological effects are unknown. Biopolymers are often categorized as environmentally uncritical due to their origin, even though supporting data are missing. To assess potential environmental effects, we considered six biopolymers and modified natural polymers in a systematic ecotoxicity screening with aquatic organisms, hypothesising that the selected polymers are not ecotoxic. We tested alginate, chitosan, cellulose fibres Jelucel® HM 200, xanthan, carboxymethyl cellulose and modified starch Emwaxy® Jel 100 in OECD TG 201 (standard), 202 (miniaturised) and 236 (adapted to include OMICs) with algae, daphnids and zebrafish embryos, respectively. A screening of transcriptomic changes in zebrafish embryos was used to identify potential toxic modes-of-action of the polymers in fish. The polymers were applied in a concentration range of 1-100 mg/L as filtrates to consider potential intrinsic effects. Additionally, the impact of particulate polymers was evaluated by exposing algae and daphnids to polymer suspensions. Physical interactions were determined in algae growth studies leading to reduction of cell growth (lowest LOEC 10 mg/L). Daphnid immobility was observed with exposure to particulate chitosan (LOEC 100 mg/L). In zebrafish embryos, predominantly Emwaxy® Jel 100 affected the expression of genes involved in metabolic and catabolic processes (607 up- and 1002 down-regulated genes). Overall, polymer filtrates usually had no significant impact within the concentration range.

A genotype-phenotype transformer to assess and explain polygenic risk

Genome-wide association studies have linked millions of genetic variants to biomedical phenotypes, but their utility has been limited by lack of mechanistic understanding and widespread epistatic interactions. Recently, Transformer models have emerged as a powerful machine learning architecture with potential to address these and other challenges. Accordingly, here we introduce the Genotype-to-Phenotype Transformer (G2PT), a framework for modeling hierarchical information flow among variants, genes, multigenic systems, and phenotypes. As proof-of-concept, we use G2PT to model the genetics of TG/HDL (triglycerides to high-density lipoprotein cholesterol), an indicator of metabolic health. G2PT predicts this trait via attention to 1,395 variants underlying at least 20 systems, including immune response and cholesterol transport, with accuracy exceeding state-of-the-art. It implicates 40 epistatic interactions, including epistasis between APOA4 and CETP in phospholipid transfer, a target pathway for cholesterol modification. This work positions hierarchical graph transformers as a next-generation approach to polygenic risk.

Saccharomyces Genome Database: advances in genome annotation, expanded biochemical pathways, and other key enhancements

Budding yeast (Saccharomyces cerevisiae) is the most extensively characterized eukaryotic model organism and has long been used to gain insight into the fundamentals of genetics, cellular biology, and the functions of specific genes and proteins. The Saccharomyces Genome Database (SGD) is a scientific resource that provides information about the genome and biology of S. cerevisiae. For more than 30 years, SGD has maintained the genetic nomenclature, chromosome maps, and functional annotation for budding yeast along with search and analysis tools to explore these data. Here, we describe recent updates at SGD, including the 2 most recent reference genome annotation updates, expanded biochemical pathway representation, changes to SGD search and data files, and other enhancements to the SGD website and user interface. These activities are part of our continuing effort to promote insights gained from yeast to enable the discovery of functional relationships between sequence and gene products in fungi and higher eukaryotes.

Transition to Piscivory Seen Through Brain Transcriptomics in a Juvenile Percid Fish: Complex Interplay of Differential Gene Transcription, Alternative Splicing, and ncRNA Activity

Pikeperch (Sander Lucioperca) belongs to main predatory fish species in freshwater bodies throughout Europe playing the key role by reducing planktivorous fish abundance. Two size classes of the young-of-the-year (YOY) pikeperch are known in Europe and North America. Our long-term fish survey elucidates late-summer size distribution of YOY pikeperch in the Lipno Reservoir (Czechia) and recognizes two distinct subcohorts: smaller pelagic planktivores heavily outnumber larger demersal piscivores. To explore molecular mechanisms accompanying the switch from planktivory to piscivory, we compared brain transcriptomes of both subcohorts and identified 148 differentially transcribed genes. The pathway enrichment analyses identified the piscivorous phase to be associated with genes involved in collagen and extracellular matrix generation with numerous Gene Ontology (GO), while the planktivorous phase was associated with genes for non-muscle-myosins (NMM) with less GO terms. Transcripts further upregulated in planktivores from the periphery of the NMM network were Pmchl, Pomcl, and Pyyb, all involved also in appetite control and producing (an)orexigenic neuropeptides. Noncoding RNAs were upregulated in transcriptomes of planktivores including three transcripts of snoRNA U85. Thirty genes mostly functionally unrelated to those differentially transcribed were alternatively spliced between the subcohorts. Our results indicate planktivores as potentially driven by voracity to initiate the switch to piscivory, while piscivores undergo a dynamic brain development. We propose a spatiotemporal spreading of juvenile development over a longer period and larger spatial scales through developmental plasticity as an adaptation to exploiting all types of resources and decreasing the intraspecific competition.

Characterization and automated classification of sentences in the biomedical literature: a case study for biocuration of gene expression and protein kinase activity

Biological knowledgebases are essential resources for biomedical researchers, providing ready access to gene function and genomic data. Professional, manual curation of knowledgebases, however, is labor-intensive and thus high-performing machine learning methods that improve biocuration efficiency are needed. Here we report on sentence-level classification to identify biocuration-relevant sentences in the full text of published references for two gene function data types: gene expression and protein kinase activity. We performed a detailed characterization of sentences from references in the WormBase bibliography and used this characterization to define three tasks for classifying sentences as either 1) fully curatable, 2) fully and partially curatable, or 3) all language-related. We evaluated various machine learning (ML) models applied to these tasks and found that GPT and BioBERT achieve the highest average performance, resulting in F1 performance scores ranging from 0.89 to 0.99 depending upon the task. Our findings demonstrate the feasibility of extracting biocuration-relevant sentences from full text. Integrating these models into professional biocuration workflows, such as those used by the Alliance of Genome Resources and the ACKnowledge community curation platform, might well facilitate efficient and accurate annotation of the biomedical literature.

UniProt: the Universal Protein Knowledgebase in 2025

The aim of the UniProt Knowledgebase (UniProtKB; https://www.uniprot.org/) is to provide users with a comprehensive, high-quality and freely accessible set of protein sequences annotated with functional information. In this publication, we describe ongoing changes to our production pipeline to limit the sequences available in UniProtKB to high-quality, non-redundant reference proteomes. We continue to manually curate the scientific literature to add the latest functional data and use machine learning techniques. We also encourage community curation to ensure key publications are not missed. We provide an update on the automatic annotation methods used by UniProtKB to predict information for unreviewed entries describing unstudied proteins. Finally, updates to the UniProt website are described, including a new tab linking protein to genomic information. In recognition of its value to the scientific community, the UniProt database has been awarded Global Core Biodata Resource status.

EnrichDO: a global weighted model for Disease Ontology enrichment analysis

BACKGROUND

Disease Ontology (DO) has been widely studied in biomedical research and clinical practice to describe the roles of genes. DO enrichment analysis is an effective means to discover associations between genes and diseases. Compared to hundreds of Gene Ontology (GO)-based enrichment analysis methods, however, DO-based methods are relatively scarce, and most current DO-based approaches are term-for-term and thus are unable to solve over-enrichment problems caused by the "true-path" rule.

RESULTS

Here, we describe a novel double-weighted model, EnrichDO, which leverages the latest annotations of the human genome with DO terms and integrates DO graph topology on a global scale. Compared to classic enrichment methods (mainly for GO) and existing DO-based enrichment tools, EnrichDO performs better in both GO and DO enrichment analysis cases. It can accurately identify more specific terms, without ignoring the truly associated parent terms, as shown in the Alzheimer's disease (AD) case (AD ranked first). Moreover, both a simulated test and a data perturbation test validate the accuracy and robustness of EnrichDO. Finally, EnrichDO is applied to other types of datasets to expand its application, including gene expression profile datasets, a host gene set of microorganisms, and hallmark gene sets. Based on the findings reported here, EnrichDO shows significant improvement via all experimental results.

CONCLUSIONS

EnrichDO provides an effective DO enrichment analysis model for gaining insight into the significance of a particular gene set in the context of disease. To increase the usability of EnrichDO, we have developed an R-based software package, which is freely available through Bioconductor (https://bioconductor.org/packages/release/bioc/html/EnrichDO.html) or at https://github.com/liangcheng-hrbmu/EnrichDO.

DNA Methylation Signatures of Cardiovascular Health Provide Insights into Diseases

Background

The association of overall cardiovascular health (CVH) with changes in DNA methylation (DNAm) has not been well characterized.

Methods

We calculated the American Heart Association's Life's Essential 8 (LE8) score to reflect CVH in five cohorts with diverse ancestry backgrounds. Epigenome-wide association studies (EWAS) for LE8 score were conducted, followed by bioinformatic analyses. DNAm loci significantly associated with LE8 score were used to calculate a CVH DNAm score. We examined the association of the CVH DNAm score with incident CVD, CVD-specific mortality, and all-cause mortality.

Results

We identified 609 CpGs associated with LE8 score at false discovery rate (FDR) < 0.05 in the discovery analysis and at Bonferroni corrected P < 0.05 in the multi-cohort replication stage. Most had low-to-moderate heterogeneity (414 CpGs [68.0%] with I2 < 0.2) in replication analysis. Pathway enrichment analyses and phenome-wide association study (PheWAS) search associated these CpGs with inflammatory or autoimmune phenotypes. We observed enrichment for phenotypes in the EWAS catalog, with 29-fold enrichment for stroke (P = 2.4e-15) and 21-fold for ischemic heart disease (P = 7.4e-38). Two-sample Mendelian randomization (MR) analysis showed significant association between 141 CpGs and ten phenotypes (261 CpG-phenotype pairs) at FDR < 0.05. For example, hypomethylation at cg20544516 (MIR33B; SREBF1) associated with lower risk of stroke (P = 8.1e-6). In multivariable prospective analyses, the CVH DNAm score was consistently associated with clinical outcomes across participating cohorts, the reduction in risk of incident CVD, CVD mortality, and all-cause mortality per standard deviation increase in the DNAm score ranged from 19% to 32%, 28% to 40%, and 27% to 45%, respectively.

Conclusions

We identified new DNAm signatures for CVH across diverse cohorts. Our analyses indicate that immune response-related pathways may be the key mechanism underpinning the association between CVH and clinical outcomes.

Effects of different straw breeding substrates on the growth of tomato seedlings and transcriptome analysis

Traditional substrate cultivation is now a routine practice in vegetable facility breeding. However, finding renewable substrates that can replace traditional substrates is urgent in today's production. In this study, we used the 'Pindstrup' substrate as control and two types of composite substrates made from fermented corn straw (i.e. 0-3 and 3-5 mm) to identify appropriate substrate conditions for tomato seedling growth under winter greenhouse conditions. Seedling growth potential related data and substrate water content related data were tested to carry out data-oriented support. Since the single physiological data cannot well explain the mechanism of tomato seedlings under winter greenhouse condition, transcriptomic analysis of tomato root and leaf tissues were conducted to provide theoretical basis. The physiological data of tomato seedlings and substrate showed that compared with 0-3 mm and Pindstrup substrate, tomato seedlings planted in 3-5 mm had stronger growth potential and stronger water retention, and were more suitable for planting tomato seedlings. Transcriptome analysis revealed a greater number of DEGs between the Pindstrup and the 3-5 mm. The genes in this group contribute to tomato growth as well as tomato stress response mechanisms, such as ABA-related genes, hormone-related genes and some TFs. The simulation network mechanism diagram adds evidence to the above conclusions. Overall, these results demonstrate the potential benefits of using the fermented corn straw of 3-5 mm for growing tomato seedlings and present a novel method of utilizing corn straw.

Saccharomyces Genome Database: Advances in Genome Annotation, Expanded Biochemical Pathways, and Other Key Enhancements

Budding yeast (Saccharomyces cerevisiae) is the most extensively characterized eukaryotic model organism and has long been used to gain insight into the fundamentals of genetics, cellular biology, and the functions of specific genes and proteins. The Saccharomyces Genome Database (SGD) is a scientific resource that provides information about the genome and biology of S. cerevisiae. For more than 30 years, SGD has maintained the genetic nomenclature, chromosome maps, and functional annotation for budding yeast along with search and analysis tools to explore these data. Here we describe recent updates at SGD, including the two most recent reference genome annotation updates, expanded biochemical pathways representation, changes to SGD search and data files, and other enhancements to the SGD website and user interface. These activities are part of our continuing effort to promote insights gained from yeast to enable the discovery of functional relationships between sequence and gene products in fungi and higher eukaryotes.

Exploration of the link between COVID-19 and gastric cancer from the perspective of bioinformatics and systems biology

Background

Coronavirus disease 2019 (COVID-19), an infectious disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has caused a global pandemic. Gastric cancer (GC) poses a great threat to people's health, which is a high-risk factor for COVID-19. Previous studies have found some associations between GC and COVID-19, whereas the underlying molecular mechanisms are not well understood.

Methods

We employed bioinformatics and systems biology to explore these links between GC and COVID-19. Gene expression profiles of COVID-19 (GSE196822) and GC (GSE179252) were obtained from the Gene Expression Omnibus (GEO) database. After identifying the shared differentially expressed genes (DEGs) for GC and COVID-19, functional annotation, protein-protein interaction (PPI) network, hub genes, transcriptional regulatory networks and candidate drugs were analyzed.

Results

We identified 209 shared DEGs between COVID-19 and GC. Functional analyses highlighted immune-related pathways as key players in both diseases. Ten hub genes (CDK1, KIF20A, TPX2, UBE2C, HJURP, CENPA, PLK1, MKI67, IFI6, IFIT2) were identified. The transcription factor/gene and miRNA/gene interaction networks identified 38 transcription factors (TFs) and 234 miRNAs. More importantly, we identified ten potential therapeutic agents, including ciclopirox, resveratrol, etoposide, methotrexate, trifluridine, enterolactone, troglitazone, calcitriol, dasatinib and deferoxamine, some of which have been reported to improve and treat GC and COVID-19.

Conclusion

This research offer valuable insights into the molecular interplay between COVID-19 and GC, potentially guiding future therapeutic strategies.

AttentionPert: accurately modeling multiplexed genetic perturbations with multi-scale effects

MOTIVATION

Genetic perturbations (e.g. knockouts, variants) have laid the foundation for our understanding of many diseases, implicating pathogenic mechanisms and indicating therapeutic targets. However, experimental assays are fundamentally limited by the number of measurable perturbations. Computational methods can fill this gap by predicting perturbation effects under novel conditions, but accurately predicting the transcriptional responses of cells to unseen perturbations remains a significant challenge.

RESULTS

We address this by developing a novel attention-based neural network, AttentionPert, which accurately predicts gene expression under multiplexed perturbations and generalizes to unseen conditions. AttentionPert integrates global and local effects in a multi-scale model, representing both the nonuniform system-wide impact of the genetic perturbation and the localized disturbance in a network of gene-gene similarities, enhancing its ability to predict nuanced transcriptional responses to both single and multi-gene perturbations. In comprehensive experiments, AttentionPert demonstrates superior performance across multiple datasets outperforming the state-of-the-art method in predicting differential gene expressions and revealing novel gene regulations. AttentionPert marks a significant improvement over current methods, particularly in handling the diversity of gene perturbations and in predicting out-of-distribution scenarios.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/BaiDing1234/AttentionPert.

Prototyping an Ontological Framework for Cellular Senescence Mechanisms: A Homeostasis Imbalance Perspective

Although cellular senescence is a key factor in organismal aging, with both positive and negative effects on individuals, its mechanisms remain largely unknown. Thus, integrating knowledge is essential to explain how cellular senescence manifests in tissue damage and age-related diseases. Here, we propose an ontological model that organizes knowledge of cellular senescence in a computer-readable form. We manually annotated and defined cellular senescence processes, molecules, anatomical structures, phenotypes, and other entities based on the Homeostasis Imbalance Process ontology (HOIP). We described the mechanisms as causal relationships of processes and modelled a homeostatic imbalance between stress and stress response in cellular senescence for a unified framework. HOIP was assessed formally, and the relationships between cellular senescence and diseases were inferred for higher-order knowledge processing. We visualized cellular senescence processes to support knowledge utilization. Our study provides a knowledge base to help elucidate mechanisms linking cellular and organismal aging.

WormBase 2024: status and transitioning to Alliance infrastructure

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

PomBase: a Global Core Biodata Resource-growth, collaboration, and sustainability

PomBase (https://www.pombase.org), the model organism database (MOD) for fission yeast, was recently awarded Global Core Biodata Resource (GCBR) status by the Global Biodata Coalition (GBC; https://globalbiodata.org/) after a rigorous selection process. In this MOD review, we present PomBase's continuing growth and improvement over the last 2 years. We describe these improvements in the context of the qualitative GCBR indicators related to scientific quality, comprehensivity, accelerating science, user stories, and collaborations with other biodata resources. This review also showcases the depth of existing connections both within the biocuration ecosystem and between PomBase and its user community.

Updates to the Alliance of Genome Resources central infrastructure

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

Comparative transcriptomics analysis reveals defense mechanisms of Manihot esculenta Crantz against Sri Lanka Cassava MosaicVirus

BACKGROUND

Cassava mosaic disease (CMD), caused by Sri Lankan cassava mosaic virus (SLCMV) infection, has been identified as a major pernicious disease in Manihot esculenta Crantz (cassava) plantations. It is widespread in Southeast Asia, especially in Thailand, which is one of the main cassava supplier countries. With the aim of restricting the spread of SLCMV, we explored the gene expression of a tolerant cassava cultivar vs. a susceptible cassava cultivar from the perspective of transcriptional regulation and the mechanisms underlying plant immunity and adaptation.

RESULTS

Transcriptomic analysis of SLCMV-infected tolerant (Kasetsart 50 [KU 50]) and susceptible (Rayong 11 [R 11]) cultivars at three infection stages-that is, at 21 days post-inoculation (dpi) (early/asymptomatic), 32 dpi (middle/recovery), and 67 dpi (late infection/late recovery)-identified 55,699 expressed genes. Differentially expressed genes (DEGs) between SLCMV-infected KU 50 and R 11 cultivars at (i) 21 dpi to 32 dpi (the early to middle stage), and (ii) 32 dpi to 67 dpi (the middle stage to late stage) were then identified and validated by real-time quantitative PCR (RT-qPCR). DEGs among different infection stages represent genes that respond to and regulate the viral infection during specific stages. The transcriptomic comparison between the tolerant and susceptible cultivars highlighted the role of gene expression regulation in tolerant and susceptible phenotypes.

CONCLUSIONS

This study identified genes involved in epigenetic modification, transcription and transcription factor activities, plant defense and oxidative stress response, gene expression, hormone- and metabolite-related pathways, and translation and translational initiation activities, particularly in KU 50 which represented the tolerant cultivar in this study.

A comprehensive overview of post-stroke depression treatment options

Nearly one-third of all stroke patients develop depression at any time after a stroke, and its presence is associated with unfavorable outcomes. This narrative review aims to provide a synopsis of possible pharmacological and non-pharmacological treatment modalities for post-stroke depression (PSD). Several studies have demonstrated the efficacy and safety of selective serotonin reuptake inhibitors in treating the symptoms of this clinical condition. The treatment of PSD has been recently enhanced by innovative approaches, such as cognitive-behavioral therapy, virtual reality, telehealth, repetitive transcranial magnetic stimulation, and non-conventional therapies, which might improve depression treatment in stroke survivors. Future high-quality randomized controlled trials are necessary to confirm this hypothesis.

Polyoxyethylene tallow amine and glyphosate exert different developmental toxicities on human pluripotent stem cells-derived heart organoid model

The early stage of heart development is highly susceptible to various environmental factors. While the use of animal models has aided in identifying numerous environmental risk factors, the variability between species and the low throughput limit their translational potential. Recently, a type of self-assembling cardiac structures, known as human heart organoids (hHOs), exhibits a remarkable biological consistency with human heart. However, the feasibility of hHOs for assessing cardiac developmental risk factors remains unexplored. Here, we focused on the cardiac developmental effects of core components of Glyphosate-based herbicides (GBHs), the most widely used herbicides, to evaluate the reliability of hHOs for the prediction of possible cardiogenesis toxicity. GBHs have been proven toxic to cardiac development based on multiple animal models, with the mechanism remaining unknown. We found that polyoxyethylene tallow amine (POEA), the most common surfactant in GBHs formulations, played a dominant role in GBHs' heart developmental toxicity. Though there were a few differences in transcriptive features, hHOs exposed to sole POEA and combined POEA and Glyphosate would suffer from both disruption of heart contraction and disturbance of commitment in cardiomyocyte isoforms. By contrast, Glyphosate only caused mild epicardial hyperplasia. This study not only sheds light on the toxic mechanism of GBHs, but also serves as a methodological demonstration, showcasing its effectiveness in recognizing and evaluating environmental risk factors, and deciphering toxic mechanisms.

TCMNPAS: a comprehensive analysis platform integrating network formulaology and network pharmacology for exploring traditional Chinese medicine

The application of network formulaology and network pharmacology has significantly advanced the scientific understanding of traditional Chinese medicine (TCM) treatment mechanisms in disease. The field of herbal biology is experiencing a surge in data generation. However, researchers are encountering challenges due to the fragmented nature of the data and the reliance on programming tools for data analysis. We have developed TCMNPAS, a comprehensive analysis platform that integrates network formularology and network pharmacology. This platform is designed to investigate in-depth the compatibility characteristics of TCM formulas and their potential molecular mechanisms. TCMNPAS incorporates multiple resources and offers a range of functions designed for automated analysis implementation, including prescription mining, molecular docking, network pharmacology analysis, and visualization. These functions enable researchers to analyze and obtain core herbs and core formulas from herbal prescription data through prescription mining. Additionally, TCMNPAS facilitates virtual screening of active compounds in TCM and its formulas through batch molecular docking, allowing for the rapid construction and analysis of networks associated with "herb-compound-target-pathway" and disease targets. Built upon the integrated analysis concept of network formulaology and network pharmacology, TCMNPAS enables quick point-and-click completion of network-based association analysis, spanning from core formula mining from clinical data to the exploration of therapeutic targets for disease treatment. TCMNPAS serves as a powerful platform for uncovering the combinatorial rules and mechanism of TCM formulas holistically. We distribute TCMNPAS within an open-source R package at GitHub ( https://github.com/yangpluszhu/tcmnpas ), and the project is freely available at http://54.223.75.62:3838/ .

Enrichment on steps, not genes, improves inference of differentially expressed pathways

Enrichment analysis is frequently used in combination with differential expression data to investigate potential commonalities amongst lists of genes and generate hypotheses for further experiments. However, current enrichment analysis approaches on pathways ignore the functional relationships between genes in a pathway, particularly OR logic that occurs when a set of proteins can each individually perform the same step in a pathway. As a result, these approaches miss pathways with large or multiple sets because of an inflation of pathway size (when measured as the total gene count) relative to the number of steps. We address this problem by enriching on step-enabling entities in pathways. We treat sets of protein-coding genes as single entities, and we also weight sets to account for the number of genes in them using the multivariate Fisher's noncentral hypergeometric distribution. We then show three examples of pathways that are recovered with this method and find that the results have significant proportions of pathways not found in gene list enrichment analysis.

The Reactome Pathway Knowledgebase 2024

The Reactome Knowledgebase (https://reactome.org), an Elixir and GCBR core biological data resource, provides manually curated molecular details of a broad range of normal and disease-related biological processes. Processes are annotated as an ordered network of molecular transformations in a single consistent data model. Reactome thus functions both as a digital archive of manually curated human biological processes and as a tool for discovering functional relationships in data such as gene expression profiles or somatic mutation catalogs from tumor cells. Here we review progress towards annotation of the entire human proteome, targeted annotation of disease-causing genetic variants of proteins and of small-molecule drugs in a pathway context, and towards supporting explicit annotation of cell- and tissue-specific pathways. Finally, we briefly discuss issues involved in making Reactome more fully interoperable with other related resources such as the Gene Ontology and maintaining the resulting community resource network.

Representation of non-coding RNA-mediated regulation of gene expression using the Gene Ontology

Regulatory non-coding RNAs (ncRNAs) are increasingly recognized as integral to the control of biological processes. This is often through the targeted regulation of mRNA expression, but this is by no means the only mechanism through which regulatory ncRNAs act. The Gene Ontology (GO) has long been used for the systematic annotation of protein-coding and ncRNA gene function, but rapid progress in the understanding of ncRNAs meant that the ontology needed to be revised to accurately reflect current knowledge. Here, a targeted effort to revise GO terms used for the annotation of regulatory ncRNAs is described, focusing on microRNAs (miRNAs), long non-coding RNAs (lncRNAs), small interfering RNAs (siRNAs) and PIWI-interacting RNAs (piRNAs). This paper provides guidance to biocurators annotating ncRNA-mediated processes using the GO and serves as background for researchers wishing to make use of the GO in their studies of ncRNAs and the biological processes they regulate.

The alliance of genome resources: transforming comparative genomics

Comparing genomic and biological characteristics across multiple species is essential to using model systems to investigate the molecular and cellular mechanisms underlying human biology and disease and to translate mechanistic insights from studies in model organisms for clinical applications. Building a scalable knowledge commons platform that supports cross-species comparison of rich, expertly curated knowledge regarding gene function, phenotype, and disease associations available for model organisms and humans is the primary mission of the Alliance of Genome Resources (the Alliance). The Alliance is a consortium of seven model organism knowledgebases (mouse, rat, yeast, nematode, zebrafish, frog, fruit fly) and the Gene Ontology resource. The Alliance uses a common set of gene ortholog assertions as the basis for comparing biological annotations across the organisms represented in the Alliance. The major types of knowledge associated with genes that are represented in the Alliance database currently include gene function, phenotypic alleles and variants, human disease associations, pathways, gene expression, and both protein-protein and genetic interactions. The Alliance has enhanced the ability of researchers to easily compare biological annotations for common data types across model organisms and human through the implementation of shared programmatic access mechanisms, data-specific web pages with a unified "look and feel", and interactive user interfaces specifically designed to support comparative biology. The modular infrastructure developed by the Alliance allows the resource to serve as an extensible "knowledge commons" capable of expanding to accommodate additional model organisms.

Updates to the Alliance of Genome Resources Central Infrastructure Alliance of Genome Resources Consortium

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively-studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, C. elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and APIs. Here we focus on developments over the last two years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse. We describe our progress towards a central persistent database to support curation, the data modeling that underpins harmonization, and progress towards a state-of-the art literature curation system with integrated Artificial Intelligence and Machine Learning (AI/ML).

Data sharing and ontology use among agricultural genetics, genomics, and breeding databases and resources of the Agbiodata Consortium

Over the last couple of decades, there has been a rapid growth in the number and scope of agricultural genetics, genomics and breeding databases and resources. The AgBioData Consortium (https://www.agbiodata.org/) currently represents 44 databases and resources (https://www.agbiodata.org/databases) covering model or crop plant and animal GGB data, ontologies, pathways, genetic variation and breeding platforms (referred to as 'databases' throughout). One of the goals of the Consortium is to facilitate FAIR (Findable, Accessible, Interoperable, and Reusable) data management and the integration of datasets which requires data sharing, along with structured vocabularies and/or ontologies. Two AgBioData working groups, focused on Data Sharing and Ontologies, respectively, conducted a Consortium-wide survey to assess the current status and future needs of the members in those areas. A total of 33 researchers responded to the survey, representing 37 databases. Results suggest that data-sharing practices by AgBioData databases are in a fairly healthy state, but it is not clear whether this is true for all metadata and data types across all databases; and that, ontology use has not substantially changed since a similar survey was conducted in 2017. Based on our evaluation of the survey results, we recommend (i) providing training for database personnel in a specific data-sharing techniques, as well as in ontology use; (ii) further study on what metadata is shared, and how well it is shared among databases; (iii) promoting an understanding of data sharing and ontologies in the stakeholder community; (iv) improving data sharing and ontologies for specific phenotypic data types and formats; and (v) lowering specific barriers to data sharing and ontology use, by identifying sustainability solutions, and the identification, promotion, or development of data standards. Combined, these improvements are likely to help AgBioData databases increase development efforts towards improved ontology use, and data sharing via programmatic means. Database URL  https://www.agbiodata.org/databases.

Semantic representation of neural circuit knowledge in Caenorhabditis elegans

In modern biology, new knowledge is generated quickly, making it challenging for researchers to efficiently acquire and synthesise new information from the large volume of primary publications. To address this problem, computational approaches that generate machine-readable representations of scientific findings in the form of knowledge graphs have been developed. These representations can integrate different types of experimental data from multiple papers and biological knowledge bases in a unifying data model, providing a complementary method to manual review for interacting with published knowledge. The Gene Ontology Consortium (GOC) has created a semantic modelling framework that extends individual functional gene annotations to structured descriptions of causal networks representing biological processes (Gene Ontology-Causal Activity Modelling, or GO-CAM). In this study, we explored whether the GO-CAM framework could represent knowledge of the causal relationships between environmental inputs, neural circuits and behavior in the model nematode C. elegans [C. elegans Neural-Circuit Causal Activity Modelling (CeN-CAM)]. We found that, given extensions to several relevant ontologies, a wide variety of author statements from the literature about the neural circuit basis of egg-laying and carbon dioxide (CO2) avoidance behaviors could be faithfully represented with CeN-CAM. Through this process, we were able to generate generic data models for several categories of experimental results. We also discuss how semantic modelling may be used to functionally annotate the C. elegans connectome. Thus, Gene Ontology-based semantic modelling has the potential to support various machine-readable representations of neurobiological knowledge.

Biochemical pathways represented by Gene Ontology-Causal Activity Models identify distinct phenotypes resulting from mutations in pathways

Gene inactivation can affect the process(es) in which that gene acts and causally downstream ones, yielding diverse mutant phenotypes. Identifying the genetic pathways resulting in a given phenotype helps us understand how individual genes interact in a functional network. Computable representations of biological pathways include detailed process descriptions in the Reactome Knowledgebase and causal activity flows between molecular functions in Gene Ontology-Causal Activity Models (GO-CAMs). A computational process has been developed to convert Reactome pathways to GO-CAMs. Laboratory mice are widely used models of normal and pathological human processes. We have converted human Reactome GO-CAMs to orthologous mouse GO-CAMs, as a resource to transfer pathway knowledge between humans and model organisms. These mouse GO-CAMs allowed us to define sets of genes that function in a causally connected way. To demonstrate that individual variant genes from connected pathways result in similar but distinguishable phenotypes, we used the genes in our pathway models to cross-query mouse phenotype annotations in the Mouse Genome Database (MGD). Using GO-CAM representations of 2 related but distinct pathways, gluconeogenesis and glycolysis, we show that individual causal paths in gene networks give rise to discrete phenotypic outcomes resulting from perturbations of glycolytic and gluconeogenic genes. The accurate and detailed descriptions of gene interactions recovered in this analysis of well-studied processes suggest that this strategy can be applied to less well-understood processes in less well-studied model systems to predict phenotypic outcomes of novel gene variants and to identify potential gene targets in altered processes.

CTD tetramers: a new online tool that computationally links curated chemicals, genes, phenotypes, and diseases to inform molecular mechanisms for environmental health

The molecular mechanisms connecting environmental exposures to adverse endpoints are often unknown, reflecting knowledge gaps. At the Comparative Toxicogenomics Database (CTD), we developed a bioinformatics approach that integrates manually curated, literature-based interactions from CTD to generate a "CGPD-tetramer": a 4-unit block of information organized as a step-wise molecular mechanism linking an initiating Chemical, an interacting Gene, a Phenotype, and a Disease outcome. Here, we describe a novel, user-friendly tool called CTD Tetramers that generates these evidence-based CGPD-tetramers for any curated chemical, gene, phenotype, or disease of interest. Tetramers offer potential solutions for the unknown underlying mechanisms and intermediary phenotypes connecting a chemical exposure to a disease. Additionally, multiple tetramers can be assembled to construct detailed modes-of-action for chemical-induced disease pathways. As well, tetramers can help inform environmental influences on adverse outcome pathways (AOPs). We demonstrate the tool's utility with relevant use cases for a variety of environmental chemicals (eg, perfluoroalkyl substances, bisphenol A), phenotypes (eg, apoptosis, spermatogenesis, inflammatory response), and diseases (eg, asthma, obesity, male infertility). Finally, we map AOP adverse outcome terms to corresponding CTD terms, allowing users to query for tetramers that can help augment AOP pathways with additional stressors, genes, and phenotypes, as well as formulate potential AOP disease networks (eg, liver cirrhosis and prostate cancer). This novel tool, as part of the complete suite of tools offered at CTD, provides users with computational datasets and their supporting evidence to potentially fill exposure knowledge gaps and develop testable hypotheses about environmental health.

A hub gene signature as a therapeutic target and biomarker for sepsis and geriatric sepsis-induced ARDS concomitant with COVID-19 infection

Background

COVID-19 and sepsis represent formidable public health challenges, characterized by incompletely elucidated molecular mechanisms. Elucidating the interplay between COVID-19 and sepsis, particularly in geriatric patients suffering from sepsis-induced acute respiratory distress syndrome (ARDS), is of paramount importance for identifying potential therapeutic interventions to mitigate hospitalization and mortality risks.

Methods

We employed bioinformatics and systems biology approaches to identify hub genes, shared pathways, molecular biomarkers, and candidate therapeutics for managing sepsis and sepsis-induced ARDS in the context of COVID-19 infection, as well as co-existing or sequentially occurring infections. We corroborated these hub genes utilizing murine sepsis-ARDS models and blood samples derived from geriatric patients afflicted by sepsis-induced ARDS.

Results

Our investigation revealed 189 differentially expressed genes (DEGs) shared among COVID-19 and sepsis datasets. We constructed a protein-protein interaction network, unearthing pivotal hub genes and modules. Notably, nine hub genes displayed significant alterations and correlations with critical inflammatory mediators of pulmonary injury in murine septic lungs. Simultaneously, 12 displayed significant changes and correlations with a neutrophil-recruiting chemokine in geriatric patients with sepsis-induced ARDS. Of these, six hub genes (CD247, CD2, CD40LG, KLRB1, LCN2, RETN) showed significant alterations across COVID-19, sepsis, and geriatric sepsis-induced ARDS. Our single-cell RNA sequencing analysis of hub genes across diverse immune cell types furnished insights into disease pathogenesis. Functional analysis underscored the interconnection between sepsis/sepsis-ARDS and COVID-19, enabling us to pinpoint potential therapeutic targets, transcription factor-gene interactions, DEG-microRNA co-regulatory networks, and prospective drug and chemical compound interactions involving hub genes.

Conclusion

Our investigation offers potential therapeutic targets/biomarkers, sheds light on the immune response in geriatric patients with sepsis-induced ARDS, emphasizes the association between sepsis/sepsis-ARDS and COVID-19, and proposes prospective alternative pathways for targeted therapeutic interventions.

Semantic Representation of Neural Circuit Knowledge in Caenorhabditis elegans

In modern biology, new knowledge is generated quickly, making it challenging for researchers to efficiently acquire and synthesise new information from the large volume of primary publications. To address this problem, computational approaches that generate machine-readable representations of scientific findings in the form of knowledge graphs have been developed. These representations can integrate different types of experimental data from multiple papers and biological knowledge bases in a unifying data model, providing a complementary method to manual review for interacting with published knowledge. The Gene Ontology Consortium (GOC) has created a semantic modelling framework that extends individual functional gene annotations to structured descriptions of causal networks representing biological processes (Gene Ontology Causal Activity Modelling, or GO-CAM). In this study, we explored whether the GO-CAM framework could represent knowledge of the causal relationships between environmental inputs, neural circuits and behavior in the model nematode C. elegans (C. elegans Neural Circuit Causal Activity Modelling (CeN-CAM)). We found that, given extensions to several relevant ontologies, a wide variety of author statements from the literature about the neural circuit basis of egg-laying and carbon dioxide (CO2) avoidance behaviors could be faithfully represented with CeN-CAM. Through this process, we were able to generate generic data models for several categories of experimental results. We also discuss how semantic modelling may be used to functionally annotate the C. elegans connectome. Thus, Gene Ontology-based semantic modelling has the potential to support various machine-readable representations of neurobiological knowledge.

Self-prediction of relations in GO facilitates its quality auditing

As applications of the gene ontology (GO) increase rapidly in the biomedical field, quality auditing of it is becoming more and more important. Existing auditing methods are mostly based on rules, observed patterns or hypotheses. In this study, we propose a machine-learning-based framework for GO to audit itself: we first predict the IS-A relations among concepts in GO, then use differences between predicted results and existing relations to uncover potential errors. Specifically, we transfer the taxonomy of GO 2020 January release into a dataset with concept pairs as items and relations between them as labels(pairs with no direct IS-A relation are labeled as ndrs). To fully obtain the representation of each pair, we integrate the embeddings for the concept name, concept definition, as well as concept node in a substring-based topological graph. We divide the dataset into 10 parts, and rotate over all the parts by choosing one part as the testing set and the remaining as the training set each time. After 10 rotations, the prediction model predicted 4,640 existing IS-A pairs as ndrs. In the GO 2022 March release, 340 of these predictions were validated, demonstrating significance with a p-value of 1.60e-46 when compared to the results of randomly selected pairs. On the other hand, the model predicted 2,840 out of 17,079 selected ndrs in GO to be IS-A's relations. After deleting those that caused redundancies and circles, 924 predicted IS-A's relations remained. Among 200 pairs randomly selected, 30 were validated as missing IS-A's by domain experts. In conclusion, this study investigates a novel way of auditing biomedical ontologies by predicting the relations in it, which was shown to be useful for discovering potential errors.

Luteolin and triptolide: Potential therapeutic compounds for post-stroke depression via protein STAT

Post stroke depression (PSD) is a common neuropsychiatric complication following stroke closely associated with the immune system. The development of medications for PSD remains to be a considerable challenge due to the unclear mechanism of PSD. Multiple researches agree that the functions of gene ontology (GO) are efficient for the investigation of disease mechanisms, and DeepPurpose (DP) is extremely valuable for the mining of new drugs. However, GO terms and DP have not yet been applied to explore the pathogenesis and drug treatment of PSD. This study aimed to interpret the mechanism of PSD and discover important drug candidates targeting risk proteins, based on immune-related risk GO functions and informatics algorithms. According to the risk genes of PSD, we identified 335 immune-related risk GO functions and 37 compounds. Based on the construction of the GO function network, we found that STAT protein may be a pivot protein in underlying the mechanism of PSD. Additionally, we also established networks of Protein-Protein Interaction as well as Gene-GO function to facilitate the evaluation of key genes. Based on DP, a total of 37 candidate compounds targeting 7 key proteins were identified with a potential for the therapy of PSD. Furthermore, we noted that the mechanisms by which luteolin and triptolide acting on STAT-related GO function might involve three crucial pathways, including specifically hsa04010 (MAPK signaling pathway), hsa04151 (PI3K-Akt signaling pathway) and hsa04060 (Cytokine-cytokine receptor interaction). Thus, this study provided fresh and powerful information for the mechanism and therapeutic strategies of PSD.

Biochemical Pathways Represented by Gene Ontology Causal Activity Models Identify Distinct Phenotypes Resulting from Mutations in Pathways

Gene inactivation can affect the process(es) in which that gene acts and causally downstream ones, yielding diverse mutant phenotypes. Identifying the genetic pathways resulting in a given phenotype helps us understand how individual genes interact in a functional network. Computable representations of biological pathways include detailed process descriptions in the Reactome Knowledgebase, and causal activity flows between molecular functions in Gene Ontology-Causal Activity Models (GO-CAMs). A computational process has been developed to convert Reactome pathways to GO-CAMs. Laboratory mice are widely used models of normal and pathological human processes. We have converted human Reactome GO-CAMs to orthologous mouse GO-CAMs, as a resource to transfer pathway knowledge between humans and model organisms. These mouse GO-CAMs allowed us to define sets of genes that function in a connected and well-defined way. To test whether individual genes from well-defined pathways result in similar and distinguishable phenotypes, we used the genes in our pathway models to cross-query mouse phenotype annotations in the Mouse Genome Database (MGD). Using GO-CAM representations of two related but distinct pathways, gluconeogenesis and glycolysis, we can identify causal paths in gene networks that give rise to discrete phenotypic outcomes for perturbations of glycolysis and gluconeogenesis. The accurate and detailed descriptions of gene interactions recovered in this analysis of well-studied processes suggest that this strategy can be applied to less well-understood processes in less well-studied model systems to predict phenotypic outcomes of novel gene variants and to identify potential gene targets in altered processes.

Causal feature selection using a knowledge graph combining structured knowledge from the biomedical literature and ontologies: A use case studying depression as a risk factor for Alzheimer's disease

BACKGROUND

Causal feature selection is essential for estimating effects from observational data. Identifying confounders is a crucial step in this process. Traditionally, researchers employ content-matter expertise and literature review to identify confounders. Uncontrolled confounding from unidentified confounders threatens validity, conditioning on intermediate variables (mediators) weakens estimates, and conditioning on common effects (colliders) induces bias. Additionally, without special treatment, erroneous conditioning on variables combining roles introduces bias. However, the vast literature is growing exponentially, making it infeasible to assimilate this knowledge. To address these challenges, we introduce a novel knowledge graph (KG) application enabling causal feature selection by combining computable literature-derived knowledge with biomedical ontologies. We present a use case of our approach specifying a causal model for estimating the total causal effect of depression on the risk of developing Alzheimer's disease (AD) from observational data.

METHODS

We extracted computable knowledge from a literature corpus using three machine reading systems and inferred missing knowledge using logical closure operations. Using a KG framework, we mapped the output to target terminologies and combined it with ontology-grounded resources. We translated epidemiological definitions of confounder, collider, and mediator into queries for searching the KG and summarized the roles played by the identified variables. We compared the results with output from a complementary method and published observational studies and examined a selection of confounding and combined role variables in-depth.

RESULTS

Our search identified 128 confounders, including 58 phenotypes, 47 drugs, 35 genes, 23 collider, and 16 mediator phenotypes. However, only 31 of the 58 confounder phenotypes were found to behave exclusively as confounders, while the remaining 27 phenotypes played other roles. Obstructive sleep apnea emerged as a potential novel confounder for depression and AD. Anemia exemplified a variable playing combined roles.

CONCLUSION

Our findings suggest combining machine reading and KG could augment human expertise for causal feature selection. However, the complexity of causal feature selection for depression with AD highlights the need for standardized field-specific databases of causal variables. Further work is needed to optimize KG search and transform the output for human consumption.

The Gene Ontology knowledgebase in 2023

The Gene Ontology (GO) knowledgebase (http://geneontology.org) is a comprehensive resource concerning the functions of genes and gene products (proteins and noncoding RNAs). GO annotations cover genes from organisms across the tree of life as well as viruses, though most gene function knowledge currently derives from experiments carried out in a relatively small number of model organisms. Here, we provide an updated overview of the GO knowledgebase, as well as the efforts of the broad, international consortium of scientists that develops, maintains, and updates the GO knowledgebase. The GO knowledgebase consists of three components: (1) the GO-a computational knowledge structure describing the functional characteristics of genes; (2) GO annotations-evidence-supported statements asserting that a specific gene product has a particular functional characteristic; and (3) GO Causal Activity Models (GO-CAMs)-mechanistic models of molecular "pathways" (GO biological processes) created by linking multiple GO annotations using defined relations. Each of these components is continually expanded, revised, and updated in response to newly published discoveries and receives extensive QA checks, reviews, and user feedback. For each of these components, we provide a description of the current contents, recent developments to keep the knowledgebase up to date with new discoveries, and guidance on how users can best make use of the data that we provide. We conclude with future directions for the project.

Proanthocyanidins attenuates ferroptosis against influenza-induced acute lung injury in mice by reducing IFN-γ

BACKGROUND

Acute lung injury (ALI) is associated with high morbidity and mortality and is partly driven promoted by ferroptosis. Proanthocyanidins (PAs) is a natural bioactive flavonoid with anti-inflammatory and antioxidant activities. PAs can also significantly protect against acute lung inflammation and ferroptosis in alveolar epithelial cells. However, it is unclear whether PAs can alleviate ALI by reducing ferroptosis. This study aimed to evaluate the protective effects of PAs and the potential mechanisms against Influenza A virus (IAV)-induced ALI.

METHODS

Mice were inoculated nasally with IAV to induce ALI. IAV-induced pulmonary inflammation and ferroptosis was tested by measuring the levels of malondialdehyde (MDA), glutathione (GSH), glutathione peroxidase 4 (GPX4), solute carrier family 7 member 11 (SLC7A11) and acyl-CoA synthetase long-chain family member (ACSL4) in lung tissue. The potential targets that PAs protect against IAV-induced ALI were determined via a systemic pharmacological analysis. The molecular mechanism of PAs in ALI treatment was investigated by assessing the level of inflammation and ferroptosis markers using Western Blot and quantitative real-time PCR.

RESULTS

Systemic pharmacological analysis suggested that PAs protect against IAV-induced pneumonia thorough TGF-β1 and its relative signaling pathway. PAs effectively alleviated histopathological lung injury, reduced inflammatory cytokines and chemokines secretion, which were increased in IAV-infected mice. Meanwhile, PAs further prevented mouse airway inflammation in ALI, concomitant with the decreased expression TGF-β1, smad2/3, p-Smad2, p-Smad3 and ferroptosis mediator IFN-γ. Furthermore，IFN-γ promotes cell lipid peroxidation and ferroptosis，PAs significantly reduced MDA and ACSL4 levels and upregulated GSH, GPX4, and SLC7A11.

CONCLUSION

Overall, PAs can attenuate ferroptosis against IAV-induced ALI via the TGF-β1/Smad2/3 pathway and is a promising novel therapeutic candidate for ALI.

A new framework for host-pathogen interaction research

COVID-19 often manifests with different outcomes in different patients, highlighting the complexity of the host-pathogen interactions involved in manifestations of the disease at the molecular and cellular levels. In this paper, we propose a set of postulates and a framework for systematically understanding complex molecular host-pathogen interaction networks. Specifically, we first propose four host-pathogen interaction (HPI) postulates as the basis for understanding molecular and cellular host-pathogen interactions and their relations to disease outcomes. These four postulates cover the evolutionary dispositions involved in HPIs, the dynamic nature of HPI outcomes, roles that HPI components may occupy leading to such outcomes, and HPI checkpoints that are critical for specific disease outcomes. Based on these postulates, an HPI Postulate and Ontology (HPIPO) framework is proposed to apply interoperable ontologies to systematically model and represent various granular details and knowledge within the scope of the HPI postulates, in a way that will support AI-ready data standardization, sharing, integration, and analysis. As a demonstration, the HPI postulates and the HPIPO framework were applied to study COVID-19 with the Coronavirus Infectious Disease Ontology (CIDO), leading to a novel approach to rational design of drug/vaccine cocktails aimed at interrupting processes occurring at critical host-coronavirus interaction checkpoints. Furthermore, the host-coronavirus protein-protein interactions (PPIs) relevant to COVID-19 were predicted and evaluated based on prior knowledge of curated PPIs and domain-domain interactions, and how such studies can be further explored with the HPI postulates and the HPIPO framework is discussed.

Screening target genes for the treatment of PCOS via analysis of single-cell sequencing data

BACKGROUND

Polycystic ovary syndrome (PCOS) is a condition of the female reproductive system and it remains imperative to identify target genes responsible for its pathogenesis and develop therapeutic drugs capable of effectively treating it.

METHODS

We performed primary screening, staging, functional analysis as well as screening of target genes and therapeutic drugs based on single cell sequencing data of 34 oocytes from the GEO database.

RESULTS

Oxidative phosphorylation played a pivotal role in the development of oocytes, insulin resistance and ovulation disorders. At the cellular level, GV and MI phases were particularly critical for the biology of pregnancy. We screened PGR, SIRT1 and ADAMTS1 as hub differentially expressed genes (DEGs) and found relevant drugs using the Drug-Gene Interaction Database. In clinical study, oral contraceptives and insulin sensitisers were found to be effective in the treatment of PCOS.

CONCLUSION

PGR, SIRT1 and ADAMTS1 were found to be down-regulated in oocytes, ovulation and female pregnancy. These 3 genes are likely biomarkers important in the treatment of PCOS. Insulin sensitiser in combination with oral contraceptive administration were found to significantly improve PCOS.Key messagesOur study used a new bioinformatics approach to find target genes for the treatment of PCOS.Our study sought to identify target genes that affect human oocyte quality by analysing single-cell sequencing data from oocytes.We testified to our data by analysing a subset of clinical data.

Potential core genes associated with COVID-19 identified via weighted gene co-expression network analysis

AIMS

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel virus belonging to the Coronaviridae family that causes coronavirus disease (COVID-19). This disease rapidly reached pandemic status, presenting a serious threat to global health. However, the detailed molecular mechanism contributing to COVID-19 has not yet been elucidated.

METHODS

The expression profiles, including the mRNA levels, of samples from patients infected with SARS-CoV-2 along with clinical data were obtained from the GSE152075 dataset in the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules, which were then implemented to evaluate the relationships between fundamental modules and clinical traits. The differentially expressed genes (DEGs), gene ontology (GO) functional enrichment, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were evaluated using R software packages.

RESULTS

A total of 377 SARS-CoV-2-infected samples and 54 normal samples with available clinical and genetic data were obtained from the GEO database. There were 1444 DEGs identified between the sample types, which were used to screen out 11 co-expression modules in the WGCNA. Six co-expression modules were significantly associated with three clinical traits (SARS-CoV-2 positivity, age, and sex). Among the DEGs in two modules significantly correlated with SARS-CoV-2 positivity, enrichment was observed in the biological process of viral infection strategies (viral translation) in the GO analysis. The KEGG signalling pathway analysis demonstrated that the DEGs in the two modules were commonly enriched in oxidative phosphorylation, ribosome, and thermogenesis pathways. Moreover, a five-core gene set (RPL35A, RPL7A, RPS15, RPS20, and RPL17) with top connectivity with other genes was identified in the SARS-CoV-2 infection modules, suggesting that these genes may be indispensable in viral transcription after infection.

CONCLUSION

The identified core genes and signalling pathways associated with SARS-CoV-2 infection can significantly supplement the current understanding of COVID-19. The five core genes encoding ribosomal proteins may be indispensable in viral protein biosynthesis after SARS-CoV-2 infection and serve as therapeutic targets for COVID-19 treatment. These findings can be used as a basis for creating a hypothetical model for future experimental studies regarding associations of SARS-CoV-2 infection with ribosomal protein function.

Bioinformatics and systems-biology analysis to determine the effects of Coronavirus disease 2019 on patients with allergic asthma

Background

The coronavirus disease (COVID-19) pandemic has posed a significant challenge for global health systems. Increasing evidence shows that asthma phenotypes and comorbidities are major risk factors for COVID-19 symptom severity. However, the molecular mechanisms underlying the association between COVID-19 and asthma are poorly understood. Therefore, we conducted bioinformatics and systems biology analysis to identify common pathways and molecular biomarkers in patients with COVID-19 and asthma, as well as potential molecular mechanisms and candidate drugs for treating patients with both COVID-19 and asthma.

Methods

Two sets of differentially expressed genes (DEGs) from the GSE171110 and GSE143192 datasets were intersected to identify common hub genes, shared pathways, and candidate drugs. In addition, murine models were utilized to explore the expression levels and associations of the hub genes in asthma and lung inflammation/injury.

Results

We discovered 157 common DEGs between the asthma and COVID-19 datasets. A protein–protein-interaction network was built using various combinatorial statistical approaches and bioinformatics tools, which revealed several hub genes and critical modules. Six of the hub genes were markedly elevated in murine asthmatic lungs and were positively associated with IL-5, IL-13 and MUC5AC, which are the key mediators of allergic asthma. Gene Ontology and pathway analysis revealed common associations between asthma and COVID-19 progression. Finally, we identified transcription factor–gene interactions, DEG–microRNA coregulatory networks, and potential drug and chemical-compound interactions using the hub genes.

Conclusion

We identified the top 15 hub genes that can be used as novel biomarkers of COVID-19 and asthma and discovered several promising candidate drugs that might be helpful for treating patients with COVID-19 and asthma.

A roadmap for the functional annotation of protein families: a community perspective

Over the last 25 years, biology has entered the genomic era and is becoming a science of 'big data'. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3-4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

A structural homology approach to identify potential cross-reactive antibody responses following SARS-CoV-2 infection

The emergence of the novel SARS-CoV-2 virus is the most important public-health issue of our time. Understanding the diverse clinical presentations of the ensuing disease, COVID-19, remains a critical unmet need. Here we present a comprehensive listing of the diverse clinical indications associated with COVID-19. We explore the theory that anti-SARS-CoV-2 antibodies could cross-react with endogenous human proteins driving some of the pathologies associated with COVID-19. We describe a novel computational approach to estimate structural homology between SARS-CoV-2 proteins and human proteins. Antibodies are more likely to interrogate 3D-structural epitopes than continuous linear epitopes. This computational workflow identified 346 human proteins containing a domain with high structural homology to a SARS-CoV-2 Wuhan strain protein. Of these, 102 proteins exhibit functions that could contribute to COVID-19 clinical pathologies. We present a testable hypothesis to delineate unexplained clinical observations vis-à-vis COVID-19 and a tool to evaluate the safety-risk profile of potential COVID-19 therapies.

Exploring automatic inconsistency detection for literature-based gene ontology annotation

Motivation

Literature-based gene ontology annotations (GOA) are biological database records that use controlled vocabulary to uniformly represent gene function information that is described in the primary literature. Assurance of the quality of GOA is crucial for supporting biological research. However, a range of different kinds of inconsistencies in between literature as evidence and annotated GO terms can be identified; these have not been systematically studied at record level. The existing manual-curation approach to GOA consistency assurance is inefficient and is unable to keep pace with the rate of updates to gene function knowledge. Automatic tools are therefore needed to assist with GOA consistency assurance. This article presents an exploration of different GOA inconsistencies and an early feasibility study of automatic inconsistency detection.

Results

We have created a reliable synthetic dataset to simulate four realistic types of GOA inconsistency in biological databases. Three automatic approaches are proposed. They provide reasonable performance on the task of distinguishing the four types of inconsistency and are directly applicable to detect inconsistencies in real-world GOA database records. Major challenges resulting from such inconsistencies in the context of several specific application settings are reported. This is the first study to introduce automatic approaches that are designed to address the challenges in current GOA quality assurance workflows. The data underlying this article are available in Github at https://github.com/jiyuc/AutoGOAConsistency.

Harmonizing model organism data in the Alliance of Genome Resources

The Alliance of Genome Resources (the Alliance) is a combined effort of 7 knowledgebase projects: Saccharomyces Genome Database, WormBase, FlyBase, Mouse Genome Database, the Zebrafish Information Network, Rat Genome Database, and the Gene Ontology Resource. The Alliance seeks to provide several benefits: better service to the various communities served by these projects; a harmonized view of data for all biomedical researchers, bioinformaticians, clinicians, and students; and a more sustainable infrastructure. The Alliance has harmonized cross-organism data to provide useful comparative views of gene function, gene expression, and human disease relevance. The basis of the comparative views is shared calls of orthology relationships and the use of common ontologies. The key types of data are alleles and variants, gene function based on gene ontology annotations, phenotypes, association to human disease, gene expression, protein-protein and genetic interactions, and participation in pathways. The information is presented on uniform gene pages that allow facile summarization of information about each gene in each of the 7 organisms covered (budding yeast, roundworm Caenorhabditis elegans, fruit fly, house mouse, zebrafish, brown rat, and human). The harmonized knowledge is freely available on the alliancegenome.org portal, as downloadable files, and by APIs. We expect other existing and emerging knowledge bases to join in the effort to provide the union of useful data and features that each knowledge base currently provides.

Missing Links Between Gene Function and Physiology in Genomics

Knowledge of biological organisms at the molecular level that has been gathered is now organized into databases, often within ontological frameworks. To enable computational comparisons of annotations across different genomes and organisms, controlled vocabularies have been essential, as is the case in the functional annotation classifications used for bacteria, such as MultiFun and the more widely used Gene Ontology. The function of individual gene products as well as the processes in which collections of them participate constitute a wealth of classes that describe the biological role of gene products in a large number of organisms in the three kingdoms of life. In this contribution, we highlight from a qualitative perspective some limitations of these frameworks and discuss challenges that need to be addressed to bridge the gap between annotation as currently captured by ontologies and databases and our understanding of the basic principles in the organization and functioning of organisms; we illustrate these challenges with some examples in bacteria. We hope that raising awareness of these issues will encourage users of Gene Ontology and similar ontologies to be careful about data interpretation and lead to improved data representation.

Evaluation of the Prognostic Relevance of Differential Claudin Gene Expression Highlights Claudin-4 as Being Suppressed by TGFβ1 Inhibitor in Colorectal Cancer

Background: Claudins (CLDNs) are a family of closely related transmembrane proteins that have been linked to oncogenic transformation and metastasis across a range of cancers, suggesting that they may be valuable diagnostic and/or prognostic biomarkers that can be used to evaluate patient outcomes. However, CLDN expression patterns associated with colorectal cancer (CRC) remain to be defined.

Methods: The mRNA levels of 21 different CLDN family genes were assessed across 20 tumor types using the Oncomine database. Correlations between these genes and patient clinical outcomes, immune cell infiltration, clinicopathological staging, lymph node metastasis, and mutational status were analyzed using the GEPIA, UALCAN, Human Protein Atlas, Tumor Immune Estimation Resource, STRING, Genenetwork, cBioportal, and DAVID databases in an effort to clarify the potential functional roles of different CLDN protein in CRC. Molecular docking analyses were used to probe potential interactions between CLDN4 and TGFβ1. Levels of CLDN4 and CLDN11 mRNA expression in clinical CRC patient samples and in the HT29 and HCT116 cell lines were assessed via qPCR. CLDN4 expression levels in these 2 cell lines were additionally assessed following TGFβ1 inhibitor treatment.

Results: These analyses revealed that COAD and READ tissues exhibited the upregulation of CLDN1, CLDN2, CLDN3, CLDN4, CLDN7, and CLDN12 as well as the downregulation of CLDN5 and CLDN11 relative to control tissues. Higher CLDN11 and CLDN14 expression as well as lower CLDN23 mRNA levels were associated with poorer overall survival (OS) outcomes. Moreover, CLDN2 and CLDN3 or CLDN11 mRNA levels were significantly associated with lymph node metastatic progression in COAD or READ lower in COAD and READ tissues. A positive correlation between the expression of CLDN11 and predicted macrophage, dendritic cell, and CD4+ T cell infiltration was identified in CRC, with CLDN12 expression further being positively correlated with CD4+ T cell infiltration whereas a negative correlation was observed between such infiltration and the expression of CLDN3 and CLDN15. A positive correlation between CLDN1, CLDN16, and neutrophil infiltration was additionally detected, whereas neutrophil levels were negatively correlated with the expression of CLDN3 and CLDN15. Molecular docking suggested that CLDN4 was able to directly bind via hydrogen bond with TGFβ1. Relative to paracancerous tissues, clinical CRC tumor tissue samples exhibited CLDN4 and CLDN11 upregulation and downregulation, respectively. LY364947 was able to suppress the expression of CLDN4 in both the HT29 and HCT116 cell lines.

Conclusion: Together, these results suggest that the expression of different CLDN family genes is closely associated with CRC tumor clinicopathological staging and immune cell infiltration. Moreover, CLDN4 expression is closely associated with TGFβ1 in CRC, suggesting that it and other CLDN family members may represent viable targets for antitumor therapeutic intervention.

WormBase in 2022-data, processes, and tools for analyzing Caenorhabditis elegans

WormBase (www.wormbase.org) is the central repository for the genetics and genomics of the nematode Caenorhabditis elegans. We provide the research community with data and tools to facilitate the use of C. elegans and related nematodes as model organisms for studying human health, development, and many aspects of fundamental biology. Throughout our 22-year history, we have continued to evolve to reflect progress and innovation in the science and technologies involved in the study of C. elegans. We strive to incorporate new data types and richer data sets, and to provide integrated displays and services that avail the knowledge generated by the published nematode genetics literature. Here, we provide a broad overview of the current state of WormBase in terms of data type, curation workflows, analysis, and tools, including exciting new advances for analysis of single-cell data, text mining and visualization, and the new community collaboration forum. Concurrently, we continue the integration and harmonization of infrastructure, processes, and tools with the Alliance of Genome Resources, of which WormBase is a founding member.

Rhea, the reaction knowledgebase in 2022

Rhea (https://www.rhea-db.org) is an expert-curated knowledgebase of biochemical reactions based on the chemical ontology ChEBI (Chemical Entities of Biological Interest) (https://www.ebi.ac.uk/chebi). In this paper, we describe a number of key developments in Rhea since our last report in the database issue of Nucleic Acids Research in 2019. These include improved reaction coverage in Rhea, the adoption of Rhea as the reference vocabulary for enzyme annotation in the UniProt knowledgebase UniProtKB (https://www.uniprot.org), the development of a new Rhea website, and the designation of Rhea as an ELIXIR Core Data Resource. We hope that these and other developments will enhance the utility of Rhea as a reference resource to study and engineer enzymes and the metabolic systems in which they function.