Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes.

Germline and somatic mutations can give rise to proteins with altered activity, including both gain and loss-of-function. The effects of these variants can be captured in disease-specific reactions and pathways that highlight the resulting changes to normal biology. A disease reaction is defined as an aberrant reaction in which a variant protein participates. A disease pathway is defined as a pathway that contains a disease reaction. Annotation of disease variants as participants of disease reactions and disease pathways can provide a standardized overview of molecular phenotypes of pathogenic variants that is amenable to computational mining and mathematical modeling. Reactome (https://reactome.org/), an open source, manually curated, peer-reviewed database of human biological pathways, in addition to providing annotations for >11 000 unique human proteins in the context of ∼15 000 wild-type reactions within more than 2000 wild-type pathways, also provides annotations for >4000 disease variants of close to 400 genes as participants of ∼800 disease reactions in the context of ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, described in wild-type reactions and pathways, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Reactome's data model enables mapping of disease variant datasets to specific disease reactions within disease pathways, providing a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity. Database URL: https://reactome.org/.

The landscape of cancer-rewired GPCR signaling axes.

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Data Valuation with Gradient Similarity.

High-quality data is crucial for accurate machine learning and actionable analytics, however, mislabeled or noisy data is a common problem in many domains. Distinguishing low- from high-quality data can be challenging, often requiring expert knowledge and considerable manual intervention. Data Valuation algorithms are a class of methods that seek to quantify the value of each sample in a dataset based on its contribution or importance to a given predictive task. These data values have shown an impressive ability to identify mislabeled observations, and filtering low-value data can boost machine learning performance. In this work, we present a simple alternative to existing methods, termed Data Valuation with Gradient Similarity (DVGS). This approach can be easily applied to any gradient descent learning algorithm, scales well to large datasets, and performs comparably or better than baseline valuation methods for tasks such as corrupted label discovery and noise quantification. We evaluate the DVGS method on tabular, image and RNA expression datasets to show the effectiveness of the method across domains. Our approach has the ability to rapidly and accurately identify low-quality data, which can reduce the need for expert knowledge and manual intervention in data cleaning tasks.

Targeting CCL2/CCR2 Signaling Overcomes MEK Inhibitor Resistance in Acute Myeloid Leukemia.

PURPOSE: Emerging evidence underscores the critical role of extrinsic factors within the microenvironment in protecting leukemia cells from therapeutic interventions, driving disease progression, and promoting drug resistance in acute myeloid leukemia (AML). This finding emphasizes the need for the identification of targeted therapies that inhibit intrinsic and extrinsic signaling to overcome drug resistance in AML. EXPERIMENTAL DESIGN: We performed a comprehensive analysis utilizing a cohort of ∼300 AML patient samples. This analysis encompassed the evaluation of secreted cytokines/growth factors, gene expression, and ex vivo drug sensitivity to small molecules. Our investigation pinpointed a notable association between elevated levels of CCL2 and diminished sensitivity to the MEK inhibitors (MEKi). We validated this association through loss-of-function and pharmacologic inhibition studies. Further, we deployed global phosphoproteomics and CRISPR/Cas9 screening to identify the mechanism of CCR2-mediated MEKi resistance in AML. RESULTS: Our multifaceted analysis unveiled that CCL2 activates multiple prosurvival pathways, including MAPK and cell-cycle regulation in MEKi-resistant cells. Employing combination strategies to simultaneously target these pathways heightened growth inhibition in AML cells. Both genetic and pharmacologic inhibition of CCR2 sensitized AML cells to trametinib, suppressing proliferation while enhancing apoptosis. These findings underscore a new role for CCL2 in MEKi resistance, offering combination therapies as an avenue to circumvent this resistance. CONCLUSIONS: Our study demonstrates a compelling rationale for translating CCL2/CCR2 axis inhibitors in combination with MEK pathway-targeting therapies, as a potent strategy for combating drug resistance in AML. This approach has the potential to enhance the efficacy of treatments to improve AML patient outcomes.

Graph Structured Neural Networks for Perturbation Biology.

Computational modeling of perturbation biology identifies relationships between molecular elements and cellular response, and an accurate understanding of these systems will support the full realization of precision medicine. Traditional deep learning, while often accurate in predicting response, is unlikely to capture the true sequence of involved molecular interactions. Our work is motivated by two assumptions: 1) Methods that encourage mechanistic prediction logic are likely to be more trustworthy, and 2) problem-specific algorithms are likely to outperform generic algorithms. We present an alternative to Graph Neural Networks (GNNs) termed Graph Structured Neural Networks (GSNN), which uses cell signaling knowledge, encoded as a graph data structure, to add inductive biases to deep learning. We apply our method to perturbation biology using the LINCS L1000 dataset and literature-curated molecular interactions. We demonstrate that GSNNs outperform baseline algorithms in several prediction tasks, including 1) perturbed expression, 2) cell viability of drug combinations, and 3) disease-specific drug prioritization. We also present a method called GSNNExplainer to explain GSNN predictions in a biologically interpretable form. This work has broad application in basic biological research and pre-clincal drug repurposing. Further refinement of these methods may produce trustworthy models of drug response suitable for use as clinical decision aids. Availability and implementation: Our implementation of the GSNN method is available at https://github.com/nathanieljevans/GSNN. All data used in this work is publicly available.

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.

Detection of neoplastic-immune hybrid cells with metastatic properties in uveal melanoma.

Background: Uveal melanoma is the most common non-cutaneous melanoma and is an intraocular malignancy affecting nearly 7,000 individuals per year worldwide. Of these, approximately 50% will progress to metastatic disease for which there are currently no effective therapies. Despite advances in molecular profiling and metastatic stratification of uveal melanoma tumors, little is known regarding their underlying biology of metastasis. Our group has identified a disseminated neoplastic cell population characterized by co-expression of immune and melanoma proteins, circulating hybrid cells (hybrids), in patients with uveal melanoma. Compared to circulating tumor cells, which lack expression of immune proteins, hybrids are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression. Methods: To ascertain mechanisms underlying enhanced hybrid cell dissemination we identified hybrid cells within primary uveal melanoma tumors using single cell RNA sequencing and evaluated their gene expression and predicted ligand-receptor interactions in relation to other melanoma and immune cells within the primary tumor. We then verified expression of upregulated hybrid pathways within patient-matched tumor and peripheral blood hybrids using cyclic immunofluorescence and quantified their protein expression relative to other non-hybrid tumor and disseminated tumor cells. Results: Among the top upregulated genes and pathways in hybrid cells were those involved in enhanced cell motility and cytoskeletal rearrangement, immune evasion, and altered cellular metabolism. In patient-matched tumor and peripheral blood, we verified gene expression by examining concordant protein expression for each pathway category: TMSB10 (cell motility), CD74 (immune evasion) and GPX1 (metabolism). Both TMSB10 and GPX1 were expressed on significantly higher numbers of disseminated hybrid cells compared to circulating tumor cells, and CD74 and GPX1 were expressed on more disseminated hybrids than tumor-resident hybrids. Lastly, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting metastasis including GAS6-AXL, CXCL12-CXCR4, LGALS9-P4HB and IGF1-IGFR1. Conclusion: These findings highlight the importance of TMSB10, GPX1 and CD74 for successful hybrid cell dissemination and survival in circulation. Our results contribute to the understanding of uveal melanoma tumor progression and interactions between tumor cells and immune cells in the tumor microenvironment that may promote metastasis.

Analysis of uveal melanoma scRNA sequencing data identifies neoplastic-immune hybrid cells that exhibit metastatic potential.

Uveal melanoma (UM) is the most common non-cutaneous melanoma and is an intraocular malignancy that affects nearly 7,000 individuals per year worldwide. Of these, nearly 50% will progress to metastatic disease for which there are currently no effective therapies. Despite advances in the molecular profiling and metastatic stratification of class 1 and 2 UM tumors, little is known regarding the underlying biology of UM metastasis. Our group has identified a disseminated tumor cell population characterized by co-expression of immune and melanoma proteins, (circulating hybrid cells (CHCs), in patients with UM. Compared to circulating tumor cells, CHCs are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression. To identify mechanisms underlying enhanced hybrid cell dissemination we sought to identify hybrid cells within a primary UM single cell RNA-seq dataset. Using rigorous doublet discrimination approaches, we identified UM hybrids and evaluated their gene expression, predicted ligand-receptor status, and cell-cell communication state in relation to other melanoma and immune cells within the primary tumor. We identified several genes and pathways upregulated in hybrid cells, including those involved in enhancing cell motility and cytoskeleton rearrangement, evading immune detection, and altering cellular metabolism. In addition, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting cancer metastasis including IGF1-IGFR1, GAS6-AXL, LGALS9-P4HB, APP-CD74 and CXCL12-CXCR4. These results contribute to our understanding of tumor progression and interactions between tumor cells and immune cells in the UM microenvironment that may promote metastasis.

Identification of Cellular Interactions in the Tumor Immune Microenvironment Underlying CD8 T Cell Exhaustion.

While immune checkpoint inhibitors show success in treating a subset of patients with certain late-stage cancers, these treatments fail in many other patients as a result of mechanisms that have yet to be fully characterized. The process of CD8 T cell exhaustion, by which T cells become dysfunctional in response to prolonged antigen exposure, has been implicated in immunotherapy resistance. Single-cell RNA sequencing (scRNA-seq) produces an abundance of data to analyze this process; however, due to the complexity of the process, contributions of other cell types to a process within a single cell type cannot be simply inferred. We constructed an analysis framework to first rank human skin tumor samples by degree of exhaustion in tumor-infiltrating CD8 T cells and then identify immune cell type-specific gene-regulatory network patterns significantly associated with T cell exhaustion. Using this framework, we further analyzed scRNA-seq data from human tumor and chronic viral infection samples to compare the T cell exhaustion process between these two contexts. In doing so, we identified transcription factor activity in the macrophages of both tissue types associated with this process. Our framework can be applied beyond the tumor immune microenvironment to any system involving cell-cell communication, facilitating insights into key biological processes that underpin the effective treatment of cancer and other complicated diseases.

ChatGPT usage in the Reactome curation process.

Appreciating the rapid advancement and ubiquity of generative AI, particularly ChatGPT, a chatbot using large language models like GPT, we endeavour to explore the potential application of ChatGPT in the data collection and annotation stages within the Reactome curation process. This exploration aimed to create an automated or semi-automated framework to mitigate the extensive manual effort traditionally required for gathering and annotating information pertaining to biological pathways, adopting a Reactome "reaction-centric" approach. In this pilot study, we used ChatGPT/GPT4 to address gaps in the pathway annotation and enrichment in parallel with the conventional manual curation process. This approach facilitated a comparative analysis, where we assessed the outputs generated by ChatGPT against manually extracted information. The primary objective of this comparison was to ascertain the efficiency of integrating ChatGPT or other large language models into the Reactome curation workflow and helping plan our annotation pipeline, ultimately improving our protein-to-pathway association in a reliable and automated or semi-automated way. In the process, we identified some promising capabilities and inherent challenges associated with the utilisation of ChatGPT/GPT4 in general and also specifically in the context of Reactome curation processes. We describe approaches and tools for refining the output given by ChatGPT/GPT4 that aid in generating more accurate and detailed output.

Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes.

Disease variant annotation in the context of biological reactions and pathways can provide a standardized overview of molecular phenotypes of pathogenic mutations that is amenable to computational mining and mathematical modeling. Reactome, an open source, manually curated, peer-reviewed database of human biological pathways, provides annotations for over 4000 disease variants of close to 400 genes in the context of ∼800 disease reactions constituting ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics (ACMG). Reactome's pathway-based, reaction-specific disease variant dataset and data model provide a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity.

The landscape of cancer rewired GPCR signaling axes.

We explored the dysregulation of GPCR ligand signaling systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes, which revealed that multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes. We showed that biosynthetic pathway enrichment from enzyme expression recapitulated pathway activity signatures from metabolomics datasets, providing valuable surrogate information for GPCRs responding to organic ligands. We found that several GPCRs signaling components were significantly associated with patient survival in a cancer type-specific fashion. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many pairs involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including e.g., muscarinic, adenosine, 5-hydroxytryptamine and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Illuminate the Functions of Dark Proteins Using the Reactome-IDG Web Portal.

Understudied or dark proteins have the potential to shed light on as-yet undiscovered molecular mechanisms that underlie phenotypes and suggest innovative therapeutic approaches for many diseases. The Reactome-IDG (Illuminating the Druggable Genome) project aims to place dark proteins in the context of manually curated, highly reliable pathways in Reactome, the most comprehensive, open-source biological pathway knowledgebase, facilitating the understanding functions and predicting therapeutic potentials of dark proteins. The Reactome-IDG web portal, deployed at https://idg.reactome.org, provides a simple, interactive web page for users to search pathways that may functionally interact with dark proteins, enabling the prediction of functions of dark proteins in the context of Reactome pathways. Enhanced visualization features implemented at the portal allow users to investigate the functional contexts for dark proteins based on tissue-specific gene or protein expression, drug-target interactions, or protein or gene pairwise relationships in the original Reactome's systems biology graph notation (SBGN) diagrams or the new simplified functional interaction (FI) network view of pathways. The protocols in this chapter describe step-by-step procedures to use the web portal to learn biological functions of dark proteins in the context of Reactome pathways. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Search for interacting pathways of a protein Support Protocol: Interacting pathway results for an annotated protein Alternate Protocol: Use individual pairwise relationships to predict interacting pathways of a protein Basic Protocol 2: Using the IDG pathway browser to study interacting pathways Basic Protocol 3: Overlaying tissue-specific expression data Basic Protocol 4: Overlaying protein/gene pairwise relationships in the pathway context Basic Protocol 5: Visualizing drug/target interactions.

Biology-inspired graph neural network encodes reactome and reveals biochemical reactions of disease.

Functional heterogeneity of healthy human tissues complicates interpretation of molecular studies, impeding precision therapeutic target identification and treatment. Considering this, we generated a graph neural network with Reactome-based architecture and trained it using 9,115 samples from Genotype-Tissue Expression (GTEx). Our graph neural network (GNN) achieves adjusted Rand index (ARI) = 0.7909, while a Resnet18 control model achieves ARI = 0.7781, on 370 held-out healthy human tissue samples from The Cancer Genome Atlas (TCGA), despite the Resnet18 using over 600 times the parameters. Our GNN also succeeds in separating 83 healthy skin samples from 95 lesional psoriasis samples, revealing that upregulation of 26S- and NUB1-mediated degradation of NEDD8, UBD, and their conjugates is central to the largest perturbed reaction network component in psoriasis. We show that our results are not discoverable using traditional differential expression and hypergeometric pathway enrichment analyses yet are supported by separate human multi-omics and small-molecule mouse studies, suggesting future molecular disease studies may benefit from similar GNN analytical approaches.

Illuminating Dark Proteins using Reactome Pathways.

Limited knowledge about a substantial portion of protein coding genes, known as "dark" proteins, hinders our understanding of their functions and potential therapeutic applications. To address this, we leveraged Reactome, the most comprehensive, open source, open-access pathway knowledgebase, to contextualize dark proteins within biological pathways. By integrating multiple resources and employing a random forest classifier trained on 106 protein/gene pairwise features, we predicted functional interactions between dark proteins and Reactome-annotated proteins. We then developed three scores to measure the interactions between dark proteins and Reactome pathways, utilizing enrichment analysis and fuzzy logic simulations. Correlation analysis of these scores with an independent single-cell RNA sequencing dataset provided supporting evidence for this approach. Furthermore, systematic natural language processing (NLP) analysis of over 22 million PubMed abstracts and manual checking of the literature associated with 20 randomly selected dark proteins reinforced the predicted interactions between proteins and pathways. To enhance the visualization and exploration of dark proteins within Reactome pathways, we developed the Reactome IDG portal, deployed at https://idg.reactome.org, a web application featuring tissue-specific protein and gene expression overlay, as well as drug interactions. Our integrated computational approach, together with the user-friendly web platform, offers a valuable resource for uncovering potential biological functions and therapeutic implications of dark proteins.

VIGET: A web portal for study of vaccine-induced host responses based on Reactome pathways and ImmPort data.

Host responses to vaccines are complex but important to investigate. To facilitate the study, we have developed a tool called Vaccine Induced Gene Expression Analysis Tool (VIGET), with the aim to provide an interactive online tool for users to efficiently and robustly analyze the host immune response gene expression data collected in the ImmPort/GEO databases. VIGET allows users to select vaccines, choose ImmPort studies, set up analysis models by choosing confounding variables and two groups of samples having different vaccination times, and then perform differential expression analysis to select genes for pathway enrichment analysis and functional interaction network construction using the Reactome's web services. VIGET provides features for users to compare results from two analyses, facilitating comparative response analysis across different demographic groups. VIGET uses the Vaccine Ontology (VO) to classify various types of vaccines such as live or inactivated flu vaccines, yellow fever vaccines, etc. To showcase the utilities of VIGET, we conducted a longitudinal analysis of immune responses to yellow fever vaccines and found an intriguing complex activity response pattern of pathways in the immune system annotated in Reactome, demonstrating that VIGET is a valuable web portal that supports effective vaccine response studies using Reactome pathways and ImmPort data.

Using the Reactome Database.

Pathway databases provide descriptions of the roles of proteins, nucleic acids, lipids, carbohydrates, and other molecular entities within their biological cellular contexts. Pathway-centric views of these roles may allow for the discovery of unexpected functional relationships in data such as gene expression profiles and somatic mutation catalogues from tumor cells. For this reason, there is a high demand for high-quality pathway databases and their associated tools. The Reactome project (a collaboration between the Ontario Institute for Cancer Research, New York University Langone Health, the European Bioinformatics Institute, and Oregon Health & Science University) is one such pathway database. Reactome collects detailed information on biological pathways and processes in humans from the primary literature. Reactome content is manually curated, expert-authored, and peer-reviewed and spans the gamut from simple intermediate metabolism to signaling pathways and complex cellular events. This information is supplemented with likely orthologous molecular reactions in mouse, rat, zebrafish, worm, and other model organisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Browsing a Reactome pathway Basic Protocol 2: Exploring Reactome annotations of disease and drugs Basic Protocol 3: Finding the pathways involving a gene or protein Alternate Protocol 1: Finding the pathways involving a gene or protein using UniProtKB (SwissProt), Ensembl, or Entrez gene identifier Alternate Protocol 2: Using advanced search Basic Protocol 4: Using the Reactome pathway analysis tool to identify statistically overrepresented pathways Basic Protocol 5: Using the Reactome pathway analysis tool to overlay expression data onto Reactome pathway diagrams Basic Protocol 6: Comparing inferred model organism and human pathways using the Species Comparison tool Basic Protocol 7: Comparing tissue-specific expression using the Tissue Distribution tool.

Circulating Neoplastic-Immune Hybrid Cells Predict Metastatic Progression in Uveal Melanoma.

Background: Uveal melanoma is an aggressive cancer with high metastatic risk. Recently, we identified a circulating cancer cell population that co-expresses neoplastic and leukocyte antigens, termed circulating hybrid cells (CHCs). In other cancers, CHCs are more numerous and better predict oncologic outcomes compared to circulating tumor cells (CTCs). We sought to investigate the potential of CHCs as a prognostic biomarker in uveal melanoma. Methods: We isolated peripheral blood monocular cells from uveal melanoma patients at the time of primary treatment and used antibodies against leukocyte and melanoma markers to identify and enumerate CHCs and CTCs by immunocytochemistry. Results: Using a multi-marker approach to capture the heterogeneous disseminated tumor cell population, detection of CHCs was highly sensitive in uveal melanoma patients regardless of disease stage. CHCs were detected in 100% of stage I-III uveal melanoma patients (entire cohort, n = 68), whereas CTCs were detected in 58.8% of patients. CHCs were detected at levels statically higher than CTCs across all stages (p = 0.05). Moreover, CHC levels, but not CTCs, predicted 3 year progression-free survival (p < 0.03) and overall survival (p < 0.04). Conclusion: CHCs are a novel and promising prognostic biomarker in uveal melanoma.

The reactome pathway knowledgebase 2022.

The Reactome Knowledgebase (https://reactome.org), an Elixir core resource, provides manually curated molecular details across a broad range of physiological and pathological biological processes in humans, including both hereditary and acquired disease processes. The 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. Recent curation work has expanded our annotations of normal and disease-associated signaling processes and of the drugs that target them, in particular infections caused by the SARS-CoV-1 and SARS-CoV-2 coronaviruses and the host response to infection. New tools support better simultaneous analysis of high-throughput data from multiple sources and the placement of understudied ('dark') proteins from analyzed datasets in the context of Reactome's manually curated pathways.

CIDO ontology updates and secondary analysis of host responses to COVID-19 infection based on ImmPort reports and literature.

BACKGROUND: With COVID-19 still in its pandemic stage, extensive research has generated increasing amounts of data and knowledge. As many studies are published within a short span of time, we often lose an integrative and comprehensive picture of host-coronavirus interaction (HCI) mechanisms. As of early April 2021, the ImmPort database has stored 7 studies (with 6 having details) that cover topics including molecular immune signatures, epitopes, and sex differences in terms of mortality in COVID-19 patients. The Coronavirus Infectious Disease Ontology (CIDO) represents basic HCI information. We hypothesize that the CIDO can be used as the platform to represent newly recorded information from ImmPort leading the reinforcement of CIDO. METHODS: The CIDO was used as the semantic platform for logically modeling and representing newly identified knowledge reported in the 6 ImmPort studies. A recursive eXtensible Ontology Development (XOD) strategy was established to support the CIDO representation and enhancement. Secondary data analysis was also performed to analyze different aspects of the HCI from these ImmPort studies and other related literature reports. RESULTS: The topics covered by the 6 ImmPort papers were identified to overlap with existing CIDO representation. SARS-CoV-2 viral S protein related HCI knowledge was emphasized for CIDO modeling, including its binding with ACE2, mutations causing different variants, and epitope homology by comparison with other coronavirus S proteins. Different types of cytokine signatures were also identified and added to CIDO. Our secondary analysis of two cohort COVID-19 studies with cytokine panel detection found that a total of 11 cytokines were up-regulated in female patients after infection and 8 cytokines in male patients. These sex-specific gene responses were newly modeled and represented in CIDO. A new DL query was generated to demonstrate the benefits of such integrative ontology representation. Furthermore, IL-10 signaling pathway was found to be statistically significant for both male patients and female patients. CONCLUSION: Using the recursive XOD strategy, six new ImmPort COVID-19 studies were systematically reviewed, the results were modeled and represented in CIDO, leading to the enhancement of CIDO. The enhanced ontology and further seconary analysis supported more comprehensive understanding of the molecular mechanism of host responses to COVID-19 infection.