A comparative molecular analysis of developing mouse forelimbs and hindlimbs using serial analysis of gene expression (SAGE).

The analysis of differentially expressed genes is a powerful approach to elucidate the genetic mechanisms underlying the morphological and evolutionary diversity among serially homologous structures, both within the same organism (e.g., hand vs. foot) and between different species (e.g., hand vs. wing). In the developing embryo, limb-specific expression of Pitx1, Tbx4, and Tbx5 regulates the determination of limb identity. However, numerous lines of evidence, including the fact that these three genes encode transcription factors, indicate that additional genes are involved in the Pitx1-Tbx hierarchy. To examine the molecular distinctions coded for by these factors, and to identify novel genes involved in the determination of limb identity, we have used Serial Analysis of Gene Expression (SAGE) to generate comprehensive gene expression profiles from intact, developing mouse forelimbs and hindlimbs. To minimize the extraction of erroneous SAGE tags from low-quality sequence data, we used a new algorithm to extract tags from -analyzed sequence data and obtained 68,406 and 68,450 SAGE tags from forelimb and hindlimb SAGE libraries, respectively. We also developed an improved method for determining the identity of SAGE tags that increases the specificity of and provides additional information about the confidence of the tag-UniGene cluster match. The most differentially expressed gene between our SAGE libraries was Pitx1. The differential expression of Tbx4, Tbx5, and several limb-specific Hox genes was also detected; however, their abundances in the SAGE libraries were low. Because numerous other tags were differentially expressed at this low level, we performed a 'virtual' subtraction with 362,344 tags from six additional nonlimb SAGE libraries to further refine this set of candidate genes. This subtraction reduced the number of candidate genes by 74%, yet preserved the previously identified regulators of limb identity. This study presents the gene expression complexity of the developing limb and identifies candidate genes involved in the regulation of limb identity. We propose that our computational tools and the overall strategy used here are broadly applicable to other SAGE-based studies in a variety of organisms. [SAGE data are all available at GEO (http://www.ncbi.nlm.nih.gov/geo/) under accession nos. GSM55 and GSM56, which correspond to the forelimb and hindlimb raw SAGE data.]

[1]  L. Wolpert Developmental Biology , 1968, Nature.

[2]  F. James Rohlf,et al.  Biometry: The Principles and Practice of Statistics in Biological Research , 1969 .

[3]  J. Bishop,et al.  The expression of three abundance classes of messenger RNA in mouse tissues , 1976, Cell.

[4]  L. Pikó,et al.  Patterns of mRNA prevalence and expression of B1 and B2 transcripts in early mouse embryos. , 1987, Development.

[5]  M. Frohman,et al.  Rapid production of full-length cDNAs from rare transcripts: amplification using a single gene-specific oligonucleotide primer. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[6]  S. Bryant,et al.  A staging system for mouse limb development. , 1989, The Journal of experimental zoology.

[7]  D. Engelke,et al.  Purification of Thermus aquaticus DNA polymerase expressed in Escherichia coli. , 1990, Analytical biochemistry.

[8]  D. Jacobs,et al.  Hox-3.6: isolation and characterization of a new murine homeobox gene located in the 5′ region of the Hox-3 cluster , 1992, Mechanisms of Development.

[9]  C. Caskey,et al.  Human and murine dystrophin mRNA transcripts are differentially expressed during skeletal muscle, heart, and brain development. , 1992, Nucleic acids research.

[10]  J. Chamberlain,et al.  A B2 repeat insertion generates alternate structures of the mouse muscle gamma-phosphorylase kinase gene. , 1993, Genomics.

[11]  W. Thilly,et al.  Specificity, efficiency, and fidelity of PCR. , 1993, PCR methods and applications.

[12]  T. Papenbrock,et al.  The murine Hoxc cluster contains five neighboring AbdB-related Hox genes that show unique spatially coordinated expression in posterior embryonic subregions , 1994, Mechanisms of Development.

[13]  D. Witte,et al.  Hoxa 11 structure, extensive antisense transcription, and function in male and female fertility. , 1995, Development.

[14]  C. Tabin,et al.  Analysis of Hox gene expression in the chick limb bud. , 1996, Development.

[15]  V. Papaioannou,et al.  Evidence of a role for T-☐ genes in the evolution of limb morphogenesis and the specification of forelimb/hindlimb identity , 1996, Mechanisms of Development.

[16]  M. Rosenfeld,et al.  P-OTX: a PIT-1-interacting homeodomain factor expressed during anterior pituitary gland development. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[17]  M. Cohn,et al.  Limbs: a model for pattern formation within the vertebrate body plan. , 1996, Trends in genetics : TIG.

[18]  David I. Wilson,et al.  Holt-Oram syndrome is caused by mutations in TBX5, a member of the Brachyury (T) gene family , 1997, Nature Genetics.

[19]  R H Hruban,et al.  Gene expression profiles in normal and cancer cells. , 1997, Science.

[20]  C. Tabin,et al.  Molecular Models for Vertebrate Limb Development , 1997, Cell.

[21]  J. Claverie,et al.  The significance of digital gene expression profiles. , 1997, Genome research.

[22]  R. Kucherlapati,et al.  Mutations in human cause limb and cardiac malformation in Holt-Oram syndrome , 1997, Nature Genetics.

[23]  Wei Zhou,et al.  Characterization of the Yeast Transcriptome , 1997, Cell.

[24]  M. Cohn,et al.  Hox9 genes and vertebrate limb specification , 1997, nature.

[25]  J. Seidman,et al.  Erratum: Mutations in human TBX5 cause limb and cardiac malformation in Holt-Oram syndrome , 1997, Nature Genetics.

[26]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[27]  K. Patel,et al.  Tbx genes and limb identity in chick embryo development. , 1998, Development.

[28]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[29]  L. Silver,et al.  Involvement of T-box genes Tbx2-Tbx5 in vertebrate limb specification and development. , 1998, Development.

[30]  M. Capecchi,et al.  The mouse Hoxc11 gene: genomic structure and expression pattern , 1998, Mechanisms of Development.

[31]  Concepción Rodríguez-Esteban,et al.  The T-box genes Tbx4 and Tbx5 regulate limb outgrowth and identity , 1999, Nature.

[32]  K. Tamura,et al.  Differential expression of Tbx4 and Tbx5 in Zebrafish Fin buds , 1999, Mechanisms of Development.

[33]  K. Yasuda,et al.  Tbx5 and Tbx4 genes determine the wing/leg identity of limb buds , 1999, Nature.

[34]  C. Lanctôt,et al.  Hindlimb patterning and mandible development require the Ptx1 gene. , 1999, Development.

[35]  L. Niswander Developmental biology: Legs to wings and back again , 1999, Nature.

[36]  S. Carroll,et al.  Selector Genes and Limb Identity in Arthropods and Vertebrates , 1999, Cell.

[37]  C. Tabin,et al.  Role of Pitx1 upstream of Tbx4 in specification of hindlimb identity. , 1999, Science.

[38]  P M Bossuyt,et al.  Genes differentially expressed in medulloblastoma and fetal brain. , 1999, Physiological genomics.

[39]  G. Landes,et al.  Analysis of human transcriptomes , 1999, Nature Genetics.

[40]  A. Kassam,et al.  Comprehensive transcript analysis in small quantities of mRNA by SAGE-lite. , 1999, Nucleic acids research.

[41]  N. Datson,et al.  MicroSAGE: a modified procedure for serial analysis of gene expression in limited amounts of tissue. , 1999, Nucleic acids research.

[42]  Concepción Rodríguez-Esteban,et al.  Role of the Bicoid-related homeodomain factor Pitx1 in specifying hindlimb morphogenesis and pituitary development. , 1999, Genes & development.

[43]  J. Buhler,et al.  Serial microanalysis of renal transcriptomes. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[44]  S. Altschul,et al.  SAGEmap: a public gene expression resource. , 2000, Genome research.

[45]  K. Takeshima,et al.  Conserved and divergent expression of T-box genes Tbx2-Tbx5 in Xenopus , 2000, Mechanisms of Development.

[46]  Yixin Wang,et al.  POWER_SAGE: comparing statistical tests for SAGE experiments , 2000, Bioinform..

[47]  E. H. Margulies,et al.  Building Arms or Legs with Molecular Models , 2000, Pediatric Research.

[48]  J. I. Izpisúa Belmonte,et al.  Perspectives on the evolutionary origin of tetrapod limbs. , 2000, The Journal of experimental zoology.

[49]  T. Papenbrock,et al.  Loss of fibula in mice overexpressing Hoxc11 , 2000, Mechanisms of Development.

[50]  P. Green,et al.  Analysis of expressed sequence tags indicates 35,000 human genes , 2000, Nature Genetics.

[51]  I. Ruvinsky,et al.  Genetic and developmental bases of serial homology in vertebrate limb evolution. , 2000, Development.

[52]  D. Gautheret,et al.  Patterns of variant polyadenylation signal usage in human genes. , 2000, Genome research.

[53]  Elliott H. Margulies,et al.  eSAGE: managing and analysing data generated with Serial Analysis of Gene Expression (SAGE) , 2000, Bioinform..

[54]  C. Fizames,et al.  Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence , 2000, Nature Genetics.

[55]  John Quackenbush,et al.  Gene Index analysis of the human genome estimates approximately 120,000 genes , 2000, Nature Genetics.

[56]  X. Estivill,et al.  The mouse brain transcriptome by SAGE: differences in gene expression between P30 brains of the partial trisomy 16 mouse model of Down syndrome (Ts65Dn) and normals. , 2000, Genome research.

[57]  F. Baas,et al.  The Human Transcriptome Map: Clustering of Highly Expressed Genes in Chromosomal Domains , 2001, Science.

[58]  E. H. Margulies,et al.  Identification and prevention of a GC content bias in SAGE libraries. , 2001, Nucleic acids research.

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

[60]  T. Liesegang The human transcriptome map: Clustering of highly expressed genes in chromosomal domains. Caron H, ∗ van Schaik B, van der Mee M, et al. Science 2001;291:1289–1292. , 2001 .

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

[62]  Ji Huang,et al.  [Serial analysis of gene expression]. , 2002, Yi chuan = Hereditas.