Background and Hypotheses: Mongolia has the highest reported incidence of - and mortality from - hepatocellular carcinoma (HCC) in the world, which is between three and seven times higher than that observed in other high-incidence populations, such as South Korea, Thailand, and China. Although its causative factors and underlying tumor biology remain unknown, risk factors such as chronic hepatitis B virus and/or hepatitis C virus infection are highly prevalent. Moreover, Mongolia has an extraordinarily high prevalence of hepatitis delta virus (HDV), more than ten times above the global estimate, which we hypothesize as a major driver of Mongolian hepatocarcinogenesis. By performing a comprehensive and integrative genomic characterization, we aimed to identify robust molecular subclasses and driver features informative of the etiology and progression of the disease. This study, the first one of this nature ever to be conducted on a Mongolian cohort, identifies similarities to previous studies performed on other populations (which reflect general characteristics of liver carcinogenesis) as well as unique features of Mongolian HCC. Study Design and Methods: Clinical information and paired tumor/non-tumor liver tissue samples were obtained from 76 HCC patients undergoing surgery at the National Cancer Center of Mongolia. Transcriptome sequencing and whole exome sequencing were performed on all tumor and non-tumor samples, followed by bioinformatics processing and quality control. Consensus clustering and regularized Cox regression analyses were performed on transcriptome data. Driver mutations and mutational signatures were assessed. These results were compared to those from 373 HCC patients of different races and ethnicities and diverse etiologies. Furthermore, molecular subclass gene signatures from multiple previous studies were assessed for validation. Results and Conclusions: Using a transcriptomics-based unsupervised approach, we identified four molecular subclasses among Mongolian HCC cases, which were associated with different patient survival outcomes. These molecular subclasses were mapped to signatures reported in other HCC cohorts, showing significant associations consistent with survival outcomes. A supervised analysis determined risk scores strongly associated with survival, consistent with the four molecular subclasses. Mutational signature analysis using recently published COSMIC and Environmental Agents Compendia identified signatures with distinct prevalence among HDV+ and HVD- groups. Signatures differentially associated with HDV+ include mutational patterns linked to alkylating agents (such as temozolomide), tobacco chewing and exposures to 1,8-Dinitropyrene and furan, whereas HDV- appeared differentially associated with aristolochic acid II. In addition to common driver mutations (TP53, CTNNB1) frequently found in pan-cancer analysis, Mongolian HCC exhibits several unique drivers (most notably GTF2IRD2B, PNRC2, and SPTA1), the latter of which may be linked to HDV infection. These results suggest the existence of new molecular mechanisms at play in Mongolian hepatocarcinogenesis.

The Genomic Landscape of Mongolian Hepatocellular Carcinoma