The gene guessing game

Methods in comparative genomics: genome correspondence, gene identification and regulatory motif discovery

In Kellis et al. (2003), we reported the genome sequences of S. paradoxus, S. mikatae, and S. bayanus and compared these three yeast species to their close relative, S. cerevisiae. Genomewide comparative analysis allowed the identification of functionally important sequences, both coding and noncoding. In this companion paper we describe the mathematical and algorithmic results underpinning the analysis of these genomes. (1) We present methods for the automatic determination of genome correspondence. The algorithms enabled the automatic identification of orthologs for more than 90% of genes and intergenic regions across the four species despite the large number of duplicated genes in the yeast genome. The remaining ambiguities in the gene correspondence revealed recent gene family expansions in regions of rapid genomic change. (2) We present methods for the identification of protein-coding genes based on their patterns of nucleotide conservation across related species. We observed the pressure to conserve the reading frame of functional proteins and developed a test for gene identification with high sensitivity and specificity. We used this test to revisit the genome of S. cerevisiae, reducing the overall gene count by 500 genes (10% of previously annotated genes) and refining the gene structure of hundreds of genes. (3) We present novel methods for the systematic de novo identification of regulatory motifs. The methods do not rely on previous knowledge of gene function and in that way differ from the current literature on computational motif discovery. Based on genomewide conservation patterns of known motifs, we developed three conservation criteria that we used to discover novel motifs. We used an enumeration approach to select strongly conserved motif cores, which we extended and collapsed into a small number of candidate regulatory motifs. These include most previously known regulatory motifs as well as several noteworthy novel motifs. The majority of discovered motifs are enriched in functionally related genes, allowing us to infer a candidate function for novel motifs. Our results demonstrate the power of comparative genomics to further our understanding of any species. Our methods are validated by the extensive experimental knowledge in yeast and will be invaluable in the study of complex genomes like that of the human.

Genes, isochores and bands in human chromosomes 21 and 22

The recently available DNA sequences from chromosomes 21 and 22 enabled us to define the relationships of different band types with isochores and with gene concentration and to compare these relationships with previous results. We showed that chromosomal bands appear as Giemsa or Reverse bands depending not on their absolute GC level, but on the composition GC level relative to those of adjacent contiguous bands. We also demonstrated that the GC-richest, and gene-richest H3+ bands are characterized by a lower DNA compaction compared with the GC-poorest, gene-poorest L1+ bands. Moreover, our results indicate that the human genome contains about 30,000 genes.

Novel approaches for identifying genes regulating lymphocyte development and function

The draft sequence of the human and mouse genomes provides an unparalleled opportunity for understanding the genetic control of immune-cell development. Strategies can begin with a gene sequence and pursue a putative immune-system function by employing mRNA-expression profiling or creating gene knockouts in embryonic stem cells. The latter can be produced by utilising the Cre/Lox system, a tetracycline operon, a gene-trap method or chemical mutagenesis. Alternatively, mutant phenotypes (derived using the mutagen ethylnitrosourea) can be traced back to gene sequences.

Assessment of the total number of human transcription units

Variation in the estimates of the number of genes encoded by the human genome (28,000-120,000) attests to the difficulty of systematically identifying human genes. Sequencing of human chromosome 22 (Chr22) provided the first comprehensive, unbiased view of an entire human chromosome, and intensive analysis of this sequence identified 545 genes and 134 pseudogenes that had similarity or identity to known proteins and/or ESTs and which were listed in the gene annotation (http://www.sanger.ac.uk/HGP/Chr22). This analysis yielded an estimate of approximately 36,000 functional expressed genes in the human genome (and 9000 pseudogenes). However, a key uncertainty in this estimate was that hundreds of additional genes beyond those annotated in the Chr22 sequence are predicted by the gene prediction program Genscan, an unknown number of which might represent additional expressed genes. To determine what fraction of these "predicted novel genes" (PNGs) represents expressed human genes, we used a sensitive RT-PCR assay to detect predicted transcripts in 17 tissues and one cell line. Our results indicate that at least 5000-9000 additional human genes which lack similarity to known genes or proteins exist in the human genome, increasing baseline gene estimates to approximately 41,000-45,000.

The gene guessing game

A recent flurry of publications and media attention has revived interest in the question of how many genes exist in the human genome. Here, I review the estimates and use genomic sequence data from human chromosomes 21 and 22 to establish my own prediction.