Xanthine Dehydrogenase (XDH): Episodic Evolution of a ``Neutral' Protein

We investigated the evolution of xanthine dehydrogenase (Xdh) in 34 species from the three multicellular kingdoms, including one plant, two fungi, and three animal phyla, two classes of vertebrates, four orders of mammals, and two orders of insects. We adopted a model-based maximum-likelihood framework of inference. After accounting for among-site rate variation and heterogeneous nucleotide composition of the sequences using the discrete gamma distribution, and using nonhomogeneous nonstationary representations of the substitution process, the rate of amino acid replacement is 30.4 x 10(-10)/site/year when Drosophila species are compared but only approximately 18 x 10(-10)/site/year when comparisons are made between mammal orders, between insect orders, or between different animal phyla and approximately 11 x 10(-10)/site/year when comparisons are made between birds and mammals, between fungi, or between the three multicellular kingdoms. To account for these observations, the rate of amino acid replacement must have been eight or more times higher in some lineages and at some times than in others. Spastic evolution of Xdh appears to be related to the particularities of the genomes in which the locus is embedded.     

Exploring metazoan evolution through dynamic and holistic changes in protein families and domains

ABSTRACT: BACKGROUND: Proteins convey the majority of biochemical and cellular activities in organisms. Over the course of evolution, proteins undergo normal sequence mutations as well as large scale mutations involving domain duplication and/or domain shuffling. These events result in the generation of new proteins and protein families. Processes that affect proteome evolution drive species diversity and adaptation. Herein, change over the course of metazoan evolution, as defined by birth/death and duplication/deletion events within protein families and domains, was examined using the proteomes of 9 metazoan and two outgroup species. RESULTS: In studying members of the three major metazoan groups, the vertebrates, arthropods, and nematodes, we found that the number of protein families increased at the majority of lineages over the course of metazoan evolution where the magnitude of these increases was greatest at the lineages leading to mammals. In contrast, the number of protein domains decreased at most lineages and at all terminal lineages. This resulted in a weak correlation between protein family birth and domain birth; however, the correlation between domain birth and domain member duplication was quite strong. These data suggest that domain birth and protein family birth occur via different mechanisms, and that domain shuffling plays a role in the formation of protein families. The ratio of protein family birth to protein domain birth (domain shuffling index) suggests that shuffling had a more demonstrable effect on protein families in nematodes and arthropods than in vertebrates. Through the contrast of high and low domain shuffling indices at the lineages of Trichinella spiralis and Gallus gallus, we propose a link between protein redundancy and evolutionary changes controlled by domain shuffling; however, the speed of adaptation among the different lineages was relatively invariant. Evaluating the functions of protein families that appeared or disappeared at the last common ancestors (LCAs) of the three metazoan clades supports a correlation with organism adaptation. Furthermore, bursts of new protein families and domains in the LCAs of metazoans and vertebrates are consistent with whole genome duplications. CONCLUSION: Metazoan speciation and adaptation were explored by birth/death and duplication/deletion events among protein families and domains. Our results provide insights into protein evolution and its bearing on metazoan evolution.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

Switch in codon bias and increased rates of amino acid substitution in the Drosophila saltans species group.

We investigated the nucleotide composition of five genes, Xdh, Adh, Sod, Per, and 28SrRNA, in nine species of Drosophila (subgenus Sophophora) and one of Scaptodrosophila. The six species of the Drosophila saltans group markedly differ from the others in GC content and codon use bias. The GC content in the third codon position, and to a lesser extent in the first position and the introns, is higher in the D. melanogaster and D. obscura groups than in the D. saltans group (in Scaptodrosophila it is intermediate but closer to the melanogaster and obscura species). Differences are greater for Xdh than for Adh, Sod, Per, and 28SrRNA, which are functionally more constrained. We infer that rapid evolution of GC content in the saltans lineage is largely due to a shift in mutation pressure, which may have been associated with diminished natural selection due to smaller effective population numbers rather than reduced recombination rates. The rate of GC content evolution impacts the rate of protein evolution and may distort phylogenetic inferences. Previous observations suggesting that GC content evolution is very limited in Drosophila may have been distorted due to the restricted number of genes and species (mostly D. melanogaster) investigated.     

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.     

A multidisciplinary approach reveals cryptic diversity in Western Palearctic Tetramorium ants (Hymenoptera: Formicidae)

Diversity of ants of the Tetramorium caespitum/impurum complex was investigated in a multidisciplinary study. Focusing on morphologically hardly distinguishable Western Palearctic samples, we demonstrate the genetic and phenotypic diversity, demarcate phylogenetic entities, and discuss the clades in terms of biogeography. Sequences of 1113bp of the mitochondrial COI gene revealed 13 lineages. COII data, worker morphometry and male genitalia morphology corroborated the COI results for seven lineages; the remaining six were disregarded because of small sample size. A comparison with published data on cuticular hydrocarbons showed correspondence. The seven entities show different distribution patterns, though some ranges overlap in Central Europe. Since no major discrepancy between the results of the different disciplines became apparent, we conclude that the seven entities within the T. caespitum/impurum complex represent seven species. Geographical evidence allows the identification of T. caespitum and T. impurum, and we therefore designate neotypes and redescribe the two species in terms of morphology and mtDNA. As the revision of about 50 taxon names would go beyond the scope of this study, we refer to the remaining five species under code names. We discuss our findings in terms of plesiomorphy and convergent evolution by visualizing the mtDNA phylogeny in morphological space.     

Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

R1 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of arthropods. One aspect of R1 evolution that has been difficult to explain is the presence of divergent lineages of R1 in the rDNA loci of the same species. Multiple lineages should compete for a limited number of insertion sites, in addition to being subject to the concerted evolution processes homogenizing the rRNA genes. The presence of multiple lineages suggests either the ability of the elements to overcome these factors and diverge within rDNA loci, or the introduction of new lineages by horizontal transmission. To address this issue, we attempted to characterize the complete set of R1 elements in the rDNA locus from five Drosophila species groups (melanogaster, obscura, testacea, quinaria, and repleta). Two major R1 lineages, A and B, that diverged about 100 MYA were found to exist in Drosophila. Elements of the A lineage were found in all 35 Drosophila species tested, while elements of the B lineage were found in only 11 species from three species groups. Phylogenetic analysis of the R1 elements, supported by comparison of their rates of nucleotide sequence substitution, revealed that both the A and the B lineages have been maintained by vertical descent. The B lineage was less stable and has undergone numerous, independent elimination events, while the A lineage has diverged into three sublineages, which were, in turn, differentially stable. We conclude that while the differential retention of multiple lineages greatly complicates its phylogenetic history, the available R1 data continue to be consistent with the strict vertical descent of these elements.     

Comparative phylogeography of Trypanosoma rangeli and Rhodnius (Hemiptera: Reduviidae) supports a long coexistence of parasite lineages and their sympatric vectors

To make reliable interpretations about evolutionary relationships between Trypanosoma rangeli lineages and their insect vectors (triatomine bugs of the genus Rhodnius) and, thus, about the determinant factors of lineage segregation within T. rangeli, we compared phylogenies of parasite isolates and vector species. Sixty-one T. rangeli isolates from invertebrate and vertebrate hosts were initially evaluated in terms of polymorphism of the spliced-leader gene (SL). Further analysis based on SL and SSUrRNA sequences from 33 selected isolates, representative of the overall phylogenetic diversity and geographical range of T. rangeli, supported four phylogenetic lineages within this species. By comparing the phylogeny of Rhodnius species with that inferred for T. rangeli isolates and through analysis of the geographical range of the isolates, we showed that there is a very significant overlap in the distribution of Rhodnius species and T. rangeli lineages. Congruence between phylogeographical analysis of both T. rangeli lineages and complexes of Rhodnius species are consistent with the hypothesis of a long coexistence of parasites and their vectors, with lineage divergence associated with sympatric species of Rhodnius apparently without association with particular vertebrate hosts. Separation of T. rangeli isolates from vectors of distinct complexes living in sympatry favours the absence of gene flow between the lineages and suggests evolution of T. rangeli lineages in independent transmission cycles, probably associated to specific Rhodnius spp. ecotopes. A polymerase chain reaction assay based on SL intergenic sequences was developed for simultaneous identification and lineage genotyping of T. rangeli in epidemiological surveys.     

Evolution of circulating wild poliovirus and of vaccine-derived poliovirus in an immunodeficient patient: a unifying model

We determined nucleotide sequences of the VP1 and 2AB genes and portions of the 2C and 3D genes of two evolving poliovirus lineages: circulating wild viruses of T geotype and Sabin vaccine-derived isolates from an immunodeficient patient. Different regions of the viral RNA were found to evolve nonsynchronously, and the rate of evolution of the 2AB region in the vaccine-derived population was not constant throughout its history. Synonymous replacements occurred not completely randomly, suggesting the need for conservation of certain rare codons (possibly to control translation elongation) and the existence of unidentified constraints in the viral RNA structure. Nevertheless the major contribution to the evolution of the two lineages came from linear accumulation of synonymous substitutions. Therefore, in agreement with current theories of viral evolution, we suggest that the majority of the mutations in both lineages were fixed as a result of successive sampling, from the heterogeneous populations, of random portions containing predominantly neutral and possibly adverse mutations. As a result of such a mode of evolution, the virus fitness may be maintained at a more or less constant level or may decrease unless more-fit variants are stochastically generated. The proposed unifying model of natural poliovirus evolution has important implications for the epidemiology of poliomyelitis.     

Structure and function of lineage-specific sequence insertions in the bacterial RNA polymerase beta' subunit

The large beta and beta' subunits of the bacterial core RNA polymerase (RNAP) are highly conserved throughout evolution. Nevertheless, large sequence insertions in beta and beta' characterize specific evolutionary lineages of bacteria. The Thermus aquaticus RNAP beta' subunit contains a 283 residue insert between conserved regions A and B that is found in only four bacterial species. The Escherichia coli RNAP beta' subunit contains a 188 residue insert in the middle of conserved region G that is found in a wide range of bacterial species. Here, we present structural studies of these two beta' insertions. We show that the inserts comprise repeats of a previously characterized fold, the sandwich-barrel hybrid motif (as predicted from previous sequence analysis) and that the inserts serve significant roles in facilitating protein/protein and/or protein/nucleic acid interactions.     

Lineage-Specific Differences in the Amino Acid Substitution Process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size.To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.     

Diversification of Neoaves: integration of molecular sequence data and fossils

Patterns of diversification and timing of evolution within Neoaves, which includes almost 95% of all bird species, are virtually unknown. On the other hand, molecular data consistently indicate a Cretaceous origin of many neoavian lineages and the fossil record seems to support an Early Tertiary diversification. Here, we present the first well-resolved molecular phylogeny for Neoaves, together with divergence time estimates calibrated with a large number of stratigraphically and phylogenetically well-documented fossils. Our study defines several well-supported clades within Neoaves. The calibration results suggest that Neoaves, after an initial split from Galloanseres in Mid-Cretaceous, diversified around or soon after the K/T boundary. Our results thus do not contradict palaeontological data and show that there is no solid molecular evidence for an extensive pre-Tertiary radiation of Neoaves.     

Determinants of rate variation in mammalian DNA sequence evolution

Attempts to analyze variation in the rates of molecular evolution among mammalian lineages have been hampered by paucity of data and by nonindependent comparisons. Using phylogenetically independent comparisons, we test three explanations for rate variation which predict correlations between rate variation and generation time, metabolic rate, and body size. Mitochondrial and nuclear genes, protein coding, rRNA, and nontranslated sequences from 61 mammal species representing 14 orders are used to compare the relative rates of sequence evolution. Correlation analyses performed on differences in genetic distance since common origin of each pair against differences in body mass, generation time, and metabolic rate reveal that substitution rate at fourfold degenerate sites in two out of three protein sequences is negatively correlated with generation time. In addition, there is a relationship between the rate of molecular evolution and body size for two nuclear-encoded sequences. No evidence is found for an effect of metabolic rate on rate of sequence evolution. Possible causes of variation in substitution rate between species are discussed.     

Comparative molecular population genetics of the Xdh locus in the cactophilic sibling species Drosophila buzzatii and D. koepferae

The Xdh (rosy) gene is one of the best studied in the Drosophila genus from an evolutionary viewpoint. Here we analyze nucleotide variation in a 1875-bp fragment of the second exon of Xdh in Argentinian populations of the cactophilic D. buzzatii and its sibling D. koepferae. The major electrophoretic alleles of D. buzzatii not only lack diagnostic amino acids in the region studied but also differ on average from each other by four to 13 amino acid changes. Our data also suggest that D. buzzatii populations belonging to different phytogeographic regions are not genetically differentiated, whereas D. koepferae exhibits a significant pattern of population structure. The Xdh region studied is twice as polymorphic in D. buzzatii as in D. koepferae. Differences in historical population size or in recombinational environment between species could account for the differences in the level of nucleotide variation. In both species, the Xdh region exhibits a great number of singletons, which significantly departs from the frequency spectrum expected under neutrality for nonsynonymous sites and also for synonymous sites in D. buzzatii. These excesses of singletons could be the signature of a recent population expansion in D. buzzatii, whereas they may be simply explained as the result of negative selection in D. koepferae.     

Nucleotide polymorphism at the xanthine dehydrogenase locus in Drosophila pseudoobscura.

Sequential polyacrylamide electrophoresis has revealed 20 allozymes of xanthine dehydrogenase (XDH) in Drosophila pseudoobscura. DNA sequence determination of seven isolates of the Xdh locus that represent six allozyme classes are presented here. Of the 5,456 sites examined, 180 are polymorphic, with 27 polymorphisms occurring at nonsynonymous, or replacement, sites. An average of nine amino acids differ between XDH allozyme classes, with 85% of the polymorphic amino acids singly represented. The level and pattern of variation observed at Xdh argue that the effective population size of the species is quite large--i.e., on the order of 2 x 10(6)--and that the populations sampled are quite ancient. In addition, as judged by two statistical tests, the levels of nucleotide polymorphism observed at Xdh are compatible with predictions from the neutral theory of molecular evolution.     

A comparative analysis of the Far Eastern grayling species Thymallus tugarinae and Thymallus grubii flavomaculatus based on the data from mitochondrial DNA cytochrome b gene sequencing

Nucleotide sequences of a fragment of mitochondrial DNA cytochrome b gene were obtained from two species of graylings that inhabit the Russian Far East, viz., Thymallus tugarinae and T. grubii flavomaculatus. A phylogenetic analysis of four Thymallus species, whose relationships have been poorly studied, was performed on these new data as well as on several sequences from the NCBI GenBank sequence database. The rate of genetic divergence between the Far Eastern graylings Thymallus tugarinae and T. g. flavomaculatus corresponded to that between European (T. thymallus) and Arctic (T. arcticus) graylings. Moreover, single nucleotide substitutions that result in alterations of the amino-acid structure of protein products (non-synonymic mutations) were revealed between sequences of the cytochrome b gene in the Far Eastern grayling species. The topology of phylogenetic trees, which was composed by means of Bayesian analysis and the maximum parsimony method, showed four independent phylogenetic lineages of graylings. In addition, the phylogenetic relationships between T. tugarinae and T. g. flavomaculatus were supported statistically. The obtained data indicated that Far Eastern graylings T. tugarinae and T. g. flavomaculatus are distinct species and confirmed that they belong to different phylogenetic lineages, as found earlier.     

AN INTEGRATIVE TEST OF THE DEAD‐END HYPOTHESIS OF SELFING EVOLUTION IN TRITICEAE (POACEAE)

Self‐fertilization is hypothesized to be an evolutionary dead end because reversion to outcrossing can rarely happen, and selfing lineages are thought to rapidly become extinct because of limited potential for adaptation and/or accumulation of deleterious mutations. We tested these two assumptions by combining morphological characters and molecular‐evolution analyses in a tribe of hermaphroditic grasses (Triticeae). First, we determined the mating system of the 19 studied species. Then, we sequenced 27 protein‐coding loci and compared base composition and substitution patterns between selfers and outcrossers. We found that the evolution of the mating system is best described by a model including outcrossing‐to‐selfing transitions only. At the molecular level, we showed that regions of low recombination exhibit signatures of relaxed selection. However, we did not detect any evidence of accumulation of nonsynonymous substitutions in selfers compared to outcrossers. Additionally, we tested for the potential deleterious effects of GC‐biased gene conversion in outcrossing species. We found that recombination and not the mating system affected substitution patterns and base composition. We suggest that, in Triticeae, although recombination patterns have remained stable, selfing lineages are of recent origin and inbreeding may have persisted for insufficient time for differences between the two mating systems to evolve.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.