zPicture: dynamic alignment and visualization tool for analyzing conservation profiles

Genomic deletion of a long-range bone enhancer misregulates sclerostin in Van Buchem disease

Mutations in distant regulatory elements can have a negative impact on human development and health, yet because of the difficulty of detecting these critical sequences, we predominantly focus on coding sequences for diagnostic purposes. We have undertaken a comparative sequence-based approach to characterize a large noncoding region deleted in patients affected by Van Buchem (VB) disease, a severe sclerosing bone dysplasia. Using BAC recombination and transgenesis, we characterized the expression of human sclerostin (SOST) from normal (SOST(wt)) or Van Buchem (SOST(vbDelta) alleles. Only the SOST(wt) allele faithfully expressed high levels of human SOST in the adult bone and had an impact on bone metabolism, consistent with the model that the VB noncoding deletion removes a SOST-specific regulatory element. By exploiting cross-species sequence comparisons with in vitro and in vivo enhancer assays, we were able to identify a candidate enhancer element that drives human SOST expression in osteoblast-like cell lines in vitro and in the skeletal anlage of the embryonic day 14.5 (E14.5) mouse embryo, and discovered a novel function for sclerostin during limb development. Our approach represents a framework for characterizing distant regulatory elements associated with abnormal human phenotypes.

Distinct transcriptional regulation and function of the human BACE2 and BACE1 genes

Amyloid beta protein (Abeta) is the principal component of neuritic plaques in Alzheimer's disease (AD). Abeta is derived from beta amyloid precursor protein (APP) by beta- and gamma-secretases. Beta-site APP cleaving enzyme 1 (BACE1) has been identified as the major beta-secretase. BACE2 is the homolog of BACE1. The BACE2 gene is on chromosome 21 and has been implicated in the pathogenesis of AD. However, the function of BACE2 in Abeta generation is controversial. Some studies have shown that BACE2 cleaved APP at the beta-site whereas other studies showed it cleaved around the alpha-secretase site. To elucidate the involvement of BACE2 in AD pathogenesis, we compared BACE2 and BACE1 gene regulation and their functions in Abeta generation. We cloned and functionally characterized the human BACE2 promoter. The BACE2 gene is controlled by a TATA-less promoter. Though Sp1 can regulate both BACE1 and BACE2 genes, comparative sequence analysis and transcription factor prediction showed little similarity between the two promoters. BACE1 increased APP cleavage at the beta-site and Abeta production whereas BACE2 did not. Overexpression of BACE2 significantly increased sAPP levels in conditioned media but markedly reduced Abeta production. Knockdown of BACE2 resulted in increased APP C83. Our data indicate that despite being homologous in amino acid sequence, BACE2 and BACE1 have distinct functions and transcriptional regulation. BACE2 is not a beta-secretase, but processes APP within the Abeta domain at a site downstream of the alpha-secretase cleavage site. Our data argue against BACE2 being involved in the formation of neuritic plaques in AD.

Evolution of transcriptional enhancers in the immunoglobulin heavy-chain gene: functional characteristics of the zebrafish Emu3' enhancer

Past studies of the channel catfish immunoglobulin heavy-chain (IgH) locus indicates that it lacks an Emu enhancer in the J(H)-Cmu1 intron but does have an enhancer, termed Emu3', in the mu-delta intergenic region. The positioning of the catfish enhancer downstream of the mu-chain exons is predicted to be unfavorable for antibody-affinity maturation in catfish, and would also have been an impediment to the evolution of class switch recombination, had it existed in early tetrapods. To determine if this downstream enhancer is a general feature of teleost fish, we have identified the location of the transcriptional enhancer in the zebrafish IgH locus. We find that zebrafish, like catfish, only have an Emu3'-like enhancer that has cross-species activity, but which is B-cell-specific in its activity. A 300-bp region of the zebrafish enhancer shares sequence homology with the core of the catfish Emu3', although there has been loss and gain of specific octamer enhancer motifs. Mutagenesis studies demonstrate that the zebrafish IgH enhancer depends on a pair of E-box motifs that are found in the enhancer core. Similarly spaced E-box motifs appear to exist in the Emu3' enhancer regions of other teleost fish, suggesting this is a common feature among fish IgH enhancers. We discuss how this distal positioning of the enhancer may influence affinity maturation in extant teleosts as well as the evolution of this process in the early tetrapods.

Improvements to GALA and dbERGE II: databases featuring genomic sequence alignment, annotation and experimental results

We describe improvements to two databases that give access to information on genomic sequence similarities, functional elements in DNA and experimental results that demonstrate those functions. GALA, the database of Genome ALignments and Annotations, is now a set of interlinked relational databases for five vertebrate species, human, chimpanzee, mouse, rat and chicken. For each species, GALA records pairwise and multiple sequence alignments, scores derived from those alignments that reflect the likelihood of being under purifying selection or being a regulatory element, and extensive annotations such as genes, gene expression patterns and transcription factor binding sites. The user interface supports simple and complex queries, including operations such as subtraction and intersections as well as clustering and finding elements in proximity to features. dbERGE II, the database of Experimental Results on Gene Expression, contains experimental data from a variety of functional assays. Both databases are now run on the DB2 database management system. Improved hardware and tuning has reduced response times and increased querying capacity, while simplified query interfaces will help direct new users through the querying process. Links are available at http://www.bx.psu.edu/.

Mulan: multiple-sequence local alignment and visualization for studying function and evolution

Multiple-sequence alignment analysis is a powerful approach for understanding phylogenetic relationships, annotating genes, and detecting functional regulatory elements. With a growing number of partly or fully sequenced vertebrate genomes, effective tools for performing multiple comparisons are required to accurately and efficiently assist biological discoveries. Here we introduce Mulan (http://mulan.dcode.org/), a novel method and a network server for comparing multiple draft and finished-quality sequences to identify functional elements conserved over evolutionary time. Mulan brings together several novel algorithms: the TBA multi-aligner program for rapid identification of local sequence conservation, and the multiTF program for detecting evolutionarily conserved transcription factor binding sites in multiple alignments. In addition, Mulan supports two-way communication with the GALA database; alignments of multiple species dynamically generated in GALA can be viewed in Mulan, and conserved transcription factor binding sites identified with Mulan/multiTF can be integrated and overlaid with extensive genome annotation data using GALA. Local multiple alignments computed by Mulan ensure reliable representation of short- and large-scale genomic rearrangements in distant organisms. Mulan allows for interactive modification of critical conservation parameters to differentially predict conserved regions in comparisons of both closely and distantly related species. We illustrate the uses and applications of the Mulan tool through multispecies comparisons of the GATA3 gene locus and the identification of elements that are conserved in a different way in avians than in other genomes, allowing speculation on the evolution of birds. Source code for the aligners and the aligner-evaluation software can be freely downloaded from http://www.bx.psu.edu/miller_lab/.

Comparative genomics

The genomes from three mammals (human, mouse, and rat), two worms, and several yeasts have been sequenced, and more genomes will be completed in the near future for comparison with those of the major model organisms. Scientists have used various methods to align and compare the sequenced genomes to address critical issues in genome function and evolution. This review covers some of the major new insights about gene content, gene regulation, and the fraction of mammalian genomes that are under purifying selection and presumed functional. We review the evolutionary processes that shape genomes, with particular attention to variation in rates within genomes and along different lineages. Internet resources for accessing and analyzing the treasure trove of sequence alignments and annotations are reviewed, and we discuss critical problems to address in new bioinformatic developments in comparative genomics.

rVISTA 2.0: evolutionary analysis of transcription factor binding sites

Identifying and characterizing the transcription factor binding site (TFBS) patterns of cis-regulatory elements represents a challenge, but holds promise to reveal the regulatory language the genome uses to dictate transcriptional dynamics. Several studies have demonstrated that regulatory modules are under positive selection and, therefore, are often conserved between related species. Using this evolutionary principle, we have created a comparative tool, rVISTA, for analyzing the regulatory potential of noncoding sequences. Our ability to experimentally identify functional noncoding sequences is extremely limited, therefore, rVISTA attempts to fill this great gap in genomic analysis by offering a powerful approach for eliminating TFBSs least likely to be biologically relevant. The rVISTA tool combines TFBS predictions, sequence comparisons and cluster analysis to identify noncoding DNA regions that are evolutionarily conserved and present in a specific configuration within genomic sequences. Here, we present the newly developed version 2.0 of the rVISTA tool, which can process alignments generated by both the zPicture and blastz alignment programs or use pre-computed pairwise alignments of several vertebrate genomes available from the ECR Browser and GALA database. The rVISTA web server is closely interconnected with the TRANSFAC database, allowing users to either search for matrices present in the TRANSFAC library collection or search for user-defined consensus sequences. The rVISTA tool is publicly available at http://rvista.dcode.org/.

ECR Browser: a tool for visualizing and accessing data from comparisons of multiple vertebrate genomes

With an increasing number of vertebrate genomes being sequenced in draft or finished form, unique opportunities for decoding the language of DNA sequence through comparative genome alignments have arisen. However, novel tools and strategies are required to accommodate this large volume of genomic information and to facilitate the transfer of predictions generated by comparative sequence alignment to researchers focused on experimental annotation of genome function. Here, we present the ECR Browser, a tool that provides easy and dynamic access to whole genome alignments of human, mouse, rat and fish sequences. This web-based tool (http://ecrbrowser.dcode.org) provides the starting point for discovery of novel genes, identification of distant gene regulatory elements and prediction of transcription factor binding sites. The genome alignment portal of the ECR Browser also permits fast and automated alignments of any user-submitted sequence to the genome of choice. The interconnection of the ECR Browser with other DNA sequence analysis tools creates a unique portal for studying and exploring vertebrate genomes.

eShadow: a tool for comparing closely related sequences

Primate sequence comparisons are difficult to interpret due to the high degree of sequence similarity shared between such closely related species. Recently, a novel method, phylogenetic shadowing, has been pioneered for predicting functional elements in the human genome through the analysis of multiple primate sequence alignments. We have expanded this theoretical approach to create a computational tool, eShadow, for the identification of elements under selective pressure in multiple sequence alignments of closely related genomes, such as in comparisons of human-to-primate or mouse-to-rat DNA. This tool integrates two different statistical methods and allows for the dynamic visualization of the resulting conservation profile. eShadow also includes a versatile optimization module capable of training the underlying Hidden Markov Model to differentially predict functional sequences. This module grants the tool high flexibility in the analysis of multiple sequence alignments and in comparing sequences with different divergence rates. Here, we describe the eShadow comparative tool and its potential uses for analyzing both multiple nucleotide and protein alignments to predict putative functional elements.