Rates Of Sequence Evolution Research Articles

Regressive evolution of an eye pigment gene in independently evolved eyeless subterranean diving beetles

Regressive evolution, the reduction or total loss of non-functional characters, is a fairly common evolutionary phenomenon in subterranean taxa. However, the genetic basis of regressive evolution is not well understood. Here we investigate the molecular evolution of the eye pigment gene cinnabar in several independently evolved lineages of subterranean water beetles using maximum likelihood analyses. We found that in eyeless lineages cinnabar has an increased rate of sequence evolution, as well as mutations leading to frame shifts and stop codons, indicative of pseudogenes. These results are consistent with the hypothesis that regressive evolution of eyes proceeds by random mutations, in the absence of selection, that ultimately lead to the loss of gene function in protein-coding genes specific to the eye pathway.

The roles of positive and negative selection in the molecular evolution of insect endosymbionts

The evolutionary rate acceleration observed in most endosymbiotic bacteria may be explained by higher mutation rates, changes in selective pressure, and increased fixation of deleterious mutations by genetic drift. Here, we explore the forces influencing molecular evolution in Blochmannia, an obligate endosymbiont of Camponotus and related ant genera. Our goals were to compare rates of sequence evolution in Blochmannia with related bacteria, to explore variation in the strength and efficacy of negative (purifying) selection, and to evaluate the effect of positive selection. For six Blochmannia pairs, plus Buchnera and related enterobacteria, estimates of sequence divergence at four genes confirm faster rates of synonymous evolution in the ant mutualist. This conclusion is based on higher d S between Blochmannia lineages despite their more recent divergence. Likewise, generally higher dN in Blochmannia indicates faster rates of nonsynonymous substitution in this group. One exception is the groEL gene, for which lower d N and d N / d S compared to Buchnera indicate exceptionally strong negative selection in Blochmannia. In addition, we explored evidence for positive selection in Blochmannia using both site-and lineage-based maximum likelihood models. These approaches confirmed heterogeneity of d N / d S among codon sites and revealed significant variation in d N / d S across Blochmannia lineages for three genes. Lineage variation affected genes independently, with no evidence of parallel changes in d N / d S across genes along a given branch. Our data also reveal instances of d N / d S greater than one; however, we do not interpret these large d N / d S ratios as evidence for positive selection. In sum, while drift may contribute to an overall rate acceleration at nonsynonymous sites in Blochmannia, variable selective pressures best explain the apparent gene-specific changes in d N / d S across lineages of this ant mutualist. In the course of this study, we reanalyzed variation at Buchnera groEL and found no evidence of positive selection that was previously reported.

Reduced Rates of Sequence Evolution of Y-Linked Satellite DNA in Rumex (Polygonaceae)

One characteristic of sex chromosomes is the accumulation of a set of different types of repetitive DNA sequences in the Y chromosomes. However, little is known about how this occurs or about how the absence of recombination affects the subsequent evolutionary fate of the repetitive sequences in the Y chromosome. Here we compare the evolutionary pathways leading to the appearance of three different families of satellite-DNA sequences within the genomes of Rumex acetosa and R. papillaris, two dioecious plant species with a complex XX/XY(1)Y(2) sex-chromosome system. We have found that two of these families, one autosomic (the RAE730 family) and one Y-linked (the RAYSI family), arose independently from the ancestral duplication of the same 120-bp repeat unit. Conversely, a comparative analysis of the three satellite-DNA families reveals no evolutionary relationships between these two and the third, RAE180, also located in the Y chromosomes. However, we have demonstrated that, regardless of the mechanisms that gave rise to these families, satellite-DNA sequences have different evolutionary fates according to their location in different types of chromosomes. Specifically, those in the Y chromosomes have evolved at half the rate of those in the autosomes, our results supporting the hypothesis that satellite DNAs in nonrecombining Y chromosomes undergo lower rates of sequence evolution and homogenization than do satellite DNAs in autosomes.

Genome evolution in Microsporidia

Microsporidia are eukaryotic intracellular parasites that evolved from fungi. However, they are highly derived and specialized, and lack several typical eukaryotic features, such as canonical mitochondria, flagella and peroxisomes. They also exhibit seemingly “primitive” traits, probably due to their extreme sequence divergence and reductive evolution, which long obscured their true phylogenetic affinity. Gene content and genomic architecture in Microsporidia are known from E. cuniculi, the only species whose complete genome has been sequenced. A genomic sequence survey of A. locustae, a second and distantly related species, carried out in our laboratory provided valuable data for many studies that improved our knowledge on the biology and evolution of Microsporidia. We present here the latest results of our survey, emphasizing comparative analyses between A. locustae and E. cuniculi in order to explore the dynamics of genome evolution and compaction in this intriguing group of parasites. We found that, in spite of their extreme rates of sequence evolution, microsporidian genome structure evolves slowly compared to other eukaryotic groups and we propose that it is a consequence of compaction. We also report the discovery and characterization of a photolyase in A. locustae that protects the spore from DNA damage. In addition, comparative analysis allows the identification of several important genes that would otherwise not be recognised due to the extreme level of sequence evolution. Finally, we compared intergenic DNA between both species aiming to understand the structure of regulatory sequences in these highly compacted genomes.

Rate variation in parasitic plants: correlated and uncorrelated patterns among plastid genes of different function

BackgroundThe analysis of synonymous and nonsynonymous rates of DNA change can help in the choice among competing explanations for rate variation, such as differences in constraint, mutation rate, or the strength of genetic drift. Nonphotosynthetic plants of the Orobanchaceae have increased rates of DNA change. In this study 38 taxa of Orobanchaceae and relatives were used and 3 plastid genes were sequenced for each taxon.ResultsPhylogenetic reconstructions of relative rates of sequence evolution for three plastid genes (rbcL, matK and rps2) show significant rate heterogeneity among lineages and among genes. Many of the non-photosynthetic plants have increases in both synonymous and nonsynonymous rates, indicating that both (1) selection is relaxed, and (2) there has been a change in the rate at which mutations are entering the population in these species. However, rate increases are not always immediate upon loss of photosynthesis. Overall there is a poor correlation of synonymous and nonsynonymous rates. There is, however, a strong correlation of synonymous rates across the 3 genes studied and the lineage-speccific pattern for each gene is strikingly similar. This indicates that the causes of synonymous rate variation are affecting the whole plastid genome in a similar way. There is a weaker correlation across genes for nonsynonymous rates. Here the picture is more complex, as could be expected if there are many causes of variation, differing from taxon to taxon and gene to gene.ConclusionsThe distinctive pattern of rate increases in Orobanchaceae has at least two causes. It is clear that there is a relaxation of constraint in many (though not all) non-photosynthetic lineages. However, there is also some force affecting synonymous sites as well. At this point, it is not possible to tell whether it is generation time, speciation rate, mutation rate, DNA repair efficiency or some combination of these factors.

Comparison of the chicken and turkey genomes reveals a higher rate of nucleotide divergence on microchromosomes than macrochromosomes

A distinctive feature of the avian genome is the large heterogeneity in the size of chromosomes, which are usually classified into a small number of macrochromosomes and numerous microchromosomes. These chromosome classes show characteristic differences in a number of interrelated features that could potentially affect the rate of sequence evolution, such as GC content, gene density, and recombination rate. We studied the effects of these factors by analyzing patterns of nucleotide substitution in two sets of chicken-turkey sequence alignments. First, in a set of 67 orthologous introns, divergence was significantly higher in microchromosomes (chromosomes 11-38; 11.7% divergence) than in both macrochromosomes (chromosomes 1-5; 9.9% divergence; P = 0.016) and intermediate-sized chromosomes (chromosomes 6-10; 9.5% divergence; P = 0.026). At least part of this difference was due to the higher incidence of CpG sites on microchromosomes. Second, using 155 orthologous coding sequences we noted a similar pattern, in which synonymous substitution rates on microchromosomes (13.1%) were significantly higher than were rates on macrochromosomes (10.3%; P = 0.024). Broadly assuming neutrality of introns and synonymous sites, or constraints on such sequences do not differ between chromosomal classes, these observations imply that microchromosomal genes are exposed to more germ line mutations than those on other chromosomes. We also find that dN/dS ratios for genes located on microchromosomes (average, 0.094) are significantly lower than those of macrochromosomes (average, 0.185; P = 0.025), suggesting that the proteins of genes on microchromosomes are under greater evolutionary constraint.

Divergence of conserved non-coding sequences: rate estimates and relative rate tests.

In many eukaryotic genomes only a small fraction of the DNA codes for proteins, but the non-protein coding DNA harbors important genetic elements directing the development and the physiology of the organisms, like promoters, enhancers, insulators, and micro-RNA genes. The molecular evolution of these genetic elements is difficult to study because their functional significance is hard to deduce from sequence information alone. Here we propose an approach to the study of the rate of evolution of functional non-coding sequences at a macro-evolutionary scale. We identify functionally important non-coding sequences as Conserved Non-Coding Nucleotide (CNCN) sequences from the comparison of two outgroup species. The CNCN sequences so identified are then compared to their homologous sequences in a pair of ingroup species, and we monitor the degree of modification these sequences suffered in the two ingroup lineages. We propose a method to test for rate differences in the modification of CNCN sequences among the two ingroup lineages, as well as a method to estimate their rate of modification. We apply this method to the full sequences of the HoxA clusters from six gnathostome species: a shark, Heterodontus francisci; a basal ray finned fish, Polypterus senegalus; the amphibian, Xenopus tropicalis; as well as three mammalian species, human, rat and mouse. The results show that the evolutionary rate of CNCN sequences is not distinguishable among the three mammalian lineages, while the Xenopus lineage has a significantly increased rate of evolution. Furthermore the estimates of the rate parameters suggest that in the stem lineage of mammals the rate of CNCN sequence evolution was more than twice the rate observed within the placental amniotes clade, suggesting a high rate of evolution of cis-regulatory elements during the origin of amniotes and mammals. We conclude that the proposed methods can be used for testing hypotheses about the rate and pattern of evolution of putative cis-regulatory elements.

Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution.

To study the evolutionary effects of polyploidy on plant gene functions, we analyzed functional genomics data for a large number of duplicated gene pairs formed by ancient polyploidy events in Arabidopsis thaliana. Genes retained in duplicate are not distributed evenly among Gene Ontology or Munich Information Center for Protein Sequences functional categories, which indicates a nonrandom process of gene loss. Genes involved in signal transduction and transcription have been preferentially retained, and those involved in DNA repair have been preferentially lost. Although the two members of each gene pair must originally have had identical transcription profiles, less than half of the pairs formed by the most recent polyploidy event still retain significantly correlated profiles. We identified several cases where groups of duplicated gene pairs have diverged in concert, forming two parallel networks, each containing one member of each gene pair. In these cases, the expression of each gene is strongly correlated with the other nonhomologous genes in its network but poorly correlated with its paralog in the other network. We also find that the rate of protein sequence evolution has been significantly asymmetric in >20% of duplicate pairs. Together, these results suggest that functional diversification of the surviving duplicated genes is a major feature of the long-term evolution of polyploids.

Testing the Phylogenetic Position of a Parasitic Plant (Cuscuta, Convolvulaceae, Asteridae): Bayesian Inference and the Parametric Bootstrap on Data Drawn from Three Genomes

Previous findings on structural rearrangements in the chloroplast genome of Cuscuta (dodder), the only parasitic genus in the morning-glory family, Convolvulaceae, were attributed to its parasitic life style, but without proper comparison to related nonparasitic members of the family. Before molecular evolutionary questions regarding genome evolution can be answered, the phylogenetic problems within the family need to be resolved. However, the phylogenetic position of parasitic angiosperms and their precise relationship to nonparasitic relatives are difficult to infer. Problems are encountered with both morphological and molecular evidence. Molecular data have been used in numerous studies to elucidate relationships of parasitic taxa, despite accelerated rates of sequence evolution. To address the question of the position of the genus Cuscuta within Convolvulaceae, we generated a new molecular data set consisting of mitochondrial (atpA) and nuclear (RPB2) genes, and analyzed these data together with an existing chloroplast data matrix (rbcL, atpB, trnL-F, and psbE-J), to which an additional chloroplast gene (rpl2) was added. This data set was analyzed with an array of phylogenetic methods, including Bayesian analysis, maximum likelihood, and maximum parsimony. Further exploration of data was done by using methods of phylogeny hypothesis testing. At least two nonparasitic lineages are shown to diverge within the Convolvulaceae before Cuscuta. However, the exact sister group of Cuscuta could not be ascertained, even though many alternatives were rejected with confidence. Caution is therefore warranted when interpreting the causes of molecular evolution in Cuscuta. Detailed comparisons with nonparasitic Convolvulaceae are necessary before firm conclusions can be reached regarding the effects of the parasitic mode of life on patterns of molecular evolution in Cuscuta.

Genome Compaction and Stability in Microsporidian Intracellular Parasites

Microsporidian genomes are extraordinary among eukaryotes for their extreme reduction: although they are similar in form to other eukaryotic genomes, they are typically smaller than many prokaryotic genomes [1, 2]. At the same time, their rates of sequence evolution are among the highest for eukaryotic organisms [3]. To explore the effects of compaction on nuclear genome evolution, we sequenced 685,000 bp of the Antonospora locustae genome (formerly Nosema locustae) and compared its organization with the recently completed genome of the human parasite Encephalitozoon cuniculi[1, 2]. Despite being very distantly related, the genomes of these two microsporidian species have retained an unexpected degree of synteny: 13% of genes are in the same context, and 30% of the genes were separated by a small number of short rearrangements. Microsporidian genomes are, therefore, paradoxically composed of rapidly evolving sequences harbored within a slowly evolving genome, although these two processes are sometimes considered to be coupled [4–7]. Microsporidian genomes show that eukaryotic genomes (like genes) do not evolve in a clock-like fashion, and genome stability may result from compaction in addition to a lack of recombination, as has been traditionally thought to occur in bacterial and organelle genomes [8–11].

Molecular phylogeny and divergence time of the Antarctic sea urchin (Sterechinus neumayeri) in relation to the South American sea urchins

Sterechinus neumayeri is an abundant regular sea urchin that lives in the shallow Antarctic waters. This organism has been used as a model system in many fields of the Antarctic biology. Yet, understanding of the evolutionary identity of the species, such as its phylogenetic relationships and divergence time, remains limited. Here, we reconstructed the molecular phylogenies of the species together with two sea urchin species in southernmost South America (Loxechinus albus and Pseudechinus magellanicus), a parechinid species (Paracentrotus lividus) and three strongylocentrotid species (Strongylocentrotus purpuratus, S. intermedius, and Hemicentrotus pulcherrimus) using mitochondrial DNA sequences of 12S rDNA-tRNA(gln) region (877nt) and cytochrome oxidase subunit I (COI, 1079nt). The rate of sequence evolution and the divergence time of the species were then estimated from the trees. The phylogenetic trees reveal that S. neumayeri is a sister group to the lineage of L. albus and P. lividus, and separated from the lineage 24–35 million years ago (m.y.a.). The divergence between S. neumayeri and L. albus coincides with the separation of Antarctica from South America, suggesting that the tectonic event must have provoked the cladogenesis of the species through vicariance.

Molecular cloning of the estrogen and progesterone receptors of the American alligator

Steroid hormones perform many essential roles in vertebrates during embryonic development, reproduction, growth, water balance, and responses to stress. The estrogens are essential for normal reproductive activity in female and male vertebrates and appear to have direct actions during sex determination in some vertebrates. To begin to understand the molecular mechanisms of estrogen action in alligators, we have isolated cDNAs encoding the estrogen receptors (ER) from the ovary. Degenerate PCR primers specific to ER were designed and used to amplify alligator ovary RNA. Two different DNA fragments (ERα and ERβ) were obtained and the full-length alligator ERα cDNA was obtained using 5 ′ and 3 ′ RACE. The inferred amino acid sequence of alligator ERα (aERα) was very similar to the chicken ERα (91% identity), although phylogenetic analyses suggested profound differences in the rate of sequence evolution for vertebrate ER sequences. We also isolated partial DNA fragments encoding ERβ and the progesterone receptor (PR) of the alligator, both of which show strong sequence similarities to avian ERβ and PR. We examined the expression levels of these three steroid receptors (ERα, ERβ, and PR) in the ovary of juvenile alligators and observed detectable levels of all three receptors. Quantitative RT-PCR showed that gonadal ERα transcript levels in juvenile alligators decreased after E 2 treatment whereas ERβ and PR transcripts were not changed. These results provide tools that will allow future studies examining the regulation and ontogenic expression of steroid receptors in alligators and expand our knowledge of vertebrate steroid receptor evolution.

The complete mitochondrial DNA sequence of the green alga Pseudendoclonium akinetum (Ulvophyceae) highlights distinctive evolutionary trends in the chlorophyta and suggests a sister-group relationship between the Ulvophyceae and Chlorophyceae.

The mitochondrial genome has undergone radical changes in both the Chlorophyta and Streptophyta, yet little is known about the dynamics of mtDNA evolution in either of these lineages. In the Chlorophyta, which comprises four of the five recognized classes of green algae (Prasinophyceae, Trebouxiophyceae, Ulvophyceae, and Chlorophyceae), the mitochondrial genome varies from 16 to 55 kb. This genome has retained a compact gene organization and a relatively complex gene repertoire ("ancestral" pattern) in the basal lineages represented by the Trebouxiophyceae and Prasinophyceae, whereas it has been reduced in size and gene complement and tends to evolve much more rapidly at the sequence level ("reduced-derived" pattern of evolution) in the Chlorophyceae and the lineage leading to the enigmatic chlorophyte Pedinomonas. To gain information about the evolutionary trends of mtDNA in the Ulvophyceae and also to gain insights into the phylogenetic relationships between ulvophytes and other chlorophytes, we have determined the mtDNA sequence of Pseudendoclonium akinetum. At 95,880 bp, Pseudendoclonium mtDNA is the largest green-algal mitochondrial genome sequenced to date and has the lowest gene density. These derived features are reminiscent of the "expanded" pattern exhibited by embryophyte mtDNAs, indicating that convergent evolution towards genome expansion has occurred independently in the Chlorophyta and Streptophyta. With 57 conserved genes, the gene repertoire of Pseudendoclonium mtDNA is slightly smaller than those of the prasinophyte Nephroselmis olivacea and the trebouxiophyte Prototheca wickerhamii. This ulvophyte mtDNA contains seven group I introns, four of which have homologs in green-algal mtDNAs displaying an "ancestral" or a "reduced-derived" pattern of evolution. Like its counterpart in the chlorophycean green alga Scenedesmus obliquus, it features numerous small, dispersed repeats in intergenic regions and introns. Its overall rate of sequence evolution appears to be accelerated to an intermediary level as compared with the rates observed in "ancestral" and "reduced-derived" mtDNAs. In agreement with the finding that Pseudendoclonium mtDNA exhibits features typical of both the "ancestral" and "reduced-derived" patterns of evolution, phylogenetic analyses of seven mtDNA-encoded proteins revealed a sister-group relationship between this ulvophyte and chlorophytes displaying "reduced-derived" mtDNAs.

A conservative test of genetic drift in the endosymbiotic bacterium Buchnera: slightly deleterious mutations in the chaperonin groEL.

The obligate endosymbiotic bacterium Buchnera aphidicola shows elevated rates of sequence evolution compared to free-living relatives, particularly at nonsynonymous sites. Because Buchnera experiences population bottlenecks during transmission to the offspring of its aphid host, it is hypothesized that genetic drift and the accumulation of slightly deleterious mutations can explain this rate increase. Recent studies of intraspecific variation in Buchnera reveal patterns consistent with this hypothesis. In this study, we examine inter- and intraspecific nucleotide variation in groEL, a highly conserved chaperonin gene that is constitutively overexpressed in Buchnera. Maximum-likelihood estimates of nonsynonymous substitution rates across Buchnera species are strikingly low at groEL compared to other loci. Despite this evidence for strong purifying selection on groEL, our intraspecific analysis of this gene documents reduced synonymous polymorphism, elevated nonsynonymous polymorphism, and an excess of rare alleles relative to the neutral expectation, as found in recent studies of other Buchnera loci. Comparisons with Escherichia coli generally show patterns predicted by their differences in N(e). The sum of these observations is not expected under relaxed or balancing selection, selective sweeps, or increased mutation rate. Rather, they further support the hypothesis that drift is an important force driving accelerated protein evolution in this obligate mutualist.

Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution.

Lineage-specific gene loss, to a large extent, accounts for the differences in gene repertoires between genomes, particularly among eukaryotes. We derived a parsimonious scenario of gene losses for eukaryotic orthologous groups (KOGs) from seven complete eukaryotic genomes. The scenario involves substantial gene loss in fungi, nematodes, and insects. Based on this evolutionary scenario and estimates of the divergence times between major eukaryotic phyla, we introduce a numerical measure, the propensity for gene loss (PGL). We explore the connection among the propensity of a gene to be lost in evolution (PGL value), protein sequence divergence, the effect of gene knockout on fitness, the number of protein-protein interactions, and expression level for the genes in KOGs. Significant correlations between PGL and each of these variables were detected. Genes that have a lower propensity to be lost in eukaryotic evolution accumulate fewer substitutions in their protein sequences and tend to be essential for the organism viability, tend to be highly expressed, and have many interaction partners. The dependence between PGL and gene dispensability and interactivity is much stronger than that for sequence evolution rate. Thus, propensity of a gene to be lost during evolution seems to be a direct reflection of its biological importance.

Phylogeny of “ Philoceanus complex” seabird lice (Phthiraptera: Ischnocera) inferred from mitochondrial DNA sequences

The Philoceanus complex is a large assemblage of lice that parasitise procellariiform seabirds (petrels, albatrosses, and their relatives). We obtained mitochondrial 12S rRNA and cytochrome oxidase I DNA sequences from 39 species from diverse hosts and localities. Resolution of deeper relationships between genera was limited, however there is evidence for two major clades, one hosted by albatrosses, the other by petrels. Based on our results, the genera hosted by albatrosses are excellent candidates for detailed analysis of cospeciation. Our results also suggest that a previous estimate of a 5-fold difference in the relative rate of sequence evolution in lice and their avian hosts is an artefact of limited taxonomic sampling.

Geobiological analysis using whole genome‐based tree building applied to the Bacteria, Archaea, and Eukarya

ABSTRACTWe constructed genomic trees based on the presence and absence of families of protein‐encoding genes observed in 55 prokaryotic and five eukaryotic genomes. There are features of the genomic trees that are not congruent with typical rRNA phylogenetic trees. In the bacteria, for example, Deinococcus radiodurans associates with the Gram‐positive bacteria, a result that is also seen in some other phylogenetic studies using whole genome data. In the Archaea, the methanogens plus Archaeoglobus form a united clade and the Euryarchaeota are divided with the two Thermoplasma genomes and Halobacterium sp. falling below the Crenarchaeota. While the former appears to be an accurate representation of methanogen‐relatedness, the misplacement of Halobacterium may be an artefact of parsimony. These results imply the last common ancestor of the Archaea was not a methanogen, leaving sulphur reduction as the most geochemically plausible metabolism for the base of the archaeal crown group. It also suggests that methanogens were not a component of the Earth's earliest biosphere and that their origin occurred sometime during the Archean. In the Eukarya, the parsimony analysis of five Eukaryotes using the Crenarchaeota as an outgroup seems to counter the Ecdysozoa hypothesis, placing Caenorhabditis elegans (Nematoda) below the common ancestor of Drosophila melanogaster (Arthropoda) and Homo sapiens (Chordata) even when efforts are made to counter the possible effects of a faster rate of sequence evolution for the C. elegans genome. Further analysis, however, suggests that the gene loss of ‘animal’ genes is highest in C. elegans and is obscuring the relationships of these organisms.

Congruent avian phylogenies inferred from mitochondrial and nuclear DNA sequences.

Recent molecular studies addressing the phylogenetic relationships of avian orders have had conflicting results. While studies using nuclear DNA sequences tend to support traditional taxonomic views, also supported by morphological data [(paleognaths (galloanseres (all other birds)))], with songbirds forming a clade within Neoaves (all other birds), analyses with complete mtDNA genomes have resulted in topologies that place songbirds as one of the earliest-diverging avian lineages. Considering that over half of the extant bird species are songbirds, these different results have very different implications for our understanding of avian evolution. We analyzed data sets comprising nearly 4 kb of mitochondrial DNA (mtDNA) (complete 12S, ND1, ND2, and cytochrome b) plus 600 bp of the nuclear gene c-mos for 15 birds that were chosen to represent all major avian clades and to minimize potential long-branch attraction problems; we used a partition-specific maximum likelihood approach. Our results show congruence with respect to the ingroup among phylogenies obtained with mtDNA and the nuclear gene c-mos, separately or combined. The data sets support a traditional avian taxonomy, with paleognaths (ratites and tinamous) occupying a basal position and with songbirds more derived and forming a monophyletic group. We also show that, for mtDNA studies, turtles may be a better outgroup for birds than crocodilians because of their slower rate of sequence evolution.

Linking dynamical and population genetic models of persistent viral infection.

This article develops a theoretical framework to link dynamical and population genetic models of persistent viral infection. This linkage is useful because, while the dynamical and population genetic theories have developed independently, the biological processes they describe are completely interrelated. Parameters of the dynamical models are important determinants of evolutionary processes such as natural selection and genetic drift. We develop analytical methods, based on coupled differential equations and Markov chain theory, to predict the accumulation of genetic diversity within the viral population as a function of dynamical parameters. These methods are first applied to the standard model of viral dynamics and then generalized to consider the infection of multiple host cell types by the viral population. Each cell type is characterized by specific parameter values. Inclusion of multiple cell types increases the likelihood of persistent infection and can increase the amount of genetic diversity within the viral population. However, the overall rate of gene sequence evolution may actually be reduced.

Molecular Phylogenies and Divergence Times of Sea Urchin Species of Strongylocentrotidae, Echinoida

Sea urchins of the family Strongylocentrotidae have been important model systems in many fields of basic biology, yet knowledge of their evolutionary identities such as the phylogenetic relationships and divergence times remains limited. Here, I inferred molecular phylogenies of seven Strongylocentrotid species (Strongylocentrotus franciscanus, S. nudus, S. purpuratus, S. intermedius, S. droebachiensis, S. pallidus, and Hemicentrotus pulcherrimus) from the analyses of mitochondrial DNA sequences of 12SrDNA (349 nt), 12SrDNA-tRNA(gln) region (862 nt), and a combined sequence of cytochrome oxidase subunit I (COI, 1080 nt) and NADH dehydrogenase subunit I (NDI, 742 nt). The rate of sequence evolution and divergence times for each species were then estimated from the trees with reference to the time of separation between Strongylocentrotidae and Parechinidae, 35 to 50 MYA. The three trees agree well with each other, and the phylogeny is summarized by ((S. franciscanus, S. nudus), (H. pulcherrimus (S. purpuratus, S. intermedius (S. droebachiensis, S. pallidus)))). It is notable that the genus Strongylocentrotus consists of two distinct clades and that H. pulcherrimus branches off within Strongylocentrotus, implying assignment of a separate, monospecific genus to this species inappropriate. The rate of sequence evolution is calibrated to be 0.24%-0.34%/Myr in 12SrDNA, 0.25%-0.36%/Myr in 12SrDNA-tRNA(gln), and 0.65%-0.93%/Myr in COI-NDI combined sequences. S. purpuratus, in particular, shows the significantly higher rate of evolution in the 12SrDNA and 12SrDNA-tRNA(gln) regions compared to other species, suggesting careful use of its sequences in comparative studies. The two clades of Strongylocentrotidae seem to have split 13-19 MYA, and H. pulcherrimus branched off 7.2-14 MYA. In the former clade, S. franciscanus and S. nudus separated 5.7-8.1 MYA. In the latter clade, S. purpuratus, S. intermedius, and the clade of S. droebachiensis and S. pallidus diverged approximately 4.6-12 MYA, and the last two closest species separated 2.1-3.1 MYA.