PseudoGeneQuest – Service for identification of different pseudogene types in the human genome
© Ortutay and Vihinen; licensee BioMed Central Ltd. 2008
Received: 18 February 2008
Accepted: 02 July 2008
Published: 02 July 2008
Pseudogenes, nonfunctional copies of genes, evolve fast due the lack of evolutionary pressures and thus appear in several different forms. PseudoGeneQuest is an online tool to search the human genome for a given query sequence and to identify different types of pseudogenes as well as novel genes and gene fragments.
The service can detect pseudogenes, that have arisen either by retrotransposition or segmental genome duplication, many of which are not listed in the public pseudogene databases. The service has a user-friendly web interface and uses a powerful computer cluster in order to perform parallel searches and provide relatively fast runtimes despite exhaustive database searches and analyses.
PseudoGeneQuest is a versatile tool for detecting novel pseudogene candidates from the human genome. The service searches human genome sequences for five types of pseudogenes and provides an output that allows easy further analysis of observations. In addition to the result file the system provides visualization of the results linked to Ensembl Genome Browser. PseudoGeneQuest service is freely available.
Pseudogenes are nonfunctional copies of genes whose transcription or translation is disrupted. They have several biological roles; for example, they may regulate gene expression and create a reservoir of diversity at a genetic or phenotypic level . The recombination between the immunoglobulin Vh gene and its pseudogenes has been proposed to drive the formation of new genes . Pseudogenes can also drive gene conversion to contribute to immunoglobulin heavy and light chain diversity .
Pseudogenes arise either by retrotransposition or duplication of genomic DNA segments . Because of these fundamentally different mechanisms, the different types of pseudogenes have characteristic structures which facilitate identification and classification.
Processed pseudogenes, which emerged via retrotransposons , lack introns, a poly A tract at the 3' end, and flanking repeats. Other pseudogenes are complete or partial duplicates of real genes with an interruption in their transcription or translation. The identification of pseudogenes from genomes can be a difficult task. Due to the missing evolutionary pressure, pseudogene sequences evolve neutrally therefore mutations accumulate at a fast pace, i.e. they diverge fast from the original sequences. Homologous sequence stretches missing functional characteristics (such as promoters, interrupted coding regions etc.) are considered pseudogenes.
Two major bioinformatic services have been developed and distributed for the analysis and annotation of pseudogenes. Pseudogene.org  is an extensive database and analysis system for pseudogenes from several genomes, including human. The database facilitates the comparison of pseudogenes in selected genomes. The data is collected with an automatic PseudoPipe system . Another service, HOPPSIGEN , is dedicated for processed pseudogenes in human and mouse. The database seems to have not been updated recently. Both these services are useful, but have some limitations. If one wants to study recently identified genes then these services are not helpful. The PseudoPipe can be downloaded, but the user should be familiar with programming. As a solution for researchers looking for all types of pseudogenes, as well as for newly identified genes, we developed an online tool that implements the methods for pseudogene identification used in our previous study for tracking human immune system related pseudogenes . At the moment PseudoGeneQuest is the only publicly available web based system, which identifies all types of pseudogenes and pseudogene fragments in the human genome using user provided query sequence as starting point. It is a unique tool for searching and identifying new human pseudogenes.
Hits for the same query that are within 60 nucleotides (nt) of each other are merged as single hits and extended by 30 nt at both ends. These candidates are aligned back to the original query sequence using fastx , and the hits with less than 40% amino acid identity are discarded. The pseudogene candidate and the query protein sequences are then aligned using tblastn. If the e-value of the alignment is greater than 10-10, the hit is discarded. If the identity is >95%, the length is >95% of the query, and the reading frame does not contain an internal stop codon, the hit is considered a new gene candidate. False positives have an e-value >10-10 or an identity <40% at the amino acid level in the fasta alignment.
Multiexon pseudogene candidates are grouped based on their repeat content as determined by Repeat Masker [18, 19]. If the repeat content is <50% of the length of the target, the sequence is identified as a duplicated multiexon pseudogene; otherwise it is classified as an interrupted processed pseudogene.
For pseudogene fragments, the alignment of a single-exon pseudogene candidate covers <70% of the query sequence. If the candidate covers >70% of the known gene, the amino acid identity is >40%, and it contains a stop codon(s) or has a frame shift(s), it is classified as a true pseudogene. In summary, the identified pseudogenes are grouped into five categories: known and newly identified pseudogenes, pseudogene fragments, duplicated pseudogenes, and interrupted processed pseudogenes. In addition, putative new genes and exons can be identified.
The PseudoGeneQuest (PGQ) service runs on a Linux cluster currently consisting of 70 processors. Run times are relatively short; however, since the PGQ is based on an excessive and genome wide analysis, results are not immediately available. Users access the system via a web page  and provide either a protein sequence in FASTA format or a protein sequence ID. Each submission has a unique search id, which is used to retrieve results from the web page when ready. The server forwards the job to the cluster, where the actual analysis is performed. The analysis generally takes 5 to 20 minutes. To limit the length of the process we request the user on the query page to specify the maximum number of the hits to be analysed during the process. PGQ offers up 500 hits to analyse, which covers the needs of most users. The user is notified by email about the results.
The results are included in the email, but a more user-friendly version with links is available online. The result file contains hits in each pseudogene category and information for any new genes or gene fragments. The location of each hit in the reference genome sequence is provided. The used human genome build and Pseudogene.org database version is also recorded at the end of the result file together with the version of the PGQ software. The human genome database is updated automatically in our system upon release of new builds.
Genome segments identified by PseudoGeneQuest
Tag in the result file
ALREADY KNOWN GENE
Hits overlapping genes annotated in the genome files.
Results overlapping with records in pseudogene.org.
REAL GENE OR EXON
The hit matches almost exactly to the query. Frequently parts of yet un-annotated or predicted genes.
Covers <70% of the length of the query. Small parts of pseudogenes.
Covers >70% of the length of the query and the reading frame is broken. Processed pseudogenes with high homology to the query sequence.
Covers >70% of the length of the query and the reading frame is intact.
PUTATIVE NEW GENE
The hit has uninterrupted reading frame. Un-annotated gene.
Multiexon hit with <50% repeat content. A recent duplication of a gene.
INTERRUPTED PROCESSED PSEUDOGENE
Multiexon hit with >50% repeat content. Old processed pseudogenes which accumulated repeats.
'Duplicated' and 'interrupted processed pseudogenes' are multiexon candidates with different repeat content. They might be duplications of real genes or other interrupted pseudogenes. To these sequences different amount of repeats has accumulated after duplication events. Thus, these findings have diverged in a different extent from their ancestor sequences. 'Duplicated pseudogenes' are copies of genomic segment s containing genes with broken gene expression. The duplication happened recently therefore there was not enough time to accumulate repeat sequences. On the other hand, if the duplication took place earlier, the original sequence can be masked by repeat sequences. If more than half of the sequence consist of repeat sequence, then we distinguish it as "interrupted processed pseudogene".
To allow users to link the PGQ results to genomic context, there is a link to connect each of the identified and categorized segments to the genomic contigs in the Entrez database. Further the genome segments be visualized using Ensembl Genome Browser ContigView . This way users can obtain additional information about the identified regions their location in chromosomes.
Results and Discussion
Originally PGQ was applied for identification of pseudogenes related to 845 human immune system related proteins . The analysis identified a total of 4816 pseudogenes related to 313 genes, along with putative new genes and pseudogene fragments. Most of the new pseudogenes were pseudogene fragments: altogether 3736 fragments for 229 genes. The expression of the pseudogenes varies substantially based on the EST data. Many of these findings were novel and not listed in the Pseudogene.org database. All the results of the earlier extensive analysis are available in the Immunome database [22, 23].
Results of test analysis.
Already known gene
Real gene or exon
Putative new gene
Interrupted processed pseudogene
For researchers interested in long list of putative pseudogenes or pseudogene fragments related to a human gene or protein not listed in Pseudogene.org, there is another tool, PseudoPipe  available from the authors of the pseudogene.org database. The system is powerful but possibly out of reach for some users if they do not have the necessary skills to install and run the program package. For them and all others, PGQ provides an alternative, which is easy to use and allows extensive analysis of the query sequence. Its results can be compared to the Pseudogene.org database, however can provide additional information because the query sequece provided by the user can be any protein sequence, not just those used for the construction of the Pseudogene.org database. The system yields results, which the users can study further with other experimental or in silico methods.
PGQ is an online tool for users who want to identify pseudogenes related to a query sequence. PGQ groups results in altogether 9 categories (Table 1). The system has a user friendly web interface for placing the queries, uses a powerful computer cluster to handle parallel searches and notifies the user via email about finished analysis. Results are provided to the email and additionally via the PGQ web page, which provides also links to the identified sequences.
Availability and requirements
PseudoGeneQuest service is available freely for the scientific community via http://bioinf.uta.fi/PseudoGeneQuest. Non-academic usage is permitted upon agreement with the authors. The users need in addition to the query sequence only a recent recent version of a web browser.
We thank the Medical Research Fund of Tampere University Hospital for financial support.
- Balakirev ES, Ayala FJ: Pseudogenes: are they "junk" or functional DNA? Annu Rev Genet 2003, 37: 123–151. 10.1146/annurev.genet.37.040103.103949View ArticlePubMed
- Vargas-Madrazo E, Almagro C, Lara-Ochoa F: Structural repertoire in VHpseudogenes of immunoglobulins: comparison with human germline genes and human amino acid sequences. J Mol Biol 1995, 246: 74–81. 10.1006/jmbi.1994.0067View ArticlePubMed
- Arakawa H, Buerstedde JM: Immunoglobulin gene conversion: insights from bursal B cells and the DT40 cell line. Dev Dyn 2004, 229: 458–464. 10.1002/dvdy.10495View ArticlePubMed
- Mighell AJ, Smith R, Robinson PA, Markham AF: Vertebrate seudogenes. FEBS Lett 2000, 468: 109–114. 10.1016/S0014-5793(00)01199-6View ArticlePubMed
- Vanin EF: Processed pseudogenes: characteristics and evolution. Annu Rev Genet 1985, 19: 253–272. 10.1146/annurev.ge.19.120185.001345View ArticlePubMed
- Karro JE, Yan Y, Zheng D, Zhang Z, Carriero N, Cayting P, Harrrison P, Gerstein M: Pseudogene.org: a comprehensive database and comparison platform for pseudogene annotation. Nucleic Acids Res 2007, 35: D55–60. 10.1093/nar/gkl851PubMed CentralView ArticlePubMed
- Zhang Z, Carriero N, Zheng D, Karro J, Harrison PM, Gerstein M: PseudoPipe: an automated pseudogene identification pipeline. Bioinformatics 2006, 22: 1437–1439. 10.1093/bioinformatics/btl116View ArticlePubMed
- Khelifi A, Duret L, Mouchiroud D: HOPPSIGEN: a database of human and mouse processed pseudogenes. Nucleic Acids Res 2005, 33: D59–66. 10.1093/nar/gni059View ArticlePubMed
- Ortutay C, Siermala M, Vihinen M: Molecular characterization of the immune system: emergence of proteins, processes, and domains. Immunogenetics 2007, 59: 333–348. 10.1007/s00251-007-0191-0View ArticlePubMed
- Harrison PM, Hegyi H, Balasubramanian S, Luscombe NM, Bertone P, Echols N, Johnson T, Gerstein M: Molecular fossils in the human genome: identification and analysis of the pseudogenes in chromosomes 21 and 22. Genome Res 2002, 12: 272–280. 10.1101/gr.207102PubMed CentralView ArticlePubMed
- Zhang Z, Harrison P, Gerstein M: Identification and analysis of over 2000 ribosomal protein pseudogenes in the human genome. Genome Res 2002, 12: 1466–1482. 10.1101/gr.331902PubMed CentralView ArticlePubMed
- Zhang Z, Harrison PM, Liu Y, Gerstein M: Millions of years of evolution preserved: a comprehensive catalog of the processed pseudogenes in the human genome. Genome Res 2003, 13: 2541–2558. 10.1101/gr.1429003PubMed CentralView ArticlePubMed
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, Lehvaslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: Perl modules for the life sciences. Genome Res 2002, 12: 1611–1618. 10.1101/gr.361602PubMed CentralView ArticlePubMed
- Barak A, Guday S, Wheeler R: The MOSIX Distributed Operating System: Load Balancing for UNIX. New York: Springer-Verlag New York, Inc. Secaucus, NJ, USA; 1993.View Article
- McGinnis S, Madden TL: BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res 2004, 32: W20–25. 10.1093/nar/gkh435PubMed CentralView ArticlePubMed
- Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci USA 1988, 85: 2444–2448. 10.1073/pnas.85.8.2444PubMed CentralView ArticlePubMed
- Smit A, Hubley R, Green P: RepeatMasker Open-3.0. 1996–2004.
- Flicek P, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Eyre T, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Holland R, Howe KL, Howe K, Johnson N, Jenkinson A, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Megy K: Ensembl 2008. Nucleic Acids Res 2008, 36: D707–714. 10.1093/nar/gkm988PubMed CentralView ArticlePubMed
- Ortutay C, Vihinen M: Immunome: a reference set of genes and proteins for systems biology of the human immune system. Cell Immunol 2006, 244(2):87–89. 10.1016/j.cellimm.2007.01.012View ArticlePubMed
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.