Fluorescence-based alternative splicing reporters for the study of epithelial plasticity in vivo

Gene editing of hematopoietic stem cells restores T-cell response in familial hemophagocytic lymphohistiocytosis

BACKGROUND

Hemophagocytic lymphohistiocytosis (HLH) is a hyperinflammatory disorder characterized by a life-threatening cytokine storm and immunopathology. Familial HLH type 3 (FHL3) accounts for approximately 30% of all inborn HLH cases worldwide. It is caused by mutations in the UNC13D gene that result in impaired degranulation of cytotoxic vesicles and hence compromised T-cell- and natural killer-cell-mediated killing. Current treatment protocols, including allogeneic hematopoietic stem cell (HSC) transplantation, still show high mortality.

OBJECTIVE

We sought to develop and evaluate a curative genome editing strategy in the preclinical FHL3 Jinx mouse model. Jinx mice harbor a cryptic splice donor site in Unc13d intron 26 and develop clinical symptoms of human FHL3 upon infection with lymphocytic choriomeningitis virus (LCMV).

METHODS

We employed clustered regularly interspaced short palindromic repeats (CRISPR)-Cas technology to delete the disease-causing mutation in HSCs and transplanted Unc13d-edited stem cells into busulfan-conditioned Jinx recipient mice. Safety studies included extensive genotyping and chromosomal aberrations analysis by single targeted linker-mediated PCR sequencing (CAST-Seq)-based off-target analyses. Cure from HLH predisposition was assessed by LCMV infection.

RESULTS

Hematopoietic cells isolated from transplanted mice revealed efficient gene editing (>95%), polyclonality of the T-cell receptor repertoire, and neither signs of off-target effects nor leukemogenesis. Unc13d transcription levels of edited and wild-type cells were comparable. While LCMV challenge resulted in acute HLH in Jinx mice transplanted with mock-edited HSCs, Jinx mice grafted with Unc13d-edited cells showed rapid virus clearance and protection from HLH.

CONCLUSIONS

Our study demonstrates that transplantation of CRISPR-Cas edited HSCs supports the development of a functional polyclonal T-cell response in the absence of genotoxicity-associated clonal outgrowth.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

A synthetic lethal screen for Snail-induced enzalutamide resistance identifies JAK/STAT signaling as a therapeutic vulnerability in prostate cancer

Despite substantial improvements in the treatment landscape of prostate cancer, the evolution of hormone therapy-resistant and metastatic prostate cancer remains a major cause of cancer-related death globally. The mainstay of treatment for advanced prostate cancer is targeting of androgen receptor signaling, including androgen deprivation therapy plus second-generation androgen receptor blockade (e.g., enzalutamide, apalutamide, darolutamide), and/or androgen synthesis inhibition (abiraterone). While these agents have significantly prolonged the lives of patients with advanced prostate cancer, is nearly universal. This therapy resistance is mediated by diverse mechanisms, including both androgen receptor-dependent mechanisms, such as androgen receptor mutations, amplifications, alternative splicing, and amplification, as well as non-androgen receptor-mediated mechanisms, such as lineage plasticity toward neuroendocrine-like or epithelial-mesenchymal transition (EMT)-like lineages. Our prior work identified the EMT transcriptional regulator Snail as critical to hormonal therapy resistance and is commonly detected in human metastatic prostate cancer. In the current study, we sought to interrogate the actionable landscape of EMT-mediated hormone therapy resistant prostate cancer to identify synthetic lethality and collateral sensitivity approaches to treating this aggressive, therapy-resistant disease state. Using a combination of high-throughput drug screens and multi-parameter phenotyping by confluence imaging, ATP production, and phenotypic plasticity reporters of EMT, we identified candidate synthetic lethalities to Snail-mediated EMT in prostate cancer. These analyses identified multiple actionable targets, such as XPO1, PI3K/mTOR, aurora kinases, c-MET, polo-like kinases, and JAK/STAT as synthetic lethalities in Snail+ prostate cancer. We validated these targets in a subsequent validation screen in an LNCaP-derived model of resistance to sequential androgen deprivation and enzalutamide. This follow-up screen provided validation of inhibitors of JAK/STAT and PI3K/mTOR as therapeutic vulnerabilities for both Snail+ and enzalutamide-resistant prostate cancer.

A ratiometric dual luciferase reporter for quantitative monitoring of pre-mRNA splicing efficiency in vivo

Precursor messenger RNA (pre-mRNA) splicing is critical for cell growth and development, and errors in RNA splicing frequently cause cellular dysfunction, abnormal gene expression, and a variety of human diseases. However, there is currently a lack of reliable systems to noninvasively monitor the mRNA splicing efficiency in cells and animals. Here, we described the design of a genetically engineered ratiometric dual luciferase reporter to continuously quantify the changes in mRNA splice variants in vivo. This reporter system is encoded within a single polypeptide but on separate exons, thus generating two distinct luciferase signals derived from spliced and unspliced mRNAs. With this reporter, the two kinds of luciferase in the same individual can minimize the influence of indirect factors on splicing, and the ratio of these two luciferase intensities represents the dynamic splicing efficiency of pre-mRNA. Our study offers a convenient and robust tool for the screening and identification of small molecules or trans-acting factors that affect the efficiency of specific splicing reactions.

U2AF2 binds IL7R exon 6 ectopically and represses its inclusion

Interleukin 7 receptor α-chain is crucial for the development and maintenance of T cells and is genetically associated with autoimmune disorders including multiple sclerosis (MS), a demyelinating disease of the CNS. Exon 6 of IL7R encodes for the transmembrane domain of the receptor and is regulated by alternative splicing: inclusion or skipping of IL7R exon 6 results in membrane-bound or soluble IL7R isoforms, respectively. We previously identified a SNP (rs6897932) in IL7R exon 6, strongly associated with MS risk and showed that the risk allele (C) increases skipping of the exon, resulting in elevated levels of sIL7R. This has important pathological consequences as elevated levels of sIL7R has been shown to exacerbate the disease in the experimental autoimmune encephalomyelitis mouse model of MS. Understanding the regulation of exon 6 splicing provides important mechanistic insights into the pathogenesis of MS. Here we report two mechanisms by which IL7R exon 6 is controlled. First, a competition between PTBP1 and U2AF2 at the polypyrimidine tract (PPT) of intron 5, and second, an unexpected U2AF2-mediated assembly of spicing factors in the exon. We noted the presence of a branchpoint sequence (BPS) (TACTAAT or TACTAAC) within exon 6, which is stronger with the C allele. We also noted that the BPS is followed by a PPT and conjectured that silencing could be mediated by the binding of U2AF2 to that tract. In support of this model, we show that evolutionary conservation of the exonic PPT correlates well with the degree of alternative splicing of exon 6 in two non-human primate species and that U2AF2 binding to this PPT recruits U2 snRNP components to the exon. These observations provide the first explanation for the stronger silencing of IL7R exon 6 with the disease associated C allele at rs6897932.

Towards control of cellular decision-making networks in the epithelial-to-mesenchymal transition

We present the epithelial-to-mesenchymal transition (EMT) from two perspectives: experimental/technological and theoretical. We review the state of the current understanding of the regulatory networks that underlie EMT in three physiological contexts: embryonic development, wound healing, and metastasis. We describe the existing experimental systems and manipulations used to better understand the molecular participants and factors that influence EMT and metastasis. We review the mathematical models of the regulatory networks involved in EMT, with a particular emphasis on the network motifs (such as coupled feedback loops) that can generate intermediate hybrid states between the epithelial and mesenchymal states. Ultimately, the understanding gained about these networks should be translated into methods to control phenotypic outcomes, especially in the context of cancer therapeutic strategies. We present emerging theories of how to drive the dynamics of a network toward a desired dynamical attractor (e.g. an epithelial cell state) and emerging synthetic biology technologies to monitor and control the state of cells.

Revealing the alternative promoter usage of SAF/MAZ gene by bichromatic fluorescent reporter construct

The large-scale identification of putative alternative promoters study shows more than 52% of human genes are regulated by alternative promoters. The human myc-associated zinc finger protein (SAF/MAZ) gene have SAF-1 and SAF-3 variants transcripted from two transcription start sites (TSSs). By using SAF/MAZ promoter as a model, we set up an approach to probe how the alternative promoters are regulated in real time. We have constructed the bichromatic fluorescent reporter driven by SAF/MAZ 5'-proximal promoter plasmids from which transactivation status of SAF-1 and SAF-3 alternative promoter could be monitored by EGFP and DsRed expression respectively. The results showed that the SAF-3 expression is regulated by alternative promoters. When the bichromatic fluorescent reporter was driven by -1692/+277 or -1401/+277 SAF/MAZ promoter the dominant expression of SAF-3 would be observed in comparison with SAF-1 expression. We also identified that Elk-1 is an inhibitory transcription factor for SAF-3 expression. The temporal diversity of SAF-1 and SAF-3 expressions can be observed via bichromatic fluorescent reporters. These imply that the bichromatic fluorescent reporter driven by alternative promoter construct might be a useful tool for decoding the temporal regulatory repertoire of alternative promoter in human genes.

Environmental influences on RNA processing: Biochemical, molecular and genetic regulators of cellular response

RNA processing has emerged as a key mechanistic step in the regulation of the cellular response to environmental perturbation. Recent work has uncovered extensive remodeling of transcriptome composition upon environmental perturbation and linked the impacts of this molecular plasticity to health and disease outcomes. These isoform changes and their underlying mechanisms are varied-involving alternative sites of transcription initiation, alternative splicing, and alternative cleavage at the 3' end of the mRNA. The mechanisms and consequences of differential RNA processing have been characterized across a range of common environmental insults, including chemical stimuli, immune stimuli, heat stress, and cancer pathogenesis. In each case, there are perturbation-specific contributions of local (cis) regulatory elements or global (trans) factors and downstream consequences. Overall, it is clear that choices in isoform usage involve a balance between the usage of specific genetic elements (i.e., splice sites, polyadenylation sites) and the timing at which certain decisions are made (i.e., transcription elongation rate). Fine-tuned cellular responses to environmental perturbation are often dependent on the genetic makeup of the cell. Genetic analyses of interindividual variation in splicing have identified genetic effects on splicing that contribute to variation in complex traits. Finally, the increase in the number of tissue types and environmental conditions analyzed for RNA processing is paralleled by the need to develop appropriate analytical tools. The combination of large datasets, novel methods and conditions explored promises to provide a much greater understanding of the role of RNA processing response in human phenotypic variation. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Interconnected feedback loops among ESRP1, HAS2, and CD44 regulate epithelial-mesenchymal plasticity in cancer

Aberrant activation of epithelial-mesenchymal transition (EMT) in carcinoma cells contributes to increased migration and invasion, metastasis, drug resistance, and tumor-initiating capacity. EMT is not always a binary process; rather, cells may exhibit a hybrid epithelial/mesenchymal (E/M) phenotype. ZEB1-a key transcription factor driving EMT-can both induce and maintain a mesenchymal phenotype. Recent studies have identified two novel autocrine feedback loops utilizing epithelial splicing regulatory protein 1 (ESRP1), hyaluronic acid synthase 2 (HAS2), and CD44 which maintain high levels of ZEB1. However, how the crosstalk between these feedback loops alters the dynamics of epithelial-hybrid-mesenchymal transition remains elusive. Here, using an integrated theoretical-experimental framework, we identify that these feedback loops can enable cells to stably maintain a hybrid E/M phenotype. Moreover, computational analysis identifies the regulation of ESRP1 as a crucial node, a prediction that is validated by experiments showing that knockdown of ESRP1 in stable hybrid E/M H1975 cells drives EMT. Finally, in multiple breast cancer datasets, high levels of ESRP1, ESRP1/HAS2, and ESRP1/ZEB1 correlate with poor prognosis, supporting the relevance of ZEB1/ESRP1 and ZEB1/HAS2 axes in tumor progression. Together, our results unravel how these interconnected feedback loops act in concert to regulate ZEB1 levels and to drive the dynamics of epithelial-hybrid-mesenchymal transition.

Construction of a tri-chromatic reporter cell line for the rapid and simple screening of splice-switching oligonucleotides targeting DMD exon 51 using high content screening

Splice-switching oligonucleotides (SSOs) that can modulate RNA splicing are used for the treatment of many genetic disorders. To enhance the efficacy of modulating splicing, it is important to optimize SSOs with regard to target sites, GC content, melting temperature (Tm value), chemistries, and lengths. Thus, in vitro assay systems that allow for the rapid and simple screening of SSOs are essential for optimizing SSO design. In this study, we established a novel tri-chromatic reporter cell line for SSO screening. This reporter cell line is designed to express three different fluorescent proteins (blue, green, and red) and was employed for high content screening (HCS, also known as high content analysis; HCA) for the evaluation of SSO-induced exon skipping by analyzing the expression levels of fluorescent proteins. The blue fluorescent protein is stably expressed throughout the cell and is useful for data normalization using cell numbers. Furthermore, both the green and red fluorescent proteins were used for monitoring the splicing patterns of target genes. Indeed, we demonstrated that this novel reporter cell line involving HCS leads to a more rapid and simple approach for the evaluation of exon skipping than widely used methods, such as RT-PCR, western blotting, and quantitative RT-PCR. Additionally, a brief screening of Locked nucleic acids (LNA)-based SSOs targeting exon 51 in DMD was performed using the reporter cell line. The LNA-based SSO cocktail shows high exon 51 skipping in a dose-dependent manner. Furthermore, the LNA-based SSO cocktails display high exon 51 skipping activities on endogenous DMD mRNA in human rhabdomyosarcoma cells.

Epithelial/mesenchymal plasticity: how have quantitative mathematical models helped improve our understanding?

Phenotypic plasticity, the ability of cells to reversibly alter their phenotypes in response to signals, presents a significant clinical challenge to treating solid tumors. Tumor cells utilize phenotypic plasticity to evade therapies, metastasize, and colonize distant organs. As a result, phenotypic plasticity can accelerate tumor progression. A well‐studied example of phenotypic plasticity is the bidirectional conversions among epithelial, mesenchymal, and hybrid epithelial/mesenchymal (E/M) phenotype(s). These conversions can alter a repertoire of cellular traits associated with multiple hallmarks of cancer, such as metabolism, immune evasion, invasion, and metastasis. To tackle the complexity and heterogeneity of these transitions, mathematical models have been developed that seek to capture the experimentally verified molecular mechanisms and act as ‘hypothesis‐generating machines’. Here, we discuss how these quantitative mathematical models have helped us explain existing experimental data, guided further experiments, and provided an improved conceptual framework for understanding how multiple intracellular and extracellular signals can drive E/M plasticity at both the single‐cell and population levels. We also discuss the implications of this plasticity in driving multiple aggressive facets of tumor progression.

EMT and MET: necessary or permissive for metastasis?

Epithelial‐to‐mesenchymal transition (EMT) and its reverse mesenchymal‐to‐epithelial transition (MET) have been suggested to play crucial roles in metastatic dissemination of carcinomas. These phenotypic transitions between states are not binary. Instead, carcinoma cells often exhibit a spectrum of epithelial/mesenchymal phenotype(s). While epithelial/mesenchymal plasticity has been observed preclinically and clinically, whether any of these phenotypic transitions are indispensable for metastatic outgrowth remains an unanswered question. Here, we focus on epithelial/mesenchymal plasticity in metastatic dissemination and propose alternative mechanisms for successful dissemination and metastases beyond the traditional EMT/MET view. We highlight multiple hypotheses that can help reconcile conflicting observations, and outline the next set of key questions that can offer valuable insights into mechanisms of metastasis in multiple tumor models.

Analysis of Alternative Pre-RNA Splicing in the Mouse Retina Using a Fluorescent Reporter

In vivo alternative splicing is controlled in a tissue and cell type specific manner. Often individual cellular components of complex tissues will express different splicing programs. Thus, when studying splicing in multicellular organisms it is critical to determine the exon inclusion levels in individual cells positioned in the context of their native tissue or organ. Here we describe how a fluorescent splicing reporter in combination with in vivo electroporation can be used to visualize alternative splicing in individual cells within mature tissues. In a test case we show how the splicing of a photoreceptor specific exon can be visualized within the mouse retina. The retina was chosen as an example of a complex tissue that is fragile and whose cells cannot be studied in culture. With minor modifications to the injection and electroporation procedure, the protocol we outline can be applied to other tissues and organs.

Fluorescent Protein-Based Quantification of Alternative Splicing of a Target Cassette Exon in Mammalian Cells

Alternative splicing is an important mechanism of regulation of gene expression and expansion of proteome complexity. Recently we developed a new fluorescence reporter for quantitative analysis of alternative splicing of a target cassette exon in live cells (Gurskaya et al., 2012). It consists of a specially designed minigene encoding red and green fluorescent proteins (Katushka and TagGFP2) and a fragment of the target gene between them. Skipping or inclusion of the alternative exon induces a frameshift; ie, alternative exon length must not be a multiple of 3. Finally, red and green fluorescence intensities of cells expressing this reporter are used to estimate the percentage of alternative (exon-skipped) and normal (exon-retained) transcripts. Here, we provide a detailed description of design and application of the fluorescence reporter of a target alternative exon splicing in mammalian cell lines.

In Vivo Visualization and Characterization of Epithelial-Mesenchymal Transition in Breast Tumors

The activation of the epithelial-to-mesenchymal transition (EMT) program is a critical step in cancer progression and metastasis, but visualization of this process at the single-cell level, especially in vivo, remains challenging. We established an in vivo approach to track the fate of tumor cells based on a novel EMT-driven fluorescent color switching breast cancer mouse model and intravital two-photon laser scanning microscopy. Specifically, the MMTV-PyMT, Rosa26-RFP-GFP, and Fsp1-Cre triple transgenic mouse model was used to monitor the conversion of RFP-positive epithelial cells to GFP-positive mesenchymal cells in mammary tumors under the control of the Fsp1 (ATL1) promoter, a gate-keeper of EMT initiation. RFP-positive cells were isolated from the tumors, sorted, and transplanted into mammary fat pads of SCID mice to monitor EMT during breast tumor formation. We found that the conversion from RFP- to GFP-positive and spindle-shaped cells was a gradual process, and that GFP-positive cells preferentially localized close to blood vessels, independent of tumor size. Furthermore, cells undergoing EMT expressed high levels of the HGF receptor, c-Met, and treatment of RFP-positive cells with the c-Met inhibitor, cabozantinib, suppressed the RFP-to-GFP conversion in vitro Moreover, administration of cabozantinib to mice with palpable RFP-positive tumors resulted in a silent EMT phenotype whereby GFP-positive cells exhibited reduced motility, leading to suppressed tumor growth. In conclusion, our imaging technique provides a novel opportunity for visualizing tumor EMT at the single-cell level and may help to reveal the intricacies underlying tumor dynamics and treatment responses. Cancer Res; 76(8); 2094-104. ©2016 AACR.

Distinct routes to metastasis: plasticity-dependent and plasticity-independent pathways

The cascade that culminates in macrometastases is thought to be mediated by phenotypic plasticity, including epithelial–mesenchymal and mesenchymal–epithelial transitions (EMT and MET). Although there is substantial support for the role of EMT in driving cancer cell invasion and dissemination, much less is known about the importance of MET in the later steps of metastatic colonization. We created novel reporters, which integrate transcriptional and post-transcriptional regulation, to test whether MET is required for metastasis in multiple in vivo cancer models. In a model of carcinosarcoma, metastasis occurred via an MET-dependent pathway; however, in two prostate carcinoma models, metastatic colonization was MET independent. Our results provide evidence for both MET-dependent and MET-independent metastatic pathways.

Modulation of alternative splicing with chemical compounds in new therapeutics for human diseases

Alternative splicing is a critical step where a limited number of human genes generate a complex and diverse proteome. Various diseases, including inherited diseases with abnormalities in the "genome code," have been found to result in an aberrant mis-spliced "transcript code" with correlation to the resulting phenotype. Chemical compound-based and nucleic acid-based strategies are trying to target this mis-spliced "transcript code". We will briefly mention about how to obtain splicing-modifying-compounds by high-throughput screening and overview of what is known about compounds that modify splicing pathways. The main focus will be on RNA-binding protein kinase inhibitors. In the main text, we will refer to diseases where splicing-modifying-compounds have been intensively investigated, with comparison to nucleic acid-based strategies. The information on their involvement in mis-splicing as well as nonsplicing events will be helpful in finding better compounds with less off-target effects for future implications in mis-splicing therapy.

Real-time imaging of the epithelial-mesenchymal transition using microRNA-200a sequence-based molecular beacon-conjugated magnetic nanoparticles

The epithelial-mesenchymal transition (EMT) plays important roles in tumor progression to metastasis. Thus, the development of an imaging probe that can monitor transient periods of the EMT process in live cells is required for a better understanding of metastatic process. Inspired by the fact that the mRNA expression levels of zinc finger E-box-binding homeobox 1 (ZEB1) increase when cells adopt mesenchyme characteristics and that microRNA-200a (miR-200a) can bind to ZEB1 mRNA, we conjugated molecular beacon (MB) mimicking mature miR-200a to magnetic nanoparticles (miR-200a-MB-MNPs) and devised an imaging method to observe transitional changes in the cells during EMT. Transforming growth factor-β1 treated epithelial cells and breast cancer cell lines representing both epithelial and mesenchymal phenotypes were used for the validation of miR-200a-MB-MNPs as an EMT imaging probe. The real-time imaging of live cells acquired with the induction of EMT revealed an increase in fluorescence signals by miR-200a-MB-MNPs, cell morphology alterations, and the loss of cell-cell adhesion. Our results suggest that miR-200a-MB-MNPs can be used as an imaging probe for the real-time monitoring of the EMT process in live cells.

Cellular migration and invasion uncoupled: increased migration is not an inexorable consequence of epithelial-to-mesenchymal transition

Metastatic dissemination requires carcinoma cells to detach from the primary tumor and invade through the basement membrane. To acquire these characteristics, epithelial tumor cells undergo epithelial-to-mesenchymal transitions (EMT), whereby cells lose polarity and E-cadherin-mediated cell-cell adhesion. Post-EMT cells have also been shown, or assumed, to be more migratory; however, there have been contradictory reports on an immortalized human mammary epithelial cell line (HMLE) that underwent EMT. In the context of carcinoma-associated EMT, it is not yet clear whether the change in migration and invasion must be positively correlated during EMT or whether enhanced migration is a necessary consequence of having undergone EMT. Here, we report that pre-EMT rat prostate cancer (PC) and HMLE cells are more migratory than their post-EMT counterparts. To determine a mechanism for increased epithelial cell migration, gene expression analysis was performed and revealed an increase in epidermal growth factor receptor (EGFR) expression in pre-EMT cells. Indeed, inhibition of EGFR in PC epithelial cells slowed migration. Importantly, while post-EMT PC and HMLE cell lines are less migratory, both remain invasive in vitro and, for PC cells, in vivo. Our study demonstrates that enhanced migration is not a phenotypic requirement of EMT, and migration and invasion can be uncoupled during carcinoma-associated EMT.

Identification of latent biomarkers in hepatocellular carcinoma by ultra-deep whole-transcriptome sequencing

There is an urgent need to identify biomarkers for hepatocellular carcinoma due to limited treatment options and the poor prognosis of this common lethal disease. Whole-transcriptome shotgun sequencing (RNA-Seq) provides new possibilities for biomarker identification. We sequenced ∼250 million pair-end reads from a pair of adjacent normal and tumor liver samples. With the aid of bioinformatics tools, we determined the transcriptome landscape and sought novel biomarkers by further empirical validations in 55 pairs of adjacent normal and tumor liver samples with various viral statuses such as HBV(+), HCV(+) and HBV(-)HCV(-). We identified a novel gene with coding regions, termed DUNQU1, which has a tissue-specific expression pattern in tumor liver samples of HCV(+) and HBV(-)HCV(-) hepatocellular carcinomas. Overexpression of DUNQU1 in Huh7 cell lines enhances the ability to form colonies in soft agar. Also, we identified three novel differentially-expressed protein-coding genes (ALG1L, SERPINA11 and TMEM82) that lack documented expression profiles in liver cancer and showed that the level of SREPINA11 is correlated with pathology stages. Moreover, we showed that the alternative splicing event of FGFR2 is associated with virus infection, tumor size, cirrhosis and tumor recurrence. The findings indicate that these new markers of hepatocellular carcinoma may be of value in improving prognosis and could have potential as new targets for developing new treatment options.