[X]uniqMAP: unique gene sequence regions in the human and mouse genomes
© Jiménez and Durbin; licensee BioMed Central Ltd. 2006
Received: 12 July 2006
Accepted: 06 October 2006
Published: 06 October 2006
Current approaches for genome-wise functional analyses, such as microarray and RNA interference studies, rely on the specificity of oligonucleotide sequences to selectively target cellular transcripts. The design of specific oligos involves the determination of unique DNA regions in the gene/transcripts of interest from the targeted organism. This process is tedious, time consuming and it does not scale up for high-throughput studies.
Taking advantage of the availability of complete genome sequence information for mouse and human, the most widely used systems for the study of mammalian genetics, we have built a database, [X]uniqMAP, that stores the precalculated unique regions for all transcripts of these two organisms. For each gene, the database discriminates between those unique regions that are shared by all transcripts and those exclusive to single transcripts. In addition, it also provides those unique regions that are shared between orthologous genes from the two organisms. The database is updated regularly to reflect changes in genome assemblies and gene builds.
Over 85% of genes have unique regions at least 19 bases long, with the majority being unique over 60% of their lengths. 14482 human genes share exactly at least a unique region with mouse genes, though such regions are typically under 40 bases long. The full data are publicly accessible online both interactively and for download. They should facilitate (i) the design of probes, primers and siRNAs for both small- and large-scale projects; and (ii) the identification of regions for the design of oligos that could be re-used to target equivalent gene/transcripts from human and mouse.
Following the completion of several whole genome sequencing projects a considerable effort has been focused on genome-wise functional analyses of a number of organisms (reviewed in ). Some of the most popular methods are the study of gene expression by microarrays and phenotypic analyses from gene knock-downs by means of RNA interference techniques [2, 3]. The success of these methods relies in the ability of reagent oligonucleotides to specifically recognise single species of transcripts within the complex mixture present in the studied cells. Therefore, when designing probes, primers and siRNAs, the sequence specificity of candidate oligonucleotides must be assessed in order to minimise potential cross-hybridisations and off-target effects [4, 5]. Although cross-reaction events have been described between siRNAs and molecules of limited sequence similarity [6, 7], the determination of specificity routinely requires the identification of oligonucleotides that are identical in sequence only to the intended target. This uniqueness assessment is usually calculated every time that a new reagent needs to be designed. However, given the availability of complete genome sequences for a number of organisms, all their unique regions could be calculated, stored and made publicly available, for example, via an online resource. In addition, this resource could also take advantage of the known contextual relationships between transcripts within a gene to categorise uniqueness at the gene and transcript levels to, for example, easily discriminate between unique regions shared by all transcripts and those exclusive to individual alternative splicing variants. This information would simplify the process of oligo design by abolishing the step to determine unique fragments, with the desired range of action, within the gene/transcript of interest. At present, no such resource is available.
The [X]uniqMAP database has been developed to store and present already pre-calculated unique regions for all EnsEMBL transcripts of the human and mouse genomes , the most widely used systems for the study of mammalian genetics. It also records those unique fragments that are shared between them, which could help to identify sequences to simultaneously target equivalent genes between these two organisms. [X]uniqMAP differs from standard genomic browsers in that it uses genome comparisons to reveal unique regions within and between organisms. These regions are presented at both gene and transcript levels. The data stored in [X]uniqMAP can be retrieved via a user-friendly web interface or as downloadable FASTA files, and it should be useful for small- as well as large-scale projects for which the identification of unique DNA regions is required.
Construction and content
[X]uniqMAP is a collection of three databases: human and mouse uniqMAP as well as XuniqMAP. The first two collect the unique DNA 19-mers for all gene/transcripts within the two genomes. XuniqMAP collects those unique fragments within human and mouse that are also shared between them. The determination of the unique regions within a genome comprises three steps: (i) the construction of a set of non-redundant (NR) sequences, keeping track of the gene structure, from all the transcripts of each gene; (ii) self-comparison of the NR sequence set and storage of those positions found to be redundant; and (iii) deduction of the unique regions from regions identified in (ii) and gene structure data from (i).
The next step was to compare the sequences of the NR set with each other in order to determine which 19-mers were shared between several genes or transcripts. Prior to this, repeats and low complexity regions were filtered out  and the sequences were fragmented into all possible overlapping 19-mers excluding those containing filtered positions. Afterwards, the 19-mers were compared against the original, unfiltered, NR set using the exonerate program , recording only those matches, in either sense, at least 18 bases long and sharing 18 identities. At the end, the genomic coordinates of the central positions of the 19-mers containing filtered positions or matching to others were stored.
Finally, the unique regions were calculated from the above values kept in the database as follows. For regions unique to a gene, the segments shared by all transcripts of that gene were located onto the exonic regions from which the filtered and redundant positions were subsequently subtracted, yielding the remaining positions as unique. An identical procedure was applied to uncover the unique regions in individual transcripts with the exception that the positions initially mapped onto the transcripts were those exclusive to them. This discrimination of gene/transcript levels in unique regions is important for the design of oligos with the desired range of action. For human, unique regions were identified in 26704 out of 31274 genes, and 23438 transcripts from a total of 36161 transcripts from genes with annotated alternative splicing variants. For mouse, 26135 (out of 27735) genes and 11350 (out of 11807) transcripts contained unique regions.
For the calculation of unique positions shared between genomes, the two sets of unique fragments within the human and mouse genomes were compared against each other, storing information only for pairs of genes that shared at least five identical 19-mers in either sense. This arbitrary cut-off resulted almost completely in one-to-one, typically orthologous, relationships between human and mouse genes. The final set contains 15104 pairs made from 14482 and 14645 genes from human and mouse, respectively. As above, the gene/transcript relationships from each organism were transferred to the shared sequences which were then labelled according to their potential range of action across both organisms (see also next section for further details).
Utility and discussion
[X]uniqMAP contains a collection of the unique regions of all gene/transcripts in human and mouse as well as those unique sequences shared between their orthologs. For high throughput studies the complete data sets can be downloaded as FASTA files from the web site for further in-house calculations. The description lines of the sequences in the FASTA files contain all the appropriate information about the genomic co-ordinates and gene/transcript coverage of the fragments (see below for a more detailed explanation). On the other hand, in the more general case of experimentalists interested in only one or just few genes, this information can be retrieved interactively from the web site by searching the database either with an identifier or a nucleotide sequence.
When a sequence is provided, this will be BLASTed  against the set of unique sequences stored in the database. Only nucleotide sequences are allowed as queries. The search parameters have been adjusted to effectively deal with queries of any length, reporting only those local matches that are highly similar to the query, at least 19 nuleotides long and contain no gaps. The results page of the sequence search comprises graphical and sequence representations of the regions of the oligonucleotide query matching to those in the database highlighting the type of unique sequences, i.e. gene versus transcript, matched as well as the presence of mismatches if any.
where the value of symbol can be either Ug or Ut, depending on whether the fragment is shared by all transcripts in a gene or unique to a single transcript, respectively. For the former case, gene_id and transcript_id will be the same and it will correspond to the EnsEMBL gene identifier.
In XuniqMAP, the display is also split into two levels, namely reference and target, corresponding to all the unique regions of the gene in the reference organism and those matched by the target organism, respectively (Figure 3b). The colour-coded scheme for the reference is the same as in uniqMAP, i.e. red for unique regions shared by all transcripts and blue for those specific to individual transcripts (Figure 3b, top). For the target, the shared unique positions are placed relative to those matched with the reference and the colours represent all the possible combination of matches that can be found between the shared sequences: (i) red when they are present in all transcripts for both genes; (ii) blue if present only in individual transcripts for both genes; and (iii) green when present in all the transcripts of one gene but only in a single transcript of the other (Figure 3b, bottom). In the sequence display, only the gene of the target sequence is shown where the central positions are coloured as in the graphical display. Links are also provided for the retrieval of the sequences in FASTA format, for each type of match. In this case, the description lines for the fragment contain the following information:
where, as in uniqMAP, the symbol nomenclature will describe whether the fragment is shared or not by all transcripts.
[X]uniqMAP is a collection of databases that stores all unique regions within human and mouse gene/transcripts as well as the sub-set of these unique regions that are identical between these two organisms. All the information is provided in the adequate genomic context for each gene, reflecting the relationships between their transcripts. This discrimination between regions shared by all transcripts of a gene from those exclusive to individual transcripts should be useful to determine the range of action desired when designing oligos. For oligo design, [X]uniqMAP information could be fed, if necessary, directly to downstream pipelines (i) to assess further the adequacy of some segments as, for example, RNAi or primers; or (ii) for extra filtering if lower similarity to other potential cross-reactive molecules as well as the precise location of the mismatched positions must be accounted for. [X]uniqMAP is based on EnsEMBL and is updated in parallel with EnsEMBL to reflect changes in genome assemblies or gene builds.
Availability and requirements
[X]uniqMAP is available online from http://www.mitocheck.org/cgi-bin/uniqMAP/Xsearch. Enquiries on technical or data access issues can be mailed to firstname.lastname@example.org.
We are grateful to the Sixth Framework European Programme, via MitoCheck, (JLJ) and the Wellcome Trust (RD) for funding.
- Lee NH: Genomic approaches for reconstructing gene networks. Pharmacogenomics. 2005, 6 (3): 245-258. 10.1517/146224220.127.116.11.PubMedView ArticleGoogle Scholar
- Eisen MB, Brown PO: DNA arrays for analysis of gene expression. Methods Enzymol. 1999, 303: 179-205.PubMedView ArticleGoogle Scholar
- Vanhecke D, Janitz M: High-throughput gene silencing using cell arrays. Oncogene. 2004, 23 (51): 8353-8358. 10.1038/sj.onc.1208027.PubMedView ArticleGoogle Scholar
- Naito Y, Yamada T, Ui-Tei K, Morishita S, Saigo K: siDirect: highly effective, target-specific siRNA design software for mammalian RNA interference. Nucleic Acids Res. 2004, W124-129. 32 Web Server
- Nielsen HB, Wernersson R, Knudsen S: Design of oligonucleotides for microarrays and perspectives for design of multi-transcriptome arrays. Nucleic Acids Res. 2003, 31 (13): 3491-3496. 10.1093/nar/gkg622.PubMedPubMed CentralView ArticleGoogle Scholar
- Jackson AL, Bartz SR, Schelter J, Kobayashi SV, Burchard J, Mao M, Li B, Cavet G, Linsley PS: Expression profiling reveals off-target gene regulation by RNAi. Nat Biotechnol. 2003, 21 (6): 635-637. 10.1038/nbt831.PubMedView ArticleGoogle Scholar
- Martinez J, Tuschl T: RISC is a 5' phosphomonoester-producing RNA endonuclease. Genes Dev. 2004, 18 (9): 975-980. 10.1101/gad.1187904.PubMedPubMed CentralView ArticleGoogle Scholar
- Birney E, Andrews D, Caccamo M, Chen Y, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T: Ensembl 2006. Nucleic Acids Res. 2006, D556-561. 10.1093/nar/gkj133. 34 Database
- Smith AFA, Green P: RepeatMasker. 1996, [http://repeatmasker.genome.washington.edu]Google Scholar
- Slater GS, Birney E: Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005, 6: 31-10.1186/1471-2105-6-31.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1006/jmbi.1990.9999.PubMedView ArticleGoogle Scholar
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.