Defining the separation landscape of topological domains for decoding consensus domain organization of the 3D genome.

Topologically associating domains (TADs) have emerged as basic structural and functional units of genome organization and have been determined by many computational methods from Hi-C contact maps. However, the TADs obtained by different methods vary greatly, which makes the accurate determination of TADs a challenging issue and hinders subsequent biological analyses about their organization and functions. Obvious inconsistencies among the TADs identified by different methods indeed make the statistical and biological properties of TADs overly depend on the chosen method rather than on the data. To this end, we use the consensus structural information captured by these methods to define the TAD separation landscape for decoding the consensus domain organization of the 3D genome. We show that the TAD separation landscape could be used to compare domain boundaries across multiple cell types for discovering conserved and divergent topological structures, decipher three types of boundary regions with diverse biological features, and identify consensus TADs (ConsTADs). We illustrate that these analyses could deepen our understanding of the relationships between the topological domains and chromatin states, gene expression, and DNA replication timing.

PhenoDriver: interpretable framework for studying personalized phenotype-associated driver genes in breast cancer.

Identifying personalized cancer driver genes and further revealing their oncogenic mechanisms is critical for understanding the mechanisms of cell transformation and aiding clinical diagnosis. Almost all existing methods primarily focus on identifying driver genes at the cohort or individual level but fail to further uncover their underlying oncogenic mechanisms. To fill this gap, we present an interpretable framework, PhenoDriver, to identify personalized cancer driver genes, elucidate their roles in cancer development and uncover the association between driver genes and clinical phenotypic alterations. By analyzing 988 breast cancer patients, we demonstrate the outstanding performance of PhenoDriver in identifying breast cancer driver genes at the cohort level compared to other state-of-the-art methods. Otherwise, our PhenoDriver can also effectively identify driver genes with both recurrent and rare mutations in individual patients. We further explore and reveal the oncogenic mechanisms of some known and unknown breast cancer driver genes (e.g. TP53, MAP3K1, HTT, etc.) identified by PhenoDriver, and construct their subnetworks for regulating clinical abnormal phenotypes. Notably, most of our findings are consistent with existing biological knowledge. Based on the personalized driver profiles, we discover two existing and one unreported breast cancer subtypes and uncover their molecular mechanisms. These results intensify our understanding for breast cancer mechanisms, guide therapeutic decisions and assist in the development of targeted anticancer therapies.

A novel heterophilic graph diffusion convolutional network for identifying cancer driver genes.

Identifying cancer driver genes plays a curial role in the development of precision oncology and cancer therapeutics. Although a plethora of methods have been developed to tackle this problem, the complex cancer mechanisms and intricate interactions between genes still make the identification of cancer driver genes challenging. In this work, we propose a novel machine learning method of heterophilic graph diffusion convolutional networks (called HGDCs) to boost cancer-driver gene identification. Specifically, HGDC first introduces graph diffusion to generate an auxiliary network for capturing the structurally similar nodes in a biomolecular network. Then, HGDC designs an improved message aggregation and propagation scheme to adapt to the heterophilic setting of biomolecular networks, alleviating the problem of driver gene features being smoothed by its neighboring dissimilar genes. Finally, HGDC uses a layer-wise attention classifier to predict the probability of one gene being a cancer driver gene. In the comparison experiments with other existing state-of-the-art methods, our HGDC achieves outstanding performance in identifying cancer driver genes. The experimental results demonstrate that HGDC not only effectively identifies well-known driver genes on different networks but also novel candidate cancer genes. Moreover, HGDC can effectively prioritize cancer driver genes for individual patients. Particularly, HGDC can identify patient-specific additional driver genes, which work together with the well-known driver genes to cooperatively promote tumorigenesis.

A social theory-enhanced graph representation learning framework for multitask prediction of drug-drug interactions.

Current machine learning-based methods have achieved inspiring predictions in the scenarios of mono-type and multi-type drug-drug interactions (DDIs), but they all ignore enhancive and depressive pharmacological changes triggered by DDIs. In addition, these pharmacological changes are asymmetric since the roles of two drugs in an interaction are different. More importantly, these pharmacological changes imply significant topological patterns among DDIs. To address the above issues, we first leverage Balance theory and Status theory in social networks to reveal the topological patterns among directed pharmacological DDIs, which are modeled as a signed and directed network. Then, we design a novel graph representation learning model named SGRL-DDI (social theory-enhanced graph representation learning for DDI) to realize the multitask prediction of DDIs. SGRL-DDI model can capture the task-joint information by integrating relation graph convolutional networks with Balance and Status patterns. Moreover, we utilize task-specific deep neural networks to perform two tasks, including the prediction of enhancive/depressive DDIs and the prediction of directed DDIs. Based on DDI entries collected from DrugBank, the superiority of our model is demonstrated by the comparison with other state-of-the-art methods. Furthermore, the ablation study verifies that Balance and Status patterns help characterize directed pharmacological DDIs, and that the joint of two tasks provides better DDI representations than individual tasks. Last, we demonstrate the practical effectiveness of our model by a version-dependent test, where 88.47 and 81.38% DDI out of newly added entries provided by the latest release of DrugBank are validated in two predicting tasks respectively.

m6Aexpress-enet: Predicting the regulatory expression m6A sites by an enet-regularization negative binomial regression model.

As the most abundant mRNA modification, m6A controls and influences many aspects of mRNA metabolism including the mRNA stability and degradation. However, the role of specific m6A sites in regulating gene expression still remains unclear. In additional, the multicollinearity problem caused by the correlation of methylation level of multiple m6A sites in each gene could influence the prediction performance. To address the above challenges, we propose an elastic-net regularized negative binomial regression model (called m6Aexpress-enet) to predict which m6A site could potentially regulate its gene expression. Comprehensive evaluations on simulated datasets demonstrate that m6Aexpress-enet could achieve the top prediction performance. Applying m6Aexpress-enet on real MeRIP-seq data from human lymphoblastoid cell lines, we have uncovered the complex regulatory pattern of predicted m6A sites and their unique enrichment pathway of the constructed co-methylation modules. m6Aexpress-enet proves itself as a powerful tool to enable biologists to discover the mechanism of m6A regulatory gene expression. Furthermore, the source code and the step-by-step implementation of m6Aexpress-enet is freely accessed at https://github.com/tengzhangs/m6Aexpress-enet.

Tri-Step Water-Assisted Strategy for Suppressing Cs4PbBr6 Phase in Printable Carbon-Based CsPbBr3 Solar Cells to Achieve High Stability.

Often deemed the "natural nemesis" of perovskites, water molecules have been largely circumvented by the majority of researchers in the field of perovskite solar cells. This has resulted in significant hurdles in investigating the beneficial impacts of water molecules on perovskite crystallization. Herein, it is found that by utilizing ethanol with minimal water content and subjecting all-inorganic perovskite to three distinct annealing temperatures within the same solvent, the residual CsBr can be effectively removed, and the formation of the Cs4PbBr6 phase can be curtailed. By selecting an optimal water content, substantial improvements are observed in the crystalline quality of CsPbBr3, the perovskite/carbon interface, and the mesoporous filling effect. The Urbach energy (Eu) is reduced from 38.96 to 35.59 meV, and the defect density decreased from 4.16 × 1014 to 3.39 × 1014 cm-3. As a result, the power conversion efficiency (PCE) improved from 7.55% in the control group to 9.37%. Under severe environmental conditions with a temperature (T) of 85 °C and a relative humidity (RH) of 40%, tracking tests over 1200 h retained 89.3% of the initial PCE. This research signifies a breakthrough in the fabrication of highly stable and efficient all-inorganic printable mesoscopic perovskite solar cells.

Competitive Formation Mechanism for Bidentate Passivation of Halogen Vacancies in Perovskite Based on 6-Chloropurine.

For metal halide perovskite solar cells, bidentate passivation (BP) is highly effective, but currently, only passivation sites rather than molecular environments are being considered. Here, the authors report an effective approach for high-performance fully printable mesoscopic perovskite solar cells (FP-PSCs) through the BP strategy using the multidentate molecule 6-chloropurine (6-CP). By utilizing density functional theory (DFT) calculations, X-ray photoelectron spectroscopy (XPS), and Fourier transform infrared spectroscopy (FTIR) characterizations, the competition mechanism is identified of BP between the chlorine atom and neighboring nitrogen atom of the imidazole and pyrimidine rings. Through BP between the chlorine atom and adjacent nitrogen atom in imidazole, the power conversion efficiency (PCE) of the pristine samples is significantly enhanced from 16.25% to 17.63% with 6-CP. The formation of BP enhances interfacial hole selectivity and charge transfer, and suppresses nonradiative recombination, improving device stability under high humidity conditions. The competition mechanism of BP between two aromatic cycles provides a path for designing molecular passivants and selecting passivation pathways to approach theoretical limits.

m 6 Aexpress-BHM: predicting m6A regulation of gene expression in multiple-groups context by a Bayesian hierarchical mixture model.

As the most abundant RNA modification, N6-methyladenosine (m6A) plays an important role in various RNA activities including gene expression and translation. With the rapid application of MeRIP-seq technology, samples of multiple groups, such as the involved multiple viral/ bacterial infection or distinct cell differentiation stages, are extracted from same experimental unit. However, our current knowledge about how the dynamic m6A regulating gene expression and the role in certain biological processes (e.g. immune response in this complex context) is largely elusive due to lack of effective tools. To address this issue, we proposed a Bayesian hierarchical mixture model (called m6Aexpress-BHM) to predict m6A regulation of gene expression (m6A-reg-exp) in multiple groups of MeRIP-seq experiment with limited samples. Comprehensive evaluations of m6Aexpress-BHM on the simulated data demonstrate its high predicting precision and robustness. Applying m6Aexpress-BHM on three real-world datasets (i.e. Flaviviridae infection, infected time-points of bacteria and differentiation stages of dendritic cells), we predicted more m6A-reg-exp genes with positive regulatory mode that significantly participate in innate immune or adaptive immune pathways, revealing the underlying mechanism of the regulatory function of m6A during immune response. In addition, we also found that m6A may influence the expression of PD-1/PD-L1 via regulating its interacted genes. These results demonstrate the power of m6Aexpress-BHM, helping us understand the m6A regulatory function in immune system.

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.

Arginine Methylation of MDH1 by CARM1 Inhibits Glutamine Metabolism and Suppresses Pancreatic Cancer.

Distinctive from their normal counterparts, cancer cells exhibit unique metabolic dependencies on glutamine to fuel anabolic processes. Specifically, pancreatic ductal adenocarcinoma (PDAC) cells rely on an unconventional metabolic pathway catalyzed by aspartate aminotransferase, malate dehydrogenase 1 (MDH1), and malic enzyme 1 to rewire glutamine metabolism and support nicotinamide adenine dinucleotide phosphate (NADPH) production. Here, we report that methylation on arginine 248 (R248) negatively regulates MDH1. Protein arginine methyltransferase 4 (PRMT4/CARM1) methylates and inhibits MDH1 by disrupting its dimerization. Knockdown of MDH1 represses mitochondria respiration and inhibits glutamine metabolism, which sensitizes PDAC cells to oxidative stress and suppresses cell proliferation. Meanwhile, re-expression of wild-type MDH1, but not its methylation-mimetic mutant, protects cells from oxidative injury and restores cell growth and clonogenic activity. Importantly, MDH1 is hypomethylated at R248 in clinical PDAC samples. Our study reveals that arginine methylation of MDH1 by CARM1 regulates cellular redox homeostasis and suppresses glutamine metabolism of pancreatic cancer.

Resilience function uncovers the critical transitions in cancer initiation.

Considerable evidence suggests that during the progression of cancer initiation, the state transition from wellness to disease is not necessarily smooth but manifests switch-like nonlinear behaviors, preventing the cancer prediction and early interventional therapy for patients. Understanding the mechanism of such wellness-to-disease transitions is a fundamental and challenging task. Despite the advances in flux theory of nonequilibrium dynamics and 'critical slowing down'-based system resilience theory, a system-level approach still lacks to fully describe this state transition. Here, we present a novel framework (called bioRFR) to quantify such wellness-to-disease transition during cancer initiation through uncovering the biological system's resilience function from gene expression data. We used bioRFR to reconstruct the biologically and dynamically significant resilience functions for cancer initiation processes (e.g. BRCA, LUSC and LUAD). The resilience functions display the similar resilience pattern with hysteresis feature but different numbers of tipping points, which implies that once the cell become cancerous, it is very difficult or even impossible to reverse to the normal state. More importantly, bioRFR can measure the severe degree of cancer patients and identify the personalized key genes that are associated with the individual system's state transition from normal to tumor in resilience perspective, indicating that bioRFR can contribute to personalized medicine and targeted cancer therapy.

Prediction of the transcription factor binding sites with meta-learning.

With the accumulation of ChIP-seq data, convolution neural network (CNN)-based methods have been proposed for predicting transcription factor binding sites (TFBSs). However, biological experimental data are noisy, and are often treated as ground truth for both training and testing. Particularly, existing classification methods ignore the false positive and false negative which are caused by the error in the peak calling stage, and therefore, they can easily overfit to biased training data. It leads to inaccurate identification and inability to reveal the rules of governing protein-DNA binding. To address this issue, we proposed a meta learning-based CNN method (namely TFBS_MLCNN or MLCNN for short) for suppressing the influence of noisy labels data and accurately recognizing TFBSs from ChIP-seq data. Guided by a small amount of unbiased meta-data, MLCNN can adaptively learn an explicit weighting function from ChIP-seq data and update the parameter of classifier simultaneously. The weighting function overcomes the influence of biased training data on classifier by assigning a weight to each sample according to its training loss. The experimental results on 424 ChIP-seq datasets show that MLCNN not only outperforms other existing state-of-the-art CNN methods, but can also detect noisy samples which are given the small weights to suppress them. The suppression ability to the noisy samples can be revealed through the visualization of samples' weights. Several case studies demonstrate that MLCNN has superior performance to others.

m6Acancer-Net: Identification of m6A-mediated cancer driver genes from gene-site heterogeneous network.

N6-methyladenosine (m6A) is the most abundant eukaryotic modification internal mRNA, which plays the crucial roles in the occurrence and development of cancer. However, current knowledge about m6A-mediated functional circuit and key genes targeted by m6A methylation in cancer is mostly elusive. Thus, here we proposed a novel network-based approach (called m6Acancer-Net) to identify m6A-mediated driver genes and their associated network in specific type of cancer, such as acute myeloid leukemia. m6A-mediated cancer driver genes are defined as genes mediated by m6A methylation, significantly mutated, and functionally interacted in cancer. m6Acancer-Net identified the m6A-mediated cancer driver genes by combining gene functional interaction network with RNA methylation, gene expression and mutation information. A cancer-specific gene-site heterogeneous network was firstly constructed by connecting the m6A site co-methylation network with the functional interaction pruned gene co-expression network generated from large scale gene expression profile of specific cancer. Then, the functional m6A-mediated genes were identified by selecting the m6A regulators as seed genes to perform the random walk with restart algorithm on the gene-site heterogeneous network. Finally, m6A-mediated cancer driver gene subnetworks were constructed by performing the heat diffusion of mutation frequency for functional m6A-mediated genes in protein-protein interaction networks. The experimental results of m6Acancer-Net on the acute myeloid leukemia (AML) and glioblastoma multiforme (GBM) data from TCGA project show that the m6A-mediated caner driver genes identified by m6Acancer-Net are targeted by m6A regulators, and mediate significant cancer-related pathways. They play crucial roles in development and prognostic stratification of cancer. Moreover, 15 m6A-mediated cancer driver genes identified in AML are validated by literatures to mediate AML progress, and 14 m6A-mediated cancer driver genes identified in GBM are validated by literatures to participate in development of GBM. m6Acancer-Net is reliable to identify the functionally significant m6A-mediated driver genes in specific cancer, and it can effectively facilitate the understanding of regulatory and therapeutic mechanism of cancer driver genes in epitranscriptome layer.

m6Aexpress-Reader: Prediction of m6A regulated expression genes by integrating m6A sites and reader binding information in specific- context.

N6-methyladenosine (m6A) is the most abundant form of mRNA modification and plays an important role in regulating gene expression. However, the mechanisms of m6A regulated gene expression in cell or condition specific, are still poorly understood. Even though, some methods are able to predict m6A regulated expression (m6A-reg-exp) genes in specific context, they don't introduce the m6A reader binding information, while this information can help to predict m6A-reg-exp genes and more clearly to explain the mechanisms of m6A-mediated gene expression process. Thus, by integrating m6A sites and reader binding information, we proposed a novel method (called m6Aexpress-Reader) to predict m6A-reg-exp genes from limited MeRIP-seq data in specific context. m6Aexpress-Reader adopts the reader binding signal strength to weight the posterior distribution of the estimated regulatory coefficients for enhancing the prediction power. By using m6Aexpress-Reader, we found the complex characteristic of m6A on gene expression regulation and the distinct regulated pattern of m6A-reg-exp genes with different reader binding. m6A readers, YTHDF2 or IGF2BP1/3 all play an important role in various cancers and the key cancer pathways. In addition, m6Aexpress-Reader reveals the distinct m6A regulated mode of reader targeted genes in cancer. m6Aexpress-Reader could be a useful tool for studying the m6A regulation on reader target genes in specific context and it can be freely accessible at: https://github.com/NWPU-903PR/m6AexpressReader.

Network controllability-based algorithm to target personalized driver genes for discovering combinatorial drugs of individual patients.

Multiple driver genes in individual patient samples may cause resistance to individual drugs in precision medicine. However, current computational methods have not studied how to fill the gap between personalized driver gene identification and combinatorial drug discovery for individual patients. Here, we developed a novel structural network controllability-based personalized driver genes and combinatorial drug identification algorithm (CPGD), aiming to identify combinatorial drugs for an individual patient by targeting personalized driver genes from network controllability perspective. On two benchmark disease datasets (i.e. breast cancer and lung cancer datasets), performance of CPGD is superior to that of other state-of-the-art driver gene-focus methods in terms of discovery rate among prior-known clinical efficacious combinatorial drugs. Especially on breast cancer dataset, CPGD evaluated synergistic effect of pairwise drug combinations by measuring synergistic effect of their corresponding personalized driver gene modules, which are affected by a given targeting personalized driver gene set of drugs. The results showed that CPGD performs better than existing synergistic combinatorial strategies in identifying clinical efficacious paired combinatorial drugs. Furthermore, CPGD enhanced cancer subtyping by computationally providing personalized side effect signatures for individual patients. In addition, CPGD identified 90 drug combinations candidates from SARS-COV2 dataset as potential drug repurposing candidates for recently spreading COVID-19.

m6A-express: uncovering complex and condition-specific m6A regulation of gene expression.

N6-methyladenosine (m6A) is the most abundant form of mRNA modification and controls many aspects of RNA metabolism including gene expression. However, the mechanisms by which m6A regulates cell- and condition-specific gene expression are still poorly understood, partly due to a lack of tools capable of identifying m6A sites that regulate gene expression under different conditions. Here we develop m6A-express, the first algorithm for predicting condition-specific m6A regulation of gene expression (m6A-reg-exp) from limited methylated RNA immunoprecipitation sequencing (MeRIP-seq) data. Comprehensive evaluations of m6A-express using simulated and real data demonstrated its high prediction specificity and sensitivity. When only a few MeRIP-seq samples may be available for the cellular or treatment conditions, m6A-express is particularly more robust than the log-linear model. Using m6A-express, we reported that m6A writers, METTL3 and METTL14, competitively regulate the transcriptional processes by mediating m6A-reg-exp of different genes in Hela cells. In contrast, METTL3 induces different m6A-reg-exp of a distinct group of genes in HepG2 cells to regulate protein functions and stress-related processes. We further uncovered unique m6A-reg-exp patterns in human brain and intestine tissues, which are enriched in organ-specific processes. This study demonstrates the effectiveness of m6A-express in predicting condition-specific m6A-reg-exp and highlights the complex, condition-specific nature of m6A-regulation of gene expression.

Radiomics model based on intravoxel incoherent motion and diffusion kurtosis imaging for predicting histopathological grade and Ki-67 expression level of soft tissue sarcomas.

BACKGROUND: Accurate identification of the histopathological grade and the Ki-67 expression level is important in clinical cases of soft tissue sarcomas (STSs). PURPOSE: To explore the feasibility of a radiomics model based on intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI) and diffusion kurtosis imaging (DKI) MRI parameter maps in predicting the histopathological grade and Ki-67 expression level of STSs. MATERIAL AND METHODS: In total, 42 patients diagnosed with STSs between May 2018 and January 2020 were selected. The MADC software in Functool of GE ADW 4.7 workstation was used to obtain standard apparent diffusion coefficient (ADC), D, D*, f, mean diffusivity, and mean kurtosis (MK). The histopathological grade and Ki-67 expression level of STSs were identified. The radiomics features of IVIM and DKI parameter maps were used as the dataset. The area under the receiver operating characteristic curve (AUC) and F1-score were calculated. RESULTS: D-SVM achieved the best diagnostic performance for histopathological grade. The AUC in the validation cohort was 0.88 (sensitivity: 0.75 [low level] and 0.83 [high level]; specificity: 0.83 [low level] and 0.75 [high level]; F1-score: 0.75 [low level] and 0.83 [high level]). MK-SVM achieved the best diagnostic performance for Ki-67 expression level. The AUC in the validation cohort was 0.83 (sensitivity: 0.83 [low level] and 0.50 [high level; specificity: 0.50 [low level] and 0.83 [high level]; F1-score: 0.77 [low level] and 0.57 [high level]). CONCLUSION: The proposed radiomics classifier could predict the pathological grade of STSs and the Ki-67 expression level in STSs.

[Effect on Danggui Shaoyao Powder on mitophagy in rat model of Alzheimer's disease based on PINK1-Parkin pathway].

This study investigated the mechanism of Danggui Shaoyao Powder(DSP) against mitophagy in rat model of Alzheimer's disease(AD) induced by streptozotocin(STZ) based on PTEN induced putative kinase 1(PINK1)-Parkin signaling pathway. The AD rat model was established by injecting STZ into the lateral ventricle, and the rats were divided into normal group, model group, DSP low-dose group(12 g·kg~(-1)·d~(-1)), DSP medium-dose group(24 g·kg~(-1)·d~(-1)), and DSP high-dose group(36 g·kg~(-1)·d~(-1)). Morris water maze test was used to detect the learning and memory function of the rats, and transmission electron microscopy and immunofluorescence were employed to detect mitophagy. The protein expression levels of PINK1, Parkin, LC3BⅠ/LC3BⅡ, and p62 were assayed by Western blot. Compared with the normal group, the model group showed a significant decrease in the learning and memory function(P<0.01), reduced protein expression of PINK1 and Parkin(P<0.05), increased protein expression of LC3BⅠ/LC3BⅡ and p62(P<0.05), and decreased occurrence of mitophagy(P<0.01). Compared with the model group, the DSP medium-and high-dose groups notably improved the learning and memory ability of AD rats, which mainly manifested as shortened escape latency, leng-thened time in target quadrants and elevated number of crossing the platform(P<0.05 or P<0.01), remarkably activated mitophagy(P<0.05), up-regulated the protein expression of PINK1 and Parkin, and down-regulated the protein expression of LC3BⅠ/LC3BⅡ and p62(P<0.05 or P<0.01). These results demonstrated that DSP might promote mitophagy mediated by PINK1-Parkin pathway to remove damaged mitochondria and improve mitochondrial function, thereby exerting a neuroprotective effect.

Prioritization of cancer driver gene with prize-collecting steiner tree by introducing an edge weighted strategy in the personalized gene interaction network.

BACKGROUND: Cancer is a heterogeneous disease in which tumor genes cooperate as well as adapt and evolve to the changing conditions for individual patients. It is a meaningful task to discover the personalized cancer driver genes that can provide diagnosis and target drug for individual patients. However, most of existing methods mainly ranks potential personalized cancer driver genes by considering the patient-specific nodes information on the gene/protein interaction network. These methods ignore the personalized edge weight information in gene interaction network, leading to false positive results. RESULTS: In this work, we presented a novel algorithm (called PDGPCS) to predict the Personalized cancer Driver Genes based on the Prize-Collecting Steiner tree model by considering the personalized edge weight information. PDGPCS first constructs the personalized weighted gene interaction network by integrating the personalized gene expression data and prior known gene/protein interaction network knowledge. Then the gene mutation data and pathway data are integrated to quantify the impact of each mutant gene on every dysregulated pathway with the prize-collecting Steiner tree model. Finally, according to the mutant gene's aggregated impact score on all dysregulated pathways, the mutant genes are ranked for prioritizing the personalized cancer driver genes. Experimental results on four TCGA cancer datasets show that PDGPCS has better performance than other personalized driver gene prediction methods. In addition, we verified that the personalized edge weight of gene interaction network can improve the prediction performance. CONCLUSIONS: PDGPCS can more accurately identify the personalized driver genes and takes a step further toward personalized medicine and treatment. The source code of PDGPCS can be freely downloaded from https://github.com/NWPU-903PR/PDGPCS .

Network control principles for identifying personalized driver genes in cancer.

To understand tumor heterogeneity in cancer, personalized driver genes (PDGs) need to be identified for unraveling the genotype-phenotype associations corresponding to particular patients. However, most of the existing driver-focus methods mainly pay attention on the cohort information rather than on individual information. Recent developing computational approaches based on network control principles are opening a new way to discover driver genes in cancer, particularly at an individual level. To provide comprehensive perspectives of network control methods on this timely topic, we first considered the cancer progression as a network control problem, in which the expected PDGs are altered genes by oncogene activation signals that can change the individual molecular network from one health state to the other disease state. Then, we reviewed the network reconstruction methods on single samples and introduced novel network control methods on single-sample networks to identify PDGs in cancer. Particularly, we gave a performance assessment of the network structure control-based PDGs identification methods on multiple cancer datasets from TCGA, for which the data and evaluation package also are publicly available. Finally, we discussed future directions for the application of network control methods to identify PDGs in cancer and diverse biological processes.