Identifying Differential Exon Splicing Using Linear Models and Correlation Coefficients

Overview

Journal BMC Bioinformatics

Publisher Biomed Central

Specialty Biology

Date 2009 Jan 22

PMID 19154578

Citations 15

Authors

Sonia H Shah

Jacqueline A Pallas

Affiliations

Soon will be listed here.

Abstract

Background: With the availability of the Affymetrix exon arrays a number of tools have been developed to enable the analysis. These however can be expensive or have several pre-installation requirements. This led us to develop an analysis workflow for analysing differential splicing using freely available software packages that are already being widely used for gene expression analysis. The workflow uses the packages in the standard installation of R and Bioconductor (BiocLite) to identify differential splicing. We use the splice index method with the LIMMA framework. The main drawback with this approach is that it relies on accurate estimates of gene expression from the probe-level data. Methods such as RMA and PLIER may misestimate when a large proportion of exons are spliced. We therefore present the novel concept of a gene correlation coefficient calculated using only the probeset expression pattern within a gene. We show that genes with lower correlation coefficients are likely to be differentially spliced.

Results: The LIMMA approach was used to identify several tissue-specific transcripts and splicing events that are supported by previous experimental studies. Filtering the data is necessary, particularly removing exons and genes that are not expressed in all samples and cross-hybridising probesets, in order to reduce the false positive rate. The LIMMA approach ranked genes containing single or few differentially spliced exons much higher than genes containing several differentially spliced exons. On the other hand we found the gene correlation coefficient approach better for identifying genes with a large number of differentially spliced exons.

Conclusion: We show that LIMMA can be used to identify differential exon splicing from Affymetrix exon array data. Though further work would be necessary to develop the use of correlation coefficients into a complete analysis approach, the preliminary results demonstrate their usefulness for identifying differentially spliced genes. The two approaches work complementary as they can potentially identify different subsets of genes (single/few spliced exons vs. large transcript structure differences).

Citing Articles

Comprehensive insights on pivotal prognostic signature involved in clear cell renal cell carcinoma microenvironment using the ESTIMATE algorithm.

Luo J, Xie Y, Zheng Y, Wang C, Qi F, Hu J Cancer Med. 2020; 9(12):4310-4323.

PMID: 32311223 PMC: 7300420. DOI: 10.1002/cam4.2983.

Unusual structure and splicing pattern of the vertebrate mitochondrial solute carrier SLC25A3 gene.

Calvello R, Cianciulli A, Panaro M J Genet. 2018; 97(1):225-233.

PMID: 29666342

A random effects model for the identification of differential splicing (REIDS) using exon and HTA arrays.

Van Moerbeke M, Kasim A, Talloen W, Reumers J, Gohlmann H, Shkedy Z BMC Bioinformatics. 2017; 18(1):273.

PMID: 28545391 PMC: 5445373. DOI: 10.1186/s12859-017-1687-8.

Algorithms for differential splicing detection using exon arrays: a comparative assessment.

Zimmermann K, Jentsch M, Rasche A, Hummel M, Leser U BMC Genomics. 2016; 16:136.

PMID: 27391904 PMC: 4391533. DOI: 10.1186/s12864-015-1322-x.

iGEMS: an integrated model for identification of alternative exon usage events.

Sood S, Szkop K, Nakhuda A, Gallagher I, Murie C, Brogan R Nucleic Acids Res. 2016; 44(11):e109.

PMID: 27095197 PMC: 4914109. DOI: 10.1093/nar/gkw263.

References

OReilly M, Marshall E, Speirs H, Brown R . WNK1, a gene within a novel blood pressure control pathway, tissue-specifically generates radically different isoforms with and without a kinase domain. J Am Soc Nephrol. 2003; 14(10):2447-56. DOI: 10.1097/01.asn.0000089830.97681.3b. View

Verissimo F, Jordan P . WNK kinases, a novel protein kinase subfamily in multi-cellular organisms. Oncogene. 2001; 20(39):5562-9. DOI: 10.1038/sj.onc.1204726. View

Nakamura N, Miyake Y, Matsushita M, Tanaka S, Inoue H, Kanazawa H . KIF1Bbeta2, capable of interacting with CHP, is localized to synaptic vesicles. J Biochem. 2002; 132(3):483-91. DOI: 10.1093/oxfordjournals.jbchem.a003246. View

Okoniewski M, Yates T, Dibben S, Miller C . An annotation infrastructure for the analysis and interpretation of Affymetrix exon array data. Genome Biol. 2007; 8(5):R79. PMC: 1929135. DOI: 10.1186/gb-2007-8-5-r79. View

Irizarry R, Bolstad B, Collin F, Cope L, Hobbs B, Speed T . Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003; 31(4):e15. PMC: 150247. DOI: 10.1093/nar/gng015. View

John G, Shang Y, Li L, Renken C, Mannella C, Selker J . The mitochondrial inner membrane protein mitofilin controls cristae morphology. Mol Biol Cell. 2005; 16(3):1543-54. PMC: 551514. DOI: 10.1091/mbc.e04-08-0697. View

Das D, Clark T, Schweitzer A, Yamamoto M, Marr H, Arribere J . A correlation with exon expression approach to identify cis-regulatory elements for tissue-specific alternative splicing. Nucleic Acids Res. 2007; 35(14):4845-57. PMC: 1950531. DOI: 10.1093/nar/gkm485. View

Delaloy C, Lu J, Houot A, Disse-Nicodeme S, Gasc J, Corvol P . Multiple promoters in the WNK1 gene: one controls expression of a kidney-specific kinase-defective isoform. Mol Cell Biol. 2003; 23(24):9208-21. PMC: 309643. DOI: 10.1128/MCB.23.24.9208-9221.2003. View

Srinivasan K, Shiue L, Hayes J, Centers R, Fitzwater S, Loewen R . Detection and measurement of alternative splicing using splicing-sensitive microarrays. Methods. 2005; 37(4):345-59. DOI: 10.1016/j.ymeth.2005.09.007. View

10.

Huizing M, DePinto V, Ruitenbeek W, Trijbels F, van den Heuvel L, Wendel U . Importance of mitochondrial transmembrane processes in human mitochondriopathies. J Bioenerg Biomembr. 1996; 28(2):109-14. DOI: 10.1007/BF02110640. View

11.

Gong T, Winnicki R, Kohrman D, Lomax M . A novel mouse kinesin of the UNC-104/KIF1 subfamily encoded by the Kif1b gene. Gene. 1999; 239(1):117-27. DOI: 10.1016/s0378-1119(99)00370-4. View

12.

Fiermonte G, Dolce V, Palmieri F . Expression in Escherichia coli, functional characterization, and tissue distribution of isoforms A and B of the phosphate carrier from bovine mitochondria. J Biol Chem. 1998; 273(35):22782-7. DOI: 10.1074/jbc.273.35.22782. View

13.

Smyth G . Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2006; 3:Article3. DOI: 10.2202/1544-6115.1027. View

14.

Clark T, Schweitzer A, Chen T, Staples M, Lu G, Wang H . Discovery of tissue-specific exons using comprehensive human exon microarrays. Genome Biol. 2007; 8(4):R64. PMC: 1896007. DOI: 10.1186/gb-2007-8-4-r64. View

15.

Xing Y, Kapur K, Wong W . Probe selection and expression index computation of Affymetrix Exon Arrays. PLoS One. 2006; 1:e88. PMC: 1762343. DOI: 10.1371/journal.pone.0000088. View

16.

Purdom E, Simpson K, Robinson M, Conboy J, Lapuk A, Speed T . FIRMA: a method for detection of alternative splicing from exon array data. Bioinformatics. 2008; 24(15):1707-14. PMC: 2638867. DOI: 10.1093/bioinformatics/btn284. View

17.

Guipponi M, Scott H, Chen H, Schebesta A, Rossier C, Antonarakis S . Two isoforms of a human intersectin (ITSN) protein are produced by brain-specific alternative splicing in a stop codon. Genomics. 1998; 53(3):369-76. DOI: 10.1006/geno.1998.5521. View

18.

Gieffers C, Korioth F, Heimann P, Ungermann C, Frey J . Mitofilin is a transmembrane protein of the inner mitochondrial membrane expressed as two isoforms. Exp Cell Res. 1997; 232(2):395-9. DOI: 10.1006/excr.1997.3539. View

19.

Robinson M, Speed T . A comparison of Affymetrix gene expression arrays. BMC Bioinformatics. 2007; 8:449. PMC: 2216046. DOI: 10.1186/1471-2105-8-449. View

20.

Mayr J, Merkel O, Kohlwein S, Gebhardt B, Bohles H, Fotschl U . Mitochondrial phosphate-carrier deficiency: a novel disorder of oxidative phosphorylation. Am J Hum Genet. 2007; 80(3):478-84. PMC: 1821108. DOI: 10.1086/511788. View