Isoform Inference From RNA-Seq Samples Based on Gene Structures on Chromosomes

Yan Ji, Jia Wei

R&D Information, AstraZeneca, Shanghai, China.

DOI: 10.4236/jbm.2013.11001                

Abstract

The emerging RNA-Seq technology makes it possible to infer splicing variants from millions of short sequence reads. Here we present a method to identify isoforms by their specific signatures on chromosomes including both exons and junctions. By applying this method to a RNA-Seq dataset of gastric cancer, we showed that our method is more accurate and sensitive than other isoform inference tools such as RSEM and Cufflinks. By constructing a network from gene list identified by our method but missed by other tools, we found that some cancer-related genes enriched in network modules have significant implications for cancer drug discovery.

Keywords

RNA-Seq; isoform inference; cancer genes; Cufflinks; RSEM

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Metzker, M.L. (2010) Sequencing technologies - the next generation. Nat Rev Genet, 11(1), 31-46.
[2]	Mortazavi, A. Williams, B. A. McCue, K. Schaeffer, L. Wold, B. (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods, 2008. 5(7): p. 621-8.
[3]	Wang, E.T., et al. (2008) Alternative isoform regulation in human tissue transcriptomes. Nature, 456(7221), 470-6.
[4]	Guttman, M., et al. (2010) Ab initio reconstruction of cell type-specific transcriptomes in mouse reveals the con-served multi-exonic structure of lincRNAs. Nat Biotechnol, 28(5): p. 503-10.
[5]	Li, B. and C.N. Dewey (2011) RSEM: accurate transcript quantification from RNA-Seq data with or without a ref-erence genome. BMC Bioinformatics, 12: p. 323.
[6]	Trapnell, C., et al. (2010) Transcript assembly and quan-tification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol, 28(5): p. 511-5.
[7]	Trapnell, C., L. Pachter, and S.L. Salzberg (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinfor-matics, 25(9): p. 1105-11.
[8]	Snowden, A.W., et al. (2003) Repression of vascular endothelial growth factor A in glioblastoma cells using engineered zinc finger transcription factors. Cancer Res, 63(24): p. 8968-76.
[9]	Ngan, E.S., et al. (2003) Overexpression of Cdc25B, an androgen receptor coactivator, in prostate cancer. Onco-gene, 22(5): p. 734-9.
[10]	Samanta, A.K., et al.(2004) Overexpression of MEKK3 confers resistance to apoptosis through activation of NFkappaB. J Biol Chem, 279(9): p. 7576-83.