Department of Zoology and Entomology, University of Queensland, Brisbane, Australia; and School of Biological Sciences, University of Sussex, Brighton, United Kingdom
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Resolving the controversy over the association between body size and rate of molecular evolution has important implications for understanding and using the molecular clock (the use of measures of DNA divergence to estimate the age of biological lineages). Clearly the rate of molecular evolution is not uniform across all lineages: for example, viruses can have a rate of molecular evolution a million times faster than mammals (Fitch 1996
), and mammalian rates are an order of magnitude faster than sharks (Martin 1999
). However, there has been a tendency to assume a local molecular clockthat closely related species will have approximately the same rate of molecular evolutionin order to use a calibration rate estimated from one lineage to date the origin of other related lineages. But if there is a general body size effect, then even closely related species could differ in substitution rate. Establishing the extent of the body size effect among vertebrate lineages is important in order to critically evaluate the common assumption of a local molecular clock.
The study of patterns of lineage-specific rates of molecular evolution has in the past been hampered by insufficient comparative data and inappropriate analyses. Available gene sequences have been concentrated in a few well-studied taxa. In particular, the most well-studied cases of lineage-specific rates have been for four mammalian species: rat, mouse, human, and cow. Although comparison of a large number of sequences from these species has served to illustrate consistent lineage-specific rates (Li et al. 1990
; Bulmer, Wolfe, and Sharp 1991
; Easteal and Collett 1994
; Yang and Nielsen 1998
), the limited number of taxa has prevented assessment of the generality of lineage-specific rates and has hampered testing of hypothetical causes of rate differences. For example, the consistently higher rates of molecular evolution in rats and mice than in humans could be caused by their shorter generation time, higher metabolic rate, lower rates of DNA repair, or all or none of these possibilities. The limited number of taxa has also led to an anecdotal approach to assessing the patterns and causes of rate variation, pointing to specific examples as evidence for or against a given hypothesis.
Furthermore, many studies seeking correlates of lineage-specific rates have been limited by failure to adequately account for phylogenetic bias in the data (e.g., Martin and Palumbi 1993
; Omland 1997
; Bleiweiss 1998
; Schmitz and Moritz 1998
). Species cannot be considered as independently derived data points in a statistical analysis because patterns of descent create hierarchies of similarity. Including multiple nonindependent comparisons could create a false association between rate and body size. But phylogenetic bias cannot wholly account for the observation of a body size effect on rate of molecular evolution because studies that have controlled for phylogenetic nonindependence have also found an inverse body size effect (Mooers and Harvey 1994
; Bromham, Rambaut, and Harvey 1996
).
Here, a phylogenetic comparative method has been used to demonstrate a significant association between life history and rate of DNA sequence evolution for a range of reptile lineages for three genes. The power of this study is limited by the available datait is difficult to find sufficient comparisons between reptile species for which appropriate sequence data, life history, and phylogenetic information are all availablebut the significant association between life history and rate of molecular evolution for the genes tested in this study suggests that this may be a general phenomenon in reptiles. As more data become available, the precise nature of the relationship, its taxonomic generality, and the mechanisms behind the pattern may be illuminated.
![]() |
Data |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
![]() |
Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
DNA sequences used in this study were taken from GenBank and aligned by eye using Se-Al (Rambaut 1996
). Three distance matrices were estimated from each alignment. Genetic distance (D) was estimated using the HKY85 model, with transition-transversion ratio estimated from the data (implemented in PAUP*; Swofford 1999
). Synonymous (DS) and nonsynonymous (DN) distances were estimated using PAML 2.0k (Yang 1999
). From each of these distance matrices, the difference in rate of molecular evolution was calculated for each comparison using the relative rates test, which uses triplets of taxa to compare the distance between two ingroup taxa and an outgroup (Sarich and Wilson 1973
; Wu and Li 1985
). When used as a clock test to detect significant departures from rate constancy, the relative rates test frequently has poor power (Scherer 1989
; Avise 1994
; Robinson et al. 1998
); for sequences less than 1,000 nucleotides, there is a substantial risk of failure to detect all cases of rate variation (Bromham et al. 2000
). Although weak as a clock test, relative rates tests provide a simple but effective tool for examining patterns of rate variation across taxa by allowing comparison of the amount of genetic change accumulated in two sister lineages since their last common ancestor from the pairwise genetic distances between taxa. The association between the difference in branch length for each comparison and the ratio of life-history variables was assessed using both parametric (Pearson correlation) and nonparametric tests (Spearman rank correlation).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
There are a number of complications of the generation time hypothesis that need clarification. Firstly, the number of germline DNA replications per organism generation varies between species, and it seems likely that it also scales with body size, potentially counteracting any generation time effect. Secondly, although mitochondria must replicate at cell division, they may also replicate throughout the lifetime of the cell at a frequency partly determined by the role and activity of the cell. This suggests that the generation time effect ought to be less obvious for mitochondrial genes, and yet they show the relationship as consistently as nuclear genes. Thirdly, the generation time effect predicts a relationship with synonymous rate (as a reflection of increase in number of mutations per unit time; Li, Tanimura, and Sharp 1987
; Ohta 1993
; Eyre-Walker and Gaut 1997
), and yet the body size effect is also a feature of total genetic distance and nonsynonymous rate (this study; Martin and Palumbi 1993
; Mooers and Harvey 1994
). Finally, generation time clearly does not provide a global explanation of lineage-specific rate variationfor example, the range of generation times in sharks overlaps with that of primates; yet shark rates are estimated to be an order of magnitude lower than primate rates (Martin and Palumbi 1993
; Martin 1995, 1999
).
Given these factors which speak against a simple and universal generation time effect, why is a relationship between body size-generation time and lineage-specific rates of molecular evolution so consistently indicated for vertebrates? The first and second complications, which concern the potential decoupling of DNA replication frequency from organismal generation time, might dampen but not negate the generation time effect. For example, the influence of mitochondrial divisions might be minimized if mitochondria are effectively quiescent in female germ cells (Allen 1995
). There are several possible explanations for the observation of a generation time effect for nonsynonymous substitutions (and total genetic distance) when the predicted effect is for synonymous change. Nearly neutral theory predicts that nonsynonymous substitutions will be approximately even over time because small animals tend to have both shorter generation time and larger population sizes, so that the effect of population size on substitution rates (faster substitution in smaller populations) could be cancelled out by the generation time effect (faster rates in small animals; Chao and Carr 1993
; Ohta 1993
). However, the magnitude of the effect of population and generation time on substitution rates is not known for real populations, and so it may be that the effect of generation time is not cancelled out by population size, leaving residual lineage-specific rate variation in the nonsynonymous substitution rate. The amino acid replacement rate might also be raised by selection, if rapid generation turnover speeds the process of selection (but see Rosenheim and Tabashnik 1993
).
The fourth complication, that generation time does not provide a global explanation for lineage-specific rates in vertebrate molecular evolution, suggests that the effect of generation time on substitution rate is overridden by other effects between widely divergent lineages, such as endothermy versus ectothermy or variation in genome size. In this way, the body size-generation time effect could be considered primarily a feature of the local molecular clock. An alternative explanation for the persistence of a generation time effect, despite complications which might reduce its influence, is that generation time may be a proxy variable for an unknown causal factor that scales with body size and its covarying life-history traits (Bromham, Rambaut, and Harvey 1996
). One such variable is metabolic rate. The metabolic rate hypothesis suggests that because small animals tend to have higher mass-specific metabolic rates than their larger relatives, they may have higher rates of production of potentially DNA-damaging metabolites, such as free oxygen radicals (Martin and Palumbi 1993
; Rand 1994
; Martin 1995, 1999
). This hypothesis is consistent with the observation that poikilotherms have slower rates of molecular evolution than endotherms. It is also consistent with higher rates of molecular evolution in mitochondrial genes than nuclear genes; mitochondria are the site of oxidative metabolism and hence might be expected to suffer the greatest effect of the production of DNA-damaging metabolites. Because metabolic rate and generation time both tend to scale with body size in vertebrates, as do many life-history traits, their effects on rate of sequence evolution may be conflated. However, two studies that explicitly compared these two hypotheses found significant correlations between rate of molecular evolution and both metabolic rate and generation time but found no evidence of an effect of metabolic rate beyond its covariation with generation time (Mooers and Harvey 1994
; Bromham, Rambaut, and Harvey 1996
). Unfortunately, insufficient metabolic rate data prevented a similar test for this study.
Clutch Size
This study revealed a positive relationship between clutch size and rate of molecular evolution for cytochrome b. This relationship could be considered consistent with a prediction of Britten (1986)
that high fecundity will speed rate of molecular evolution by increasing the number of times the genome is copied per generation. This hypothesis shares its basic premise with the generation time effect and male-driven evolutionincreasing the number of DNA replications per unit time adds more replication errors to the substitution rate.
The relationship between fecundity and number of DNA replications is not simple, as it varies between species and between sexes within a species. Production of more offspring will only increase the substitution rate by the mechanism described if increase in gamete production is achieved by increasing the average number of cell divisions taken to produce each gamete. For example, this may be true of mammalian spermatogenesis (gametes produced by a chain of cell divisions such that later sperm have undergone more cell divisions than earlier sperm), but it is less likely to be true for mammalian oogenesis. Higher fecundity may contribute to faster substitution rates in eusocial taxa than in their close nonsocial relatives (Schmitz and Moritz 1998
; L. Bromham and R. Leys, unpublished data). Because oogenesis is continuous in Hymenoptera, eusocial species that produce a large number of sterile workers before their reproductive offspring must have a much higher number of DNA replications per generation than a nonsocial species that produces reproductive offspring directly (L. Bromham and R. Leys, unpublished data). However, it seems likely that oogenesis in reptiles resembles the process in mammals, although the germ cell development differs among reptile groups (Hubert 1985
). The positive relationship between rate and clutch size may warrant further investigation as more sequence and life-history data become available.
Comparative Method and Molecular Rates
Studies of determinants of lineage-specific rates of molecular evolution have been limited by the perceived lack of a comparative method for testing correlates of rate variation (Slowinski and Arbogast 1999
). The method described here and in earlier studies (Mooers and Harvey 1994
; Bromham, Rambaut, and Harvey 1996
) provides a general statistical framework for testing for correlates of rate of molecular evolution, for assessing the generality of relationships across taxa, as well as for comparing hypotheses. This approach is limited by the number of taxa with available DNA sequences and life-history data and by the phylogenetic resolution needed to choose independent comparisons, but growth of sequence databases will allow this test to be extended to more taxa and to more genes.
Recently, Slowinski and Arbogast (1999)
suggested that the inverse relationship between body size and rate of molecular evolution revealed by Martin and Palumbi (1993)
was an artifact of their analysis for two reasons: (1) failure to account for phylogenetic bias, and (2) bias toward underestimation of substitution rate in larger-bodied taxa. The comparative method employed here overcomes both these problems. Selection of phylogenetically independent pairs ensures that any differences between members of a pair must have arisen since their last common ancestor, and this will have occurred independently of any other pair. Such differences provide phylogenetically independent data for use in statistical tests. The comparative relative rates approach also overcomes the problem of biases in the calculation of absolute rates from fossil-based dates of divergence, which could lead to consistent underestimation of rates for larger-bodied lineages with deeper divergences (Slowinski and Arbogast 1999
). Because each data point is the comparison between two lineages of the same age, any age bias in estimation of genetic distance must apply equally to both members of a pair. Bias in estimating distance could arise if the lineage with the faster rate saturates more rapidly, but this would tend to reduce the observed difference in branch length, dampening any observed relationship with size.
Implications for the Molecular Clock
The observation of a correlation between body size and rate of molecular evolution for several genes for a range of reptilian taxa echoes the body size-generation time effect observed for mammals and birds and suggests a general trend in vertebrate molecular evolution. Clearly, the data are currently limited, and the relationship must be tested for a greater number of taxa and genes before general conclusions can be drawn. Nevertheless, these results suggest that the inverse body size relationship should be explored further because it has two important implications: it provides a window on mechanisms of substitution, and it cautions reliance on a local molecular clock without rigorous testing.
The cause of the inverse body size relationship remains uncertain. The influence of number of DNA replications per unit time on substitution rate is supported by a number of observations and thus provides a plausible explanation for the influence of body size on vertebrate substitution rates. Although both generation time and metabolic rate may influence rate of molecular evolution, it is clear that neither provides a global explanation of lineage-specific substitution rates in vertebrates. It seems likely that many factors contribute to different rates of molecular evolution. Because body size correlates with a large suite of life-history traits, it is important that other traits that might influence rate of DNA change are also considered.
Whatever the cause, the observation of an inverse body size relationship for some taxa and genes suggests caution in assuming a local molecular clock for reptile species without rigorous testing. It should be acknowledged that testing for rate constancy is problematic because clock tests designed to detect lineage-specific rate variation generally have low power for most sequences used in phylogenetic analysis (Bromham et al. 2000
). This is critical for the accuracy of molecular date estimates as failure to detect sequences with lineage-specific rate variation can result in consistent overestimation of molecular date estimates (Bromham et al. 2000
). A life-history trend in rate of molecular evolution suggests that rates of molecular evolution must evolve over trees as life-history traits evolve, posing a dynamic picture of rate of molecular evolution, which complicates molecular phylogenetic inference (Sanderson 1997
).
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: molecular clock
molecular evolution
generation time effect
substitution rate
relative rates
comparative method
Address for correspondence and reprints: Lindell Bromham, School of Biological Sciences, University of Sussex, Brighton BN1 9QG, United Kingdom. l.d.bromham{at}sussex.ac.uk
.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Allen J., 1995 Separate sexes and the mitochondrial theory of ageing J. Theor. Biol 180:135-140[ISI]
Avise J. C., 1994 Molecular markers, natural history and evolution Chapman & Hall, New York
Bleiweiss R., 1998 Relative-rate tests and biological causes of molecular evolution in hummingbirds Mol. Biol. Evol 15:481-491
Branch B., 1998 Field guide to snakes and other reptiles of Southern Africa Ralph Curtis Books, Florida
Britten R. J., 1986 Rates of DNA sequence evolution differ between taxonomic groups Science 231:1393-1398[ISI][Medline]
Bromham L., A. Rambaut, P. H. Harvey, 1996 Determinants of rate variation in mammalian DNA sequence evolution J. Mol. Evol 43:610-621[ISI][Medline]
Bromham L., A. Rambaut, M. D. Hendy, D. Penny, 2000 The power of relative rates tests depends on the data J. Mol. Evol 50:296-301[ISI][Medline]
Bulmer M., K. Wolfe, P. M. Sharp, 1991 Synonymous nucleotide substitution rate in mammalian genes: implications for the molecular clock and the relationship of mammalian orders Proc. Natl. Acad. Sci. USA 88:5974-5978[Abstract]
Campbell B. N., 1997 Hic Sunt Serpentesmolecular phylogenetics and the Boidae (Serpentes: Booidea) Doctoral dissertation, Queen's University, Kingston, Ontario, Canada
Chang B. H.-J., L. C. Shimmin, S.-K. Shyue, D. Hewett-Emmett, W.-H. Li, 1994 Weak male-driven molecular evolution in rodents Proc. Natl. Acad. Sci. USA 91:827-831
Chao L., D. E. Carr, 1993 The molecular clock and the relationship between population size and generation time Evolution 47:688-690[ISI]
Cree A., C. H. Daugherty, J. M. Hay, 1995 Reproduction of a rare New Zealand reptile, the Tuatara Sphenodon punctatus, on rat-free and rat-inhabited islands Conserv. Biol 9:373-383[ISI]
Donnellan S. C., M. N. Hutchinson, K. M. Saint, 1999 Molecular evidence for the phylogeny of Australian gekkonid lizards Biol. J. Linn. Soc 67:97-118[ISI]
Easteal S., C. Collett, 1994 Consistent variation in amino-acid substitution rate, despite uniformity of mutation rate: protein evolution in mammals is not neutral Mol. Biol. Evol 11:643-647[Abstract]
Ehrhart L. M., 1995 A review of sea turtle reproduction Pp. 2938 in K. A. Bjorndal, ed. Biology and conservation of sea turtles. Smithsonian Institution Press, Washington
Ellegren H., A. K. Fridolfsson, 1997 Male-driven evolution of DNA sequences in birds Nature Genet 17:182-184[ISI][Medline]
Ernst C. H., 1992 Venomous reptiles of North America Smithsonian Institution Press, Washington
Ernst C. H., R. W. Barbour, 1989 Turtles of the world Smithsonian Institution Press, Washington
Eyre-Walker A., B. S. Gaut, 1997 Correlated rates of synonymous site evolution across plant genomes Mol. Biol. Evol 14:455-460[Abstract]
Felsenstein J., 1985 Phylogenies and the comparative method Am. Nat 125:1-15[ISI]
Fitch W. M., 1996 The variety of human virus evolution Mol. Phylogenet. Evol 5:247-258[ISI][Medline]
Greer A. E., 1989 The biology and evolution of Australian lizards Surrey Beatty and Sons, Chipping Norton, Australia
. 1997 The biology and evolution of Australian snakes Surrey Beatty and Sons, Chipping Norton, Australia
Harless M., H. Morlock, 1979 Turtle perspectives and research John Wiley, New York
Harlow P., R. Shine, 1992 Food habits and the reproductive biology of the pacific island boas J. Herpetol 26:60-66[ISI]
Harvey P. H., M. Pagel, 1991 The comparative method in evolutionary biology Oxford University Press, Oxford
Harvey P. H., A. Purvis, 1991 Comparative methods for explaining adaptations Nature 351:619-624[ISI][Medline]
Heatwole H., 1999 Sea snakes University of New South Wales (UNSW) Press, Sydney
Heatwole H., J. Taylor, 1987 Ecology of reptiles Surrey Beatty & Sons, Chipping Norton, Australia
Heppell S. S., 1998 Application of life-history theory and population model analysis to turtle conservation Copeia 1998:367375
Hubert J., 1985 The origin and development of oocytes Pp. 4174 in C. Gans, F. Billet, and P. F. A. Maderson, eds. Biology of the Reptilia: development A. John Wiley & Sons, New York
Kennett R., 1999 Reproduction of two species of freshwater turtle, Chelodina rugosa and Elseya dentata, from the wet-dry tropics of northern Australia J. Zool. Lond. 247:457473
Keogh J. S., R. Shine, S. Donnellan, 1998 Phylogenetic relationships of terrestrial Australo-Papuan elapid snakes (subfamily Hydrophiinae) based on cytochrome b and 16S rRNA sequences Mol. Phylogenet. Evol 10:67-81[ISI][Medline]
Kraus F., W. M. Brown, 1998 Phylogenetic relationships of colubroid snakes based on mitochondrial DNA sequences Zool. J. Linn. Soc 122:455-487[ISI]
Kraus F., D. G. Mink, W. M. Brown, 1996 Crotaline intergeneric relationships based on mitochondrial DNA sequence data Copeia 4:763773
Li W.-H., M. Gouy, P. Sharp, C. O'Huigin, Y.-W. Yang, 1990 Molecular phylogeny of Rodentia, Lagomorpha, Artiodactyla and Carnivora and molecular clocks Proc. Natl. Acad. Sci. USA 87:6703-6707[Abstract]
Li W.-H., M. Tanimura, P. M. Sharp, 1987 An evaluation of the molecular clock hypothesis using mammalian DNA sequences J. Mol. Evol 25:330-342[ISI][Medline]
Martin A. P., 1995 Metabolic rate and directional nucleotide substitution in animal mitochondrial DNA Mol. Biol. Evol 12:1124-1131[Abstract]
. 1999 Substitution rates of organelle and nuclear genes in sharks: implicating metabolic rate (again) Mol. Biol. Evol 16:996-1002[Abstract]
Martin A. P., S. R. Palumbi, 1993 Body size, metabolic rate, generation time and the molecular clock Proc. Natl. Acad. Sci. USA 90:4087-4091[Abstract]
Miller J. D., 1997 Reproduction in sea turtles Pp. 5182 in P. L. Lutz and J. A. Musick, eds. The biology of sea turtles. CRC Press, Boca Raton
Mooers A. Ø., P. H. Harvey, 1994 Metabolic rate, generation time and the rate of molecular evolution in birds Mol. Phylogenet. Evol 3:344-350[Medline]
Murphy J. C., R. W. Henderson, 1997 Tales of giant snakes: a historical natural history of anacondas and pythons Krieger Publishing Company, Florida
Ohta T., 1993 An examination of the generation time effect on molecular evolution Proc. Natl. Acad. Sci. USA 90:10676-10680[Abstract]
Omland K. E., 1997 Correlated rates of molecular and morphological evolution Evolution 51:1381-1393[ISI]
Pritchard P. C. H., P. Trebbau, 1984 Turtles of Venezuela Society for the Study of Amphibians and Reptiles, Athens, Ga.
Rambaut A. E., 1996 Se-Al sequence alignment editor University of Oxford, http://evolve.zoo.ox.ac.uk/software/Se-Al, Oxford, U.K
Rand D. M., 1994 Thermal habit, metabolic rate and the evolution of mitochondrial DNA Trends Ecol. Evol 9:125-131[ISI]
Robinson M., M. Gouy, C. Gautier, D. Mouchirod, 1998 Sensitivity of relative rates tests to taxonomic sampling Mol. Biol. Evol 15:1091-1098[Abstract]
Rosenheim J., B. Tabashnik, 1993 Generation time and evolution Nature 365:791-792
Rossman D. A., N. B. Ford, R. A. Seigel, 1996 The garter snakes: evolution and ecology University of Oklahoma Press, Norman
Saint K. M., C. C. Austin, S. C. Donnellan, M. N. Hutchinson, 1998 c-mos, a nuclear marker useful for squamate phylogenetic analysis Mol. Phylogenet. Evol 10:259-263[ISI][Medline]
Sanderson M. J., 1997 A nonparametric approach to estimating divergence times in the absence of rate constancy J. Mol. Evol 14:1218-1231
Sarich V. M., A. C. Wilson, 1973 Generation time and genomic evolution in primates Science 179:1144-1147[ISI][Medline]
Scherer S., 1989 The relative-rate test of the molecular clock hypothesis: a note of caution Mol. Biol. Evol 6:436-441
Schmitz J., R. F. A. Moritz, 1998 Sociality and the rate of rDNA sequence evolution in wasps (Vespidae) and honeybees (Apis) J. Mol. Evol 47:606-612[ISI][Medline]
Shaffer H. B., P. Meylan, M. L. McKnight, 1997 Tests of turtle phylogeny: molecular, morphological, and paleontological approaches Syst. Biol 46:235-268[ISI][Medline]
Shimmin L. C., B. H.-J. Chang, W. H. Li, 1993 Male-driven evolution of DNA sequences Nature 362:745-747[ISI][Medline]
. 1994 Contrasting rates of nucleotide substitution in the X-linked and Y-linked zinc-finger genes J. Mol. Evol 39:569-578[ISI][Medline]
Shine R., 1994 Allometric patterns in the ecology of Australian snakes Copeia 1994:851867
Shine R., J. B. Iverson, 1995 Patterns of survival growth and maturation in turtles Oikos 72:343-348[ISI]
Shine R., R. A. Seigel, 1996 A neglected life-history trait: clutch-size variance in snakes J. Zool. Lond 239:209-223
Slowinski J. B., B. S. Arbogast, 1999 Is the rate of molecular evolution inversely related to body size? Syst. Biol 48:396-399[ISI][Medline]
Swofford D. L., 1999 PAUP*: phylogenetic analysis using parsimony (*and other methods). Version 4 Sinauer Associates, Sunderland, Mass
Whitfield Gibbons J., 1990 Life history and ecology of the slider turtle Smithsonian Institution Press, Washington
Wu C.-I., W.-H. Li, 1985 Evidence for higher rates of nucleotide substitutions in rodents than in man Proc. Natl. Acad. Sci. USA 82:1741-1745[Abstract]
Wynn A. H., G. R. Zug, 1985 Observations on the reproductive biology of Candoia carinata (Serpentes, Boidae) Snake 17:15-24
Yang Z., 1999 Phylogenetic analysis by maximum likelihood (PAML) University College, http://abacus.gene.ucl.ac.uk/software/paml.html, London
Yang Z., R. Nielsen, 1998 Synonymous and nonsynonymous rate variation in nuclear genes of mammals J. Mol. Evol 46:409-418[ISI][Medline]