RNA-Seq analyses generate comprehensive transcriptomic landscape and reveal complex transcript patterns in hepatocellular carcinoma

RNA Sequencing Analysis of Patients with Chronic Hepatitis B Treated Using PEGylated Interferon

Purpose

Worldwide, chronic hepatitis B virus (CHB) infection is a public health concern, ultimately leading to liver cirrhosis and hepatocellular carcinoma. Currently, patients with CHB can be treated using polyethylene glycol (PEG)ylated interferon (PEG-IFN) antiviral therapy, which has both immune modulatory and antiviral properties. This study aimed to reveal the mechanism underlying the effect of PEG-IFN therapy, to rationally optimize this therapeutic option.

Patients and Methods

Ten patients with CHB who were positive for the hepatitis B virus e antigen (HBeAg) and were receiving PEG-IFN treatment were enrolled. Clinical and virological parameters were monitored during 48 weeks of treatment. In addition, peripheral blood mononuclear cells (PBMCs) were collected from the 10 patients at 0, 24, and 36 weeks. RNA sequencing technology was used to analyze the RNA expression profile in the PBMC samples.

Results

Following PEG-IFN treatment, we identified 217 differentially expressed genes (DEGs), most of which were upregulated. Gene ontology enrichment analysis of the DEGs revealed that they were enriched in 29 clusters, mainly associated with "antiviral defense", "innate immunity", "immunity", "defense response to virus", "response to virus", "type I interferon signaling pathway", "negative regulation of viral genome replication", "innate immune response", and "RNA-binding".

Conclusion

After PEG-IFN treatment, a certain mRNA expression profile was observed in patients with CHB, providing further mechanistic insights into the antiviral effect of this therapy.

Transcriptome profiling and analysis of patients with esophageal squamous cell carcinoma from Kazakhstan

Esophageal squamous cell carcinoma (ESCC) is the predominant subtype of esophageal cancer in Central Asia, often diagnosed at advanced stages. Understanding population-specific patterns of ESCC is crucial for tailored treatments. This study aimed to unravel ESCC's genetic basis in Kazakhstani patients and identify potential biomarkers for early diagnosis and targeted therapies. ESCC patients from Kazakhstan were studied. We analyzed histological subtypes and conducted in-depth transcriptome sequencing. Differential gene expression analysis was performed, and significantly dysregulated pathways were identified using KEGG pathway analysis (p-value < 0.05). Protein-protein interaction networks were constructed to elucidate key modules and their functions. Among Kazakhstani patients, ESCC with moderate dysplasia was the most prevalent subtype. We identified 42 significantly upregulated and two significantly downregulated KEGG pathways, highlighting molecular mechanisms driving ESCC pathogenesis. Immune-related pathways, such as viral protein interaction with cytokines, rheumatoid arthritis, and oxidative phosphorylation, were elevated, suggesting immune system involvement. Conversely, downregulated pathways were associated with extracellular matrix degradation, crucial in cancer invasion and metastasis. Protein-protein interaction network analysis revealed four distinct modules with specific functions, implicating pathways in esophageal cancer development. High-throughput transcriptome sequencing elucidated critical molecular pathways underlying esophageal carcinogenesis in Kazakhstani patients. Insights into dysregulated pathways offer potential for early diagnosis and precision treatment strategies for ESCC. Understanding population-specific patterns is essential for personalized approaches to ESCC management.

Prioritization of potential pharmacological targets for the development of anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway: the reconstruction and analysis of associative gene networks help

Hepatocellular carcinoma (HCC) is a common severe type of liver cancer characterized by an extremely aggressive course and low survival rates. It is known that disruptions in the regulation of apoptosis activation are some of the key features inherent in most cancer cells, which determines the pharmacological induction of apoptosis as an important strategy for cancer therapy. The computer design of chemical compounds capable of specifically regulating the external signaling pathway of apoptosis induction represents a promising approach for creating new effective ways of therapy for liver cancer and other oncological diseases. However, at present, most of the studies are devoted to pharmacological effects on the internal (mitochondrial) apoptosis pathway. In contrast, the external pathway induced via cell death receptors remains out of focus. Aberrant gene methylation, along with hepatitis C virus (HCV) infection, are important risk factors for the development of hepatocellular carcinoma. The reconstruction of gene networks describing the molecular mechanisms of interaction of aberrantly methylated genes with key participants of the extrinsic apoptosis pathway and their regulation by HCV proteins can provide important information when searching for pharmacological targets. In the present study, 13 criteria were proposed for prioritizing potential pharmacological targets for developing anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway. The criteria are based on indicators of the structural and functional organization of reconstructed gene networks of hepatocarcinoma, the extrinsic apoptosis pathway, and regulatory pathways of virus-extrinsic apoptosis pathway interaction and aberrant gene methylation-extrinsic apoptosis pathway interaction using ANDSystem. The list of the top 100 gene targets ranked according to the prioritization rating was statistically significantly (p-value = 0.0002) enriched for known pharmacological targets approved by the FDA, indicating the correctness of the prioritization method. Among the promising potential pharmacological targets, six highly ranked genes (JUN, IL10, STAT3, MYC, TLR4, and KHDRBS1) are likely to deserve close attention.

Screening differential circular RNA expression profiles and the potential role of hsa_circ_0085465 in liver cancer

Aims

This study aimed to screen the circular RNAs (circRNAs) that are differentially expressed between liver cancer and paired paracarcinoma tissues and then elucidate their role in cancer progression.

Materials and Methods

High-throughput sequencing of cancer and paired paracarcinoma tissues was followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the parental genes of the differentially expressed circRNAs, which were also verified via real-time quantitative polymerase chain reaction analysis of the tissues. In addition, the function of selected circRNAs was determined using 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4- sulfophenyl)-2H-tetrazolium (MTS) and transwell assays.

Results

Total 218 and 121 circRNAs were differentially upregulated and downregulated, respectively; these were mainly enriched with GO and KEGG terms related to biological functions. From five representatives of the differentially expressed circRNAs, we selected hsa_circ_0085465 for further analysis, discovering that its overexpression promoted the proliferation, migration, and invasion of 97 L cells.

Conclusion

Taken together, our results suggest that hsa_circ_0085465 is relevant to liver cancer progression.

Transcriptome analysis creates a new era of precision medicine for managing recurrent hepatocellular carcinoma

The high incidence of hepatocellular carcinoma (HCC) recurrence negatively impacts outcomes of patients treated with curative intent despite advances in surgical techniques and other locoregional liver-targeting therapies. Over the past few decades, the emergence of transcriptome analysis tools, including real-time quantitative reverse transcription PCR, microarrays, and RNA sequencing, has not only largely contributed to our knowledge about the pathogenesis of recurrent HCC but also led to the development of outcome prediction models based on differentially expressed gene signatures. In recent years, the single-cell RNA sequencing technique has revolutionized our ability to study the complicated crosstalk between cancer cells and the immune environment, which may benefit further investigations on the role of different immune cells in HCC recurrence and the identification of potential therapeutic targets. In the present article, we summarized the major findings yielded with these transcriptome methods within the framework of a causal model consisting of three domains: primary cancer cells; carcinogenic stimuli; and tumor microenvironment. We provided a comprehensive review of the insights that transcriptome analyses have provided into diagnostics, surveillance, and treatment of HCC recurrence.

Comprehensive analysis of a TPX2-related TRHDE-AS1/PKIA ceRNA network involving prognostic signatures in Hepatitis B virus-infected hepatocellular carcinoma

Hepatitis B virus (HBV) infection is a main carcinogenic factor of hepatocellular carcinoma (HCC). TPX2 microtubule nucleation factor is recently recommended as a novel prognostic biomarker in HBV-infected HCC tissues. This study aimed to explore a TPX2-related ceRNA regulatory network in HBV-infected HCC and the potential impact on HCC prognosis. We comprehensively identified 541 differential expressed lncRNAs (DElncRNAs), 37 DEmiRNAs and 439 DEmRNAs from HBV-related TCGA-HCC cohorts in TPX2low and TPX2high groups. Based on their RNA-RNA interaction and expression analysis, four DElncRNAs (TRHDE-AS1, DLX6-AS1, SNHG14, HOXA11-AS), four DEmiRNAs (miR-23b, miR-320a, miR-589, miR-126) and five DEmRNAs (PKIA, PCDHA2, SHCBP1, PRSS16, KIF18A) in HCC tumor vs normal groups were subjected to the hub regulatory networks analysis and further prognostic value analysis. Importantly, the TRHDE-AS1/miR-23b/PKIA ceRNA network was associated with HCC prognosis. Furthermore, cellular location analysis and base-base interaction analysis indicated that the cytoplasmic lncRNA TRHDE-AS1 was regarded as a ceRNA to sponging miR-23b and then regulating PKIA. Interestingly, correlation analysis suggested the expression correlation between TRHDE-AS1 and PKIA in HCC. Finally, we further performed the methylation and immune infiltration analysis to explore the functional process of PKIA in HCC. We proposed a ceRNA regulatory network may help elucidate the mechanism by which TPX2 contributes to the prognosis of HBV-related HCC.

RNA sequencing analysis of hepatocellular carcinoma identified oxidative phosphorylation as a major pathologic feature

Dysregulation of expression of functional genes and pathways plays critical roles in the etiology and progression of hepatocellular carcinoma (HCC). Next generation-based RNA sequencing (RNA-seq) offers unparalleled power to comprehensively characterize HCC at the whole transcriptome level. In this study, 17 fresh-frozen HCC samples with paired non-neoplastic liver tissue from Caucasian patients undergoing liver resection or transplantation were used for RNA-seq analysis. Pairwise differential expression analysis of the RNA-seq data was performed to identify genes, pathways, and functional terms differentially regulated in HCC versus normal tissues. At a false discovery rate (FDR) of 0.10, 13% (n = 4335) of transcripts were up-regulated and 19% (n = 6454) of transcripts were down-regulated in HCC versus non-neoplastic tissue. Eighty-five Kyoto Encyclopedia of Genes and Genomes pathways were differentially regulated (FDR, <0.10), with almost all pathways (n = 83) being up-regulated in HCC versus non-neoplastic tissue. Among the top up-regulated pathways was oxidative phosphorylation (hsa00190; FDR, 1.12E-15), which was confirmed by Database for Annotation, Visualization, and Integrated Discovery (DAVID) gene set enrichment analysis. Consistent with potential oxidative stress due to activated oxidative phosphorylation, DNA damage-related signals (e.g., the up-regulated hsa03420 nucleotide excision repair [FDR, 1.14E-04] and hsa03410 base excision repair [FDR, 2.71E-04] pathways) were observed. Among down-regulated genes (FDR, <0.10), functional terms related to cellular structures (e.g., cell membrane [FDR, 3.05E-21] and cell junction [FDR, 2.41E-07], were highly enriched, suggesting compromised formation of cellular structure in HCC at the transcriptome level. Interestingly, the olfactory transduction (hsa04740; FDR, 1.53E-07) pathway was observed to be down-regulated in HCC versus non-neoplastic tissue, suggesting impaired liver chemosensory functions in HCC. Our findings suggest oxidative phosphorylation and the associated DNA damage may be the major driving pathologic feature in HCC.

Alternative splicing and liver disease

Alternative splicing produces complex and dynamic changes in the protein isoforms that are necessary for the proper biological functioning of the metabolic pathways involved in liver development and hepatocyte homeostasis. Changes in the physiological state of alternatively spliced forms are increasingly linked to liver pathologies. This may occur when the expression or function of the set of proteins controlling the alternative splicing processes are altered by external effectors such as oxidative stress and other environmental variations. Studies addressing these modifications reveal a complex interplay between the expression levels of different proteins that regulate the alternative splicing process as well as the changes in alternative splicing. This interplay results in a cascade of different protein isoforms that correlate with the progression of non-alcoholic fatty liver disease, hepatocellular carcinoma, and alcoholic liver disease. However, research on the detailed molecular mechanism underlying the production of these isoforms is needed. It is imperative to identify the physiological processes affected by the differentially spliced isoforms and confirm their role on the onset and maintenance of the pathology. This is required to design potential therapeutic approaches targeting the key splicing changes to revert the pathological condition as well as identify prognostic markers. In this review, we describe the complexity of the splicing process through an example to encourage researchers to go down this path. Subsequently, rather than a catalog of splicing events we have hand-picked and discuss a few selected studies of specific liver pathologies and suggested ways to focus research on these areas.

Advances in Multi-Omics Applications in HBV-Associated Hepatocellular Carcinoma

Hepatitis B virus (HBV) specifically infects liver cells, leading to progressive liver cirrhosis and significantly increasing the risk of hepatocellular carcinoma (HCC). The maturity of sequencing technology, improvement in bioinformatics data analysis and progress of omics technologies had improved research efficiency. The occurrence and progression of HCC are affected by multisystem and multilevel pathological changes. With the application of single-omics technologies, including genomics, transcriptomics, metabolomics and proteomics in tissue and body fluid samples, and even the novel development of multi-omics analysis on a single-cell platform, HBV-associated HCC changes can be better analyzed. The review summarizes the application of single omics and combined analysis of multi-omics data in HBV-associated HCC and proposes the importance of multi-omics analysis in the type of HCC, which provide the possibility for the precise diagnosis and therapy of HBV-associated HCC.

Quantitative analysis of the splice variants expressed by the major hepatitis B virus genotypes

Hepatitis B virus (HBV) is a major human pathogen that causes liver diseases. The main HBV RNAs are unspliced transcripts that encode the key viral proteins. Recent studies have shown that some of the HBV spliced transcript isoforms are predictive of liver cancer, yet the roles of these spliced transcripts remain elusive. Furthermore, there are nine major HBV genotypes common in different regions of the world, these genotypes may express different spliced transcript isoforms. To systematically study the HBV splice variants, we transfected human hepatoma cells, Huh7, with four HBV genotypes (A2, B2, C2 and D3), followed by deep RNA-sequencing. We found that 13-28 % of HBV RNAs were splice variants, which were reproducibly detected across independent biological replicates. These comprised 6 novel and 10 previously identified splice variants. In particular, a novel, singly spliced transcript was detected in genotypes A2 and D3 at high levels. The biological relevance of these splice variants was supported by their identification in HBV-positive liver biopsy and serum samples, and in HBV-infected primary human hepatocytes. Interestingly the levels of HBV splice variants varied across the genotypes, but the spliced pregenomic RNA SP1 and SP9 were the two most abundant splice variants. Counterintuitively, these singly spliced SP1 and SP9 variants had a suboptimal 5' splice site, supporting the idea that splicing of HBV RNAs is tightly controlled by the viral post-transcriptional regulatory RNA element.

Screening and Identification of Potential Biomarkers in Hepatitis B Virus-Related Hepatocellular Carcinoma by Bioinformatics Analysis

Hepatocellular carcinoma (HCC) is one of the most lethal cancers globally. Hepatitis B virus (HBV) infection might cause chronic hepatitis and cirrhosis, leading to HCC. To screen prognostic genes and therapeutic targets for HCC by bioinformatics analysis and determine the mechanisms underlying HBV-related HCC, three high-throughput RNA-seq based raw datasets, namely GSE25599, GSE77509, and GSE94660, were obtained from the Gene Expression Omnibus database, and one RNA-seq raw dataset was acquired from The Cancer Genome Atlas (TCGA). Overall, 103 genes were up-regulated and 127 were down-regulated. A protein–protein interaction (PPI) network was established using Cytoscape software, and 12 pivotal genes were selected as hub genes. The 230 differentially expressed genes and 12 hub genes were subjected to functional and pathway enrichment analyses, and the results suggested that cell cycle, nuclear division, mitotic nuclear division, oocyte meiosis, retinol metabolism, and p53 signaling-related pathways play important roles in HBV-related HCC progression. Further, among the 12 hub genes, kinesin family member 11 (KIF11), TPX2 microtubule nucleation factor (TPX2), kinesin family member 20A (KIF20A), and cyclin B2 (CCNB2) were identified as independent prognostic genes by survival analysis and univariate and multivariate Cox regression analysis. These four genes showed higher expression levels in HCC than in normal tissue samples, as identified upon analyses with Oncomine. In addition, in comparison with normal tissues, the expression levels of KIF11, TPX2, KIF20A, and CCNB2 were higher in HBV-related HCC than in HCV-related HCC tissues. In conclusion, our results suggest that KIF11, TPX2, KIF20A, and CCNB2 might be involved in the carcinogenesis and development of HBV-related HCC. They can thus be used as independent prognostic genes and novel biomarkers for the diagnosis of HBV-related HCC and development of pertinent therapeutic strategies.

Analysis of Differentially Expressed Genes in a Chinese Cohort of Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is a leading malignancy in China with both high incidence and mortality. Towards improving outcomes, clinically-relevant biomarkers are urgently needed for use as prognostic and treatment targets. Herein we applied RNA-seq for deep sequencing of ten matched pairs of ESCC and adjacent non-cancerous tissues (NT) from Chinese patients. Transcriptomic data mapped to approximately 64% of all annotated genes with 2,047 and 708 unigenes being differentially up-regulated and down-regulated, respectively, between ESCCs and NT samples (p<0.05). Dividing cases by pathological grade revealed significant differentially expressed genes (DEG) between ESCC and NT in both low and high differentiation cases (p<0.05) whereas gene expression differences were not significantly different between high and low differentiation ESCC tissues (p=0.053). Moreover, the majority of ESCC and NT tissues formed clusters in principal component analyses. The veracity of the DEG list was validated in a larger cohort of 45 patient samples, with down-regulated CLIC3, up-regulated CLIC4 and unchanged expression of CLIC2 confirmed in ESCC using quantitative PCR and Western blotting. Our data reveal both previously identified ESCC biomarkers along with novel candidates and represent a ready resource of DEGs in ESCC for further investigation.

Tumor-infiltrating immune cells in hepatocellular carcinoma: Tregs is correlated with poor overall survival

Systematic interrogation of tumor-infiltrating immune cells (TIICs) is key to the prediction of clinical outcome and development of immunotherapies. However, little is known about the TIICs of hepatocellular carcinoma (HCC) and its impact on the prognosis of patients and potential for immunotherapy. We applied CIBERSORT of 1090 tumors to infer the infiltration of 22 subsets of TIICs using gene expression data. Unsupervised clustering analysis by 22 TIICs revealed 4 clusters of tumors, mainly defined by macrophages and T cells, with distinct prognosis and associations with immune checkpoint molecules, including PD-1, CD274, CTLA-4, LAG-3 and IFNG. We found tumors with decreased number of M1 macrophages or increased regulatory T cells were associated with poor prognosis. Based on the multivariate Cox analysis, a nomogram was also established for clinical application. In conclusion, composition of the TIICs in HCC was quite different, which is an important determinant of prognosis with great potential to identify candidates for immunotherapy.

AllEnricher: a comprehensive gene set function enrichment tool for both model and non-model species

BACKGROUND

Function genomic studies will generally result in lists of genes that may provide clues for exploring biological questions and discovering unanticipated functions, based on differential gene expression analysis, differential epigenomic analysis or co-expression network analysis. While tools have been developed to identify biological functions that are enriched in the genes sets, there remains a need for comprehensive tools that identify functional enrichment of genes for both model and non-model species from a different function classification perspective.

RESULTS

We developed AllEnricher, a tool that calculates gene set function enrichment, with user-defined updatable libraries backing up for both model and non-model species as well as providing comprehensive functional interpretation from multiple dimensions, including GO, KEGG, Reactome, DO and DisGeNET.

CONCLUSIONS

AllEnricher incorporates up to date information from different public resources and provides a comprehensive resolution for biologists to make sense out of specific gene sets, making it an advanced open-source tool for gene set function analysis.

Multiple novel hepatocellular carcinoma signature genes are commonly controlled by the master pluripotency factor OCT4

BACKGROUND

Worldwide, hepatocellular carcinoma (HCC) is a common solid tumor with a poor prognosis. HCC is often due to hepatitis B virus (HBV) infection. As yet, efficacious HCC treatment regimens for late-stage HCC patients are lacking. Therefore, the identification of more specific and sensitive biomarkers for its early diagnosis and treatment remains an urgent need.

METHODS

Total RNAs from paired HBV-derived HCC tumors and adjacent peritumor tissues (APTs) were subjected to RNA sequencing (RNA-seq), and differentially expressed genes (DEGs) between HCC tumors and APTs were selected and verified.

RESULTS

We identified 166 DEGs and found that eight top-ranked and verified DEGs (TK1, CTTN, CEP72, TRIP13, FTH1, FLAD1, CHRM2, AMBP) all contained putative OCT4 binding motifs in their promoter regions. TK1, TRIP13 and OCT4 were found to exhibit concurrent higher expression levels in HCC tumors than in APTs. The mRNA levels of TK1, TRIP13 and OCT4 in a cohort of 384 HCC samples from the TCGA database were all found to be negatively correlated with patient overall survival, relapse-free survival and progression-free survival, underscoring the HCC biomarker status of TK1 and TRIP13 on one hand, and implicating their association with OCT4 on the other hand. Furthermore, OCT4 proteins were found to bind to the promoters of both genes in vitro and in vivo. Knocking out OCT4 in HCC-derived cell lines reduced the expression of TK1 and TRIP13 and significantly decreased their tumorigenicity.

CONCLUSIONS

Using RNA-seq, we identified several novel HCC signature genes that may serve as biomarkers for its diagnosis and prognosis. Their common transcriptional regulation by OCT4 suggests key roles in the development of HCC, and indicates that OCT4 may serve as a potential therapeutic target.

Inflammation associated with noise-induced hearing loss

Inflammation is a complex biological response to harmful stimuli including infection, tissue damage, and toxins. Thus, it is not surprising that cochlear damage by noise includes an inflammatory component. One mechanism by which inflammation is generated by tissue damage is the activation of damage-associated molecular patterns (DAMPs). Many of the cellular receptors for DAMPS, including Toll-like receptors, NOD-like receptors, and DNA receptors, are also receptors for pathogens, and function in the innate immune system. DAMP receptors are known to be expressed by cochlear cells, and binding of molecules released by damaged cells to these receptors result in the activation of cell stress pathways. This leads to the generation of pro-inflammatory cytokines and chemokines that recruit pro-inflammatory leukocytes. Extensive evidence indicates pro-inflammatory cytokines including TNF alpha and interleukin 1 beta, and chemokines including CCL2, are induced in the cochlea after noise exposure. The recruitment of macrophages into the cochlea has also been demonstrated. These provide substrates for noise damage to be enhanced by inflammation. Evidence is provided by the effectiveness of anti-inflammatory drugs in ameliorating noise-induced hearing loss. Involvement of inflammation provides a wide variety of additional anti-inflammatory and pro-resolution agents as potential pharmacological interventions in noise-induced hearing loss.

Detection of Novel Fusion Transcript VTI1A-CFAP46 in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is now the second-highest cause of cancer death worldwide. Recent studies have discovered a wide range of somatic mutations in HCC. These mutations involve various vital signaling pathways such as: Wnt/β-Catenin, p53, telome-rase reverse transcriptase (TERT), chromatin remodeling, RAS/MAPK signaling, and oxidative stress. However, fusion transcripts have not been broadly explored in HCC.

Methods

To identify novel fusion transcripts in HCC, in the first phase of our study, we performed targeted RNA sequencing (in HCC and paired non-HCC tissues) on 6 patients with a diagnosis of HCC undergoing liver transplantation.

Results

As a result of these studies, we discovered the novel fusion transcript, VTI1A-CFAP46. In the second phase of our study, we measured the expression of wild-type VTI1A in 21 HCC specimens, which showed that 10 of 21 exhibited upregulation of wild-type VTI1A in their tumors. VTI1A (Vesicle Transport via Interaction with t-SNARE homolog 1A) is a member of the Soluble N-ethylmaleimide-Sensitive Factor (NSF) attachment protein receptor (SNARE) gene family, which is essential for membrane trafficking and function in endocytosis, autophagy, and Golgi transport. Notably, it is known that autophagy is involved in HCC.

Conclusions

The link between novel fusion transcript VTI1A-CFAP46 and autophagy as a potential therapeutic target in HCC patients deserves further investigation. Moreover, this study shows that fusion transcripts are worthy of additional exploration in HCC.

Comprehensive analysis of transcriptome profiles in hepatocellular carcinoma

Background

Hepatocellular carcinoma is the second most deadly cancer with late presentation and limited treatment options, highlighting an urgent need to better understand HCC to facilitate the identification of early-stage biomarkers and uncover therapeutic targets for the development of novel therapies for HCC.

Methods

Deep transcriptome sequencing of tumor and paired non-tumor liver tissues was performed to comprehensively evaluate the profiles of both the host and HBV transcripts in HCC patients. Differential gene expression patterns and the dys-regulated genes associated with clinical outcomes were analyzed. Somatic mutations were identified from the sequencing data and the deleterious mutations were predicted. Lastly, human-HBV chimeric transcripts were identified, and their distribution, potential function and expression association were analyzed.

Results

Expression profiling identified the significantly upregulated TP73 as a nodal molecule modulating expression of apoptotic genes. Approximately 2.5% of dysregulated genes significantly correlated with HCC clinical characteristics. Of the 110 identified genes, those involved in post-translational modification, cell division and/or transcriptional regulation were upregulated, while those involved in redox reactions were downregulated in tumors of patients with poor prognosis. Mutation signature analysis identified that somatic mutations in HCC tumors were mainly non-synonymous, frequently affecting genes in the micro-environment and cancer pathways. Recurrent mutations occur mainly in ribosomal genes. The most frequently mutated genes were generally associated with a poorer clinical prognosis. Lastly, transcriptome sequencing suggest that HBV replication in the tumors of HCC patients is rare. HBV-human fusion transcripts are a common observation, with favored HBV and host insertion sites being the HBx C-terminus and gene introns (in tumors) and introns/intergenic-regions (in non-tumors), respectively. HBV-fused genes in tumors were mainly involved in RNA binding while those in non-tumors tissues varied widely. These observations suggest that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts may occur during tumorigenesis.

Conclusions

Transcriptome sequencing of HCC patients reveals key cancer molecules and clinically relevant pathways deregulated/mutated in HCC patients and suggests that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts likely occur during the process of tumorigenesis.

Ancestry patterns inferred from massive RNA-seq data

There is a growing body of evidence suggesting that patterns of gene expression vary within and between human populations. However, the impact of this variation in human diseases has been poorly explored, in part owing to the lack of a standardized protocol to estimate biogeographical ancestry from gene expression studies. Here we examine several studies that provide new solid evidence indicating that the ancestral background of individuals impacts gene expression patterns. Next, we test a procedure to infer genetic ancestry from RNA-seq data in 25 data sets where information on ethnicity was reported. Genome data of reference continental populations retrieved from The 1000 Genomes Project were used for comparisons. Remarkably, only eight out of 25 data sets passed FastQC default filters. We demonstrate that, for these eight population sets, the ancestral background of donors could be inferred very efficiently, even in data sets including samples with complex patterns of admixture (e.g., American-admixed populations). For most of the gene expression data sets of suboptimal quality, ancestral inference yielded odd patterns. The present study thus brings a cautionary note for gene expression studies highlighting the importance to control for the potential confounding effect of ancestral genetic background.

Dissecting the expression landscape of cytochromes P450 in hepatocellular carcinoma: towards novel molecular biomarkers

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related deaths around the world. Recent advances in genomic technologies have allowed the identification of various molecular signatures in HCC tissues. For instance, differential gene expression levels of various cytochrome P450 genes (CYP450) have been reported in studies performed on limited numbers of HCC tissue samples, or focused on a small subset on CYP450s. In the present study, we monitored the expression landscape of all the members of the CYP450 family (57 genes) in more than 200 HCC tissues using RNA-Seq data from The Cancer Genome Atlas. Using stringent statistical filters and data from paired tissues, we identified significantly dysregulated CYP450 genes in HCC. Moreover, the expression level of selected CYP450s was validated by qPCR on cDNA samples from an independent cohort. Threshold values (sensitivity and specificity) based on dysregulated gene expression were also determined to allow for confident identification of HCC tissues. Finally, a global look at expression levels of the 57 members of the CYP450 family across ten different cancer types revealed specific expression signatures. Overall, this study provides useful information on the transcriptomic landscape of CYP450 genes in HCC and on new potential HCC biomarkers.

Observational Study of PD-L1, TGF-β, and Immune Cell Infiltrates in Hepatocellular Carcinoma

Introduction: Hepatocellular carcinoma (HCC) typically develops in cirrhotic livers, with increased programed death ligand 1 (PD-L1) and transforming growth factor beta (TGF-β) activity implicated in immunosuppression.

Methods: In an observational study of HCC liver samples, we determined the incidence of PD-L1 and immune cell (IC) infiltrates, and signs of TGF-β activity. HCCs were characterized by the incidence and distribution of PD-L1⁺ cells, and CD8⁺, CD68⁺, and FoxP3⁺ infiltrating ICs in HCC and surrounding liver. Gene expression signatures (GESs) associated with TGF-β activity and ICs were evaluated by RNAseq.

Results: In non-neoplastic cirrhotic and non-cirrhotic liver, PD-L1 occurred on sinusoidal lining cells (mostly Kupffer cells), endothelial cells and ICs. In HCC, PD-L1⁺ tumor cells were rare. Most PD-L1⁺ cells were identified as ICs. CD8⁺, CD68⁺, and FoxP3⁺ ICs were associated with HCC, particularly in the invasive margin. CD8⁺ cell incidence correlated with PD-L1⁺ cells, consistent with PD-L1 being upregulated in response to pre-existing cytotoxic T-lymphocyte activity. TGFB1 mRNA levels and TGF-β activation GES correlated with the strength of the tumor-associated macrophage GES.

Conclusion: Inhibition of PD-L1⁺ ICs and TGF-β activity and their respective immunomodulatory pathways may contribute to antitumor effects in HCC.

Profiling of LINE-1-Related Genes in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a prime public health concern that accounts for most of the primary liver malignancies in humans. The most common etiological factor of HCC is hepatitis B virus (HBV). Despite recent advances in treatment strategies, there has been little success in improving the survival of HCC patients. To develop a novel therapeutic approach, evaluation of a working hypothesis based on different viewpoints might be important. Long interspersed element 1 (L1) retrotransposons have been suggested to play a role in HCC. However, the molecular machineries that can modulate L1 biology in HBV-related HCC have not been well-evaluated. Here, we summarize the profiles of expression and/or activation status of L1-related genes in HBV-related HCC, and HBV- and HCC-related genes that may impact L1-mediated tumorigenesis. L1 restriction factors appear to be suppressed by HBV infection. Since some of the L1 restriction factors also limit HBV, these factors may be exhausted in HBV-infected cells, which causes de-suppression of L1. Several HBV- and HCC-related genes that interact with L1 can affect oncogenic processes. Thus, L1 may be a novel prime therapeutic target for HBV-related HCC. Studies in this area will provide insights into HCC and other types of cancers.

Integrated analysis of lncRNA-associated ceRNA network reveals potential biomarkers for the prognosis of hepatitis B virus-related hepatocellular carcinoma

Background

There is evidence that abnormal expression of lncRNAs is associated with hepatitis B virus (HBV) infection-induced hepatocellular carcinoma (HCC). However, the mechanisms remain not fully elucidated. The study aimed to identify novel lncRNAs and explore their underlying mechanisms based on the ceRNA hypothesis.

Methods

The RNA and miRNA expression profiling in 20 tumor and matched adjacent tissues from HBV–HCC patients were retrieved from the Gene Expression Omnibus database under accession numbers GSE77509 and GSE76903, respectively. Differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs) were identified using the EdgeR package. Protein–protein interaction (PPI) network was constructed for DEGs followed by module analysis. The ceRNA network was constructed based on interaction relationships between miRNAs and mRNAs/lncRNAs. The functions of DEGs were predicted using DAVID and BinGO databases. The prognosis values (overall survival [OS] and recurrence-free survival [RFS]) of ceRNA network genes were determined using The Cancer Genome Atlas (TCGA) data with Cox regression analysis and Kaplan–Meier method.

Results

The present study screened 643 DELs, 83 DEMs, and 1,187 DEGs. PPI network analysis demonstrated that CDK1 and CCNE1 were hub genes and extracted in functionally related modules. E2F2, CDK1, and CCNE1 were significantly enriched into cell cycle pathway. FAM182B-miR-125b-5p-E2F2 and LINC00346-miR-10a-5p-CDK1/CCNE1 ceRNA axes were obtained by constructing the ceRNA network. Patients with high expressions of DELs and DEGs in the above ceRNA axes had poor OS, while patients with the high expression of DEMs possessed excellent OS. CDK1 was also an RFS-related biomarker, with its high expression predicting poor RFS. The upregulation of LINC00346 and CDK1 but the downregulation of miR-10a-5p in HCC was validated in other microarray datasets and TCGA database.

Conclusion

The LINC00346-miR-10a-5p-CDK1 axis may be an important mechanism for HBV-related HCC, and genes in this ceRNA axis may be potential prognostic biomarkers and therapeutic targets.

Transcriptome-Wide Analysis Reveals the Landscape of Aberrant Alternative Splicing Events in Liver Cancer

Alternative splicing (AS) is assumed to be a pivotal determinant for the generation of diverse transcriptional variants in cancer. However, the comprehensive dysregulation of AS and the prospective biological and clinical relevance in hepatocellular carcinoma (HCC) remain obscure. Here, we identified and depicted the AS landscape in HCC by performing reference-based assembly of sequencing reads from over 600 RNA sequencing (RNA-seq) libraries. We detected various differentially spliced ASEs across patients covering not only protein-coding genes, but also considerable numbers of noncoding genes. Strikingly, alternative transcription initiation was found to frequently occur in HCC. These differential ASEs were highly related to "cancer hallmarks" and involved in metabolism-related pathways in particular. In addition, 243 differential ASEs were identified as risk predictors for HCC patient survival. The isoform switch of metabolism-related gene UGP2 (UDP-glucose pyrophosphorylase 2) might play an essential role in HCC. We further constructed regulatory networks between RNA-binding protein (RBP) genes and the corresponding ASEs. Further analysis demonstrated that the regulated networks were enriched in a variety of metabolism-related pathways. Conclusion: Differential ASEs are prevalent in HCC, where alternative transcription initiation was found to frequently occur. We found that genes having differential ASEs were significantly enriched in metabolism-related pathways. The expression variations, binding relations, and even mutations of RBP genes largely influenced differential ASEs in HCC.

Disulfide bridge formation influences ligand recognition by the ATAD2 bromodomain

The ATPase family, AAA domain-containing protein 2 (ATAD2) has a C-terminal bromodomain, which functions as a chromatin reader domain recognizing acetylated lysine on the histone tails within the nucleosome. ATAD2 is overexpressed in many cancers and its expression is correlated with poor patient outcomes, making it an attractive therapeutic target and potential biomarker. We solved the crystal structure of the ATAD2 bromodomain and found that it contains a disulfide bridge near the base of the acetyllysine binding pocket (Cys1057-Cys1079). Site-directed mutagenesis revealed that removal of a free C-terminal cysteine (C1101) residue greatly improved the solubility of the ATAD2 bromodomain in vitro. Isothermal titration calorimetry experiments in combination with the Ellman's assay demonstrated that formation of an intramolecular disulfide bridge negatively impacts the ligand binding affinities and alters the thermodynamic parameters of the ATAD2 bromodomain interaction with a histone H4K5ac peptide as well as a small molecule bromodomain ligand. Molecular dynamics simulations indicate that the formation of the disulfide bridge in the ATAD2 bromodomain does not alter the structure of the folded state or flexibility of the acetyllysine binding pocket. However, consideration of this unique structural feature should be taken into account when examining ligand-binding affinity, or in the design of new bromodomain inhibitor compounds that interact with this acetyllysine reader module.

Identification of key genes and pathways by bioinformatics analysis with TCGA RNA sequencing data in hepatocellular carcinoma

Improved insight into the molecular characteristics of hepatocellular carcinoma (HCC) is required to predict prognosis and to develop a new rationale for targeted therapeutic strategy. Bioinformatics methods, including functional enrichment and network analysis combined with survival analysis, are required to process a large volume of data to obtain further information on differentially expressed genes (DEGs). The RNA sequencing data related to HCC in The Cancer Genome Atlas (TCGA) database were analyzed to screen DEGs, which were separately submitted to perform gene enrichment analysis to identify gene sets and signaling pathways, and to construct a protein-protein interaction (PPI) network. Subsequently, hub genes were selected by the core level in the network, and the top hub genes were focused on gene expression analysis and survival analysis. A total of 610 DEGs were identified, including 444 upregulated and 166 downregulated genes. The upregulated DEGs were significantly enriched in the Gene Ontology analysis (GO): Cell division and in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway: Cell cycle, whereas the downregulated DEGs were enriched in GO: Negative regulation of growth and in the KEGG pathway: Retinol metabolism, with significant differences. Cyclin-dependent kinase (CDK)1 was selected as the top hub gene by the PPI network, which exhibited a similar expression trend with the data from the Gene Expression Omnibus (GEO) database. Survival analysis revealed a significantly negative correlation between CDK1 expression level and overall survival in the TCGA group (P<0.01) and the GEO group (P<0.01). Therefore, high-throughput TCGA data analysis appears to be an effective method for screening tumor molecular markers, and high expression of CDK1 is a prognostic factor for HCC.

Epigenome-wide association study in hepatocellular carcinoma: Identification of stochastic epigenetic mutations through an innovative statistical approach

Hepatocellular carcinoma (HCC) results from accumulation of both genetic and epigenetic alterations. We investigated the genome-wide DNA methylation profile in 69 pairs of HCC and adjacent non-cancerous liver tissues using the Infinium HumanMethylation 450K BeadChip array. An innovative analytical approach has been adopted to identify Stochastic Epigenetic Mutations (SEMs) in HCC.HCC and peritumoral tissues showed a different epigenetic profile, mainly characterized by loss of DNA methylation in HCC. Total number of SEMs was significantly higher in HCC tumor (median: 77,370) than in peritumoral (median: 5,656) tissues and correlated with tumor grade. A significant positive association emerged between SEMs measured in peritumoral tissue and hepatitis B and/or C virus infection status. A restricted number of SEMs resulted to be shared by more than 90% of HCC tumor samples and never present in peritumoral tissue. This analysis allowed the identification of four epigenetically regulated candidate genes (AJAP1, ADARB2, PTPRN2, SDK1), potentially involved in the pathogenesis of HCC.In conclusion, HCC showed a methylation profile completely deregulated and very far from adjacent non-cancerous liver tissues. The SEM analysis provided valuable clues for further investigations in understanding the process of tumorigenesis in HCC.

ATAD2 in cancer: a pharmacologically challenging but tractable target

INTRODUCTION

ATAD2 protein is an emerging oncogene that has strongly been linked to the etiology of multiple advanced human cancers. Therapeutically, despite the fact that genetic suppression/knockdown studies have validated it as a compelling drug target for future therapeutic development, recent druggability assessment data suggest that direct targeting of ATAD2's bromodomain (BRD) may be a very challenging task. ATAD2's BRD has been predicted as a 'difficult to drug' or 'least druggable' target due to the concern that its binding pocket, and the areas around it, seem to be unfeasible for ligand binding. Areas covered: In this review, after shedding light on the multifaceted roles of ATAD2 in normal physiology as well as in cancer-etiology, we discuss technical challenges rendered by ATAD2's BRD active site and the recent drug discovery efforts to find small molecule inhibitors against it. Expert opinion: The identification of a novel low-nanomolar semi-permeable chemical probe against ATAD2's BRD by recent drug discovery campaign has demonstrated it to be a pharmacologically tractable target. Nevertheless, the development of high quality bioavailable inhibitors against ATAD2 is still a pending task. Moreover, ATAD2 may also potentially be utilized as a promising target for future development of RNAi-based therapy to treat cancers.

Cholesterol overload in the liver aggravates oxidative stress-mediated DNA damage and accelerates hepatocarcinogenesis

Primary liver cancers represent the second leading cause of cancer-related deaths worldwide. Diverse etiological factors include chronic viral hepatitis, aflatoxin and alcohol exposure as well as aberrant liver lipid overload. Cholesterol has been identified as a key inducer of metabolic impairment, oxidative stress and promoter of cellular dysfunction. The aim of this work was to address the oxidative stress-mediated DNA damage induced by cholesterol overload, and its role in the development of hepatocellular carcinoma.

C57BL/6 male mice were fed with a high cholesterol diet, followed by a single dose of N-diethylnitrosamine (DEN, 10 μg/g, ip). Reactive oxygen species generation, DNA oxidation, antioxidant and DNA repair proteins were analyzed at different time points. Diet-induced cholesterol overload caused enhanced oxidative DNA damage in the liver and was associated with a decrease in key DNA repair genes as early as 7 days. Interestingly, we found a cell survival response, induced by cholesterol, judged by a decrement in Bax to Bcl2 ratio. Importantly, N-acetyl-cysteine supplementation significantly prevented DNA oxidation damage. Furthermore, at 8 months after DEN administration, tumor growth was significantly enhanced in mice under cholesterol diet in comparison to control animals. Together, these results suggest that cholesterol overload exerts an oxidative stress-mediated effects and promotes the development of liver cancer.

Molecular mechanisms of pathogenesis in hepatocellular carcinoma revealed by RNA‑sequencing

The present study aimed to explore the underlying molecular mechanisms of hepatocellular carcinoma (HCC). RNA‑sequencing profiles GSM629264 and GSM629265, from the GSE25599 data set, were downloaded from the Gene Expression Omnibus database and processed by quality evaluation. GSM629264 and GSM629265 were from HCC and adjacent non‑cancerous tissues, respectively. TopHat software was used for alignment analysis, followed by the detection of novel splicing sites. In addition, the Cufflinks software package was used to analyze gene expressions, and the Cuffdiff program was used to screen for differently expressed genes (DEGs) and differentially expressed splicing variants. Gene ontology functional enrichment and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of DEGs were also performed. Transcription factors (TFs) and microRNAs (miRNAs) that regulate DEGs were identified, and a protein‑protein interaction (PPI) network was constructed. The hub node in the PPI network was obtained, and the TFs and miRNAs that regulated the hub node were further predicted. The quality of the sequencing data met the standards for analysis, and the clean reads were ~65%. Most sequencing reads mapped into coding sequence exons (CDS_exons), whereas other reads mapped into exon 3' untranslated regions (UTR_Exons), 5'UTR_Exons and Introns. Upregulated and downregulated DEGs between HCC and adjacent non‑cancerous tissues were screened. Genes of differentially expressed splicing variants were identified, including vesicle‑associated membrane protein 4, phosphatidylinositol glycan anchor biosynthesis class C, protein disulfide isomerase family A member 4 and growth arrest specific 5. Screened DEGs were enriched in the complement pathway. In the PPI network, ubiquitin C (UBC) was the hub node. UBC was predicted to be regulated by several TFs, including specificity protein 1 (SP1), FBJ murine osteosarcoma viral oncogene homolog (FOS), proto‑oncogene c‑JUN (JUN), FOS‑like antigen 2 (FOSL2) and SWI/SNF‑related, matrix‑associated, actin‑dependent regulator of chromatin, subfamily A, member 4 (SMARCA4), and several miRNAs, including miR‑30 and miR‑181. Results from the present study demonstrated that UBC, SP1, FOS, JUN, FOSL2, SMARCA4, miR‑30 and miR‑181 may participate in the development of HCC.

An RNA-based signature enables high specificity detection of circulating tumor cells in hepatocellular carcinoma

Circulating tumor cells (CTCs) are shed into the bloodstream by invasive cancers, but the difficulty inherent in identifying these rare cells by microscopy has precluded their routine use in monitoring or screening for cancer. We recently described a high-throughput microfluidic CTC-iChip, which efficiently depletes hematopoietic cells from blood specimens and enriches for CTCs with well-preserved RNA. Application of RNA-based digital PCR to detect CTC-derived signatures may thus enable highly accurate tissue lineage-based cancer detection in blood specimens. As proof of principle, we examined hepatocellular carcinoma (HCC), a cancer that is derived from liver cells bearing a unique gene expression profile. After identifying a digital signature of 10 liver-specific transcripts, we used a cross-validated logistic regression model to identify the presence of HCC-derived CTCs in nine of 16 (56%) untreated patients with HCC versus one of 31 (3%) patients with nonmalignant liver disease at risk for developing HCC (P < 0.0001). Positive CTC scores declined in treated patients: Nine of 32 (28%) patients receiving therapy and only one of 15 (7%) patients who had undergone curative-intent ablation, surgery, or liver transplantation were positive. RNA-based digital CTC scoring was not correlated with the standard HCC serum protein marker alpha fetoprotein (P = 0.57). Modeling the sequential use of these two orthogonal markers for liver cancer screening in patients with high-risk cirrhosis generates positive and negative predictive values of 80% and 86%, respectively. Thus, digital RNA quantitation constitutes a sensitive and specific CTC readout, enabling high-throughput clinical applications, such as noninvasive screening for HCC in populations where viral hepatitis and cirrhosis are prevalent.

NEK2 serves as a prognostic biomarker for hepatocellular carcinoma

Never in mitosis gene A (NIMA)-related kinase 2 (NEK2) is a microtubule-associated protein that regulates spindle assembly in human cells and is overexpressed in various malignancies. However, the role of NEK2 in hepatocellular carcinoma (HCC) remains undetermined. We performed RNA-seq of the HCC cell line SMMC-7721 and the normal liver cell line HL-7702 using the Ion Proton System. NEK2 expression was detected using quantitative reverse transcription polymerase chain reaction in two cell lines and 5 matched HCC and adjacent non-tumorous liver tissues. The correlation between survival and NEK2 expression was analyzed in 359 patients with HCC using RNASeqV2 data available from The Cancer Genome Atlas (TCGA) website (https://tcga-data.nci.nih.gov/tcga/). The expression of NEK2, phospho-AKT and MMP-2 was evaluated by immunohistochemistry in 63 cases of HCC and matched adjacent non-tumorous liver tissues. Relationships between protein expression and clinicopathological parameters were assessed, and the correlations between NEK2 with phospho-AKT and MMP-2 expressions were evaluated. A total of 610 differentially expressed genes (DEGs) were revealed in the transcriptome comparison, 297 of which were upregulated and 313 were downregulated in HCC. NEK2, as the most obviously different DEG in cells and tissues from the RNA-seq data, was listed as an HCC candidate biomarker for further verification. NEK2 was overexpressed in HCC cells and tissues (P=0.002, P=0.013) and HCC patients with a high expression of NEK2 had a poor prognosis (P=0.0145). Clinical analysis indicated that the overexpression of NEK2 in HCC was significantly correlated with diolame complete (P<0.001), tumor nodule number (P=0.012) and recurrence (P=0.004). NEK2 expression was positively correlated with the expression of phospho-AKT (r=0.883, P<0.01) and MMP-2 (r=0.781, P<0.01). Overexpression of NEK2 was associated with clinicopathological characteristics and poor patient outcomes, suggesting that NEK2 serves as a prognostic biomarker for HCC. Alteration of NEK2 protein levels may contribute to invasion and metastasis of HCC, which may occur through activation of AKT signaling and promotion of MMP-2 expression.

Alternative RNA Splicing in the Pathogenesis of Liver Disease

Non-alcoholic fatty liver disease (NAFLD) is becoming increasingly prevalent due to the worldwide obesity epidemic and currently affects one-third of adults or about one billion people worldwide. NAFLD is predicted to affect over 50% of the world's population by the end of the next decade. It is the most common form of liver disease and is associated with increased risk for progression to a more severe form non-alcoholic steatohepatitis, as well as insulin resistance, type 2 diabetes mellitus, cirrhosis, and eventually hepatocellular carcinoma. This review article will focus on the role of alternative splicing in normal liver physiology and dysregulation in liver disease.

Circadian Homeostasis of Liver Metabolism Suppresses Hepatocarcinogenesis

Chronic jet lag induces spontaneous hepatocellular carcinoma (HCC) in wild-type mice following a mechanism very similar to that observed in obese humans. The process initiates with non-alcoholic fatty liver disease (NAFLD) that progresses to steatohepatitis and fibrosis before HCC detection. This pathophysiological pathway is driven by jet-lag-induced genome-wide gene deregulation and global liver metabolic dysfunction, with nuclear receptor-controlled cholesterol/bile acid and xenobiotic metabolism among the top deregulated pathways. Ablation of farnesoid X receptor dramatically increases enterohepatic bile acid levels and jet-lag-induced HCC, while loss of constitutive androstane receptor (CAR), a well-known liver tumor promoter that mediates toxic bile acid signaling, inhibits NAFLD-induced hepatocarcinogenesis. Circadian disruption activates CAR by promoting cholestasis, peripheral clock disruption, and sympathetic dysfunction.

Molecular Pathogenesis of Hepatocellular Carcinoma

The pathogenesis of hepatocellular carcinoma (HCC) is a multistep process involving the progressive accumulation of molecular alterations pinpointing different molecular and cellular events. The next-generation sequencing technology is facilitating the global and systematic evaluation of molecular landscapes in HCC. There is emerging evidence supporting the importance of cancer metabolism and tumor microenvironment in providing a favorable and supportive niche to expedite HCC development. Moreover, recent studies have identified distinct surface markers of cancer stem cell (CSC) in HCC, and they also put forward the profound involvement of altered signaling pathways and epigenetic modifications in CSCs, in addition to the concomitant drug resistance and metastasis. Taken together, multiple key genetic and non-genetic factors, as well as liver CSCs, result in the development and progression of HCC.

Clinical and Prognostic Significance of ANCCA in Squamous Cell Lung Carcinoma Patients

BACKGROUND AND AIMS

AAA+ nuclear coregulatory cancer associated (ANCCA) has been demonstrated as playing important roles in diverse biological processes including tumorigenesis. However, the clinical and prognostic significance of ANCCA in squamous cell lung carcinoma (SCLC) patients is still unknown. The aim of this study is to identify the role of ANCCA in SCLC patients.

METHODS

ANCCA mRNA and protein expressions were detected in SCLC tissues and cell lines by real-time PCR and Western blot. We examined the ANCCA protein expression in 152 SCLC samples by immunohistochemistry and analyzed the association between the expression of ANCCA protein and clinicopathological characteristics of SCLC patients.

RESULTS

ANCCA mRNA and protein expression are increased in SCLC tissues and cell lines. Moreover, ANCCA protein overexpression was associated with differentiated degree, clinical stage, lymph node metastasis, and distant metastasis. In uni- and multivariate analyses, ANCCA protein overexpression was an independent poor prognostic factor for SCLC patients.

CONCLUSIONS

ANCCA act as a potential biomarker for therapeutic strategy and prognostic prediction for SCLC.

Transcriptome-Wide Analysis of Hepatitis B Virus-Mediated Changes to Normal Hepatocyte Gene Expression

Globally, a chronic hepatitis B virus (HBV) infection remains the leading cause of primary liver cancer. The mechanisms leading to the development of HBV-associated liver cancer remain incompletely understood. In part, this is because studies have been limited by the lack of effective model systems that are both readily available and mimic the cellular environment of a normal hepatocyte. Additionally, many studies have focused on single, specific factors or pathways that may be affected by HBV, without addressing cell physiology as a whole. Here, we apply RNA-seq technology to investigate transcriptome-wide, HBV-mediated changes in gene expression to identify single factors and pathways as well as networks of genes and pathways that are affected in the context of HBV replication. Importantly, these studies were conducted in an ex vivo model of cultured primary hepatocytes, allowing for the transcriptomic characterization of this model system and an investigation of early HBV-mediated effects in a biologically relevant context. We analyzed differential gene expression within the context of time-mediated gene-expression changes and show that in the context of HBV replication a number of genes and cellular pathways are altered, including those associated with metabolism, cell cycle regulation, and lipid biosynthesis. Multiple analysis pipelines, as well as qRT-PCR and an independent, replicate RNA-seq analysis, were used to identify and confirm differentially expressed genes. HBV-mediated alterations to the transcriptome that we identified likely represent early changes to hepatocytes following an HBV infection, suggesting potential targets for early therapeutic intervention. Overall, these studies have produced a valuable resource that can be used to expand our understanding of the complex network of host-virus interactions and the impact of HBV-mediated changes to normal hepatocyte physiology on viral replication.

Identification of important long non-coding RNAs and highly recurrent aberrant alternative splicing events in hepatocellular carcinoma through integrative analysis of multiple RNA-Seq datasets

Hepatocellular carcinoma (HCC) is an aggressive and deadly cancer. The molecular pathogenesis of the disease remains poorly understood. To better understand HCC biology and explore potential biomarkers and therapeutic targets, we investigated the whole transcriptome of HCC. Considering the genetic heterogeneity of HCC, four datasets from four studies consisting of 15 pairs of HCC and adjacent normal samples were analyzed. We observed that the number of lncRNAs expressed in each HCC sample was consistently greater than the adjacent normal sample. Moreover, 15 lncRNAs were identified expressed in five to seven HCC tissues but were not detected in any adjacent normal tissue. Differential expression analysis detected 35 up- and 80 down-regulated lncRNAs in HCC samples compared with adjacent normal samples. In addition, five differentially expressed lncRNAs were predicted to play a role in oxidation and reduction process. With regard to splicing alterations, we identified nine highly recurrent differential splicing events belonging to eight genes USO1, RPS24, CCDC50, THNSL2, NUMB, FN1 (two events), SLC39A14 and NR1I3. Of them, splicing alterations of SLC39A14 and NR1I3 were reported for the association with HCC for the first time. The splicing dysregulation in HCC may be influenced by three splicing factors ESRP2, CELF2 and SRSF5 which were significantly down-regulated in HCC samples. This study revealed uncharacterized aspects of HCC transcriptome and identified important lncRNAs and splicing isoforms with the potential to serve as biomarkers and therapeutic targets for the disease.

CAGE profiling of ncRNAs in hepatocellular carcinoma reveals widespread activation of retroviral LTR promoters in virus-induced tumors

An increasing number of noncoding RNAs (ncRNAs) have been implicated in various human diseases including cancer; however, the ncRNA transcriptome of hepatocellular carcinoma (HCC) is largely unexplored. We used CAGE to map transcription start sites across various types of human and mouse HCCs with emphasis on ncRNAs distant from protein-coding genes. Here, we report that retroviral LTR promoters, expressed in healthy tissues such as testis and placenta but not liver, are widely activated in liver tumors. Despite HCC heterogeneity, a subset of LTR-derived ncRNAs were more than 10-fold up-regulated in the vast majority of samples. HCCs with a high LTR activity mostly had a viral etiology, were less differentiated, and showed higher risk of recurrence. ChIP-seq data show that MYC and MAX are associated with ncRNA deregulation. Globally, CAGE enabled us to build a mammalian promoter map for HCC, which uncovers a new layer of complexity in HCC genomics.

Cellular transcriptomics: gelsolin negatively regulates the expression of apoptosis-associated genes and inhibits apoptosis in hepatocarcinoma cells

Gelsolin (GSN), which is a Ca(2+)-dependent actin filament severing and capping protein, plays a critical role in the cancer progress and has the potential for providing a novel thread for cancer therapy. In current study, we demonstrate the roles of GSN on anti-apoptosis of hepatocarcinoma cells by transcriptome RNA-seq method. Then flow cytometry (FCM), in-cell immunoblotting and transmission electron microscopy (TEM) were used to examine the GSN regulatory cell apoptosis. The results revealed GSN significantly suppresses apoptosis-associated functional categories through down-regulating apoptosis-associated genes in 5 apoptosis terms and 6 relevant KEGG pathways. FCM showed a significant lower apoptotic rate in GSN-SMMC7721 (P<0.05). In-cell immunoblotting detected discrepant expression of the apoptosis factors among GSN expressed/shRNA transfectants (P<0.05). TEM observed the discernible apoptosis morphology. Above results suggest a negative relationship between GSN expression and hepatocarcinoma cell apoptosis. GSN overexpression suppresses apoptosis while down-regulated GSN promotes apoptosis. The possible mechanism could be associated with the regulation of GSN on the apoptosis-associated pathways and the apoptosis factors caspase 3 and bcl-2.

Gene expression and splicing alterations analyzed by high throughput RNA sequencing of chronic lymphocytic leukemia specimens

Background

To determine differentially expressed and spliced RNA transcripts in chronic lymphocytic leukemia specimens a high throughput RNA-sequencing (HTS RNA-seq) analysis was performed.

Methods

Ten CLL specimens and five normal peripheral blood CD19+ B cells were analyzed by HTS RNA-seq. The library preparation was performed with Illumina TrueSeq RNA kit and analyzed by Illumina HiSeq 2000 sequencing system.

Results

An average of 48.5 million reads for B cells, and 50.6 million reads for CLL specimens were obtained with 10396 and 10448 assembled transcripts for normal B cells and primary CLL specimens respectively. With the Cuffdiff analysis, 2091 differentially expressed genes (DEG) between B cells and CLL specimens based on FPKM (fragments per kilobase of transcript per million reads and false discovery rate, FDR q < 0.05, fold change >2) were identified. Expression of selected DEGs (n = 32) with up regulated and down regulated expression in CLL from RNA-seq data were also analyzed by qRT-PCR in a test cohort of CLL specimens. Even though there was a variation in fold expression of DEG genes between RNA-seq and qRT-PCR; more than 90 % of analyzed genes were validated by qRT-PCR analysis. Analysis of RNA-seq data for splicing alterations in CLL and B cells was performed by Multivariate Analysis of Transcript Splicing (MATS analysis). Skipped exon was the most frequent splicing alteration in CLL specimens with 128 significant events (P-value <0.05, minimum inclusion level difference >0.1).

Conclusion

The RNA-seq analysis of CLL specimens identifies novel DEG and alternatively spliced genes that are potential prognostic markers and therapeutic targets. High level of validation by qRT-PCR for a number of DEG genes supports the accuracy of this analysis. Global comparison of transcriptomes of B cells, IGVH non-mutated CLL (U-CLL) and mutated CLL specimens (M-CLL) with multidimensional scaling analysis was able to segregate CLL and B cell transcriptomes but the M-CLL and U-CLL transcriptomes were indistinguishable. The analysis of HTS RNA-seq data to identify alternative splicing events and other genetic abnormalities specific to CLL is an added advantage of RNA-seq that is not feasible with other genome wide analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1708-9) contains supplementary material, which is available to authorized users.

SNP rs2057482 in HIF1A gene predicts clinical outcome of aggressive hepatocellular carcinoma patients after surgery

Hypoxia-inducible factor 1α (HIF-1α) plays an important role in tumor growth and metastasis. Genetic variations of HIF1A gene have been shown to influence the developing risk and prognosis in many types of human malignancies. However, their association with clinical outcomes of hepatocellular carcinoma (HCC) patients remains unclear. To investigate the predictive role of single nucleotide polymorphisms (SNPs) in HIF1A gene in HCC patients’ outcomes, we genotyped three functional SNPs (rs2057482, rs1957757 and rs2301113) in HIF1A gene and assessed their associations with clinicopathological parameters and prognosis of 492 surgical HCC patients. The patients with variant alleles (CT+TT) of SNP rs2057482 had a significantly lower recurrence risk when compared with patients with the CC genotype. In stratified analysis, the protective effect of rs2057482 CT+TT genotype was more evident in patients with adverse strata, compared with patients with favorable strata. Additionally, strong joint predictive effect between rs2057482 genotypes and AFP level, stage or differentiation were observed. Functional assay also indicated the significant effect of rs2057482 on gene expression. In conclusion, SNP rs2057482 in HIF1A gene is significantly associated with clinical outcomes of Chinese HCC patients after surgery, especially in those with aggressive status, which warrants further validation in other patient populations.

ATAD2 as a Poor Prognostic Marker for Hepatocellular Carcinoma after Curative Resection

Purpose

Cancer cells frequently express genes that are specifically or preferentially expressed in male germ cells under normal conditions. The ATPase family AAA domain-containing 2 (ATAD2) is one such and works as an important cofactor for MYC-dependent transcription. In hepatocellular carcinoma (HCC), ATAD2 has been identified as a candidate driver gene located within the amplified 8q24 locus. However, the prognostic significance of ATAD2 protein expression in HCC remains uncertain.

Materials and Methods

We investigated ATAD2 protein expression by immunohistochemistry in tumor tissue from 182 HCC patients who underwent curative resection. Associations of ATAD2 expression with clinicopathologic variables or prognosis of HCC patients were analyzed.

Results

ATAD2 expression was observed in 119 (65.4%) of the 182 HCCs and tended to be independent predictor of early recurrence (p=0.059). ATAD2 expression showed an unfavorable influence on recurrence-free survival (RFS) (p < 0.001). Subgroup analysis among patients with tumor size ≤ 5.0 cm (n=109), patients at Barcelona Clinic Liver Cancer stage 0 or A (n=92), and patients with α-fetoprotein ≤ 20 ng/mL (n=61), the ATAD2-positive groups unfavorably influenced RFS (p=0.008, p=0.009, and p=0.013, respectively). In addition, ATAD2 expression was an independent predictor of shorter RFS (p=0.002). ATAD2 expression showed an unfavorable influence on disease-specific survival (p=0.001), but was not an independent predictor of shorter disease-specific survival (p=0.109).

Conclusion

ATAD2 protein expression may be a potential predictor of RFS in HCC patients after curative resection and ATAD2 may have prognostic value in patients with early stage HCC or normal serum α-fetoprotein level.

Employment of digital gene expression profiling to identify potential pathogenic and therapeutic targets of fulminant hepatic failure

Background

The dysregulated cytokine metabolism and activity are crucial to the development of fulminant hepatic failure (FHF), and many different cytokines have been identified. However, the precise gene expression profile and their interactions association with FHF are yet to be further elucidated.

Methods

In this study, we detected the digital gene expression profile (DGEP) by high-throughput sequencing in normal and FHF mouse liver, and the candidate genes and potential targets for FHF therapy were verified. And the FHF mouse model was induced by D-Galactosamine (GalN)/lipopolysaccharide (LPS).

Results

Totally 12727 genes were detected, and 3551 differentially expressed genes (DEGs) were obtained from RNA-seq data in FHF mouse liver. In FHF mouse liver, many of those DEGs were identified as differentially expressed in metabolic process, biosynthetic process, response to stimulus and response to stress, etc. Similarly, pathway enrichment analysis in FHF mouse liver showed that many significantly DEGs were also enriched in metabolic pathways, apoptosis, chemokine signaling pathways, etc. Considering the important role of nuclear factor-kappa B (NF-κB) in metabolic regulation and delicate balance between cell survival and death, several DEGs involved in NF-κB pathway were selected for experimental validation. As compared to normal control, NF-κBp65 and its inhibitory protein IκBα were both significantly increased, and NF-κB targeted genes including tumor necrosis factor α(TNFα), inducible nitric oxide synthase (iNOS), interleukin-1β, chemokines CCL3 and CCL4 were also increased in hepatic tissues of FHF. In addition, after NF-κB was successfully pre-blocked, there were significant alteration of hepatic pathological damage and mortality of FHF mouse model.

Conclusions

This study provides the globe gene expression profile of FHF mouse liver, and demonstrates the possibility of NF-κB gene as a potential therapeutic target for FHF.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-015-0380-9) contains supplementary material, which is available to authorized users.

Liver cancer oncogenomics: opportunities and dilemmas for clinical applications

Primary liver cancers are among the most rapidly evolving malignant tumors worldwide. An underlying chronic inflammatory liver disease, which precedes liver cancer development for several decades and frequently creates a pro-oncogenic microenvironment, impairs progress in therapeutic approaches. Molecular heterogeneity of liver cancer is potentiated by a crosstalk between epithelial tumor and stromal cells that complicate translational efforts to unravel molecular mechanisms of hepatocarcinogenesis with a drugable intend. Next-generation sequencing has greatly advanced our understanding of cancer development. With regards to liver cancer, the unprecedented coverage of next-generation sequencing has created a detailed map of genetic alterations and identified key somatic changes such as CTNNB1 and TP53 as well as several previously unrecognized recurrent disease-causing alterations that could contribute to new therapeutic approaches. Importantly, these investigations indicate that a classical oncogene addiction cannot be assumed for primary liver cancer. Therefore, hepatocarcinogenesis can be considered a paradigm suitable for individualized medicine.

Examining the effects of nuclear GPCRs on gene expression using isolated nuclei

The measurement of changes in the transcriptome is a common end point for various pathologic and pharmacologic studies. In recent years, with the discovery of a host of potential pharmacologic targets located directly on the nuclear membrane, the need to assess their potential control over the transcriptome has arisen. Here we present techniques for assessing changes in gene expression in isolated nuclei in response to stimulation by endogenous GPCRs on the nuclear membrane.