Analysis of super-enhancer using machine learning and its application to medical biology.

The analysis of super-enhancers (SEs) has recently attracted attention in elucidating the molecular mechanisms of cancer and other diseases. SEs are genomic structures that strongly induce gene expression and have been reported to contribute to the overexpression of oncogenes. Because the analysis of SEs and integrated analysis with other data are performed using large amounts of genome-wide data, artificial intelligence technology, with machine learning at its core, has recently begun to be utilized. In promoting precision medicine, it is important to consider information from SEs in addition to genomic data; therefore, machine learning technology is expected to be introduced appropriately in terms of building a robust analysis platform with a high generalization performance. In this review, we explain the history and principles of SE, and the results of SE analysis using state-of-the-art machine learning and integrated analysis with other data are presented to provide a comprehensive understanding of the current status of SE analysis in the field of medical biology. Additionally, we compared the accuracy between existing machine learning methods on the benchmark dataset and attempted to explore the kind of data preprocessing and integration work needed to make the existing algorithms work on the benchmark dataset. Furthermore, we discuss the issues and future directions of current SE analysis.

Mechanism of ERBB2 gene overexpression by the formation of super-enhancer with genomic structural abnormalities in lung adenocarcinoma without clinically actionable genetic alterations.

BACKGROUND: In an extensive genomic analysis of lung adenocarcinomas (LUADs), driver mutations have been recognized as potential targets for molecular therapy. However, there remain cases where target genes are not identified. Super-enhancers and structural variants are frequently identified in several hundred loci per case. Despite this, most cancer research has approached the analysis of these data sets separately, without merging and comparing the data, and there are no examples of integrated analysis in LUAD. METHODS: We performed an integrated analysis of super-enhancers and structural variants in a cohort of 174 LUAD cases that lacked clinically actionable genetic alterations. To achieve this, we conducted both WGS and H3K27Ac ChIP-seq analyses using samples with driver gene mutations and those without, allowing for a comprehensive investigation of the potential roles of super-enhancer in LUAD cases. RESULTS: We demonstrate that most genes situated in these overlapped regions were associated with known and previously unknown driver genes and aberrant expression resulting from the formation of super-enhancers accompanied by genomic structural abnormalities. Hi-C and long-read sequencing data further corroborated this insight. When we employed CRISPR-Cas9 to induce structural abnormalities that mimicked cases with outlier ERBB2 gene expression, we observed an elevation in ERBB2 expression. These abnormalities are associated with a higher risk of recurrence after surgery, irrespective of the presence or absence of driver mutations. CONCLUSIONS: Our findings suggest that aberrant gene expression linked to structural polymorphisms can significantly impact personalized cancer treatment by facilitating the identification of driver mutations and prognostic factors, contributing to a more comprehensive understanding of LUAD pathogenesis.

Histone methyltransferase SUV420H1/KMT5B contributes to poor prognosis in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has a high rate of recurrence and poor prognosis, even after curative surgery. Multikinase inhibitors have been applied for HCC patients, but their effect has been restricted. This study aims to clarify the clinical impact of SUV420H1/KMT5B, one of the methyltransferases for histone H4 at lysine 20, and elucidate the novel mechanisms of HCC progression. We retrospectively investigated SUV420H1 expression using HCC clinical tissue samples employing immunohistochemical analysis (n = 350). We then performed loss-of-function analysis of SUV420H1 with cell cycle analysis, migration assay, invasion assay and RNA sequence for Gene Ontology (GO) pathway analysis in vitro, and animal experiments with xenograft mice in vivo. The SUV420H1-high-score group (n = 154) had significantly poorer prognosis for both 5-year overall and 2-year/5-year disease-free survival than the SUV420H1-low-score group (n = 196) (p < 0.001 and p < 0.05, respectively). The SUV420H1-high-score group had pathologically larger tumor size, more tumors, poorer differentiation, and more positive vascular invasion than the SUV420H1-low-score group. Multivariate analysis demonstrated that SUV420H1 high score was the poorest independent factor for overall survival. SUV420H1 knockdown could suppress cell cycle from G1 to S phase and cell invasion. GO pathway analysis showed that SUV420H1 contributed to cell proliferation, cell invasion, and/or metastasis. Overexpression of SUV420H1 clinically contributed to poor prognosis in HCC, and the inhibition of SUV420H1 could repress tumor progression and invasion both in vitro and in vivo; thus, further analyses of SUV420H1 are necessary for the discovery of future molecularly targeted drugs.

Clinical features and impact of p53 status on sporadic mismatch repair deficiency and Lynch syndrome in uterine cancer.

The clinical features of sporadic mismatch repair deficiency (MMRd) and Lynch syndrome (LS) in Japanese patients with endometrial cancer (EC) were examined by evaluating the prevalence and prognostic factors of LS and sporadic MMRd in patients with EC. Targeted sequencing of five LS susceptibility genes (MLH1, MSH2, MSH6, PMS2, and EPCAM) was carried out in 443 patients with EC who were pathologically diagnosed with EC at the National Cancer Center Hospital between 2011 and 2018. Pathogenic variants in these genes were detected in 16 patients (3.7%). Immunohistochemistry for MMR proteins was undertaken in 337 of the 433 (77.9%) EC patients, and 91 patients (27.0%) showed absent expression of at least one MMR protein. The 13 cases of LS with MMR protein loss (93.8%) showed a favorable prognosis with a 5-year overall survival (OS) rate of 100%, although there was no statistically significant difference between this group and the sporadic MMRd group (p = 0.27). In the MMRd without LS group, the 5-year OS rate was significantly worse in seven patients with an aberrant p53 expression pattern than in those with p53 WT (53.6% vs. 93.9%, log-rank test; p = 0.0016). These results suggest that p53 abnormalities and pathogenic germline variants in MMR genes could be potential biomarkers for the molecular classification of EC with MMRd.

Application of non-negative matrix factorization in oncology: one approach for establishing precision medicine.

The increase in the expectations of artificial intelligence (AI) technology has led to machine learning technology being actively used in the medical field. Non-negative matrix factorization (NMF) is a machine learning technique used for image analysis, speech recognition, and language processing; recently, it is being applied to medical research. Precision medicine, wherein important information is extracted from large-scale medical data to provide optimal medical care for every individual, is considered important in medical policies globally, and the application of machine learning techniques to this end is being handled in several ways. NMF is also introduced differently because of the characteristics of its algorithms. In this review, the importance of NMF in the field of medicine, with a focus on the field of oncology, is described by explaining the mathematical science of NMF and the characteristics of the algorithm, providing examples of how NMF can be used to establish precision medicine, and presenting the challenges of NMF. Finally, the direction regarding the effective use of NMF in the field of oncology is also discussed.

A comprehensive appraisal of HER2 heterogeneity in HER2-amplified and HER2-low colorectal cancer.

BACKGROUND: This study aimed to elucidate the clinicopathological and molecular features of HER2-amplified and HER2-low colorectal cancers (CRCs). We also characterised HER2 expression statuses in CRCs focusing on their intratumoral heterogeneity and alterations in metastatic lesions to establish practical HER2 status assessment. METHODS: We evaluated 1009 CRCs for HER2 expression and HER2 amplification by immunohistochemistry and FISH, respectively, and correlated the results to clinicopathological and molecular data. For HER2-positive tumours, HER2 expression in metastatic lesions was also assessed. RESULTS: Twenty-five HER2-amplified (2.5%) and 46 HER2-low tumours (4.6%) were identified. HER2-amplified tumours consistently lacked a mucinous component and HER2-low tumours tended to be in the right colon, but no other clinicopathological features were noted. KRAS, NRAS or BRAF mutations were detected in only two HER2-amplified tumours (8%), whereas 23 HER2-low tumours (50%) had one of these mutations. Most HER2-amplified and HER2-low tumours showed a homogeneous or mosaic HER2 expression pattern and a clustered heterogeneous expression pattern was rather rare. HER2 expression was maintained in most metastatic lesions in both HER2-amplified (93%) and HER2-low tumours (81%). CONCLUSIONS: These results suggest that biopsy-based assessment of primary lesions is appropriate for the identification of CRC patients eligible for systemic HER2-targeted therapy.

Utility of molecular subtypes and genetic alterations for evaluating clinical outcomes in 1029 patients with endometrial cancer.

BACKGROUND: We investigated the utility of a molecular classifier tool and genetic alterations for predicting prognosis in Japanese patients with endometrial cancer. METHODS: A total of 1029 patients with endometrial cancer from two independent cohorts were classified into four molecular subtype groups. The primary and secondary endpoints were relapse-free survival (RFS) and overall survival (OS), respectively. RESULTS: Among the 265 patients who underwent initial surgery, classified according to immunohistochemistry, patients with DNA polymerase epsilon exonuclease domain mutation had an excellent prognosis (RFS and OS), patients with no specific molecular profile (NSMP) and mismatch repair protein deficiency had an intermediate prognosis, and those with protein 53 abnormal expression (p53abn) had the worst prognosis (P < 0.001). In the NSMP group, mutant KRAS and wild-type ARID1A were associated with significantly poorer 5-year RFS (41.2%) than other genomic characteristics (P < 0.001). The distribution of the subtypes differed significantly between patients with recurrence/progression and classified by sequencing (n = 764) and patients who underwent initial surgery (P < 0.001). Among patients with recurrence/progression, 51.4% had the opportunity to receive molecular targeted therapy. CONCLUSIONS: A molecular classifier is a useful tool for determining prognosis and eligibility for molecularly targeted therapy in patients with endometrial cancer.

Tumor microenvironment-mediated immune profiles and efficacy of anti-PD-L1 antibody plus chemotherapy stratified by DLL3 expression in small-cell lung cancer.

BACKGROUND: Delta-like ligand 3 (DLL3) is a therapeutic target in small-cell lung cancer (SCLC). However, how DLL3 expression status affects the tumor microenvironment (TME) and clinical outcomes in SCLC remains unclear. METHODS: This retrospective study included patients with postoperative limited-stage (LS)-SCLC and extensive-stage (ES)-SCLC treated with platinum and etoposide (PE) plus anti-programmed cell death ligand 1 (PD-L1) antibody. We investigated the relationship of DLL3 expression with TME, mutation status, tumor neoantigens, and immunochemotherapy. RESULTS: In the LS-SCLC cohort (n = 59), whole-exome sequencing revealed that DLL3High cases had significantly more neoantigens (P = 0.004) and a significantly higher rate of the signature SBS4 associated with smoking (P = 0.02) than DLL3Low cases. Transcriptome analysis in the LS-SCLC cohort revealed that DLL3High cases had significantly suppressed immune-related pathways and dendritic cell (DC) function. SCLC with DLL3High had significantly lower proportions of T cells, macrophages, and DCs than those with DLL3Low. In the ES-SCLC cohort (n = 30), the progression-free survival associated with PE plus anti-PD-L1 antibody was significantly worse in DLL3High cases than in DLL3Low cases (4.7 vs. 7.4 months, P = 0.01). CONCLUSIONS: Although SCLC with DLL3High had a higher neoantigen load, these tumors were resistant to immunochemotherapy due to suppressed tumor immunity by inhibiting antigen-presenting functions.

Clinical impact of genetic alterations of CTNNB1 in patients with grade 3 endometrial endometrioid carcinoma.

To identify prognostic factors in patients with grade 3 (high-grade) endometrial endometrioid carcinoma, we evaluated the spectrum of genomic alterations and examined whether previously reported molecular subtypes of endometrial carcinoma were adapted to clinical outcome prediction. Seventy-five Japanese patients with grade 3 endometrial endometrioid carcinoma, who underwent a potentially curative resection procedure between 1997 and 2018 at the National Cancer Center Hospital, were included. We classified the patients into four risk groups of the disease based on the Proactive Molecular Risk Classifier for Endometrial Cancer. Genomic alterations in PTEN, ARID1A, TP53, and PIK3CA were detected in more than 30% of the patients. Overall survival and recurrence-free survival of patients with genomic alterations in CTNNB1 were poorer than those of patients with wild-type CTNNB1 (p = 0.006 and p = 0.004, respectively). Compared with that of alterations prevalent in Caucasians, the frequency of genomic alterations in POLE and TP53 was higher in our study than in The Cancer Genome Atlas dataset (p = 0.01 and p = 0.01, respectively). The tendency for recurrence-free survival in the POLE exonuclease domain mutation group was better than that in the TP53 mutation and mismatch repair-deficient groups (p = 0.08 and p = 0.07, respectively), consistent with the Proactive Molecular Risk Classifier for Endometrial Cancer risk classifier definition. The CTNNB1 mutation is a potential novel biomarker for the prognosis of patients with grade 3 endometrial endometrioid carcinoma, and prognosis classification using Proactive Molecular Risk Classifier for Endometrial Cancer may help screen Japanese patients with the disease.

Clinicopathological and molecular characteristics of RSPO fusion-positive colorectal cancer.

BACKGROUND: RSPO fusions that lead to WNT pathway activation are potential therapeutic targets in colorectal cancer (CRC), but their clinicopathological significance remains unclear. METHODS: We screened 1019 CRCs for RSPO fusions using multiplex reverse transcription-PCR. The RSPO fusion-positive tumours were subjected to whole-exome sequencing (WES). RESULTS: Our analysis identified 29 CRCs with RSPO fusions (2.8%), consisting of five with an EIF3E-RSPO2 fusion and 24 with PTPRK-RSPO3 fusions. The patients were 17 women and 12 men. Thirteen tumours (45%) were right-sided. Histologically, approximately half of the tumours (13/29, 45%) had a focal or extensive mucinous component that was significantly more frequent than the RSPO fusion-negative tumours (13%; P = 8.1 × 10-7). Four tumours (14%) were mismatch repair-deficient. WES identified KRAS, BRAF, and NRAS mutations in a total of 27 tumours (93%). In contrast, pathogenic mutations in major WNT pathway genes, such as APC, CTNNB1 and RNF43, were absent. RSPO fusion status did not have a statistically significant influence on the overall or recurrence-free survival. These clinicopathological and genetic features were also confirmed in a pooled analysis of previous studies. CONCLUSION: RSPO fusion-positive CRCs constitute a rare subgroup of CRCs with several characteristic clinicopathological and genetic features.

Epithelial Expression of YAP and TAZ Is Sequentially Required in Lung Development.

TAZ (transcriptional coactivator with PDZ-binding motif) and YAP (Yes-associated protein) are key molecules of the Hippo pathway. Recent studies revealed that these molecules are essential in lung development; however, the precise signaling cascade involving these molecules and the differences in their roles during lung development remain unknown. We aimed to investigate YAP and TAZ functions using lung epithelium-specific Taz and Yap conditional knockout mice. We generated lung epithelium-specific Taz and Yap conditional knockout mice and investigated the functions of YAP and TAZ in lung development. Selective TAZ deficiency in mouse lung epithelial cells resulted in abnormal alveolarization, which mimics lung emphysema, in adults, whereas YAP deficiency caused disruption of bronchial morphogenesis during the embryonic stage. We report that TAZ and YAP are sequentially expressed in the lung and that this could explain their different phenotypes. Furthermore, we report that YAP stimulates Shh (Sonic hedgehog) expression and regulates the FGF (fibroblast growth factor)-SHH feedback loop, thereby contributing to normal bronchial morphogenesis. We also found that TGF-ß (transforming growth factor-ß) stimulation induced Shh expression in the lung epithelial cells, and both TAZ and YAP are essential in this novel pathway. Our results provide a novel insight into the molecular mechanisms underlying lung development and contribute to a better understanding of the characteristics of TAZ and YAP.

Genomic alterations in STK11 can predict clinical outcomes in cervical cancer patients.

OBJECTIVE: Cervical cancer is the fourth most common cause of cancer-related deaths in Asian women, due to its poor prognosis. This study aimed to decipher genomic alteration profiles of a cohort of Japanese cervical cancer patients to understand why certain patients benefited from molecular targeted therapies and their prognostic significance. METHODS: During 2008-2018, 154 cervical cancer patients underwent a potentially curative resection procedure at the National Cancer Center Hospital. Genomic DNA samples were analyzed using Ion AmpliSeq™ Cancer Hotspot Panel v2. Alterations in the copy number of PIK3CA, ERBB2, PTEN, and STK11 were detected using the TaqMan assay. HPV-positive results were confirmed by genomic testing and in situ hybridization assay. RESULTS: The frequency of genomic alterations in PIK3CA (36%), STK11 (16%), PTEN (11%), TP53 (11%), and KRAS (8%) was >5%. KRAS mutations were preferentially detected in patients with adenocarcinomas, and the frequency of PIK3CA mutations in patients with squamous cell carcinomas was higher than that in patients with other histological cancer types. HPV-positive results were observed in 139/154 (90.3%) patients, and TP53 mutants were detected in HPV-negative specimens. In this study, the overall survival of patients with genomic alterations in STK11 was worse than in patients with wild-type STK11 (hazard ratio = 10.6, P = 0.0079) and TCGA dataset (hazard ratio = 2.46, P = 0.029). CONCLUSIONS: More than one-third of Japanese cervical cancer patients exhibit mutations targeted by molecular targeted therapies. We have proposed the prognostic value of STK11 genomic alterations.

Histone methyltransferase SMYD2 selective inhibitor LLY-507 in combination with poly ADP ribose polymerase inhibitor has therapeutic potential against high-grade serous ovarian carcinomas.

Dysfunction of histone methylation is known to be related to cancer progression. The histone methyltransferase SMYD2 methylates histone protein H3 and non-histone proteins, including poly ADP ribose polymerase 1 (PARP1). There have been reports of SMYD2 overexpression in several types of cancers. However, there are no reports regarding its role in high-grade serous ovarian carcinomas (HGSOCs). Therefore, we investigated the expression profile and conducted functional analysis on SMYD2 in HGSOC cells. In addition, we verified whether SMYD2 inhibition increases the susceptibility of HGSOC cells to PARP inhibitors. We analyzed the expression of histone methyltransferase SMYD2 by quantitative real-time polymerase chain reaction and immunohistochemistry using HGSOC clinical tissues (n = 35). We performed functional analyses, including cell proliferation assay, cell cycle analysis, and immunoblotting, after treatment with SMYD2 siRNAs and SMYD2 selective inhibitor LLY-507 in HGSOC cells. We also performed colony-formation assay after combination treatment with LLY-507 and PARP inhibitor olaparib in HGSOC cells. The expression profiles of SMYD2 showed significant overexpression of SMYD2 in HGSOC clinical tissues. The knockdown or inhibition of SMYD2 by siRNAs or LLY-507, respectively, suppressed cell growth by increasing the proportion of apoptotic cells. LLY-507 showed additive effect with olaparib in the colony-formation assay. These findings suggest that LLY-507 can be used alone or in combination with a PARP inhibitor for the treatment of patients with HGSOC.

The histone methyltransferase WHSC1 is regulated by EZH2 and is important for ovarian clear cell carcinoma cell proliferation.

BACKGROUND: Wolf-Hirschhorn syndrome candidate gene-1 (WHSC1), a histone methyltransferase, has been found to be upregulated and its expression to be correlated with expression of enhancer of zeste homolog 2 (EZH2) in several cancers. In this study, we evaluated the role of WHSC1 and its therapeutic significance in ovarian clear cell carcinoma (OCCC). METHODS: First, we analyzed WHSC1 expression by quantitative PCR and immunohistochemistry using 23 clinical OCCC specimens. Second, the involvement of WHSC1 in OCCC cell proliferation was evaluated by MTT assays after siRNA-mediated WHSC1 knockdown. We also performed flow cytometry (FACS) to address the effect of WHSC1 on cell cycle. To examine the functional relationship between EZH2 and WHSC1, we knocked down EZH2 using siRNAs and checked the expression levels of WHSC1 and its histone mark H3K36m2 in OCCC cell lines. Finally, we checked WHSC1 expression after treatment with the selective inhibitor, GSK126. RESULTS: Both quantitative PCR and immunohistochemical analysis revealed that WHSC1 was significantly overexpressed in OCCC tissues compared with that in normal ovarian tissues. MTT assay revealed that knockdown of WHSC1 suppressed cell proliferation, and H3K36me2 levels were found to be decreased in immunoblotting. FACS revealed that WHSC1 knockdown affected the cell cycle. We also confirmed that WHSC1 expression was suppressed by EZH2 knockdown or inhibition, indicating that EZH2 is upstream of WHSC1 in OCCC cells. CONCLUSIONS: WHSC1 overexpression induced cell growth and its expression is, at least in part, regulated by EZH2. Further functional analysis will reveal whether WHSC1 is a promising therapeutic target for OCCC.

Effects of SMYD2-mediated EML4-ALK methylation on the signaling pathway and growth in non-small-cell lung cancer cells.

A specific subtype of non-small-cell lung cancer (NSCLC) characterized with an EML4-ALK fusion gene, which drives constitutive oncogenic activation of anaplastic lymphoma kinase (ALK), shows a good clinical response to ALK inhibitors. We have reported multiple examples implying the biological significance of methylation on non-histone proteins including oncogenic kinases in human carcinogenesis. Through the process to search substrates for various methyltransferases using an in vitro methyltransferase assay, we found that a lysine methyltransferase, SET and MYND domain-containing 2 (SMYD2), could methylate lysine residues 1451, 1455, and 1610 in ALK protein. Knockdown of SMYD2 as well as treatment with a SMYD2 inhibitor in two NSCLC cell lines with an EML4-ALK gene significantly attenuated the phosphorylation levels of the EML4-ALK protein. Substitutions of each of these three lysine residues to an alanine partially or almost completely diminished in vitro methylation of ALK. In addition, we found that exogenous introduction of EML4-ALK protein with the substitution of lysine 1610 to an alanine in these two cell lines reduced the phosphorylation levels of AKT, one of the downstream oncogenic molecules in the EML4-ALK pathway, and suppressed the growth of the two cell lines. We further showed that the combination of a SMYD2 inhibitor and an ALK inhibitor additively suppressed the growth of these two NSCLC cells, compared with single-agent treatment. Our results shed light on a novel mechanism that modulates the kinase activity of the ALK fused gene product and imply that SMYD2-mediated ALK methylation might be a promising target for development of a novel class of treatment for tumors with the ALK fused gene.

Dysregulation of protein methyltransferases in human cancer: An emerging target class for anticancer therapy.

Protein methylation is one of the important post-translational modifications. Although its biological and physiological functions were unknown for a long time, we and others have characterized a number of protein methyltransferases, which have unveiled the critical functions of protein methylation in various cellular processes, in particular, in epigenetic regulation. In addition, it had been believed that protein methylation is an irreversible phenomenon, but through identification of a variety of protein demethylases, protein methylation is now considered to be dynamically regulated similar to protein phosphorylation. A large amount of evidence indicated that protein methylation has a pivotal role in post-translational modification of histone proteins as well as non-histone proteins and is involved in various processes of cancer development and progression. As dysregulation of this modification has been observed frequently in various types of cancer, small-molecule inhibitors targeting protein methyltransferases and demethylases have been actively developed as anticancer drugs; clinical trials for some of these drugs have already begun. In this review, we discuss the biological and physiological importance of protein methylation in human cancer, especially focusing on the significance of protein methyltransferases as emerging targets for anticancer therapy.

Recognition of modification status on a histone H3 tail by linked histone reader modules of the epigenetic regulator UHRF1.

Multiple covalent modifications on a histone tail are often recognized by linked histone reader modules. UHRF1 [ubiquitin-like, containing plant homeodomain (PHD) and really interesting new gene (RING) finger domains 1], an essential factor for maintenance of DNA methylation, contains linked two-histone reader modules, a tandem Tudor domain and a PHD finger, tethered by a 17-aa linker, and has been implicated to link histone modifications and DNA methylation. Here, we present the crystal structure of the linked histone reader modules of UHRF1 in complex with the amino-terminal tail of histone H3. Our structural and biochemical data provide the basis for combinatorial readout of unmodified Arg-2 (H3-R2) and methylated Lys-9 (H3-K9) by the tandem tudor domain and the PHD finger. The structure reveals that the intermodule linker plays an essential role in the formation of a histone H3-binding hole between the reader modules by making extended contacts with the tandem tudor domain. The histone H3 tail fits into the hole by adopting a compact fold harboring a central helix, which allows both of the reader modules to simultaneously recognize the modification states at H3-R2 and H3-K9. Our data also suggest that phosphorylation of a linker residue can modulate the relative position of the reader modules, thereby altering the histone H3-binding mode. This finding implies that the linker region plays a role as a functional switch of UHRF1 involved in multiple regulatory pathways such as maintenance of DNA methylation and transcriptional repression.