Integrated genomic/epigenomic analysis stratifies subtypes of clear cell ovarian carcinoma, highlighting their cellular origin.

The cellular origin of clear cell ovarian carcinoma (CCOC), a major histological subtype of ovarian carcinoma remains elusive. Here, we explored the candidate cellular origin and identify molecular subtypes using integrated genomic/epigenomic analysis. We performed whole exome-sequencing, microarray, and DNA methylation array in 78 CCOC samples according to the original diagnosis. The findings revealed that ARID1A and/or PIK3CA mutations were mutually exclusive with DNA repair related genes, including TP53, BRCA1, and ATM. Clustering of CCOC and other ovarian carcinomas (n = 270) with normal tissues from the fallopian tube, ovarian surface epithelium, endometrial epithelium, and pelvic peritoneum mesothelium (PPM) in a methylation array showed that major CCOC subtypes (with ARID1A and/or PIK3CA mutations) were associated with the PPM-lile cluster (n = 64). This cluster was sub-divided into three clusters: (1) mismatch repair (MMR) deficient with tumor mutational burden-high (n = 2), (2) alteration of ARID1A (n = 51), and (3) ARID1A wild-type (n = 11). The remaining samples (n = 14) were subdivided into (4) ovarian surface epithelium-like (n = 11) and (5) fallopian tube-like (considered as high-grade serous histotype; n = 3). Among these, subtypes (1-3) and others (4 and 5) were found to be associated with immunoreactive signatures and epithelial-mesenchymal transition, respectively. These results contribute to the stratification of CCOC into biological subtypes.

Molecular analysis of non-small cell lung cancer using a dual-targeted DNA and RNA comprehensive genomic profiling panel.

BACKGROUND: Comprehensive cancer genomic profiling tests have recently been used clinically to guide optimal treatment. Currently approved tests use DNA from tissue or plasma samples to analyze a few hundred genes. RNA panels complement DNA panels to detect fusion and exon skipping. METHODS: Between April 2017 and March 2022, we analyzed non-small cell lung cancer samples using Todai OncoPanel, a matched tumor/normal pair panel targeting both DNA and RNA. Publicly available genomic data was downloaded from the Center for Cancer Genomics and Advanced Therapeutics database on 2022/11/3. RESULTS: Sixty non-small cell lung cancer (NSCLC) samples were analyzed. With the DNA panel, 32 samples (53%) had TP53 loss-of-function mutations. Among adenocarcinoma, 17 (33%) had EGFR activating mutations, and 6 (12%) had ERBB2 activating mutations. One BRCA1 and one BRCA2 pathogenic germline variant were also detected. With the RNA panel, 11 fusion genes were detected, all in adenocarcinoma. Specifically, EML4-ALK and KIF5B-RET were detected from one sample each, and 9 others were all novel fusions with unknown pathogenicity. In addition, 4 of 60 (7%) NSCLC samples harbored MET exon 14 skipping. Analysis of the Center for Cancer Genomics and Advanced Therapeutics database found 37 MET exon 14 splice site mutations in 1514 NSCLC samples (2%, p = 0.039). CONCLUSIONS: Analysis of NSCLC with Todai OncoPanel detected many druggable targets. Its RNA panel may detect MET exon 14 skipping with high sensitivity.

Region-specific DNA hydroxymethylation along the malignant progression of IDH-mutant gliomas.

The majority of low-grade isocitrate dehydrogenase-mutant (IDHmt) gliomas undergo malignant progression (MP), but their underlying mechanism remains unclear. IDHmt gliomas exhibit global DNA methylation, and our previous report suggested that MP could be partly attributed to passive demethylation caused by accelerated cell cycles. However, during MP, there is also active demethylation mediated by ten-eleven translocation, such as DNA hydroxymethylation. Hydroxymethylation is reported to potentially contribute to gene expression regulation, but its role in MP remains under investigation. Therefore, we conducted a comprehensive analysis of hydroxymethylation during MP of IDHmt astrocytoma. Five primary/malignantly progressed IDHmt astrocytoma pairs were analyzed with oxidative bisulfite and the Infinium EPIC methylation array, detecting 5-hydroxymethyl cytosine at over 850,000 locations for region-specific hydroxymethylation assessment. Notably, we observed significant sharing of hydroxymethylated genomic regions during MP across the samples. Hydroxymethylated CpGs were enriched in open sea and intergenic regions (p < 0.001), and genes undergoing hydroxymethylation were significantly associated with cancer-related signaling pathways. RNA sequencing data integration identified 91 genes with significant positive/negative hydroxymethylation-expression correlations. Functional analysis suggested that positively correlated genes are involved in cell-cycle promotion, while negatively correlated ones are associated with antineoplastic functions. Analyses of The Cancer Genome Atlas clinical data on glioma were in line with these findings. Motif-enrichment analysis suggested the potential involvement of the transcription factor KLF4 in hydroxymethylation-based gene regulation. Our findings shed light on the significance of region-specific DNA hydroxymethylation in glioma MP and suggest its potential role in cancer-related gene expression and IDHmt glioma malignancy.

Analysis of quality metrics in comprehensive cancer genomic profiling using a dual DNA-RNA panel.

Background: The nucleic acid quality from formalin-fixed paraffin-embedded (FFPE) tumor vary among samples, resulting in substantial variability in the quality of comprehensive cancer genomic profiling tests. The objective of the study is to investigate how nucleic acid quality affects sequencing quality. We also examined the variations in nucleic acid quality among different hospitals or cancer types. Methods: Three nucleic acid quality metrics (ddCq, Q-value, and DV200) and five sequencing quality metrics (on-target rate, mean depth, coverage uniformity, target exon coverage, and coverage of the housekeeping gene) were examined using 585 samples from the Todai OncoPanel, a dual DNA-RNA panel. Results: In the DNA panel, ddCq served as an indicator of sequencing depth and Q-value reflected the uniformity of sequencing across different regions. It was essential to have favorable values not only for ddCq but also for Q-value to obtain ideal sequencing results. For the RNA panel, DV200 proved to be a valuable metric for assessing the coverage of the housekeeping genes. Significant inter-hospital differences were observed for DNA quality (ddCq and Q-value), but not for RNA quality (DV200). Differences were also observed among cancer types, with Q-value being the lowest in lung and the highest in cervix, while DV200 was the highest in lung and the lowest in bowel. Conclusions: We demonstrated distinct characteristics and high predictive performances of ddCq, Q-value, and DV200. Variations were observed in the nucleic acid quality across hospitals and cancer types. Further study is warranted on preanalytical factors in comprehensive cancer genomic profiling tests.

Transcriptomic and epigenetic dissection of spinal ependymoma (SP-EPN) identifies clinically relevant subtypes enriched for tumors with and without NF2 mutation.

Ependymomas encompass multiple clinically relevant tumor types based on localization and molecular profiles. Tumors of the methylation class "spinal ependymoma" (SP-EPN) represent the most common intramedullary neoplasms in children and adults. However, their developmental origin is ill-defined, molecular data are scarce, and the potential heterogeneity within SP-EPN remains unexplored. The only known recurrent genetic events in SP-EPN are loss of chromosome 22q and NF2 mutations, but neither types and frequency of these alterations nor their clinical relevance have been described in a large, epigenetically defined series. Transcriptomic (n = 72), epigenetic (n = 225), genetic (n = 134), and clinical data (n = 112) were integrated for a detailed molecular overview on SP-EPN. Additionally, we mapped SP-EPN transcriptomes to developmental atlases of the developing and adult spinal cord to uncover potential developmental origins of these tumors. The integration of transcriptomic ependymoma data with single-cell atlases of the spinal cord revealed that SP-EPN display the highest similarities to mature adult ependymal cells. Unsupervised hierarchical clustering of transcriptomic data together with integrated analysis of methylation profiles identified two molecular SP-EPN subtypes. Subtype A tumors primarily carried previously known germline or sporadic NF2 mutations together with 22q loss (bi-allelic NF2 loss), resulting in decreased NF2 expression. Furthermore, they more often presented as multilocular disease and demonstrated a significantly reduced progression-free survival as compared to SP-EP subtype B. In contrast, subtype B predominantly contained samples without NF2 mutation detected in sequencing together with 22q loss (monoallelic NF2 loss). These tumors showed regular NF2 expression but more extensive global copy number alterations. Based on integrated molecular profiling of a large multi-center cohort, we identified two distinct SP-EPN subtypes with important implications for genetic counseling, patient surveillance, and drug development priorities.

Case report and literature review: exploration of molecular therapeutic targets in recurrent malignant meningioma through comprehensive genetic analysis with Todai OncoPanel.

Background: Despite accumulating research on the molecular characteristics of meningiomas, no definitive molecularly targeted therapy for these tumors has been established to date. Molecular mechanisms underlying meningioma progression also remain unclear. Comprehensive genetic testing approaches can reveal actionable gene aberrations in meningiomas. However, there is still limited information on whether profiling the molecular status of subsequent recurrent meningiomas could influence the choice of molecular-targeted therapies. Case presentation: We report a case of meningioma with malignant progression and multiple recurrences. We performed matched tumor pair analysis using the Todai OncoPanel to investigate the possibility of additional standard treatments. The loss of several chromosomal regions, including NF2 and CDKN2A, which is associated with aggressive meningiomas, was considered a significant driver event for malignant progression. Using additional matched tumor pair analysis, mutations in TRAF7, ARID1A, and ERBB3 were identified as subclonal driver events at the time of recurrence. No genetic aberrations were found for which evidence-based targeted therapy was applicable. We also reviewed previous reports of molecular therapies in meningioma to discuss issues with the current molecular testing approach. Conclusion: Gene panel testing platforms such as the Todai OncoPanel represent a powerful approach to elucidate actionable genetic alterations in various types of tumors, although their use is still limited to the diagnosis and prediction of prognosis in meningiomas. To enable targeted molecular therapy informed by gene-panel testing, further studies including matched tumor pair analyses are required to understand the molecular characteristics of meningiomas and develop treatments based on genetic abnormalities.

Gut insulin action protects from hepatocarcinogenesis in diabetic mice comorbid with nonalcoholic steatohepatitis.

Diabetes is known to increase the risk of nonalcoholic steatohepatitis (NASH) and hepatocellular carcinoma (HCC). Here we treat male STAM (STelic Animal Model) mice, which develop diabetes, NASH and HCC associated with dysbiosis upon low-dose streptozotocin and high-fat diet (HFD), with insulin or phlorizin. Although both treatments ameliorate hyperglycemia and NASH, insulin treatment alone lead to suppression of HCC accompanied by improvement of dysbiosis and restoration of antimicrobial peptide production. There are some similarities in changes of microflora from insulin-treated patients comorbid with diabetes and NASH. Insulin treatment, however, fails to suppress HCC in the male STAM mice lacking insulin receptor specifically in intestinal epithelial cells (ieIRKO), which show dysbiosis and impaired gut barrier function. Furthermore, male ieIRKO mice are prone to develop HCC merely on HFD. These data suggest that impaired gut insulin signaling increases the risk of HCC, which can be countered by restoration of insulin action in diabetes.

Association of systemic adverse reaction patterns with long-term dynamics of humoral and cellular immunity after coronavirus disease 2019 third vaccination.

The objective of this study was to clarify the impact of adverse reactions on immune dynamics. We investigated the pattern of systemic adverse reactions after the second and third coronavirus disease 2019 (COVID-19) vaccinations and their relationship with immunoglobulin G against severe acute respiratory syndrome coronavirus 2 spike 1 protein titers, neutralizing antibody levels, peak cellular responses, and the rate of decrease after the third vaccination in a large-scale community-based cohort in Japan. Participants who received a third vaccination with BNT162b2 (Pfizer/BioNTech) or mRNA-1273 (Moderna), had two blood samples, had not had COVID-19, and had information on adverse reactions after the second and third vaccinations (n = 2198) were enrolled. We collected data on sex, age, adverse reactions, comorbidities, and daily medicine using a questionnaire survey. Participants with many systemic adverse reactions after the second and third vaccinations had significantly higher humoral and cellular immunity in the peak phase. Participants with multiple systemic adverse reactions after the third vaccination had small changes in the geometric values of humoral immunity and had the largest geometric mean of cellar immunity in the decay phase. Systemic adverse reactions after the third vaccination helped achieve high peak values and maintain humoral and cellular immunity. This information may help promote uptake of a third vaccination, even among those who hesitate due to adverse reactions.

Multiancestry genomic and transcriptomic analysis of gastric cancer.

Gastric cancer is among the most common malignancies worldwide, characterized by geographical, epidemiological and histological heterogeneity. Here, we report an extensive, multiancestral landscape of driver events in gastric cancer, involving 1,335 cases. Seventy-seven significantly mutated genes (SMGs) were identified, including ARHGAP5 and TRIM49C. We also identified subtype-specific drivers, including PIGR and SOX9, which were enriched in the diffuse subtype of the disease. SMGs also varied according to Epstein-Barr virus infection status and ancestry. Non-protein-truncating CDH1 mutations, which are characterized by in-frame splicing alterations, targeted localized extracellular domains and uniquely occurred in sporadic diffuse-type cases. In patients with gastric cancer with East Asian ancestry, our data suggested a link between alcohol consumption or metabolism and the development of RHOA mutations. Moreover, mutations with potential roles in immune evasion were identified. Overall, these data provide comprehensive insights into the molecular landscape of gastric cancer across various subtypes and ancestries.

Loss of viral genome with altered immune microenvironment during tumour progression of Epstein-Barr virus-associated gastric carcinoma.

Epstein-Barr virus (EBV) is one of the major drivers of gastric carcinogenesis. EBV infection is established before tumour initiation and is generally maintained throughout tumour development; however, the significance of EBV in tumour maintenance and progression remains to be elucidated. Here, we report eight cases of EBV-associated gastric carcinoma (EBVaGC) with intratumoural heterogenous expression of EBV-encoded small RNA (EBER), a highly expressed latent gene of EBV, and demonstrate clinicopathological characteristics of these rare cases. By performing detailed histological assessment of EBER-positive and -negative components of each case, detection of EBV genome in tumour cells by fluorescence in situ hybridisation, TP73 methylation analysis, whole exome sequencing, and targeted gene panel sequencing, we identified tumours in two patients to be collision tumours of different origins. In the other six patients, some genetic/epigenetic alterations were shared between EBER-positive and -negative components, suggesting that EBV was eliminated from tumour cells during progression. Interestingly, in both tumour types, programmed death ligand 1 and intratumoural infiltration of CD8+ T lymphocytes were lower in EBER-negative than in EBER-positive components, suggesting an immunogenic role of EBV. To the best of our knowledge, this study is the first to demonstrate the detailed histological features and genetic/epigenetic alterations in EBVaGC with heterogenous EBER expression; the loss of EBV may benefit tumour progression and immune evasion and might be clinically important for selecting treatment strategies for such cancers. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Clinical utility of Todai OncoPanel in the setting of approved comprehensive cancer genomic profiling tests in Japan.

Comprehensive cancer genome profiling (CGP) has been nationally reimbursed in Japan since June 2019. Less than 10% of the patients have been reported to undergo recommended treatment. Todai OncoPanel (TOP) is a dual DNA-RNA panel as well as a paired tumor-normal matched test. Two hundred patients underwent TOP as part of Advanced Medical Care B with approval from the Ministry of Health, Labour and Welfare between September 2018 and December 2019. Tests were carried out in patients with cancers without standard treatment or when patients had already undergone standard treatment. Data from DNA and RNA panels were analyzed in 198 and 191 patients, respectively. The percentage of patients who were given therapeutic or diagnostic recommendations was 61% (120/198). One hundred and four samples (53%) harbored gene alterations that were detected with the DNA panel and had potential treatment implications, and 14 samples (7%) had a high tumor mutational burden. Twenty-two samples (11.1%) harbored 30 fusion transcripts or MET exon 14 skipping that were detected by the RNA panel. Of those 30 transcripts, 6 had treatment implications and 4 had diagnostic implications. Thirteen patients (7%) were found to have pathogenic or likely pathogenic germline variants and genetic counseling was recommended. Overall, 12 patients (6%) received recommended treatment. In summary, patients benefited from both TOP DNA and RNA panels while following the same indication as the approved CGP tests. (UMIN000033647).

Multi-omics analysis defines highly refractory RAS burdened immature subgroup of infant acute lymphoblastic leukemia.

KMT2A-rearranged infant acute lymphoblastic leukemia (ALL) represents the most refractory type of childhood leukemia. To uncover the molecular heterogeneity of this disease, we perform RNA sequencing, methylation array analysis, whole exome and targeted deep sequencing on 84 infants with KMT2A-rearranged leukemia. Our multi-omics clustering followed by single-sample and single-cell inference of hematopoietic differentiation establishes five robust integrative clusters (ICs) with different master transcription factors, fusion partners and corresponding stages of B-lymphopoietic and early hemato-endothelial development: IRX-type differentiated (IC1), IRX-type undifferentiated (IC2), HOXA-type MLLT1 (IC3), HOXA-type MLLT3 (IC4), and HOXA-type AFF1 (IC5). Importantly, our deep mutational analysis reveals that the number of RAS pathway mutations predicts prognosis and that the most refractory subgroup of IC2 possesses 100% frequency and the heaviest burden of RAS pathway mutations. Our findings highlight the previously under-appreciated intra- and inter-patient heterogeneity of KMT2A-rearranged infant ALL and provide a rationale for the future development of genomics-guided risk stratification and individualized therapy.

The serological diversity of serum IgG/IgA/IgM against SARS-CoV-2 nucleoprotein, spike, and receptor-binding domain and neutralizing antibodies in patients with COVID-19 in Japan.

Background: We compared the temporal changes of immunoglobulin M (IgM), IgG, and IgA antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleoprotein (N), spike 1 subunit (S1), and receptor-binding domain (RBD), and neutralizing antibodies (NAbs) against SARS-CoV-2 in patients with coronavirus disease 2019 (COVID-19) to understand the humoral immunity in COVID-19 patients for developing drugs and vaccines for COVID-19. Methods: A total of five confirmed COVID-19 cases in Nissan Tamagawa Hospital in early August 2020 were recruited in this study. Using a fully automated chemiluminescence immunoassay analyzer, we measured the levels of IgG, IgA, and IgM against SARS-CoV-2 N, S1, and RBD and NAbs against SARS-CoV-2 in COVID-19 patients' sera acquired multiple times in individuals from 0 to 76 days after symptom onset. Results: IgG levels against SARS-CoV-2 structural proteins increased over time in all cases but IgM and IgA levels against SARS-CoV-2 showed different increasing trends among individuals in the early stage. In particular, we observed IgA increasing before IgG and IgM in some cases. The NAb levels were more than cut-off value in 4/5 COVID-19 patients some of whose antibodies against RBD did not exceed the cut-off value in the early stage. Furthermore, NAb levels against SARS-CoV-2 increased and kept above cut-off value more than around 70 days after symptom onset in all cases. Conclusion: Our findings indicate COVID-19 patients should be examined for IgG, IgA, and IgM against SARS-CoV-2 structural proteins and NAbs against SARS-CoV-2 to analyze the diversity of patients' immune mechanisms.

Tumor mutational burden measurement using comprehensive genomic profiling assay.

BACKGROUND: Tumors with a high number of mutations in the genome, or tumor mutational burden, are presumed to be more likely to respond to immune checkpoint inhibitors. However, the optimal method to calculate tumor mutational burden using comprehensive genomic profiling assays is unknown. METHODS: Todai OncoPanel is a dual panel of a deoxyribonucleic acid panel and a ribonucleic acid panel. Todai OncoPanel deoxyribonucleic acid panel version 6 is an improvement over version 3 with increased number of targeted genes and limited targeting of intronic regions. We calculated tumor mutational burden measured by Todai OncoPanel deoxyribonucleic acid panel versions 3 and 6 using three different calculation methods: all mutations within the targeted region (target tumor mutational burden), all mutations within the coding region (all coding tumor mutational burden) and non-synonymous mutations (non-synonymous coding tumor mutational burden). We then compared them with whole exosome sequencing tumor mutational burden. In addition, 16 lung cancer patients whose samples were analyzed using Todai OncoPanel deoxyribonucleic acid version 3 were treated with anti-PD-1 or PD-L1 antibody monotherapy. RESULTS: When compared with whole exosome sequencing tumor mutational burden as the standard, tumor mutational burden measured by Todai OncoPanel deoxyribonucleic acid version 3 resulted in accuracy of 71% for all three calculation methods. In version 6, accuracy was 96% for target tumor mutational burden and all coding tumor mutational burden and 91% for non-synonymous coding tumor mutational burden. Patients with either partial response or stable disease had higher non-synonymous coding tumor mutational burden (6.7/Mb vs. 1.6/Mb, P = 0.02) and higher PD-L1 expression (40% vs. 3%, P = 0.01) and a trend toward higher target tumor mutational burden (9.2/Mb vs. 2.4/Mb, P = 0.09) compared with patients with progressive disease. CONCLUSIONS: Increase in targeted gene number and limiting intronic regions improved tumor mutational burden measurement by Todai OncoPanel when compared with whole exosome sequencing tumor mutational burden. Target tumor mutational burden may be the method of choice to measure tumor mutational burden.

Genetic diagnosis of pseudomyxoma peritonei originating from mucinous borderline tumor inside an ovarian teratoma.

BACKGROUND: Pseudomyxoma peritonei is a rare disease condition mainly caused by primary mucinous tumors from the appendix and rarely from the ovary, such as when mucinous ovarian tumors arise from within a teratoma. Molecular analyses of pseudomyxoma from the appendix showed that KRAS and GNAS pathogenic variants are common genetic features of pseudomyxoma peritonei. However, the origin of the tumors is difficult to be identified via genetic variants alone. This study presents a case of pseudomyxoma peritonei of ovarian origin, which was diagnosed by comprehensive genomic profiling with ploidy analysis in a series of primary, recurrent, and autopsy tumor specimens. CASE PRESENTATION: A 40-year-old woman was diagnosed with Stage IC2 mucinous ovarian tumor of borderline malignancy with mature cystic teratoma, upon clinical pathology. Immunohistochemical analysis suggested that the mucinous tumor was derived from the intestinal component of an ovarian teratoma. Three years later, intraperitoneal recurrence was detected, which subsequently progressed to pseudomyxoma peritonei. Genomic analysis detected KRAS (G12D), GNAS (R201C), and FBXW7 (R367*) variants in the primary tumor. In addition, the tumor showed aneuploidy with loss of heterozygosity (LOH) in all its chromosomes, which suggested that the primary ovarian tumor was derived from germ cells. Existence of one Barr body suggested the existence of uniparental disomy of the tumors throughout the genome, instead of a haploid genotype. All three pathogenic variants remained positive in the initial recurrent tumor, as well as in the paired DNA from the whole blood in pseudomyxoma peritonei. The pathogenic variant of KRAS (G12D) was also identified in the autopsy specimen of the appendix by droplet digital polymerase chain reaction. CONCLUSIONS: This study pathologically and genetically confirmed that the primary ovarian borderline tumor was derived from the intestinal component of an ovarian teratoma, and that the subsequent pseudomyxoma peritonei progressed from the primary ovarian tumor. Integrative genomic analysis was useful to identify cellular origin of tumors, as well as to precisely interpret the process of disease progression.

Composite Clonal Analysis Reveals Transition of NSCLC Subtypes Through Accumulation of Gene Mutations: A Case Report.

We analyzed an EGFR-mutated lung cancer with a pathologic diagnosis of combined large cell neuroendocrine carcinoma with mixed adenocarcinoma subtypes. Targeted next-generation sequencing of each component suggested that mutations in RB1, TP53, and SMAD4 and apparent loss of heterozygosity of TP53 and SMAD4 accompanied the transition of different adenocarcinoma subtypes. Additional gene mutations including PTEN, MST1R, and PIK3CA were noted during transdifferentiation from acinar adenocarcinoma to large cell neuroendocrine carcinoma. Combined DNA and RNA analysis using Todai OncoPanel revealed that transdifferentiation to different pathologic subtypes occurred in a single tumor through the accumulation of gene mutations.

Genetic and epigenetic basis of hepatoblastoma diversity.

Hepatoblastoma (HB) is the most common pediatric liver malignancy; however, hereditary predisposition and acquired molecular aberrations related to HB clinicopathological diversity are not well understood. Here, we perform an integrative genomic profiling of 163 pediatric liver tumors (154 HBs and nine hepatocellular carcinomas) based on the data acquired from a cohort study (JPLT-2). The total number of somatic mutations is precious low (0.52/Mb on exonic regions) but correlated with age at diagnosis. Telomerase reverse transcriptase (TERT) promoter mutations are prevalent in the tween HBs, selective in the transitional liver cell tumor (TLCT, > 8 years old). DNA methylation profiling reveals that classical HBs are characterized by the specific hypomethylated enhancers, which are enriched with binding sites for ASCL2, a regulatory transcription factor for definitive endoderm in Wnt-pathway. Prolonged upregulation of ASCL2, as well as fetal-liver-like methylation patterns of IGF2 promoters, suggests their "cell of origin" derived from the premature hepatoblast, similar to intestinal epithelial cells, which are highly proliferative. Systematic molecular profiling of HB is a promising approach for understanding the epigenetic drivers of hepatoblast carcinogenesis and deriving clues for risk stratification.

Accumulation of Molecular Aberrations Distinctive to Hepatocellular Carcinoma Progression.

Cancer develops through the accumulation of genetic and epigenetic aberrations. To identify sequential molecular alterations that occur during the development of hepatocellular carcinoma (HCC), we compared 52 early and 108 overt HCC samples by genome sequencing. Gene mutations in the p53/RB1 pathway, WNT pathway, MLL protein family, SWI/SNF complexes, and AKT/PI3K pathway were common in HCC. In the early phase of all entities, TERT was the most frequently upregulated gene owing to diverse mechanisms. Despite frequent somatic mutations in driver genes, including CTNNB1 and TP53, early HCC was a separate molecular entity from overt HCC, as each had a distinct expression profile. Notably, WNT target genes were not activated in early HCC regardless of CTNNB1 mutation status because ß-catenin did not translocate into the nucleus due to the E-cadherin/ß-catenin complex at the membrane. Conversely, WNT targets were definitively upregulated in overt HCC, with CTNNB1 mutation associated with downregulation of CDH1 and hypomethylation of CpG islands in target genes. Similarly, cell-cycle genes downstream of the p53/RB pathway were upregulated only in overt HCC, with TP53 or RB1 gene mutations associated with chromosomal deletion of 4q or 16q. HCC was epigenetically distinguished into four subclasses: normal-like methylation, global-hypomethylation (favorable prognosis), stem-like methylation (poor prognosis), and CpG island methylation. These methylation statuses were globally maintained through HCC progression. Collectively, these data show that as HCC progresses, additional molecular events exclusive of driver gene mutations cooperatively contribute to transcriptional activation of downstream targets according to methylation status. SIGNIFICANCE: In addition to driver gene mutations in the WNT and p53 pathways, further molecular events are required for aberrant transcriptional activation of these pathways as HCC progresses.