Conservation of core gene expression in vertebrate tissues
© Chan et al.; licensee BioMed Central Ltd. 2009
Received: 23 January 2009
Accepted: 18 March 2009
Published: 16 April 2009
Vertebrates share the same general body plan and organs, possess related sets of genes, and rely on similar physiological mechanisms, yet show great diversity in morphology, habitat and behavior. Alteration of gene regulation is thought to be a major mechanism in phenotypic variation and evolution, but relatively little is known about the broad patterns of conservation in gene expression in non-mammalian vertebrates.
We measured expression of all known and predicted genes across twenty tissues in chicken, frog and pufferfish. By combining the results with human and mouse data and considering only ten common tissues, we have found evidence of conserved expression for more than a third of unique orthologous genes. We find that, on average, transcription factor gene expression is neither more nor less conserved than that of other genes. Strikingly, conservation of expression correlates poorly with the amount of conserved nonexonic sequence, even using a sequence alignment technique that accounts for non-collinearity in conserved elements. Many genes show conserved human/fish expression despite having almost no nonexonic conserved primary sequence.
There are clearly strong evolutionary constraints on tissue-specific gene expression. A major challenge will be to understand the precise mechanisms by which many gene expression patterns remain similar despite extensive cis-regulatory restructuring.
Vertebrates all share a body plan, gene number and gene catalog [1–4] inherited from a common progenitor, but so far it has been unclear to what degree gene expression is conserved. King and Wilson  initially posited that phenotypic differences among primates are mainly due to adaptive changes in gene regulation, rather than to changes in protein-coding sequence or function, and this idea has accumulated supporting evidence in recent years [6–12]. Recent work has indicated that gene expression evolves in a fashion similar to other traits, where in the absence of selection, random mutations introduce variants within a population [11, 13–19]. Changes negatively affecting fitness are probably eliminated by purifying selection: core cellular processes seem to be coexpressed from yeast to human , and conservation of the expression of individual genes in specific tissues has been observed across distantly related vertebrates [21–24], perhaps reflecting requirements for patterning and development as well as conserved functions of organs, tissues and cell types. Conversely, changes that benefit fitness (for example, under new ecological pressures) may become fixed: changes in gene expression are believed to underlie many differences in morphology, physiology and behavior and, indeed, subtle differences in gene regulation can result in spatial and temporal alterations in transcript levels, with phenotypic consequences at the cell, tissue and organismal levels [5, 25]. The degree to which stabilizing selection constrains directional selection and neutral drift across the full vertebrate subphylum is, to our knowledge, unknown.
Comparative genomic analyses provide a perspective on the evolution of both cis- and trans-regulatory mechanisms, and they are often used as a starting point for the identification of regulatory mechanisms. One estimate, using collinear multiple-genome alignments, suggested that roughly a million sequence elements are conserved in vertebrates (particularly among mammals, which represent the majority of sequenced vertebrates) [26–29], with most being nonexonic , and a series of studies have demonstrated the cis-regulatory potential of the most highly conserved nonexonic elements (for example, [27, 29, 30]). Another study  found that only 29% of nonexonic mammalian conserved bases are evident in chicken, and that nearly all aligning sequence in fish overlaps exons, raising the possibility that gene regulatory mechanisms may be very different among vertebrate clades. Absence of conserved sequence does not imply lack of regulatory conservation, however, as many known cis-regulatory elements seem to undergo rapid turnover [32, 33], and there are examples in which orthologous genes have similar expression patterns despite apparent lack of sequence conservation in regulatory regions . As further evidence of pervasive regulatory restructuring in vertebrate evolution, an analysis  that accounted for shuffling (non-collinearity) of locally conserved sequences suggested that the number of conserved elements may be several fold higher than collinear alignments detect, particularly between distant vertebrate relatives, such as mammals and fish.
Trans-acting factors (transcription factors or TFs) also show examples of striking conservation, such as among the homeotic factors, and diversifying selection . Studies comparing expression patterns between human and chimpanzee liver found that TF genes were enriched among the genes with greatest human-specific increase in expression levels [37, 38], supporting arguments for alteration of trans-regulatory architecture as a driving evolutionary mechanism . On the other hand, in the Drosophila developmental transition, expression of transcription factor genes is more evolutionarily stable than expression of their targets, on average . The fact that enhancers will often function similarly in fish and mammals, even when the enhancer itself is not conserved, indicates that mechanisms underlying cell-specific and developmental expression are likely to be widely conserved across vertebrates [41, 42].
Global trends in conservation of gene expression, conservation of cis-regulatory sequence and relationships between the two are not completely understood [13, 39, 41], partly because the cis-regulatory 'lexicon' (that is, how TF binding sites combine to form enhancers) remains mostly unknown, testing individual enhancers is tedious and expensive, and many vertebrates are not amenable to genetic experimentation. These issues are of both academic and practical consequence: in addition to our curiosity about the origin and distinctive characteristics of the human species, primary sequence conservation is widely used to identify regulatory mechanisms. We reasoned that expression profiling data from species spanning much greater phylogenetic distance than humans and mice, and thus having greater opportunity for both neutral drift and positive selection, would allow assessment of the degree of conservation of tissue gene expression among all vertebrates, and a comparison of the conservation of expression to the conservation of nonexonic primary sequence. Here, we describe a survey of gene expression in adult tissues and organs in the main vertebrate clades: mammals, avians/reptiles, amphibians and fish. Our analyses demonstrate that core tissue-specific gene expression patterns are conserved across all major vertebrate lineages, but that the correspondence between conservation of expression and amount of conserved nonexonic sequence is weak overall, at least at a level that is detectable by current alignment approaches.
Tissue-specific gene expression is broadly correlated across vertebrates
To examine gene expression in a broad range of vertebrates, we collected a compendium of gene expression datasets, consisting of previously published datasets for human  and mouse , and newly generated datasets containing 20 tissues each from chicken (Gallus gallus), frog (Xenopus tropicalis) and pufferfish (Tetraodon nigroviridis). Details of the experiments are found in the Materials and methods; lists of tissues are found in Additional data file 1. Clustering analyses of each dataset separately (Additional data file 2) shows that prominent tissue-specific expression patterns are found in all vertebrates.
Thousands of individual tissue-specific gene expression events are conserved across all vertebrate clades
Conservation of expression does not correlate with proportion or amount of conserved nonexonic sequence
We next asked what gene properties correlate with conservation of expression among the 3,074 measured unique orthologs. We considered the following gene properties: those that are contained in our data, that is, median expression level and Shannon entropy as a measure of tissue specificity and preferential expression in individual tissues; GO annotations; and sequence properties, that is, length of gene, size of encoded protein, presence of a DNA-binding domain (for known and predicted TFs), sequence conservation of encoded protein (pairwise BLASTP bit score) and amount of conserved nonexonic sequence (measured in several ways) (Additional data files 7 and 8; see Materials and methods for details).
Several observations emerged from this analysis. First, the genes with the highest expression similarity between species are most often genes expressed in a highly tissue-specific manner in tissues with specialized functions. Although the Pearson correlation is heavily influenced by extreme values, thus giving higher weight to tissue-specific pairs, most of these high scoring genes were also classified as conserved by our binary measure. Among the 50 genes with highest median pairwise Pearson correlation of expression are structural components of the eye lens, liver-synthesized proteins involved in the complement system and blood coagulation, and neurotransmitter receptors and transporters. This observation is supported by the GO categories enriched among genes with high expression similarity, such as synaptic transmission (GO:0007268), visual perception (GO:0007601), wound healing (GO:0042060) and muscle development (GO:0007517) (Wilcoxon-Mann-Whitney test (WMW) p-values 1.55 × 10-4, 2.36 × 10-3, 2.24 × 10-3 and 4.98 × 10-5, respectively; Additional data file 8). In contrast, we did not find any evidence that the expression of TFs (228 of the 3,074 measured orthologs) is more or less conserved than that of non-TFs, in contrast to previous reports of both higher  and lower  rates of evolution of TF expression. A slightly lower proportion of TFs did seem to show conservation events relative to non-TFs using the binary measure, but this difference is due to the fact that TFs are expressed in fewer tissues: the difference is not seen when comparing TFs and non-TFs with similar overall expression levels (data not shown).
We reasoned that pervasive shuffling might obscure most of the cis-regulatory elements, particularly in pufferfish. In order to address this possibility, we developed a technique similar to that of Sanges et al.  to detect shuffled conserved sequence elements (SCEs), which may be non-collinear, across the five species (see Materials and methods for details). Among the total 4,898 1-1-1-1-1 orthologs, we identified 491,028, 457,074, 79,001, 54,134 and 11,731 SCEs in human, mouse, chicken, frog and pufferfish with median lengths of 164, 80, 68, 68 and 65 nucleotides, respectively. These SCEs showed good overlap with those in  (75.5% of the sequences in  within regions we aligned were identified as SCEs in our analysis) and they were calibrated to minimize false positives (see Materials and methods). However, we still did not observe a strong relationship between the degree of conservation and the proportion or number of aligned bases in each species (median Spearman correlation: -0.062 and 0.042 for binary and Pearson models, respectively, versus proportion of aligned nonexonic bases in each species; Figure 6c,d; similar correlations are obtained with number of aligned non-exonic bases).
Because TF binding sites are degenerate, it is conceivable that these genes have a high number of conserved TF binding sites, despite their lack of primary sequence conservation. To examine this possibility we used Enhancer Element Locator (EEL)  to align TF binding sites defined by 138 motif models downloaded from the JASPAR database . Over all 4,804 aligned human/pufferfish ortholog pairs, the number of genes that scored highly using EEL was only slightly higher with real ortholog pairs than with randomly assigned orthologs with similar amounts of nonexonic associated sequence in both genomes (p = 0.24, Kolmogorov-Smirnov test; see Materials and methods and Additional data file 9) and there is almost no correlation between EEL score and conservation of expression (EEL score against median versus pufferfish normalized intensity Pearson r = 0.022). We conclude that the regulatory architecture of the vast majority of genes has diverged beyond recognition by any current approaches, despite the apparently very similar regulatory output in many cases, and the likelihood that at least some orthologous TFs are functioning in the same tissues.
Our data provide a resource of large-scale gene expression data in tissues of three non-mammalian vertebrates and demonstrate that there is conservation of core vertebrate tissue gene expression. Our analysis almost certainly underestimates the proportion of genes with conserved expression patterns, because our analysis focused on only ten large adult tissues in mature animals in captivity. Nonetheless, our results already provide an index of several thousand highly conserved tissue gene expression events and a picture of core gene expression in major tissues and organs of the common progenitor of all vertebrates, which most likely resembled the expression patterns shown in Figure 4.
In our analysis, some biological processes emerged as more highly conserved than others. Genes involved in the more conserved processes on the whole tend to be preferentially expressed in tissues with a limited number of cell types (brain, eye, liver, heart and muscle) that carry out specialized functions particular to those tissues. This finding is consistent with the notion that mechanisms underlying important biological processes should be conserved across taxa . It is likely that the conservation of gene expression in these tissues extends beyond the base of vertebrates; coexpression of neuronal genes, for example, is observed as far as nematodes . Genes expressed in tissues subject to greater environmental influence (such as intestine, stomach and spleen) may be more likely to take on new roles and diverge in expression as means of adaptation. We find gene expression similarity in the testis across vertebrates to be relatively low compared with other tissues, in support of earlier observations of accelerated evolution of testis gene expression in Drosophila  and primates , and consistent with mating competition and speciation. We also note that relatively low conservation of gene expression in the kidney is consistent with its divergent function; in teleosts, the kidney is a major lymphoid organ .
Our finding that the correlation between the amount of conserved nonexonic sequence and conservation of gene expression is low underscores the apparent malleability of the cis-regulatory 'lexicon' [32, 35, 57]. It is easy to rationalize cases in which there is conserved sequence but little or no conservation of expression in our study, because it is possible that the majority of conserved sequence identified is either not cis-regulatory, is functioning as cis-regulatory in a context that we have not measured and/or regulates neighboring genes. What is most striking is that many genes with the highest conserved tissue-specific expression have almost no nonexonic conserved sequence outside mammals. Divergence in trans-regulatory architecture does not provide a satisfying explanation for this observation: although there are examples in which the binding specificity of TFs evolves [58, 59], as a general rule the individual monomeric TF sequence specificities seem to be unchanged over large evolutionary distances , and DNA-contacting residues are often the most conserved . The fact that conservation of TF expression is comparable to that of other genes with similar levels of tissue specificity also supports the notion that many ancestral vertebrate regulatory mechanisms are still in use. Moreover, enhancers can be functional across large evolutionary distances (for example, human reporters in zebrafish), even when the enhancers are not conserved, or are below the threshold of detection by current alignment techniques .
Understanding the mechanisms underlying conservation of expression patterns is a major challenge in our understanding of evolution, and of genome function and gene regulation: even within a single genome, it is difficult to find cis-regulatory modules shared by coexpressed genes, indicating that there are many ways to achieve the same expression output. We propose that our catalog of conserved (and non-conserved) expression will be useful to test ideas regarding enhancer definition. In particular, we predict that the small size of pufferfish genes and knowledge of the expression of TFs in each tissue may facilitate searches for enhancers on the basis of density and arrangement of TF binding sites, rather than primary sequence alignment [34, 62]. These and other techniques require experimental data for training and testing, and we now provide such data for tens of thousands of genes across the vertebrate lineage, including several thousand unique orthologs.
Materials and methods
Data tables including microarray data, probe sequences, the 1-1-1-1-1 ortholog list, Phastcons information, UCE information, GO annotations and novel non-collinear alignments are found in the Additional data files and also on our project website .
Five chickens (three young females, one in-lay adult female and one rooster) were obtained from Sunnybrook Health Sciences Centre, Toronto, Canada. The bursa of Fabricius and thymus were dissected from the young females. The oviduct and ovaries were dissected from the in-lay adult female and the testis dissected from the rooster. Approximately 75 adult female and five adult male frogs were obtained from Nasco (Atkinson, USA), housed in 20 L aquariums at 20°C and fed a diet of Aquamax Grower 600 trout food (Purina Mills, Gray Summit, USA). Exactly 100 green spotted pufferfish, of unknown age and sex, described as caught wild in Malaysia, were purchased from Aquarium Services Warehouse Outlets (Thornhill, Canada) and housed in a single 120 gallon (453 L) aquarium tank at room temperature. Their markings, size (about 6–7 cm) and behavior matched descriptions of Tetraodon nigroviridis. The fish were fed frozen brine shrimp twice a day.
Tissue collection and mRNA isolation
Animal handling and euthanization procedures were performed according to protocols approved by the University of Toronto. Tissue processing and mRNA extraction were as described previously . Tissues collected are listed in Additional file 1. All fish and frogs were anaesthetized by adding 35 ml of a 1:10 clove oil (Rougier Pharma, Mirabel, Canada):ethanol mixture to 1.5 L tank water containing a single animal. Following cessation of movement (15–20 s), they were euthanized by decapitation. Chickens were euthanized by barbiturate injection followed by decapitation.
Microarray design, processing and data collection
Oligonucleotide microarrays (60-mer) were custom-designed by us and manufactured by Agilent Technologies to survey the expression of a set of known and predicted mRNAs compiled from Ensembl  version 37 and JGI (for X. tropicalis, v3.0) corresponding to 24,877, 25,594 and 25,937 known and predicted loci in Tetraodon, chicken and X. tropicalis, respectively, using 41,533, 41,534 and 41,523 probes. cDNA synthesis, labeling reactions and microarray processing procedures were performed as described previously . Each tissue was assayed in duplicate, in fluor-reversed pairs with biological replicates and averaged (average Pearson correlation between replicate arrays = 0.89). Using a previously described procedure to define expression , 72%, 77% and 89% of genes on the arrays were detected in at least one of the tissues profiled in chicken, frog and pufferfish, respectively, with 40%, 33% and 52% of genes detected in each tissue, on average (Additional data file 10), comparable to previously reported figures of 52% and 59% in human and mouse [44, 65]. Our percentages may be lower because many predicted genes were included on the arrays. Microarray-based expression data for 55 mouse tissues  and 60 human tissues  were downloaded and mapped to Ensembl v37 identifiers by aligning probe sequences to Ensembl transcript sequences using BLAT . Expression similarity was calculated using the relative expression levels or 'ratios' (median-subtracted, variance stabilization-normalized asinh intensities) as well as using the absolute values (normalized intensities) in each tissue, as indicated in the Results and figures. All new microarray data have been uploaded to the Gene Expression Omnibus (GEO) database, accession numbers [GEO:GSE12974, GEO:GSE12975, GEO:GSE12976], and are also found on our project website .
Determination of gene orthology relationships
Inparanoid  was used to analyze each possible pairwise all-versus-all protein BLAST  comparison between all known proteins for each of the five vertebrate species in version 37 of the Ensembl database to delineate pairwise gene orthology relationships. These relationships were then assembled into unique 1-1-1-1-1 ortholog groups across the five species using custom Perl scripts in an approach analogous to that of Alexeyenko et al. .
Gene Ontology analysis
GO annotations were downloaded from Ensembl BioMart  for each species. Annotations for chicken, frog and pufferfish (Tetraodon) were further supplemented by mapping the corresponding mouse annotations by any type of orthology as defined by Ensembl. Annotations were up-propagated and terms with few or too many annotated genes were removed as described previously . All GO WMW analyses performed across the set of unique ortho-logs were done with human annotations.
Definition of expression conservation events in each gene and tissue using the binary measure
We split each set of orthologs in each tissue in each species on the basis of their measured normalized expression intensities according to the following thresholds: top 1/2, top 1/3, top 1/4, top 1/5 and top 1/6. Because the human and mouse datasets were designed independently from those of the other three species, there were orthologous genes with missing measurements. In order to facilitate comparison between as many unique orthologs between the five species as possible, we applied variance stabilizing normalization  to the human and mouse orthologs in order to make them comparable, under the assumption that gene expression is mostly conserved between the mammals relative to the rest of the vertebrate phylogeny. We combined the human and mouse ortholog data by averaging to obtain a set of 3,074 unique orthologs with measured microarray data across ten common tissues.
Calculating Shannon entropy as a measure of tissue specificity
Shannon entropy, which measures the degree of overall tissue specificity of a gene, was calculated as described by Schug et al. . Briefly, the relative expression of a gene g in a tissue t relative to its expression given in N tissues is defined as pt|g= wg, t/Σ1 ≤ t ≤ Nwg, t, where wg, tis the expression level of the gene g in tissue t. The Shannon entropy of a gene's expression distribution is then calculated as H g = Σ1 ≤ t ≤ N-pt|glog2pt|g. This value is expressed in bits and ranges from zero to log2(10) genes expressed in a single tissue and uniformly expressed in all the common tissues examined, respectively.
Definition of overlap between nonexonic regions associated with human 1-1-1-1-1 orthologs and Phastcons elements and UCEs
Locations of Phastcons elements  and UCEs  were downloaded from the hg17 version of the UCSC genome browser  and from the supplementary website  of Bejerano et al. . The number of Phastcons elements and UCEs (and the number of bases) that overlapped all nonexonic sequences, including intergenic sequences, up to 50 kb upstream and downstream of each human 1-1-1-1-1 ortholog was tabulated using custom Perl scripts. The proportion of Phastcons and UCE bases covered in non-coding regions was calculated as the number of bases out of 50 kb of flanking bases upstream and downstream of each gene with coding regions masked out.
Definition and calculation of other gene features and statistical comparisons to expression conservation
A measure of protein sequence conservation between two species was derived by performing pairwise BLASTP  between the protein-coding sequences (downloaded from Ensembl EnsMart ) of an orthologous gene pair and retaining the bit score. The median bit score was taken as a measure of protein sequence conservation over all species. Average and maximum expression level and the expression rank within tissues were calculated in the expected manner for each gene across the ten common tissues. The number of bases in an aligned conserved element (aligned as described below) was obtained by summing the number of bases in each species within a five-way gapped alignment between all species. The total number of aligned bases in an aligned conserved element is the sum of counts in each species. TF genes within the set of unique orthologs were defined by the presence of a DNA-binding domain in the mouse protein sequence in the Pfam database . Our list of TF genes is given in Additional data file 11. Conservation of expression was measured by either our binary or Pearson measures, both of which yield a real value (an integer between zero and six in the case of the binary method). With the exception of GO annotations and TF identities, all measures for each gene are compiled in Additional data file 7; TFs are listed separately in Additional data file 11. Comparisons of properties with real values were made by Spearman p-value, and these are also given in Additional data file 7 and shown graphically in Additional data file 12. Comparisons of categorical properties were made by WMW p-value and are given in Additional data file 8.
Multiple sequence alignment algorithm
All repeat-masked intronic, 3' untranslated region and intergenic non-coding sequence upstream and downstream of protein-coding sequences in orthologous groups that we have identified using the Inparanoid algorithm  were downloaded from Ensembl v37 . Using LAGAN , we built a multiple global alignment across all of the genomic sequences within each set of unique orthologs to identify conserved non-coding elements in all species. We used up to 50 kb of upstream and downstream sequence, up to the point that the transcript of another annotated gene was encountered.
We initially applied a conservation cutoff of 55% in a 50 nucleotide window to search for conserved regions between the human and the orthologous genomes. After removing conserved elements that were annotated as exonic, we extracted each sequence element in the alignment that was conserved only in a subset of the genomes and built the most parsimonious ancestral reconstruction of each of the sequences using Fitch's algorithm, treating the gap character as a fifth symbol. This was then used to search against the other orthologous genome(s) using the CHAOS aligner  with very sensitive parameters. The Smith-Waterman threshold was varied between 40 and 60 (60 is the default conservation setting for the BLASTZ  alignment program; this setting is both relatively sensitive and very specific, and thus almost all hits above this threshold will be real). This homology filtering step was used to identify non-collinear conserved sequences that may have changed position and orientation relative to exons over evolutionary time. In contrast to the method used by Sanges et al. , we tried using both the mouse and the chicken sequence as the homology filter. By using the chicken as the base organism we were able to significantly lower the false positive rate (only one decoy (i.e. permuted gene identity) sequence with a score above 45), while our true positive rate was unchanged (all of the human/frog and human/pufferfish conserved regions recovered with mouse as the base were also recovered with pufferfish (Tetraodon) as the base).
Through the procedure described above, our approach takes advantage of the conservation of order of the conserved elements when no rearrangements have taken place, and the flexibility of aligning less conserved regions that have been shuffled around in the genome. By running Fitch's algorithm on the aligned sequences (unlike the Sanges et al.  approach of simply using the mouse sequence as input to CHAOS) we increased the sensitivity of our alignment technique. In particular, it was easier to align Tetraodon sequence to the common ancestor of human and chicken genomes than to either genome individually.
Enhancer Element Locator
We scanned the nonexonic sequence associated (as defined for the multiple sequence alignment algorithm) with 4,804 human/fish ortholog pairs for conserved TF binding sites by applying EEL  using the default parameters to perform a local pairwise TF binding site alignment. We did not attempt to align the 94 (of 4,898) ortholog pairs for which one of the orthologs appears in the intron of another gene. The motif models input into EEL were the 138 models returned by JASPAR webserver  on 27 February 2009. The score of only the best alignment from each orthologous gene group was captured and was used to construct a distribution of 4,804 gene group alignment scores. We also constructed a shuffled (negative control) distribution of EEL scores by attempting to align the TF binding sites between the human non-coding sequences in each of the ortholog pairs and those of six non-orthologous Tetraodon genes. These genes were selected among the other 4,803 non-orthologous paired Tetraodon genes to be the six with the most similar amount of non-coding sequence as the proper ortholog, under the constraint that three of the genes had to have less non-coding sequence and three had to have more non-coding sequence. As evidence that the EEL analysis is detecting some very limited degree of conservation, we note that the distribution of EEL scores is slightly higher in real than in randomly assigned orthologs starting around an EEL score of 175; in particular, we estimate that 1.3% (63 of 4,084) of all EEL scores are non-random by subtracting the proportion of random scores above 175 (39.8%) from the proportion of real scores above 175 (41.1%).
Additional data files
The following additional data files are available. Additional data file 1: Tissues analyzed in this study. The tissues at the top, highlighted in color, are those considered to be among the ten common tissue types. Those with identical coloring were combined (by averaging normalized intensities) for the analysis of conservation of gene expression among the ten common tissues. Additional data file 2: Microarray gene expression data obtained in this study. Clustergrams show the microarray datasets in chicken, frog and pufferfish, displayed as relative expression ratio (see Materials and methods) of each gene within each of the 20 tissues profiled. Rows and columns were ordered independently for each dataset, and high-level branches broken and rearranged to obtain a diagonal appearance as described in . Additional data file 3: Dendrogram of correlations among ten common tissues, using 1 – Pearson correlation and average linkage over 3,074 genes. Additional data file 4: Gene Ontology categories that tend to be expressed highly in each of the ten common tissues. Selected GO biological process categories enriched amongst genes highly expressed within each of the ten common tissues in each species are shown. The tissue and GO category order were manually arranged in the heat map. (A full matrix of WMW scores is given in Additional data file 13.) Additional data file 5: Binary matrix of genes classified as having fully conserved expression events, based on ranked microarray spot intensity, at five different thresholds (1/6, 1/5, 1/4, 1/3, 1/2). Additional data file 6: Cumulative distributions summarizing pairwise comparisons of conservation of gene expression using the Pearson correlation measure. The cumulative distributions show the proportion of all 3,074 genes with Pearson r (normalized intensities) below the value shown on the horizontal axis, for real orthologs (green) and randomly matched genes (blue). Additional data file 7: Feature matrix used to compare conservation of expression measures to other attributes of individual genes, with Spearman correlations and p-values. Additional data file 8: WMW p-values for categorical gene attributes, with ranks determined by relative conservation of gene expression by median Pearson correlation for each species against Tetraodon. Additional data file 9: Cumulative distribution of EEL scores for real and permuted orthology between human and pufferfish. Additional data file 10: Breakdown of the proportion of all genes in each species that are expressed within each tissue. Additional data file 11: List of genes classified as TFs on the basis of containing a known DNA binding domain. Additional data file 12: Clustergrams showing Spearman correlations and p-values for comparisons of gene expression conservation versus other gene properties. Additional data file 13: WMW enrichment p-values of genes associated with GO biological process annotations expressed within each tissue of each species (full matrix used to create Additional data file 4).
We thank Rick Winterbottom, Norm White, Lloyd Berger, Usha Bhargava and Michael Rennie for their assistance with the animal husbandry and dissections and Peter Roy, Martha Bulyk and Anthony Phillipakis for helpful conversations. This work was supported by grants from the Howard Hughes Medical Institute, the Canadian Foundation for Innovation, and the Ontario Research Fund to TRH, NSERC Discovery grants to QDM and MB, an NCIC grant to AW, and NIH grant R01-GM1080 to MB.
- International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.View Article
- International Human Genome Sequencing Consortium: Finishing the euchromatic sequence of the human genome. Nature. 2004, 431: 931-945. 10.1038/nature03001.View Article
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biémont C, Skalli Z, Cattolico L, Poulain J, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate protokaryotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.PubMedView Article
- Mouse Genome Sequencing Consortium, Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, Antonarakis SE, Attwood J, Baertsch R, Bailey J, Barlow K, Beck S, Berry E, Birren B, Bloom T, Bork P, Botcherby M, Bray N, Brent MR, Brown DG, Brown SD, Bult C, Burton J, Butler J, Campbell RD, et al: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.View Article
- King MC, Wilson AC: Evolution at two levels in humans and chimpanzees. Science. 1975, 188: 107-116. 10.1126/science.1090005.PubMedView Article
- Cooper TF, Rozen DE, Lenski RE: Parallel changes in gene expression after 20,000 generations of evolution in Escherichia coli. Proc Natl Acad Sci USA. 2003, 100: 1072-1077. 10.1073/pnas.0334340100.PubMedPubMed CentralView Article
- Ferea TL, Botstein D, Brown PO, Rosenzweig RF: Systematic changes in gene expression patterns following adaptive evolution in yeast. Proc Natl Acad Sci USA. 1999, 96: 9721-9726. 10.1073/pnas.96.17.9721.PubMedPubMed CentralView Article
- Gompel N, Prud'homme B, Wittkopp PJ, Kassner VA, Carroll SB: Chance caught on the wing: cis-regulatory evolution and the origin of pigment patterns in Drosophila. Nature. 2005, 433: 481-487. 10.1038/nature03235.PubMedView Article
- Miller CT, Beleza S, Pollen AA, Schluter D, Kittles RA, Shriver MD, Kingsley DM: cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans. Cell. 2007, 131: 1179-1189. 10.1016/j.cell.2007.10.055.PubMedPubMed CentralView Article
- Prud'homme B, Gompel N, Rokas A, Kassner VA, Williams TM, Yeh SD, True JR, Carroll SB: Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene. Nature. 2006, 440: 1050-1053. 10.1038/nature04597.PubMedView Article
- Shapiro MD, Marks ME, Peichel CL, Blackman BK, Nereng KS, Jonsson B, Schluter D, Kingsley DM: Genetic and developmental basis of evolutionary pelvic reduction in threespine sticklebacks. Nature. 2004, 428: 717-723. 10.1038/nature02415.PubMedView Article
- Wittkopp PJ, Haerum BK, Clark AG: Evolutionary changes in cis and trans gene regulation. Nature. 2004, 430: 85-88. 10.1038/nature02698.PubMedView Article
- Wray GA: Transcriptional regulation and the evolution of development. Int J Dev Biol. 2003, 47: 675-684.PubMed
- Hammock EA, Young LJ: Microsatellite instability generates diversity in brain and sociobehavioral traits. Science. 2005, 308: 1630-1634. 10.1126/science.1111427.PubMedView Article
- Lim MM, Wang Z, Olazabal DE, Ren X, Terwilliger EF, Young LJ: Enhanced partner preference in a promiscuous species by manipulating the expression of a single gene. Nature. 2004, 429: 754-757. 10.1038/nature02539.PubMedView Article
- Anand S, Wang WC, Powell DR, Bolanowski SA, Zhang J, Ledje C, Pawashe AB, Amemiya CT, Shashikant CS: Divergence of Hoxc8 early enhancer parallels diverged axial morphologies between mammals and fishes. Proc Natl Acad Sci USA. 2003, 100: 15666-15669. 10.1073/pnas.2535667100.PubMedPubMed CentralView Article
- Khaitovich P, Weiss G, Lachmann M, Hellmann I, Enard W, Muetzel B, Wirkner U, Ansorge W, Paabo S: A neutral model of transcriptome evolution. PLoS Biol. 2004, 2: E132-10.1371/journal.pbio.0020132.PubMedPubMed CentralView Article
- Zhang Y, Sturgill D, Parisi M, Kumar S, Oliver B: Constraint and turnover in sex-biased gene expression in the genus Drosophila. Nature. 2007, 450: 233-237. 10.1038/nature06323.PubMedPubMed CentralView Article
- Yanai I, Graur D, Ophir R: Incongruent expression profiles between human and mouse orthologous genes suggest widespread neutral evolution of transcription control. OMICS. 2004, 8: 15-24. 10.1089/153623104773547462.PubMedView Article
- ENCODE Project Consortium, Birney E, Stamatoyannopoulos JA, Dutta A, Guigó R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, Thurman RE, Kuehn MS, Taylor CM, Neph S, Koch CM, Asthana S, Malhotra A, Adzhubei I, Greenbaum JA, Andrews RM, Flicek P, Boyle PJ, Cao H, Carter NP, Clelland GK, Davis S, Day N, Dhami P, Dillon SC, Dorschner MO, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447: 799-816. 10.1038/nature05874.View Article
- Yan W, Ma L, Zilinski CA, Matzuk MM: Identification and characterization of evolutionarily conserved pufferfish, zebrafish, and frog orthologs of GASZ. Biol Reprod. 2004, 70: 1619-1625. 10.1095/biolreprod.103.024778.PubMedView Article
- Chervenak AP, Basu P, Shin M, Redmond LC, Sheng G, Lloyd JA: Identification, characterization, and expression pattern of the chicken EKLF gene. Dev Dyn. 2006, 235: 1933-1940. 10.1002/dvdy.20829.PubMedView Article
- Shi X, Bosenko DV, Zinkevich NS, Foley S, Hyde DR, Semina EV, Vihtelic TS: Zebrafish pitx3 is necessary for normal lens and retinal development. Mech Dev. 2005, 122: 513-527. 10.1016/j.mod.2004.11.012.PubMedView Article
- Marza E, Barthe C, Andre M, Villeneuve L, Helou C, Babin PJ: Developmental expression and nutritional regulation of a zebrafish gene homologous to mammalian microsomal triglyceride transfer protein large subunit. Dev Dyn. 2005, 232: 506-518. 10.1002/dvdy.20251.PubMedView Article
- Wray GA: The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007, 8: 206-216. 10.1038/nrg2063.PubMedView Article
- Bejerano G, Pheasant M, Makunin I, Stephen S, Kent WJ, Mattick JS, Haussler D: Ultraconserved elements in the human genome. Science. 2004, 304: 1321-1325. 10.1126/science.1098119.PubMedView Article
- Pennacchio LA, Ahituv N, Moses AM, Prabhakar S, Nobrega MA, Shoukry M, Minovitsky S, Dubchak I, Holt A, Lewis KD, Plajzer-Frick I, Akiyama J, De Val S, Afzal V, Black BL, Couronne O, Eisen MB, Visel A, Rubin EM: In vivo enhancer analysis of human conserved non-coding sequences. Nature. 2006, 444: 499-502. 10.1038/nature05295.PubMedView Article
- Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, Weinstock GM, Wilson RK, Gibbs RA, Kent WJ, Miller W, Haussler D: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005, 15: 1034-1050. 10.1101/gr.3715005.PubMedPubMed CentralView Article
- Woolfe A, Goodson M, Goode DK, Snell P, McEwen GK, Vavouri T, Smith SF, North P, Callaway H, Kelly K, Walter K, Abnizova I, Gilks W, Edwards YJ, Cooke JE, Elgar G: Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 2005, 3: e7-10.1371/journal.pbio.0030007.PubMedPubMed CentralView Article
- Prabhakar S, Poulin F, Shoukry M, Afzal V, Rubin EM, Couronne O, Pennacchio LA: Close sequence comparisons are sufficient to identify human cis-regulatory elements. Genome Res. 2006, 16: 855-863. 10.1101/gr.4717506.PubMedPubMed CentralView Article
- Thomas JW, Touchman JW, Blakesley RW, Bouffard GG, Beckstrom-Sternberg SM, Margulies EH, Blanchette M, Siepel AC, Thomas PJ, McDowell JC, Maskeri B, Hansen NF, Schwartz MS, Weber RJ, Kent WJ, Karolchik D, Bruen TC, Bevan R, Cutler DJ, Schwartz S, Elnitski L, Idol JR, Prasad AB, Lee-Lin SQ, Maduro VV, Summers TJ, Portnoy ME, Dietrich NL, Akhter N, Ayele K, et al: Comparative analyses of multi-species sequences from targeted genomic regions. Nature. 2003, 424: 788-793. 10.1038/nature01858.PubMedView Article
- Dermitzakis ET, Clark AG: Evolution of transcription factor binding sites in mammalian gene regulatory regions: conservation and turnover. Mol Biol Evol. 2002, 19: 1114-1121.PubMedView Article
- Odom DT, Dowell RD, Jacobsen ES, Gordon W, Danford TW, MacIsaac KD, Rolfe PA, Conboy CM, Gifford DK, Fraenkel E: Tissue-specific transcriptional regulation has diverged significantly between human and mouse. Nat Genet. 2007, 39: 730-732. 10.1038/ng2047.PubMedPubMed CentralView Article
- Oda-Ishii I, Bertrand V, Matsuo I, Lemaire P, Saiga H: Making very similar embryos with divergent genomes: conservation of regulatory mechanisms of Otx between the ascidians Halocynthia roretzi and Ciona intestinalis. Development. 2005, 132: 1663-1674. 10.1242/dev.01707.PubMedView Article
- Sanges R, Kalmar E, Claudiani P, D'Amato M, Muller F, Stupka E: Shuffling of cis-regulatory elements is a pervasive feature of the vertebrate lineage. Genome Biol. 2006, 7: R56-10.1186/gb-2006-7-7-r56.PubMedPubMed CentralView Article
- Huntley S, Baggott DM, Hamilton AT, Tran-Gyamfi M, Yang S, Kim J, Gordon L, Branscomb E, Stubbs L: A comprehensive catalog of human KRAB-associated zinc finger genes: insights into 27 the evolutionary history of a large family of transcriptional repressors. Genome Res. 2006, 16: 669-677. 10.1101/gr.4842106.PubMedPubMed CentralView Article
- Blekhman R, Oshlack A, Chabot AE, Smyth GK, Gilad Y: Gene regulation in primates evolves under tissue-specific selection pressures. PLoS Genet. 2008, 4: e1000271-10.1371/journal.pgen.1000271.PubMedPubMed CentralView Article
- Gilad Y, Oshlack A, Smyth GK, Speed TP, White KP: Expression profiling in primates reveals a rapid evolution of human transcription factors. Nature. 2006, 440: 242-245. 10.1038/nature04559.PubMedView Article
- Wagner GP, Lynch VJ: The gene regulatory logic of transcription factor evolution. Trends Ecol Evol. 2008, 23: 377-385. 10.1016/j.tree.2008.03.006.PubMedView Article
- Rifkin SA, Kim J, White KP: Evolution of gene expression in the Drosophila melanogaster subgroup. Nat Genet. 2003, 33: 138-144. 10.1038/ng1086.PubMedView Article
- Venkatesh B, Yap WH: Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements. Bioessays. 2005, 27: 100-107. 10.1002/bies.20134.PubMedView Article
- Fisher S, Grice EA, Vinton RM, Bessling SL, McCallion AS: Conservation of RET regulatory function from human to zebrafish without sequence similarity. Science. 2006, 312: 276-279. 10.1126/science.1124070.PubMedView Article
- Schadt EE, Edwards SW, GuhaThakurta D, Holder D, Ying L, Svetnik V, Leonardson A, Hart KW, Russell A, Li G, Cavet G, Castle J, McDonagh P, Kan Z, Chen R, Kasarskis A, Margarint M, Caceres RM, Johnson JM, Armour CD, Garrett-Engele PW, Tsinoremas NF, Shoemaker DD: A comprehensive transcript index of the human genome generated using microarrays and computational approaches. Genome Biol. 2004, 5: R73-10.1186/gb-2004-5-10-r73.PubMedPubMed CentralView Article
- Zhang W, Morris QD, Chang R, Shai O, Bakowski MA, Mitsakakis N, Mohammad N, Robinson MD, Zirngibl R, Somogyi E, Laurin N, Eftekharpour E, Sat E, Grigull J, Pan Q, Peng WT, Krogan N, Greenblatt J, Fehlings M, Kooy van der D, Aubin J, Bruneau BG, Rossant J, Blencowe BJ, Frey BJ, Hughes TR: The functional landscape of mouse gene expression. J Biol. 2004, 3: 21-10.1186/jbiol16.PubMedPubMed CentralView Article
- Chung WY, Albert R, Albert I, Nekrutenko A, Makova KD: Rapid and asymmetric divergence of duplicate genes in the human gene coexpression network. BMC Bioinformatics. 2006, 7: 46-10.1186/1471-2105-7-46.PubMedPubMed CentralView Article
- Gu X, Zhang Z, Huang W: Rapid evolution of expression and regulatory divergences after yeast gene duplication. Proc Natl Acad Sci USA. 2005, 102: 707-712. 10.1073/pnas.0409186102.PubMedPubMed CentralView Article
- Remm M, Storm CE, Sonnhammer EL: Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001, 314: 1041-1052. 10.1006/jmbi.2000.5197.PubMedView Article
- Hedges SB: The origin and evolution of model organisms. Nat Rev Genet. 2002, 3: 838-849. 10.1038/nrg929.PubMedView Article
- Khaitovich P, Hellmann I, Enard W, Nowick K, Leinweber M, Franz H, Weiss G, Lachmann M, Paabo S: Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees. Science. 2005, 309: 1850-1854. 10.1126/science.1108296.PubMedView Article
- Khaitovich P, Enard W, Lachmann M, Paabo S: Evolution of primate gene expression. Nat Rev Genet. 2006, 7: 693-702. 10.1038/nrg1940.PubMedView Article
- Hallikas O, Palin K, Sinjushina N, Rautiainen R, Partanen J, Ukkonen E, Taipale J: Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity. Cell. 2006, 124: 47-59. 10.1016/j.cell.2005.10.042.PubMedView Article
- Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32: D91-D94. 10.1093/nar/gkh012.PubMedPubMed CentralView Article
- Whitehead A, Crawford DL: Variation within and among species in gene expression: raw material for evolution. Mol Ecol. 2006, 15: 1197-1211. 10.1111/j.1365-294X.2006.02868.x.PubMedView Article
- Stuart JM, Segal E, Koller D, Kim SK: A gene-coexpression network for global discovery of conserved genetic modules. Science. 2003, 302: 249-255. 10.1126/science.1087447.PubMedView Article
- Meiklejohn CD, Parsch J, Ranz JM, Hartl DL: Rapid evolution of male-biased gene expression in Drosophila. Proc Natl Acad Sci USA. 2003, 100: 9894-9899. 10.1073/pnas.1630690100.PubMedPubMed CentralView Article
- Zapata A, Diez B, Cejalvo T, Gutierrez-de Frias C, Cortes A: Ontogeny of the immune system of fish. Fish Shellfish Immunol. 2006, 20: 126-136. 10.1016/j.fsi.2004.09.005.PubMedView Article
- Wilson MD, Barbosa-Morais NL, Schmidt D, Conboy CM, Vanes L, Tybulewicz VL, Fisher EM, Tavare S, Odom DT: Species-specific transcription in mice carrying human chromosome 21. Science. 2008, 322: 434-438. 10.1126/science.1160930.PubMedPubMed CentralView Article
- Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, Tanenbaum DM, White TJ, Sninsky JJ, Hernandez RD, Civello D, Adams MD, Cargill M, Clark AG: Natural selection on protein-coding genes in the human genome. Nature. 2005, 437: 1153-1157. 10.1038/nature04240.PubMedView Article
- Lopez-Bigas N, De S, Teichmann SA: Functional protein divergence in the evolution of Homo sapiens. Genome Biol. 2008, 9: R33-10.1186/gb-2008-9-2-r33.PubMedPubMed CentralView Article
- Berger MF, Badis G, Gehrke AR, Talukder S, Philippakis AA, Pena-Castillo L, Alleyne TM, Mnaimneh S, Botvinnik OB, Chan ET, Khalid F, Zhang W, Newburger D, Jaeger SA, Morris QD, Bulyk ML, Hughes TR: Variation in homeodomain DNA binding revealed by high-resolution analysis of sequence preferences. Cell. 2008, 133: 1266-1276. 10.1016/j.cell.2008.05.024.PubMedPubMed CentralView Article
- Luscombe NM, Thornton JM: Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity. J Mol Biol. 2002, 320: 991-1009. 10.1016/S0022-2836(02)00571-5.PubMedView Article
- Segal E, Raveh-Sadka T, Schroeder M, Unnerstall U, Gaul U: Predicting expression patterns from regulatory sequence in Drosophila segmentation. Nature. 2008, 451: 535-540. 10.1038/nature06496.PubMedView Article
- Conservation of Core Gene Expression in Vertebrate Tissues: Supplementary Data. [http://hugheslab.ccbr.utoronto.ca/supplementary-data/vertebrate_expression]
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Herrero J, Holland R, Howe K, Howe K, Johnson N, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Melsopp C, Megy K, et al: Ensembl 2007. Nucleic Acids Res. 2007, 35: D610-D617. 10.1093/nar/gkl996.PubMedPubMed CentralView Article
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMedPubMed CentralView Article
- Kent WJ: BLAT – the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMedPubMed CentralView Article
- O'Brien KP, Remm M, Sonnhammer EL: Inparanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res. 2005, 33: D476-D480. 10.1093/nar/gki107.PubMedPubMed CentralView Article
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView Article
- Alexeyenko A, Tamas I, Liu G, Sonnhammer EL: Automatic clustering of orthologs and inparalogs shared by multiple proteomes. Bioinformatics. 2006, 22: e9-e15. 10.1093/bioinformatics/btl213.PubMedView Article
- Kasprzyk A, Keefe D, Smedley D, London D, Spooner W, Melsopp C, Hammond M, Rocca-Serra P, Cox T, Birney E: EnsMart: a generic system for fast and flexible access to biological data. Genome Res. 2004, 14: 160-169. 10.1101/gr.1645104.PubMedPubMed CentralView Article
- Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18 (Suppl 1): S96-S104.PubMedView Article
- Schug J, Schuller WP, Kappen C, Salbaum JM, Bucan M, Stoeckert CJ: Promoter features related to tissue specificity as measured by Shannon entropy. Genome Biol. 2005, 6: R33-10.1186/gb-2005-6-4-r33.PubMedPubMed CentralView Article
- UCSC Genome Bioinformatics. [http://genome.ucsc.edu]
- Electronic Supplement: Ultraconserved Elements in the Human Genome. [http://www.soe.ucsc.edu/~jill/ultra.html]
- Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2008, 36: D281-D288. 10.1093/nar/gkm960.PubMedPubMed CentralView Article
- Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A, Batzoglou S: LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003, 13: 721-731. 10.1101/gr.926603.PubMedPubMed CentralView Article
- Brudno M, Malde S, Poliakov A, Do CB, Couronne O, Dubchak I, Batzoglou S: Glocal alignment: finding rearrangements during alignment. Bioinformatics. 2003, 19 (Suppl 1): i54-i62. 10.1093/bioinformatics/btg1005.PubMedView Article
- Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, Hardison RC, Haussler D, Miller W: Human-mouse alignments with BLASTZ. Genome Res. 2003, 13: 103-107. 10.1101/gr.809403.PubMedPubMed CentralView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.