Open Access

t2prhd: a tool to study the patterns of repeat evolution

  • Botond Sipos1Email author,
  • Kálmán Somogyi1,
  • István Andó1 and
  • Zsolt Pénzes1, 2
BMC Bioinformatics20089:27

DOI: 10.1186/1471-2105-9-27

Received: 22 August 2007

Accepted: 18 January 2008

Published: 18 January 2008

Abstract

Background

The models developed to characterize the evolution of multigene families (such as the birth-and-death and the concerted models) have also been applied on the level of sequence repeats inside a gene/protein. Phylogenetic reconstruction is the method of choice to study the evolution of gene families and also sequence repeats in the light of these models. The characterization of the gene family evolution in view of the evolutionary models is done by the evaluation of the clustering of the sequences with the originating loci in mind. As the locus represents positional information, it is straightforward that in the case of the repeats the exact position in the sequence should be used, as the simple numbering according to repeat order can be misleading.

Results

We have developed a novel rapid visual approach to study repeat evolution, that takes into account the exact repeat position in a sequence. The "pairwise repeat homology diagram" visualizes sequence repeats detected by a profile HMM in a pair of sequences and highlights their homology relations inferred by a phylogenetic tree. The method is implemented in a Perl script (t2prhd) available for downloading at http://t2prhd.sourceforge.net and is also accessible as an online tool at http://t2prhd.brc.hu. The power of the method is demonstrated on the EGF-like and fibronectin-III-like (Fn-III) domain repeats of three selected mammalian Tenascin sequences.

Conclusion

Although pairwise repeat homology diagrams do not carry all the information provided by the phylogenetic tree, they allow a rapid and intuitive assessment of repeat evolution. We believe, that t2prhd is a helpful tool with which to study the pattern of repeat evolution. This method can be particularly useful in cases of large datasets (such as large gene families), as the command line interface makes it possible to automate the generation of pairwise repeat homology diagrams with the aid of scripts.

Background

The conceptual models developed to explain the evolution of multigene families [1, 2] have also been applied to sequence repeats inside a gene/protein. Studies of repeat evolution [39] have often revealed complex patterns: some repeats may evolve in concert, subject to homogenization, while other repeats may maintain their sequence identity, which is more consistent with birth-and-death and/or divergent evolution. Phylogenetic reconstruction is an established approach with which to study the mode of evolution of multigene families [1014] and repeats [39]. In the analysis of closely related species, the association of genes on the same locus (orthologues) on phylogenetic trees suggests independent evolution of the respective members, as expected for the birth-and-death model. On the contrary, the association of paralogues from the same species suggests homogenization, which is consistent with the concerted model. The application of the same logic to repeats poses a number of problems. Studies on the evolution of multigene families by phylogenetic methods make use of a prior concept of homology defined by the locus. As the term "locus" denotes genomic position, it follows that the exact repeat position in the full sequence (starting and ending positions) should be used to identify them. This approach has seldom been used in the previous studies, probably because it complicates the interpretation of the phylogenetic trees. Instead, a simple numbering according to repeat order is often applied, although this can be misleading: if a repeat is not detected by the selected method (for example the profile Hidden Markov Model (HMM) method) in one of the compared sequences, the ordinal numbers of the following repeats will be shifted by one and hence the analysis will become laborious.

Results and Discussion

To simplify the analyses of repeat evolution, a novel rapid visual method has been developed to highlight the relationships of repeats detected in a pair of sequences by the profile HMM method. The primary aim of the application is the analysis of tandem amino acid sequence repeats but DNA sequences can also be used with a profile HMM trained with DNA sequences. When DNA sequences are used, the script does not consider the reverse complement of the sequences, so the usefulness of this application in a genomic context (for example in the case of transportable elements) might be limited, though the implementation can be easily extended to handle this kind of analysis.

The "pairwise repeat homology diagram" visualizes repeats in a simple scheme, together with their homology relations (orthology and paralogy) inferred by a phylogenetic tree, providing an intuitive way to analyse the patterns of repeat evolution. This visualization also facilitates a survey of the sequence structure (size, linker sequences and non-repetitive regions). Regions in which consecutive repeats are connected with those in the other sequence forming a ladder-like pattern are most likely to have evolved independently following the birth-and-death and/or divergent process. On the contrary, in regions where repeats have only internal or no identified connections, repeats possibly evolve in concert. The patterns demonstrated by internal connections can reveal units of concerted evolution or recent internal duplications. It should be taken into account, however, that if there are clades formed by more than two identical sequences, the branching pattern of such clades, and hence the identified relation, is arbitrary and uninformative in this kind of analysis.

Implementation

The method is implemented in a Perl script with command line interface (t2prhd, standing for "tree to pairwise repeat homology diagram"). The repeats are first identified with hmmsearch from the HMMER package [15] and a profile HMM specified by the user. Raw search results are parsed by using BioPerl[16] modules and the extracted repeats are aligned with hmmalign. The resulting alignment is converted into Fasta and sequencial PHYLIP formats and a phylogenetic tree is built by using CLUSTAL W[17] (with default parameters) or by using PhyML[18]. The script reads in the resulting tree as a Bio::Tree::TreeI object. After getting the list of all leaf nodes it finds the "sister leaf nodes" (leaf node pairs having the same ancestor) by an algorithm that for n leaf nodes (repeats) needs in the best case n 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGUbGBaeaacqaIYaGmaaaaaa@2EC8@ and in the worst case n iterations. So, the asymptotic upper bound for the time complexity of this algorithm is O(n) and the asymptotic lower bound is Ω( n 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGUbGBaeaacqaIYaGmaaaaaa@2EC8@ ). These sister leaf nodes have a most recent common ancestor as indicated by the tree and accordingly we regard them as unambiguously identified homologues.

The script creates an SVG (Scalable Vector Graphics) file containing the diagram by using the XML::Writer module and saves the outputs and logs produced by the external applications. Optionally, it can also generate a LaTeX output. The identified homology relations are represented by lines or arcs connecting the respective repeats: blue lines are drawn between repeats from different sequences (orthology), and brown arcs in cases of internal relations (paralogy). The colour intensity of the connecting lines is a function of the patristic distance between the respective leaf nodes (d) divided by the total tree length (T): ( 1 d T ) w MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaWaaeWaaeaacqaIXaqmcqGHsisljuaGdaWcaaqaaiabdsgaKbqaaiabdsfaubaaaOGaayjkaiaawMcaamaaCaaaleqabaGaem4DaChaaaaa@3407@ . The default value of the "colour gradient parameter" w is 1 (linear colour scale), but by setting this parameter the colour scale can be tuned so that one can discriminate between close distance values. To make the interpretation of the colour scale easier a legend is drawn.

The generated SVG file can be viewed by Firefox 1.5 [19] or higher as an example and can be rasterized (and also viewed) using the Batik toolkit [20] or any other image editor capable to handle the SVG format. The generated LaTeX file should be processed by pdflatex, the pgf/TikZ, fancyhdr, xcolor and fullpage packages are required. The script has a manual page embedded in POD format.

The script is also accessible as an online tool [21] (with CLUSTAL W back-end only) through a web interface created by using the Pise form generator [22]. The real power of this visual method is manifested in studies with large data sets, where the analysis of numerous or large trees would be highly laborious. In these cases, it is very advantageous that the command line interface enables the use of the scripts to automate the diagram generation.

Example

Tenascins are extracellular matrix glycoproteins containing regions of repeated EGF-like and fibronectin-III-like (Fn-III) domains. The evolution of these repeated domains in mammalian Tenascins has been studied in detail by means of phylogenetic and other methods [8]. To illustrate the power of the method, we generated pairwise repeat homology diagrams with three selected protein sequences. The protein sequences corresponding to the DNA sequences studied by Hughes [8] were used (abbreviations and GenBank accession numbers in parentheses): human Tenascin X (TXH, [GenBank: AAB47488.1]), human Tenascin C (TCH, [GenBank: CAA39628.1]) and murine Tenascin X (TXM, [GenBank: AAB82015.1]). The profile HMM files were downloaded from Pfam (EGF-like: PF00008.17, Fn-III: PF00041.11).

Hughes [8] concluded that EGF domain repeats underwent homogenization within each Tenascin gene after duplication, but remained conserved after the divergence of rodents and primates. The same conclusions can be drawn after the evaluation of the diagrams generated with paralogous (Figure 1A and 1B) and orthologous (Figure 1C) sequence pairs. Our diagrams are also in accordance with the conclusions of Hughes regarding the evolution of Fn-III type domain repeats. These repeats can be divided into three categories. The last three C-terminal repeats demonstrate conservation since the duplication of the Tenascin X and Tenascin C genes (Figure 1A', B' and 1C'). Other repeats became homogenized within each gene subsequent to gene duplication, but have remained conserved following the divergence of primates and rodents (Figure 1C'). The repeats of the third category have evolved in a concerted fashion in rodent and primate lineages since their divergence (Figure 1C').
https://static-content.springer.com/image/art%3A10.1186%2F1471-2105-9-27/MediaObjects/12859_2007_Article_2012_Fig1_HTML.jpg
Figure 1

Pairwise repeat homology diagrams of selected Tenascin proteins (SVG outputs). The diagrams were produced with the default PhyML parameters of the script (WAG+Γ+I, 4 gamma categories, gamma parameter and proportion of invariable sites estimated by ML, BIONJ starting tree) and with colour gradient paramater w set to 2. The repeats (EGF on A, B and C; Fn-III on A', B' and C') are indicated by red rectangles in the protein sequence schemes (N-termini on top). The identified homology relations are represented by lines or arcs connecting the respective repeats (blue lines between repeats of different protein sequences (orthology), and brown arcs in case of internal relations (paralogy)). The colour scale bar demonstrates the colour intensities as a function of the patristic distance between the respective clades divided by the total tree length. The lines represent the values between 0 and 1 by units of 0.1. The sequence positions are shown on the right of each scheme. Abbreviations: TCH: Homo sapiens Tenascin C; TXH: Homo sapiens Tenascin X and TXM: Mus musculus Tenascin X.

Conclusion

Although pairwise repeat homology diagrams do not carry all the information about the phylogenetic tree on which they are based, by visualizing the exact positions of the repeats and the homology relations, they permit a rapid and intuitive assessment of the patterns of repeat evolution (compare Figure 1C' with Figure 2). These features make t2prhd a powerful tool, especially in cases of massive datasets, as in studies of repeat evolution in large gene families.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2105-9-27/MediaObjects/12859_2007_Article_2012_Fig2_HTML.jpg
Figure 2

Phylogenetic tree of Fn-III repeats of Tenascin X proteins of Homo sapiens and Mus musculus. Maximum Likelihood tree produced during the analysis of human versus murine Tenascin X sequences with Fn-III profile HMM (diagram on Figure 1C'). Repeats are identified by starting and ending positions. Abbreviations: see Figure 1.

Availability and requirements

  • Project name: t2prhd

  • Project home page: http://t2prhd.sourceforge.net

  • Online access: http://t2prhd.brc.hu

  • Operating system(s): OS Independent (Written in an interpreted language)

  • Programming language: Perl

  • Other requirements: Perl version 5.8.8 or higher with the standard modules, BioPerl modules (version 1.4.0 or higher), XML::Writer module, HMMER package 2.3.2 or higher, CLUSTAL W version 1.83 and/or PhyML version 2.4.4 or higher

  • License: GNU General Public License

  • Any restrictions to use by non-academics: none

Declarations

Acknowledgements

Joint doctoral fellowship of The National Centre for Scholarships Abroad of the Romanian Ministry of Education and Research and the Hungarian Scholarship Board to BS; FEBS Long-term Fellowship to KS; Bolyai Research Scholarship of the Hungarian Academy of Sciences to ZP. This work was supported by the Hungarian Scientific Research Fund grant (NI60442). Further support for the IT infrastructure was provided by the OMAI DT-D/4/005/2006 and BIO-00126 funds of the Hungarian Ministry of Education. We would like to thank for the comments and suggestions to the anonymous reviewers, for useful discussion on the source code to Szilárd Páll and for the assistance in online hosting to Ernő Homolya.

Authors’ Affiliations

(1)
Institute of Genetics, Biological Research Center of the Hungarian Academy of Sciences
(2)
Department of Ecology, University of Szeged

References

  1. Nei M, Rooneyi AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet 2005, 39: 121–52. 10.1146/annurev.genet.39.073003.112240PubMed CentralView ArticlePubMedGoogle Scholar
  2. Liao D: Concerted evolution: molecular mechanism and biological implications. Am J Hum Genet 1999, 64: 24–30. 10.1086/302221PubMed CentralView ArticlePubMedGoogle Scholar
  3. Paoletti M, Saupe SJ, Clavé C: Genesis of a fungal non-self recognition repertoire. PLoS ONE 2007, 2: e283. 10.1371/journal.pone.0000283PubMed CentralView ArticlePubMedGoogle Scholar
  4. Carmon A, M W, Hassan J, Baron M, Macintyre R: Concerted evolution within the Drosophila dumpy gene. Genetics 2007, 176: 309–25. 10.1534/genetics.106.060897PubMed CentralView ArticlePubMedGoogle Scholar
  5. Johannesson H, Townsend JP, Hung CY, Cole GT, Taylor JW: Concerted evolution in the repeats of an immunomodulating cell surface protein, SOWgp, of the human pathogenic fungi Coccidioides immitis and C. posadasii . Genetics 2005, 171: 109–17. 10.1534/genetics.105.040923PubMed CentralView ArticlePubMedGoogle Scholar
  6. Galindo BE, Vacquier VD, Swanson WJ: Positive selection in the egg receptor for abalone sperm lysin. Proc Natl Acad Sci USA 2003, 100: 4639–43. 10.1073/pnas.0830022100PubMed CentralView ArticlePubMedGoogle Scholar
  7. Meeds T, Lockard E, Livingston BT: Special evolutionary properties of genes encoding a protein with a simple amino acid repeat. J Mol Evol 2001, 53: 180–90. 10.1007/s002390010208View ArticlePubMedGoogle Scholar
  8. Hughes AL: Concerted evolution of exons and introns in the MHC-linked tenascin-X gene of mammals. Mol Biol Evol 1999, 16: 1558–67.View ArticlePubMedGoogle Scholar
  9. Swanson WJ, Vacquier VD: Concerted evolution in an egg receptor for a rapidly evolving abalone sperm protein. Science 1998, 281: 710–2. 10.1126/science.281.5377.710View ArticlePubMedGoogle Scholar
  10. Cvitanich C, Salcido M, Judelson HS: Concerted evolution of a tandemly arrayed family of mating-specific genes in Phytophthora analyzed through inter- and intraspecific comparisons. Mol Genet Genomics 2006, 275: 169–84. 10.1007/s00438-005-0074-8View ArticlePubMedGoogle Scholar
  11. Frankel N, Carrari F, Hasson E, Iusem ND: Evolutionary history of the Asr gene family. Gene 2006, 378: 74–83. 10.1016/j.gene.2006.05.010View ArticlePubMedGoogle Scholar
  12. Kwakkenbosi MJ, Matmati M, Madsen O, Pouwels W, Wang Y, Bontrop RE, Heidt PJ, Hoek RM, Hamann J: An unusual mode of concerted evolution of the EGF-TM7 receptor chimera EMR2. FASEB J 2006, 20: 2582–4. 10.1096/fj.06-6500fjeView ArticleGoogle Scholar
  13. Quesada H, Ramos-Onsins SE, Aguadé M: Birth-and-death evolution of the Cecropin multigene family in Drosophila . J Mol Evol 2005, 60: 1–11. 10.1007/s00239-004-0053-4View ArticlePubMedGoogle Scholar
  14. Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci USA 1997, 94: 7799–806. 10.1073/pnas.94.15.7799PubMed CentralView ArticlePubMedGoogle Scholar
  15. HMMER: biosequence analysis using profile hiddenMarkov models[http://hmmer.janelia.org]
  16. Stajich JE, et al.: The Bioperl toolkit: Perl modules for the life sciences. Genome Res 2002, 12: 1611–8. 10.1101/gr.361602PubMed CentralView ArticlePubMedGoogle Scholar
  17. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position specific gap penalties and weight matrix choice. Nucl Acids Res 1994, 22: 4673–80. 10.1093/nar/22.22.4673PubMed CentralView ArticlePubMedGoogle Scholar
  18. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52: 696–704. 10.1080/10635150390235520View ArticlePubMedGoogle Scholar
  19. Firefox[http://www.mozilla.com/en-US/firefox]
  20. Batik SVG Toolkit[http://xmlgraphics.apache.org/batik/]
  21. t2prhd at brc.hu[http://t2prhd.brc.hu]
  22. Letondal C: A Web interface generator for molecular biology programs in Unix. Bioinformatics 2001, 17: 73–82. 10.1093/bioinformatics/17.1.73View ArticlePubMedGoogle Scholar

Copyright

© Sipos et al; licensee BioMed Central Ltd. 2008

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.