ABSTRACT
Molecular heterogeneity of hepatobiliary tumor including intertumoral and intratumoral disparity always leads to drug resistance. Here, seven hepatobiliary tumor organoids are generated to explore heterogeneity and evolution via single-cell RNA sequencing. HCC272 with high status of epithelia-mesenchymal transition proves broad-spectrum drug resistance. By examining the expression pattern of cancer stem cells markers (e.g., PROM1, CD44, and EPCAM), it is found that CD44 positive population may render drug resistance in HCC272. UMAP and pseudo-time analysis identify the intratumoral heterogeneity and distinct evolutionary trajectories, of which catenin beta-1 (CTNNB1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and nuclear paraspeckle assembly transcript 1 (NEAT1) advantage expression clusters are commonly shared across hepatobiliary organoids. CellphoneDB analysis further implies that metabolism advantage organoids with enrichment of hypoxia signal upregulate NEAT1 expression in CD44 subgroup and mediate drug resistance that relies on Jak-STAT pathway. Moreover, metabolism advantage clusters shared in several organoids have similar characteristic genes (GAPDH, NDRG1 (N-Myc downstream regulated 1), ALDOA, and CA9). The combination of GAPDH and NDRG1 is an independent risk factor and predictor for patient survival. This study delineates heterogeneity of hepatobiliary tumor organoids and proposes that the collaboration of intratumoral heterogenic subpopulations renders malignant phenotypes and drug resistance.