Whole-transcriptome analysis of hepatocellular carcinoma

Glycogene Expression Profiling of Hepatic Cells by RNA-Seq Analysis for Glyco-Biomarker Identification

Glycans are primarily generated by “glycogenes,” which consist of more than 200 genes for glycosynthesis, including sugar-nucleotide synthases, sugar-nucleotide transporters, and glycosyltransferases. Measuring the expression level of glycogenes is one of the approaches to analyze the glycomes of particular biological and clinical samples. To develop an effective strategy for identifying the glycosylated biomarkers, we performed transcriptome analyses using quantitative real-time polymerase chain reaction (qRT-PCR) arrays and RNA sequencing (RNA-Seq). First, we measured and analyzed the transcriptome from the primary culture of human liver cells and hepatocarcinoma cells using RNA-Seq. This analysis revealed similar but distinctive expression profiles of glycogenes among hepatic cells as indicated by the qRT-PCR arrays, which determined a copy number of 186 glycogenes. Both data sets indicated that altered expression of glycosyltransferases affect the glycosylation of particular glycoproteins, which is consistent with the mass analysis data. Moreover, RNA-Seq analysis can uncover mutations in glycogenes and search differently expressed genes out of more than 50,000 distinct human gene transcripts including candidate biomarkers that were previously reported for hepatocarcinoma cells. Identification of candidate glyco-biomarkers from the expression profile of the glycogenes and proteins from liver cancer tissues available from public database emphasized the possibility that even though the expression level of biomarkers might not be altered, the expression of the glycogenes modifying biomarkers, generating glyco-biomarkers, might be different. Pathway analysis revealed that ~20% of the glycogenes exhibited different expression levels in normal and cancer cells. Thus, transcriptome analyses using both qRT-PCR array and RNA-Seq in combination with glycome and glycoproteome analyses can be advantageous to identify “glyco-biomarkers” by reinforcing information at the expression levels of both glycogenes and proteins.

Transcriptome sequencing identified hub genes for hepatocellular carcinoma by weighted-gene co-expression analysis

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide, and it remains a challenge to understand the genetic mechanisms underlying hepatocarcinogenesis. A global gene network of differential expression profiles in HCC has yet to be fully characterized. In the present study, we performed transcriptome sequencing (mRNA and lncRNA) in liver cancer and cirrhotic tissues of nine HCC patients. We identified differentially expressed genes (DEGs) and constructed a weighted gene co-expression network for the DEGs. In total, 755 DEGs (747 mRNA and eight lncRNA) were identified, and several co-expression modules were significantly associated with HCC clinical traits, including tumor location, tumor grade, and the α-fetoprotein (AFP) level. Of note, we identified 15 hub genes in the module associated with AFP level, and three (SPX, AFP and ADGRE1) of four hub genes were validated in an independent HCC cohort (n=78). Identification of hub genes for HCC clinical traits has implications for further understanding of the molecular genetic basis of HCC.

Combretastatin A-1 phosphate, a microtubule inhibitor, acts on both hepatocellular carcinoma cells and tumor-associated macrophages by inhibiting the Wnt/β-catenin pathway

Combretastatin A-1 phosphate (CA1P) is a microtubule polymerization inhibitor that binds to the colchicine-binding site of tubulin. We demonstrated that CA1P has outstanding anti-cancer activity against hepatocellular carcinoma (HCC) in vitro and in vivo. As determined by fluorescence staining and western blots (WBs), CA1P induced reactive oxygen species (ROS) accumulation and apoptosis in HepG2 cells with a down-regulation of Mcl-1. Additional studies indicated that CA1P induced microtubule depolymerization-mediated AKT inactivation, which resulted in GSK-3β activation, Wnt/β-Catenin pathway inhibition, and Mcl-1 down-regulation. The induction of HepG2 cell apoptosis by CA1P was prevented by a GSK-3β-specific inhibitor. Furthermore, immunohistochemistry studies on hepatocellular carcinoma mouse models showed that CA1P had activity against tumor-associated macrophages (TAMs). CA1P induced TAM apoptosis in vitro through the same mechanism observed with HepG2 cells, and it eliminated TAMs in the tumor microenvironment (TME) in vivo. In TME, the expression of TGF-β and TNF-α was also altered. The adoptive transfer of macrophages partly rescued the growth of tumor inhibited by CA1P. These findings indicate that CA1P has great potential to impact both cancer cells and the microenvironment, and our results should accelerate the application of CA1P for HCC therapy in clinic.

Novel fusion transcripts in bladder cancer identified by RNA-seq

Urothelial carcinoma (UC) is the most common type of bladder cancer and is the second most frequently diagnosed genitourinary tumor. The identification of fusion genes in bladder cancer might provide new perspectives for its classification and significance. In this study, we present a thorough search on three UC samples for novel fusion transcripts in bladder cancer using high-throughput RNA sequencing. We used stringent requirements for 819 fusion candidates and nominated 10 candidate fusion transcripts. Among them four novel fusion genes SEPT9/CYHR, IGF1R/TTC23, SYT8/TNNI2 and CASZ1/DFFA were validated and characterized in 48 formalin-fixed paraffin-embedded (FFPE) specimens of bladder cancer. Chromosomal rearrangements of regions 17q25, 15q26.3 and 1p36.22 resulting in the fusion transcripts SEPT9/CYHR, IGF1R/TTC23 and CASZ1/DFFA, appeared to be rare or unique events because they were not detected in the 48 UC samples. In contrast, the SYT8/TNNI2 fusion transcript resulting from transcription-induced chimerism by read-through mechanisms was a rather common and tumor-specific event occurring in 37.5% (18/48) of the UC specimens. Further investigation of functional and clinical relevance of novel fusion genes remains to be elucidated to reveal their role in bladder carcinogenesis.