The Network of Cancer Genes (NCG): a comprehensive catalogue of known and candidate cancer genes from cancer sequencing screens

A survey of k-mer methods and applications in bioinformatics

The rapid progression of genomics and proteomics has been driven by the advent of advanced sequencing technologies, large, diverse, and readily available omics datasets, and the evolution of computational data processing capabilities. The vast amount of data generated by these advancements necessitates efficient algorithms to extract meaningful information. K-mers serve as a valuable tool when working with large sequencing datasets, offering several advantages in computational speed and memory efficiency and carrying the potential for intrinsic biological functionality. This review provides an overview of the methods, applications, and significance of k-mers in genomic and proteomic data analyses, as well as the utility of absent sequences, including nullomers and nullpeptides, in disease detection, vaccine development, therapeutics, and forensic science. Therefore, the review highlights the pivotal role of k-mers in addressing current genomic and proteomic problems and underscores their potential for future breakthroughs in research.

scSID: A lightweight algorithm for identifying rare cell types by capturing differential expression from single-cell sequencing data

Single-cell RNA sequencing (scRNA-seq) is currently an important technology for identifying cell types and studying diseases at the genetic level. Identifying rare cell types is biologically important as one of the downstream data analyses of single-cell RNA sequencing. Although rare cell identification methods have been developed, most of these suffer from insufficient mining of intercellular similarities, low scalability, and being time-consuming. In this paper, we propose a single-cell similarity division algorithm (scSID) for identifying rare cells. It takes cell-to-cell similarity into consideration by analyzing both inter-cluster and intra-cluster similarities, and discovers rare cell types based on the similarity differences. We show that scSID outperforms other existing methods by benchmarking it on different experimental datasets. Application of scSID to multiple datasets, including 68K PBMC and intestine, highlights its exceptional scalability and remarkable ability to identify rare cell populations.

Clinical use of tigecycline may contribute to the widespread dissemination of carbapenem-resistant hypervirulent Klebsiella pneumoniae strains

The emergence of carbapenem-resistant hypervirulent Klebsiella pneumoniae (CR-hvKP) poses grave threats to human health. These strains increased dramatically in clinical settings in China in the past few years but not in other parts of the world. Four isogenic K. pneumoniae strains, including classical K. pneumoniae, carbapenem-resistant K. pneumoniae (CRKP), hypervirulent K. pneumoniae (hvKP) and CR-hvKP, were created and subjected to phenotypic characterization, competition assays, mouse sepsis model and rat colonization tests to investigate the mechanisms underlying the widespread nature of CR-hvKP in China. Acquisition of virulence plasmid led to reduced fitness and abolishment of colonization in the gastrointestinal tract, which may explain why hvKP is not clinically prevalent after its emergence for a long time. However, tigecycline treatment facilitated the colonization of hvKP and CR-hvKP and reduced the population of Lactobacillus spp. in animal gut microbiome. Feeding with Lactobacillus spp. could significantly reduce the colonization of hvKP and CR-hvKP in the animal gastrointestinal tract. Our data implied that the clinical use of tigecycline to treat carbapenem-resistant K. pneumoniae infections facilitated the high spread of CR-hvKP in clinical settings in China and demonstrated that Lactobacillus spp. was a potential candidate for anticolonization strategy against CR-hvKP.

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans - a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes' cellular and molecular biology - which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus - a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates.

Disrupted gut harmony in attention-deficit/hyperactivity disorder: Dysbiosis and decreased short-chain fatty acids

Background

Attention-Deficit Hyperactivity Disorder (ADHD) is a prevalent neurodevelopmental disorder with complex genetic and environmental underpinnings. Emerging evidence suggests a significant role of gut microbiota in ADHD pathophysiology. This study investigates variations in gut microbiota composition and Short-Chain Fatty Acid (SCFA) profiles between children and adolescents with ADHD and healthy controls.

Methods

The study included 42 ADHD patients and 31 healthy controls, aged 6-18 years. Fecal samples were analyzed for microbial composition using 16S rRNA gene sequencing and for SCFA profiles through gas chromatography-mass spectrometry (GC-MS). The study assessed both α and β diversity of gut microbiota and quantified various SCFAs to compare between the groups.

Results

ADHD subjects demonstrated significantly reduced gut microbiota diversity, as indicated by lower α-diversity indices (Shannon index, Observed species, Faith PD index) and a trend towards significance in β-diversity (Weighted UniFrac). Notably, the ADHD group exhibited significantly lower levels of key SCFAs, including acetic, propionic, isobutyric, isovaleric, and valeric acids, highlighting a distinct microbial and metabolic profile in these individuals.

Conclusion

This study uncovers significant alterations in gut microbiota and SCFA profiles in children with ADHD, compared to healthy controls. The observed changes in SCFAs, known for their associations with other behavioral and neurologic pathologies, and for their role in neural signaling. These findings offer a metabolite fingerprint that could potentially lead to novel diagnostic and treatment approaches for ADHD, emphasizing the importance of gut microbiota in the disorder's pathogenesis and management.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

A Functional Survey of the Regulatory Landscape of Estrogen Receptor-Positive Breast Cancer Evolution

Only a handful of somatic alterations have been linked to endocrine therapy resistance in hormone-dependent breast cancer, potentially explaining ∼40% of relapses. If other mechanisms underlie the evolution of hormone-dependent breast cancer under adjuvant therapy is currently unknown. In this work, we employ functional genomics to dissect the contribution of cis-regulatory elements (CRE) to cancer evolution by focusing on 12 megabases of noncoding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Parallel epigenetic perturbation (CRISPRi) in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies a limited set of noncoding changes potentially involved in therapy resistance. Overall, our data uncover how endocrine therapies trigger the emergence of transient features which could ultimately be exploited to hinder the adaptive process. Significance: This study shows that cells adapting to endocrine therapies undergo changes in the usage or regulatory regions. Dormant cells are less vulnerable to regulatory perturbation but gain transient dependencies which can be exploited to decrease the formation of dormant persisters.

ECD-CDGI: An efficient energy-constrained diffusion model for cancer driver gene identification

The identification of cancer driver genes (CDGs) poses challenges due to the intricate interdependencies among genes and the influence of measurement errors and noise. We propose a novel energy-constrained diffusion (ECD)-based model for identifying CDGs, termed ECD-CDGI. This model is the first to design an ECD-Attention encoder by combining the ECD technique with an attention mechanism. ECD-Attention encoder excels at generating robust gene representations that reveal the complex interdependencies among genes while reducing the impact of data noise. We concatenate topological embedding extracted from gene-gene networks through graph transformers to these gene representations. We conduct extensive experiments across three testing scenarios. Extensive experiments show that the ECD-CDGI model possesses the ability to not only be proficient in identifying known CDGs but also efficiently uncover unknown potential CDGs. Furthermore, compared to the GNN-based approach, the ECD-CDGI model exhibits fewer constraints by existing gene-gene networks, thereby enhancing its capability to identify CDGs. Additionally, ECD-CDGI is open-source and freely available. We have also launched the model as a complimentary online tool specifically crafted to expedite research efforts focused on CDGs identification.

RAS-ON inhibition overcomes clinical resistance to KRAS G12C-OFF covalent blockade

Selective KRASG12C inhibitors have been developed to covalently lock the oncogene in the inactive GDP-bound state. Two of these molecules, sotorasib and adagrasib, are approved for the treatment of adult patients with KRASG12C-mutated previously treated advanced non-small cell lung cancer. Drug treatment imposes selective pressures leading to the outgrowth of drug-resistant variants. Mass sequencing from patients' biopsies identified a number of acquired KRAS mutations -both in cis and in trans- in resistant tumors. We demonstrate here that disease progression in vivo can also occur due to adaptive mechanisms and increased KRAS-GTP loading. Using the preclinical tool tri-complex KRASG12C-selective covalent inhibitor, RMC-4998 (also known as RM-029), that targets the active GTP-bound (ON) state of the oncogene, we provide a proof-of-concept that the clinical stage KRASG12C(ON) inhibitor RMC-6291 alone or in combination with KRASG12C(OFF) drugs can be an alternative potential therapeutic strategy to circumvent resistance due to increased KRAS-GTP loading.

Systematic discovery of gene-environment interactions underlying the human plasma proteome in UK Biobank

Understanding how gene-environment interactions (GEIs) influence the circulating proteome could aid in biomarker discovery and validation. The presence of GEIs can be inferred from single nucleotide polymorphisms that associate with phenotypic variability - termed variance quantitative trait loci (vQTLs). Here, vQTL association studies are performed on plasma levels of 1463 proteins in 52,363 UK Biobank participants. A set of 677 independent vQTLs are identified across 568 proteins. They include 67 variants that lack conventional additive main effects on protein levels. Over 1100 GEIs are identified between 101 proteins and 153 environmental exposures. GEI analyses uncover possible mechanisms that explain why 13/67 vQTL-only sites lack corresponding main effects. Additional analyses also highlight how age, sex, epistatic interactions and statistical artefacts may underscore associations between genetic variation and variance heterogeneity. This study establishes the most comprehensive database yet of vQTLs and GEIs for the human proteome.

Unraveling druggable cancer-driving proteins and targeted drugs using artificial intelligence and multi-omics analyses

The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .

Identifying cooperating cancer driver genes in individual patients through hypergraph random walk

OBJECTIVE

Identifying cancer driver genes, especially rare or patient-specific cancer driver genes, is a primary goal in cancer therapy. Although researchers have proposed some methods to tackle this problem, these methods mostly identify cancer driver genes at single gene level, overlooking the cooperative relationship among cancer driver genes. Identifying cooperating cancer driver genes in individual patients is pivotal for understanding cancer etiology and advancing the development of personalized therapies.

METHODS

Here, we propose a novel Personalized Cooperating cancer Driver Genes (PCoDG) method by using hypergraph random walk to identify the cancer driver genes that cooperatively drive individual patient cancer progression. By leveraging the powerful ability of hypergraph in representing multi-way relationships, PCoDG first employs the personalized hypergraph to depict the complex interactions among mutated genes and differentially expressed genes of an individual patient. Then, a hypergraph random walk algorithm based on hyperedge similarity is utilized to calculate the importance scores of mutated genes, integrating these scores with signaling pathway data to identify the cooperating cancer driver genes in individual patients.

RESULTS

The experimental results on three TCGA cancer datasets (i.e., BRCA, LUAD, and COADREAD) demonstrate the effectiveness of PCoDG in identifying personalized cooperating cancer driver genes. These genes identified by PCoDG not only offer valuable insights into patient stratification correlating with clinical outcomes, but also provide an useful reference resource for tailoring personalized treatments.

CONCLUSION

We propose a novel method that can effectively identify cooperating cancer driver genes for individual patients, thereby deepening our understanding of the cooperative relationship among personalized cancer driver genes and advancing the development of precision oncology.

Multiview representation learning for identification of novel cancer genes and their causative biological mechanisms

Tumorigenesis arises from the dysfunction of cancer genes, leading to uncontrolled cell proliferation through various mechanisms. Establishing a complete cancer gene catalogue will make precision oncology possible. Although existing methods based on graph neural networks (GNN) are effective in identifying cancer genes, they fall short in effectively integrating data from multiple views and interpreting predictive outcomes. To address these shortcomings, an interpretable representation learning framework IMVRL-GCN is proposed to capture both shared and specific representations from multiview data, offering significant insights into the identification of cancer genes. Experimental results demonstrate that IMVRL-GCN outperforms state-of-the-art cancer gene identification methods and several baselines. Furthermore, IMVRL-GCN is employed to identify a total of 74 high-confidence novel cancer genes, and multiview data analysis highlights the pivotal roles of shared, mutation-specific, and structure-specific representations in discriminating distinctive cancer genes. Exploration of the mechanisms behind their discriminative capabilities suggests that shared representations are strongly associated with gene functions, while mutation-specific and structure-specific representations are linked to mutagenic propensity and functional synergy, respectively. Finally, our in-depth analyses of these candidates suggest potential insights for individualized treatments: afatinib could counteract many mutation-driven risks, and targeting interactions with cancer gene SRC is a reasonable strategy to mitigate interaction-induced risks for NR3C1, RXRA, HNF4A, and SP1.

Using clusterProfiler to characterize multiomics data

With the advent of multiomics, software capable of multidimensional enrichment analysis has become increasingly crucial for uncovering gene set variations in biological processes and disease pathways. This is essential for elucidating disease mechanisms and identifying potential therapeutic targets. clusterProfiler stands out for its comprehensive utilization of databases and advanced visualization features. Importantly, clusterProfiler supports various biological knowledge, including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, through performing over-representation and gene set enrichment analyses. A key feature is that clusterProfiler allows users to choose from various graphical outputs to visualize results, enhancing interpretability. This protocol describes innovative ways in which clusterProfiler has been used for integrating metabolomics and metagenomics analyses, identifying and characterizing transcription factors under stress conditions, and annotating cells in single-cell studies. In all cases, the computational steps can be completed within ~2 min. clusterProfiler is released through the Bioconductor project and can be accessed via https://bioconductor.org/packages/clusterProfiler/ .

CGMega: explainable graph neural network framework with attention mechanisms for cancer gene module dissection

Cancer is rarely the straightforward consequence of an abnormality in a single gene, but rather reflects a complex interplay of many genes, represented as gene modules. Here, we leverage the recent advances of model-agnostic interpretation approach and develop CGMega, an explainable and graph attention-based deep learning framework to perform cancer gene module dissection. CGMega outperforms current approaches in cancer gene prediction, and it provides a promising approach to integrate multi-omics information. We apply CGMega to breast cancer cell line and acute myeloid leukemia (AML) patients, and we uncover the high-order gene module formed by ErbB family and tumor factors NRG1, PPM1A and DLG2. We identify 396 candidate AML genes, and observe the enrichment of either known AML genes or candidate AML genes in a single gene module. We also identify patient-specific AML genes and associated gene modules. Together, these results indicate that CGMega can be used to dissect cancer gene modules, and provide high-order mechanistic insights into cancer development and heterogeneity.

TOTEM: a multi-cancer detection and localization approach using circulating tumor DNA methylation markers

BACKGROUND

Detection of cancer and identification of tumor origin at an early stage improve the survival and prognosis of patients. Herein, we proposed a plasma cfDNA-based approach called TOTEM to detect and trace the cancer signal origin (CSO) through methylation markers.

METHODS

We performed enzymatic conversion-based targeted methylation sequencing on plasma cfDNA samples collected from a clinical cohort of 500 healthy controls and 733 cancer patients with seven types of cancer (breast, colorectum, esophagus, stomach, liver, lung, and pancreas) and randomly divided these samples into a training cohort and a testing cohort. An independent validation cohort of 143 healthy controls, 79 liver cancer patients and 100 stomach cancer patients were recruited to validate the generalizability of our approach.

RESULTS

A total of 57 multi-cancer diagnostic markers and 873 CSO markers were selected for model development. The binary diagnostic model achieved an area under the curve (AUC) of 0.907, 0.908 and 0.868 in the training, testing and independent validation cohorts, respectively. With a training specificity of 98%, the specificities in the testing and independent validation cohorts were 100% and 98.6%, respectively. Overall sensitivity across all cancer stages was 65.5%, 67.3% and 55.9% in the training, testing and independent validation cohorts, respectively. Early-stage (I and II) sensitivity was 50.3% and 45.7% in the training and testing cohorts, respectively. For cancer patients correctly identified by the binary classifier, the top 1 and top 2 CSO accuracies were 77.7% and 86.5% in the testing cohort (n = 148) and 76.0% and 84.0% in the independent validation cohort (n = 100). Notably, performance was maintained with only 21 diagnostic and 214 CSO markers, achieving a training AUC of 0.865, a testing AUC of 0.866, and an integrated top 2 accuracy of 83.1% in the testing cohort.

CONCLUSIONS

TOTEM demonstrates promising potential for accurate multi-cancer detection and localization by profiling plasma methylation markers. The real-world clinical performance of our approach needs to be investigated in a much larger prospective cohort.

Epigenetic-based differentiation therapy for Acute Myeloid Leukemia

Despite the development of novel therapies for acute myeloid leukemia, outcomes remain poor for most patients, and therapeutic improvements are an urgent unmet need. Although treatment regimens promoting differentiation have succeeded in the treatment of acute promyelocytic leukemia, their role in other acute myeloid leukemia subtypes needs to be explored. Here we identify and characterize two lysine deacetylase inhibitors, CM-444 and CM-1758, exhibiting the capacity to promote myeloid differentiation in all acute myeloid leukemia subtypes at low non-cytotoxic doses, unlike other commercial histone deacetylase inhibitors. Analyzing the acetylome after CM-444 and CM-1758 treatment reveals modulation of non-histone proteins involved in the enhancer-promoter chromatin regulatory complex, including bromodomain proteins. This acetylation is essential for enhancing the expression of key transcription factors directly involved in the differentiation therapy induced by CM-444/CM-1758 in acute myeloid leukemia. In summary, these compounds may represent effective differentiation-based therapeutic agents across acute myeloid leukemia subtypes with a potential mechanism for the treatment of acute myeloid leukemia.

Identifying new cancer genes based on the integration of annotated gene sets via hypergraph neural networks

MOTIVATION

Identifying cancer genes remains a significant challenge in cancer genomics research. Annotated gene sets encode functional associations among multiple genes, and cancer genes have been shown to cluster in hallmark signaling pathways and biological processes. The knowledge of annotated gene sets is critical for discovering cancer genes but remains to be fully exploited.

RESULTS

Here, we present the DIsease-Specific Hypergraph neural network (DISHyper), a hypergraph-based computational method that integrates the knowledge from multiple types of annotated gene sets to predict cancer genes. First, our benchmark results demonstrate that DISHyper outperforms the existing state-of-the-art methods and highlight the advantages of employing hypergraphs for representing annotated gene sets. Second, we validate the accuracy of DISHyper-predicted cancer genes using functional validation results and multiple independent functional genomics data. Third, our model predicts 44 novel cancer genes, and subsequent analysis shows their significant associations with multiple types of cancers. Overall, our study provides a new perspective for discovering cancer genes and reveals previously undiscovered cancer genes.

AVAILABILITY AND IMPLEMENTATION

DISHyper is freely available for download at https://github.com/genemine/DISHyper.

Prognostic significance of alpha-2-macrglobulin and low-density lipoprotein receptor-related protein-1 in various cancers

Cancer is the leading cause of death worldwide. The World Health Organization (WHO) estimates that 10 million fatalities occurred in 2023. Breast cancer (BC) ranked first among malignancies with 2.26 million cases, lung cancer (LC) second with 2.21 million cases, and colon and rectum cancers (CC, CRC) third with 1.93 million cases. These results highlight the importance of investigating novel cancer prognoses and anti-cancer markers. In this study, we investigated the potential effects of alpha-2 macroglobulin and its receptor, LRP1, on the outcomes of breast, lung, and colorectal malignancies. Immunohistochemical staining was used to analyze the expression patterns of A2M and LRP1 in 545 cases of invasive ductal breast carcinoma (IDC) and 51 cases of mastopathies/fibrocystic breast disease (FBD); 256 cases of non-small cell lung carcinomas (NSCLCs) and 45 cases of non-malignant lung tissue (NMLT); and 108 cases of CRC and 25 cases of non-malignant colorectal tissue (NMCT). A2M and LRP1 expression levels were also investigated in breast (MCF-7, BT-474, SK-BR-3, T47D, MDA-MB-231, and MDA-MB-231/BO2), lung (NCI-H1703, NCI-H522, and A549), and colon (LS 180, Caco-2, HT-29, and LoVo) cancer cell lines. Based on our findings, A2M and LRP1 exhibited various expression patterns in the examined malignancies, which were related to one another. Additionally, the stroma of lung and colorectal cancer has increased levels of A2M/LRP1 areas, which explains the significance of the stroma in the development and maintenance of tumor homeostasis. A2M expression was shown to be downregulated in all types of malignancies under study and was positively linked with an increase in cell line aggressiveness. Although more invasive cells had higher levels of A2M expression, an IHC analysis showed the opposite results. This might be because exogenous alpha-2-macroglobulin is present, which has an inhibitory effect on several cancerous enzymes and receptor-dependent signaling pathways. Additionally, siRNA-induced suppression of the transcripts for A2M and LRPP1 revealed their connection, which provides fresh information on the function of the LRP1 receptor in A2M recurrence in cancer. Further studies on different forms of cancer may corroborate the fact that both A2M and LRP1 have high potential as innovative therapeutic agents.

Detecting genetic gain and loss events in terms of protein domain: Method and implementation

Continuous gain and loss of genes are the primary driving forces of bacterial evolution and environmental adaptation. Studying bacterial evolution in terms of protein domain, which is the fundamental function and evolutionary unit of proteins, can provide a more comprehensive understanding of bacterial differentiation and phenotypic adaptation processes. Therefore, we proposed a phylogenetic tree-based method for detecting genetic gain and loss events in terms of protein domains. Specifically, the method focuses on a single domain to trace its evolution process or on multiple domains to investigate their co-evolution principles. This novel method was validated using 122 Shigella isolates. We found that the loss of a significant number of domains was likely the main driving force behind the evolution of Shigella, which could reduce energy expenditure and preserve only the most essential functions. Additionally, we observed that simultaneously gained and lost domains were often functionally related, which can facilitate and accelerate phenotypic evolutionary adaptation to the environment. All results obtained using our method agree with those of previous studies, which validates our proposed method.

CancerGATE: Prediction of cancer-driver genes using graph attention autoencoders

Discovery of the cancer type specific-driver genes is important for understanding the molecular mechanisms of each cancer type and for providing proper treatment. Recently, graph deep learning methods became widely used in finding cancer-driver genes. However, previous methods had limited performance in individual cancer types due to a small number of cancer-driver genes used in training and biases toward the cancer-driver genes used in training the models. Here, we introduce a novel pipeline, CancerGATE that predicts the cancer-driver genes using graph attention autoencoder (GATE) to learn in a self-supervised manner and can be applied to each of the cancer types. CancerGATE utilizes biological network topology and multi-omics data from 15 types of cancer of 20,079 samples from the cancer genome atlas (TCGA). Attention coefficients calculated in the model are used to prioritize cancer-driver genes by comparing coefficients of cancer and normal contexts. CancerGATE shows a higher AUPRC with a difference ranging from 1.5 % to 36.5 % compared to the previous graph deep learning models in each cancer type. We also show that CancerGATE is free from the bias toward cancer-driver genes used in training, revealing mechanisms of the cancer-driver genes in specific cancer types. Finally, we propose novel cancer-driver gene candidates that could be therapeutic targets for specific cancer types.

Sea cucumber (Holothuria glaberrima) intestinal microbiome dataset from Puerto Rico, generated by shotgun sequencing

The sea cucumber (H. glaberrima) is a species found in the shallow waters near coral reefs and seagrass beds in Puerto Rico. To characterize the microbial taxonomic composition and functional profiles present in the sea cucumber, total DNA was obtained from their intestinal system, fosmid libraries constructed, and subsequent sequencing was performed. The diversity profile displayed that the most predominant domain was Bacteria (76.56 %), followed by Viruses (23.24 %) and Archaea (0.04 %). Within the 11 phyla identified, the most abundant was Proteobacteria (73.16 %), followed by Terrabacteria group (3.20 %) and Fibrobacterota, Chlorobiota, Bacteroidota (FCB) superphylum (1.02 %). The most abundant species were Porvidencia rettgeri (21.77 %), Pseudomonas stutzeri (14.78 %), and Alcaligenes faecalis (5.00 %). The functional profile revealed that the most abundant functions are related to transporters, MISC (miscellaneous information systems), organic nitrogen, energy, and carbon utilization. The data collected in this project on the diversity and functional profiles of the intestinal system of the H. glaberrima provided a detailed view of its microbial ecology. These findings may motivate comparative studies aimed at understanding the role of the microbiome in intestinal regeneration.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Generalizable transcriptome-based tumor malignant level evaluation and molecular subtyping towards precision oncology

In cancer treatment, therapeutic strategies that integrate tumor-specific characteristics (i.e., precision oncology) are widely implemented to provide clinical benefits for cancer patients. Here, through in-depth integration of tumor transcriptome and patients' prognoses across cancers, we investigated dysregulated and prognosis-associated genes and catalogued such important genes in a cancer type-dependent manner. Utilizing the expression matrices of these genes, we built models to quantitatively evaluate the malignant levels of tumors across cancers, which could add value to the clinical staging system for improved prediction of patients' survival. Furthermore, we performed a transcriptome-based molecular subtyping on hepatocellular carcinoma, which revealed three subtypes with significantly diversified clinical outcomes, mutation landscapes, immune microenvironment, and dysregulated pathways. As tumor transcriptome was commonly profiled in clinical practice with low experimental complexity and cost, this work proposed easy-to-perform approaches for practical clinical promotion towards better healthcare and precision oncology of cancer patients.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

GLIMS: A two-stage gradual-learning method for cancer genes prediction using multi-omics data and co-splicing network

Identifying cancer genes is vital for cancer diagnosis and treatment. However, because of the complexity of cancer occurrence and limited cancer genes knowledge, it is hard to identify cancer genes accurately using only a few omics data, and the overall performance of existing methods is being called for further improvement. Here, we introduce a two-stage gradual-learning strategy GLIMS to predict cancer genes using integrative features from multi-omics data. Firstly, it uses a semi-supervised hierarchical graph neural network to predict the initial candidate cancer genes by integrating multi-omics data and protein-protein interaction (PPI) network. Then, it uses an unsupervised approach to further optimize the initial prediction by integrating the co-splicing network in post-transcriptional regulation, which plays an important role in cancer development. Systematic experiments on multi-omics cancer data demonstrated that GLIMS outperforms the state-of-the-art methods for the identification of cancer genes and it could be a useful tool to help advance cancer analysis.

The plant endomicrobiome: Structure and strategies to produce stress resilient future crop

Plants have a microbiome, a diverse community of microorganisms, including bacteria, fungi, and viruses, living inside and on their tissues. Versatile endophytic microorganisms inhabited in every plant part without causing disease and develop endophytic microbiome or endo-microbiome. Plant endo-microbiome are drawn by the nutrient rich micro-environment, and in turn some microbes mutualistically endorse and protect plant from adverse environmental stresses. Plant endo-microbiome interact within well-designed host equilibrium containing xylem, phloem, nutrients, phytohormones, metabolites and shift according to environmental and nutritional change. Plant endo-microbiome regulate and respond to environmental variations, pathogens, herbivores by producing stress regulators, organic acids, secondary metabolites, stress hormones as well as unknown substances and signalling molecules. Endomicrobiome efficiently synthesizes multiple bioactive compounds, stress phytohormones with high competence. The technological innovation as next generation genomics biology and high-throughput multiomics techniques stepping stones on the illumination of critical endo-microbiome communities and functional characterization that aid in improving plant physiology, biochemistry and immunity interplay for best crop productivity. This review article contains deeper insight in endomicrobiome related research work in last years, recruitment, niche development, nutrient dynamics, stress removal mechanisms, bioactive services in plant health development, community architecture and communication, and immunity interplay in producing stress resilient future crop.

High-quality chromosome-level genome assembly of Nicotiana benthamiana

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

Empowering prediction of miRNA-mRNA interactions in species with limited training data through transfer learning

MicroRNAs (miRNAs) play a crucial role in mRNA regulation. Identifying functionally important mRNA targets of a specific miRNA is essential for uncovering its biological function and assisting miRNA-based drug development. Datasets of high-throughput direct bona fide miRNA-target interactions (MTIs) exist only for a few model organisms, prompting the need for computational prediction. However, the scarcity of data poses a challenge in training accurate machine learning models for MTI prediction. In this study, we explored the potential of transfer learning technique (with ANN and XGB) to address the limited data challenge by leveraging the similarities in interaction rules between species. Furthermore, we introduced a novel approach called TransferSHAP for estimating the feature importance of transfer learning in tabular dataset tasks. We demonstrated that transfer learning improves MTI prediction accuracy for species with limited datasets and identified the specific interaction features the models employed to transfer information across different species.

Polyploidy Promotes Hypertranscription, Apoptosis Resistance, and Ciliogenesis in Cancer Cells and Mesenchymal Stem Cells of Various Origins: Comparative Transcriptome In Silico Study

Mesenchymal stem cells (MSC) attract an increasing amount of attention due to their unique therapeutic properties. Yet, MSC can undergo undesirable genetic and epigenetic changes during their propagation in vitro. In this study, we investigated whether polyploidy can compromise MSC oncological safety and therapeutic properties. For this purpose, we compared the impact of polyploidy on the transcriptome of cancer cells and MSC of various origins (bone marrow, placenta, and heart). First, we identified genes that are consistently ploidy-induced or ploidy-repressed through all comparisons. Then, we selected the master regulators using the protein interaction enrichment analysis (PIEA). The obtained ploidy-related gene signatures were verified using the data gained from polyploid and diploid populations of early cardiomyocytes (CARD) originating from iPSC. The multistep bioinformatic analysis applied to the cancer cells, MSC, and CARD indicated that polyploidy plays a pivotal role in driving the cell into hypertranscription. It was evident from the upregulation of gene modules implicated in housekeeping functions, stemness, unicellularity, DNA repair, and chromatin opening by means of histone acetylation operating via DNA damage associated with the NUA4/TIP60 complex. These features were complemented by the activation of the pathways implicated in centrosome maintenance and ciliogenesis and by the impairment of the pathways related to apoptosis, the circadian clock, and immunity. Overall, our findings suggest that, although polyploidy does not induce oncologic transformation of MSC, it might compromise their therapeutic properties because of global epigenetic changes and alterations in fundamental biological processes. The obtained results can contribute to the development and implementation of approaches enhancing the therapeutic properties of MSC by removing polyploid cells from the cell population.

Distinct genomic contexts predict gene presence-absence variation in different pathotypes of Magnaporthe oryzae

Fungi use the accessory gene content of their pangenomes to adapt to their environments. While gene presence-absence variation contributes to shaping accessory gene reservoirs, the genomic contexts that shape these events remain unclear. Since pangenome studies are typically species-wide and do not analyze different populations separately, it is yet to be uncovered whether presence-absence variation patterns and mechanisms are consistent across populations. Fungal plant pathogens are useful models for studying presence-absence variation because they rely on it to adapt to their hosts, and members of a species often infect distinct hosts. We analyzed gene presence-absence variation in the blast fungus, Magnaporthe oryzae (syn. Pyricularia oryzae), and found that presence-absence variation genes involved in host-pathogen and microbe-microbe interactions may drive the adaptation of the fungus to its environment. We then analyzed genomic and epigenomic features of presence-absence variation and observed that proximity to transposable elements, gene GC content, gene length, expression level in the host, and histone H3K27me3 marks were different between presence-absence variation genes and conserved genes. We used these features to construct a model that was able to predict whether a gene is likely to experience presence-absence variation with high precision (86.06%) and recall (92.88%) in M. oryzae. Finally, we found that presence-absence variation genes in the rice and wheat pathotypes of M. oryzae differed in their number and their genomic context. Our results suggest that genomic and epigenomic features of gene presence-absence variation can be used to better understand and predict fungal pangenome evolution. We also show that substantial intra-species variation can exist in these features.

Medieval DNA from Soqotra points to Eurasian origins of an isolated population at the crossroads of Africa and Arabia

Soqotra, an island situated at the mouth of the Gulf of Aden in the northwest Indian Ocean between Africa and Arabia, is home to ~60,000 people subsisting through fishing and semi-nomadic pastoralism who speak a Modern South Arabian language. Most of what is known about Soqotri history derives from writings of foreign travellers who provided little detail about local people, and the geographic origins and genetic affinities of early Soqotri people has not yet been investigated directly. Here we report genome-wide data from 39 individuals who lived between ~650 and 1750 CE at six locations across the island and document strong genetic connections between Soqotra and the similarly isolated Hadramawt region of coastal South Arabia that likely reflects a source for the peopling of Soqotra. Medieval Soqotri can be modelled as deriving ~86% of their ancestry from a population such as that found in the Hadramawt today, with the remaining ~14% best proxied by an Iranian-related source with up to 2% ancestry from the Indian sub-continent, possibly reflecting genetic exchanges that occurred along with archaeologically documented trade from these regions. In contrast to all other genotyped populations of the Arabian Peninsula, genome-level analysis of the medieval Soqotri is consistent with no sub-Saharan African admixture dating to the Holocene. The deep ancestry of people from medieval Soqotra and the Hadramawt is also unique in deriving less from early Holocene Levantine farmers and more from groups such as Late Pleistocene hunter-gatherers from the Levant (Natufians) than other mainland Arabians. This attests to migrations by early farmers having less impact in southernmost Arabia and Soqotra and provides compelling evidence that there has not been complete population replacement between the Pleistocene and Holocene throughout the Arabian Peninsula. Medieval Soqotra harboured a small population that showed qualitatively different marriage practices from modern Soqotri, with first-cousin unions occurring significantly less frequently than today.

Exploring risk factors and molecular targets in leukemia patients with COVID-19: a bioinformatics analysis of differential gene expression

The molecular mechanism of COVID-19's pathogenic effects in leukemia patients is still poorly known. Our study investigated the possible disease mechanism of COVID-19 and its associated risk factors in patients with leukemia utilizing differential gene expression analysis. We also employed network-based approaches to identify molecular targets that could potentially diagnose and treat COVID-19-infected leukemia patients. Our study demonstrated a shared set of 60 genes that are expressed differentially among patients with leukemia and COVID-19. Most of these genes are expressed in blood and bone marrow tissues and are predominantly implicated in the pathogenesis of different hematologic malignancies, increasingly imperiling COVID-19 morbidity and mortality among the affected patients. Additionally, we also found that COVID-19 may influence the expression of several cancer-associated genes in leukemia patients, such as CCR7, LEF1, and 13 candidate cancer-driver genes. Furthermore, our findings reveal that COVID-19 may predispose leukemia patients to altered blood homeostasis, increase the risk of COVID-19-related liver injury, and deteriorate leukemia-associated injury and patient prognosis. Our findings imply that molecular signatures, like transcription factors, proteins such as TOP21, and 25 different microRNAs, may be potential targets for diagnosing and treating COVID-19-infected leukemia patients. Nevertheless, additional experimental studies will contribute to further validating the study's findings.

Personalized Driver Gene Prediction Using Graph Convolutional Networks with Conditional Random Fields

Cancer is a complex and evolutionary disease mainly driven by the accumulation of genetic variations in genes. Identifying cancer driver genes is important. However, most related studies have focused on the population level. Cancer is a disease with high heterogeneity. Thus, the discovery of driver genes at the individual level is becoming more valuable but is a great challenge. Although there have been some computational methods proposed to tackle this challenge, few can cover all patient samples well, and there is still room for performance improvement. In this study, to identify individual-level driver genes more efficiently, we propose the PDGCN method. PDGCN integrates multiple types of data features, including mutation, expression, methylation, copy number data, and system-level gene features, along with network structural features extracted using Node2vec in order to construct a sample-gene interaction network. Prediction is performed using a graphical convolutional neural network model with a conditional random field layer, which is able to better combine the network structural features with biological attribute features. Experiments on the ACC (Adrenocortical Cancer) and KICH (Kidney Chromophobe) datasets from TCGA (The Cancer Genome Atlas) demonstrated that the method performs better compared to other similar methods. It can identify not only frequently mutated driver genes, but also rare candidate driver genes and novel biomarker genes. The results of the survival and enrichment analyses of these detected genes demonstrate that the method can identify important driver genes at the individual level.

Detecting and understanding meaningful cancerous mutations based on computational models of mRNA splicing

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.

CanVaxKB: a web-based cancer vaccine knowledgebase

Cancer vaccines have been increasingly studied and developed to prevent or treat various types of cancers. To systematically survey and analyze different reported cancer vaccines, we developed CanVaxKB (https://violinet.org/canvaxkb), the first web-based cancer vaccine knowledgebase that compiles over 670 therapeutic or preventive cancer vaccines that have been experimentally verified to be effective at various stages. Vaccine construction and host response data are also included. These cancer vaccines are developed against various cancer types such as melanoma, hematological cancer, and prostate cancer. CanVaxKB has stored 263 genes or proteins that serve as cancer vaccine antigen genes, which we have collectively termed 'canvaxgens'. Top three mostly used canvaxgens are PMEL, MLANA and CTAG1B, often targeting multiple cancer types. A total of 193 canvaxgens are also reported in cancer-related ONGene, Network of Cancer Genes and/or Sanger Cancer Gene Consensus databases. Enriched functional annotations and clusters of canvaxgens were identified and analyzed. User-friendly web interfaces are searchable for querying and comparing cancer vaccines. CanVaxKB cancer vaccines are also semantically represented by the community-based Vaccine Ontology to support data exchange. Overall, CanVaxKB is a timely and vital cancer vaccine source that facilitates efficient collection and analysis, further helping researchers and physicians to better understand cancer mechanisms.

An atlas of cell-type-specific interactome networks across 44 human tumor types

BACKGROUND

Biological processes are controlled by groups of genes acting in concert. Investigating gene-gene interactions within different cell types can help researchers understand the regulatory mechanisms behind human complex diseases, such as tumors.

METHODS

We collected extensive single-cell RNA-seq data from tumors, involving 563 patients with 44 different tumor types. Through our analysis, we identified various cell types in tumors and created an atlas of different immune cell subsets across different tumor types. Using the SCINET method, we reconstructed interactome networks specific to different cell types. Diverse functional data was then integrated to gain biological insights into the networks, including somatic mutation patterns and gene functional annotation. Additionally, genes with prognostic relevance within the networks were also identified. We also examined cell-cell communications to investigate how gene interactions modulate cell-cell interactions.

RESULTS

We developed a data portal called CellNetdb for researchers to study cell-type-specific interactome networks. Our findings indicate that these networks can be used to identify genes with topological specificity in different cell types. We also found that prognostic genes can deconvolved into cell types through analyzing network connectivity. Additionally, we identified commonalities and differences in cell-type-specific networks across different tumor types. Our results suggest that these networks can be used to prioritize risk genes.

CONCLUSIONS

This study presented CellNetdb, a comprehensive repository featuring an atlas of cell-type-specific interactome networks across 44 human tumor types. The findings underscore the utility of these networks in delineating the intricacies of tumor microenvironments and advancing the understanding of molecular mechanisms underpinning human tumors.

magpie: A power evaluation method for differential RNA methylation analysis in N6-methyladenosine sequencing

Recently, novel biotechnologies to quantify RNA modifications became an increasingly popular choice for researchers who study epitranscriptome. When studying RNA methylations such as N6-methyladenosine (m6A), researchers need to make several decisions in its experimental design, especially the sample size and a proper statistical power. Due to the complexity and high-throughput nature of m6A sequencing measurements, methods for power calculation and study design are still currently unavailable. In this work, we propose a statistical power assessment tool, magpie, for power calculation and experimental design for epitranscriptome studies using m6A sequencing data. Our simulation-based power assessment tool will borrow information from real pilot data, and inspect various influential factors including sample size, sequencing depth, effect size, and basal expression ranges. We integrate two modules in magpie: (i) a flexible and realistic simulator module to synthesize m6A sequencing data based on real data; and (ii) a power assessment module to examine a set of comprehensive evaluation metrics.

Exploring gene-patient association to identify personalized cancer driver genes by linear neighborhood propagation

BACKGROUND

Driver genes play a vital role in the development of cancer. Identifying driver genes is critical for diagnosing and understanding cancer. However, challenges remain in identifying personalized driver genes due to tumor heterogeneity of cancer. Although many computational methods have been developed to solve this problem, few efforts have been undertaken to explore gene-patient associations to identify personalized driver genes.

RESULTS

Here we propose a method called LPDriver to identify personalized cancer driver genes by employing linear neighborhood propagation model on individual genetic data. LPDriver builds personalized gene network based on the genetic data of individual patients, extracts the gene-patient associations from the bipartite graph of the personalized gene network and utilizes a linear neighborhood propagation model to mine gene-patient associations to detect personalized driver genes. The experimental results demonstrate that as compared to the existing methods, our method shows competitive performance and can predict cancer driver genes in a more accurate way. Furthermore, these results also show that besides revealing novel driver genes that have been reported to be related with cancer, LPDriver is also able to identify personalized cancer driver genes for individual patients by their network characteristics even if the mutation data of genes are hidden.

CONCLUSIONS

LPDriver can provide an effective approach to predict personalized cancer driver genes, which could promote the diagnosis and treatment of cancer. The source code and data are freely available at https://github.com/hyr0771/LPDriver .

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

Robustness evaluations of pathway activity inference methods on gene expression data

BACKGROUND

With the exponential growth of high-throughput technologies, multiple pathway analysis methods have been proposed to estimate pathway activities from gene expression profiles. These pathway activity inference methods can be divided into two main categories: non-Topology-Based (non-TB) and Pathway Topology-Based (PTB) methods. Although some review and survey articles discussed the topic from different aspects, there is a lack of systematic assessment and comparisons on the robustness of these approaches.

RESULTS

Thus, this study presents comprehensive robustness evaluations of seven widely used pathway activity inference methods using six cancer datasets based on two assessments. The first assessment seeks to investigate the robustness of pathway activity in pathway activity inference methods, while the second assessment aims to assess the robustness of risk-active pathways and genes predicted by these methods. The mean reproducibility power and total number of identified informative pathways and genes were evaluated. Based on the first assessment, the mean reproducibility power of pathway activity inference methods generally decreased as the number of pathway selections increased. Entropy-based Directed Random Walk (e-DRW) distinctly outperformed other methods in exhibiting the greatest reproducibility power across all cancer datasets. On the other hand, the second assessment shows that no methods provide satisfactory results across datasets.

CONCLUSION

However, PTB methods generally appear to perform better in producing greater reproducibility power and identifying potential cancer markers compared to non-TB methods.

Pathway analysis through mutual information

MOTIVATION

In pathway analysis, we aim to establish a connection between the activity of a particular biological pathway and a difference in phenotype. There are many available methods to perform pathway analysis, many of them rely on an upstream differential expression analysis, and many model the relations between the abundances of the analytes in a pathway as linear relationships.

RESULTS

Here, we propose a new method for pathway analysis, MIPath, that relies on information theoretical principles and, therefore, does not model the association between pathway activity and phenotype, resulting in relatively few assumptions. For this, we construct a graph of the data points for each pathway using a nearest-neighbor approach and score the association between the structure of this graph and the phenotype of these same samples using Mutual Information while adjusting for the effects of random chance in each score. The initial nearest neighbor approach evades individual gene-level comparisons, hence making the method scalable and less vulnerable to missing values. These properties make our method particularly useful for single-cell data. We benchmarked our method on several single-cell datasets, comparing it to established and new methods, and found that it produces robust, reproducible, and meaningful scores.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/statisticalbiotechnology/mipath, or through Python Package Index as "mipathway."

PRADclass: Hybrid Gleason Grade-Informed Computational Strategy Identifies Consensus Biomarker Features Predictive of Aggressive Prostate Adenocarcinoma

BACKGROUND

Prostate adenocarcinoma (PRAD) is a common cancer diagnosis among men globally, yet large gaps in our knowledge persist with respect to the molecular bases of its progression and aggression. It is mostly indolent and slow-growing, but aggressive prostate cancers need to be recognized early for optimising treatment, with a view to reducing mortality.

METHODS

Based on TCGA transcriptomic data pertaining to PRAD and the associated clinical metadata, we determined the sample Gleason grade, and used it to execute: (i) Gleason-grade wise linear modeling, followed by five contrasts against controls and ten contrasts between grades; and (ii) Gleason-grade wise network modeling via weighted gene correlation network analysis (WGCNA). Candidate biomarkers were obtained from the above analysis and the consensus found. The consensus biomarkers were used as the feature space to train ML models for classifying a sample as benign, indolent or aggressive.

RESULTS

The statistical modeling yielded 77 Gleason grade-salient genes while the WGCNA algorithm yielded 1003 trait-specific key genes in grade-wise significant modules. Consensus analysis of the two approaches identified two genes in Grade-1 (SLC43A1 and PHGR1), 26 genes in Grade-4 (including LOC100128675, PPP1R3C, NECAB1, UBXN10, SERPINA5, CLU, RASL12, DGKG, FHL1, NCAM1, and CEND1), and seven genes in Grade-5 (CBX2, DPYS, FAM72B, SHCBP1, TMEM132A, TPX2, UBE2C). A RandomForest model trained and optimized on these 35 biomarkers for the ternary classification problem yielded a balanced accuracy ∼ 86% on external validation.

CONCLUSIONS

The consensus of multiple parallel computational strategies has unmasked candidate Gleason grade-specific biomarkers. PRADclass, a validated AI model featurizing these biomarkers achieved good performance, and could be trialed to predict the differentiation of prostate cancers. PRADclass is available for academic use at: https://apalania.shinyapps.io/pradclass (online) and https://github.com/apalania/pradclass (command-line interface).

KombOver: Efficient k-core and K-truss based characterization of perturbations within the human gut microbiome

The microbes present in the human gastrointestinal tract are regularly linked to human health and disease outcomes. Thanks to technological and methodological advances in recent years, metagenomic sequencing data, and computational methods designed to analyze metagenomic data, have contributed to improved understanding of the link between the human gut microbiome and disease. However, while numerous methods have been recently developed to extract quantitative and qualitative results from host-associated microbiome data, improved computational tools are still needed to track microbiome dynamics with short-read sequencing data. Previously we have proposed KOMB as a de novo tool for identifying copy number variations in metagenomes for characterizing microbial genome dynamics in response to perturbations. In this work, we present KombOver (KO), which includes four key contributions with respect to our previous work: (i) it scales to large microbiome study cohorts, (ii) it includes both k-core and K-truss based analysis, (iii) we provide the foundation of a theoretical understanding of the relation between various graph-based metagenome representations, and (iv) we provide an improved user experience with easier-to-run code and more descriptive outputs/results. To highlight the aforementioned benefits, we applied KO to nearly 1000 human microbiome samples, requiring less than 10 minutes and 10 GB RAM per sample to process these data. Furthermore, we highlight how graph-based approaches such as k-core and K-truss can be informative for pinpointing microbial community dynamics within a myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) cohort. KO is open source and available for download/use at: https://github.com/treangenlab/komb.

A Interacting Model: How TRIM21 Orchestrates with Proteins in Intracellular Immunity

Tripartite motif-containing protein 21 (TRIM21), identified as both a cytosolic E3 ubiquitin ligase and FcR (Fragment crystallizable receptor), primarily interacts with proteins via its PRY/SPRY domains and promotes their proteasomal degradation to regulate intracellular immunity. But how TRIM21 involves in intracellular immunity still lacks systematical understanding. Herein, it is probed into the TRIM21-related literature and raises an interacting model about how TRIM21 orchestrates proteins in cytosol. In this novel model, TRIM21 generally interacts with miscellaneous protein in intracellular immunity in two ways: For one, TRIM21 solely plays as an E3, ubiquitylating a glut of proteins that contain specific interferon-regulatory factor, nuclear transcription factor kappaB, virus sensors and others, and involving inflammatory responses. For another, TRIM21 serves as both E3 and specific FcR that detects antibody-complexes and facilitates antibody destroying target proteins. Correspondingly delineated as Fc-independent signaling and Fc-dependent signaling in this review, how TRIM21's interactions contribute to intracellular immunity, expecting to provide a systematical understanding of this important protein and invest enlightenment for further research on the pathogenesis of related diseases and its prospective application is elaborated.

Dysregulated microRNAs in prostate cancer: In silico prediction and in vitro validation

Objectives

MicroRNAs, which are micro-coordinators of gene expression, have been recently investigated as a potential treatment for cancer. The study used computational techniques to identify microRNAs that could target a set of genes simultaneously. Due to their multi-target-directed nature, microRNAs have the potential to impact multiple key pathways and their pathogenic cross-talk.

Materials and Methods

We identified microRNAs that target a prostate cancer-associated gene set using integrated bioinformatics analyses and experimental validation. The candidate gene set included genes targeted by clinically approved prostate cancer medications. We used STRING, GO, and KEGG web tools to confirm gene-gene interactions and their clinical significance. Then, we employed integrated predicted and validated bioinformatics approaches to retrieve hsa-miR-124-3p, 16-5p, and 27a-3p as the top three relevant microRNAs. KEGG and DIANA-miRPath showed the related pathways for the candidate genes and microRNAs.

Results

The Real-time PCR results showed that miR-16-5p simultaneously down-regulated all genes significantly except for PIK3CA/CB in LNCaP; miR-27a-3p simultaneously down-regulated all genes significantly, excluding MET in LNCaP and PIK3CA in PC-3; and miR-124-3p could not down-regulate significantly PIK3CB, MET, and FGFR4 in LNCaP and FGFR4 in PC-3. Finally, we used a cell cycle assay to show significant G0/G1 arrest by transfecting miR-124-3p in LNCaP and miR-16-5p in both cell lines.

Conclusion

Our findings suggest that this novel approach may have therapeutic benefits and these predicted microRNAs could effectively target the candidate genes.

Computational analysis of protein-protein interactions of cancer drivers in renal cell carcinoma

Renal cell carcinoma (RCC) is the most common type of kidney cancer with rising cases in recent years. Extensive research has identified various cancer driver proteins associated with different subtypes of RCC. Most RCC drivers are encoded by tumor suppressor genes and exhibit enrichment in functional categories such as protein degradation, chromatin remodeling, and transcription. To further our understanding of RCC, we utilized powerful deep-learning methods based on AlphaFold to predict protein-protein interactions (PPIs) involving RCC drivers. We predicted high-confidence complexes formed by various RCC drivers, including TCEB1, KMT2C/D and KDM6A of the COMPASS-related complexes, TSC1 of the MTOR pathway, and TRRAP. These predictions provide valuable structural insights into the interaction interfaces, some of which are promising targets for cancer drug design, such as the NRF2-MAFK interface. Cancer somatic missense mutations from large datasets of genome sequencing of RCCs were mapped to the interfaces of predicted and experimental structures of PPIs involving RCC drivers, and their effects on the binding affinity were evaluated. We observed more than 100 cancer somatic mutations affecting the binding affinity of complexes formed by key RCC drivers such as VHL and TCEB1. These findings emphasize the importance of these mutations in RCC pathogenesis and potentially offer new avenues for targeted therapies.

Genome-wide identification of the OMT gene family in Cucumis melo L. and expression analysis under abiotic and biotic stress

Background

O-methyltransferase (OMT)-mediated O-methylation is a frequent modification that occurs during natural product biosynthesis, and it increases the diversity and stability of secondary metabolites. However, detailed genome-wide identification and expression analyses of OMT gene family members have not been performed in melons. In this study, we aimed to perform the genome-wide identification of OMT gene family members in melon to identify and clarify their actions during stress.

Methods

Genome-wide identification of OMT gene family members was performed using data from the melon genome database. The Cucumis melo OMT genes (CmOMTs) were then compared with the genes from two representative monocotyledons and three representative dicotyledons. The basic information, cis-regulatory elements in the promoter, predicted 3-D-structures, and GO enrichment results of the 21 CmOMTs were analyzed.

Results

In our study, 21 CmOMTs (named CmOMT1-21) were obtained by analyzing the melon genome. These genes were located on six chromosomes and divided into three groups composed of nine, six, and six CmOMTs based on phylogenetic analysis. Gene structure and motif descriptions were similar within the same classes. Each CmOMT gene contains at least one cis-acting element associated with hormone transport regulation. Analysis of cis-acting elements illustrated the potential role of CmOMTs in developmental regulation and adaptations to various abiotic and biotic stresses. The RNA-seq and quantitative real-time PCR (qRT-PCR) results indicated that NaCl stress significantly induced CmOMT6/9/14/18 and chilling and high temperature and humidity (HTH) stresses significantly upregulated CmOMT14/18. Furthermore, the expression pattern of CmOMT18 may be associated with Fusarium oxysporum f. sp. melonis race 1.2 (FOM1.2) and powdery mildew resistance. Our study tentatively explored the biological functions of CmOMT genes in various stress regulation pathways and provided a conceptual basis for further detailed studies of the molecular mechanisms.

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.

Childhood environment influences epigenetic age and methylation concordance of a CpG clock locus in British-Bangladeshi migrants

Migration from one location to another often comes with a change in environmental conditions. Here, we analysed features of DNA methylation in young, adult British-Bangladeshi women who experienced different environments during their childhoods: a) migrants, who grew up in Bangladesh with exposure to comparatively higher pathogen loads and poorer health care, and b) second-generation British-Bangladeshis, born to Bangladeshi parents, who grew up in the UK. We used buccal DNA to estimate DNA methylation-based age (DNAm age) from 14 migrants and 11 second-generation migrants, aged 18-35 years. 'AgeAccel,' a measure of DNAm age, independent of chronological age, showed that the group of women who spent their childhood in Bangladesh had higher AgeAccel (P = 0.028), compared to their UK peers. Since epigenetic clocks have been proposed to be associated with maintenance processes of epigenetic systems, we evaluated the preference for concordant DNA methylation at the luteinizing hormone/choriogonadotropin receptor (LHCGR/LHR) locus, which harbours one of the CpGs contributing to Horvath's epigenetic clock. Measurements on both strands of individual, double-stranded DNA molecules indicate higher stability of DNA methylation states at this LHCGR/LHR locus in samples of women who grew up in Bangladesh. Together, our two independent analytical approaches imply that childhood environments may induce subtle changes that are detectable long after exposure occurred, which might reflect altered activity of the epigenetic maintenance system or a difference in the proportion of cell types in buccal tissue. This exploratory work supports our earlier findings that adverse childhood environments lead to phenotypic life history trade-offs.