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

Abstract

Nasopharyngeal carcinoma (NPC) is a prevalent malignancy in Southeast Asia among the Chinese population. Aberrant regulation of transcripts has been implicated in many types of cancers including NPC. Herein, we characterized mRNA and miRNA transcriptomes by RNA sequencing (RNASeq) of NPC model systems. Matched total mRNA and small RNA of undifferentiated Epstein–Barr virus (EBV)-positive NPC xenograft X666 and its derived cell line C666, well-differentiated NPC cell line HK1, and the immortalized nasopharyngeal epithelial cell line NP460 were sequenced by Solexa technology. We found 2812 genes and 149 miRNAs (human and EBV) to be differentially expressed in NP460, HK1, C666 and X666 with RNASeq; 533 miRNA–mRNA target pairs were inversely regulated in the three NPC cell lines compared to NP460. Integrated mRNA/miRNA expression profiling and pathway analysis show extracellular matrix organization, Beta-1 integrin cell surface interactions, and the PI3K/AKT, EGFR, ErbB, and Wnt pathways were potentially deregulated in NPC. Real-time quantitative PCR was performed on selected mRNA/miRNAs in order to validate their expression. Transcript sequence variants such as short insertions and deletions (INDEL), single nucleotide variant (SNV), and isomiRs were characterized in the NPC model systems. A novel TP53 transcript variant was identified in NP460, HK1, and C666. Detection of three previously reported novel EBV-encoded BART miRNAs and their isomiRs were also observed. Meta-analysis of a model system to a clinical system aids the choice of different cell lines in NPC studies. This comprehensive characterization of mRNA and miRNA transcriptomes in NPC cell lines and the xenograft provides insights on miRNA regulation of mRNA and valuable resources on transcript variation and regulation in NPC, which are potentially useful for mechanistic and preclinical studies.

Highlights

• •

  Using RNASeq we characterized the mRNA and miRNA transcriptomes in NPC and NP models.

• •

  2812 Genes and 149 miRNAs (human and EBV) were differentially expressed in NPC vs NP models.

• •

  533 miRNA–mRNA target pairs were inversely regulated in HK1, C666, and X666 vs NP460.

• •

  ECM, β1 integrin, PI3K/AKT, EGFR, ErbB, and Wnt pathways appeared to be deregulated in NPC.

• •

  A novel TP53 mutation was identified in NP460, HK1, and C666.