The Reactome pathway knowledgebase

Potential molecular patterns for tuberculosis susceptibility in diabetic patients with poor glycaemic control: a pilot study

Type 2 diabetes (DM2) is an increasingly prevalent disease that challenges tuberculosis (TB) control strategies worldwide. It is significant that DM2 patients with poor glycemic control (PDM2) are prone to developing tuberculosis. Furthermore, elucidating the molecular mechanisms that govern this susceptibility is imperative to address this problem. Therefore, a pilot transcriptomic study was performed. Human blood samples from healthy controls (CTRL, HbA1c < 6.5%), tuberculosis (TB), comorbidity TB-DM2, DM2 (HbA1c 6.5-8.9%), and PDM2 (HbA1c > 10%) groups (n = 4 each) were analyzed by differential expression using microarrays. We use a network strategy to identify potential molecular patterns linking the differentially expressed genes (DEGs) specific for TB-DM2 and PDM2 (p-value < 0.05, fold change > 2). We define OSM, PRKCD, and SOCS3 as key regulatory genes (KRGs) that modulate the immune system and related pathways. RT-qPCR assays confirmed upregulation of OSM, PRKCD, and SOCS3 genes (p < 0.05) in TB-DM2 patients (n = 18) compared to CTRL, DM2, PDM2, or TB groups (n = 17, 19, 15, and 9, respectively). Furthermore, OSM, PRKCD, and SOCS3 were associated with PDM2 susceptibility pathways toward TB-DM2 and formed a putative protein-protein interaction confirmed in STRING. Our results reveal potential molecular patterns where OSM, PRKCD, and SOCS3 are KRGs underlying the compromised immune response and susceptibility of patients with PDM2 to develop tuberculosis. Therefore, this work paved the way for fundamental research of new molecular targets in TB-DM2. Addressing their cellular implications, and the impact on the diagnosis, treatment, and clinical management of TB-DM2 could help improve the strategy to end tuberculosis for this vulnerable population.

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/.

RCPA: An Open-Source R Package for Data Processing, Differential Analysis, Consensus Pathway Analysis, and Visualization

Identifying impacted pathways is important because it provides insights into the biology underlying conditions beyond the detection of differentially expressed genes. Because of the importance of such analysis, more than 100 pathway analysis methods have been developed thus far. Despite the availability of many methods, it is challenging for biomedical researchers to learn and properly perform pathway analysis. First, the sheer number of methods makes it challenging to learn and choose the correct method for a given experiment. Second, computational methods require users to be savvy with coding syntax, and comfortable with command-line environments, areas that are unfamiliar to most life scientists. Third, as learning tools and computational methods are typically implemented only for a few species (i.e., human and some model organisms), it is difficult to perform pathway analysis on other species that are not included in many of the current pathway analysis tools. Finally, existing pathway tools do not allow researchers to combine, compare, and contrast the results of different methods and experiments for both hypothesis testing and analysis purposes. To address these challenges, we developed an open-source R package for Consensus Pathway Analysis (RCPA) that allows researchers to conveniently: (1) download and process data from NCBI GEO; (2) perform differential analysis using established techniques developed for both microarray and sequencing data; (3) perform both gene set enrichment, as well as topology-based pathway analysis using different methods that seek to answer different research hypotheses; (4) combine methods and datasets to find consensus results; and (5) visualize analysis results and explore significantly impacted pathways across multiple analyses. This protocol provides many example code snippets with detailed explanations and supports the analysis of more than 1000 species, two pathway databases, three differential analysis techniques, eight pathway analysis tools, six meta-analysis methods, and two consensus analysis techniques. The package is freely available on the CRAN repository. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Processing Affymetrix microarrays Basic Protocol 2: Processing Agilent microarrays Support Protocol: Processing RNA sequencing (RNA-Seq) data Basic Protocol 3: Differential analysis of microarray data (Affymetrix and Agilent) Basic Protocol 4: Differential analysis of RNA-Seq data Basic Protocol 5: Gene set enrichment analysis Basic Protocol 6: Topology-based (TB) pathway analysis Basic Protocol 7: Data integration and visualization.

MRNDR: Multihead Attention-Based Recommendation Network for Drug Repurposing

As is well-known, the process of developing new drugs is extremely expensive, whereas drug repurposing represents a promising approach to augment the efficiency of new drug development. While this method can indeed spare us from expensive drug toxicity and safety experiments, it still demands a substantial amount of time to carry out precise efficacy experiments for specific diseases, thereby consuming a significant quantity of resources. Therefore, if we can prescreen potential other indications for selected drugs, it could result in substantial cost savings. In light of this, this paper introduces a drug repurposing recommendation model called MRNDR, which stands for Multi-head attention-based Recommendation Network for Drug Repurposing. This model serves as a prediction tool for drug-disease relationships, leveraging the multihead self-attention mechanism that demonstrates robust generalization capabilities. These capabilities stem not only from our extensive million-level training data set, BioRE (Biology Recommended Entity data), but also from the utilization of the WRDS (Weighted Representation Distance Score) algorithm proposed by us. The MRNDR model has achieved new state-of-the-art results on the GP-KG public data set, with an MRR (Mean Reciprocal Rank) score of 0.308 and a Hits@10 score of 0.628. This represents significant improvements of 4.7% (MRR) and 18.1% (Hits@10) over the current best-performing models. Additionally, to further validate the practical utility of the model, we examined results recommended by MRNDR that were not present in the training data set. Some of these recommendations have undergone clinical trials, as evidenced by their presence on ClinicalTrials.gov and the China Clinical Trials Center, indirectly confirming the applicability of MRNDR. The MRNDR model can predict the reusability of candidate drugs, reducing the need for manual expert assessments and enabling efficient drug repurposing.

Genetic associations of protein-coding variants in venous thromboembolism

Previous genetic studies of venous thromboembolism (VTE) have been largely limited to common variants, leaving the genetic determinants relatively incomplete. We performed an exome-wide association study of VTE among 14,723 cases and 334,315 controls. Fourteen known and four novel genes (SRSF6, PHPT1, CGN, and MAP3K2) were identified through protein-coding variants, with broad replication in the FinnGen cohort. Most genes we discovered exhibited the potential to predict future VTE events in longitudinal analysis. Notably, we provide evidence for the additive contribution of rare coding variants to known genome-wide polygenic risk in shaping VTE risk. The identified genes were enriched in pathways affecting coagulation and platelet activation, along with liver-specific expression. The pleiotropic effects of these genes indicated the potential involvement of coagulation factors, blood cell traits, liver function, and immunometabolic processes in VTE pathogenesis. In conclusion, our study unveils the valuable contribution of protein-coding variants in VTE etiology and sheds new light on its risk stratification.

Alterations to the Hidradenitis Suppurativa Serum Proteome with Spleen Tyrosine Kinase Antagonism: Proteomic Results from a Phase 2 Clinical Trial

Hidradenitis suppurativa is a disease in great need of novel therapies. Given the heterogeneous nature of the disease and the variable response to therapies, biomarkers are essential to predict response to therapies and increase our understanding of disease pathogenesis. Our recent phase 2 clinical trial of spleen tyrosine kinase antagonism using fostamatinib in hidradenitis suppurativa demonstrated a 75% clinical response, with the greatest benefit in individuals with elevated serum inflammation and IgG. In this study, we present results of an in-depth serum proteomic analysis in this patient cohort identifying downregulation of IL-12B as well as B-cell-associated proteins CCL19 and CCL20 and IFN-γ-mediated proteins CXCL10 and CX3CL1. Clinical responders demonstrated greater reduction in serum IL-17A, IL-6, IL-8, and CX3CL1 compared with clinical nonresponders. Baseline levels of CCL28 were associated with clinical response to fostamatinib therapy at week 12. Overall, this suggests that fostamatinib, by targeting B-cell receptor and Fc receptor activity in B cells, monocytes, and macrophages, has a significant molecular impact on the inflammatory serum proteome of hidradenitis suppurativa. In addition, potential therapeutic biomarkers may aid in patient selection for targeted therapy.

Ranking of cell clusters in a single-cell RNA-sequencing analysis framework using prior knowledge

Prioritization or ranking of different cell types in a single-cell RNA sequencing (scRNA-seq) framework can be performed in a variety of ways, some of these include: i) obtaining an indication of the proportion of cell types between the different conditions under study, ii) counting the number of differentially expressed genes (DEGs) between cell types and conditions in the experiment or, iii) prioritizing cell types based on prior knowledge about the conditions under study (i.e., a specific disease). These methods have drawbacks and limitations thus novel methods for improving cell ranking are required. Here we present a novel methodology that exploits prior knowledge in combination with expert-user information to accentuate cell types from a scRNA-seq analysis that yield the most biologically meaningful results with respect to a disease under study. Our methodology allows for ranking and prioritization of cell types based on how well their expression profiles relate to the molecular mechanisms and drugs associated with a disease. Molecular mechanisms, as well as drugs, are incorporated as prior knowledge in a standardized, structured manner. Cell types are then ranked/prioritized based on how well results from data-driven analysis of scRNA-seq data match the predefined prior knowledge. In additional cell-cell communication perturbations between disease and control networks are used to further prioritize/rank cell types. Our methodology has substantial advantages to more traditional cell ranking techniques and provides an informative complementary methodology that utilizes prior knowledge in a rapid and automated manner, that has previously not been attempted by other studies. The current methodology is also implemented as an R package entitled Single Cell Ranking Analysis Toolkit (scRANK) and is available for download and installation via GitHub (https://github.com/aoulas/scRANK).

Reconstruction Set Test (RESET): A computationally efficient method for single sample gene set testing based on randomized reduced rank reconstruction error

We have developed a new, and analytically novel, single sample gene set testing method called Reconstruction Set Test (RESET). RESET quantifies gene set importance based on the ability of set genes to reconstruct values for all measured genes. RESET is realized using a computationally efficient randomized reduced rank reconstruction algorithm (available via the RESET R package on CRAN) that can effectively detect patterns of differential abundance and differential correlation for self-contained and competitive scenarios. As demonstrated using real and simulated scRNA-seq data, RESET provides superior performance at a lower computational cost relative to other single sample approaches.

Inhibition of Acyl-CoA Synthetase Long Chain Isozymes Decreases Multiple Myeloma Cell Proliferation and Causes Mitochondrial Dysfunction

Multiple myeloma (MM) is an incurable cancer of plasma cells with a 5-year survival rate of 59%. Dysregulation of fatty acid (FA) metabolism is associated with MM development and progression; however, the underlying mechanisms remain unclear. Acyl-CoA synthetase long-chain family members (ACSLs) convert free long-chain fatty acids into fatty acyl-CoA esters and play key roles in catabolic and anabolic fatty acid metabolism. The Cancer Dependency Map data suggested that ACSL3 and ACSL4 were among the top 25% Hallmark Fatty Acid Metabolism genes that support MM fitness. Here, we show that inhibition of ACSLs in human myeloma cell lines using the pharmacological inhibitor Triascin C (TriC) causes apoptosis and decreases proliferation in a dose- and time-dependent manner. RNA-seq of MM.1S cells treated with TriC for 24 h showed a significant enrichment in apoptosis, ferroptosis, and ER stress. Proteomics of MM.1S cells treated with TriC for 48 h revealed that mitochondrial dysfunction and oxidative phosphorylation were significantly enriched pathways of interest, consistent with our observations of decreased mitochondrial membrane potential and increased mitochondrial superoxide levels. Interestingly, MM.1S cells treated with TriC for 24 h also showed decreased mitochondrial ATP production rates and overall lower cellular respiration.

Mitf regulates gene expression networks implicated in B cell homeostasis, germinal center responses, and tolerance

Introduction

The microphthalmia transcription factor Mitf has been shown to regulate B cell activation and tolerance. However, the underlying B cell-specific mechanisms responsible, and those that distinguish Mitf from closely related Mitf/TFE (MiT) transcription factors Tfe3, Tfeb, and Tfec, remain obscure.

Methods

Two complementary mouse models of Mitf and MiT deficiency were used: the Mitfmi-vga9/mi-vga9 systemic loss-of-function mutation, and B-cell specific MiT family inactivation via transgenic expression of a trans-dominant negative (TDN) protein (TDN-B). These models were employed to identify MiT family candidate target genes and pathways.

Results

Both models displayed spontaneous splenomegaly coincident with elevated plasma cell numbers, autoantibody titers, and proteinuria. These abnormalities appeared dependent on T helper cells, but independent of other non-B cell intrinsic effects of systemic Mitf inactivation. MiT inactivation in B cells augmented aspects of lupus-like autoimmune disease on the C57BL/6-Faslpr/lpr background. In both models, RNAseq of ex vivo resting B cells showed transcriptional upregulation of genes that control cell cycle, germinal center responses, and plasma cell differentiation. Among the genes strongly upregulated in both models were Socs6, Isp53 (Baiap1), S1pR2, and IgG2b/c. Mitf null B cells, but not TDN-B cells, showed evidence of type I interferon dysregulation.

Discussion

These studies clarify Mitf's role as 1) a key regulator of a B cell intrinsic germinal center program that influences self-tolerance through novel target genes, and 2) a regulator of systemic inflammatory processes that can impact the B cell microenvironment. This distinction of Mitf's function from that of related MiT transcription factors advances our understanding of B cell regulation and autoimmunity.

Gene set correlation enrichment analysis for interpreting and annotating gene expression profiles

Pathway analysis, including nontopology-based (non-TB) and topology-based (TB) methods, is widely used to interpret the biological phenomena underlying differences in expression data between two phenotypes. By considering dependencies and interactions between genes, TB methods usually perform better than non-TB methods in identifying pathways that include closely relevant or directly causative genes for a given phenotype. However, most TB methods may be limited by incomplete pathway data used as the reference network or by difficulties in selecting appropriate reference networks for different research topics. Here, we propose a gene set correlation enrichment analysis method, Gscore, based on an expression dataset-derived coexpression network to examine whether a differentially expressed gene (DEG) list (or each of its DEGs) is associated with a known gene set. Gscore is better able to identify target pathways in 89 human disease expression datasets than eight other state-of-the-art methods and offers insight into how disease-wide and pathway-wide associations reflect clinical outcomes. When applied to RNA-seq data from COVID-19-related cells and patient samples, Gscore provided a means for studying how DEGs are implicated in COVID-19-related pathways. In summary, Gscore offers a powerful analytical approach for annotating individual DEGs, DEG lists, and genome-wide expression profiles based on existing biological knowledge.

Protein Profiling of Placental Extracellular Vesicles in Gestational Diabetes Mellitus

Throughout pregnancy, some degree of insulin resistance is necessary to divert glucose towards the developing foetus. In gestational diabetes mellitus (GDM), insulin resistance is exacerbated in combination with insulin deficiency, causing new-onset maternal hyperglycaemia. The rapid reversal of insulin resistance following delivery strongly implicates the placenta in GDM pathogenesis. In this case-control study, we investigated the proteomic cargo of human syncytiotrophoblast-derived extracellular vesicles (STBEVs), which facilitate maternal-fetal signalling during pregnancy, in a UK-based cohort comprising patients with a gestational age of 38-40 weeks. Medium/large (m/l) and small (s) STBEVs were isolated from GDM (n = 4) and normal (n = 5) placentae using ex vivo dual-lobe perfusion and subjected to mass spectrometry. Bioinformatics were used to identify differentially carried proteins and mechanistic pathways. In m/lSTBEVs, 56 proteins were differently expressed while in sSTBEVs, no proteins reached statistical difference. Differences were also observed in the proteomic cargo between m/lSTBEVs and sSTBEVs, indicating that the two subtypes of STBEVs may have divergent modes of action and downstream effects. In silico functional enrichment analysis of differentially expressed proteins in m/lSTBEVs from GDM and normal pregnancy found positive regulation of cytoskeleton organisation as the most significantly enriched biological process. This work presents the first comparison of two populations of STBEVs' protein cargos (m/l and sSTBEVs) from GDM and normal pregnancy isolated using placenta perfusion. Further investigation of differentially expressed proteins may contribute to an understanding of GDM pathogenesis and the development of novel diagnostic and therapeutic tools.

Evolutionary proteogenomic landscape from pre-invasive to invasive lung adenocarcinoma

Lung adenocarcinoma follows a stepwise progression from pre-invasive to invasive. However, there remains a knowledge gap regarding molecular events from pre-invasive to invasive. Here, we conduct a comprehensive proteogenomic analysis comprising whole-exon sequencing, RNA sequencing, and proteomic and phosphoproteomic profiling on 98 pre-invasive and 99 invasive lung adenocarcinomas. The deletion of chr4q12 contributes to the progression from pre-invasive to invasive adenocarcinoma by downregulating SPATA18, thus suppressing mitophagy and promoting cell invasion. Proteomics reveals diverse enriched pathways in normal lung tissues and pre-invasive and invasive adenocarcinoma. Proteomic analyses identify three proteomic subtypes, which represent different stages of tumor progression. We also illustrate the molecular characterization of four immune clusters, including endothelial cells, B cells, DCs, and immune depression subtype. In conclusion, this comprehensive proteogenomic study characterizes the molecular architecture and hallmarks from pre-invasive to invasive lung adenocarcinoma, guiding the way to a deeper understanding of the tumorigenesis and progression of this disease.

Identification of potential biological processes and key genes in diabetes-related stroke through weighted gene co-expression network analysis

BACKGROUND

Type 2 diabetes mellitus (T2DM) is an established risk factor for acute ischemic stroke (AIS). Although there are reports on the correlation of diabetes and stroke, data on its pathogenesis is limited. This study aimed to explore the underlying biological mechanisms and promising intervention targets of diabetes-related stroke.

METHODS

Diabetes-related datasets (GSE38642 and GSE44035) and stroke-related datasets (GSE16561 and GSE22255) were obtained from the Gene Expression omnibus (GEO) database. The key modules for stroke and diabetes were identified by weight gene co-expression network analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) analyses were employed in the key module. Genes in stroke- and diabetes-related key modules were intersected to obtain common genes for T2DM-related stroke. In order to discover the key genes in T2DM-related stroke, the Cytoscape and protein-protein interaction (PPI) network were constructed. The key genes were functionally annotated in the Reactome database.

RESULTS

By intersecting the diabetes- and stroke-related crucial modules, 24 common genes for T2DM-related stroke were identified. Metascape showed that neutrophil extracellular trap formation was primarily enriched. The hub gene was granulin precursor (GRN), which had the highest connectivity among the common genes. In addition, functional enrichment analysis indicated that GRN was involved in neutrophil degranulation, thus regulating neutrophil extracellular trap formation.

CONCLUSIONS

This study firstly revealed that neutrophil extracellular trap formation may represent the common biological processes of diabetes and stroke, and GRN may be potential intervention targets for T2DM-related stroke.

Tissue-adjusted pathway analysis of cancer (TPAC): A novel approach for quantifying tumor-specific gene set dysregulation relative to normal tissue

We describe a novel single sample gene set testing method for cancer transcriptomics data named tissue-adjusted pathway analysis of cancer (TPAC). The TPAC method leverages information about the normal tissue-specificity of human genes to compute a robust multivariate distance score that quantifies gene set dysregulation in each profiled tumor. Because the null distribution of the TPAC scores has an accurate gamma approximation, both population and sample-level inference is supported. As we demonstrate through an analysis of gene expression data for 21 solid human cancers from The Cancer Genome Atlas (TCGA) and associated normal tissue expression data from the Human Protein Atlas (HPA), TPAC gene set scores are more strongly associated with patient prognosis than the scores generated by existing single sample gene set testing methods.

HTNpedia: A Knowledge Base for Hypertension Research

BACKGROUND

Hypertension is notably a serious public health concern due to its high prevalence and strong association with cardiovascular disease and renal failure. It is reported to be the fourth leading disease that leads to death worldwide.

OBJECTIVES

Currently, there is no active operational knowledge base or database for hypertension or cardiovascular illness.

METHODS

The primary data source was retrieved from the research outputs obtained from our laboratory team working on hypertension research. We have presented a preliminary dataset and external links to the public repository for detailed analysis to readers.

RESULTS

As a result, HTNpedia was created to provide information regarding hypertension-related proteins and genes.

CONCLUSION

The complete webpage is accessible via www.mkarthikeyan.bioinfoau.org/HTNpedia.

Large-scale functional inference for skin-expressing lncRNAs using expression and sequence information

Long noncoding RNAs (lncRNAs) regulate the expression of protein-coding genes and have been shown to play important roles in inflammatory skin diseases. However, we still have limited understanding of the functional impact of lncRNAs in skin, partly due to their tissue specificity and lower expression levels compared with protein-coding genes. We compiled a comprehensive list of 18,517 lncRNAs from different sources and studied their expression profiles in 834 RNA-Seq samples from multiple inflammatory skin conditions and cytokine-stimulated keratinocytes. Applying a balanced random forest to predict involvement in biological functions, we achieved a median AUROC of 0.79 in 10-fold cross-validation, identifying significant DNA binding domains (DBDs) for 39 lncRNAs. G18244, a skin-expressing lncRNA predicted for IL-4/IL-13 signaling in keratinocytes, was highly correlated in expression with F13A1, a protein-coding gene involved in macrophage regulation, and we further identified a significant DBD in F13A1 for G18244. Reflecting clinical implications, AC090198.1 (predicted for IL-17 pathway) and AC005332.6 (predicted for IFN-γ pathway) had significant negative correlation with the SCORAD metric for atopic dermatitis. We also utilized single-cell RNA and spatial sequencing data to validate cell type specificity. Our research demonstrates lncRNAs have important immunological roles and can help prioritize their impact on inflammatory skin diseases.

Germline gene fusions across species reveal the chromosomal instability regions and cancer susceptibility

The canine transmissible venereal tumor (CTVT) is a clonal cell-mediated cancer with a long evolutionary history and extensive karyotype rearrangements in its genome. However, little is known about its genetic similarity to human tumors. Here, using multi-omics data we identified 11 germline gene fusions (GGFs) in CTVT, which showed higher genetic susceptibility than others. Additionally, we illustrate a mechanism of a complex gene fusion of three gene segments (HSD17B4-DMXL1-TNFAIP8) that we refer to "greedy fusion". Our findings also provided evidence that expressions of GGFs are downregulated during the tumor regressive phase, which is associated with DNA methylation level. This study presents a comprehensive landscape of gene fusions (GFs) in CTVT, which offers a valuable genetic resource for exploring potential genetic mechanisms underlying the development of cancers in both dogs and humans.

Advances in computational methods for identifying cancer driver genes

Cancer driver genes (CDGs) are crucial in cancer prevention, diagnosis and treatment. This study employed computational methods for identifying CDGs, categorizing them into four groups. The major frameworks for each of these four categories were summarized. Additionally, we systematically gathered data from public databases and biological networks, and we elaborated on computational methods for identifying CDGs using the aforementioned databases. Further, we summarized the algorithms, mainly involving statistics and machine learning, used for identifying CDGs. Notably, the performances of nine typical identification methods for eight types of cancer were compared to analyze the applicability areas of these methods. Finally, we discussed the challenges and prospects associated with methods for identifying CDGs. The present study revealed that the network-based algorithms and machine learning-based methods demonstrated superior performance.

Enzyme Databases in the Era of Omics and Artificial Intelligence

Enzyme research is important for the development of various scientific fields such as medicine and biotechnology. Enzyme databases facilitate this research by providing a wide range of information relevant to research planning and data analysis. Over the years, various databases that cover different aspects of enzyme biology (e.g., kinetic parameters, enzyme occurrence, and reaction mechanisms) have been developed. Most of the databases are curated manually, which improves reliability of the information; however, such curation cannot keep pace with the exponential growth in published data. Lack of data standardization is another obstacle for data extraction and analysis. Improving machine readability of databases is especially important in the light of recent advances in deep learning algorithms that require big training datasets. This review provides information regarding the current state of enzyme databases, especially in relation to the ever-increasing amount of generated research data and recent advancements in artificial intelligence algorithms. Furthermore, it describes several enzyme databases, providing the reader with necessary information for their use.

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.

Unveiling the protein signature of the human osteosarcoma 3AB-OS cancer stem cell line

In cancer research today, one of the major challenges is the eradication of cancer stem cells (CSCs) within the tumor mass. These cells play a crucial role in initiating, growing, and maintaining the tumor. Evidence has demonstrated the presence and significance of CSCs in the development and progression of osteosarcoma (OS). However, our understanding of the specific markers for OS stem cells remains limited. In this study, we aim to identify distinct biomarkers for this cell population by conducting a proteomic analysis comparing OS stem cells to their non-stem counterparts. Our investigation focuses on a particular cell line called 3AB-OS, which exhibits stem-like characteristics, and its differentiated parental cell line, MG63. Through this research, we discovered 63 proteins exclusively expressed in 3AB-OS cells. Applying an in silico bioinformatics approach, we determined that the majority of these proteins are associated with RNA metabolism. Additionally, we identified a potential correlation between the insulin-like growth factor-binding proteins (IGF2BPs) signaling pathway and the tumorigenic and stemness features observed in 3AB-OS cells.

Transcriptional Regulation of Airway Epithelial Cell Differentiation: Insights into the Notch Pathway and Beyond

The airway epithelium is a critical component of the respiratory system, serving as a barrier against inhaled pathogens and toxins. It is composed of various cell types, each with specific functions essential to proper airway function. Chronic respiratory diseases can disrupt the cellular composition of the airway epithelium, leading to a decrease in multiciliated cells (MCCs) and an increase in secretory cells (SCs). Basal cells (BCs) have been identified as the primary stem cells in the airway epithelium, capable of self-renewal and differentiation into MCCs and SCs. This review emphasizes the role of transcription factors in the differentiation process from BCs to MCCs and SCs. Recent advancements in single-cell RNA sequencing (scRNAseq) techniques have provided insights into the cellular composition of the airway epithelium, revealing specialized and rare cell types, including neuroendocrine cells, tuft cells, and ionocytes. Understanding the cellular composition and differentiation processes within the airway epithelium is crucial for developing targeted therapies for respiratory diseases. Additionally, the maintenance of BC populations and the involvement of Notch signaling in BC self-renewal and differentiation are discussed. Further research in these areas could provide valuable insights into the mechanisms underlying airway epithelial homeostasis and disease pathogenesis.

Proteogenomics of different urothelial bladder cancer stages reveals distinct molecular features for papillary cancer and carcinoma in situ

The progression of urothelial bladder cancer (UC) is a complicated multi-step process. We perform a comprehensive multi-omics analysis of 448 samples from 190 UC patients, covering the whole spectrum of disease stages and grades. Proteogenomic integration analysis indicates the mutations of HRAS regulated mTOR signaling to form urothelial papilloma rather than papillary urothelial cancer (PUC). DNA damage is a key signaling pathway in the progression of carcinoma in situ (CIS) and related to APOBEC signature. Glucolipid metabolism increase and lower immune cell infiltration are associated with PUC compared to CIS. Proteomic analysis distinguishes the origins of invasive tumors (PUC-derived and CIS-derived), related to distinct clinical prognosis and molecular features. Additionally, loss of RBPMS, associated with CIS-derived tumors, is validated to increase the activity of AP-1 and promote metastasis. This study reveals the characteristics of two distinct branches (PUC and CIS) of UC progression and may eventually benefit clinical practice.

Human blood vessel organoids reveal a critical role for CTGF in maintaining microvascular integrity

The microvasculature plays a key role in tissue perfusion and exchange of gases and metabolites. In this study we use human blood vessel organoids (BVOs) as a model of the microvasculature. BVOs fully recapitulate key features of the human microvasculature, including the reliance of mature endothelial cells on glycolytic metabolism, as concluded from metabolic flux assays and mass spectrometry-based metabolomics using stable tracing of 13C-glucose. Pharmacological targeting of PFKFB3, an activator of glycolysis, using two chemical inhibitors results in rapid BVO restructuring, vessel regression with reduced pericyte coverage. PFKFB3 mutant BVOs also display similar structural remodelling. Proteomic analysis of the BVO secretome reveal remodelling of the extracellular matrix and differential expression of paracrine mediators such as CTGF. Treatment with recombinant CTGF recovers microvessel structure. In this work we demonstrate that BVOs rapidly undergo restructuring in response to metabolic changes and identify CTGF as a critical paracrine regulator of microvascular integrity.

Uniformly shaped harmonization combines human transcriptomic data from different platforms while retaining their biological properties and differential gene expression patterns

Introduction: Co-normalization of RNA profiles obtained using different experimental platforms and protocols opens avenue for comprehensive comparison of relevant features like differentially expressed genes associated with disease. Currently, most of bioinformatic tools enable normalization in a flexible format that depends on the individual datasets under analysis. Thus, the output data of such normalizations will be poorly compatible with each other. Recently we proposed a new approach to gene expression data normalization termed Shambhala which returns harmonized data in a uniform shape, where every expression profile is transformed into a pre-defined universal format. We previously showed that following shambhalization of human RNA profiles, overall tissue-specific clustering features are strongly retained while platform-specific clustering is dramatically reduced. Methods: Here, we tested Shambhala performance in retention of fold-change gene expression features and other functional characteristics of gene clusters such as pathway activation levels and predicted cancer drug activity scores. Results: Using 6,793 cancer and 11,135 normal tissue gene expression profiles from the literature and experimental datasets, we applied twelve performance criteria for different versions of Shambhala and other methods of transcriptomic harmonization with flexible output data format. Such criteria dealt with the biological type classifiers, hierarchical clustering, correlation/regression properties, stability of drug efficiency scores, and data quality for using machine learning classifiers. Discussion: Shambhala-2 harmonizer demonstrated the best results with the close to 1 correlation and linear regression coefficients for the comparison of training vs validation datasets and more than two times lesser instability for calculation of drug efficiency scores compared to other methods.

A co-formulation of interferons alpha2b and gamma distinctively targets cell cycle in the glioblastoma-derived cell line U-87MG

BACKGROUND

HeberFERON is a co-formulation of α2b and γ interferons, based on their synergism, which has shown its clinical superiority over individual interferons in basal cell carcinomas. In glioblastoma (GBM), HeberFERON has displayed promising preclinical and clinical results. This led us to design a microarray experiment aimed at identifying the molecular mechanisms involved in the distinctive effect of HeberFERON compared to the individual interferons in U-87MG model.

METHODS

Transcriptional expression profiling including a control (untreated) and three groups receiving α2b-interferon, γ-interferon and HeberFERON was performed using an Illumina HT-12 microarray platform. Unsupervised methods for gene and sample grouping, identification of differentially expressed genes, functional enrichment and network analysis computational biology methods were applied to identify distinctive transcription patterns of HeberFERON. Validation of most representative genes was performed by qPCR. For the cell cycle analysis of cells treated with HeberFERON for 24 h, 48 and 72 h we used flow cytometry.

RESULTS

The three treatments show different behavior based on the gene expression profiles. The enrichment analysis identified several mitotic cell cycle related events, in particular from prometaphase to anaphase, which are exclusively targeted by HeberFERON. The FOXM1 transcription factor network that is involved in several cell cycle phases and is highly expressed in GBMs, is significantly down regulated. Flow cytometry experiments corroborated the action of HeberFERON on the cell cycle in a dose and time dependent manner with a clear cellular arrest as of 24 h post-treatment. Despite the fact that p53 was not down-regulated, several genes involved in its regulatory activity were functionally enriched. Network analysis also revealed a strong relationship of p53 with genes targeted by HeberFERON. We propose a mechanistic model to explain this distinctive action, based on the simultaneous activation of PKR and ATF3, p53 phosphorylation changes, as well as its reduced MDM2 mediated ubiquitination and export from the nucleus to the cytoplasm. PLK1, AURKB, BIRC5 and CCNB1 genes, all regulated by FOXM1, also play central roles in this model. These and other interactions could explain a G2/M arrest and the effect of HeberFERON on the proliferation of U-87MG.

CONCLUSIONS

We proposed molecular mechanisms underlying the distinctive behavior of HeberFERON compared to the treatments with the individual interferons in U-87MG model, where cell cycle related events were highly relevant.

Algorithmically Reconstructed Molecular Pathways as the New Generation of Prognostic Molecular Biomarkers in Human Solid Cancers

Individual gene expression and molecular pathway activation profiles were shown to be effective biomarkers in many cancers. Here, we used the human interactome model to algorithmically build 7470 molecular pathways centered around individual gene products. We assessed their associations with tumor type and survival in comparison with the previous generation of molecular pathway biomarkers (3022 "classical" pathways) and with the RNA transcripts or proteomic profiles of individual genes, for 8141 and 1117 samples, respectively. For all analytes in RNA and proteomic data, respectively, we found a total of 7441 and 7343 potential biomarker associations for gene-centric pathways, 3020 and 2950 for classical pathways, and 24,349 and 6742 for individual genes. Overall, the percentage of RNA biomarkers was statistically significantly higher for both types of pathways than for individual genes (p < 0.05). In turn, both types of pathways showed comparable performance. The percentage of cancer-type-specific biomarkers was comparable between proteomic and transcriptomic levels, but the proportion of survival biomarkers was dramatically lower for proteomic data. Thus, we conclude that pathway activation level is the advanced type of biomarker for RNA and proteomic data, and momentary algorithmic computer building of pathways is a new credible alternative to time-consuming hypothesis-driven manual pathway curation and reconstruction.

Dysbindin Domain-Containing 1 in Prostate Cancer: New Insights into Bioinformatic Validation of Molecular and Immunological Features

Prostate cancer (PCa) is one of the most prevalent cancers in men, yet its pathogenic pathways remain poorly understood. Transcriptomics and high-throughput sequencing can help uncover cancer diagnostic targets and understand biological circuits. Using prostate adenocarcinoma (PRAD) datasets of various web-based applications (GEPIA, UALCAN, cBioPortal, SR Plot, hTFtarget, Genome Browser, and MetaCore), we found that upregulated dysbindin domain-containing 1 (DBNDD1) expression in primary prostate tumors was strongly correlated with pathways involving the cell cycle, mitotic in KEGG, WIKI, and REACTOME database, and transcription factor-binding sites with the DBNDD1 gene in prostate samples. DBNDD1 gene expression was influenced by sample type, cancer stage, and promoter methylation levels of different cancers, such as PRAD, liver hepatocellular carcinoma (LIHC), and lung adenocarcinoma (LUAD). Regulation of glycogen synthase kinase (GSK)-3β in bipolar disorder and ATP/ITP/GTP/XTP/TTP/CTP/UTP metabolic pathways was closely correlated with the DBNDD1 gene and its co-expressed genes in PCa. DBNDD1 gene expression was positively associated with immune infiltration of B cells, Myeloid-derived suppressor cell (MDSC), M2 macrophages, andneutrophil, whereas negatively correlated with CD8+ T cells, T follicular helper cells, M1 macrophages, and NK cells in PCa. These findings suggest that DBNDD1 may serve as a viable prognostic marker not only for early-stage PCa but also for immunotherapies.

A synthetic lethal approach to drug targeting of G-quadruplexes based on CX-5461

DNA G-quadruplex (G4) structures are enriched at human genome loci critical for cancer development, such as in oncogene promoters, telomeres, and rDNA. Medicinal chemistry approaches to developing drugs that target G4 structures date back to over 20 years ago. Small-molecule drugs were designed to target and stabilize G4 structures, thereby blocking replication and transcription, resulting in cancer cell death. CX-3543 (Quarfloxin) was the first G4-targeting drug to enter clinical trials in 2005; however, because of the lack of efficacy, it was withdrawn from Phase 2 clinical trials. Efficacy problems also occurred in the clinical trial of patients with advanced hematologic malignancies using CX-5461 (Pidnarulex), another G4-stabilizing drug. Only after the discovery of synthetic lethal (SL) interactions between Pidnarulex and the BRCA1/2-mediated homologous recombination (HR) pathway in 2017, promising clinical efficacy was achieved. In this case, Pidnarulex was used in a clinical trial to treat solid tumors deficient in BRCA2 and PALB2. The history of the development of Pidnarulex highlights the importance of SL in identifying cancer patients responsive to G4-targeting drugs. In order to identify additional cancer patients responsive to Pidnarulex, several genetic interaction screens have been performed with Pidnarulex and other G4-targeting drugs using human cancer cell lines or C. elegans. Screening results confirmed the synthetic lethal interaction between G4 stabilizers and HR genes and also uncovered other novel genetic interactions, including genes in other DNA damage repair pathways and genes in transcription, epigenetic, and RNA processing deficiencies. In addition to patient identification, synthetic lethality is also important for the design of drug combination therapy for G4-targeting drugs in order to achieve better clinical outcomes.

The salivary proteome in relation to oral mucositis in autologous hematopoietic stem cell transplantation recipients: a labelled and label-free proteomics approach

BACKGROUND

Oral mucositis is a frequently seen complication in the first weeks after hematopoietic stem cell transplantation recipients which can severely affects patients quality of life. In this study, a labelled and label-free proteomics approach were used to identify differences between the salivary proteomes of autologous hematopoietic stem cell transplantation (ASCT) recipients developing ulcerative oral mucositis (ULC-OM; WHO score ≥ 2) or not (NON-OM).

METHODS

In the TMT-labelled analysis we pooled saliva samples from 5 ULC-OM patients at each of 5 timepoints: baseline, 1, 2, 3 weeks and 3 months after ASCT and compared these with pooled samples from 5 NON-OM patients. For the label-free analysis we analyzed saliva samples from 9 ULC-OM and 10 NON-OM patients at 6 different timepoints (including 12 months after ASCT) with Data-Independent Acquisition (DIA). As spectral library, all samples were grouped (ULC-OM vs NON-OM) and analyzed with Data Dependent Analysis (DDA). PCA plots and a volcano plot were generated in RStudio and differently regulated proteins were analyzed using GO analysis with g:Profiler.

RESULTS

A different clustering of ULC-OM pools was found at baseline, weeks 2 and 3 after ASCT with TMT-labelled analysis. Using label-free analysis, week 1-3 samples clustered distinctly from the other timepoints. Unique and up-regulated proteins in the NON-OM group (DDA analysis) were involved in immune system-related processes, while those proteins in the ULC-OM group were intracellular proteins indicating cell lysis.

CONCLUSIONS

The salivary proteome in ASCT recipients has a tissue protective or tissue-damage signature, that corresponded with the absence or presence of ulcerative oral mucositis, respectively.

TRIAL REGISTRATION

The study is registered in the national trial register (NTR5760; automatically added to the International Clinical Trial Registry Platform).

Genetic differences according to onset age and lung function in asthma: A cluster analysis

BACKGROUND

The extent of differences between genetic risks associated with various asthma subtypes is still unknown. To better understand the heterogeneity of asthma, we employed an unsupervised method to identify genetic variants specifically associated with asthma subtypes. Our goal was to gain insight into the genetic basis of asthma.

METHODS

In this study, we utilized the UK Biobank dataset to select asthma patients (All asthma, n = 50,517) and controls (n = 283,410). We excluded 14,431 individuals who had no information on predicted values of forced expiratory volume in one second percent (FEV1%) and onset age, resulting in a final total of 36,086 asthma cases. We conducted k-means clustering based on asthma onset age and predicted FEV1% using these samples (n = 36,086). Cluster-specific genome-wide association studies were then performed, and heritability was estimated via linkage disequilibrium score regression. To further investigate the pathophysiology, we conducted eQTL analysis with GTEx and gene-set enrichment analysis with FUMA.

RESULTS

Clustering resulted in four distinct clusters: early onset asthma^normalLF (early onset with normal lung function, n = 8172), early onset asthma^reducedLF (early onset with reduced lung function, n = 8925), late-onset asthma^normalLF (late-onset with normal lung function, n = 12,481), and late-onset asthma^reducedLF (late-onset with reduced lung function, n = 6508). Our GWASs in four clusters and in All asthma sample identified 5 novel loci, 14 novel signals, and 51 cluster-specific signals. Among clusters, early onset asthma^normalLF and late-onset asthma^reducedLF were the least correlated (r_g  = 0.37). Early onset asthma^reducedLF showed the highest heritability explained by common variants (h²  = 0.212) and was associated with the largest number of variants (71 single nucleotide polymorphisms). Further, the pathway analysis conducted through eQTL and gene-set enrichment analysis showed that the worsening of symptoms in early onset asthma correlated with lymphocyte activation, pathogen recognition, cytokine receptor activation, and lymphocyte differentiation.

CONCLUSIONS

Our findings suggest that early onset asthma^reducedLF was the most genetically predisposed cluster, and that asthma clusters with reduced lung function were genetically distinct from clusters with normal lung function. Our study revealed the genetic variation between clusters that were segmented based on onset age and lung function, providing an important clue for the genetic mechanism of asthma heterogeneity.

Gene and pathway based burden analyses in familial lymphoid cancer cases: Rare variants in immune pathway genes

Genome-wide association studies have revealed common genetic variants with small effect sizes associated with diverse lymphoid cancers. Family studies have uncovered rare variants with high effect sizes. However, these variants explain only a portion of the heritability of these cancers. Some of the missing heritability may be attributable to rare variants with small effect sizes. We aim to identify rare germline variants associated with familial lymphoid cancers using exome sequencing. One case per family was selected from 39 lymphoid cancer families based on early onset of disease or rarity of subtype. Control data was from Non-Finnish Europeans in gnomAD exomes (N = 56,885) or ExAC (N = 33,370). Gene and pathway-based burden tests for rare variants were performed using TRAPD. Five putatively pathogenic germline variants were found in four genes: INTU, PEX7, EHHADH, and ASXL1. Pathway-based association tests identified the innate and adaptive immune systems, peroxisomal pathway and olfactory receptor pathway as associated with lymphoid cancers in familial cases. Our results suggest that rare inherited defects in the genes involved in immune system and peroxisomal pathway may predispose individuals to lymphoid cancers.

Systems-level analyses of protein-protein interaction network dysfunctions via epichaperomics identify cancer-specific mechanisms of stress adaptation

Systems-level assessments of protein-protein interaction (PPI) network dysfunctions are currently out-of-reach because approaches enabling proteome-wide identification, analysis, and modulation of context-specific PPI changes in native (unengineered) cells and tissues are lacking. Herein, we take advantage of chemical binders of maladaptive scaffolding structures termed epichaperomes and develop an epichaperome-based 'omics platform, epichaperomics, to identify PPI alterations in disease. We provide multiple lines of evidence, at both biochemical and functional levels, demonstrating the importance of these probes to identify and study PPI network dysfunctions and provide mechanistically and therapeutically relevant proteome-wide insights. As proof-of-principle, we derive systems-level insight into PPI dysfunctions of cancer cells which enabled the discovery of a context-dependent mechanism by which cancer cells enhance the fitness of mitotic protein networks. Importantly, our systems levels analyses support the use of epichaperome chemical binders as therapeutic strategies aimed at normalizing PPI networks.

A multi-scale map of protein assemblies in the DNA damage response

The DNA damage response (DDR) ensures error-free DNA replication and transcription and is disrupted in numerous diseases. An ongoing challenge is to determine the proteins orchestrating DDR and their organization into complexes, including constitutive interactions and those responding to genomic insult. Here, we use multi-conditional network analysis to systematically map DDR assemblies at multiple scales. Affinity purifications of 21 DDR proteins, with/without genotoxin exposure, are combined with multi-omics data to reveal a hierarchical organization of 605 proteins into 109 assemblies. The map captures canonical repair mechanisms and proposes new DDR-associated proteins extending to stress, transport, and chromatin functions. We find that protein assemblies closely align with genetic dependencies in processing specific genotoxins and that proteins in multiple assemblies typically act in multiple genotoxin responses. Follow-up by DDR functional readouts newly implicates 12 assembly members in double-strand-break repair. The DNA damage response assemblies map is available for interactive visualization and query (ccmi.org/ddram/).

Signaling mechanisms in renal compensatory hypertrophy revealed by multi-omics

Loss of a kidney results in compensatory growth of the remaining kidney, a phenomenon of considerable clinical importance. However, the mechanisms involved are largely unknown. Here, we use a multi-omic approach in a unilateral nephrectomy model in male mice to identify signaling processes associated with renal compensatory hypertrophy, demonstrating that the lipid-activated transcription factor peroxisome proliferator-activated receptor alpha (PPARα) is an important determinant of proximal tubule cell size and is a likely mediator of compensatory proximal tubule hypertrophy.

Bioinformatic analysis of gene expression data reveals Src family protein tyrosine kinases as key players in androgenetic alopecia

Introduction

Androgenetic alopecia (AGA) is a common progressive scalp hair loss disorder that leads to baldness. This study aimed to identify core genes and pathways involved in premature AGA through an in-silico approach.

Methods

Gene expression data (GSE90594) from vertex scalps of men with premature AGA and men without pattern hair loss was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between the bald and haired samples were identified using the limma package in R. Gene ontology and Reactome pathway enrichment analyses were conducted separately for the up-regulated and down-regulated genes. The DEGs were annotated with the AGA risk loci, and motif analysis in the promoters of the DEGs was also carried out. STRING Protein-protein interaction (PPI) and Reactome Functional Interaction (FI) networks were constructed using the DEGs, and the networks were analyzed to identify hub genes that play could play crucial roles in AGA pathogenesis.

Results and discussion

The in-silico study revealed that genes involved in the structural makeup of the skin epidermis, hair follicle development, and hair cycle are down-regulated, while genes associated with the innate and adaptive immune systems, cytokine signaling, and interferon signaling pathways are up-regulated in the balding scalps of AGA. The PPI and FI network analyses identified 25 hub genes namely CTNNB1, EGF, GNAI3, NRAS, BTK, ESR1, HCK, ITGB7, LCK, LCP2, LYN, PDGFRB, PIK3CD, PTPN6, RAC2, SPI1, STAT3, STAT5A, VAV1, PSMB8, HLA-A, HLA-F, HLA-E, IRF4, and ITGAM that play crucial roles in AGA pathogenesis. The study also implicates that Src family tyrosine kinase genes such as LCK, and LYN in the up-regulation of the inflammatory process in the balding scalps of AGA highlighting their potential as therapeutic targets for future investigations.

Early feeding leads to molecular maturation of the gut mucosal immune system in suckling piglets

Introduction

Diet-microbiota-host interactions are increasingly studied to comprehend their implications in host metabolism and overall health. Keeping in mind the importance of early life programming in shaping intestinal mucosal development, the pre-weaning period can be utilised to understand these interactions in suckling piglets. The objective of this study was to investigate the consequences of early life feeding on the time-resolved mucosal transcriptional program as well as mucosal morphology.

Methods

A customised fibrous feed was provided to piglets (early-fed or EF group; 7 litters) from five days of age until weaning (29 days of age) in addition to sow's milk, whereas control piglets (CON; 6 litters) suckled mother's milk only. Rectal swabs, intestinal content, and mucosal tissues (jejunum, colon) were obtained pre- and post-weaning for microbiota analysis (16S amplicon sequencing) and host transcriptome analysis (RNA sequencing).

Results

Early feeding accelerated both microbiota colonisation as well as host transcriptome, towards a more "mature state", with a more pronounced response in colon compared to jejunum. Early feeding elicited the largest impact on the colon transcriptome just before weaning (compared to post-weaning time-points), exemplified by the modulation of genes involved in cholesterol and energy metabolism and immune response. The transcriptional impact of early feeding persisted during the first days post-weaning and was highlighted by a stronger mucosal response to the weaning stress, via pronounced activation of barrier repair reactions, which is a combination of immune activation, epithelial migration and "wound-repair" like processes, compared to the CON piglets.

Discussion

Our study demonstrates the potential of early life nutrition in neonatal piglets as a means to support their intestinal development during the suckling period, and to improve adaptation during the weaning transition.

Benchmarking of protein interaction databases for integration with manually reconstructed signalling network models

Protein interaction databases are critical resources for network bioinformatics and integrating molecular experimental data. Interaction databases may also enable construction of predictive computational models of biological networks, although their fidelity for this purpose is not clear. Here, we benchmark protein interaction databases X2K, Reactome, Pathway Commons, Omnipath and Signor for their ability to recover manually curated edges from three logic-based network models of cardiac hypertrophy, mechano-signalling and fibrosis. Pathway Commons performed best at recovering interactions from manually reconstructed hypertrophy (137 of 193 interactions, 71%), mechano-signalling (85 of 125 interactions, 68%) and fibroblast networks (98 of 142 interactions, 69%). While protein interaction databases successfully recovered central, well-conserved pathways, they performed worse at recovering tissue-specific and transcriptional regulation. This highlights a knowledge gap where manual curation is critical. Finally, we tested the ability of Signor and Pathway Commons to identify new edges that improve model predictions, revealing important roles of protein kinase C autophosphorylation and Ca2+ /calmodulin-dependent protein kinase II phosphorylation of CREB in cardiomyocyte hypertrophy. This study provides a platform for benchmarking protein interaction databases for their utility in network model construction, as well as providing new insights into cardiac hypertrophy signalling. KEY POINTS: Protein interaction databases are used to recover signalling interactions from previously developed network models. The five protein interaction databases benchmarked recovered well-conserved pathways, but did poorly at recovering tissue-specific pathways and transcriptional regulation, indicating the importance of manual curation. We identify new signalling interactions not previously used in the network models, including a role for Ca2+ /calmodulin-dependent protein kinase II phosphorylation of CREB in cardiomyocyte hypertrophy.

Unsupervised clustering of venous thromboembolism patients by clinical features at presentation identifies novel endotypes that improve prognostic stratification

BACKGROUND

Individuals with acute venous thromboembolism (VTE) constitute a heterogeneous group of patients with diverse clinical characteristics and outcome.

OBJECTIVES

To identify endotypes of individuals with acute VTE based on clinical characteristics at presentation through unsupervised cluster analysis and to evaluate their molecular proteomic profile and clinical outcome.

METHODS

Data from 591 individuals from the Genotyping and Molecular phenotyping of Venous thromboembolism (GMP-VTE) project were explored. Hierarchical clustering was applied to 58 variables to define VTE endotypes. Clinical characteristics, three-year incidence of thromboembolic events or death, and acute-phase plasma proteomics were assessed.

RESULTS

Four endotypes were identified, exhibiting different patterns of clinical characteristics and clinical course. Endotype 1 (n = 300), comprising older individuals with comorbidities, had the highest incidence of thromboembolic events or death (HR [95 % CI]: 3.76 [1.96-7.19]), followed by endotype 4 (n = 127) (HR [95 % CI]: 2.55 [1.26-5.16]), characterised by men with history of VTE and provoking risk factors, and endotype 3 (n = 57) (HR [95 % CI]: 1.57 [0.63-3.87]), composed of young women with provoking risk factors, vs. reference endotype 2 (n = 107). The reference endotype was constituted by individuals diagnosed with PE without comorbidities, who had the lowest incidence of the investigated endpoint. Differentially expressed proteins associated with the endotypes were related to distinct biological processes, supporting differences in molecular pathophysiology. The endotypes had superior prognostic ability compared to existing risk stratifications such as provoked vs unprovoked VTE and D-dimer levels.

CONCLUSION

Four endotypes of VTE were identified by unsupervised phenotype-based clustering that diverge in clinical outcome and plasmatic protein signature. This approach might support the future development of individualized treatment in VTE.

Reconstruction Set Test (RESET): a computationally efficient method for single sample gene set testing based on randomized reduced rank reconstruction error

We have developed a new, and analytically novel, single sample gene set testing method called Reconstruction Set Test (RESET). RESET quantifies gene set importance at both the sample-level and for the entire dataset based on the ability of set genes to reconstruct values for all measured genes. RESET addresses four important limitations of current techniques: 1) existing single sample methods are designed to detect mean differences and struggle to identify differential correlation patterns, 2) computationally efficient techniques are self-contained methods and cannot directly detect competitive scenarios where set genes differ from non-set genes in the same sample, 3) the scores generated by current methods can only be accurately compared across samples for a single set and not between sets, and 4) the computational performance of even the fastest existing methods be significant on very large datasets. RESET is realized using a computationally efficient randomized reduced rank reconstruction algorithm (available via the RESET R package on CRAN) that can effectively detect patterns of differential abundance and differential correlation for self-contained and competitive scenarios. As demonstrated using real and simulated scRNA-seq data, RESET provides superior accuracy at a lower computational cost relative to other single sample approaches.

Integrated host/microbe metagenomics enables accurate lower respiratory tract infection diagnosis in critically ill children

BACKGROUNDLower respiratory tract infection (LRTI) is a leading cause of death in children worldwide. LRTI diagnosis is challenging because noninfectious respiratory illnesses appear clinically similar and because existing microbiologic tests are often falsely negative or detect incidentally carried microbes, resulting in antimicrobial overuse and adverse outcomes. Lower airway metagenomics has the potential to detect host and microbial signatures of LRTI. Whether it can be applied at scale and in a pediatric population to enable improved diagnosis and treatment remains unclear.METHODSWe used tracheal aspirate RNA-Seq to profile host gene expression and respiratory microbiota in 261 children with acute respiratory failure. We developed a gene expression classifier for LRTI by training on patients with an established diagnosis of LRTI (n = 117) or of noninfectious respiratory failure (n = 50). We then developed a classifier that integrates the host LRTI probability, abundance of respiratory viruses, and dominance in the lung microbiome of bacteria/fungi considered pathogenic by a rules-based algorithm.RESULTSThe host classifier achieved a median AUC of 0.967 by cross-validation, driven by activation markers of T cells, alveolar macrophages, and the interferon response. The integrated classifier achieved a median AUC of 0.986 and increased the confidence of patient classifications. When applied to patients with an uncertain diagnosis (n = 94), the integrated classifier indicated LRTI in 52% of cases and nominated likely causal pathogens in 98% of those.CONCLUSIONLower airway metagenomics enables accurate LRTI diagnosis and pathogen identification in a heterogeneous cohort of critically ill children through integration of host, pathogen, and microbiome features.FUNDINGSupport for this study was provided by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and the National Heart, Lung, and Blood Institute (UG1HD083171, 1R01HL124103, UG1HD049983, UG01HD049934, UG1HD083170, UG1HD050096, UG1HD63108, UG1HD083116, UG1HD083166, UG1HD049981, K23HL138461, and 5R01HL155418) as well as by the Chan Zuckerberg Biohub.

Comprehensive proteogenomic characterization of early duodenal cancer reveals the carcinogenesis tracks of different subtypes

The subtypes of duodenal cancer (DC) are complicated and the carcinogenesis process is not well characterized. We present comprehensive characterization of 438 samples from 156 DC patients, covering 2 major and 5 rare subtypes. Proteogenomics reveals LYN amplification at the chromosome 8q gain functioned in the transmit from intraepithelial neoplasia phase to infiltration tumor phase via MAPK signaling, and illustrates the DST mutation improves mTOR signaling in the duodenal adenocarcinoma stage. Proteome-based analysis elucidates stage-specific molecular characterizations and carcinogenesis tracks, and defines the cancer-driving waves of the adenocarcinoma and Brunner's gland subtypes. The drug-targetable alanyl-tRNA synthetase (AARS1) in the high tumor mutation burden/immune infiltration is significantly enhanced in DC progression, and catalyzes the lysine-alanylation of poly-ADP-ribose polymerases (PARP1), which decreases the apoptosis of cancer cells, eventually promoting cell proliferation and tumorigenesis. We assess the proteogenomic landscape of early DC, and provide insights into the molecular features corresponding therapeutic targets.

Integrative proteogenomic characterization of early esophageal cancer

Esophageal squamous cell carcinoma (ESCC) is malignant while the carcinogenesis is still unclear. Here, we perform a comprehensive multi-omics analysis of 786 trace-tumor-samples from 154 ESCC patients, covering 9 histopathological stages and 3 phases. Proteogenomics elucidates cancer-driving waves in ESCC progression, and reveals the molecular characterization of alcohol drinking habit associated signatures. We discover chromosome 3q gain functions in the transmit from nontumor to intraepithelial neoplasia phases, and find TP53 mutation enhances DNA replication in intraepithelial neoplasia phase. The mutations of AKAP9 and MCAF1 upregulate glycolysis and Wnt signaling, respectively, in advanced-stage ESCC phase. Six major tracks related to different clinical features during ESCC progression are identified, which is validated by an independent cohort with another 256 samples. Hyperphosphorylated phosphoglycerate kinase 1 (PGK1, S203) is considered as a drug target in ESCC progression. This study provides insight into the understanding of ESCC molecular mechanism and the development of therapeutic targets.

NaRnEA: An Information Theoretic Framework for Gene Set Analysis

Gene sets are being increasingly leveraged to make high-level biological inferences from transcriptomic data; however, existing gene set analysis methods rely on overly conservative, heuristic approaches for quantifying the statistical significance of gene set enrichment. We created Nonparametric analytical-Rank-based Enrichment Analysis (NaRnEA) to facilitate accurate and robust gene set analysis with an optimal null model derived using the information theoretic Principle of Maximum Entropy. By measuring the differential activity of ~2500 transcriptional regulatory proteins based on the differential expression of each protein's transcriptional targets between primary tumors and normal tissue samples in three cohorts from The Cancer Genome Atlas (TCGA), we demonstrate that NaRnEA critically improves in two widely used gene set analysis methods: Gene Set Enrichment Analysis (GSEA) and analytical-Rank-based Enrichment Analysis (aREA). We show that the NaRnEA-inferred differential protein activity is significantly correlated with differential protein abundance inferred from independent, phenotype-matched mass spectrometry data in the Clinical Proteomic Tumor Analysis Consortium (CPTAC), confirming the statistical and biological accuracy of our approach. Additionally, our analysis crucially demonstrates that the sample-shuffling empirical null models leveraged by GSEA and aREA for gene set analysis are overly conservative, a shortcoming that is avoided by the newly developed Maximum Entropy analytical null model employed by NaRnEA.

The formation and function of the neutrophil phagosome

Neutrophils are the most abundant circulating leukocyte and are crucial to the initial innate immune response to infection. One of their key pathogen-eliminating mechanisms is phagocytosis, the process of particle engulfment into a vacuole-like structure called the phagosome. The antimicrobial activity of the phagocytic process results from a collaboration of multiple systems and mechanisms within this organelle, where a complex interplay of ion fluxes, pH, reactive oxygen species, and antimicrobial proteins creates a dynamic antimicrobial environment. This complexity, combined with the difficulties of studying neutrophils ex vivo, has led to gaps in our knowledge of how the neutrophil phagosome optimizes pathogen killing. In particular, controversy has arisen regarding the relative contribution and integration of nicotinamide adenine dinucleotide phosphate (NADPH) oxidase-derived antimicrobial agents and granule-delivered antimicrobial proteins. Clinical syndromes arising from dysfunction in these systems in humans allow useful insight into these mechanisms, but their redundancy and synergy add to the complexity. In this article, we review the current knowledge regarding the formation and function of the neutrophil phagosome, examine new insights into the phagosomal environment that have been permitted by technological advances in recent years, and discuss aspects of the phagocytic process that are still under debate.

Analysis of genetically determined gene expression suggests role of inflammatory processes in exfoliation syndrome

BACKGROUND

Exfoliation syndrome (XFS) is an age-related systemic disorder characterized by excessive production and progressive accumulation of abnormal extracellular material, with pathognomonic ocular manifestations. It is the most common cause of secondary glaucoma, resulting in widespread global blindness. The largest global meta-analysis of XFS in 123,457 multi-ethnic individuals from 24 countries identified seven loci with the strongest association signal in chr15q22-25 region near LOXL1. Expression analysis have so far correlated coding and a few non-coding variants in the region with LOXL1 expression levels, but functional effects of these variants is unclear. We hypothesize that analysis of the contribution of the genetically determined component of gene expression to XFS risk can provide a powerful method to elucidate potential roles of additional genes and clarify biology that underlie XFS.

RESULTS

Transcriptomic Wide Association Studies (TWAS) using PrediXcan models trained in 48 GTEx tissues leveraging on results from the multi-ethnic and European ancestry GWAS were performed. To eliminate the possibility of false-positive results due to Linkage Disequilibrium (LD) contamination, we i) performed PrediXcan analysis in reduced models removing variants in LD with LOXL1 missense variants associated with XFS, and variants in LOXL1 models in both multiethnic and European ancestry individuals, ii) conducted conditional analysis of the significant signals in European ancestry individuals, and iii) filtered signals based on correlated gene expression, LD and shared eQTLs, iv) conducted expression validation analysis in human iris tissues. We observed twenty-eight genes in chr15q22-25 region that showed statistically significant associations, which were whittled down to ten genes after statistical validations. In experimental analysis, mRNA transcript levels for ARID3B, CD276, LOXL1, NEO1, SCAMP2, and UBL7 were significantly decreased in iris tissues from XFS patients compared to control samples. TWAS genes for XFS were significantly enriched for genes associated with inflammatory conditions. We also observed a higher incidence of XFS comorbidity with inflammatory and connective tissue diseases.

CONCLUSION

Our results implicate a role for connective tissues and inflammation pathways in the etiology of XFS. Targeting the inflammatory pathway may be a potential therapeutic option to reduce progression in XFS.

Multilingual translation for zero-shot biomedical classification using BioTranslator

Existing annotation paradigms rely on controlled vocabularies, where each data instance is classified into one term from a predefined set of controlled vocabularies. This paradigm restricts the analysis to concepts that are known and well-characterized. Here, we present the novel multilingual translation method BioTranslator to address this problem. BioTranslator takes a user-written textual description of a new concept and then translates this description to a non-text biological data instance. The key idea of BioTranslator is to develop a multilingual translation framework, where multiple modalities of biological data are all translated to text. We demonstrate how BioTranslator enables the identification of novel cell types using only a textual description and how BioTranslator can be further generalized to protein function prediction and drug target identification. Our tool frees scientists from limiting their analyses within predefined controlled vocabularies, enabling them to interact with biological data using free text.

Transcriptome-Based Traits of Radioresistant Sublines of Non-Small Cell Lung Cancer Cells

Radioresistance is a major obstacle for the successful therapy of many cancers, including non-small cell lung cancer (NSCLC). To elucidate the mechanism of radioresistance of NSCLC cells and to identify key molecules conferring radioresistance, the radioresistant subclones of p53 wild-type A549 and p53-deficient H1299 cell cultures were established. The transcriptional changes between parental and radioresistant NSCLC cells were investigated by RNA-seq. In total, expression levels of 36,596 genes were measured. Changes in the activation of intracellular molecular pathways of cells surviving irradiation relative to parental cells were quantified using the Oncobox bioinformatics platform. Following 30 rounds of 2 Gy irradiation, a total of 322 genes were differentially expressed between p53 wild-type radioresistant A549IR and parental A549 cells. For the p53-deficient (H1299) NSCLC cells, the parental and irradiated populations differed in the expression of 1628 genes and 1616 pathways. The expression of genes associated with radioresistance reflects the complex biological processes involved in clinical cancer cell eradication and might serve as a potential biomarker and therapeutic target for NSCLC treatment.

Data synthesis and adversarial networks: A review and meta-analysis in cancer imaging

Despite technological and medical advances, the detection, interpretation, and treatment of cancer based on imaging data continue to pose significant challenges. These include inter-observer variability, class imbalance, dataset shifts, inter- and intra-tumour heterogeneity, malignancy determination, and treatment effect uncertainty. Given the recent advancements in image synthesis, Generative Adversarial Networks (GANs), and adversarial training, we assess the potential of these technologies to address a number of key challenges of cancer imaging. We categorise these challenges into (a) data scarcity and imbalance, (b) data access and privacy, (c) data annotation and segmentation, (d) cancer detection and diagnosis, and (e) tumour profiling, treatment planning and monitoring. Based on our analysis of 164 publications that apply adversarial training techniques in the context of cancer imaging, we highlight multiple underexplored solutions with research potential. We further contribute the Synthesis Study Trustworthiness Test (SynTRUST), a meta-analysis framework for assessing the validation rigour of medical image synthesis studies. SynTRUST is based on 26 concrete measures of thoroughness, reproducibility, usefulness, scalability, and tenability. Based on SynTRUST, we analyse 16 of the most promising cancer imaging challenge solutions and observe a high validation rigour in general, but also several desirable improvements. With this work, we strive to bridge the gap between the needs of the clinical cancer imaging community and the current and prospective research on data synthesis and adversarial networks in the artificial intelligence community.