Gene structure-based splice variant deconvolution using a microarry platform

MOTIVATION Alternative splicing allows a single gene to generate multiple mRNAs, which can be translated into functionally and structurally diverse proteins. One gene can have multiple variants coexisting at different concentrations. Estimating the relative abundance of each variant is important for the study of underlying biological function. Microarrays are standard tools that measure gene expression. But most design and analysis has not accounted for splice variants. Thus splice variant-specific chip designs and analysis algorithms are needed for accurate gene expression profiling. RESULTS Inspired by Li and Wong (2001), we developed a gene structure-based algorithm to determine the relative abundance of known splice variants. Probe intensities are modeled across multiple experiments using gene structures as constraints. Model parameters are obtained through a maximum likelihood estimation (MLE) process/framework. The algorithm produces the relative concentration of each variant, as well as an affinity term associated with each probe. Validation of the algorithm is performed by a set of controlled spike experiments as well as endogenous tissue samples using a human splice variant array.

[1]  Teresa A. Webster,et al.  Probe selection for high-density oligonucleotide arrays , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Masaru Tomita,et al.  Inferring alternative splicing patterns in mouse from a full-length cDNA library and microarray data. , 2002, Genome research.

[3]  S. P. Fodor,et al.  Large-Scale Transcriptional Activity in Chromosomes 21 and 22 , 2002, Science.

[4]  Tyson A. Clark,et al.  Genomewide Analysis of mRNA Processing in Yeast Using Splicing-Specific Microarrays , 2002, Science.

[5]  Xiang-Dong Fu,et al.  Profiling alternative splicing on fiber-optic arrays , 2002, Nature Biotechnology.

[6]  D. Haussler,et al.  Assembly of the working draft of the human genome with GigAssembler. , 2001, Genome research.

[7]  T. Jatkoe,et al.  Predicting splice variant from DNA chip expression data. , 2001, Genome research.

[8]  M B Eisen,et al.  Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Timothy B. Stockwell,et al.  The Sequence of the Human Genome , 2001, Science.

[10]  J. V. Moran,et al.  Initial sequencing and analysis of the human genome. , 2001, Nature.

[11]  R. Stoughton,et al.  Experimental annotation of the human genome using microarray technology , 2001, Nature.

[12]  C. Li,et al.  Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[13]  International Human Genome Sequencing Consortium.  Initial sequencing and analysis of the human genome , 2001, Nature.

[14]  Kevin Burrage,et al.  ISIS, the intron information system, reveals the high frequency of alternative splicing in the human genome , 2000, Nature Genetics.

[15]  Arya M. Sharma,et al.  Association of a human G-protein β3 subunit variant with hypertension , 1998, Nature Genetics.

[16]  U. Francke,et al.  Silent mutation induces exon skipping of fibrillin-1 gene in Marfan syndrome , 1997, Nature Genetics.

[17]  K. Klinger,et al.  Alternative splicing of exon 3 of the human growth hormone receptor is the result of an unusual genetic polymorphism. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Y. D'Aubenton-Carafa,et al.  Tissue-specific splicing of two mutually exclusive exons of the chicken beta-tropomyosin pre-mRNA: positive and negative regulations. , 1996, Biochimie.

[19]  D. Helfman,et al.  Alternatively spliced exons of the beta tropomyosin gene exhibit different affinities for F-actin and effects with nonmuscle caldesmon. , 1995, Journal of cell science.

[20]  K. Beisel,et al.  Identification of novel alternatively spliced isoforms of the tropomyosin-encoding gene, TMnm, in the rat cochlea. , 1994, Gene.

[21]  S. Bernstein,et al.  Genetic and biochemical analysis of alternative RNA splicing. , 1994, Advances in genetics.

[22]  C. Lin,et al.  Human fibroblast tropomyosin isoforms: characterization of cDNA clones and analysis of tropomyosin isoform expression in human tissues and in normal and transformed cells. , 1993, Cell motility and the cytoskeleton.

[23]  D. Helfman,et al.  The molecular basis for tropomyosin isoform diversity , 1991, BioEssays : news and reviews in molecular, cellular and developmental biology.

[24]  E. Brody,et al.  RNA secondary structure repression of a muscle-specific exon in HeLa cell nuclear extracts. , 1991, Science.

[25]  J. G. Patton,et al.  Alternative splicing in the control of gene expression. , 1989, Annual review of genetics.

[26]  J. Sutcliffe,et al.  Alternative mRNA splicing: the Shaker gene. , 1988, Trends in genetics : TIG.

[27]  J. Liautard,et al.  Complete nucleotide sequence of the adult skeletal isoform of human skeletal muscle β-tropomyosin , 1988 .

[28]  D. Mccormick.  Sequence the Human Genome , 1986, Bio/Technology.

[29]  D. Helfman,et al.  Nonmuscle and muscle tropomyosin isoforms are expressed from a single gene by alternative RNA splicing and polyadenylation , 1986, Molecular and cellular biology.

[30]  F. Walsh,et al.  A muscle-type tropomyosin in human fibroblasts: evidence for expression by an alternative RNA splicing mechanism. , 1985, Proceedings of the National Academy of Sciences of the United States of America.