|
|
|
|
Published online before print
February 14, 2005, 10.1101/gr.3266405 Genome Res. 15:343-351, 2005 ©2005 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/05 $5.00
Complex genomic rearrangements lead to novel primate gene function1 European Molecular Biology Laboratory, 69012 Heidelberg, Germany 2 Max-Delbrueck-Centrum, D-13092 Berlin, Germany
Orthologous genes that maintain a single-copy status in a broad range of species may indicate a selection against gene duplication. If this is the case, then duplicates of such genes that do survive may have escaped the dosage control by rapid and sizable changes in their function. To test this hypothesis and to develop a strategy for the identification of novel gene functions, we have analyzed 22 primate-specific intrachromosomal duplications of genes with a single-copy ortholog in all other completely sequenced metazoans. When comparing this set to genes not exposed to the single-copy status constraint, we observed a higher tendency of the former to modify their gene structure, often through complex genomic rearrangements. The analysis of the most dramatic of these duplications, affecting 10% of human Chromosome 2, enabled a detailed reconstruction of the events leading to the appearance of a novel gene family. The eight members of this family originated from the highly conserved nucleoporin RanBP2 by several genetic rearrangements such as segmental duplications, inversions, translocations, exon loss, and domain accretion. We have experimentally verified that at least one of the newly formed proteins has a cellular localization different from RanBP2's, and we show that positive selection did act on specific domains during evolution.
Gene duplication is known to play a leading role in evolution for the creation of novel gene function (Ohno 1970
The maintenance of a single-copy status across diverse metazoans is likely to be an indicator for such a selection against duplication, considering the high frequency of duplications in eukaryotes (Lynch and Conery 2000 In order to apply this strategy to human, we searched for primate-specific intrachromosomal duplications of genes with single orthologs in all other metazoans. The detailed study of the most dramatic of such duplications revealed that the break from a well-conserved genetic stability is, indeed, linked to the birth of new gene function and is caused by complex genetic mechanisms. Using comparative genomics, we could reconstruct the genomic events leading to the emergence of a novel human gene function.
Detection of primate-specific gene duplications We systematically searched for single-copy orthologs in all completely sequenced metazoan genomes except human, and then identified their multicopy orthologs located in the same human chromosome (inparalogs) (Table 1). We restricted the analysis to intrachromosomal duplications because these events are expected to be more recent when compared to interchromosomal translocations (D. Torrents, M. Suyama, and P. Bork, unpubl.), and hence more likely to have occurred after the rodent-primate split. Furthermore, the inclusion of orthologs located in different chromosomes would have affected the specificity of our analysis, owing to the possible inclusion of processed pseudogenes. In total we detected 22 metazoan single-copy orthologs with at least one additional paralog in human. Of those we excluded four identical copies that might represent copy number polymorphisms (Sebat et al. 2004
In order to refine the duplication events and to predict the putative functional role of the duplicated genes, we reconstructed the functional relationships between the human copies and their orthologs in the other species. We thus compared the score of the global alignment within the human paralogs and between them and the corresponding mouse ortholog. Higher scores between one of the human copies and the mouse ortholog than within the human paralogs should reveal divergent paralogs in terms of rearrangements in the gene structure. Thus, depending on this score, we defined two categories for the classification of the observed duplications. When the pairwise alignment score is higher within the human paralogs, it is likely that the duplications occurred recently and were not affected by genomic rearrangements. These gene copies might not have been subjected to functional selection, might differ only in their expression patterns (Gu et al. 2004 Of the 18 human paralogous families recorded here, six (33%) fall into this group, while in the remaining 12 duplications (67%) one of the human genes aligns globally better with the mouse ortholog than with the other paralogs. To test whether our data set is, indeed, enriched in genes with a modified structure, we performed a comparable analysis on a reference set obtained by relaxing the constraint of a single-copy status across metazoa. For this purpose, we required a single ortholog in rodents but allowed complete freedom elsewhere, meaning gene absence or presence of multiple paralogs in other metazoans. In this case we detected 215 duplications, and in 55% of them the human copies align better with the mouse ortholog than with the other paralogs. It should be noted that the reference set is likely to include ancient duplications that either underwent rodent-specific loss or whose rodent paralogs were overlooked in the database. This might inflate the final count by increasing the number of human paralogs more similar to the mouse orthologs than to the other copies. Despite this possible bias, we found an enrichment in divergent human paralogs in the data set of genes that have originated from metazoan single-copy orthologs. This supports our hypothesis that there are considerable functional differences among this specific set of novel genes. The duplicated human genes with a modified structure are all candidates for the rapid evolution of new functions through rapid genomic rearrangements. As a case study to test the validity of our strategy, we analyzed in detail the most dramatic of these duplications, interspersed in more than 10% of human Chromosome 2.
A novel gene family descended from RanBP2 and with a domain contribution from GCC2 Single-copy orthologs of RanBP2 are unambiguously detectable in all fully sequenced animal genomes, but not in other eukaryotes such as plants and fungi (Fig. 1A). Although there are slight differences in the domain architectures among the species, the overall domain organization of the protein is well conserved (Fig. 1A). The N-terminal leucine-rich region is followed by RanGTP-binding domains (RanBP1 homologous domains or RBD) related to the one present in RanBP1, a variable number of zinc-finger motifs, and a C terminus with homology to cyclophilin A (Fig. 1A).
In addition to the human RanBP2 ortholog, we detected eight partial copies located in regions that have arisen from intrachromosomal segmental duplications (Figs. 1B, 2A). Despite the fact that human Chromosome 2 appears to be relatively poor in segmental duplications (Bailey et al. 2002 , Imh1p, and p230/golgin-245) domain (Munro and Nichols 1999
Genomic rearrangements led to the formation of the new gene Although the current coverage and quality of the chimp genome do not allow a final assessment, we were able to detect at least four partial genes of the RGP family in addition to the complete RanBP2 ortholog. This is an indication that at least some of the RGP duplications occurred before the human-chimp split. A further confirmation derives from the recent publication of hominoid-specific gene duplications detected by cDNA-array-based comparative genomic hybridization, in which the cDNA corresponding to RanBP2L1 isoform 2 cohybridizes with the genomic DNA of several hominoids (Fortna et al. 2004 200 kb already contains the gene structure of the newly formed RGP family, which evolved from RanBP2 and GCC2 by successive rearrangements. Although it is not possible to assess the temporal order of the events, the ancestor locus underwent duplication, inversion, partial deletion of the long RanBP2 exon 20, and acquisition of the 3'-end of the GCC2 gene coding for the GRIP domain (see Fig. 2B legend for more details). The resulting progenitor locus contains the core coding and noncoding regions common to all the duplicated copies. We performed a similar comparative analysis using the surrounding genomic regions of the eight RGP copies to reconstruct the entire duplication scenario and to assess the parental relationships between them (Fig. 2B). The reliability of this reconstruction is confirmed by the agreement in topology between the trees obtained by using both the gene order approach (Fig. 2B) and the cDNA sequences of RanBP2 and RGPs (Fig. 1B).
The region surrounding the junction site of each of the duplicated segments is enriched in DNA repeats (55%) when compared to a random control (39%) (Supplemental Fig. S1a; Supplemental Table S1; see the legend to the figure for the procedure used to detect the content in DNA repeats). The repeats that mainly contribute to the enrichment are the Alu elements with a peak localized at the junction site (Supplemental Fig. S1b). These results are in agreement with the proposed mechanism for segmental duplications involving Alu-mediated recombinations (Babcock et al. 2003
Alternative splicing of the RGP genes
Positive selection drives the evolution of the RGP family
Five of the seven residues of exon 20 that are potentially under positive selection belong to the RBD domains of the RGP proteins (underlined in Table 2). Four of them correspond to residues directly involved in the binding of RanBP2 to Ran (Vetter et al. 1999
Distinct cellular localization of the RGP proteins
In this study we report primate-specific duplications of genes that have maintained a single-copy status in all other metazoans sequenced so far. Our hypothesis is that the escape from an evolutionarily conserved single-copy status in one specific lineage is connected to the acquisition of novel molecular functions. In addition, we expect that this functional divergence is acquired through massive genetic rearrangements over a relatively short period of time so as to avoid negative dosage effects. Indeed, in the set of 22 recent duplications with a single-copy ortholog in all other metazoans, we observe an enrichment of genes with a modified structure, compared to human duplications not exposed to the same evolutionary constraints. Further confirming the validity of our strategy, we were also able to detect previously reported acquisitions of primate-specific functions by gene duplication. An example is the TRE2 gene (Table 1), which derived from the chimeric fusion of the duplicates of two parental genes (USP32 and TBC1D3) and was recently described as a novel hominoid-specific gene (Paulding et al. 2003 Lineage-specific duplications with novel gene structures not only hint at new function, but also allow the tracing of the evolutionary events leading to the actual genomic arrangements. This is particularly true for primate-specific duplications as they are recent and more likely to be still contained in larger segments. Thus, the comparison of the paralogous and orthologous regions permits the identification of the ancestor locus and allows a detailed reconstruction of the evolutionary mechanism.
Both the discovery of new gene function and the reconstruction of the underlying mechanism were applicable to the most dramatic of the primate-specific duplications reported here, which involves the nucleoporin RanBP2. Diverse genetic events were apparently required over a short time period to allow the emergence of the RGP gene family (Figs. 1, 2). Moreover, within the duplicated segments, RGPs represent the vast majority of the surviving coding regions, indicating that they drive the evolution of the entire region. Finally, the cytoplasmic localization of one of the RGP variants, originally termed RanBP2L1 isoform 2 (Nothwang et al. 1998 Why did RGP appear in primates and nowhere else? There are two possible answers to this question depending on the evolutionary scenario considered. Under the more neutral one, the high number of complex rearrangements required to escape the dosage imbalance renders this event so unlikely that it would only happen in a single lineage. The alternative, and more selective scenario, requires that the positive contribution of the newly formed genes to lineage-specific traits outweighs the negative impact of the dosage imbalance. In this case, the new genes that we detected are likely to have evolved in primates because of their specific genomic and functional context. As is usual in evolutionary biology, it is difficult to determine which one of these scenarios is correct; however, if the selective one is true, a thorough functional characterization of the new genes reported here would reveal the functionally most relevant areas for primate evolution.
Detection of primate-specific gene duplications of single-copy genes The complete protein sets of Anopheles gambiae (BDGP build 2a), Caenorhabditis briggsae (Cb25.agp8), Caenorhabditis elegans (Wormbase build 102), Drosophila melanogaster (BDGP build 3a), Fugu rubripes (Fugu build 2.0), Homo sapiens (NCBI build 34), Mus musculus (NCBI build 30), and Rattus norvegicus (BGDP build 3.1) were downloaded from the Ensembl Web site (http://www.ensembl.org, Dec 20th 2003), and reduced to one translation at each locus (longest transcript). All-against-all Smith-Waterman searches and a previously described algorithm (Tatusov et al. 1997
Calculation of the global alignment score and analysis of the primate-specific duplications
RanBP2 orthology assignment, sequence analysis, and gene structure definition
Analysis of the genomic DNA sequence and detection of the duplicated fragments
Evolutionary analysis
The Ka/Ks ratios were measured using the ML method implemented in PAML (http://abacus.gene.ucl.ac.uk/software/paml.html/; Yang 1997
Amplification of cDNAs corresponding to the RGP gene family
Localization of RGP proteins HeLa cells were grown on coverslips in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum, 100 U/mL penicillin, and 100 µg/mL streptomycin. Transfections were performed in six-well plates with Lipofectamine Plus reagent (Invitrogen), according to the manufacturer's instructions. Then, 24 h after transfection, cells were washed in PBS, fixed for 10 min in 3.7% paraformaldehyde in PBS, washed again in PBS, and stained with a polyclonal anti-calnexin antibody (Santa Cruz). Cells were then mounted in Fluoromount G (Southern Biotechnology). Images were acquired using a Zeiss LSM510 FCS confocal microscope.
We thank Elena Conti for the RBD structure analysis, Sean Hooper for his help with the Ka/Ks calculations, and David Torrents for helpful discussions.
3 Corresponding author. E-mail peer.bork{at}embl.de; fax 49 6221 387 517. [Supplemental material is available online at www.genome.org. and at http://www.bork.embl.de/~ciccarel/RGP_add_data.html.] Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.3266405. Article published online ahead of print in February 2005.
Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W., and Lipman, D.J. 1997. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25: 3389-3402.
Babcock, M., Pavlicek, A., Spiteri, E., Kashork, C.D., Ioshikhes, I., Shaffer, L.G., Jurka, J., and Morrow, B.E. 2003. Shuffling of genes within low-copy repeats on 22q11 (LCR22) by Alu-mediated recombination events during evolution. Genome Res. 13: 2519-2532.
Bailey, J.A., Gu, Z., Clark, R.A., Reinert, K., Samonte, R.V., Schwartz, S., Adams, M.D., Myers, E.W., Li, P.W., and Eichler, E.E. 2002. Recent segmental duplications in the human genome. Science 297: 1003-1007. Bailey, J.A., Liu, G., and Eichler, E.E. 2003. An Alu transposition model for the origin and expansion of human segmental duplications. Am. J. Hum. Genet. 73: 823-834.[CrossRef][Medline]
Birney, E., Thompson, J.D., and Gibson, T.J. 1996. PairWise and SearchWise: Finding the optimal alignment in a simultaneous comparison of a protein profile against all DNA translation frames. Nucleic Acids Res. 24: 2730-2739. Cai, Y., Gao, Y., Sheng, Q., Miao, S., Cui, X., Wang, L., Zong, S., and Koide, S.S. 2002. Characterization and potential function of a novel testis-specific nucleoporin BS-63. Mol. Reprod. Dev. 61: 126-134.[CrossRef][Medline]
Courseaux, A. and Nahon, J.L. 2001. Birth of two chimeric genes in the Hominidae lineage. Science 291: 1293-1297. Eichler, E.E. 2001. Recent duplication, domain accretion and the dynamic mutation of the human genome. Trends Genet. 17: 661-669.[CrossRef][Medline]
Fan, Y., Linardopoulou, E., Friedman, C., Williams, E., and Trask, B.J. 2002a. Genomic structure and evolution of the ancestral chromosome fusion site in 2q13-2q14.1 and paralogous regions on other human chromosomes. Genome Res. 12: 1651-1662.
Fan, Y., Newman, T., Linardopoulou, E., and Trask, B.J. 2002b. Gene content and function of the ancestral chromosome fusion site in human chromosome 2q13-2q14.1 and paralogous regions. Genome Res. 12: 1663-1672. Fortna, A., Kim, Y., MacLaren, E., Marshall, K., Hahn, G., Meltesen, L., Brenton, M., Hink, R., Burgers, S., Hernandez-Boussard, T., et al. 2004. Lineage-specific gene duplication and loss in human and great ape evolution. PLoS Biol. 2: E207.[CrossRef][Medline] Gu, Z., Rifkin, S.A., White, K.P., and Li, W.H. 2004. Duplicate genes increase gene expression diversity within and between species. Nat. Genet. 36: 577-579.[CrossRef][Medline]
Guy, J., Hearn, T., Crosier, M., Mudge, J., Viggiano, L., Koczan, D., Thiesen, H.J., Bailey, J.A., Horvath, J.E., Eichler, E.E., et al. 2003. Genomic sequence and transcriptional profile of the boundary between pericentromeric satellites and genes on human chromosome arm 10p. Genome Res. 13: 159-172.
Hooper, S.D. and Berg, O.G. 2003. On the nature of gene innovation: Duplication patterns in microbial genomes. Mol. Biol. Evol. 20: 945-954. Johnson, M.E., Viggiano, L., Bailey, J.A., Abdul-Rauf, M., Goodwin, G., Rocchi, M., and Eichler, E.E. 2001. Positive selection of a gene family during the emergence of humans and African apes. Nature 413: 514-519.[CrossRef][Medline]
Kent, W.J. 2002. BLATThe BLAST-like alignment tool. Genome Res. 12: 656-664. Kimura. 1983. The neutral theory of molecular evolution. Cambridge University Press, Cambridge, UK.
Kumar, S., Tamura, K., Jakobsen, I.B., and Nei, M. 2001. MEGA2: Molecular evolutionary genetics analysis software. Bioinformatics 17: 1244-1245.
Letunic, I., Copley, R.R., Schmidt, S., Ciccarelli, F.D., Doerks, T., Schultz, J., Ponting, C.P., and Bork, P. 2004. SMART 4.0: Towards genomic data integration. Nucleic Acids Res. 32: D142-D144. Li, W.H., Gu, Z., Wang, H., and Nekrutenko, A. 2001. Evolutionary analyses of the human genome. Nature 409: 847-849.[CrossRef][Medline] Loftus, B.J., Kim, U.J., Sneddon, V.P., Kalush, F., Brandon, R., Fuhrmann, J., Mason, T., Crosby, M.L., Barnstead, M., Cronin, L., et al. 1999. Genome duplications and other features in 12 Mb of DNA sequence from human chromosome 16p and 16q. Genomics 60: 295-308.[CrossRef][Medline] Long, M., Betran, E., Thornton, K., and Wang, W. 2003. The origin of new genes: Glimpses from the young and old. Nat. Rev. Genet. 4: 865-875.[Medline]
Luke, M.R., Kjer-Nielsen, L., Brown, D.L., Stow, J.L., and Gleeson, P.A. 2003. GRIP domain-mediated targeting of two new coiled-coil proteins, GCC88 and GCC185, to subcompartments of the trans-Golgi network. J. Biol. Chem. 278: 4216-4226. Lupski, J.R., Wise, C.A., Kuwano, A., Pentao, L., Parke, J.T., Glaze, D.G., Ledbetter, D.H., Greenberg, F., and Patel, P.I. 1992. Gene dosage is a mechanism for Charcot-Marie-Tooth disease type 1A. Nat. Genet. 1: 29-33.[CrossRef][Medline]
Lynch, M. and Conery, J.S. 2000. The evolutionary fate and consequences of duplicate genes. Science 290: 1151-1155.
Maston, G.A. and Ruvolo, M. 2002. Chorionic gonadotropin has a recent origin within primates and an evolutionary history of selection. Mol. Biol. Evol. 19: 320-335. McConville, M.J., Ilgoutz, S.C., Teasdale, R.D., Foth, B.J., Matthews, A., Mullin, K.A., and Gleeson, P.A. 2002. Targeting of the GRIP domain to the trans-Golgi network is conserved from protists to animals. Eur. J. Cell Biol. 81: 485-495.[CrossRef][Medline] Munro, S. and Nichols, B.J. 1999. The GRIP domainA novel Golgi-targeting domain found in several coiled-coil proteins. Curr. Biol. 9: 377-380.[CrossRef][Medline] Nadeau, J.H. and Sankoff, D. 1997. Comparable rates of gene loss and functional divergence after genome duplications early in vertebrate evolution. Genetics 147: 1259-1266.[Abstract]
Nielsen, R. and Yang, Z. 1998. Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148: 929-936. Nothwang, H.G., Rensing, C., Kubler, M., Denich, D., Brandl, B., Stubanus, M., Haaf, T., Kurnit, D., and Hildebrandt, F. 1998. Identification of a novel Ran binding protein 2 related gene (RANBP2L1) and detection of a gene cluster on human chromosome 2q11-q12. Genomics 47: 383-392.[CrossRef][Medline] Ohno, S. 1970. Evolution by gene duplication. Springer-Verlag, Berlin-Heidelberg-New York. Papp, B., Pal, C., and Hurst, L.D. 2003. Dosage sensitivity and the evolution of gene families in yeast. Nature 424: 194-197.[CrossRef][Medline]
Paulding, C.A., Ruvolo, M., and Haber, D.A. 2003. The Tre2 (USP6) oncogene is a hominoid-specific gene. Proc. Natl. Acad. Sci. 100: 2507-2511. Prince, V.E. and Pickett, F.B. 2002. Splitting pairs: The diverging fates of duplicated genes. Nat. Rev. Genet. 3: 827-837.[CrossRef][Medline] Samonte, R.V. and Eichler, E.E. 2002. Segmental duplications and the evolution of the primate genome. Nat. Rev. Genet. 3: 65-72.[Medline]
Sebat, J., Lakshmi, B., Troge, J., Alexander, J., Young, J., Lundin, P., Maner, S., Massa, H., Walker, M., Chi, M., et al. 2004. Large-scale copy number polymorphism in the human genome. Science 305: 525-528. She, X., Jiang, Z., Clark, R.A., Liu, G., Cheng, Z., Tuzun, E., Church, D.M., Sutton, G., Halpern, A.L., and Eichler, E.E. 2004. Shotgun sequence assembly and recent segmental duplications within the human genome. Nature 431: 927-930.[CrossRef][Medline]
Suyama, M., Torrents, D., and Bork, P. 2004. BLAST2GENE: A comprehensive conversion of BLAST output into independent genes and gene fragments. Bioinformatics 20: 1968-1970.
Tatusov, R.L., Koonin, E.V., and Lipman, D.J. 1997. A genomic perspective on protein families. Science 278: 631-637. Tatusova, T.A. and Madden, T.L. 1999. BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol. Lett. 174: 247-250.[CrossRef][Medline]
Thompson, J., Gibson, T., Plewniak, F., Jeanmougin, F., and Higgins, D. 1997. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25: 4876-4882.
Thomson, T.M., Lozano, J.J., Loukili, N., Carrio, R., Serras, F., Cormand, B., Valeri, M., Diaz, V.M., Abril, J., Burset, M., et al. 2000. Fusion of the human gene for the polyubiquitination coeffector UEV1 with Kua, a newly identified gene. Genome Res. 10: 1743-1756. Vetter, I.R., Nowak, C., Nishimoto, T., Kuhlmann, J., and Wittinghofer, A. 1999. Structure of a Ran-binding domain complexed with Ran bound to a GTP analogue: Implications for nuclear transport. Nature 398: 39-46.[CrossRef][Medline] Walsh, J.B. 1995. How often do duplicated genes evolve new functions? Genetics 139: 421-428.[Abstract]
Wu, J., Matunis, M.J., Kraemer, D., Blobel, G., and Coutavas, E. 1995. Nup358, a cytoplasmically exposed nucleoporin with peptide repeats, Ran-GTP binding sites, zinc fingers, a cyclophilin A homologous domain, and a leucine-rich region. J. Biol. Chem. 270: 14209-14213. Wu, M., Lu, L., Hong, W., and Song, H. 2004. Structural basis for recruitment of GRIP domain golgin-245 by small GTPase Arl1. Nat. Struct. Mol. Biol. 11: 86-94.[CrossRef][Medline]
Yang, Z. 1997. PAML: A program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13: 555-556.
Yang, Z. and Nielsen, R. 2002. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol. Biol. Evol. 19: 908-917.
Yang, Z., Nielsen, R., Goldman, N., and Pedersen, A.M. 2000. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155: 431-449.
Yang, J., Lusk, R., and Li, W.-H. 2003. Organismal complexity, protein complexity, and gene duplicability. Proc. Natl. Acad. Sci. 100: 15661-15665. Yokoyama, N., Hayashi, N., Seki, T., Pante, N., Ohba, T., Nishii, K., Kuma, K., Hayashida, T., Miyata, T., Aebi, U., et al. 1995. A giant nucleopore protein that binds Ran/TC4. Nature 376: 184-188.[CrossRef][Medline]
Yoshino, A., Bieler, B.M., Harper, D.C., Cowan, D.A., Sutterwala, S., Gay, D.M., Cole, N.B., McCaffery, J.M., and Marks, M.S. 2003. A role for GRIP domain proteins and/or their ligands in structure and function of the trans Golgi network. J. Cell Sci. 116: 4441-4454.
Yunis, J.J. and Prakash, O. 1982. The origin of man: A chromosomal pictorial legacy. Science 215: 1525-1530.
Zdobnov, E.M., von Mering, C., Letunic, I., Torrents, D., Suyama, M., Copley, R.R., Christophides, G.K., Thomasova, D., Holt, R.A., Subramanian, G.M., et al. 2002. Comparative genome and proteome analysis of Anopheles gambiae and Drosophila melanogaster. Science 298: 149-159. Zhang, J., Zhang, Y.P., and Rosenberg, H.F. 2002. Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat. Genet. 30: 411-415.[CrossRef][Medline]
http://genome.ucsc.edu/; University of California San Francisco genome browser. http://www.ensembl.org/; Ensembl. http://www.bork.embl.de/~ciccarel/RGP_add_data.html; Supplemental material to this paper. http://abacus.gene.ucl.ac.uk/software/paml.html/; PAML. http://smart.embl-heidelberg.de/; SMART database. http://www.megasoftware.net/; MEGA. http://www.ncbi.nlm.nih.gov/; National Center for Biotechnology Information.
Received September 20, 2004; accepted in revised format December 15, 2004. This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||