Display options
Share it on

BioData Min. 2016 Apr 05;9:13. doi: 10.1186/s13040-016-0090-8. eCollection 2016.

Distinguishing highly similar gene isoforms with a clustering-based bioinformatics analysis of PacBio single-molecule long reads.

BioData mining

Ma Liang, Castle Raley, Xin Zheng, Geetha Kutty, Emile Gogineni, Brad T Sherman, Qiang Sun, Xiongfong Chen, Thomas Skelly, Kristine Jones, Robert Stephens, Bin Zhou, William Lau, Calvin Johnson, Tomozumi Imamichi, Minkang Jiang, Robin Dewar, Richard A Lempicki, Bao Tran, Joseph A Kovacs, Da Wei Huang

Affiliations

  1. Critical Care Medicine Department, Clinical Center, Frederick, MD USA.
  2. Leidos BioMedical Research, Inc., Frederick National Laboratory for Cancer Research, NIH, Frederick, MD USA.
  3. Center of Information Technology, National Institutes of Health (NIH), Bethesda, MD USA.
  4. Current Affiliation: National Cancer Institute, NIH, Bethesda, MD USA.

PMID: 27051465 PMCID: PMC4820869 DOI: 10.1186/s13040-016-0090-8

Abstract

BACKGROUND: Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs.

RESULTS: We present a clustering-based analysis pipeline integrated with PacBio sequencing data for profiling highly similar gene isoforms. This approach was first evaluated in comparison to de novo assembly of 454 reads using a benchmark admixture containing 10 known, cloned msg genes encoding the major surface glycoprotein of Pneumocystis jirovecii. All 10 msg isoforms were successfully reconstructed with the expected length (~1.5 kb) and correct sequence by the new approach, while 454 reads could not be correctly assembled using various assembly programs. When using an additional benchmark admixture containing 22 known P. jirovecii msg isoforms, this approach accurately reconstructed all but 4 these isoforms in their full-length (~3 kb); these 4 isoforms were present in low concentrations in the admixture. Finally, when applied to the original clinical sample from which the 22 known msg isoforms were cloned, this approach successfully identified not only all known isoforms accurately (~3 kb each) but also 48 novel isoforms.

CONCLUSIONS: PacBio sequencing integrated with the clustering-based analysis pipeline achieves high-throughput and high-resolution discrimination of highly similar sequences, and can serve as a new approach for genome-wide characterization of gene isoforms and other highly repetitive sequences.

Keywords: Bioinformatics analysis; Gene isoforms; Major surface glycoprotein; NGS; PacBio; Pneumocystis; Repetitive sequences; Uclust

References

  1. JAMA. 2009 Jun 24;301(24):2578-85 - PubMed
  2. BMC Genomics. 2012 Sep 13;13:475 - PubMed
  3. MBio. 2012 Dec 26;4(1):e00428-12 - PubMed
  4. Genome Biol. 2013 Jul 03;14(7):405 - PubMed
  5. BMC Bioinformatics. 2004 Aug 19;5:113 - PubMed
  6. Nat Rev Genet. 2010 Jan;11(1):31-46 - PubMed
  7. Nat Rev Genet. 2011 Nov 29;13(1):36-46 - PubMed
  8. J Eukaryot Microbiol. 2007 Jan-Feb;54(1):8-13 - PubMed
  9. Infect Immun. 1998 Sep;66(9):4268-73 - PubMed
  10. BMC Genomics. 2009 Aug 07;10:367 - PubMed
  11. Genome Biol. 2011 Aug 11;12(8):R73 - PubMed
  12. Microb Inform Exp. 2014 Jan 15;4(1):1 - PubMed
  13. J Infect Dis. 2008 Sep 1;198(5):741-9 - PubMed
  14. Nature. 2002 Oct 3;419(6906):498-511 - PubMed
  15. Science. 2009 Jan 2;323(5910):133-8 - PubMed
  16. Clin Infect Dis. 2012 May;54(10):1437-44 - PubMed
  17. Genetics. 2005 Aug;170(4):1589-600 - PubMed
  18. Sci Rep. 2013 Oct 03;3:2837 - PubMed
  19. Brief Bioinform. 2014 May;15(3):431-42 - PubMed
  20. Nat Rev Microbiol. 2009 Jul;7(7):493-503 - PubMed
  21. Biochim Biophys Acta. 2014 Oct;1842(10):1932-1941 - PubMed
  22. Nature. 2008 Nov 27;456(7221):470-6 - PubMed
  23. Nat Commun. 2016 Feb 22;7:10740 - PubMed
  24. Bioinformatics. 2010 Oct 1;26(19):2460-1 - PubMed
  25. BMC Genomics. 2014 Feb 04;15:96 - PubMed
  26. Nat Biotechnol. 2011 Nov 13;30(1):99-104 - PubMed
  27. Nat Genet. 2008 Dec;40(12):1413-5 - PubMed
  28. J Eukaryot Microbiol. 2006;53 Suppl 1:S89-91 - PubMed
  29. J Biomed Biotechnol. 2010;2010:753675 - PubMed
  30. Science. 2007 Sep 28;317(5846):1921-6 - PubMed
  31. Nat Biotechnol. 2013 Nov;31(11):1009-14 - PubMed
  32. J Data Mining Genomics Proteomics. 2013 Jul 31;4(3):null - PubMed
  33. Genomics. 2010 Jun;95(6):315-27 - PubMed
  34. Science. 2005 Jul 15;309(5733):416-22 - PubMed

Publication Types

Grant support