On the rate of DNA sequence evolution inDrosophila

Summary Analysis of the rate of nucleotide substitution at silent sites inDrosophila genes reveals three main points. First, the silent rate varies (by a factor of two) among nuclear genes; it is inversely related to the degree of codon usage bias, and so selection among synonymous codons appears to constrain the rate of silent substitution in some genes. Second, mitochondrial genes may have evolved only as fast as nuclear genes with weak codon usage bias (and two times faster than nuclear genes with high codon usage bias); this is quite different from the situation in mammals where mitochondrial genes evolve approximately 5–10 times faster than nuclear genes. Third, the absolute rate of substitution at silent sites in nuclear genes inDrosophila is about three times hihger than the average silent rate in mammals.     

Dating the early evolution of plants: detection and molecular clock analyses of orthologs

Orthologs generally are under selective pressure against loss of function, while paralogs usually accumulate mutations and finally die or deviate in terms of function or regulation. Most ortholog detection methods contaminate the resulting datasets with a substantial amount of paralogs. Therefore we aimed to implement a straightforward method that allows the detection of ortholog clusters with a reduced amount of paralogs from completely sequenced genomes. The described cross-species expansion of the reciprocal best BLAST hit method is a time-effective method for ortholog detection, which results in 68% truly orthologous clusters and the procedure specifically enriches single-copy orthologs. The detection of true orthologs can provide a phylogenetic toolkit to better understand evolutionary processes. In a study across six photosynthetic eukaryotes, nuclear genes of putative mitochondrial origin were shown to be over-represented among single copy orthologs. These orthologs are involved in fundamental biological processes like amino acid metabolism or translation. Molecular clock analyses based on this dataset yielded divergence time estimates for the red/green algae (1,142 MYA), green algae/land plant (725 MYA), mosses/seed plant (496 MYA), gymno-/angiosperm (385 MYA) and monocotyledons/core eudicotyledons (301 MYA) divergence times. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mitochondrial DNA evolution in theDrosophila nasuta subgroup of species

Summary TheDrosophila nasuta group consists of about 12 closely related species distributed throughout the Indo-Pacific region. They are of great interest because of their evolutionary idiosyncrasies including little morphological differentiation, the ability to intercross in the laboratory often producing fertile offspring, and substantial chromosomal evolution. Studies of metric traits, reproductive isolation, and chromosomal and enzyme polymorphisms have failed to resolve the phylogeny of the species. We report the results of a survey of the mitochondrial DNA (mtDNA) restriction patterns of the species. The phylogeny obtained is consistent with other available information and suggests thatD. albomicans may represent the ancestral lineage of the group. The amount of polymorphism in local populations (=1.0% per site) is within the typical range observed in other animals, includingDrosophila. The degree of differentiation between species is, however, low: the origin of the group is tentatively dated about 6–8 million years ago. This study confirms the usefulness of mtDNA restriction patterns for ascertaining the phylogeny of closely related species.     

Lack of correlation between body mass and metabolic rate in Drosophila melanogaster

We examined the association between body mass and metabolic rate in Drosophila melanogaster under a variety of conditions. These included comparisons of body mass and metabolic rate in flies from different laboratory lines measured at different ages, over different metabolic sampling periods, and comparisons using wet versus dry mass data. In addition, the relationship between body mass and metabolic rate was determined for flies recently collected from wild populations. In no case was there a significant correlation between body mass and metabolic rate. These results indicate that care must be taken when attempting to account for the effects of body mass on metabolic rate. Expressing such data in mass-specific units may be an inappropriate method of attempting to control for the effects of differences in body mass.     

Network models for sequence evolution

We introduce a general class of models for sequence evolution that includes network phylogenies. Networks, a generalization of strictly tree-like phylogenies, are proposed to model situations where multiple lineages contribute to the observed sequences. An algorithm to compute the probability distribution of binary character-state configurations is presented and statistical inference for this model is developed in a likelihood framework. A stepwise procedure based on likelihood ratios is used to explore the space of models. Starting with a star phylogeny, new splits (nontrivial bipartitions of the sequence set) are successively added to the model until no significant change in the likelihood is observed. A novel feature of our approach is that the new splits are not necessarily constrained to be consistent with a treelike mode of evolution. The fraction of invariable sites is estimated by maximum likelihood simultaneously with other model parameters and is essential to obtain a good fit to the data. The effect of finite sequence length on the inference methods is discussed. Finally, we provide an illustrative example using aligned VPl genes from the foot and mouth disease viruses (FMDV). The different serotypes of the FMDV exhibit a range of treelike and network evolutionary relationships.Correspondence to: A. von Haeseler     

Mitochondrial DNA evolution in primates: Transition rate has been extremely low in the lemur

Summary Based on mitochondrial DNA (mt-DNA) sequence data from a wide range of primate species, branching order in the evolution of primates was inferred by the maximum likelihood method of Felsenstein without assuming rate constancy among lineages. Bootstrap probabilities for being the maximum likelihood tree topology among alternatives were estimated without performing a maximum likelihood estimation for each resampled data set. Variation in the evolutionary rate among lineages was examined for the maximum likelihood tree by a method developed by Kishino and Hasegawa. From these analyses it appears that the transition rate of mtDNA evolution in the lemur has been extremely low, only about 1/10 that in other primate lines, whereas the transversion rate does not differ significantly from that of other primates. Furthermore, the transition rate in catarrhines, except the gibbon, is higher than those in the tarsier and in platyrrhines, and the transition rate in the gibbon is lower than those in other catarrhines. Branching dates in primate evolution were estimated by a molecular clock analysis of mtDNA, taking into account the rate of variation among different lines, and the results were compared with those estimated from nuclear DNA. Under the most likely model, where the evolutionary rate of mtDNA has been unifrom within a great apes/human calde, human/chimpanzee clustering is preferred to the alternative branching orders among human, chimpanzee, and gorilla.     

Partitioning the variation in mammalian substitution rates

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.     

Mitochondrial DNA variation in human metabolic rate and energy expenditure

The role of climate in driving selection of mtDNA as Homo sapiens migrated out of Africa into Eurasia remains controversial. We evaluated the role of mtDNA variation in resting metabolic rate (RMR) and total energy expenditure (TEE) among 294 older, community-dwelling African and European American adults from the Health, Aging and Body Composition Study. Common African haplogroups L0, L2 and L3 had significantly lower RMRs than European haplogroups H, JT and UK with haplogroup L1 RMR being intermediate to these groups. This study links mitochondrial haplogroups with ancestry-associated differences in metabolic rate and energy expenditure.     

Nutritional ecology of dimorphic herbivores: digestion and passage rates in Nubian ibex

We compared forage digestion and passage rates among three groups of Nubian ibex (Capra ibex nubiana) — mature males, non-lactating females, and lactating females — to test hypotheses relating intraspecific digestive ability to body mass and reproduction costs. We hypothesized that large males (60 kg) would exhibit longer forage retention times and more complete digestion of fermentable cell walls than adult females (23 kg). We tested these predictions by measuring digestion and retention of a grass hay and an alfalfa hay, forages that exhibited contrasting rates and extents of cell wall digestion. Consistent with predictions, males retained both forages longer than non-lactating females. However, by substantially increasing gut fill, lactating females increased both intake and retention time with respect to non-lactating females. Contrary to predictions, all three groups digested the grass (66% digestible) and alfalfa hay (63%) equally well. Alfalfa cell wall was less digestible than that of grass hay (60% vs 69% digestible), and retention time of alfalfa was consistently, but not statistically significantly, shorter. Fiber digestion was not correlated with retention time, emphasizing the ability of behavioral processes to modify digestion rate. We postulate that females achieved their greater digestion rate by masticating forages much more thoroughly than males.     

The origin of allometric scaling laws in biology

The empirical rules relating metabolic rate and body size are described in terms of (i) a scaling exponent, which refers to the ratio of the fractional change in metabolic rate to a change in body size, (ii) a proportionality constant, which describes the rate of energy expenditure in an organism of unit mass. This article integrates the chemiosmotic theory of energy transduction with the methods of quantum statistics to propose a molecular mechanism which, in sharp contrast to competing models, explains both the variation in scaling exponents and the taxon-specific differences in proportionality constants. The new model is universal in the sense that it applies to unicellular organisms, plants and animals.     

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.     

Primate brains,the ‘island rule’ and the evolution of Homo floresiensis

The taxonomic status of the small bodied hominin, Homo floresiensis, remains controversial. One contentious aspect of the debate concerns the small brain size estimated for specimen LB1 (Liang Bua 1). Based on intraspecific mammalian allometric relationships between brain and body size, it has been argued that the brain of LB1 is too small for its body mass and is therefore likely to be pathological. The relevance and general applicability of these scaling rules has, however, been challenged, and it is not known whether highly encephalized primates adapt to insular habitats in a consistent manner. Here, an analysis of brain and body size evolution in seven extant insular primates reveals that although insular primates follow the ‘island rule’, having consistently reduced body masses compared with their mainland relatives, neither brain mass nor relative brain size follow similar patterns, contrary to expectations that energetic constraints will favour decreased relative brain size. Brain:body scaling relationships previously used to assess the plausibility of dwarfism in H. floresiensis tend to underestimate body masses of insular primates. In contrast, under a number of phylogenetic scenarios, the evolution of brain and body mass in H. floresiensis is consistent with patterns observed in other insular primates.     

Molecular evolution and phylogeny of the human AIDS viruses LAV,HTLV-III,and ARV

Summary A phylogenetic tree for the human lymphadenopathy-associated virus (LAV), the human T-cell lymphotrophic virus type III (HTLV-III), and the acquired immune deficiency syndrome (AIDS)-associated retrovirus (ARV) has been constructed from comparisons of the amino acid sequences of their gag proteins. A method is proposed for estimating the divergence times among these AIDS viruses and the rates of nucleotide substitution for their RNA genomes. The analysis indicates that the LAV and HTLV-III strains diverged from one another after 1977 and that their common ancestor diverged from the ARV virus no more than 10 years earlier. Hence, the evolutionary diversity among strains of the AIDS viruses apparently has been generated within the last 20 years. It is estimated that the genome of the AIDS virus has a nucleotide substitution rate on the order of 10^–3 per site per year, with the rate in the second half of the genome being double that in the first half.     

Utility of nuclear DNA intron markers at lower taxonomic levels: phylogenetic resolution among nine Tragelaphus spp

Phylogenetic relationships among the nine spiral-horn antelope species of the African bovid tribe Tragelaphini are controversial. In particular, mitochondrial DNA sequencing studies are not congruent with previous morphological investigations. To test the utility of nuclear DNA intron markers at lower taxonomic levels and to provide additional data pertinent to tragelaphid evolution, we sequenced four nuclear DNA segments (MGF, PRKCI, SPTBN, and THY) and combined these data with mitochondrial DNA sequences from three genes (cytochrome b, 12S rRNA, and 16S rRNA). Our molecular supermatrix comprised 4682 characters which were analyzed independently and in combination. Parsimony and model based phylogenetic analyses of the combined nuclear DNA data are congruent with those derived from the analysis of mitochondrial gene sequences. The corroboration between nuclear and mtDNA gene trees reject the possibility that genetic processes such as lineage sorting, gene duplication/deletion and hybrid speciation account for the conflict evident in the previously published phylogenies. It suggests rather that the morphological characters used to delimit the Tragelaphid species are subject to convergent evolution. Divergence times among species, calculated using a relaxed Bayesian molecular clock, are consistent with hypotheses proposing that climatic oscillations and their impact on habitats were the major forces driving speciation in the tribe Tragelaphini.     

Tempo and mode of sequence evolution in mitochondrial DNA of HawaiianDrosophila

Summary Sequence comparisons were made for up to 667 bp of DNA cloned from 14 kinds of HawaiianDrosophila and five other dipteran species. These sequences include parts of the genes for NADH dehydrogenase (subunits 1, 2, and 5) and rRNA (from the large ribosomal subunit). Because the times of divergence among these species are known approximately, the sequence comparisons give insight into the evolutionary dynamics of this molecule. Transitions account for nearly all of the differences between sequences that have diverged by less than 2%; for these sequences the mean rate of divergence appears to be about 2%/Myr. In comparisons involving greater divergence times and greater sequence divergence, relatively more of the sequence differences are due to transversions. Specifically, the fraction of these differences that are counted as transversions rises from an initial value of less than 0.1 to a plateau value of nearly 0.6. The time required to reach half of the plateau value, about 10 Myr, is similar to that for mammalian mtDNA. The mtDNAs of flies and mammals are also alike in the shape of the curve relating the percentage of positions at which there are differences in protein-coding regions to the time of divergence. For both groups of animals, the curve has a steep initial slope ascribable to fast accumulation of synonymous substitutions and a shallow final slope resulting from the slow accumulation of substitutions causing amino acid replacements. However, the percentage of all sites that can experience a high rate of substitution appears to be only about 8% for fly mtDNA compared to about 20% for mammalian mtDNA. The low percentage of hypervariable sites may be a consequence of a functional constraint associated with the low content of guanine and cytosine in fly mtDNA.     

Interhomologue sequence variation of alpha satellite DNA from human chromosome 17: Evidence for concerted evolution along haplotypic lineages

Alpha satellite DNA is a family of tandemly repeated DNA found at the centromeres of all primate chromosomes. Different human chromosomes 17 in the population are characterized by distinct alpha satellite haplotypes, distinguished by the presence of variant repeat forms that have precise monomeric deletions. Pairwise comparisons of sequence diversity between variant repeat units from each haplotype show that they are closely related in sequence. Direct sequencing of PCR-amplified alpha satellite reveals heterogeneous positions between the repeat units on a chromosome as two bands at the same position on a sequencing ladder. No variation was detected in the sequence and location of these heterogeneous positions between chromosomes 17 from the same haplotype, but distinct patterns of variation were detected between chromosomes from different haplotypes. Subsequent sequence analysis of individual repeats from each haplotype confirmed the presence of extensive haplotype-specific sequence variation. Phylogenetic inference yielded a tree that suggests these chromosome 17 repeat units evolve principally along haplotypic lineages. These studies allow insight into the relative rates and/or timing of genetic turnover processes that lead to the homogenization of tandem DNA families. Correspondence to: H.F. Willard     

Restriction-site variation of PCR-amplified chloroplast DNA regions and its implication for the evolution and taxonomy of Actinidia

 Twenty six restriction sites from five PCR-amplified chloroplast DNA sequences (rbcL, psbA, rpoB, and two spacers flanking the trnL gene) were mapped and analysed in 20 Actinidia taxa, encompassing all four sections into which the genus is divided. At least three species out of the 20 examined have been found to have originated through natural interspecific hybridisation on the basis of the discrepancy between morphological and biochemical traits and the cpDNA profiles of pairs of species. A widely reticulate evolution has therefore been postulated in Actinidia. Wagner and weighted parsimony analysis produced consensus trees that did not match the traditional taxonomy based on morphological characters. The molecular data clearly showed that some taxa, such as A. rufa and A. kolomikta, occupy a wrong position and most, if not all, of the traditional groups represented by sections and series are weakly supported, since they appear as polyphyletic. A. chinensis and A. deliciosa were confirmed to be very closely related. Since chloroplast DNA is paternally inherited in Actinidia, A. chinensis is a paternal progenitor, if not the only one, of A. deliciosa, the domesticated kiwifruit. Received: 18 August 1997 / Accepted: 6 October 1997     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.     

Phenological resonance and quantum life history

The principle of 'quantum life history' is proposed here as a complementary viewpoint to current modeling of body size and life history evolution which usually considers a 'fast-slow continuum' of covarying life history traits. This principle emphasizes the discrete (and primary) nature of development time caused by the effect of phenological resonance (the compliance of development time with periodicities of earth rotation). The body mass, in turn, complies with development time, which generates body mass attractors. This principle is illustrated with mammals as exemplary group. The adaptive radiation of Cenozoic mammals is supposed to proceed as a competition-driven diversification of body sizes and development times around the strongest (year-long) resonant mode of development time corresponding to body mass of about 1 kg. Mammals with this body mass are shown here to have a largest genome size and a lowest (body mass-corrected) basal metabolic rate. This extends the previously reported negative relation between genome size and metabolic rate to the realm of nonlinearity, and suggests that selection against the accumulation of non-coding DNA in the genome is relaxed in mammals with this body mass.