Integrated mRNA and microRNA transcriptome sequencing characterizes sequence variants and mRNA-microRNA regulatory network in nasopharyngeal carcinoma model systems

Viral Encoded miRNAs in Tumorigenesis: Theranostic Opportunities in Precision Oncology

About 15% of all human cancers have a viral etiology. Although progress has been made, understanding the viral oncogenesis and associated molecular mechanisms remain complex. The discovery of cellular miRNAs has led to major breakthroughs. Interestingly, viruses have also been discovered to encode their own miRNAs. These viral, small, non-coding miRNAs are also known as viral-miRNAs (v-miRNAs). Although the function of v-miRNAs largely remains to be elucidated, their role in tumorigenesis cannot be ignored. V-miRNAs have also been shown to exploit the cellular machinery to benefit viral replication and survival. Although the discovery of Hepatitis C virus (HCV), and its viral miRNAs, is a work in progress, the existence of HPV-, EBV-, HBV-, MCPyV- and KSHV-encoded miRNA has been documented. V-miRNAs have been shown to target host factors to advance tumorigenesis, evade and suppress the immune system, and deregulate both the cell cycle and the apoptotic machinery. Although the exact mechanisms of v-miRNAs-induced tumorigenesis are still unclear, v-miRNAs are active role-players in tumorigenesis, viral latency and cell transformation. Furthermore, v-miRNAs can function as posttranscriptional gene regulators of both viral and host genes. Thus, it has been proposed that v-miRNAs may serve as diagnostic biomarkers and therapeutic targets for cancers with a viral etiology. Although significant challenges exist in their clinical application, emerging reports demonstrate their potent role in precision medicine. This review will focus on the roles of HPV-, HCV-, EBV-, HBV-, MCPyV-, and KSHV-produced v-miRNAs in tumorigenesis, as effectors in immune evasion, as diagnostic biomarkers and as novel anti-cancer therapeutic targets. Finally, it will discuss the challenges and opportunities associated with v-miRNAs theranostics in precision oncology.

General Features and Novel Gene Signatures That Identify Epstein-Barr Virus-Associated Epithelial Cancers

Epstein-Barr virus (EBV) is associated with various types of human malignancies, including nasopharyngeal carcinoma (NPC), EBV-associated gastric carcinoma (EBVaGC), and oral squamous cell carcinoma (OSCC). The present study aimed to identify gene signatures and common signaling pathways that can be used to predict the prognosis of EBV-associated epithelial cancers (EBVaCAs) by performing an integrated bioinformatics analysis of cell lines and tumor tissues. We identified 12 differentially expressed genes (DEGs) in the EBVaCA cell lines. Among them, only four DEGs, including BAMBI, SLC26A9, SGPP2, and TMC8, were significantly upregulated. However, SLC26A9 and TMC8, but not BAMBI and SGPP2, were significantly upregulated in EBV-positive tumor tissues compared to EBV-negative tumor tissues. Next, we identified IL6/JAK/STAT3 and TNF-α/NF-κB signaling pathways as common hallmarks of EBVaCAs. The expression of key genes related to the two hallmarks was upregulated in both EBV-infected cell lines and EBV-positive tumor tissues. These results suggest that SLC26A9 and TMC8 might be gene signatures that can effectively predict the prognosis of EBVaCAs and provide new insights into the molecular mechanisms of EBV-driven epithelial cancers.

Chemical capping improves template switching and enhances sequencing of small RNAs

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5′ cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5′ cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA—a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

ΔNp63α is a super enhancer-enriched master factor controlling the basal-to-luminal differentiation transcriptional program and gene regulatory networks in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) originates via malignant transformation of the pseudostratified nasopharyngeal epithelium, composed of basal and luminal cells. Super enhancers (SEs) are large clusters of cis-elements involved in the regulation of gene expression through epigenetic regulatory mechanisms. In this study, we demonstrated that basal cell-specific proteins are highly expressed, whereas luminal cell proteins are downregulated in NPC, implying a perturbation of basal-to-luminal differentiation during NPC development. We characterized NPC cell models according to different molecular signatures associated with their differentiation status and found that distinct SE landscapes are tightly associated with basal or luminal-like molecular signatures in NPC cells. Furthermore, the transcription of ΔNP63α, a prominent isoform of TP63, was found to be driven by SEs in NPC cells. Data from chromatin immunoprecipitation (ChIP)-sequencing showed that ΔNP63α largely occupied regions of SEs associated with basal cell-specific genes. Silencing of ΔNP63α led to a loss of H3K27ac occupancy at basal-type SEs and triggered a basal-to-luminal gene expression signature switch, suggesting that ΔNP63α is a master factor contributing to the perturbation of luminal differentiation. Integrative transcriptomics analysis also revealed that ΔNP63α acts as a core factor involved in the dysregulation of gene expression in NPC. Furthermore, ΔNP63α enhanced EGF-stimulated NF-κB activation in NPC cells by activating SE-mediated EGFR transcription. Finally, depletion of ΔNP63α in NPC cells induced robust growth inhibition of NPC cells in vitro and in vivo. Our data revealed that ΔNP63α-dependent SE reprogramming contributes to the blockade of luminal differentiation and uncontrolled proliferation in NPC.

Targeting the signaling in Epstein-Barr virus-associated diseases: mechanism, regulation, and clinical study

Epstein–Barr virus-associated diseases are important global health concerns. As a group I carcinogen, EBV accounts for 1.5% of human malignances, including both epithelial- and lymphatic-originated tumors. Moreover, EBV plays an etiological and pathogenic role in a number of non-neoplastic diseases, and is even involved in multiple autoimmune diseases (SADs). In this review, we summarize and discuss some recent exciting discoveries in EBV research area, which including DNA methylation alterations, metabolic reprogramming, the changes of mitochondria and ubiquitin-proteasome system (UPS), oxidative stress and EBV lytic reactivation, variations in non-coding RNA (ncRNA), radiochemotherapy and immunotherapy. Understanding and learning from this advancement will further confirm the far-reaching and future value of therapeutic strategies in EBV-associated diseases.

EBV-LMP1 induces APOBEC3s and mitochondrial DNA hypermutation in nasopharyngeal cancer

An Epstein-Barr virus (EBV)-encoded latent membrane protein 1 (LMP1) is a principal oncogene that plays a pivotal role in EBV-associated malignant tumors including nasopharyngeal cancer (NPC). Recent genomic landscape studies revealed that NPC also contained many genomic mutations, suggesting the role of LMP1 as a driver gene for the induction of these genomic mutations. Nonetheless, its exact mechanism has not been investigated. In this study, we report that LMP1 alters the expression profile of APOBEC3s(A3s), host deaminases that introduce consecutive C-to-U mutations (hypermutation). In vitro, LMP1 induces APOBEC3B (A3B) and 3F(A3F), in a nasopharyngeal cell line, AdAH. Overexpression of LMP1, A3B, or A3F induces mtDNA hypermutation, which is also detectable from NPC specimens. Expression of LMP1 and A3B in NPC was correlated with neck metastasis. These results provide evidence as to which LMP1 induces A3s and mtDNA hypermutation, and how LMP1 facilitates metastasis is also discussed.

Enrichment and identification of differentially expressed genes in hepatocellular carcinoma stem-like cells

Cancer stem cells are considered to be tumor-initiating cells. To explain the initiation or progression of hepatocellular carcinoma (HCC), we previously established a culture system that may enrich hepatic cancer stem-like cells (HCSCs). However, the regulatory mechanisms by which HCSCs acquire stem cell properties remain unclear. In the present study, three pairs of HCSCs and case-matched human HCC cells were analyzed by high-throughput screening, and novel biomarkers and pathways for the regulation of HCSCs were identified. The results led to the identification and stratification of 406 differentially expressed genes (DEGs), among which 73 GO terms were found to be significantly associated with DEGs in HCSCs, and only complement and coagulation cascade pathways were identified during the development of HCSCs. By combining the results of the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, it was revealed that 7 genes were downregulated in the complement and coagulation cascade pathways, and 7 miRNAs were predicted to target several downregulated genes involved in these pathways. The results may contribute toward hepatic cancer stem cell studies and novel drug research for HCC treatment.

Transcriptomics and solid tumors: The next frontier in precision cancer medicine

Transcriptomics, which encompasses assessments of alternative splicing and alternative polyadenylation, identification of fusion transcripts, explorations of noncoding RNAs, transcript annotation, and discovery of novel transcripts, is a valuable tool for understanding cancer mechanisms and identifying biomarkers. Recent advances in high-throughput technologies have enabled large-scale gene expression profiling. Importantly, RNA expression profiling of tumor tissue has been successfully used to determine clinically actionable molecular alterations. The WINTHER precision medicine clinical trial was the first prospective trial in diverse solid malignancies that assessed both genomics and transcriptomics to match treatments to specific molecular alterations. The use of transcriptome analysis in WINTHER and other trials increased the number of targetable -omic changes compared to genomic profiling alone. Other applications of transcriptomics involve the evaluation of tumor and circulating noncoding RNAs as predictive and prognostic biomarkers, the improvement of risk stratification by the use of prognostic and predictive multigene assays, the identification of fusion transcripts that drive tumors, and an improved understanding of the impact of DNA changes as some genomic alterations are silenced at the RNA level. Finally, RNA sequencing and gene expression analysis have been incorporated into clinical trials to identify markers predicting response to immunotherapy. Many issues regarding the complexity of the analysis, its reproducibility and variability, and the interpretation of the results still need to be addressed. The integration of transcriptomics with genomics, proteomics, epigenetics, and tumor immune profiling will improve biomarker discovery and our understanding of disease mechanisms and, thereby, accelerate the implementation of precision oncology.

Role of miR-96/EVI1/miR-449a Axis in the Nasopharyngeal Carcinoma Cell Migration and Tumor Sphere Formation

The Wnt signaling pathway is one of the major signaling pathways used by cancer stem cells (CSC). Ecotropic Viral Integration Site 1 (EVI1) has recently been shown to regulate oncogenic development of tumor cells by interacting with multiple signaling pathways, including the Wnt signaling. In the present study, we found that the Wnt modulator ICG-001 could inhibit the expression of EVI1 in nasopharyngeal carcinoma (NPC) cells. Results from loss-of-function and gain-of-function studies revealed that EVI1 expression positively regulated both NPC cell migration and growth of CSC-enriched tumor spheres. Subsequent studies indicated ICG-001 inhibited EVI1 expression via upregulated expression of miR-96. Results from EVI1 3′UTR luciferase reporter assay confirmed that EVI1 is a direct target of miR-96. Further mechanistic studies revealed that ICG-001, overexpression of miR-96, or knockdown of EVI1 expression could restore the expression of miR-449a. The suppressive effect of miR-449a on the cell migration and tumor sphere formation was confirmed in NPC cells. Taken together, the miR-96/EVI1/miR-449a axis is a novel pathway involved in ICG-001-mediated inhibition of NPC cell migration and growth of the tumor spheres.

Differential expression analysis comparing EBV uninfected to infected human cell lines identifies induced non-micro small non-coding RNAs

Epstein–Barr virus (EBV) is a ubiquitous human herpes virus, which is implicated in cancer and various autoimmune diseases. This study profiles non-micro small non-coding RNA expression changes induced by latent EBV infection. Using small RNA-Seq, 346 non-micro small RNAs were identified as being significantly differentially expressed between EBV(+) BJAB-B1 and EBV(−) BJAB cell lines. Select small RNA expression changes were experimentally validated in the BJAB-B1 cell line as well as the EBV-infected Raji and Jijoye cell lines. This latter analysis recapitulated the previously identified induction of vault RNA1, while also finding novel evidence for the deregulation of several tRNAs and a snoRNA.

Hsa_circ_0011385 accelerates the progression of thyroid cancer by targeting miR-361-3p

Background

Thyroid cancer is an endocrine malignancy that is growing in incidence worldwide. Despite progress in diagnostics and treatment of thyroid cancer, prognosis remains poor. Emerging research has shown that circular RNAs (circRNAs) have crucial regulatory roles in cancers. However, the possible functions and mechanisms of hsa_circ_0011385 remain undetermined.

Materials and methods

Expression levels of hsa_circ_0011385 and miR-361-3p were evaluated by qRT-PCR assay. The interaction between hsa_circ_0011385 and miR-361-3p was verified by dual-luciferase reporter assay. Effects of hsa_circ_0011385 or miR-361-3p on cell viability, proliferation, cell cycle, apoptosis, migration and invasion were confirmed by cell counting kit-8 (CCK-8), carboxyfluoresceinsuccinimidyl ester (CFSE), flow cytometry, and Transwell assays in vitro. The effect of hsa_circ_0011385 on thyroid cancer progression was also determined by in vivo tumor formation assay. Target genes of miR-361-3p were predicted by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and the expression of apoptosis- and metastasis-related proteins were assessed by Western blot assay.

Results

Hsa_circ_0011385 upregulated in thyroid cancer; hsa_circ_0011385 knockdown inhibited thyroid cancer cell proliferation, migration and invasion, and promoted cell cycle arrest and apoptosis. In addition, hsa_circ_0011385 could negatively regulate miR-361-3p by functioning as a sponge. Hsa_circ_0011385 promoted thyroid cancer cell proliferation, migration and invasion and suppressed cell cycle arrest and apoptosis by targeting miR-361-3p in vitro. We also found that hsa_circ_0011385 knockdown dramatically inhibited thyroid cancer growth in vivo. Furthermore, hsa_circ_0011385 regulated expression of apoptosis and metastasis-related proteins in thyroid cancer.

Conclusions

Hsa_circ_0011385facilitated thyroid cancer cell proliferation, invasion and migration, and inhibited thyroid cancer cell cycle arrest and apoptosis by targeting miR-361-3p, suggesting that the hsa_circ_0011385/miR-361-3p axis might be a promising therapeutic target for thyroid cancer.

Single cell RNA-seq reveals the landscape of tumor and infiltrating immune cells in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is one of the most malignant tumors in Southern China and southeast Asia, which is characterized by a dense lymphocyte infiltration and a poor prognosis. The emergence of single-cell sequencing represents a powerful tool to resolve tumor heterogeneity and delineate the complex communication among the tumor cells with neighboring stromal and immune cells in the tumor microenvironment (TME). Here, we performed single cell RNA-seq and analyzed tumor cells together with the infiltrating immune cells from three NPC tumor tissues. In our study, the malignant cells display the intra- and inter-tumoral heterogeneity among the individual patients. Analysis of the immune cells reveal the heterogeneous composition of the distinct immune cells and the various functional states of T cells in NPC tumors. Additionally, coupled with the reconstruct of the T cell receptor (TCR) sequences from immune cells full-length single-cell sequence data, we identify the diverse T cell clonotypes and expansion distribution in individual tumors. Overall, we firstly reveal the landscape of tumor and infiltrating immune cells in nasopharyngeal cancer. These results provide deeper insights on the mechanisms of tumor clearance by immune cells in the surrounding microenvironment, which will be helpful in improving the targeted and immune therapies for NPC.

Transcriptomic analysis of microRNAs-mRNAs regulating innate immune response of zebrafish larvae against Vibrio parahaemolyticus infection

In recent years, microRNAs (miRNAs) have been shown to play important roles in immunity. Analyses of the functions of miRNAs and their targets are useful in understanding the regulation of the immune response. To understand the relationships between miRNAs and their targets during infection, we used zebrafish as an infection model in which to characterize the miRNA and mRNA transcriptomes of zebrafish larvae infected with Vibrio parahaemolyticus. We identified the differentially expressed miRNAs and mRNAs. Overall, 37 known zebrafish miRNAs were differentially expressed in the infection group and 107 predicted target genes of 26 miRNAs were differentially expressed in the mRNA transcriptome. These targets with specific Gene Ontology (GO) terms, such as peripheral nervous system neuron axonogenesis, organophosphate metabolic process, heme binding, protein binding, tetrapyrrole binding, protein dimerization activity, and aromatase activity, which regulate nerve conduction, energy metabolism, hematopoiesis, and protein synthesis. They were also associated with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways such as phototransduction, tryptophan metabolism, notch signaling, and purine metabolism. Our findings indicate that miRNAs regulate the innate immune response via complex networks, and zebrafish (Danio rerio, dre)-miR-205-3p, dre-miR-141-5p, dre-miR-200a-5p, dre-miR-92a-2-5p, dre-miR-192, and dre-miR-1788 may play important roles in the innate immune response by regulating target genes.

The emerging role of Epstein-Barr virus encoded microRNAs in nasopharyngeal carcinoma

Epstein-Barr virus (EBV) is an oncogenic herpes virus that is closely associated with the initiation and development of nasopharyngeal carcinoma (NPC), lymphoma and other malignant tumors. EBV encodes 44 mature miRNAs that regulate viral and host cell gene expression and plays a variety of roles in biological functions and the development of cancer. In this review, we summarized the biological functions and molecular mechanisms of Epstein-Barr virus-encoded microRNAs (EBV miRNAs) in tumor immune evasion, proliferation, anti-apoptosis, invasion, metastasis and as a potential biomarker for NPC diagnosis and prognosis. The knowledge generated by EBV miRNAs can be used for EBV miRNA-based precision cancer treatments in the near future.

MicroRNA-Related Polymorphisms in Infectious Diseases-Tiny Changes With a Huge Impact on Viral Infections and Potential Clinical Applications

MicroRNAs (miRNAs) are single-stranded sequences of non-coding RNA with approximately 22 nucleotides that act posttranscriptionally on gene expression. miRNAs are important gene regulators in physiological contexts, but they also impact the pathogenesis of various diseases. The role of miRNAs in viral infections has been explored by different authors in both population-based as well as in functional studies. However, the effect of miRNA polymorphisms on the susceptibility to viral infections and on the clinical course of these diseases is still an emerging topic. Thus, this review will compile and organize the findings described in studies that evaluated the effects of genetic variations on miRNA genes and on their binding sites, in the context of human viral diseases. In addition to discussing the basic aspects of miRNAs biology, we will cover the studies that investigated miRNA polymorphisms in infections caused by hepatitis B virus, hepatitis C virus, human immunodeficiency virus, Epstein–Barr virus, and human papillomavirus. Finally, emerging topics concerning the importance of miRNA genetic variants will be presented, focusing on the context of viral infectious diseases.

Epigenetic Modifications as Biomarkers of Tumor Development, Therapy Response, and Recurrence across the Cancer Care Continuum

Aberrant epigenetic modifications are an early event in carcinogenesis, with the epigenetic landscape continuing to change during tumor progression and metastasis—these observations suggest that specific epigenetic modifications could be used as diagnostic and prognostic biomarkers for many cancer types. DNA methylation, post-translational histone modifications, and non-coding RNAs are all dysregulated in cancer and are detectable to various degrees in liquid biopsies such as sputum, urine, stool, and blood. Here, we will focus on the application of liquid biopsies, as opposed to tissue biopsies, because of their potential as non-invasive diagnostic tools and possible use in monitoring therapy response and progression to metastatic disease. This includes a discussion of septin-9 (SEPT9) DNA hypermethylation for detecting colorectal cancer, which is by far the most developed epigenetic biomarker assay. Despite their potential as prognostic and diagnostic biomarkers, technical issues such as inconsistent methodology between studies, overall low yield of epigenetic material in samples, and the need for improved histone and non-coding RNA purification methods are limiting the use of epigenetic biomarkers. Once these technical limitations are overcome, epigenetic biomarkers could be used to monitor cancer development, disease progression, therapeutic response, and recurrence across the entire cancer care continuum.

Role of Epstein-Barr Virus in the Development of Nasopharyngeal Carcinoma

Southern China experiences larger extent of total cancer pathologies, of which nasopharyngeal carcinoma has the highest incidence under otorhinolaryngeal malignant carcinomas. Risk factor of nasopharyngeal carcinoma varies from hereditary causes to virus infection, among which Epstein-Barr virus (EBV) infection is the mostly investigated. The study into mechanism of EBV in occurrence, development and prognosis of nasopharyngeal carcinoma has been studied for several decades. The pathophysiology in making of EBV into a cancerogen includes proteins as latent membrane protein 1 (LMPs) and nucleic acids as micro-RNAs. In this paper, we reviewed till date studies focusing on relationship between EBV and nasopharyngeal carcinoma.

Large-scale SNP screenings identify markers linked with GCRV resistant traits through transcriptomes of individuals and cell lines in Ctenopharyngodon idella

Grass carp (Ctenopharyngodon idella) is an important economic species in freshwater aquaculture and its industry has been confined due to variety degeneration and frequent diseases. Marker-assisted selection is a feasible method for selective breeding of new varieties. Transcriptome data have greatly facilitated high-throughput single nucleotide polymorphism (SNP) marker discovery and phenotype association study. In this study, we gained a total of 25,981 and 5,775 high quality SNPs in two transcriptomes from individuals and cell lines, respectively. Comparative transcriptome analysis identified 413 and 832 grass carp reovirus (GCRV)-resistant-association SNPs as well as 1,381 and 1,606 GCRV-susceptible-association SNPs in individuals and cell lines, respectively. Integrated analysis indicated 22 genes with single SNP share common resistant/susceptible traits in two transcriptomes. Furthermore, we infected grass carp with GCRV, genotyping and association analyses were performed, and 9 in 22 SNPs were confirmed by PCR-RFLP. Meanwhile, mRNA expression profiles of 6 genes containing confirmed SNPs were examined by qRT-PCR. The results demonstrated that mRNA expressions were significant differences in resistant/susceptible individuals and cell lines. The present study develops an important strategy for high throughput screening of phenotype association genetic markers and the results will serve in grass carp breeding for GCRV resistance.

Identification of miRNA-mRNA Modules in Colorectal Cancer Using Rough Hypercuboid Based Supervised Clustering

Differences in the expression profiles of miRNAs and mRNAs have been reported in colorectal cancer. Nevertheless, information on important miRNA-mRNA regulatory modules in colorectal cancer is still lacking. In this regard, this study presents an application of the RH-SAC algorithm on miRNA and mRNA expression data for identification of potential miRNA-mRNA modules. First, a set of miRNA rules was generated using the RH-SAC algorithm. The mRNA targets of the selected miRNAs were identified using the miRTarBase database. Next, the expression values of target mRNAs were used to generate mRNA rules using the RH-SAC. Then all miRNA-mRNA rules have been integrated for generating networks. The RH-SAC algorithm unlike other existing methods selects a group of co-expressed miRNAs and mRNAs that are also differentially expressed. In total 17 miRNAs and 141 mRNAs were selected. The enrichment analysis of selected mRNAs revealed that our method selected mRNAs that are significantly associated with colorectal cancer. We identified novel miRNA/mRNA interactions in colorectal cancer. Through experiment, we could confirm that one of our discovered miRNAs, hsa-miR-93-5p, was significantly up-regulated in 75.8% CRC in comparison to their corresponding non-tumor samples. It could have the potential to examine colorectal cancer subtype specific unique miRNA/mRNA interactions.

LncRNAs: Bridging environmental sensing and gene expression

The survival of all organisms is dependent on complex, coordinated responses to environmental cues. Non-coding RNAs have been identified as major players in regulation of gene expression, with recent evidence supporting roles for long non-coding (lnc)RNAs in both transcriptional and post-transcriptional control. Evidence from our laboratory shows that lncRNAs have the ability to form hybridized structures called R-loops with specific DNA target sequences in S. cerevisiae, thereby modulating gene expression. In this Point of View, we provide an overview of the nature of lncRNA-mediated control of gene expression in the context of our studies using the GAL gene cluster as a model for controlling the timing of transcription.

SPINK6 Promotes Metastasis of Nasopharyngeal Carcinoma via Binding and Activation of Epithelial Growth Factor Receptor

Nasopharyngeal carcinoma has the highest rate of metastasis among head and neck cancers, and distant metastasis is the major reason for treatment failure. The underlying molecular mechanisms of nasopharyngeal carcinoma metastasis are not fully understood. Here, we report the identification of serine protease inhibitor Kazal-type 6 (SPINK6) as a functional regulator of nasopharyngeal carcinoma metastasis via EGFR signaling. SPINK6 mRNA was upregulated in tumor and highly metastatic nasopharyngeal carcinoma cells. Immunohistochemical staining of 534 nasopharyngeal carcinomas revealed elevated SPINK6 expression as an independent unfavorable prognostic factor for overall, disease-free, and distant metastasis-free survival. Ectopic SPINK6 expression promoted in vitro migration and invasion as well as in vivo lymph node metastasis and liver metastasis of nasopharyngeal carcinoma cells, whereas silencing SPINK6 exhibited opposing effects. SPINK6 enhanced epithelial-mesenchymal transition by activating EGFR and the downstream AKT pathway. Inhibition of EGFR with a neutralizing antibody or erlotinib reversed SPINK6-induced nasopharyngeal carcinoma cell migration and invasion. Erlotinib also inhibited SPINK6-induced metastasis in vivo Notably, SPINK6 bound to the EGFR extracellular domain independent of serine protease-inhibitory activity. Overall, our results identified a novel EGFR-activating mechanism in which SPINK6 has a critical role in promoting nasopharyngeal carcinoma metastasis, with possible implications as a prognostic indicator in nasopharyngeal carcinoma patients. Cancer Res; 77(2); 579-89. ©2016 AACR.

A Multi-Step miRNA-mRNA Regulatory Network Construction Approach Identifies Gene Signatures Associated with Endometrioid Endometrial Carcinoma

We aimed to identify endometrioid endometrial carcinoma (EEC)-related gene signatures using a multi-step miRNA-mRNA regulatory network construction approach. Pathway analysis showed that 61 genes were enriched on many carcinoma-related pathways. Among the 14 highest scoring gene signatures, six genes had been previously shown to be endometrial carcinoma. By qRT-PCR and next generation sequencing, we found that a gene signature (CPEB1) was significantly down-regulated in EEC tissues, which may be caused by hsa-miR-183-5p up-regulation. In addition, our literature surveys suggested that CPEB1 may play an important role in EEC pathogenesis by regulating the EMT/p53 pathway. The miRNA-mRNA network is worthy of further investigation with respect to the regulatory mechanisms of miRNAs in EEC. CPEB1 appeared to be a tumor suppressor in EEC. Our results provided valuable guidance for the functional study at the cellular level, as well as the EEC mouse models.

Valuable lessons-learned in transcriptomics experimentation

We have collected several valuable lessons that will help improve transcriptomics experimentation. These lessons relate to experiment design, execution, and analysis. The cautions, but also the pointers, may help biologists avoid common pitfalls in transcriptomics experimentation and achieve better results with their transcriptome studies.

Frequent BCOR aberrations in extranodal NK/T-Cell lymphoma, nasal type

Extranodal natural killer/T cell lymphoma (ENKTL) is a rare subtype of lymphoma. Recurrent mutations in the JAK-STAT pathway, recently reported in ENKTL cases, are interesting in terms of both pathogenesis and inhibitor therapy. However, the frequencies of these mutations are low and variable among reports, and other pathognomonic mutations in ENKTL remain to be elucidated. In the present study, targeted capture sequencing of 602 cancer-related genes from 25 frozen ENKTL samples was performed, 11 of which were matched to normal samples. Several recurrent somatic mutations involving BCOR (32%), TP53 (16%), DDX3X (12%), FAT4 (8%), NRAS (8%), MLL3 (12%), and MIR17HG (8%) were identified. The pattern of BCOR aberrations (1 nonsense and 5 frame-shift mutations, a mutation leading to a splicing error, and gene loss) suggested that loss of function of BCOR was the functionally important outcome of such changes. The literature was reviewed and the public data on BCOR aberrations was reanalyzed and it was found that the aberrations were frequently found in myeloid neoplasms, but, interestingly, were highly specific to ENKTL among lymphoid malignancies. Given the high frequency and pattern of aberration, BCOR is likely to play an important role in ENKTL pathogenesis as a tumor suppressor gene.

Targeted deletion of miR-132/-212 impairs memory and alters the hippocampal transcriptome

miR-132 and miR-212 are structurally related microRNAs that have been found to exert powerful modulatory effects within the central nervous system (CNS). Notably, these microRNAs are tandomly processed from the same noncoding transcript, and share a common seed sequence: thus it has been difficult to assess the distinct contribution of each microRNA to gene expression within the CNS. Here, we employed a combination of conditional knockout and transgenic mouse models to examine the contribution of the miR-132/-212 gene locus to learning and memory, and then to assess the distinct effects that each microRNA has on hippocampal gene expression. Using a conditional deletion approach, we show that miR-132/-212 double-knockout mice exhibit significant cognitive deficits in spatial memory, recognition memory, and in tests of novel object recognition. Next, we utilized transgenic miR-132 and miR-212 overexpression mouse lines and the miR-132/-212 double-knockout line to explore the distinct effects of these two miRNAs on the transcriptional profile of the hippocampus. Illumina sequencing revealed that miR-132/-212 deletion increased the expression of 1138 genes; Venn analysis showed that 96 of these genes were also downregulated in mice overexpressing miR-132. Of the 58 genes that were decreased in animals overexpressing miR-212, only four of them were also increased in the knockout line. Functional gene ontology analysis of downregulated genes revealed significant enrichment of genes related to synaptic transmission, neuronal proliferation, and morphogenesis, processes known for their roles in learning, and memory formation. These data, coupled with previous studies, firmly establish a role for the miR-132/-212 gene locus as a key regulator of cognitive capacity. Further, although miR-132 and miR-212 share a seed sequence, these data indicate that these miRNAs do not exhibit strongly overlapping mRNA targeting profiles, thus indicating that these two genes may function in a complex, nonredundant manner to shape the transcriptional profile of the CNS. The dysregulation of miR-132/-212 expression could contribute to signaling mechanisms that are involved in an array of cognitive disorders.

The interplay of host genetic factors and Epstein-Barr virus in the development of nasopharyngeal carcinoma

The interplay between host cell genetics and Epstein-Barr virus (EBV) infection contributes to the development of nasopharyngeal carcinoma (NPC). Understanding the host genetic and epigenetic alterations and the influence of EBV on cell signaling and host gene regulation will aid in understanding the molecular pathogenesis of NPC and provide useful biomarkers and targets for diagnosis and therapy. In this review, we provide an update of the oncogenes and tumor suppressor genes associated with NPC, as well as genes associated with NPC risk including those involved in carcinogen detoxification and DNA repair. We also describe the importance of host genetics that govern the human leukocyte antigen (HLA) complex and immune responses, and we describe the impact of EBV infection on host cell signaling changes and epigenetic regulation of gene expression. High-power genomic sequencing approaches are needed to elucidate the genetic basis for inherited susceptibility to NPC and to identify the genes and pathways driving its molecular pathogenesis.

Host Gene Expression Is Regulated by Two Types of Noncoding RNAs Transcribed from the Epstein-Barr Virus BamHI A Rightward Transcript Region

In Epstein-Barr virus-infected epithelial cancers, the alternatively spliced BamHI A rightward transcripts (BARTs) are the most abundant viral polyadenylated RNA. The BART introns form the template for the production of 44 microRNAs (miRNAs), and the spliced and polyadenylated exons form nuclear non-protein-coding RNAs. Analysis of host cell transcription by RNA-seq during latency in AGS cells identified a large number of reproducibly changed genes. Genes that were downregulated were enriched for BART miRNA targets. Bioinformatics analysis predicted activation of the myc pathway and downregulation of XBP1 as likely mediators of the host transcriptional changes. Effects on XBP1 activity were not detected in these cells; however, myc activation was confirmed through use of a myc-responsive luciferase reporter. To identify potential regulatory properties of the spliced, polyadenylated BART RNAs, a full-length cDNA clone of one of the BART isoforms was obtained and expressed in the Epstein-Barr virus (EBV)-negative AGS cells. The BART cDNA transcript remained primarily nuclear yet induced considerable and consistent changes in cellular transcription, as profiled by RNA-seq. These transcriptional changes significantly overlapped the transcriptional changes induced during latent EBV infection of these same cells, where the BARTs are exclusively nuclear and do not encode proteins. These data suggest that the nuclear BART RNAs are functional long noncoding RNAs (lncRNAs). The abundant expression of multiple forms of noncoding RNAs that contribute to growth regulation without expression of immunogenic proteins would be an important mechanism for viral oncogenesis in the presence of a functional immune system.

IMPORTANCE

Infection with Epstein-Barr virus (EBV) is nearly ubiquitous in the human population; however, it does contribute to the formation of multiple types of cancer. In immunocompromised patients, EBV causes multiple types of lymphomas by expressing viral oncogenes that promote growth and survival of infected B lymphocytes. EBV-positive gastric carcinoma does not require immune suppression, and the viral oncoproteins that are frequent targets for an immunological response are not expressed. This study demonstrates using transcriptional analysis that the expression of various classes of viral non-protein-coding RNAs likely contribute to the considerable changes in the host transcriptional profile in the AGS gastric cancer cell line. This is the first report to show that the highly expressed polyadenylated BamHI A rightward transcripts (BART) viral transcript in gastric carcinoma is in fact a functional viral long noncoding RNA. These studies provide new insight into how EBV can promote transformation in the absence of viral protein expression.

Integrated analysis of miRNA, gene, and pathway regulatory networks in hepatic cancer stem cells

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide. HCC has a poor prognosis associated with tumor recurrence and drug resistance, which has been attributed to the existence of hepatic cancer stem cells (HCSCs). However, the characteristics and regulatory mechanisms of HCSCs remain unclear. We therefore established a novel system to enrich HCSCs and we demonstrate that these HCSCs exhibit cancer stem cell properties.

METHODS

We used miRNA and mRNA high-throughput sequencing data sets to determine molecular signatures and regulatory mechanisms in HCSCs. Paired miRNA and gene deep sequencing data in HCSCs versus HCC cells were used to identify candidate biomarkers of HCSCs. Using network analysis, we studied the relationship between miRNA and gene biomarkers, and KEGG pathway enrichment analysis was performed to study the function of candidate biomarkers.

RESULTS

We identified 9 up- and 9 down-regulated miRNAs and 115 up- and 402 down-regulated genes in HCSCs compared with HCC cells. A miRNA-gene network was constructed using 651 miRNA-gene interactions (between 7 up-regulated miRNAs and 274 down-regulated genes), and 103 miRNA-gene interactions (between 9 down-regulated miRNAs and 62 up-regulated genes). Pathway enrichment analysis identified five tumor invasion- and metastasis-related pathways and MAPK signaling associated with HCSCs. We further discovered two novel pathways that likely play a role in the regulation of HCSCs.

CONCLUSIONS

We identified a molecular expression signature and pathway regulatory mechanisms in HCSCs with potential diagnostic and therapeutic value.

Omics-based identification of biomarkers for nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a head and neck cancer that is highly found in distinct geographic areas, such as Southeast Asia. The management of NPC remains burdensome as the prognosis is poor due to the late presentation of the disease and the complex nature of NPC pathogenesis. Therefore, it is necessary to find effective molecular markers for early detection and therapeutic measure of NPC. In this paper, the discovery of molecular biomarker for NPC through the emerging omics technologies including genomics, miRNA-omics, transcriptomics, proteomics, and metabolomics will be extensively reviewed. These markers have been shown to play roles in various cellular pathways in NPC progression. The knowledge on their function will help us understand in more detail the complexity in tumor biology, leading to the better strategies for early detection, outcome prediction, detection of disease recurrence, and therapeutic approach.

Candidate pathways and genes for nasopharyngeal carcinoma based on bioinformatics study

PURPOSE

To reveal the potential microRNAs (miRNAs), genes, pathways and regulatory network involved in the process of nasopharyngeal carcinoma (NPC) by using the method of bioinformatics.

METHODS

Gene expression profiles GSE12452 (31 NPC and 10 normal samples) and GSE53819 (18 NPC and 18 normal samples), as well as miRNA expression profiles GSE32960 (312 NPC and 18 normal samples) and GSE36682 (62 NPC and 6 normal samples) were obtained from Gene Expression Omnibus database. The differentially expressed genes (DEGs) and miRNAs (DEmiRNAs) between NPC and normal samples were identified by using t-test based on MATLAB software (FDR < 0.01), followed by pathway enrichment analysis based on DAVID software (P-value < 0.1). Then, DEmiRNA-DEG regulatory network was constructed.

RESULTS

A total of 1254 DEGs and 107 DEmiRNAs were identified, respectively. Then, 16 pathways (including cell cycle) and 32 pathways (including pathways in cancer) were enriched by DEGs and target genes of DEmiRNAs, respectively. Furthermore, DEmiRNA-DEG regulatory network was constructed, containing 12 DEmiRNAs (including has-miR-615-3P) and 180 DEGs (including MCM4 and CCNE2).

CONCLUSION

has-miR-615-3p might take part in the pathogenetic process of NPC through regulating MCM4 which is enriched in cell cycle. The DEmiRNAs identified in the present study might serve as new biomarkers for NPC.

Making an impact: FEBS Open Bio celebrates its 3rd anniversary

As FEBS celebrates its 50th anniversary, it also marks the launch 3 years ago of its newest journal, FEBS Open Bio. This online-only open access journal is now indexed in Web of Science. Papers originally submitted to The FEBS Journal can be transferred to FEBS Open Bio together with any reviewer reports, offering authors a fast-track alternative for publication of their work.

Current transcriptomics in pig immunity research

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.