Background Checkpoint-blockade immunotherapy targeting programmed cell death protein 1 (PD-1) has recently shown promising efficacy in hepatocellular carcinoma (HCC). However, the factors affecting and predicting the response to anti-PD-1 immunotherapy in HCC are still unclear. Herein, we report the dynamic variation characteristics and specificities of the gut microbiome during anti-PD-1 immunotherapy in HCC using metagenomic sequencing. Results Fecal samples from patients responding to immunotherapy showed higher taxa richness and more gene counts than those of non-responders. For dynamic analysis during anti-PD-1 immunotherapy, the dissimilarity of beta diversity became prominent across patients as early as Week 6. In non-responders, Proteobacteria increased from Week 3, and became predominant at Week 12. Twenty responder-enriched species, including Akkermansia muciniphila and Ruminococcaceae spp., were further identified. The related functional genes and metabolic pathway analysis, such as carbohydrate metabolism and methanogenesis, verified the potential bioactivities of responder-enriched species. Conclusions Gut microbiome may have a critical impact on the responses of HCC patients treated with anti-PD-1 immunotherapy. The dynamic variation characteristics of the gut microbiome may provide early predictions of the outcomes of immunotherapy in HCC, which is critical for disease-monitoring and treatment decision-making.

Synchronous multifocal tumors are common in the hepatobiliary and pancreatic system but because of similarities in their histological features, oncologists have difficulty in identifying their precise tissue clonal origin through routine histopathological methods. To address this problem and assist in more precise diagnosis, we developed a computational approach for tissue origin diagnosis based on naive Bayes algorithm (TOD‐Bayes) using ubiquitous RNA‐Seq data. Massive tissue‐specific RNA‐Seq data sets were first obtained from The Cancer Genome Atlas (TCGA) and ∼1,000 feature genes were used to train and validate the TOD‐Bayes algorithm. The accuracy of the model was >95% based on tenfold cross validation by the data from TCGA. A total of 18 clinical cancer samples (including six negative controls) with definitive tissue origin were subsequently used for external validation and 17 of the 18 samples were classified correctly in our study (94.4%). Furthermore, we included as cases studies seven tumor samples, taken from two individuals who suffered from synchronous multifocal tumors across tissues, where the efforts to make a definitive primary cancer diagnosis by traditional diagnostic methods had failed. Using our TOD‐Bayes analysis, the two clinical test cases were successfully diagnosed as pancreatic cancer (PC) and cholangiocarcinoma (CC), respectively, in agreement with their clinical outcomes. Based on our findings, we believe that the TOD‐Bayes algorithm is a powerful novel methodology to accurately identify the tissue origin of synchronous multifocal tumors of unknown primary cancers using RNA‐Seq data and an important step toward more precision‐based medicine in cancer diagnosis and treatment.

Background: Drug resistance is a major obstacle to effective cancer therapy. In order to detect the change in tumor genomic states under drug selection pressure, we use next-generation sequencing technology to investigate the underlying potential mechanisms of drug resistance. Methods: In our study, we presented a bladder cancer patient who had been a bona fide responder to first-line gemcitabine plus cisplatin regimen and second-line pazopanib (tyrosine kinase inhibitor (TKI) for FGFR3-TACC3 fusion) but finally had disease progression as an ideal case for showing genomic alteration during drug resistance. We applied whole-exome sequencing and ultra-deep target sequencing to the patient pre- and post- pazopanib resistance. Protein-protein interaction (PPI) network and Gene Ontology (GO) analyses were used to analysis protein interactions and genomic alterations. Patient-derived xenograft (PDX) model was built to test drug sensitivity. Results: Twelve mutations scattered in 12 genes were identified by WES pre- pazopanib resistance, while 63 mutations in 50 genes arose post- pazopanib resistance. PPI network showed proteins from multiple epigenetic regulator families were involved post- pazopanib resistance, including subunits of chromatin remodeler SWI/SNF complex ARID1A/1B and SMARCA4, histone acetylation writers CREBBP, histone methylation writer NSD1 and erasers KDM6A/5A. GO enrichment analysis showed pazopanib resistance genes were prominently tagged for chromatin modification, transcription, as well as gland development, leaving genes with the best adaptive FGFR TKI-coping mechanisms. In addition, significantly elevated tumor mutational burden suggested possible utility of immunotherapy. Intriguingly, PDX model suggested that, sensitivity to original chemotherapy regimen (cisplatin) was restored in patient tumor post-pazopanib. Conclusions: Epigenetic regulation may play a role in acquired TKI resistance. Our study traced the complete tumor genomic variation course from chemo-resistant but TKI-sensitive to TKI-resistant but chemo-(re) sensitive, revealing the potential complex dynamic drug-driven mechanisms of resistance.