51
|
Kochetov AV. Alternative translation start sites and hidden coding potential of eukaryotic mRNAs. Bioessays 2008; 30:683-91. [PMID: 18536038 DOI: 10.1002/bies.20771] [Citation(s) in RCA: 136] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
It is widely suggested that a eukaryotic mRNA typically contains one translation start site and encodes a single functional protein product. However, according to current points of view on translation initiation mechanisms, eukaryotic ribosomes can recognize several alternative translation start sites and the number of experimentally verified examples of alternative translation is growing rapidly. Also, the frequent occurrence of alternative translation events and their functional significance are supported by the results of computational evaluations. The functional role of alternative translation and its contribution to eukaryotic proteome complexity are discussed.
Collapse
|
52
|
Tran MK, Schultz CJ, Baumann U. Conserved upstream open reading frames in higher plants. BMC Genomics 2008; 9:361. [PMID: 18667093 PMCID: PMC2527020 DOI: 10.1186/1471-2164-9-361] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2008] [Accepted: 07/31/2008] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Upstream open reading frames (uORFs) can down-regulate the translation of the main open reading frame (mORF) through two broad mechanisms: ribosomal stalling and reducing reinitiation efficiency. In distantly related plants, such as rice and Arabidopsis, it has been found that conserved uORFs are rare in these transcriptomes with approximately 100 loci. It is unclear how prevalent conserved uORFs are in closely related plants. RESULTS We used a homology-based approach to identify conserved uORFs in five cereals (monocots) that could potentially regulate translation. Our approach used a modified reciprocal best hit method to identify putative orthologous sequences that were then analysed by a comparative R-nomics program called uORFSCAN to find conserved uORFs. CONCLUSION This research identified new genes that may be controlled at the level of translation by conserved uORFs. We report that conserved uORFs are rare (<150 loci contain them) in cereal transcriptomes, are generally short (less than 100 nt), highly conserved (50% median amino acid sequence similarity), position independent in their 5'-UTRs, and their start codon context and the usage of rare codons for translation does not appear to be important.
Collapse
Affiliation(s)
- Michael K Tran
- Australian Centre for Plant Functional Genomics PMB 1 Glen Osmond SA 5064, Australia.
| | | | | |
Collapse
|
53
|
Karagyozov L, Godfrey R, Böhmer SA, Petermann A, Hölters S, Ostman A, Böhmer FD. The structure of the 5'-end of the protein-tyrosine phosphatase PTPRJ mRNA reveals a novel mechanism for translation attenuation. Nucleic Acids Res 2008; 36:4443-53. [PMID: 18603590 PMCID: PMC2490741 DOI: 10.1093/nar/gkn391] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Analysis of the human protein-tyrosine phosphatase (PTP) PTPRJ mRNA detected three in-frame AUGs at the 5′-end (starting at nt +14, +191 and +356) with no intervening stop codons. This tandem AUG arrangement is conserved between humans and the mouse and is unique among the genes of the classical PTPs. Until now it was assumed that the principal open reading frame (ORF) starts at AUG356. Our experiments showed that: (i) translation of the mRNA synthesized under the PTPRJ promoter starts predominantly at AUG191, leading to the generation of a 55 amino acid sequence preceding the signal peptide; (ii) the longer form is being likewise correctly processed into mature PTPRJ; (iii) the translation of the region between AUG191 and AUG356 inhibits the overall expression, a feature which depends on the sequence of the encoded peptide. Specifically, a sequence of 13 amino acids containing multiple arginine residues (RRTGWRRRRRRRR) confers the inhibition. In the absence of uORF these previously unrecognized characteristics of the 5′-end of the mRNA present a novel mechanism to suppress, and potentially to regulate translation.
Collapse
Affiliation(s)
- Luchezar Karagyozov
- Institute of Molecular Cell Biology, Center for Molecular Biomedicine, Friedrich-Schiller-University Jena, Jena, Germany
| | | | | | | | | | | | | |
Collapse
|
54
|
Hayden CA, Bosco G. Comparative genomic analysis of novel conserved peptide upstream open reading frames in Drosophila melanogaster and other dipteran species. BMC Genomics 2008; 9:61. [PMID: 18237443 PMCID: PMC2276209 DOI: 10.1186/1471-2164-9-61] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2007] [Accepted: 02/01/2008] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Upstream open reading frames (uORFs) are elements found in the 5'-region of an mRNA transcript, capable of regulating protein production of the largest, or major ORF (mORF), and impacting organismal development and growth in fungi, plants, and animals. In Drosophila, approximately 40% of transcripts contain upstream start codons (uAUGs) but there is little evidence that these are translated and affect their associated mORF. RESULTS Analyzing 19,389 Drosophila melanogaster transcript annotations and 666,153 dipteran EST sequences we have identified 44 putative conserved peptide uORFs (CPuORFs) in Drosophila melanogaster that show evidence of negative selection, and therefore are likely to be translated. Transcripts with CPuORFs constitute approximately 0.3% of the total number of transcripts, a similar frequency to the Arabidopsis genome, and have a mean length of 70 amino acids, much larger than the mean length of plant CPuORFs (40 amino acids). There is a statistically significant clustering of CPuORFs at cytological band 57 (p = 10-5), a phenomenon that has never been described for uORFs. Based on GO term and Interpro domain analyses, genes in the uORF dataset show a higher frequency of ORFs implicated in mitochondrial import than the genome-wide frequency (p < 0.01) as well as methyltransferases (p < 0.02). CONCLUSION Based on these data, it is clear that Drosophila contain putative CPuORFs at frequencies similar to those found in plants. They are distinguished, however, by the type of mORF they tend to associate with, Drosophila CPuORFs preferentially occurring in transcripts encoding mitochondrial proteins and methyltransferases. This provides a basis for the study of CPuORFs and their putative regulatory role in mitochondrial function and disease.
Collapse
Affiliation(s)
- Celine A Hayden
- Department of Molecular and Cellular Biology, University of Arizona, Tucson, AZ 85721, USA.
| | | |
Collapse
|
55
|
Puyaubert J, Denis L, Alban C. Dual targeting of Arabidopsis holocarboxylase synthetase1: a small upstream open reading frame regulates translation initiation and protein targeting. PLANT PHYSIOLOGY 2008; 146:478-91. [PMID: 18156294 PMCID: PMC2245827 DOI: 10.1104/pp.107.111534] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2007] [Accepted: 12/13/2007] [Indexed: 05/24/2023]
Abstract
Protein biotinylation is an original and very specific posttranslational modification, compartmented in plants, between mitochondria, plastids, and the cytosol. This reaction modifies and activates few carboxylases committed in key metabolisms and is catalyzed by holocarboxylase synthetase (HCS). The molecular bases of this complex compartmentalization and the relative function of each of the HCS genes, HCS1 and HCS2, identified in Arabidopsis (Arabidopsis thaliana) are mainly unknown. Here, we showed by reverse genetics that the HCS1 gene is essential for plant viability, whereas disruption of the HCS2 gene in Arabidopsis does not lead to any obvious phenotype when plants are grown under standard conditions. These findings strongly suggest that HCS1 is the only protein responsible for HCS activity in Arabidopsis cells, including the cytosolic, mitochondrial, and plastidial compartments. A closer study of HCS1 gene expression enabled us to propose an original mechanism to account for this multiplicity of localizations. Located in the HCS1 messenger RNA 5'-untranslated region, an upstream open reading frame regulates the translation initiation of HCS1 and the subsequent targeting of HCS1 protein. Moreover, an exquisitely precise alternative splicing of HCS1 messenger RNA can regulate the presence and absence of this upstream open reading frame. The existence of these complex and interdependent mechanisms creates a rich molecular platform where different parameters and factors could control HCS targeting and hence biotin metabolism.
Collapse
Affiliation(s)
- Juliette Puyaubert
- CNRS (UMR 5168)/CEA/Université Joseph Fourier/INRA (UMR 1200), CEA-Grenoble, Institut de Recherche en Technologies et Sciences pour le Vivant, Grenoble cedex 9, France
| | | | | |
Collapse
|
56
|
Rosenstiel P, Huse K, Franke A, Hampe J, Reichwald K, Platzer C, Roberts RG, Mathew CG, Platzer M, Schreiber S. Functional characterization of two novel 5' untranslated exons reveals a complex regulation of NOD2 protein expression. BMC Genomics 2007; 8:472. [PMID: 18096043 PMCID: PMC2228316 DOI: 10.1186/1471-2164-8-472] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2007] [Accepted: 12/20/2007] [Indexed: 02/06/2023] Open
Abstract
Background NOD2 is an innate immune receptor for the bacterial cell wall component muramyl-dipeptide. Mutations in the leucine-rich repeat region of NOD2, which lead to an impaired recognition of muramyl-dipeptide, have been associated with Crohn disease, a human chronic inflammatory bowel disease. Tissue specific constitutive and inducible expression patterns of NOD2 have been described that result from complex regulatory events for which the molecular mechanisms are not yet fully understood. Results We have identified two novel exons of the NOD2 gene (designated exon 1a and 1b), which are spliced to the canonical exon 2 and constitute the 5' untranslated region of two alternative transcript isoforms (i.e. exon 1a/1b/2 and exon 1a/2). The two novel transcripts are abundantly expressed and seem to comprise the majority of NOD2 transcripts under physiological conditions. We confirm the expression of the previously known canonical first exon (designated exon 1c) of the gene in unstimulated mononuclear cells. The inclusion of the second alternative exon 1b, which harbours three short upstream open reading frames (uORFs), is downregulated upon stimulation with TNF-α or under pro-inflammatory conditions in the inflamed intestinal mucosa in vivo. Using the different 5' UTR splice forms fused to a firefly luciferase (LUC) reporter we demonstrate a rapamycin-sensitive inhibitory effect of the uORFs on translation efficacy. Conclusion The differential usage of two alternative promoters in the NOD2 gene leads to tissue-specific and context-dependent NOD2 transcript isoform patterns. We demonstrate for the first time that context-dependent alternative splicing is linked to uORF-mediated translational repression. The results suggest complex parallel control mechanisms that independently regulate NOD2 expression in the context of inflammatory signaling.
Collapse
Affiliation(s)
- Philip Rosenstiel
- Institute of Clinical Molecular Biology, University Hospital Schleswig-Holstein, Campus Kiel, Germany.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
57
|
Nakagawa S, Niimura Y, Gojobori T, Tanaka H, Miura KI. Diversity of preferred nucleotide sequences around the translation initiation codon in eukaryote genomes. Nucleic Acids Res 2007; 36:861-71. [PMID: 18086709 PMCID: PMC2241899 DOI: 10.1093/nar/gkm1102] [Citation(s) in RCA: 139] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Understanding regulatory mechanisms of protein synthesis in eukaryotes is essential for the accurate annotation of genome sequences. Kozak reported that the nucleotide sequence GCCGCC(A/G)CCAUGG (AUG is the initiation codon) was frequently observed in vertebrate genes and that this ‘consensus’ sequence enhanced translation initiation. However, later studies using invertebrate, fungal and plant genes reported different ‘consensus’ sequences. In this study, we conducted extensive comparative analyses of nucleotide sequences around the initiation codon by using genomic data from 47 eukaryote species including animals, fungi, plants and protists. The analyses revealed that preferred nucleotide sequences are quite diverse among different species, but differences between patterns of nucleotide bias roughly reflect the evolutionary relationships of the species. We also found strong biases of A/G at position −3, A/C at position −2 and C at position +5 that were commonly observed in all species examined. Genes with higher expression levels showed stronger signals, suggesting that these nucleotides are responsible for the regulation of translation initiation. The diversity of preferred nucleotide sequences around the initiation codon might be explained by differences in relative contributions from two distinct patterns, GCCGCCAUG and AAAAAAAUG, which implies the presence of multiple molecular mechanisms for controlling translation initiation.
Collapse
Affiliation(s)
- So Nakagawa
- Department of Systems Biology, School of Biomedical Science, Department of Bioinformatics, Medical Research Institute, Tokyo Medical and Dental University, Yushima, Tokyo, Japan
| | | | | | | | | |
Collapse
|
58
|
Cvijović M, Dalevi D, Bilsland E, Kemp GJL, Sunnerhagen P. Identification of putative regulatory upstream ORFs in the yeast genome using heuristics and evolutionary conservation. BMC Bioinformatics 2007; 8:295. [PMID: 17686169 PMCID: PMC1964767 DOI: 10.1186/1471-2105-8-295] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2007] [Accepted: 08/08/2007] [Indexed: 11/30/2022] Open
Abstract
Background The translational efficiency of an mRNA can be modulated by upstream open reading frames (uORFs) present in certain genes. A uORF can attenuate translation of the main ORF by interfering with translational reinitiation at the main start codon. uORFs also occur by chance in the genome, in which case they do not have a regulatory role. Since the sequence determinants for functional uORFs are not understood, it is difficult to discriminate functional from spurious uORFs by sequence analysis. Results We have used comparative genomics to identify novel uORFs in yeast with a high likelihood of having a translational regulatory role. We examined uORFs, previously shown to play a role in regulation of translation in Saccharomyces cerevisiae, for evolutionary conservation within seven Saccharomyces species. Inspection of the set of conserved uORFs yielded the following three characteristics useful for discrimination of functional from spurious uORFs: a length between 4 and 6 codons, a distance from the start of the main ORF between 50 and 150 nucleotides, and finally a lack of overlap with, and clear separation from, neighbouring uORFs. These derived rules are inherently associated with uORFs with properties similar to the GCN4 locus, and may not detect most uORFs of other types. uORFs with high scores based on these rules showed a much higher evolutionary conservation than randomly selected uORFs. In a genome-wide scan in S. cerevisiae, we found 34 conserved uORFs from 32 genes that we predict to be functional; subsequent analysis showed the majority of these to be located within transcripts. A total of 252 genes were found containing conserved uORFs with properties indicative of a functional role; all but 7 are novel. Functional content analysis of this set identified an overrepresentation of genes involved in transcriptional control and development. Conclusion Evolutionary conservation of uORFs in yeasts can be traced up to 100 million years of separation. The conserved uORFs have certain characteristics with respect to length, distance from each other and from the main start codon, and folding energy of the sequence. These newly found characteristics can be used to facilitate detection of other conserved uORFs.
Collapse
Affiliation(s)
- Marija Cvijović
- Department of Cell and Molecular Biology, Lundberg Laboratory, Göteborg University, PO Box 462 SE-405 30 Göteborg, Sweden
- Max-Planck Institute for Molecular Genetics, Ihnestraße 63, D-14195 Berlin, Germany
| | - Daniel Dalevi
- Department of Computer Science and Engineering, Chalmers University of Technology, SE-412 96 Göteborg, Sweden
| | - Elizabeth Bilsland
- Department of Cell and Molecular Biology, Lundberg Laboratory, Göteborg University, PO Box 462 SE-405 30 Göteborg, Sweden
- Biochemistry Department, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK
| | - Graham JL Kemp
- Department of Computer Science and Engineering, Chalmers University of Technology, SE-412 96 Göteborg, Sweden
| | - Per Sunnerhagen
- Department of Cell and Molecular Biology, Lundberg Laboratory, Göteborg University, PO Box 462 SE-405 30 Göteborg, Sweden
| |
Collapse
|
59
|
Matsui M, Yachie N, Okada Y, Saito R, Tomita M. Bioinformatic analysis of post-transcriptional regulation by uORF in human and mouse. FEBS Lett 2007; 581:4184-8. [PMID: 17692847 DOI: 10.1016/j.febslet.2007.07.057] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2007] [Revised: 07/03/2007] [Accepted: 07/24/2007] [Indexed: 01/22/2023]
Abstract
RNA decay is thought to exert an important influence on gene expression by maintaining a steady-state level of transcripts and/or by eliminating aberrant transcripts. However, the sequence elements which control such processes have not been determined. Upstream open reading frames (uORFs) in the transcripts of several genes are reported to control translational initiation by stalling ribosomes and thereby promote RNA decay. We therefore performed bioinformatic analysis of the tissue-wide expression profiles and mRNA half-life of transcripts containing uORFs in humans and mice to assess the relationship between RNA decay and the presence of uORFs in transcripts. The expression levels of transcripts containing uORF were markedly lower than those not containing uORF. Moreover, the half-life of the uORF-containing transcripts was also shorter. These results suggest that uORFs are sequence elements that down-regulate RNA transcripts via RNA decay mechanisms.
Collapse
Affiliation(s)
- Motomu Matsui
- Institute for Advanced Biosciences, Keio University, Tsuruoka 997-0035, Japan
| | | | | | | | | |
Collapse
|
60
|
Hayden CA, Jorgensen RA. Identification of novel conserved peptide uORF homology groups in Arabidopsis and rice reveals ancient eukaryotic origin of select groups and preferential association with transcription factor-encoding genes. BMC Biol 2007; 5:32. [PMID: 17663791 PMCID: PMC2075485 DOI: 10.1186/1741-7007-5-32] [Citation(s) in RCA: 113] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2007] [Accepted: 07/30/2007] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Upstream open reading frames (uORFs) can mediate translational control over the largest, or major ORF (mORF) in response to starvation, polyamine concentrations, and sucrose concentrations. One plant uORF with conserved peptide sequences has been shown to exert this control in an amino acid sequence-dependent manner but generally it is not clear what kinds of genes are regulated, or how extensively this mechanism is invoked in a given genome. RESULTS By comparing full-length cDNA sequences from Arabidopsis and rice we identified 26 distinct homology groups of conserved peptide uORFs, only three of which have been reported previously. Pairwise Ka/Ks analysis showed that purifying selection had acted on nearly all conserved peptide uORFs and their associated mORFs. Functions of predicted mORF proteins could be inferred for 16 homology groups and many of these proteins appear to have a regulatory function, including 6 transcription factors, 5 signal transduction factors, 3 developmental signal molecules, a homolog of translation initiation factor eIF5, and a RING finger protein. Transcription factors are clearly overrepresented in this data set when compared to the frequency calculated for the entire genome (p = 1.2 x 10(-7)). Duplicate gene pairs arising from a whole genome duplication (ohnologs) with a conserved uORF are much more likely to have been retained in Arabidopsis (Arabidopsis thaliana) than are ohnologs of other genes (39% vs 14% of ancestral genes, p = 5 x 10(-3)). Two uORF groups were found in animals, indicating an ancient origin of these putative regulatory elements. CONCLUSION Conservation of uORF amino acid sequence, association with homologous mORFs over long evolutionary time periods, preferential retention after whole genome duplications, and preferential association with mORFs coding for transcription factors suggest that the conserved peptide uORFs identified in this study are strong candidates for translational controllers of regulatory genes.
Collapse
Affiliation(s)
- Celine A Hayden
- Department of Plant Sciences, University of Arizona, Tucson, AZ 85721-0036, USA
| | - Richard A Jorgensen
- Department of Plant Sciences, University of Arizona, Tucson, AZ 85721-0036, USA
| |
Collapse
|
61
|
Roy SW, Penny D, Neafsey DE. Evolutionary conservation of UTR intron boundaries in Cryptococcus. Mol Biol Evol 2007; 24:1140-8. [PMID: 17374879 DOI: 10.1093/molbev/msm045] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Despite significant progress, the general functional and evolutionary significance of the untranslated regions (UTRs) of eukaryotic transcripts remain mysterious. Particularly mysterious is the common occurrence of spliceosomal introns in transcript UTRs because UTR splicing is not necessary for restoration of transcript coding sequence. In general, it is not known to what extent such splicing performs an important function or merely represents spliceosomal "noise." We conducted the first analysis of evolutionary conservation of UTR splicing. Among 4 species from Cryptococcus neoformans species complex, we find high levels of conservation of UTR intron boundary sequences, strongly suggesting that UTR intron splicing is conserved by purifying selection. We estimate that 50-90% of splice boundaries are maintained by selection. Donor site sequences are more highly conserved than acceptor sequences, and splicing boundaries are more conserved in 5' UTRs than in 3' UTRs. In addition, we report a variety of differences between patterns of UTR splicing in Cryptococcus and corresponding patterns in animals and plants. These results focus attention on the functional roles of eukaryotic UTRs and deepen the mystery of UTR intron splicing.
Collapse
Affiliation(s)
- Scott William Roy
- Allan Wilson Centre for Molecular Ecology and Evolution, Massey University, Palmerston North, New Zealand.
| | | | | |
Collapse
|
62
|
Sunnerhagen P. Cytoplasmatic post-transcriptional regulation and intracellular signalling. Mol Genet Genomics 2007; 277:341-55. [PMID: 17333280 DOI: 10.1007/s00438-007-0221-5] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2006] [Accepted: 02/06/2007] [Indexed: 02/05/2023]
Abstract
Studies of intracellular signalling have traditionally focused on regulation at the levels of initiation of transcription on one hand, and post-translational regulation on the other. More recently, it is becoming apparent that the post-transcriptional level of gene expression is also subject to regulation by signalling pathways. The emphasis in this review is on short-term regulation of mRNAs at the levels of degradation and frequency of translation. Interplay between the mRNA translation and degradation machineries and mainly the TOR, stress-induced MAP kinase (SAPK), and DNA damage checkpoint pathways is discussed. Since a large fraction of the molecular mechanisms has been dissected using molecular genetics methods in yeast, most of the examples in this review are from budding and fission yeast. Some parallels are drawn to plant and animal cells. This review is intended for those more familiar with intracellular signalling, and who realise that post-transcriptional regulation may be an underemphasised level of signalling output.
Collapse
Affiliation(s)
- Per Sunnerhagen
- Department of Cell and Molecular Biology, Lundberg Laboratory, Göteborg University, P.O. Box 462, 405 30, Göteborg, Sweden.
| |
Collapse
|
63
|
Abstract
Background The rate of evolution varies spatially along genomes and temporally in time. The presence of evolutionary rate variation is an informative signal that often marks functional regions of genomes and historical selection events. There exist many tests for temporal rate variation, or heterotachy, that start by partitioning sampled sequences into two or more groups and testing rate homogeneity among the groups. I develop a Bayesian method to infer phylogenetic trees with a divergence point, or dramatic temporal shifts in selection pressure that affect many nucleotide sites simultaneously, located at an unknown position in the tree. Results Simulation demonstrates that the method is most able to detect divergence points when rate variation and the number of affected sites is high, but not beyond biologically relevant values. The method is applied to two viral data sets. A divergence point is identified separating the B and C subtypes, two genetically distinct variants of HIV that have spread into different human populations with the AIDS epidemic. In contrast, no strong signal of temporal rate variation is found in a sample of F and H genotypes, two genetic variants of HBV that have likely evolved with humans during their immigration and expansion into the Americas. Conclusion Temporal shifts in evolutionary rate of sufficient magnitude are detectable in the history of sampled sequences. The ability to detect such divergence points without the need to specify a prior hypothesis about the location or timing of the divergence point should help scientists identify historically important selection events and decipher mechanisms of evolution.
Collapse
Affiliation(s)
- Karin S Dorman
- Department of Statistics, and the Program in Bioinformatics and Computational Biology, Iowa State University, Ames, IA, USA.
| |
Collapse
|
64
|
|