expHRD: an individualized, transcriptome-based prediction model for homologous recombination deficiency assessment in cancer.

BACKGROUND: Homologous recombination deficiency (HRD) stands as a clinical indicator for discerning responsive outcomes to platinum-based chemotherapy and poly ADP-ribose polymerase (PARP) inhibitors. One of the conventional approaches to HRD prognostication has generally centered on identifying deleterious mutations within the BRCA1/2 genes, along with quantifying the genomic scars, such as Genomic Instability Score (GIS) estimation with scarHRD. However, the scarHRD method has limitations in scenarios involving tumors bereft of corresponding germline data. Although several RNA-seq-based HRD prediction algorithms have been developed, they mainly support cohort-wise classification, thereby yielding HRD status without furnishing an analogous quantitative metric akin to scarHRD. This study introduces the expHRD method, which operates as a novel transcriptome-based framework tailored to n-of-1-style HRD scoring. RESULTS: The prediction model has been established using the elastic net regression method in the Cancer Genome Atlas (TCGA) pan-cancer training set. The bootstrap technique derived the HRD geneset for applying the expHRD calculation. The expHRD demonstrated a notable correlation with scarHRD and superior performance in predicting HRD-high samples. We also performed intra- and extra-cohort evaluations for clinical feasibility in the TCGA-OV and the Genomic Data Commons (GDC) ovarian cancer cohort, respectively. The innovative web service designed for ease of use is poised to extend the realms of HRD prediction across diverse malignancies, with ovarian cancer standing as an emblematic example. CONCLUSIONS: Our novel approach leverages the transcriptome data, enabling the prediction of HRD status with remarkable precision. This innovative method addresses the challenges associated with limited available data, opening new avenues for utilizing transcriptomics to inform clinical decisions.

Cell-cell communication network-based interpretable machine learning predicts cancer patient response to immune checkpoint inhibitors.

Immune checkpoint inhibitors (ICIs) have revolutionized cancer treatment. However, only some patients respond to ICIs, and current biomarkers for ICI efficacy have limited performance. Here, we devised an interpretable machine learning (ML) model trained using patient-specific cell-cell communication networks (CCNs) decoded from the patient's bulk tumor transcriptome. The model could (i) predict ICI efficacy for patients across four cancer types (median AUROC: 0.79) and (ii) identify key communication pathways with crucial players responsible for patient response or resistance to ICIs by analyzing more than 700 ICI-treated patient samples from 11 cohorts. The model prioritized chemotaxis communication of immune-related cells and growth factor communication of structural cells as the key biological processes underlying response and resistance to ICIs, respectively. We confirmed the key communication pathways and players at the single-cell level in patients with melanoma. Our network-based ML approach can be used to expand ICIs' clinical benefits in cancer patients.

Drug approval prediction based on the discrepancy in gene perturbation effects between cells and humans.

BACKGROUND: Poor translation between in vitro and clinical studies due to the cells/humans discrepancy in drug target perturbation effects leads to safety failures in clinical trials, thus increasing drug development costs and reducing patients' life quality. Therefore, developing a predictive model for drug approval considering the cells/humans discrepancy is needed to reduce drug attrition rates in clinical trials. METHODS: Our machine learning framework predicts drug approval in clinical trials based on the cells/humans discrepancy in drug target perturbation effects. To evaluate the discrepancy to predict drug approval (1404 approved and 1070 unapproved drugs), we analysed CRISPR-Cas9 knockout and loss-of-function mutation rate-based gene perturbation effects on cells and humans, respectively. To validate the risk of drug targets with the cells/humans discrepancy, we examined the targets of failed and withdrawn drugs due to safety problems. FINDINGS: Drug approvals in clinical trials were correlated with the cells/humans discrepancy in gene perturbation effects. Genes tolerant to perturbation effects on cells but intolerant to those on humans were associated with failed drug targets. Furthermore, genes with the cells/humans discrepancy were related to drugs withdrawn due to severe side effects. Motivated by previous studies assessing drug safety through chemical properties, we improved drug approval prediction by integrating chemical information with the cells/humans discrepancy. INTERPRETATION: The cells/humans discrepancy in gene perturbation effects facilitates drug approval prediction and explains drug safety failures in clinical trials. FUNDING: S.K. received grants from the Korean National Research Foundation (2021R1A2B5B01001903 and 2020R1A6A1A03047902) and IITP (2019-0-01906, Artificial Intelligence Graduate School Program, POSTECH).

Information about immune cell proportions and tumor stage improves the prediction of recurrence in patients with colorectal cancer.

Predicting cancer recurrence is essential to improving the clinical outcomes of patients with colorectal cancer (CRC). Although tumor stage information has been used as a guideline to predict CRC recurrence, patients with the same stage show different clinical outcomes. Therefore, there is a need to develop a method to identify additional features for CRC recurrence prediction. Here, we developed a network-integrated multiomics (NIMO) approach to select appropriate transcriptome signatures for better CRC recurrence prediction by comparing the methylation signatures of immune cells. We validated the performance of the CRC recurrence prediction based on two independent retrospective cohorts of 114 and 110 patients. Moreover, to confirm that the prediction was improved, we used both NIMO-based immune cell proportions and TNM (tumor, node, metastasis) stage data. This work demonstrates the importance of (1) using both immune cell composition and TNM stage data and (2) identifying robust immune cell marker genes to improve CRC recurrence prediction.

An evolution-based machine learning to identify cancer type-specific driver mutations.

Identifying cancer type-specific driver mutations is crucial for illuminating distinct pathologic mechanisms across various tumors and providing opportunities of patient-specific treatment. However, although many computational methods were developed to predict driver mutations in a type-specific manner, the methods still have room to improve. Here, we devise a novel feature based on sequence co-evolution analysis to identify cancer type-specific driver mutations and construct a machine learning (ML) model with state-of-the-art performance. Specifically, relying on 28 000 tumor samples across 66 cancer types, our ML framework outperformed current leading methods of detecting cancer driver mutations. Interestingly, the cancer mutations identified by sequence co-evolution feature are frequently observed in interfaces mediating tissue-specific protein-protein interactions that are known to associate with shaping tissue-specific oncogenesis. Moreover, we provide pre-calculated potential oncogenicity on available human proteins with prediction scores of all possible residue alterations through user-friendly website (http://sbi.postech.ac.kr/w/cancerCE). This work will facilitate the identification of cancer type-specific driver mutations in newly sequenced tumor samples.

Development of bioinformatics and multi-omics analyses in organoids.

Pre-clinical models are critical in gaining mechanistic and biological insights into disease progression. Recently, patient-derived organoid models have been developed to facilitate our understanding of disease development and to improve the discovery of therapeutic options by faithfully recapitulating in vivo tissues or organs. As technological developments of organoid models are rapidly growing, computational methods are gaining attention in organoid researchers to improve the ability to systematically analyze experimental results. In this review, we summarize the recent advances in organoid models to recapitulate human diseases and computational advancements to analyze experimental results from organoids. [BMB Reports 2023; 56(1): 43-48].

Deconvolution of bulk tumors into distinct immune cell states predicts colorectal cancer recurrence.

Predicting colorectal cancer recurrence after tumor resection is crucial because it promotes the administration of proper subsequent treatment or management to improve the clinical outcomes of patients. Several clinical or molecular factors, including tumor stage, metastasis, and microsatellite instability status, have been used to assess the risk of recurrence, although their predictive ability is limited. Here, we predicted colorectal cancer recurrence based on cellular deconvolution of bulk tumors into two distinct immune cell states: cancer-associated (tumor-infiltrating immune cell-like) and noncancer-associated (peripheral blood mononuclear cell-like). Prediction model performed significantly better when immune cells were deconvoluted into two states rather than a single state, suggesting that the difference in cancer recurrence was better explained by distinct states of immune cells. It indicates the importance of distinguishing immune cell states using cellular deconvolution to improve the prediction of colorectal cancer recurrence.

Enzyme activity engineering based on sequence co-evolution analysis.

The utility of engineering enzyme activity is expanding with the development of biotechnology. Conventional methods have limited applicability as they require high-throughput screening or three-dimensional structures to direct target residues of activity control. An alternative method uses sequence evolution of natural selection. A repertoire of mutations was selected for fine-tuning enzyme activities to adapt to varying environments during the evolution. Here, we devised a strategy called sequence co-evolutionary analysis to control the efficiency of enzyme reactions (SCANEER), which scans the evolution of protein sequences and direct mutation strategy to improve enzyme activity. We hypothesized that amino acid pairs for various enzyme activity were encoded in the evolutionary history of protein sequences, whereas loss-of-function mutations were avoided since those are depleted during the evolution. SCANEER successfully predicted the enzyme activities of beta-lactamase and aminoglycoside 3'-phosphotransferase. SCANEER was further experimentally validated to control the activities of three different enzymes of great interest in chemical production: cis-aconitate decarboxylase, α-ketoglutaric semialdehyde dehydrogenase, and inositol oxygenase. Activity-enhancing mutations that improve substrate-binding affinity or turnover rate were found at sites distal from known active sites or ligand-binding pockets. We provide SCANEER to control desired enzyme activity through a user-friendly webserver.

The implication of holocytochrome c synthase mutation in Korean familial hypoplastic amelogenesis imperfecta.

OBJECTIVES: This study aimed to comprehensively characterise genetic variants of amelogenesis imperfecta in a single Korean family through whole-exome sequencing and bioinformatics analysis. MATERIAL AND METHODS: Thirty-one individuals of a Korean family, 9 of whom were affected and 22 unaffected by amelogenesis imperfecta, were enrolled. Whole-exome sequencing was performed on 12 saliva samples, including samples from 8 affected and 4 unaffected individuals. The possible candidate genes associated with the disease were screened by segregation analysis and variant filtering. In silico mutation impact analysis was then performed on the filtered variants based on sequence conservation and protein structure. RESULTS: Whole-exome sequencing data revealed an X-linked dominant, heterozygous genomic missense mutation in the mitochondrial gene holocytochrome c synthase (HCCS). We also found that HCCS is potentially related to the role of mitochondria in amelogenesis. The HCCS variant was expected to be deleterious in both evolution-based and large population-based analyses. Further, the variant was predicted to have a negative effect on catalytic function of HCCS by in silico analysis of protein structure. In addition, HCCS had significant association with amelogenesis in literature mining analysis. CONCLUSIONS: These findings suggest new evidence for the relationship between amelogenesis and mitochondria function, which could be implicated in the pathogenesis of amelogenesis imperfecta. CLINICAL RELEVANCE: The discovery of HCCS mutations and a deeper understanding of the pathogenesis of amelogenesis imperfecta could lead to finding solutions for the fundamental treatment of this disease. Furthermore, it enables dental practitioners to establish predictable prosthetic treatment plans at an early stage by early detection of amelogenesis imperfecta through personalised medicine.

Evolutionary rewiring of regulatory networks contributes to phenotypic differences between human and mouse orthologous genes.

Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).

Two-Level Blockchain System for Digital Crime Evidence Management.

Digital evidence, such as evidence from CCTV and event data recorders, is highly valuable in criminal investigations, and is used as definitive evidence in trials. However, there are risks when digital evidence obtained during the investigation of a case is managed through a physical hard disk drive until it is submitted to the court. Previous studies have focused on the integrated management of digital evidence in a centralized system, but if a centralized system server is attacked, major operations and investigation information may be leaked. Therefore, there is a need to reliably manage digital evidence and investigation information using blockchain technology in a distributed system environment. However, when large amounts of data-such as evidence videos-are stored in a blockchain, the data that must be processed only within one block before being created increase, causing performance degradation. Therefore, we propose a two-level blockchain system that separates digital evidence into hot and cold blockchains. In the criminal investigation process, information that frequently changes is stored in the hot blockchain, and unchanging data such as videos are stored in the cold blockchain. To evaluate the system, we measured the storage and inquiry processing performance of digital crime evidence videos according to the different capacities in the two-level blockchain system.

Network-based machine learning in colorectal and bladder organoid models predicts anti-cancer drug efficacy in patients.

Cancer patient classification using predictive biomarkers for anti-cancer drug responses is essential for improving therapeutic outcomes. However, current machine-learning-based predictions of drug response often fail to identify robust translational biomarkers from preclinical models. Here, we present a machine-learning framework to identify robust drug biomarkers by taking advantage of network-based analyses using pharmacogenomic data derived from three-dimensional organoid culture models. The biomarkers identified by our approach accurately predict the drug responses of 114 colorectal cancer patients treated with 5-fluorouracil and 77 bladder cancer patients treated with cisplatin. We further confirm our biomarkers using external transcriptomic datasets of drug-sensitive and -resistant isogenic cancer cell lines. Finally, concordance analysis between the transcriptomic biomarkers and independent somatic mutation-based biomarkers further validate our method. This work presents a method to predict cancer patient drug responses using pharmacogenomic data derived from organoid models by combining the application of gene modules and network-based approaches.

Capicua restricts cancer stem cell-like properties in breast cancer cells.

Cancer stem cells (CSCs) play a central role in cancer initiation, progression, therapeutic resistance, and recurrence in patients. Here we present Capicua (CIC), a developmental transcriptional repressor, as a suppressor of CSC properties in breast cancer cells. CIC deficiency critically enhances CSC self-renewal and multiple CSC subpopulations of breast cancer cells without altering their growth rate or invasiveness. Loss of CIC relieves repression of ETV4 and ETV5 expression, consequently promoting self-renewal capability, EpCAM+/CD44+/CD24low/- expression, and ALDH activity. In xenograft models, CIC deficiency significantly increases CSC frequency and drives tumor initiation through derepression of ETV4. Consistent with the experimental data, the CD44high/CD24low CSC-like feature is inversely correlated with CIC levels in breast cancer patients. We also identify SOX2 as a downstream target gene of CIC that partly promotes CSC properties. Taken together, our study demonstrates that CIC suppresses breast cancer formation via restricting cancer stemness and proposes CIC as a potential regulator of stem cell maintenance.

Capicua suppresses colorectal cancer progression via repression of ETV4 expression.

BACKGROUND: Although major driver gene mutations have been identified, the complex molecular heterogeneity of colorectal cancer (CRC) remains unclear. Capicua (CIC) functions as a tumor suppressor in various types of cancers; however, its role in CRC progression has not been examined. METHODS: Databases for gene expression profile in CRC patient samples were used to evaluate the association of the levels of CIC and Polyoma enhancer activator 3 (PEA3) group genes (ETS translocation variant 1 (ETV1), ETV4, and ETV5), the best-characterized CIC targets in terms of CIC functions, with clinicopathological features of CRC. CIC and ETV4 protein levels were also examined in CRC patient tissue samples. Gain- and loss-of function experiments in cell lines and mouse xenograft models were performed to investigate regulatory functions of CIC and ETV4 in CRC cell growth and invasion. qRT-PCR and western blot analyses were performed to verify the CIC regulation of ETV4 expression in CRC cells. Rescue experiments were conducted using siRNA against ETV4 and CIC-deficient CRC cell lines. RESULTS: CIC expression was decreased in the tissue samples of CRC patients. Cell invasion, migration, and proliferation were enhanced in CIC-deficient CRC cells and suppressed in CIC-overexpressing cells. Among PEA3 group genes, ETV4 levels were most dramatically upregulated and inversely correlated with the CIC levels in CRC patient samples. Furthermore, derepression of ETV4 was more prominent in CIC-deficient CRC cells, when compared with that observed for ETV1 and ETV5. The enhanced cell proliferative and invasive capabilities in CIC-deficient CRC cells were completely recovered by knockdown of ETV4. CONCLUSION: Collectively, the CIC-ETV4 axis is not only a key module that controls CRC progression but also a novel therapeutic and/or diagnostic target for CRC.

Domain-mediated interactions for protein subfamily identification.

Within a protein family, proteins with the same domain often exhibit different cellular functions, despite the shared evolutionary history and molecular function of the domain. We hypothesized that domain-mediated interactions (DMIs) may categorize a protein family into subfamilies because the diversified functions of a single domain often depend on interacting partners of domains. Here we systematically identified DMI subfamilies, in which proteins share domains with DMI partners, as well as with various functional and physical interaction networks in individual species. In humans, DMI subfamily members are associated with similar diseases, including cancers, and are frequently co-associated with the same diseases. DMI information relates to the functional and evolutionary subdivisions of human kinases. In yeast, DMI subfamilies contain proteins with similar phenotypic outcomes from specific chemical treatments. Therefore, the systematic investigation here provides insights into the diverse functions of subfamilies derived from a protein family with a link-centric approach and suggests a useful resource for annotating the functions and phenotypic outcomes of proteins.

Link clustering explains non-central and contextually essential genes in protein interaction networks.

Recent studies have shown that many essential genes (EGs) change their essentiality across various contexts. Finding contextual EGs in pathogenic conditions may facilitate the identification of therapeutic targets. We propose link clustering as an indicator of contextual EGs that are non-central in protein-protein interaction (PPI) networks. In various human and yeast PPI networks, we found that 29-47% of EGs were better characterized by link clustering than by centrality. Importantly, non-central EGs were prone to change their essentiality across different human cell lines and between species. Compared with central EGs and non-EGs, non-central EGs had intermediate levels of expression and evolutionary conservation. In addition, non-central EGs exhibited a significant impact on communities at lower hierarchical levels, suggesting that link clustering is associated with contextual essentiality, as it depicts locally important nodes in network structures.

Evolutionary coupling analysis identifies the impact of disease-associated variants at less-conserved sites.

Genome-wide association studies have discovered a large number of genetic variants in human patients with the disease. Thus, predicting the impact of these variants is important for sorting disease-associated variants (DVs) from neutral variants. Current methods to predict the mutational impacts depend on evolutionary conservation at the mutation site, which is determined using homologous sequences and based on the assumption that variants at well-conserved sites have high impacts. However, many DVs at less-conserved but functionally important sites cannot be predicted by the current methods. Here, we present a method to find DVs at less-conserved sites by predicting the mutational impacts using evolutionary coupling analysis. Functionally important and evolutionarily coupled sites often have compensatory variants on cooperative sites to avoid loss of function. We found that our method identified known intolerant variants in a diverse group of proteins. Furthermore, at less-conserved sites, we identified DVs that were not identified using conservation-based methods. These newly identified DVs were frequently found at protein interaction interfaces, where species-specific mutations often alter interaction specificity. This work presents a means to identify less-conserved DVs and provides insight into the relationship between evolutionarily coupled sites and human DVs.

Homogeneous photocatalytic Fe3+/Fe2+ redox cycle for simultaneous Cr(VI) reduction and organic pollutant oxidation: Roles of hydroxyl radical and degradation intermediates.

The sustained oxidation of aqueous organic pollutants using hydroxyl radicals (HO) generated in the UV-irradiated solution of ferric ions was investigated in the presence of Cr(VI). The synergistic effect of simultaneous 4-chlorophenol (4-CP) oxidation and Cr(VI) reduction is explained in terms of the various roles of OH radical, degradation intermediates, and Fe3+/Fe2+ redox cycle. The photolysis of FeIII(OH)2+ generates OH radical which degrades the organic substrate. The reduction of Cr(VI) was inhibited by the OH radical-induced re-oxidation of Cr(III) in the absence of 4-CP. The complete removal of Cr(VI) was achieved only in the presence of phenolic substrates which not only reacts with OH radical (hence inhibiting the reoxidation of Cr(III)) but also generates reducing intermediates which effectively reduce Cr(VI). Fe2+ also converted Cr(VI) to Cr(III) with regenerating Fe3+, which makes the overall process photocatalytic. The photocatalytic activity for the simultaneous removal of 4-CP and Cr(VI) was largely maintained up to five cycles. Such simultaneous and synergic photoactivity was also observed for other phenolic compounds (4-bromophenol, 4-nitrophenol, phenol). The simultaneous and synergic removal of phenolic compounds and Cr(VI) can be enabled through the redox couple of Fe3+/Fe2+ working as a homogeneous photocatalyst.

Divergence of Noncoding Regulatory Elements Explains Gene-Phenotype Differences between Human and Mouse Orthologous Genes.

Mice have been widely used as a model organism to investigate human gene-phenotype relationships based on a conjecture that orthologous genes generally perform similar functions and are associated with similar phenotypes. However, phenotypes associated with orthologous genes often turn out to be quite different between human and mouse. Herein, we devised a method to quantitatively compare phenotypes annotations associated with mouse models and human. Using semantic similarity comparisons, we identified orthologous genes with different phenotype annotations, of which the similarity score is on a par with that of random gene pairs. Analysis of sequence evolution and transcriptomic changes revealed that orthologous genes with phenotypic differences are correlated with changes in noncoding regulatory elements and tissue-specific expression profiles rather than changes in protein-coding sequences. To map accurate gene-phenotype relationships using model organisms, we propose that careful consideration of the evolutionary divergence of noncoding regulatory elements and transcriptomic profiles is essential.

Capicua suppresses hepatocellular carcinoma progression by controlling the ETV4-MMP1 axis.

Hepatocellular carcinoma (HCC) is developed by multiple steps accompanying progressive alterations of gene expression, which leads to increased cell proliferation and malignancy. Although environmental factors and intracellular signaling pathways that are critical for HCC progression have been identified, gene expression changes and the related genetic factors contributing to HCC pathogenesis are still insufficiently understood. In this study, we identify a transcriptional repressor, Capicua (CIC), as a suppressor of HCC progression and a potential therapeutic target. Expression of CIC is posttranscriptionally reduced in HCC cells. CIC levels are correlated with survival rates in patients with HCC. CIC overexpression suppresses HCC cell proliferation and invasion, whereas loss of CIC exerts opposite effects in vivo as well as in vitro. Levels of polyoma enhancer activator 3 (PEA3) group genes, the best-known CIC target genes, are correlated with lethality in patients with HCC. Among the PEA3 group genes, ETS translocation variant 4 (ETV4) is the most significantly up-regulated in CIC-deficient HCC cells, consequently promoting HCC progression. Furthermore, it induces expression of matrix metalloproteinase 1 (MMP1), the MMP gene highly relevant to HCC progression, in HCC cells; and knockdown of MMP1 completely blocks the CIC deficiency-induced HCC cell proliferation and invasion. CONCLUSION: Our study demonstrates that the CIC-ETV4-MMP1 axis is a regulatory module controlling HCC progression. (Hepatology 2018;67:2287-2301).