Clinicopathological-Associated Regulatory Network of Deregulated circRNAs in Hepatocellular Carcinoma

Simple Summary

Here, we present a novel strategy to identify key signatures of clinically-relevant co-expressed circRNA-mRNA networks in pertinent cancer-pathways that modulate the prognosis of HCC patients, by integrating clinicopathological features, circRNA and mRNA expression profiles. Five master circRNAs were identified and experimentally demonstrated to upregulate proliferate and promote transformation. Through further integration with miRNA-expression profiles, clinically-relevant competing-endogenous-RNA (ceRNA) networks of circRNA-miRNA-mRNAs were constructed. The most up-regulated nodal-circRNA, circGPC3 was experimentally demonstrated to up-regulate cell-cycle, migration and invasion. circGPC3 was found to act as a sponge of miR-378a-3p to regulate ASPM expression and modulate cell transformation. These 5 nodal circRNAs has potential to be good prognostic biomarkers with good prognostic performance. circGPC3 has great potential to be a promising non-invasive prognostic biomarker for early HCC. We have thus demonstrated the robustness of bioinformatically-predicted master circRNAs in clinically-relevant, circRNA-mRNA networks, underscoring the important roles that these identified deregulated key/master circRNAs play in HCC.

Abstract

Hepatocellular carcinoma (HCC) is one of the most common and lethal cancers worldwide. Here, we present a novel strategy to identify key circRNA signatures of clinically relevant co-expressed circRNA-mRNA networks in pertinent cancer-pathways that modulate prognosis of HCC patients, by integrating clinic-pathological features, circRNA and mRNA expression profiles. Through further integration with miRNA expression profiles, clinically relevant competing-endogenous-RNA (ceRNA) networks of circRNA-miRNA-mRNAs were constructed. At least five clinically relevant nodal-circRNAs, co-expressed with numerous genes, were identified from the circRNA-mRNA networks. These nodal circRNAs upregulated proliferation (except circRaly) and transformation in cells. The most upregulated nodal-circRNA, circGPC3, associated with higher-grade tumors and co-expressed with 33 genes, competes with 11 mRNAs for two shared miRNAs. circGPC3 was experimentally demonstrated to upregulate cell-cycle and migration/invasion in both transformed and non-transformed liver cell-lines. circGPC3 was further shown to act as a sponge of miR-378a-3p to regulate APSM (Abnormal spindle-like microcephaly associated) expression and modulate cell transformation. This study identifies 5 key nodal master circRNAs in a clinically relevant circRNA-centric network that are significantly associated with poorer prognosis of HCC patients and promotes tumorigenesis in cell-lines. The identification and characterization of these key circRNAs in clinically relevant circRNA-mRNA and ceRNA networks may facilitate the design of novel strategies targeting these important regulators for better HCC prognosis.

SwissMTB: establishing comprehensive molecular cancer diagnostics in Swiss clinics

Background

Molecular precision oncology is an emerging practice to improve cancer therapy by decreasing the risk of choosing treatments that lack efficacy or cause adverse events. However, the challenges of integrating molecular profiling into routine clinical care are manifold. From a computational perspective these include the importance of a short analysis turnaround time, the interpretation of complex drug-gene and gene-gene interactions, and the necessity of standardized high-quality workflows. In addition, difficulties faced when integrating molecular diagnostics into clinical practice are ethical concerns, legal requirements, and limited availability of treatment options beyond standard of care as well as the overall lack of awareness of their existence.

Methods

To the best of our knowledge, we are the first group in Switzerland that established a workflow for personalized diagnostics based on comprehensive high-throughput sequencing of tumors at the clinic. Our workflow, named SwissMTB (Swiss Molecular Tumor Board), links genetic tumor alterations and gene expression to therapeutic options and clinical trial opportunities. The resulting treatment recommendations are summarized in a clinical report and discussed in a molecular tumor board at the clinic to support therapy decisions.

Results

Here we present results from an observational pilot study including 22 late-stage cancer patients. In this study we were able to identify actionable variants and corresponding therapies for 19 patients. Half of the patients were analyzed retrospectively. In two patients we identified resistance-associated variants explaining lack of therapy response. For five out of eleven patients analyzed before treatment the SwissMTB diagnostic influenced treatment decision.

Conclusions

SwissMTB enables the analysis and clinical interpretation of large numbers of potentially actionable molecular targets. Thus, our workflow paves the way towards a more frequent use of comprehensive molecular diagnostics in Swiss hospitals.

Electronic supplementary material

The online version of this article (10.1186/s12911-018-0680-0) contains supplementary material, which is available to authorized users.

NGS-pipe: a flexible, easily extendable and highly configurable framework for NGS analysis

Motivation

Next-generation sequencing is now an established method in genomics, and massive amounts of sequencing data are being generated on a regular basis. Analysis of the sequencing data is typically performed by lab-specific in-house solutions, but the agreement of results from different facilities is often small. General standards for quality control, reproducibility and documentation are missing.

Results

We developed NGS-pipe, a flexible, transparent and easy-to-use framework for the design of pipelines to analyze whole-exome, whole-genome and transcriptome sequencing data. NGS-pipe facilitates the harmonization of genomic data analysis by supporting quality control, documentation, reproducibility, parallelization and easy adaptation to other NGS experiments.

Availability and implementation

https://github.com/cbg-ethz/NGS-pipe

Genomic variant annotation workflow for clinical applications

Annotation and interpretation of DNA aberrations identified through next-generation sequencing is becoming an increasingly important task. Even more so in the context of data analysis pipelines for medical applications, where genomic aberrations are associated with phenotypic and clinical features. Here we describe a workflow to identify potential gene targets in aberrated genes or pathways and their corresponding drugs. To this end, we provide the R/Bioconductor package rDGIdb, an R wrapper to query the drug-gene interaction database (DGIdb). DGIdb accumulates drug-gene interaction data from 15 different resources and allows filtering on different levels. The rDGIdb package makes these resources and tools available to R users. Moreover, rDGIdb queries can be automated through incorporation of the rDGIdb package into NGS sequencing pipelines.

Genomic variant annotation workflow for clinical applications

Annotation and interpretation of DNA aberrations identified through next-generation sequencing is becoming an increasingly important task. Even more so in the context of data analysis pipelines for medical applications, where genomic aberrations are associated with phenotypic and clinical features. Here we describe a workflow to identify potential gene targets in aberrated genes or pathways and their corresponding drugs. To this end, we provide the R/Bioconductor package rDGIdb, an R wrapper to query the drug-gene interaction database (DGIdb). DGIdb accumulates drug-gene interaction data from 15 different source databases and allows filtering on different levels. The rDGIdb package makes these resources and tools available to R users. Moreover, DGIdb queries can be automated through incorporation of the rDGIdb package into NGS sequencing pipelines.

Global re-wiring of p53 transcription regulation by the hepatitis B virus X protein

BACKGROUND

The tumour suppressor p53 is a central player in transcription regulation and cell fate determination. By interacting with p53 and altering its sequence-specific binding to the response elements, the hepatitis B virus X protein (HBx) was reported to re-direct p53 regulation of some genes.

RESULTS

Coupling massively parallel deep sequencing with p53 chromatin immunoprecipitation, we demonstrate that HBx modulates global p53 site selection and that this was strongly influenced by altered interaction with transcription co-factors/co-regulators as well as post-translational modifications. Specifically, HBx attenuated p53-TBP-RB1 transcription complex recruitment and interaction and this was associated with hyper-phosphorylation of p53 at serine 315 by HBx. Concurrently, HBx enhanced p53 DNA occupancy to other response elements either alone by displacing specific transcription factors such as CEBPB and NFkB1, or in complex with distinct interacting co-factors Sp1, JUN and E2F1. Importantly, re-wiring of p53 transcription regulation by HBx was linked to the deregulation of genes involved in cell proliferation and death, suggesting a role of HBx in errant cell fate determination mediated by altered p53 site selection of target genes.

CONCLUSIONS

Our study thus presents first evidence of global modes of p53 transcription alteration by HBx and provides new insights to understand and potentially curtail the viral oncoprotein.

Differentially Expressed miRNAs in Hepatocellular Carcinoma Target Genes in the Genetic Information Processing and Metabolism Pathways

To date, studies of the roles of microRNAs (miRNAs) in hepatocellular carcinoma (HCC) have either focused on specific individual miRNAs and a small number of suspected targets or simply reported a list of differentially expressed miRNAs based on expression profiling. Here, we seek a more in-depth understanding of the roles of miRNAs and their targets in HCC by integrating the miRNA and messenger RNA (mRNA) expression profiles of tumorous and adjacent non-tumorous liver tissues of 100 HCC patients. We assessed the levels of 829 mature miRNAs, of which 32 were significantly differentially expressed. Statistical analysis indicates that six of these miRNAs regulate a significant proportion of their in silico predicted target mRNAs. Three of these miRNAs (miR-26a, miR-122, and miR-130a) were down-regulated in HCC, and their up-regulated gene targets are primarily associated with aberrant cell proliferation that involves DNA replication, transcription and nucleotide metabolism. The other three miRNAs (miR-21, miR-93, and miR-221) were up-regulated in HCC, and their down-regulated gene targets are primarily involved in metabolism and immune system processes. We further found evidence for a coordinated miRNA-induced regulation of important cellular processes, a finding to be considered when designing therapeutic applications based on miRNAs.

Methylation profiles reveal distinct subgroup of hepatocellular carcinoma patients with poor prognosis

Hepatocellular Carcinoma (HCC) is one of the leading causes of cancer-associated mortality worldwide. However, the role of epigenetic changes such as aberrant DNA methylation in hepatocarcinogenesis remains largely unclear. In this study, we examined the methylation profiles of 59 HCC patients. Using consensus hierarchical clustering with feature selection, we identified three tumor subgroups based on their methylation profiles and correlated these subgroups with clinicopathological parameters. Interestingly, one tumor subgroup is different from the other 2 subgroups and the methylation profile of this subgroup is the most distinctly different from the non-tumorous liver tissues. Significantly, this subgroup of patients was found to be associated with poor overall as well as disease-free survival. To further understand the pathways modulated by the deregulation of methylation in HCC patients, we integrated data from both the methylation as well as the gene expression profiles of these 59 HCC patients. In these patients, while 4416 CpG sites were differentially methylated between the tumors compared to the adjacent non-tumorous tissues, only 536 of these CpG sites were associated with differences in the expression of their associated genes. Pathway analysis revealed that forty-four percent of the most significant upstream regulators of these 536 genes were involved in inflammation-related NFκB pathway. These data suggest that inflammation via the NFκB pathway play an important role in modulating gene expression of HCC patients through methylation. Overall, our analysis provides an understanding on aberrant methylation profile in HCC patients.

Population estimation of human embryonic stem cell cultures

Traditionally, the population of human embryonic stem cell (hESC) culture is estimated through haemacytometer counts, which include harvesting the cells and manually analyzing a fraction of an entire population. Obviously, through this highly invasive method, it is not possible to preserve any spatial information on the cell population. The goal of this study is to identify a fast and consistent method for in situ automated hESC population estimation to quantitatively estimate the cell growth. Therefore, cell cultures were fixed, stained, and their nuclei imaged through high-resolution microscopy, and the images were processed with different image analysis techniques. The proposed method first identifies signal and background by computing an image specific threshold for image segmentation. By applying a morphological operator (watershed), we split most physically overlapping nuclei, leading to a pixel area distribution of isolated signal areas on the image. On the basis of this distribution, we derive a nucleus area model, describing the distribution of the area of cell debris, single nuclei, and small groups of connected nuclei. Through the model, we can give a quantitative estimation of the population. The focus of this study is on low-density human embryonic stem cell populations; hence cultures were measured at days 2-3 after seeding. Compared with manual cell counts, the automatic method achieved higher accuracy with <6% error.

The types and frequencies of improved nonmusculoskeletal symptoms reported after chiropractic spinal manipulative therapy

OBJECTIVE

To investigate the frequency and types of improved nonmusculoskeletal symptoms reported after chiropractic spinal manipulative therapy.

DESIGN

Retrospective information obtained by chiropractors through standardized interview of patients on return visit within 2 weeks of previous treatment.

SETTING

The private practice of 87 Swedish chiropractors (response rate 81%).

SUBJECTS

Twenty consecutive (presumably naive) patients per chiropractor (1504 valid questionnaires returned, 86% of optimal number of replies).

INTERVENTION

Spinal manipulation with or without additional therapy provided by chiropractors.

MAIN OUTCOME MEASURES

Self-reported improved nonmusculoskeletal symptoms (reactions).

RESULTS

At least I reaction was reported after the previous treatment in 21% to 25% of cases. Of these responses, 26% were related to the airway passages (usually reported as "easier to breathe"), 25% were related to the digestive system (mostly reported as "improved function"), 14% were classified under eyes/vision (usually reported as "improved vision"), and 14% under heart/ circulation (about half of these reported as "improved circulation"). The number of spinal areas treated was positively associated with the number of reactions.

CONCLUSION

A minority of chiropractic patients report having positive nonmusculoskeletal reactions after spinal manipulative therapy but such reports cluster predominantly around specific symptoms. It would be interesting to find out if these can be verified objectively and, if so, to investigate if they are caused by the treatment or if they are signs of natural variations in human physiology.