Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996     

Synonymous and Nonsynonymous Substitutions in Genes from Gramineae: Intragenic Correlations

In this work, we have investigated the relationships between synonymous and nonsynonymous rates and base composition in coding sequences from Gramineae to analyze the factors underlying the variation in substitutional rates. We have shown that in these genes the rates of nucleotide divergence, both synonymous and nonsynonymous, are, to some extent, dependent on each other and on the base composition. In the first place, the variation in nonsynonymous rate is related to the GC level at the second codon position (the higher the GC₂ level, the higher the amino acid replacement rate). The correlation is especially strong with T₂, the coefficients being significant in the three data sets analyzed. This correlation between nonsynonymous rate and base composition at the second codon position is also detectable at the intragenic level, which implies that the factors that tend to increase the intergenic variance in nonsynonymous rates also affect the intragenic variance. On the other hand, we have shown that the synonymous rate is strongly correlated with the GC₃ level. This correlation is observed both across genes and at the intragenic level. Similarly, the nonsynonymous rate is also affected at the intragenic level by GC₃ level, like the silent rate. In fact, synonymous and nonsynonymous rates exhibit a parallel behavior in relation to GC₃ level, indicating that the intragenic patterns of both silent and amino acid divergence rates are influenced in a similar way by the intragenic variation of GC₃. This result, taken together with the fact that the number of genes displaying intragenic correlation coefficients between synonymous and nonsynonymous rates is not very high, but higher than random expectation (in the three data sets analyzed), strongly suggests that the processes of silent and amino acid replacement divergence are, at least in part, driven by common evolutionary forces in genes from Gramineae. Received: 2 July 1998 / Accepted: 18 April 1999     

Sensitivity of Patterns of Molecular Evolution to Alterations in Methodology: A Critique of Hughes and Yeager

Employing a set of 43 othologous mouse and rat genes, Hughes and Yeager (J. Mol. Evol. 45:125–130, 1997) reported (1) no correlation between synonymous and nonsynonymous rates of nucleotide substitution, (2) a positive correlation between intronic GC contents (GC _i) and intronic substitution rates (K _i), (3) that the average K _i value was very similar to the average K _s value, and (4) that the compositional correlation between the rat and the mouse genes is stronger at the third codon position (GC₃) than at the first and second codon positions (GC₁₂). We have examined the robustness of these results to alterations in substitution rate estimation protocol, alignment protocol, and statistical procedure. We find that a significant correlation between K _a and K _s is observed either if a rank correlation statistic is used instead of regression analysis, if one outlier is excluded from the analysis, or if a regression weighted by gene size is employed. The correlation between K _i and GC _i we find to be sensitive to changes in alignment protocol and disappears on the use of weighted means. The finding that K _s and K _i are approximately the same is dependent on the method for estimating K _s values. Finally, the variance around the regression line of rat GC₃ versus mouse GC₃ we find to be significantly higher than that in GC₁₂. The source of the discrepancy between this and Hughes and Yeager's result is unclear. The variance around the line for GC₄ is higher still, as might be expected. Using a methodology that may be considered preferable to that of Hughes and Yeager, we find that all four of their results are contradicted. More importantly this analysis reinforces the need for caution in assembling and analyzing data sets, as the degree of sensitivity to what many might consider minor methodological alterations is unexpected. Received: 2 February 1998 / Accepted: 23 March 1998     

Structural and Evolutionary Relationships Among Chitinases of Flowering Plants

The analysis of nuclear-encoded chitinase sequences from various angiosperms has allowed the categorization of the chitinases into discrete classes. Nucleotide sequences of their catalytic domains were compared in this study to investigate the evolutionary relationships between chitinase classes. The functionally distinct class III chitinases appear to be more closely related to fungal enzymes involved in morphogenesis than to other plant chitinases. The ordering of other plant chitinases into additional classes mainly relied on the presence of auxiliary domains—namely, a chitin-binding domain and a carboxy-terminal extension—flanking the main catalytic domain. The results of our phylogenetic analyses showed that classes I and IV form discrete and well-supported monophyletic groups derived from a common ancestral sequence that predates the divergence of dicots and monocots. In contrast, other sequences included in classes I* and II, lacking one or both types of auxiliary domains, were nested within class I sequences, indicating that they have a polyphyletic origin. According to phylogenetic analyses and the calculation of evolutionary rates, these chitinases probably arose from different class I lineages by relatively recent deletion events. The occurrence of such evolutionary trends in cultivated plants and their potential involvement in host–pathogen interactions are discussed. Received: 5 July 1996 / Accepted: 9 January 1997     

Molecular evolution of the Metazoan protein kinase C multigene family

Protein kinases C (PKCs) comprise closely related Ser/Thr kinases, ubiquitously present in animal tissues; they respond to second messengers, e.g., Ca²⁺ and/or diacylglycerol, to express their activities. Two PKCs have been sequenced from Geodia cydonium, a member of the lowest multicellular animals, the sponges (Porifera). One sponge G. cydonium PKC, GCPKC1, belongs to the ``novel' (Ca<sup>2+</sup>-independent) PKC (nPKC) subfamily while the second one, GCPKC2, has the hallmarks of the ``conventional' (Ca²⁺-dependent) PKC (cPKC) subfamily. The alignment of the Ser/Thr catalytic kinase domains, of the predicted aa sequences for these cDNAs with respective segments from previously reported sequences, revealed highest homology to PKCs from animals but also distant relationships to Ser/Thr kinases from protozoa, plants, and bacteria. However, a comparison of the complete structures of the sponge PKCs, which are—already—identical to those of nPKCs and cPKCs from higher metazoa, with the structures of protozoan, plant, and bacterial Ser/Thr kinases indicates that the metazoan PKCs have to be distinguished from the nonmetazoan enzymes. These data indicate that metazoan PKCs have a universal common ancestor which they share with the nonmetazoan Ser/Thr kinases with respect to the kinase domain, but they differ from them in overall structural composition. Received: 10 January 1996 / Accepted: 12 March 1996     

Rates of Molecular Evolution in RNA Viruses: A Quantitative Phylogenetic Analysis

The study of rates of nucleotide substitution in RNA viruses is central to our understanding of their evolution. Herein we report a comprehensive analysis of substitution rates in 50 RNA viruses using a recently developed maximum likelihood phylogenetic method. This analysis revealed a significant relationship between genetic divergence and isolation time for an extensive array of RNA viruses, although more rate variation was usually present among lineages than would be expected under the constraints of a molecular clock. Despite the lack of a molecular clock, the range of statistically significant variation in overall substitution rates was surprisingly narrow for those viruses where a significant relationship between genetic divergence and time was found, as was the case when synonymous sites were considered alone, where the molecular clock was rejected less frequently. An analysis of the ecological and genetic factors that might explain this rate variation revealed some evidence of significantly lower substitution rates in vector-borne viruses, as well as a weak correlation between rate and genome length. Finally, a simulation study revealed that our maximum likelihood estimates of substitution rates are valid, even if the molecular clock is rejected, provided that sufficiently large data sets are analyzed. Received: 23 February 2001 / Accepted: 3 July 2001     

The compositional evolution of the murid genome

Murid rodents show much less variation in isochore base composition than do most other mammals, a difference which has been referred to as the murid shift. We have investigated the murid shift by asking (1) whether the murid shift is ongoing and (2) whether there is any evidence of selection or biased gene conversion affecting base composition in the present-day mouse genome. By estimating the ancestral base composition of protein-coding genes in murids we can confirm that the murid shift is ongoing. Tests using nongenic polymorphism data fail to reject the hypothesis that base composition is due to mutation bias alone. However, the patterns of compositional change suggested by the polymorphism and divergence data differ, suggesting the possibility of two murid shifts.     

Molecular Clock Calibrations and Metazoan Divergence Dates

It has recently been argued that living metazoans diverged over 800 million years ago, based on evidence from 22 nuclear genes for such a deep divergence between vertebrates and arthropods (Gu 1998). Two ``internal' calibration points were used. However, only one fossil divergence date (the mammal–bird split) was directly used to calibrate the molecular clock. The second calibration point (the primate–rodent split) was based on molecular estimates that were ultimately also calibrated by the same mammal–bird split. However, the first tetrapods that can be assigned with confidence to either the mammal (synapsid) lineage or the bird (diapsid) lineage are approximately 288 million years old, while the first mammals that can be assigned with confidence to either the primate or the rodent lineages are 65 million years old, or 85 million years old if ferungulates are part of the primate lineage and zhelestids are accepted as ferungulate relatives. Recalibration of the protein data using these fossil dates indicates that metazoans diverged between 791 and 528 million years ago, a result broadly consistent with the palaeontological documentation of the ``Cambrian explosion.' The third, ``external' calibration point (the metazoan–fungal divergence) was similarly problematic, since it was based on a controversial molecular study (which in turn used fossil dates including the mammal–bird split); direct use of fossils for this calibration point gives the absurd dating of 455 million years for metazoan divergences. Similar calibration problems affect another recent study (Wang et al. 1999), which proposes divergences for metazoans of 1000 million years or more: recalibrations of their clock again yields much more recent dates, some consistent with a ``Cambrian explosion' scenario. Molecular clock studies have persuasively argued for the imperfection of the fossil record but have rarely acknowledged that their inferences are also directly based on this same record. Received: 26 January 1999 / Accepted: 14 April 1999     

The Power of Relative Rates Tests Depends on the Data

One of the most useful features of molecular phylogenetic analyses is the potential for estimating dates of divergence of evolutionary lineages from the DNA of extant species. But lineage-specific variation in rate of molecular evolution complicates molecular dating, because a calibration rate estimated from one lineage may not be an accurate representation of the rate in other lineages. Many molecular dating studies use a ``clock test' to identify and exclude sequences that vary in rate between lineages. However, these clock tests should not be relied upon without a critical examination of their effectiveness at removing rate variable sequences from any given data set, particularly with regard to the sequence length and number of variable sites. As an illustration of this problem we present a power test of a frequently employed triplet relative rates test. We conclude that (1) relative rates tests are unlikely to detect moderate levels of lineage-specific rate variation (where one lineage has a rate of molecular evolution 1.5 to 4.0 times the other) for most commonly used sequences in molecular dating analyses, and (2) this lack of power is likely to result in substantial error in the estimation of dates of divergence. As an example, we show that the well-studied rate difference between murid rodents and great apes will not be detected for many of the sequences used to date the divergence between these two lineages and that this failure to detect rate variation is likely to result in consistent overestimation the date of the rodent–primate split. Received: 9 June 1999 / Accepted: 22 October 1999     

Molecular Evolutionary Analysis of the Thiamine-Diphosphate-Dependent Enzyme,Transketolase

Members of the transketolase group of thiamine-diphosphate-dependent enzymes from 17 different organisms including mammals, yeast, bacteria, and plants have been used for phylogenetic reconstruction. Alignment of the amino acid and DNA sequences for 21 transketolase enzymes and one putative transketolase reveals a number of highly conserved regions and invariant residues that are of predicted importance for enzyme activity, based on the crystal structure of yeast transketolase. One particular sequence of 36 residues has some similarities to the nucleotide-binding motif and we designate it as the transketolase motif. We report further evidence that the recP protein from Streptococcus pneumoniae might be a transketolase and we list a number of invariant residues which might be involved in substrate binding. Phylogenies derived from the nucleotide and the amino acid sequences by various methods show a conventional clustering for mammalian, plant, and gram-negative bacterial transketolases. The branching order of the gram-positive bacteria could not be inferred reliably. The formaldehyde transketolase (sometimes known as dihydroxyacetone synthase) of the yeast Hansenula polymorpha appears to be orthologous to the mammalian enzymes but paralogous to the other yeast transketolases. The occurrence of more than one transketolase gene in some organisms is consistent with several gene duplications. The high degree of similarity in functionally important residues and the fact that the same kinetic mechanism is applicable to all characterized transketolase enzymes is consistent with the proposition that they are all derived from one common ancestral gene. Transketolase appears to be an ancient enzyme that has evolved slowly and might serve as a model for a molecular clock, at least within the mammalian clade. Received: 13 September 1995 / Accepted: 14 November 1996     

Implications of Fluctuations in Substitution Rates: Impact on the Uncertainty of Branch Lengths and on Relative-Rate Tests

Many tests of the lineage dependence of substitution rates, computations of the error of evolutionary distances, and simulations of molecular evolution assume that the rate of evolution is constant in time within each lineage descended from a common ancestor. However, estimates of the index of dispersion of numbers of mammalian substitutions suggest that the rate has time-dependent variations consistent with a fractal-Gaussian-rate Poisson process, which assumes common descent without assuming rate constancy. While this model does not affect certain relative-rate tests, it substantially increases the uncertainty of branch lengths. Thus, fluctuations in the rate of substitution cannot be neglected in calculations that rely on evolutionary distances, such as the confidence intervals of divergence times and certain phylogenetic reconstructions. The fractal-Gaussian-rate Poisson process is compared and contrasted with previous models of molecular evolution, including other Poisson processes, the fractal renewal process, a Lévy-stable process, a fractional-difference process, and a log-Brownian process. The fractal models are more compatible with mammalian data than the nonfractal models considered, and they may also be better supported by Darwinian theory. Although the fractal-Gaussian-rate Poisson process has not been proven to have better agreement with data or theory than the other fractal models, its Gaussian nature simplifies the exploration of its impact on evolutionary distance errors and relative-rate tests. Received: 29 September 1999 / Accepted: 20 January 2000     

The Effects of Social Structure, Geographical Structure, and Population Size on the Evolution of Mitochondrial DNA: II. Molecular Clocks and the Lineage Sorting Period

Evolutionary geneticists have increasingly used sequence variation in mitochondrial DNA (mtDNA) as a source of historical information. However, conclusions based on these data remain tentative because a sufficiently clear understanding of the evolutionary dynamics of mtDNA has yet to be developed. In this paper we present the results of computer simulations designed to illustrate the effects of social structure, geographical structure, and population size on the rate of nucleotide substitution and lineage sorting of mtDNA. The model is based in part on the social structure of macaque monkeys. Simulated populations of females were divided into 25 social groups; the animals in each were distributed in a hierarchy of four dominance rank categories. The probabilities for offspring survivorship were varied among dominance ranks to reflect the fitness consequences of social structure. Population size was varied across runs from 100 to 300 females. The pattern of female migration was also varied to mimic either the island model or the stepping-stone model. All these variables are shown to affect the lineage sorting period (LSP), and certain combinations of parameter values can cause the retention of mtDNA polymorphisms for a very long time. In addition, the simulations exhibited a negative relationship between the LSP and substitution rate over a modest and realistic range of LSP values. An important implication of these results is that estimates of time since isolation based on the assumption of a constant molecular clock may be biased and unreliable.     

Phylogeny of Organisms Investigated by the Base-Pair Changes in the Stem Regions of Small and Large Ribosomal Subunit RNAs

In order to obtain the evolutionary distance data that are as purely additive as possible, we have developed a novel method for evaluating the evolutionary distances from the base-pair changes in stem regions of ribosomal RNAs (rRNAs). The application of this method to small-subunit (SSU) and large-subunit (LSU) rRNAs provides the distance data, with which both the unweighted pair group method of analysis and the neighbor-joining method give almost the same tree topology of most organisms except for some Protoctista, thermophilic bacteria, parasitic organisms, and endosymbionts. Although the evolutionary distances calculated with LSU rRNAs are somewhat longer than those with SSU rRNAs, the difference, probably due to a slight difference in functional constraint, is substantially decreased when the distances are converted into the divergence times of organisms by the measure of the time scale estimated in each type of rRNAs. The divergence times of main branches agree fairly well with the geological record of organisms, at least after the appearance of oxygen-releasing photosynthesis, although the divergence times of Eukaryota, Archaebacteria, and Eubacteria are somewhat overestimated in comparison with the geological record of Earth formation. This result is explained by considering that the mutation rate is determined by the accumulation of misrepairs for DNA damage caused by radiation and that the effect of radiation had been stronger before the oxygen molecules became abundant in the atmosphere of the Earth. Received: 23 October 1997 / Accepted: 12 August 1998     

The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Phylogenetics of Perissodactyla and Tests of the Molecular Clock

Two mitochondrial genes, the protein-coding cytochrome c oxidase subunit II (COII) gene and a portion of the 12S rRNA gene, were used for phylogenetic investigation of the mammalian order Perissodactyla. The primary objective of the study was to utilize the extensive fossil record of perissodactyls for calibrating molecular clocks and comparing estimates of divergence times using both genes and two fossil calibration points. Secondary objectives included clarification of previously unresolved relationships within Tapiridae and comparison of the results of separate and combined analyses of two genes. Analyses included several perissodactyl lineages representing all three families (Tapiridae, Equidae, and Rhinocerotidae), most extant genera, all four species of tapirs, two to four species of rhinoceros, and two species of Equus. The application of a relatively recent fossil calibration point and a relatively ancient calibration point produced greatly different estimates of evolutionary rates and divergence times for both genes, even though a relative rates test did not find significant rate differences among taxa. A likelihood-ratio test, however, rejected a molecular clock for both genes. Neither calibration point produced estimates of divergence times consistent with paleontological evidence over a range of perissodactyl radiations. The combined analysis of both genes produces a well-resolved phylogeny with Perissodactyla that conforms to traditional views of interfamilial relationships and supports monophyly of neotropical tapirs. Combining the data sets increases support for most nodes but decreases the support for a neotropical tapir clade because the COII and 12S rRNA data sets are in conflict for tapir relationships. Received: 6 January 1999 / Accepted: 2 August 1999     

Molecular Characterization and Evolution of the Amylase Multigene Family of Drosophila ananassae

Drosophila ananassae is known to produce numerous alpha-amylase variants. We have cloned seven different Amy genes in an African strain homozygous for the AMY1,2,3,4 electrophoretic pattern. These genes are organized as two main clusters: the first one contains three intronless copies on the 2L chromosome arm, two of which are tandemly arranged. The other cluster, on the 3L arm, contains two intron-bearing copies. The amylase variants AMY1 and AMY2 have been assigned to the intronless cluster, and AMY3 and AMY4 to the second one. The divergence of coding sequences between clusters is moderate (6.1% in amino acids), but the flanking regions are very different, which could explain their differential regulation. Within each cluster, coding and noncoding regions are conserved. Two very divergent genes were also cloned, both on chromosome 3L, but very distant from each other and from the other genes. One is the Amyrel homologous (41% divergent), the second one, Amyc1 (21.6% divergent) is unknown outside the D. ananassae subgroup. These two genes have unknown functions. Received: 30 May 2000 / Accepted: 17 July 2000     

Curiously Modern DNA for a ``250 Million-Year-Old' Bacterium

Studies of ancient DNA have attracted considerable attention in scientific journals and the popular press. Several of the more extreme claims for ancient DNA have been questioned on biochemical grounds (i.e., DNA surviving longer than expected) and evolutionary grounds (i.e., nucleotide substitution patterns not matching theoretical expectations for ancient DNA). A recent letter to Nature from Vreeland et al. (2000), however, tops all others with respect to age and condition of the specimen. These researchers extracted and cultured a bacterium from an inclusion body from what they claim is a 250 million-year (Myr)-old salt crystal. If substantiated, this observation could fundamentally alter views about bacterial physiology, ecology and evolution. Here we report on molecular evolutionary analyses of the 16S rDNA from this specimen. We find that 2-9-3 differs from a modern halophile, Salibacillus marismortui, by just 3 unambiguous bp in 16S rDNA, versus the ∼59 bp that would be expected if these bacteria evolved at the same rate as other bacteria. We show, using a Poisson distribution, that unless it can be shown that S. marismortui evolves 5 to 10 times more slowly than other bacteria for which 16S rDNA substitution rates have been established, Vreeland et al.'s claim would be rejected at the 0.05 level. Also, a molecular clock test and a relative rates test fail to substantiate Vreeland et al.'s claim that strain 2-9-3 is a 250-Myr-old bacterium. The report of Vreeland et al. thus falls into a long series of suspect ancient DNA studies. Received: 12 April 2001 / Accepted: 9 June 2001     

Molecular evolution of hemoglobins of antarctic fishes (Notothenioidei)

Amino acid sequences of α- and β-chains of human hemoglobin and of hemoglobins of coelacanth and 24 teleost fish species, including 11 antarctic and two temperate Notothenioidei, were analyzed using maximum parsimony. Trees were derived for the α- and β-chains separately and for tandemly arranged sequences, using the human and coelacanth sequences as outgroups in all analyses. The topologies of the trees of the α-and β-chains are highly congruent and indicate a specific pattern of gene duplications and gene expression of teleost hemoglobins which has not yet been investigated into more detail. The Notothenioid fish generally contain a single major hemoglobin and often a second minor component. The α- and β-chains of the major components form a monophyletic group in all investigated trees, with the nonantarctic Pseudaphritis as their sister taxon. The minor chains also are a monophyletic group and form an unresolved cluster with the major chains and the hemoglobins of tuna and red gurnard. The Notothenioid families Nototheniidae and Bathydraconidae appear to be paraphyletic. Received: 26 March 1997 / Accepted: 7 May 1997     

Genomic organization of the extracellular coding region of the human FGFR4 and FLT4 genes: Evolution of the genes encoding receptor tyrosine kinases with immunoglobulin-like domains

Receptor tyrosine kinases with five, seven, and three Ig-like domains in their extracellular region are grouped in subclasses IIIa, IIIb, and IIIc, respectively. Here, we describe the genomic organization of the extracellular coding region of the human FGFR4 (IIIc) and FLT4 (IIIb) genes and compare it to that of the human FGFR1(IIIc), KIT, and FMS (IIIa). The results show that while genes belonging to the same subclass have an identical exon/intron structure in their extracellular coding region—as they do in their intracellular coding region—genes of related subclasses only have a similar exon/intron structure. These results strongly support the hypothesis that the genes of the three subclasses evolved from a common ancestor by duplications involving entire genes, already in pieces. Hypotheses on the origin of introns and on the difference in the number of extracellular Ig-like domains in the three gene subclasses are discussed. Received: 19 August 1996 / Accepted: 2 January 1997