title: puma 3.0: improved uncertainty propagation methods for gene and transcript expression analysis.
abstract: BACKGROUND
  Microarrays have been a popular tool for gene expression profiling at genome-scale for over a decade due to the low cost, short turn-around time, excellent quantitative accuracy and ease of data generation. The Bioconductor package puma incorporates a suite of analysis methods for determining uncertainties from Affymetrix GeneChip data and propagating these uncertainties to downstream analysis. As isoform level expression profiling receives more and more interest within genomics in recent years, exon microarray technology offers an important tool to quantify expression level of the majority of exons and enables the possibility of measuring isoform level expression. However, puma does not include methods for the analysis of exon array data. Moreover, the current expression summarisation method for Affymetrix 3' GeneChip data suffers from instability for low expression genes. For the downstream analysis, the method for differential expression detection is computationally intensive and the original expression clustering method does not consider the variance across the replicated technical and biological measurements. It is therefore necessary to develop improved uncertainty propagation methods for gene and transcript expression analysis.
  
  
  RESULTS
  We extend the previously developed Bioconductor package puma with a new method especially designed for GeneChip Exon arrays and a set of improved downstream approaches. The improvements include: (i) a new gamma model for exon arrays which calculates isoform and gene expression measurements and a level of uncertainty associated with the estimates, using the multi-mappings between probes, isoforms and genes, (ii) a variant of the existing approach for the probe-level analysis of Affymetrix 3' GeneChip data to produce more stable gene expression estimates, (iii) an improved method for detecting differential expression which is computationally more efficient than the existing approach in the package and (iv) an improved method for robust model-based clustering of gene expression, which takes technical and biological replicate information into consideration.
  
  
  CONCLUSIONS
  With the extensions and improvements, the puma package is now applicable to the analysis of both Affymetrix 3' GeneChips and Exon arrays for gene and isoform expression estimation. It propagates the uncertainty of expression measurements into more efficient and comprehensive downstream analysis at both gene and isoform level. Downstream methods are also applicable to other expression quantification platforms, such as RNA-Seq, when uncertainty information is available from expression measurements. puma is available through Bioconductor and can be found at http://www.bioconductor.org.
authors: College of Computer Science and Technology, Nanjing University of Aeronautics and Astronautics, 29 Yudao St., Nanjing 210016, China. xuejun.liu@nuaa.edu.cn
authors: Xuejun
authors: X
authors: Liu
authors: Zhenzhu
authors: Z
authors: Gao
authors: Li
authors: L
authors: Zhang
authors: Magnus
authors: M
authors: Rattray
mesh_terms: D002294:Carcinoma, Squamous Cell / Q000235:genetics / Q000378:metabolism; D016000:Cluster Analysis; D005091:Exons; D020869:Gene Expression Profiling / Q000379:methods*; D023281:Genomics / Q000379:methods; D006258:Head and Neck Neoplasms / Q000235:genetics / Q000378:metabolism; D006801:Humans; D020411:Oligonucleotide Array Sequence Analysis; D059368:RNA Isoforms / Q000378:metabolism; D012984:Software; D000077195:Squamous Cell Carcinoma of Head and Neck; D035501:Uncertainty
chemical_list: D059368:RNA Isoforms
references: Nucleic Acids Res. 2011 Oct;39(18):e123
references: 21745820
references: Stat Appl Genet Mol Biol. 2010;9:Article42
references: 21194414
references: BMC Bioinformatics. 2010 Feb 18;11:94
references: 20167110
references: Biostatistics. 2003 Apr;4(2):249-64
references: 12925520
references: Bioinformatics. 2008 Aug 1;24(15):1707-14
references: 18573797
references: Nat Methods. 2010 Dec;7(12):1009-15
references: 21057496
references: Bioinformatics. 2011 Jul 1;27(13):i383-91
references: 21685096
references: Bioinformatics. 2005 Oct 1;21(19):3748-54
references: 16091409
references: Nat Biotechnol. 2006 Sep;24(9):1151-61
references: 16964229
references: Bioinformatics. 2006 Sep 1;22(17):2107-13
references: 16820429
references: Bioinformatics. 2012 Jul 1;28(13):1721-8
references: 22563066
references: Nucleic Acids Res. 2010 Jan;38(1):e4
references: 19854940
references: BMC Bioinformatics. 2007 Mar 21;8:98
references: 17376221
references: BMC Bioinformatics. 2008 Mar 26;9:164
references: 18366762
references: BMC Bioinformatics. 2010 Apr 29;11:221
references: 20429936
references: Science. 2008 Jun 6;320(5881):1344-9
references: 18451266
references: Nat Biotechnol. 2006 Sep;24(9):1115-22
references: 16964225
references: BMC Bioinformatics. 2009 Jul 09;10:211
references: 19589155
references: Proc Natl Acad Sci U S A. 2001 Jan 2;98(1):31-6
references: 11134512
references: RNA. 2008 Aug;14(8):1470-9
references: 18566192
references: Bioinformatics. 2005 Sep 15;21(18):3637-44
references: 16020470