Vertebrate genome sequencing: building a backbone for comparative genomics

Newly evolved introns in human retrogenes provide novel insights into their evolutionary roles

BACKGROUND

Retrogenes generally do not contain introns. However, in some instances, retrogenes may recruit internal exonic sequences as introns, which is known as intronization. A retrogene that undergoes intronization is a good model with which to investigate the origin of introns. Nevertheless, previously, only two cases in vertebrates have been reported.

RESULTS

In this study, we systematically screened the human (Homo sapiens) genome for retrogenes that evolved introns and analyzed their patterns in structure, expression and origin. In total, we identified nine intron-containing retrogenes. Alignment of pairs of retrogenes and their parents indicated that, in addition to intronization (five cases), retrogenes also may have gained introns by insertion of external sequences into the genes (one case) or reversal of the orientation of transcription (three cases). Interestingly, many intronizations were promoted not by base substitutions but by cryptic splice sites, which were silent in the parental genes but active in the retrogenes. We also observed that the majority of introns generated by intronization did not involve frameshifts.

CONCLUSIONS

Intron gains in retrogenes are not as rare as previously thought. Furthermore, diverse mechanisms may lead to intron creation in retrogenes. The activation of cryptic splice sites in the intronization of retrogenes may be triggered by the change of gene structure after retroposition. A high percentage of non-frameshift introns in retrogenes may be because non-frameshift introns do not dramatically affect host proteins. Introns generated by intronization in human retrogenes are generally young, which is consistent with previous findings for Caenorhabditis elegans. Our results provide novel insights into the evolutionary role of introns.

Changes of globin expression in the Japanese medaka (Oryzias latipes) in response to acute and chronic hypoxia

Fishes live in an aquatic environment with low or temporally changing O(2) availability. Variations in O(2) levels require many anatomical, behavioral, physiological, and biochemical adaptations that ensure the uptake of an adequate amount of O(2). Some fish species are comparatively well adapted to tolerate low O(2) partial pressure (hypoxia). The Japanese ricefish medaka (Oryzias latipes) is an important model organism for biomedical research that shows remarkable tolerance towards hypoxia. We have investigated the regulation and role of globins under hypoxia. We applied four different regimes of chronic hypoxia (24 and 48 h at PO(2) = 2 or 4 kPa) as well as acute hypoxia (2 h at PO(2) = 0.5 kPa) to adult male medaka. Changes of mRNA levels of seven globin genes (adult hemoglobin α and β, myoglobin, neuroglobin, cytoglobin 1 and 2, globin X), three hypoxia-response genes (lactate dehydrogenase b, phosphoglycerate kinase, adrenomedullin 1) and two putative reference genes (cyclophilin, acidic ribosomal phosphoprotein P0) were monitored by means of quantitative real-time reverse-transcription PCR. We observed strong upregulation of myoglobin, which is also expressed in the medaka brain, as previously demonstrated for carp, goldfish and zebrafish. The hemoglobin chains were found upregulated, whereas earlier studies found down-regulation of hemoglobin in hypoxic zebrafish. By contrast, neuroglobin mRNA was not affected by hypoxia in medaka, but had been found upregulated in zebrafish. Globin X is induced in medaka brain, but down-regulated in zebrafish. Thus, the patterns of hypoxia response of globins are strikingly different in various fish species, which can be interpreted as indication for different roles of the various globins in hypoxia response and for alternative metabolic strategies of fish species in coping with O(2) deprivation.

Chromosomal redistribution of male-biased genes in mammalian evolution with two bursts of gene gain on the X chromosome

Mammalian X chromosomes evolved under various mechanisms including sexual antagonism, the faster-X process, and meiotic sex chromosome inactivation (MSCI). These forces may contribute to nonrandom chromosomal distribution of sex-biased genes. In order to understand the evolution of gene content on the X chromosome and autosome under these forces, we dated human and mouse protein-coding genes and miRNA genes on the vertebrate phylogenetic tree. We found that the X chromosome recently acquired a burst of young male-biased genes, which is consistent with fixation of recessive male-beneficial alleles by sexual antagonism. For genes originating earlier, however, this pattern diminishes and finally reverses with an overrepresentation of the oldest male-biased genes on autosomes. MSCI contributes to this dynamic since it silences X-linked old genes but not X-linked young genes. This demasculinization process seems to be associated with feminization of the X chromosome with more X-linked old genes expressed in ovaries. Moreover, we detected another burst of gene originations after the split of eutherian mammals and opossum, and these genes were quickly incorporated into transcriptional networks of multiple tissues. Preexisting X-linked genes also show significantly higher protein-level evolution during this period compared to autosomal genes, suggesting positive selection accompanied the early evolution of mammalian X chromosomes. These two findings cast new light on the evolutionary history of the mammalian X chromosome in terms of gene gain, sequence, and expressional evolution.

Evolutionary conservation of zinc finger transcription factor binding sites in promoters of genes co-expressed with WT1 in prostate cancer

BACKGROUND

Gene expression analyses have led to a better understanding of growth control of prostate cancer cells. We and others have identified the presence of several zinc finger transcription factors in the neoplastic prostate, suggesting a potential role for these genes in the regulation of the prostate cancer transcriptome. One of the transcription factors (TFs) identified in the prostate cancer epithelial cells was the Wilms tumor gene (WT1). To rapidly identify coordinately expressed prostate cancer growth control genes that may be regulated by WT1, we used an in silico approach.

RESULTS

Evolutionary conserved transcription factor binding sites (TFBS) recognized by WT1, EGR1, SP1, SP2, AP2 and GATA1 were identified in the promoters of 24 differentially expressed prostate cancer genes from eight mammalian species. To test the relationship between sequence conservation and function, chromatin of LNCaP prostate cancer and kidney 293 cells were tested for TF binding using chromatin immunoprecipitation (ChIP). Multiple putative TFBS in gene promoters of placental mammals were found to be shared with those in human gene promoters and some were conserved between genomes that diverged about 170 million years ago (i.e., primates and marsupials), therefore implicating these sites as candidate binding sites. Among those genes coordinately expressed with WT1 was the kallikrein-related peptidase 3 (KLK3) gene commonly known as the prostate specific antigen (PSA) gene. This analysis located several potential WT1 TFBS in the PSA gene promoter and led to the rapid identification of a novel putative binding site confirmed in vivo by ChIP. Conversely for two prostate growth control genes, androgen receptor (AR) and vascular endothelial growth factor (VEGF), known to be transcriptionally regulated by WT1, regulatory sequence conservation was observed and TF binding in vivo was confirmed by ChIP.

CONCLUSION

Overall, this targeted approach rapidly identified important candidate WT1-binding elements in genes coordinately expressed with WT1 in prostate cancer cells, thus enabling a more focused functional analysis of the most likely target genes in prostate cancer progression. Identifying these genes will help to better understand how gene regulation is altered in these tumor cells.

Construction, alignment and analysis of twelve framework physical maps that represent the ten genome types of the genus Oryza

Bacterial artificial chromosome (BAC) fingerprint and end-sequenced physical maps representing the ten genome types of Oryza are presented

We describe the establishment and analysis of a genus-wide comparative framework composed of 12 bacterial artificial chromosome fingerprint and end-sequenced physical maps representing the 10 genome types of Oryza aligned to the O. sativa ssp. japonica reference genome sequence. Over 932 Mb of end sequence was analyzed for repeats, simple sequence repeats, miRNA and single nucleotide variations, providing the most extensive analysis of Oryza sequence to date.

An evaluation of information content as a metric for the inference of putative conserved noncoding regions in DNA sequences using a genetic algorithms approach

In previous work, we presented GAMI, an approach to motif inference that uses a genetic algorithms search. GAMI is designed specifically to find putative conserved regulatory motifs in noncoding regions of divergent species, and is designed to allow for analysis of long nucleotide sequences. In this work, we compare GAMI's performance when run with its original fitness function (a simple count of the number of matches) and when run with information content, as well as several variations on these metrics. Results indicate that information content does not identify highly conserved regions, and thus is not the appropriate metric for this task, while variations on information content as well as the original metric succeed in identifying putative conserved regions.

Identification of novel light-induced genes in the suprachiasmatic nucleus

Background

The transmission of information about the photic environment to the circadian clock involves a complex array of neurotransmitters, receptors, and second messenger systems. Exposure of an animal to light during the subjective night initiates rapid transcription of a number of immediate-early genes in the suprachiasmatic nucleus of the hypothalamus. Some of these genes have known roles in entraining the circadian clock, while others have unknown functions. Using laser capture microscopy, microarray analysis, and quantitative real-time PCR, we performed a comprehensive screen for changes in gene expression immediately following a 30 minute light pulse in suprachiasmatic nucleus of mice.

Results

The results of the microarray screen successfully identified previously known light-induced genes as well as several novel genes that may be important in the circadian clock. Newly identified light-induced genes include early growth response 2, proviral integration site 3, growth-arrest and DNA-damage-inducible 45 beta, and TCDD-inducible poly(ADP-ribose) polymerase. Comparative analysis of promoter sequences revealed the presence of evolutionarily conserved CRE and associated TATA box elements in most of the light-induced genes, while other core clock genes generally lack this combination of promoter elements.

Conclusion

The photic signalling cascade in the suprachiasmatic nucleus activates an array of immediate-early genes, most of which have unknown functions in the circadian clock. Detected evolutionary conservation of CRE and TATA box elements in promoters of light-induced genes suggest that the functional role of these elements has likely remained the same over evolutionary time across mammalian orders.

Comparative genomics in vertebrate evolution and development

The vast quantities of publicly available DNA sequencing data and genome resources are enabling biologists to investigate age-old problems in biology that were not addressable previously. In this review, we discuss how comparative genomics is practiced and how the data can be used to make biological inferences with respect to vertebrate evolution and development. Examples are taken from the well-known HOX clusters, which are always a high-priority target for genomic analyses due to their inferred role in the evolution of metazoans. In addition, we briefly discuss the application of genomic approaches to problems in comparative endocrinology.

Recapitulation of human betaB1-crystallin promoter activity in transgenic zebrafish

Development of the eye is morphologically similar among vertebrates, indicating that the underlying mechanism regulating the process may have been highly conserved during evolution. Herein we analyzed the promoter of the human betaB1-crytallin gene in zebrafish by transgenic experiments. To delineate the evolutionarily conserved regulatory elements, we performed serial deletion assays in the promoter region. The results demonstrated that the -90/+61-bp upstream proximal promoter region is sufficient to confer lens-tissue specificity to the human betaB1-crystallin gene in transgenic zebrafish. Through phylogenetic sequence comparisons and an electrophoretic mobility shift assay (EMSA), a highly conserved cis-element of a six-base pair sequence TG(A/C)TGA, the consensus sequence for the Maf protein binding site, within the proximal promoter region was revealed. Further, a site-mutational assay showed that this element is crucial for promoter activity. These data suggest that the fundamental transcriptional regulatory mechanism of the betaB1-crystallin gene has been well conserved between humans and zebrafish, and plausibly among all vertebrates, during evolution.

Marine organism cell biology and regulatory sequence discoveryin comparative functional genomics

The use of bioinformatics to integrate phenotypic and genomic data from mammalian models is well established as a means of understanding human biology and disease. Beyond direct biomedical applications of these approaches in predicting structure–function relationships between coding sequences and protein activities, comparative studies also promote understanding of molecular evolution and the relationship between genomic sequence and morphological and physiological specialization. Recently recognized is the potential of comparative studies to identify functionally significant regulatory regions and to generate experimentally testable hypotheses that contribute to understanding mechanisms that regulate gene expression, including transcriptional activity, alternative splicing and transcript stability. Functional tests of hypotheses generated by computational approaches require experimentally tractable in vitro systems, including cell cultures. Comparative sequence analysis strategies that use genomic sequences from a variety of evolutionarily diverse organisms are critical for identifying conserved regulatory motifs in the 5′-upstream, 3′-downstream and introns of genes. Genomic sequences and gene orthologues in the first aquatic vertebrate and protovertebrate organisms to be fully sequenced (Fugu rubripes, Ciona intestinalis, Tetraodon nigroviridis, Danio rerio) as well as in the elasmobranchs, spiny dogfish shark (Squalus acanthias) and little skate (Raja erinacea), and marine invertebrate models such as the sea urchin (Strongylocentrotus purpuratus) are valuable in the prediction of putative genomic regulatory regions. Cell cultures have been derived for these and other model species. Data and tools resulting from these kinds of studies will contribute to understanding transcriptional regulation of biomedically important genes and provide new avenues for medical therapeutics and disease prevention.

Targeted, haplotype-resolved resequencing of long segments of the human genome

Currently, challenges exist to acquire long-range (hundreds of kilobase pairs) phase-discriminated sequence across substantial numbers of individuals. We have developed a straightforward method for isolating and characterizing specific genomic regions in a haplospecific manner. Real-time PCR is carried out to STS content map and genotype pools of fosmid clones arrayed in 384-well microtiter plates. Single-nucleotide polymorphisms, microsatellite markers, and insertion-deletion polymorphisms are used to differentiate the target region into haplotype-specific tiling paths. DNA of clones from these tiling paths is retrieved from the library and either sequenced by standard shotgun methods or amplified in vitro and sequenced by a primer-based, directed method. This approach provides convenient access to complete, haplotype-resolved resequencing data from multiple individuals across tens to hundreds of thousands of basepairs. We illustrate its implementation with a detailed example of more than 400 kbp from the human CFTR region, across 15 individuals, and summarize our experience applying it to many other human loci.

Characterizing the chromosomes of the Australian model marsupial Macropus eugenii (tammar wallaby)

Marsupials occupy a phylogenetic middle ground that is very valuable in genome comparisons of mammal and other vertebrate species. For this reason, whole genome sequencing is being undertaken for two distantly related marsupial species, including the model kangaroo species Macropus eugenii (the tammar wallaby). As a first step towards the molecular characterization of the tammar genome, we present a detailed description of the tammar karyotype, report the development of a set of molecular anchor markers and summarize the comparative mapping data for this species.

A 6-Mb contig-based comparative gene and linkage map of the rat schwannoma tumor suppressor region at 10q32.3

Frequent genetic aberrations of malignant schwannomas induced by the alkylating agent N-ethyl-N-nitrosourea in hybrids from inbred BD rat strains include allelic imbalances of the telomeric 20 Mb of chromosome 5 (Dis-2) and of the telomeric 5 Mb of chromosome 10q32 (Dis-1) in 59 and 94% of the tumors, respectively. The Dis-1 minimal loss of heterozygosity consensus region extends from D10Rat4 to the telomere and harbors a putative tumor suppressor gene(s). We constructed a 6-Mb BAC/PAC contig containing more than 70 known genes, 18 mapped microsatellites, and further ESTs/reference RNAs. A continuous block of strongly conserved synteny with mouse chromosome 11E2 and human chromosome 17q25.3 was found. Combining the sequence information from the rat and closely related syntenic regions of different mammalian species produces nearly complete gene maps as a basis for a positional candidate approach and gives insight into mammalian genomic evolution.

Annotation of cis-regulatory elements by identification, subclassification, and functional assessment of multispecies conserved sequences

An important step toward improving the annotation of the human genome is to identify cis-acting regulatory elements from primary DNA sequence. One approach is to compare sequences from multiple, divergent species. This approach distinguishes multispecies conserved sequences (MCS) in noncoding regions from more rapidly evolving neutral DNA. Here, we have analyzed a region of approximately 238kb containing the human alpha globin cluster that was sequenced and/or annotated across the syntenic region in 22 species spanning 500 million years of evolution. Using a variety of bioinformatic approaches and correlating the results with many aspects of chromosome structure and function in this region, we were able to identify and evaluate the importance of 24 individual MCSs. This approach sensitively and accurately identified previously characterized regulatory elements but also discovered unidentified promoters, exons, splicing, and transcriptional regulatory elements. Together, these studies demonstrate an integrated approach by which to identify, subclassify, and predict the potential importance of MCSs.

Identification of neuronal enhancers of the proopiomelanocortin gene by transgenic mouse analysis and phylogenetic footprinting

The proopiomelanocortin (POMC) gene is expressed in the pituitary and arcuate neurons of the hypothalamus. POMC arcuate neurons play a central role in the control of energy homeostasis, and rare loss-of-function mutations in POMC cause obesity. Moreover, POMC is the prime candidate gene within a highly significant quantitative trait locus on chromosome 2 associated with obesity traits in several human populations. Here, we identify two phylogenetically conserved neuronal POMC enhancers designated nPE1 (600 bp) and nPE2 (150 bp) located approximately 10 to 12 kb upstream of mammalian POMC transcriptional units. We show that mouse or human genomic regions containing these enhancers are able to direct reporter gene expression to POMC hypothalamic neurons, but not the pituitary of transgenic mice. Conversely, deletion of nPE1 and nPE2 in the context of the entire transcriptional unit of POMC abolishes transgene expression in the hypothalamus without affecting pituitary expression. Our results indicate that the nPEs are necessary and sufficient for hypothalamic POMC expression and that POMC expression in the brain and pituitary is controlled by independent sets of enhancers. Our study advances the understanding of the molecular nature of hypothalamic POMC neurons and will be useful to determine whether polymorphisms in POMC regulatory regions play a role in the predisposition to obesity.

Comparative sequencing provides insights about the structure and conservation of marsupial and monotreme genomes

Sequencing and comparative analyses of genomes from multiple vertebrates are providing insights about the genetic basis for biological diversity. To date, these efforts largely have focused on eutherian mammals, chicken, and fish. In this article, we describe the generation and study of genomic sequences from noneutherian mammals, a group of species occupying unusual phylogenetic positions. A large sequence data set (totaling >5 Mb) was generated for the same orthologous region in three marsupial (North American opossum, South American opossum, and Australian tammar wallaby) and one monotreme (platypus) genomes. These ancient mammalian genomes are characterized by unusual architectural features with respect to G + C and repeat content, as well as compression relative to human. Approximately 14% and 34% of the human sequence forms alignments with the orthologous sequence from platypus and the marsupials, respectively; these numbers are distinctly lower than that observed with nonprimate eutherian mammals (45-70%). The alignable sequences between human and each marsupial species are not completely overlapping (only 80% common to all three species) nor are the platypus-alignable sequences completely contained within the marsupial-alignable sequences. Phylogenetic analysis of synonymous coding positions reveals that platypus has a notably long branch length, with the human-platypus substitution rate being on average 55% greater than that seen with human-marsupial pairs. Finally, analyses of the major mammalian lineages reveal distinct patterns with respect to the common presence of evolutionarily conserved vertebrate sequences. Our results confirm that genomic sequence from noneutherian mammals can contribute uniquely to unraveling the functional and evolutionary histories of the mammalian genome.

Mammalian phylogenomics comes of age

The relatively new field of phylogenomics is beginning to reveal the potential of genomic data for evolutionary studies. As the cost of whole genome sequencing falls, anticipation of complete genome sequences from divergent species, reflecting the major lineages of modern mammals, is no longer a distant dream. In this article, we describe how comparative genomic data from mammals is progressing to resolve long-standing phylogenetic controversies, to refine dogma on how chromosomes evolve and to guide annotation of human and other vertebrate genomes.

MAP2: multiple alignment of syntenic genomic sequences

We describe a multiple alignment program named MAP2 based on a generalized pairwise global alignment algorithm for handling long, different intergenic and intragenic regions in genomic sequences. The MAP2 program produces an ordered list of local multiple alignments of similar regions among sequences, where different regions between local alignments are indicated by reporting only similar regions. We propose two similarity measures for the evaluation of the performance of MAP2 and existing multiple alignment programs. Experimental results produced by MAP2 on four real sets of orthologous genomic sequences show that MAP2 rarely missed a block of transitively similar regions and that MAP2 never produced a block of regions that are not transitively similar. Experimental results by MAP2 on six simulated data sets show that MAP2 found the boundaries between similar and different regions precisely. This feature is useful for finding conserved functional elements in genomic sequences. The MAP2 program is freely available in source code form at http://bioinformatics.iastate.edu/aat/sas.html for academic use.

Animal phylogenomics: multiple interspecific genome comparisons

The utility of DNA sequence information for phylogenetics and phylogeography is now well known. Rather than attempt to summarize studies addressing this well-demonstrated utility, this chapter focuses on fundamental approaches and techniques that implement the collection of DNA sequence data for comparative phylogenetic purposes in a genomic context (phylogenomics). Whole genome sequencing approaches have changed the way we think about phylogenetics and have opened the way for new perspectives on "old" phylogenetics concerns. Some of these concerns are which gene regions to use and how much sequence information is needed for robust phylogenetic inference. Whole genome sequences of a few animal model organisms have gone a long way to implement approaches to better understand these important phylogenetic concerns. This chapter also addresses how genomics has made it more important for a clear understanding of orthology of gene regions in comparative biology. Finally, genome-enabled technologies that are affecting comparative biology are also discussed.

A globin gene of ancient evolutionary origin in lower vertebrates: evidence for two distinct globin families in animals

Hemoglobin, myoglobin, neuroglobin, and cytoglobin are four types of vertebrate globins with distinct tissue distributions and functions. Here, we report the identification of a fifth and novel globin gene from fish and amphibians, which has apparently been lost in the evolution of higher vertebrates (Amniota). Because its function is presently unknown, we tentatively call it globin X (GbX). Globin X sequences were obtained from three fish species, the zebrafish Danio rerio, the goldfish Carassius auratus, and the pufferfish Tetraodon nigroviridis, and the clawed frog Silurana tropicalis. Globin X sequences are distinct from vertebrate hemoglobins, myoglobins, neuroglobins, and cytoglobins. Globin X displays the highest identity scores with neuroglobin (approximately 26% to 35%), although it is not a neuronal protein, as revealed by RT-PCR experiments on goldfish RNA from various tissues. The distal ligand-binding and the proximal heme-binding histidines (E7 and F8), as well as the conserved phenylalanine CD1 are present in the globin X sequences, but because of extensions at the N-terminal and C-terminal, the globin X proteins are longer than the typical eight alpha-helical globins and comprise about 200 amino acids. In addition to the conserved globin introns at helix positions B12.2 and G7.0, the globin X genes contain two introns in E10.2 and H10.0. The intron in E10.2 is shifted by 1 bp in respect to the vertebrate neuroglobin gene (E11.0), providing possible evidence for an intron sliding event. Phylogenetic analyses confirm an ancient evolutionary relationship of globin X with neuroglobin and suggest the existence of two distinct globin types in the last common ancestor of Protostomia and Deuterostomia.

A 1463 gene cattle-human comparative map with anchor points defined by human genome sequence coordinates

A second-generation 5000 rad radiation hybrid (RH) map of the cattle genome was constructed primarily using cattle ESTs that were targeted to gaps in the existing cattle-human comparative map, as well as to sparsely populated map intervals. A total of 870 targeted markers were added, bringing the number of markers mapped on the RH(5000) panel to 1913. Of these, 1463 have significant BLASTN hits (E < e(-5)) against the human genome sequence. A cattle-human comparative map was created using human genome sequence coordinates of the paired orthologs. One-hundred and ninety-five conserved segments (defined by two or more genes) were identified between the cattle and human genomes, of which 31 are newly discovered and 34 were extended singletons on the first-generation map. The new map represents an improvement of 20% genome-wide comparative coverage compared with the first-generation map. Analysis of gene content within human genome regions where there are gaps in the comparative map revealed gaps with both significantly greater and significantly lower gene content. The new, more detailed cattle-human comparative map provides an improved resource for the analysis of mammalian chromosome evolution, the identification of candidate genes for economically important traits, and for proper alignment of sequence contigs on cattle chromosomes.

Modeling human hematopoietic and cardiovascular diseases in zebrafish

Zebrafish have emerged as a useful vertebrate model system in which unbiased large-scale screens have revealed hundreds of mutations affecting vertebrate development. Many zebrafish mutants closely resemble known human disorders, thus providing intriguing prospects for uncovering the genetic basis of human diseases and for the development of pharmacologic agents that inhibit or correct the progression of developmental disorders. The rapid pace of advances in genomic sequencing and map construction, in addition to morpholino targeting and transgenic techniques, have facilitated the identification and analysis of genes associated with zebrafish mutants, thus promoting the development of zebrafish as a model for human disorders. This review aims to illustrate how the zebrafish has been used to identify unknown genes, to assign function to known genes, and to delineate genetic pathways, all contributing valuable leads toward understanding human pathophysiology.

Comparative genomics: methods and applications

Interpreting the functional content of a given genomic sequence is one of the central challenges of biology today. Perhaps the most promising approach to this problem is based on the comparative method of classic biology in the modern guise of sequence comparison. For instance, protein-coding regions tend to be conserved between species. Hence, a simple method for distinguishing a functional exon from the chance absence of stop codons is to investigate its homologue from closely related species. Predicting regulatory elements is even more difficult than exon prediction, but again, comparisons pinpointing conserved sequence motifs upstream of translation start sites are helping to unravel gene regulatory networks. In addition to interspecific studies, intraspecific sequence comparison yields insights into the evolutionary forces that have acted on a species in the past. Of particular interest here is the identification of selection events such as selective sweeps. Both intra- and interspecific sequence comparisons are based on a variety of computational methods, including alignment, phylogenetic reconstruction, and coalescent theory. This article surveys the biology and the central computational ideas applied in recent comparative genomics projects. We argue that the most fruitful method of understanding the functional content of genomes is to study them in the context of related genomic sequences. In particular, such a study may reveal selection, a fundamental pointer to biological relevance.

Interspecies comparison of neuroglobin, cytoglobin and myoglobin: Sequence evolution and candidate regulatory elements

Neuroglobin and cytoglobin are two novel members of the vertebrate globin family. Their physiological role is poorly understood, although both proteins bind oxygen reversibly and may be involved in cellular oxygen homeostasis. Here we investigate the selective constraints on coding and non-coding sequences of the neuroglobin and cytoglobin genes in human, mouse, rat and fish. Neuroglobin and cytoglobin are highly conserved, displaying very low levels of non-synonymous nucleotide substitutions. An oxygen supply function predicts distinct modes of gene regulation, involving hypoxia-responsive transcription factors. To detect conserved candidate regulatory elements, we compared the neuroglobin and cytoglobin genes in mammals and fish. The myoglobin gene was included to test if it also contains hypoxia-responsive regulatory elements. Long conserved non-coding sequences, indicative of gene-regulatory elements, were found in the cytoglobin and myoglobin, but not in the neuroglobin gene. Sequence comparison and experimental data allowed us to delimit upstream regions of the neuroglobin and cytoglobin genes that contain the putative promoters, defining candidate regulatory regions for functional tests. The neuroglobin and the myoglobin genes both lack conserved hypoxia-responsive elements (HREs) for transcriptional activation, but contain conserved hypoxia-inducible mRNA stabilization signals in their 3' untranslated regions. The cytoglobin gene, in contrast, harbors both conserved HREs and mRNA stabilization sites, strongly suggestive of an oxygen-dependent regulation.

Projector: automatic contig mapping for gap closure purposes

Projector was designed for automatic positioning of contigs from an unfinished prokaryotic genome onto a template genome of a closely related strain or species. Projector mapped 84 contigs of Lactococcus lactis MG1363 (corresponding to 81% of the assembly nucleotides) against the genome of L.lactis IL1403. Ninety three percent of subsequent gap closure PCRs were successful. Moreover, a significant improvement in the N50 and N80 values (describing the assembly quality) was observed after the use of Projector. Because increasing numbers of bacterial genomes are being sequenced, Projector provides an efficient method to close a significant number of remaining gaps in the late stages of a genome sequencing project.

Seven evolutionarily conserved human rhodopsin G protein-coupled receptors lacking close relatives

We report seven new members of the superfamily of human G protein-coupled receptors (GPCRs) found by searches in the human genome databases, termed GPR100, GPR119, GPR120, GPR135, GPR136, GPR141, and GPR142. We also report 16 orthologues of these receptors in mouse, rat, fugu (pufferfish) and zebrafish. Phylogenetic analysis shows that these are additional members of the family of rhodopsin-type GPCRs. GPR100 shows similarity with the orphan receptor SALPR. Remarkably, the other receptors do not have any close relative among other known human rhodopsin-like GPCRs. Most of these orphan receptors are highly conserved through several vertebrate species and are present in single copies. Analysis of expressed sequence tag (EST) sequences indicated individual expression patterns, such as for GPR135, which was found in a wide variety of tissues including eye, brain, cervix, stomach and testis. Several ESTs for GPR141 were found in marrow and cancer cells, while the other receptors seem to have more restricted expression patterns.

A cattle-human comparative map built with cattle BAC-ends and human genome sequence

As a step toward the goal of adding the cattle genome to those available for multispecies comparative genome analysis, 40,224 cattle BAC clones were end-sequenced, yielding 60,547 sequences (BAC end sequences, BESs) after trimming with an average read length of 515 bp. Cattle BACs were anchored to the human and mouse genome sequences by BLASTN search, revealing 29.4% and 10.1% significant hits (E < e-5), respectively. More than 60% of all cattle BES hits in both the human and mouse genomes are located within known genes. In order to confirm in silico predictions of orthologyand their relative position on cattle chromosomes, 84 cattle BESs with similarity to sequences on HSA11 were mapped using a cattle-hamster radiation hybrid (RH) panel. Resulting RH maps of BTA15 and BTA29 cover approximately 85% of HSA11 sequence, revealing a complex patchwork shuffling of segments not explained by a simple translocation followed by internal rearrangements. Overlay of the mouse conserved syntenies onto HSA11 revealed that segmental boundaries appear to be conserved in all three species. The BAC clone-based comparative map provides a foundation for the evolutionary analysis of mammalian karyotypes and for sequencing of the cattle genome.

MultiPipMaker and supporting tools: Alignments and analysis of multiple genomic DNA sequences

Analysis of multiple sequence alignments can generate important, testable hypotheses about the phylogenetic history and cellular function of genomic sequences. We describe the MultiPipMaker server, which aligns multiple, long genomic DNA sequences quickly and with good sensitivity (available at http://bio.cse.psu.edu/ since May 2001). Alignments are computed between a contiguous reference sequence and one or more secondary sequences, which can be finished or draft sequence. The outputs include a stacked set of percent identity plots, called a MultiPip, comparing the reference sequence with subsequent sequences, and a nucleotide-level multiple alignment. New tools are provided to search MultiPipMaker output for conserved matches to a user-specified pattern and for conserved matches to position weight matrices that describe transcription factor binding sites (singly and in clusters). We illustrate the use of MultiPipMaker to identify candidate regulatory regions in WNT2 and then demonstrate by transfection assays that they are functional. Analysis of the alignments also confirms the phylogenetic inference that horses are more closely related to cats than to cows.

Molecular evolution meets the genomics revolution

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.